Open Access

Repurposing non‑invasive prenatal testing data: Population study of single nucleotide variants associated with colorectal cancer and Lynch syndrome

  • Authors:
    • Natalia Forgacova
    • Juraj Gazdarica
    • Jaroslav Budis
    • Jan Radvanszky
    • Tomas Szemes
  • View Affiliations

  • Published online on: September 13, 2021     https://doi.org/10.3892/ol.2021.13040
  • Article Number: 779
  • Copyright: © Forgacova et al. This is an open access article distributed under the terms of Creative Commons Attribution License.

Metrics: Total Views: 0 (Spandidos Publications: | PMC Statistics: )
Total PDF Downloads: 0 (Spandidos Publications: | PMC Statistics: )


Abstract

In our previous work, genomic data generated through non‑invasive prenatal testing (NIPT) based on low‑coverage massively parallel whole‑genome sequencing of total plasma DNA of pregnant women in Slovakia was described as a valuable source of population specific data. In the present study, these data were used to determine the population allele frequency of common risk variants located in genes associated with colorectal cancer (CRC) and Lynch syndrome (LS). Allele frequencies of identified variants were compared with six world populations to detect significant differences between populations. Finally, variants were interpreted, functional consequences were searched for and clinical significance of variants was investigated using publicly available databases. Although the present study did not identify any pathogenic variants associated with CRC or LS in the Slovak population using NIPT data, significant differences were observed in the allelic frequency of risk CRC variants previously reported in genome‑wide association studies and common variants located in genes associated with LS. As Slovakia is one of the leading countries with the highest incidence of CRC among male patients in the world, there is a need for studies dedicated to investigating the cause of such a high incidence of CRC in Slovakia. The present study also assumed that extensive cross‑country data aggregation of NIPT results would represent an unprecedented source of information concerning human genome variation in cancer research.

Introduction

Colorectal cancer (CRC) is the third most common cancer worldwide and the second most common cancer in Europe, causing an estimated 9.4% of all cancer deaths in Europe. CRC is also a serious societal problem in Slovakia, with an incidence rate of 15,7% and a mortality rate of 15,6% [(1,2) GLOBOCAN 2020 data]. Many risk factors and causes are associated with the likelihood of developing CRC, but the main reason is still not fully understood. The considerable geographical variability suggests that CRC is a complex polygenic disease caused by genetic and environmental factors and their interactions. Age, sex, lifestyle, and dietary habits (3,4), including meat and alcohol consumption (5,6), tobacco smoking (79), obesity and lack of physical activity (4,10,11) play a major role in the pathogenesis of CRC. Other well-known risk factors may also be inflammatory bowel diseases, acromegaly, renal transplantation with long-term immunosuppression, diabetes mellitus and insulin resistance, cholecystectomy, or androgen deprivation therapy (3,4).

Beside these, inherited susceptibility plays a significant role in the etiology of CRC because it can be responsible for about 35% of all cases of colorectal cancer. However, high-penetrance germline variants in known genes (APC, BRCA2, KRAS, NTS, SMAD4, POLE, BRAF, BMPR1A, POLD1, STK11, MUTYH and DNA mismatch repair genes), which are associated with severe hereditary syndromes, such as familial adenomatous polyposis and Lynch syndrome (also called hereditary non-polyposis colorectal cancer), account for only 5–7% of total CRC cases (4,12). Therefore, the remaining unknown heritability is probably explained by the interaction of common, low-penetrance variants identified through genome-wide association studies (GWAS). GWAS in ethnic/racial minority populations offers the opportunity to uncover genetic susceptibility factors and discover new genomic regions and loci that contribute risk for CRC development. Since 2007, more than 100 common risk variants have been successfully identified in GWAS, which have helped to elucidate the etiology of CRC (1318).

Lynch syndrome (LS) is an autosomal dominant hereditary cancer syndrome that accounts for approximately 3% of all colorectal cancer cases (19). From a clinical point of view, 10–82% (20) of LS cases are associated with a lifetime risk of developing CRC, unless the risk is significantly lower in other types of cancer (21,22). LS is caused by pathogenic germline mutations in a class of genes called DNA mismatch repair (MMR) genes, mainly MLH1, located in 3p22.2 chromosome, and MSH2, located in 2p21 chromosome (23), which represent 70–85% of cases of LS (24). Mutations found in MSH6 (2p16.3), PMS2 (7p22.1) (25) and MLH3 genes have lower incidence (26). Molecular investigations have also shown that MSH3 (27) and germline 3′deletions of the EPCAM gene, which lead to epigenetic silencing of MSH2 (28), are also implicated in the pathogenesis of LS. As a consequence of MMR pathway inactivation and loss of expression of MMR proteins, DNA replication errors accumulate typically resulting in microsatellite instability (MSI), which is generally detected in LS patients' tumor tissues (29). The diagnosis of LS involves three main steps, identification of patients and their familial history that meet the Amsterdam or Bethesda guidelines, presence of MSI in tumors and immunohistochemical analysis (IHC) of MMR protein expression. A definitive diagnosis of LS must be confirmed by detecting the germline mutations in MMR genes (30).

Non-invasive prenatal testing (NIPT) based on low-coverage massively parallel whole-genome sequencing of plasma DNA from pregnant women generates a large amount of data that provides the resources to investigate human genetic variations in the population. In our previous studies, we described the re-use of the data from NIPT for genome-scale population specific frequency determination of small DNA variants (31) and CNVs (32). Since pregnant women represent a relatively standard sample of the local female population, we assumed this NIPT data could also be used in the population study of CRC, the most common cancer in Slovakia. Some research concerning on genomic analysis of plasma from NIPT has also demonstrated NIPT data's efficiency and utility for viral genetic studies (33), genetic profiling of Vietnamese population (34) or detection of CNV aberrations (32,35).

The main aim of our study was a detailed analysis of common variants (MAF >0,05) that showed evidence of association with CRC in GWAS datasets and characterization of population variability from data generated by NIPT. We assumed that the genetic factors, mainly the increased specific population frequency of CRC and LS variants could be responsible for the high incidence of CRC in Slovakia. To test this hypothesis, allele frequencies of risk CRC variants identified in the Slovak population were compared with allele frequencies of risk CRC variants in 6 worldwide populations. As LS is among the most common hereditary CRC syndromes, the aim of our study was also to analyze population allele frequencies and describe clinical impacts of relevant variants located in known LS predisposing genes. To our knowledge, this was the first population study of CRC using NIPT data conducted exclusively in the Slovak population.

Materials and methods

Data source

The laboratory procedure used, to generate the NIPT data, were as follows: DNA from plasma of peripheral maternal blood was isolated for NIPT analysis from 1,501 pregnant women after obtaining a written informed consent consistent with the Helsinki declaration from the subjects. The population cohort consisted from women in reproductive age between 17–48 years with a median of 35 years. Genomic information from a sample consisted of maternal and fetal DNA fragments. Each included individual agreed to use their genomic data in an anonymized form for general biomedical research. The NIPT study (study ID 35900_2015) was approved by the Ethical Committee of the Bratislava Self-Governing Region (Sabinovska ul.16, 820 05 Bratislava) on 30th April of 2015 under the decision ID 03899_2015. Blood samples were collected to EDTA tubes and plasma was separated in dual centrifugation procedure. DNA was isolated from 700 µl of plasma using DNA Blood Mini kit (Qiagen) according to standard protocol. Sequencing libraries were prepared from each sample using TruSeq Nano kit HT (Illumina) following standard protocol with omission of DNA fragmentation step. Each sample was normalized to 4 nm library and the final concentration of libraries was 2,8 pM. Individual barcode labelled libraries were pooled and sequenced using low-coverage whole-genome sequencing on an Illumina NextSeq500 platform (Illumina) by performing paired end sequencing of 2×35 bases (36).

Data analysis

The datasets generated and/or analyzed during the current study are available in the DSpace repository, https://dspace.uniba.sk/xmlui/handle/123456789/27 (31).

Analyses of common variants previously reported to be risk variants for CRC

We combined genotype data from all previously reported GWAS studies available online (https://www.gwascentral.org/) for the years 2007–2020, specifically 66 GWAS studies of CRC risk variants that included individuals with European, Asian and African American ancestry. Using data from these GWAS datasets, we identified 116 risk variants associated with CRC, which were then merged with our data of identified variants from NIPT. Risk variants that were not found in NIPT data were excluded from the analysis. All identified variants in the Slovak population used for further analyses were common (MAFs >0.05). Subsequently, allele frequencies of CRC risk variants for each population (East Asian, South Asian, African, American, Finnish European and non-Finnish European) were extracted from the gnomAD database available online (v3.0, downloaded from http://gnomad.broadinstitute.org/downloads) and compared with our frequencies determined for the Slovak population from NIPT data. Allele frequency in each population and allele frequency differences were plotted using boxplots. Outliers of boxplots that represent variants with highly different frequencies between Slovak and non-Finnish populations were annotated via published literature and studies [in dbSNP (https://www.ncbi.nlm.nih.gov/snp/) and GWAS (https://www.gwascentral.org/)]. To assess the relations between allele frequency of CRC risk variants in each population, we also used Principal Component Analysis (PCA) using matplotlib.pyplot library, which reduces the dimension of the data to a graphically interpretable 2D or 3D dimension. Consequently, we obtained information on which populations have similar or different allele frequencies of the identified CRC risk variants.

Analyses of variants located in genes associated with LS

After analyzing variants associated with CRC, we focused on the study of variants associated with LS. First, we filtered out a group of variants located in 7 genes known to be associated with LS (MLH1, PMS2, MSH6, MLH3, MSH2, TGFBR2, EPCAM). The genomic locations of genes were determined by the GeneCards database (https://www.genecards.org/). From the dataset of identified variants in LS associated genes, we excluded variants in low complexity genomic regions (soft masked in the reference FASTA file), eliminating the variants that could represent sequencing artifacts or repetitive regions. All variants that were used for further analysis were annotated using Ensembl Variant Effect Predictor (VEP, version 101_GRCh38). In our dataset, based on ClinVar database annotation of the most common types of pathogenic and likely pathogenic variants associated with LS (https://www.ncbi.nlm.nih.gov/clinvar), we selected variants including frameshift, missense, nonsense, splice site, non-coding and UTR variants. After this filtering, allele frequencies for both groups of variants (all variants identified in LS genes and selected types of variants) for each population (East Asian, South Asian, African, American, Finnish European and non-Finnish European) were extracted from the gnomAD database (v3.0, downloaded from http://gnomad.broadinstitute.org/downloads) and compared with our frequencies determined for the Slovak population from NIPT data. Allele frequency in each population and allele frequency differences were plotted using boxplots and PCA analysis using matplotlib.pyplot library. Outliers of boxplots representing variants with allele frequency differences more than 10% were annotated via published literature and studies [in dbSNP (https://www.ncbi.nlm.nih.gov/snp/) and GWAS (https://www.gwascentral.org/].

Results

Analyses of common variants previously reported to be risk variants for CRC

In the analysis of the 66 GWAS studies that included all identified risk variants associated with colorectal carcinogenesis from 2007–2020, we investigated 116 variants. There were 25 independent CRC risk variants in Asian population, 62 risk variants found in European population, 27 risk variants in both European and Asian population and 2 risk variants located in African American population that were previously reported in GWAS (Table I).

Table I.

Identification of 116 risk variants associated with colorectal cancer from 66 genome-wide association studies between 2007 and 2020.

Table I.

Identification of 116 risk variants associated with colorectal cancer from 66 genome-wide association studies between 2007 and 2020.

First author, yearrs_IDChrPOSPOPGenePMID(Refs.)
Wang et al, 2017rs725250519q1333,084,158AFR AGPATCH128295283(46)
Wang et al, 2017rs5684893619q13.346,321,507AFR ASYMPK28295283(46)
Lu et al, 2019rs75426651p31.362,673,037ASNL1TD130529582(13)
Law et al, 2019rs121435411p32.355,247,852ASNTTC22   31089142(43)
Lu et al, 2019rs2013952361q44245,181,421ASNEFCAB230529582(13)
Lu et al, 2019rs76065622p16.348,686,695ASNPPP1R2130529582(13)
Lu et al, 2019rs1135695143q22.2133,748,789ASNSLCO2A130529582(13)
Lu et al, 2019rs126590175q23.2125,988,175ASNALDH7A1, PHAX30529582(13)
Law et al, 2019rs6399335q31.1134,467,751ASNC5orf66, LOC10537918831089142(43)
Jia et al, 2013rs6471615q31.1134,499,092ASNPITX123263487(58)
Law et al, 2019rs69337906p21.141,672,769ASNTFEB31089142(43)
Zeng et al, 2016rs47116896p21.141,692,812ASNTFEB26965516(14)
Schmit et al, 2018rs69063596p21.3135,528,378ASNFKBP529917119(18)
Lu et al, 2019rs38300416p21.3232,191,339ASNNOTCH430529582(13)
Lu et al, 2019rs658428310q24.2101,290,301ASNNKX2-330529582(13)
Lu et al, 2019rs7796913212p11.2131,594,813ASNDENND5B30529582(13)
Zeng et al, 2016rs1106443712p13.316,982,162ASNSPSB226965516(14)
Lu et al, 2019rs273098512q1243,130,624ASNPRICKLE130529582(13)
Lu et al, 2019rs188645013q22.173,986,628ASNKLF5, KLF1230529582(13)
Lu et al, 2019rs434175416q23.280,039,621ASNWWOX, MAF30529582(13)
Lu et al, 2019rs107864317p1210,707,241ASNPIRT30529582(13)
Law et al, 2019rs7397558817p13.3816,741ASNNXN31089142(43)
Law et al, 2019rs979788519q13.241,873,001ASNTMEM9131089142(43)
Law et al, 2019rs605528620p12.37,718,045ASNNone31089142(43)
Jia et al, 2008rs242327920p12.37,812,350ASNHAO123263487(58)
Law et al, 2019rs217959320q13.1242,660,286ASNTOX231089142(43)
Lu et al, 2019rs1383120q13.3257,475,191ASNGNAS30529582(13)
Law et al, 2019rs617767191p34.338,461,319EURNone31089142(43)
Peters et al, 2013rs109112511q25.3183,112,059EURLAMC123266556(59)
Whiffin et al, 2014rs109112511q25.3183,112,059EURLAMC124737748(60)
Houlston et al, 2010rs66877581q41222,164,948EURNone20972440(61)
Houlston et al, 2010rs66911701q41222,045,446EURDUSP1020972440(61)
Law et al, 2019rs116924352q11.298,275,354EURACTR1B31089142(43)
Law et al, 2019rs118930632q33.1199,601,925EURLOC10537383131089142(43)
Law et al, 2019rs75934222q33.1200,131,695EURNone31089142(43)
Orlando et al, 2016rs9921572q35219,154,781EURTMBIM127005424(62)
Law et al, 2019rs98318613p21.153,088,285EURNone31089142(43)
Law et al, 2019rs126359463q13.2112,916,918EURNone31089142(43)
Houlston et al, 2010rs109365993q26.2169,774,313EURMYNN20972440(61)
Schmit et al, 2018rs13708214q22.294,943,383EURNone29917119(18)
Law et al, 2019rs170352894q24106,048,291EURNone31089142(43)
Law et al, 2019rs756868614q31.21145,621,328EURHHIP31089142(43)
Schmit et al, 2014rs355092824q32.2163,333,405EURFSTL525023989(63)
Schmit et al, 2018rs587917125p13.140,281,797EURPTGER429917119(18)
Schmit et al, 2018rs27359405p15.331,296,486EURTERT29917119(18)
Peters et al, 2012rs28536685p15.331299910EURTERT21761138(64)
Schmit et al, 2018rs624049686p12.155,714,314EURBMP529917119(18)
Dunlop et al, 2012rs13213116p21.236,622,900EURCDKN1A22634755(48)
Law et al, 2019rs92717706p21.3232,594,248EURLOC10798744931089142(43)
Law et al, 2019rs31310436p21.3330,758,466EURHCG2031089142(43)
Law et al, 2019rs20706996p24.112,292,772EUREDN131089142(43)
Law et al, 2019rs69288646q21105,966,894EURNone31089142(43)
Law et al, 2019rs109518787p12.346,926,695EURNone31089142(43)
Law et al, 2019rs38010817p12.347,511,161EURTNS331089142(43)
Tomlinson et al, 2008rs168927668q23.3117,630,683EUREIF3H18372905(64)
Law et al, 2019rs14128349p21.322,110,131EURCDKN2B-AS131089142(43)
Schmit et al, 2018rs1099486010q11.2352,645,424EURA1CF29917119(18)
Al Tassan et al, 2015rs1090484910p1316,955,267EURCUBN25990418(15)
Law et al, 2019rs445016811p15.410,286,755EURSBF231089142(43)
Dunlop et al, 2012rs382499911q13.474,345,550EURPOLD322634755(48)
Tenesa et al, 2008rs380284211q23.1111,171,709EURCOLCA218372901(66)
Houlston et al, 2010rs713670212q13.1250,880,216EURLARP420972440(61)
Law et al, 2019rs739837512q13.357,540,848EURLRP1   31089142(43)
Schmit et al, 2018rs7201372612q24.21115,890,835EURMED13L29917119(18)
Schmit et al, 2018rs1016198013q13.234,093,518EURSTARD1329917119(18)
Law et al, 2019rs1242760013q13.337,460,648EURSMAD931089142(43)
Law et al, 2019rs4559703513q22.173,649,152EURKLF531089142(43)
Law et al, 2019rs133088913q22.378,609,615EURLINC0044631089142(43)
Law et al, 2019rs799393413q34111,074,915EURCOL4A231089142(43)
Tomlinson et al, 2011rs195763614q22.254,560,018EURBMP421655089(67)
Houlston et al, 2008rs444423514q22.254,410,919EURBMP419011631(68)
Tomlinson et al, 2011rs444423514q22.254,410,919EURBMP421655089(67)
Tomlinson et al, 2008rs1163271515q13.333,004,247EURNone18372905(65)
Tomlinson et al, 2008rs1696968115q13.332,993,111EURSCG518372905(65)
Tomlinson et al, 2011rs1696968115q13.332,993,111EURSCG521655089(67)
Tomlinson et al, 2008rs477958415q13.332,994,756EURCRAC118372905(65)
Law et al, 2019rs477631615q22.3167,007,813EURSMAD631089142(43)
Law et al, 2019rs1015251815q2368,177,162EURNone31089142(43)
Law et al, 2019rs749513215q26.191,172,901EURCRTC331089142(43)
Houlston et al, 2008rs992921816q22.168,820,946EURCDH119011631(68)
Law et al, 2019rs6133691816q23.280,007,266EURNone31089142(43)
Schmit et al, 2018rs269683916q24.186,340,448EURFOXF129917119(18)
Broderick et al, 2007rs493982718q2146,453,463EURSMAD717934461(69)
Tenesa et al, 2008rs493982718q2146,453,463EURSMAD718372901(66)
Law et al, 2019rs28524519p13.1116,420,817EURNone31089142(43)
Houlston et al, 2008rs1041121019q13.1133,532,300EURRHPN219011631(68)
Law et al, 2019rs1297927819q13.3349,218,602EURMAMSTR31089142(43)
Tomlinson et al, 2011rs481380220p12.36,699,595EURBMP221655089(67)
Peters et al, 2012rs481380220p12.36,699,595EURBMP221761138(64)
Houlston et al, 2008rs96125320p12.36,404,281EURBMP219011631(68)
Schmit et al, 2018rs229544420q11.2233,173,883EURPIGU29917119(18)
Schmit et al, 2018rs181050220q13.1349,057,488EURPTPN129917119(18)
Law et al, 2019rs378708920q13.3362,316,630EURRTEL131089142(42)
Houlston et al, 2010rs492538620q13.3360,921,044EURLAMA520972440(61)
Schumacher et al, 2015rs81248133p14.143,476,841EUR, ASNLRIG126151821(70)
Schumacher et al, 2015rs353603283p22.140,924,962EUR, ASNCTNNB126151821(70)
Lu et al, 2019rs14765706p22.129,809,860EUR, ASNHLA-G30529582(13)
Tomlinson et al, 2008rs24501158q23.3117,624,093EUR, ASNEIF3H18372905(67)
Tomlinson et al, 2008rs64696568q23.3117,647,788EUR, ASNEIF3H18372905(67)
Haiman et al, 2007rs69832678q24.21128,413,305EUR, ASNPOU5F1B17618282(71)
Tomlinson et al, 2007rs69832678q24.21128,413,305EUR, ASNPOU5F1B17618284(72)
Hutter et al, 2010rs69832678q24.21128,413,305EUR, ASNPOU5F1B21129217(73)
Cui et al, 2011rs69832678q24.21128,413,305EUR, ASNPOU5F1B21242260(4)
Law et al, 2019rs1225514110q25.2114,294,892EUR, ASNVTI1A   31089142(43)
Zhang et al, 2014rs70401710q22.380,819,132EUR, ASNZMIZ124836286(16)
Whiffin et al, 2014rs103520910q24.2101,345,366EUR, ASNSLC25A2824737748(60)
Zeng et al, 2016rs491968710q24.32104,595,248EUR, ASNCYP17A126965516(14)
Zhang et al, 2014rs1119617210q25.2114,726,843EUR, ASNTCF7L224836286(16)
Wang et al, 2014rs1224100810q25.2114,280,702EUR, ASNVTI1A25105248(75)
Zhang et al, 2014rs153511q12.261,597,972EUR, ASNFADS224836286(16)
Zhang et al, 2014rs17455011q12.261,571,478EUR, ASNFADS124836286(16)
Zhang et al, 2014rs424621511q12.261,564,299EUR, ASNFEN124836286(16)
Zhang et al, 2014rs17453711q12.261,552,680EUR, ASNMYRF24836286(16)
Law et al, 2019rs1084943812p13.316,412,036EUR, ASNNone   31089142(43)
Zhang et al, 2014rs1084943212p13.316,385,727EUR, ASNPLEKHG624836286(16)
Peters et al, 2013rs321781012p13.324,388,271EUR, ASNCCND223266556(59)
Whiffin et al, 2014rs321781012p13.324,388,271EUR, ASNCCND224737748(60)
Jia et al, 2013rs1077421412p13.324,368,352EUR, ASNCCND223263487(58)
Zhang et al, 2014rs1260352617p13.3800,593EUR, ASNNXN24836286(16)
Zhang et al, 2014rs722963918q21.146,450,976EUR, ASNSMAD724836286(16)
Zhang et al, 2014rs180046919q13.241,860,296EUR, ASNTMEM9124836286(16)
Zhang et al, 2014rs224171419q13.241,869,392EUR, ASNB9D224836286(16)
Schumacher et al, 2015rs60668252020q13.13897,353EUR, ASNPREX126498495(70)
Schumacher et al, 2015rs606682520q13.1347,340,117EUR, ASNPREX126151821(70)
Dunlop et al, 2012rs5934683Xp22.29,751,474EUR, ASNSHROOM222634755(48)

[i] Chr, chromosome; POS, position; POP, population; AFR A, African American; ASN, Asian; EUR, European.

After merging all identified variants from GWAS (116 risk variants) with our NIPT data, we identified 106 common risk CRC variants (Table SI), while 10 risk variants that were not called in the Slovak population were excluded from further analysis. The allele frequencies of 106 variants identified in our population sample (Slovak population) and the allele frequencies of variants for 6 world populations (East Asian, South Asian, African, American, Finnish European and non-Finnish European) obtained by gnomAD database (Table SII) are shown in graphical comparison by Boxplots (Fig. 1) and principal component analysis (PCA) (Fig. 2). As shown in Fig. 1, the MAF ranged from 0.0–0.963109 in 6 world populations that is comparable with the Slovak population (0.0521–0.931). The median allele frequency for the Slovak population reached the value of 0.4072, which is closest to the value of the median of the American population (MED=0.3991) and Finnish population (MED=0.4166). PCA placed our sample set most closely to the two gnomAD population sample sets, i.e., to the Finnish and non-Finnish European population.

Next, we compared known allele frequencies of 106 CRC risk variants in our sample set from the Slovak population to allele frequencies of CRC variants in six world populations. The final findings of allele frequency differences are shown in Fig. 3. The median allele frequency for comparing the Slovak population and non-Finnish European population reached the value of 0.002285. Together, we identified 14 outliers in Fig. 3 (3 of 14 variants reached a similar value and overlapped, so they are not clearly visible in Fig. 3). Since the same variant rs4246215 was identified in 4 different population comparisons (Slovak-American, Slovak-Finnish European, Slovak-non-Finnish European, and also identified in Slovak-East Asian population comparison) and the rs3131043 variant in 2 population comparisons (Slovak-Finnish European and Slovak-non-Finnish European population comparison), we identified a total of 10 variants, whose difference in population allele frequency was more than 10%. Table II includes annotation information about these variants by dbSNP NCBI, ClinVar database, and population comparison in which they were identified.

Table II.

Outliers identified in boxplots that show allele frequency differences of Slovak and the other six world populations for 106 risk colorectal cancer variants identified from genome-wide association studies.

Table II.

Outliers identified in boxplots that show allele frequency differences of Slovak and the other six world populations for 106 risk colorectal cancer variants identified from genome-wide association studies.

rs_IDPopulation comparisonChrPOSVariant typeGeneConsequenceClinical significance
rs5934683Slovak-East AsianchrX9783434SNVGPR143Intron VariantNot Reported in ClinVar
rs7252505Slovak-Africanchr1933084158SNVGPATCH1Intron VariantNot Reported in ClinVar
rs4779584Slovak-East Asianchr1532702555SNVNoneNoneNot Reported in ClinVar
rs174550 Slovak-Americanchr1161804006SNVFADS1Intron VariantNot Reported in ClinVar
rs4246215 Slovak-Americanchr1161796827SNVFEN13′ UTR VariantNot Reported in ClinVar
Slovak-European (Finnish)chr1161796827SNVFEN13′ UTR VariantNot Reported in ClinVar
Slovak-European (non-Finnish)chr1161796827SNVFEN13′ UTR VariantNot Reported in ClinVar
Slovak-East Asianchr1161796827SNVFEN13′ UTR VariantNot Reported in ClinVar
rs10904849Slovak-European (non-Finnish)chr1016955267SNVCUBNIntron VariantNot Reported in ClinVar
rs6928864Slovak-Africanchr6105519019SNVNoneNoneNot Reported in ClinVar
rs3131043Slovak-European (non-Finnish)chr630790689SNVHCG20Intron VariantNot Reported in ClinVar
Slovak-European (Finnish)chr630790689SNVHCG20Intron VariantNot Reported in ClinVar
rs12659017Slovak-East Asianchr5126652483SNVNoneNoneNot Reported in ClinVar
rs397775554Slovak-European (non-Finnish)chr5 40281696-40281704IndelNoneNoneNot Reported in ClinVar

[i] Chr, chromosome; POS, position; SNV, single nucleotide variant.

Analyses of variants located in genes associated with LS

In the analysis of LS, we identified 1212 variants in our sample set from NIPT that were located in genes known to be associated with LS, i.e., MLH1, PMS2, MSH6, TGFBR2, MLH3, MSH2 and EPCAM. After excluding variants from low complexity regions, we obtained 648 variants that were finally annotated by VEP and used for further analysis. The allele frequencies of 648 variants identified in our population sample (Slovak population) and the allele frequencies of variants for 6 world populations (East Asian, South Asian, African, American, Finnish European and non-Finnish European) obtained by gnomAD database (Table SIII) are shown in graphical comparison by Boxplots (Fig. 4) and principal component analysis (PCA) (Fig. 5). As shown in Fig. 4, the MAF ranged from 0.0–1.0 in 6 world populations. In the Slovak population, all variants were with MAF>0.05 (0.0502–1.0). The median allele frequency for the Slovak population reached the value of 0.2204, which is closest to the value of the median of the South Asian population (MED=0.221274). PCA placed our sample set most closely to the non-Finnish European population (Fig. 5).

In the next step, to identify variants having significantly different frequencies, we compared known allele frequencies of 648 variants located in genes associated with LS identified in our sample set from the Slovak population to allele frequencies of these variants in six gnomAD world populations. The final findings of allele frequency differences are shown in Fig. 6. The median allele frequency for the comparison of the Slovak population and non-Finnish European population reached the value of −0.01093. By comparing the allele frequency of variants of the Slovak and non-Finnish populations, we identified a total of 64 outliers. Most outliers were found in the MSH2 gene, others in MSH6, TGFBR2, PMS2, MLH1 and EPCAM. We did not identify any outlying variant in the MLH3 gene. The variation type of all outliers was ‘intronic variant’ and the clinical significance of all outliers was not reported in ClinVar. All this annotation information, also including chromosome position, variants ID, reference and alternative allele, genes, is available in Table SIV.

Our analysis also included allele frequency comparison of selected variants from 648 variants identified in our sample set of NIPT data. We focused on frameshift, missense, nonsense, splice site, non-coding and UTR variants, annotated in the ClinVar database as the most common types of pathogenic and likely pathogenic variants associated with LS. However, from these selected types of variants, we found only UTR and non-coding variants in our dataset of 648 variants. Other types of variants (downstream, upstream, and intron) were excluded from further analysis. Finally, we selected 18 variants, 10 UTR and 8 non-coding variants (all selected variants with annotation information by VEP and ClinVar are available in Table III). We compared known allele frequencies of these 18 selected variants identified in our population sample (Slovak population) to the six gnomAD world populations. The final findings of allele frequency differences are shown in Fig. 7. The median allele frequency for the comparison of the Slovak population and non-Finnish European population reached the value of 0.014215, which is closest to the value of the median allele frequency comparison of the South Asian and Slovak population (MED=0.014186). By comparing the allele frequency of variants of the Slovak and six gnomAD world population, we identified a total of 4 outliers-rs10951973, rs10951972 (identified in Slovak-American population comparison), rs6791557 (in Slovak-American and Slovak-non-Finnish European population comparison) and rs9852378 (in Slovak-South Asian population comparison). All outliers were non-coding variants, rs10951973 and rs10951972 located in the PMS2 and rs6791557 located in the TGFBR2 were not reported in ClinVar. The rs9852378 SNP, detected in the MLH1, was reported as benign by ClinVar.

Table III.

Identification of 18 selected variants (UTR and non-coding) from all 648 variants in genes associated with Lynch syndrome from non-invasive prenatal testing data in the Slovak population.

Table III.

Identification of 18 selected variants (UTR and non-coding) from all 648 variants in genes associated with Lynch syndrome from non-invasive prenatal testing data in the Slovak population.

rs_IDChrPOSVariant typeGeneConsequenceClinical significance
rs11901645chr247510079SNVMSH23′ UTR variantNot reported in ClinVar
rs11891189chr247510259SNVMSH23′ UTR variantNot reported in ClinVar
rs876936chr247513059SNVMSH23′ UTR variantNot reported in ClinVar
rs72872839chr247565070SNVMSH23′ UTR variantNot reported in ClinVar
rs17036769chr247633503SNVMSH23′ UTR variantNot reported in ClinVar
rs2969774chr247661684SNVMSH23′ UTR variantNot reported in ClinVar
rs2952372chr247661919SNVMSH23′ UTR variantNot reported in ClinVar
rs2969773chr247662141SNVMSH23′ UTR variantNot reported in ClinVar
rs2705765chr247662389SNVMSH23′ UTR variantNot reported in ClinVar
rs12328344chr247662542SNVMSH23′ UTR variantNot reported in ClinVar
rs10427209chr247709297SNVMSH6Non-coding transcript exon variantNot reported in ClinVar
rs10427344chr247709476SNVMSH6Non-coding transcript exon variantNot reported in ClinVar
rs3136240chr247784947SNVMSH6Non-coding transcript exon variantNot reported in ClinVar
rs6791557chr330614676SNVTGFBR2Non-coding transcript exon variantNot reported in ClinVar
rs1817338chr330631239SNVTGFBR2Non-coding transcript exon variantNot reported in ClinVar
rs9852378chr336997280SNVMLH1Non-coding transcript exon variantBenign in ClinVar
rs10951972chr76002187SNVPMS2Non-coding transcript exon variantNot reported in ClinVar
rs10951973chr76002205SNVPMS2Non-coding transcript exon variantNot reported in ClinVar

[i] Chr, chromosome; POS, position; SNV, single nucleotide variant.

Finally, we analyzed allele frequencies of pathogenic and likely pathogenic variants associated with Lynch syndrome annotated in the ClinVar database. From 229 SNPs with pathogenic and likely pathogenic clinical significance, only 15 have non-zero AF records in the gnomAD database. As shown in Fig. 8, all found AF are significantly below 5% (Table IV and Fig. 8).

Table IV.

Identification of 15 pathogenic and likely pathogenic variants associated with Lynch syndrome with non-zero allele frequency in gnomAD database.

Table IV.

Identification of 15 pathogenic and likely pathogenic variants associated with Lynch syndrome with non-zero allele frequency in gnomAD database.

Allele frequency in population

rs_IDChrPOSREF AlleleALT AlleleAfricanAmericanEast AsianEuropean (Finnish)European (non-Finnish)South Asian
rs63750615247403333GT00000 3.28×10−5
rs1194793421247414417AGA0000 2.75×10−50
rs63750636247476492CT0000 1.55×10−50
rs63749873247795903CG0000 3.10×10−50
rs587783056247799684GTTG 4.76×10−500000
rs63751017247800714CT0000 3.10×10−50
rs876660943247806359GT0000 1.55×10−50
rs63751221337001045CT 2.38×10−500000
rs58777933875977589GA 4.86×10−5000 1.58×10−50
rs26760816175982885CT 4.80×10−500 9.61×10−5 1.55×10−50
rs6375142275986838GA 2.38×10−500000
rs6375025075986933AAT0000 9.30×10−50
rs20064058575992018GA0000 1.55×10−50
rs26760815475995572ACTGTA0000 1.55×10−50
rs6375087176002590GA0 7.33×10−50000

[i] Chr, chromosome; POS, position; REF, reference; ALT, alternative.

Discussion

Population genetic studies currently have a huge impact on the study of genomics (37). The detection of risk variants in a population and identifying their genetic relationships have advanced our understanding of the human genome's variability and led to the elucidation of many factors that influence cancer risk. In recent years, NGS technologies have played a key role in colorectal cancer research and have become a useful tool for cancer diagnostics and screening (3841). Due to the high incidence of colorectal cancer in the Slovak population, it is crucial to determine the possible causes of the high incidence of this disease in Slovakia.

Non-invasive prenatal testing of common fetal chromosomal aberrations, using low-coverage massively parallel whole-genome sequencing of maternal plasma cfDNA of pregnant women, has become the fastest low-cost genomic DNA test that is rapidly implemented in clinical practice. Currently, more than 3 million NIPT tests are carried out worldwide each year, and the large amount of data generated during NIPT provides the resources to investigate human genetic variations in the population (31). In our study, we analyzed low-coverage massively parallel whole-genome sequencing data of total plasma DNA from pregnant women generated for NIPT screening to characterize the variants in genes associated with CRC and LS in the Slovak population. To our knowledge, the present study is the first population analysis of CRC and LS variants worldwide and also in the Slovak population using NIPT data. We illustrate the utility of these genomic data for clinical genetics and population studies.

Over the past two decades, GWAS offer the opportunity to uncover genetic susceptibility factors for CRC and provide insights into the biological basis of CRC etiology. These studies have demonstrated that only a fraction of CRC heritability is explained by known risk-conferring genetic variation, whereas the remaining genetic risk of CRC may be accounted for by a combination of high-prevalence and low-penetrance of common genetic variants. To date, a large number of common genetic variants have been identified by the GWAS approach, which has intimately connected to the onset of CRC (13,18,4246).

By pooling GWAS data of risk variants associated with colorectal carcinogenesis from 2007–2020 and data variants in our population sample from NIPT, we have identified 106 common risk CRC variants. When we compared allele frequencies of these variants to allele frequencies in six gnomAD world population, finally 13 common risk variants were found that showed statistically significant differences in population allele frequencies-rs5934683, rs7252505, rs4779584, rs1535, rs174550, rs4246215, rs11196172, rs10904849, rs6928864, rs3131043, rs1476570, rs12659017, rs397775554.

The SNP rs5934683 is located on chromosome Xp22.2 between two genes, GPR143 (G protein-coupled receptor 143), which is expressed by melanocytes and retinal pigment epithelium and SHROOM2 (shroom family member 2), a human homolog of the Xenopus laevis APX gene that has important functions in cell morphogenesis including endothelial and epithelial tissue development (44). Missense mutations in this gene have been detected in large-scale screens for recurring mutations in cancer cell lines. Both GPR143 and SHROOM2 play a role in melanosome biogenesis and retinal pigmentation. It is known that abnormal retinal pigmentation, similar to the congenital hypertrophy of retinal pigment epithelium lesions, are typical of the familial adenomatous polyposis syndrome (FAP), one of the inherited syndromes of CRC (47). The relationship between Xp22.2 and CRC risk represents the first evidence for the role of X-chromosome variation in predisposition to non-sex-specific cancer (48).

The SNP rs7252505, located in the 19q13 risk locus, is in an intron of the gene GPATCH1 (G-patch domain containing 1). Although GPATCH1 is expressed in the colon, little is known about its function other than the fact that it contains a G-patch domain, a domain typically associated with RNA processing. One study found that rs7252505 was associated with CRC in African Americans (46,49).

Intergenic variant rs4779584 in chromosomal region 15q13.3 lies between SCG5 and GREM1, and the association between this SNP to CRC has been identified in several GWAS studies (13,50).

The rs4246215 polymorphism is located in the FEN1 in the long arm of chromosome 11 (11q12.2). The association between this SNP and the potential risk of different types of cancers, including esophageal, lung, gastrointestinal, gallbladder, breast cancer in Chinese and Iran populations, glioma and childhood leukemia, has been previously studied. The rs4246215 variant was also associated with colorectal cancer in East Asians and the Chinese population (51,52).

To identify variants that may predispose to LS and may cause the high incidence of CRC in Slovakia, we used NIPT data, including variants with at least 5% AF and coverage at least 100 reads per variant. To verify the reliability of the found variants using NIPT, we selected 15 variants with AF below 5% and validated them using Sanger sequencing. For this reason, it is not possible to find rare variants with AF under 5%. Initially, we selected gene variants known to be associated with LS and we focused on their population AF in gnomAD database and as well as on pathogenicity as reported in public database ClinVar. No publications are available for all variants showing statistically significant differences in population allele frequencies and selected 18 variants. The rs9852378 SNP was reported as benign by ClinVar, and other variants were not reported in ClinVar.

Our study has several key shortcomings. None of the variants identified in this study are pathogenic or likely pathogenic due to their extremely low frequency in the general population (Fig. 8). From the total number of pathogenic or likely pathogenic variants annotated in the ClinVar database, we could determine the population frequency of only 15 variants even when using the gnomAD database (Table IV). Second, the sample size was relatively small and it is strongly biased towards females. We assume that even larger sample sets will also not offer opportunities to detect such low frequencies of LS variants in the population using NIPT data. Third, a substantial portion of identified variants was removed from analyses due to technical limitations, mainly because of their location in low complexity regions. Although these could be technical artifacts (53), they could also be real variants having biological effects that are yet generally hardly determinable, but likely existing (54). Moreover, colorectal cancer is a disease caused by a combination of multiple genes and environmental factors. To assess the relationship between the variants identified in population and CRC development, it is very important in future research to study the interaction between genes and also the environment on the colorectal cancer risk. Although suitable for the determination of general population frequencies of independent variants, NIPT data are unsuitable for calculations (such as polygenic risk score determinations) based on exact combinations of these variants in individuals, which may de facto determine the risk of individuals to develop certain diseases.

The underlying mechanism for a high incidence of CRC in the Slovak population is still unclear at the moment; however, it is possible that genetic factors, like the most common inherited syndrome-LS, play a crucial role in colorectal etiology. We have performed a literature search in PubMed focused on population studies of CRC and LS in Slovakia from 2010–2020 using next-generation sequencing. In the Slovak population, only a few population studies of risk variants have been conducted to elucidate the etiology of CRC (5557). In general, little is known about risk variants associated with CRC or LS in the Slovak population.

Identifying mutations associated with CRC in populations with high mortality rate, such as the Slovak population, is important to reduce the incidence of this multifactor disorder. The findings from these studies suggest a lack of understanding of the mechanism of many risk variants of CRC. Due to study limitations, we could not identify any pathogenic variants associated with LS in the Slovak population using NIPT data. On the other hand, NIPT data is not a major obstacle to better results, as pathogenic variants have extremely low frequencies in the general population. Even in most cases, the frequencies are not known. However, we identified several promising common risk variants associated with CRC previously reported in GWAS studies that represent variants with highly different frequencies between Slovak and non-Finnish populations in boxplots. Since NIPT expands rapidly to millions of individuals each year, the reuse of these data reduces the cost of large-scale population studies and likely provides an acceptable background for information about genomic variation. Finally, future population studies on larger sample sets with various types of mutations are needed to reveal new mechanisms of pathogenicity and links to new biological pathways, which may be useful in designing preventive strategies and treatment of CRC.

Supplementary Material

Supporting Data

Acknowledgements

Not applicable.

Funding

The present study was supported by the PANGAIA project H2020-MSCA-RISE-2019 (grant no. 872539) funded under H2020-EU.1.3.3. Programme, the OP Integrated Infrastructure for the project ‘Long term strategic research and development focused on the occurrence of Lynch syndrome in the Slovak population and possibilities of prevention of tumors associated with this syndrome’ (grant no. 313011V578) co-financed by the European Regional Development Fund (ERDF), the Operational Program Integrated Infrastructure (grant no. 313011F988) co-financed by the ERDF, and the Scientific Grant Agency of the Ministry of Education, Science, Research and Sport of the Slovak Republic and the Slovak Academy of Sciences (grant no. 1/0305/19).

Availability of data and materials

The datasets generated and/or analyzed during the current study are available in the DSpace repository, https://dspace.uniba.sk/xmlui/handle/123456789/27.

Authors' contributions

NF and JG performed data analysis. NF was responsible for the literature search and manuscript writing. JG, JB and JR were responsible for designing the study and supervising the work. TS conceived the idea of the project and TS, JB and JR performed proofreading of the manuscript. JG and JB confirm the authenticity of all the raw data. All authors have read and approved the final manuscript.

Ethics approval and consent to participate

The non-invasive prenatal testing study (study ID. 35900_2015) was approved by the Ethical Committee of the Bratislava Self-Governing Region (Sabinovska ul.16, 820 05 Bratislava) on 30th April 2015 (approval no. 03899_2015).

Patient consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

References

1 

Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA and Jemal A: Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 68:394–424. 2018. View Article : Google Scholar : PubMed/NCBI

2 

Global Cancer Observatory, . International Agency for Research on Cancer. (Lyon, France). 2020.https://gco.iarc.fr/today/data/factsheets/populations/703-slovakia-fact-sheets.pdfNovember 9–2020

3 

Thanikachalam K and Khan G: Colorectal cancer and nutrition. Nutrients. 11:1642019. View Article : Google Scholar : PubMed/NCBI

4 

Rawla P, Sunkara T and Barsouk A: Epidemiology of colorectal cancer: Incidence, mortality, survival, and risk factors. Prz Gastroenterol. 14:89–103. 2019.PubMed/NCBI

5 

Cai S, Li Y, Ding Y, Chen K and Jin M: Alcohol drinking and the risk of colorectal cancer death: A meta-analysis. Eur J Cancer Prev. 23:532–539. 2014. View Article : Google Scholar : PubMed/NCBI

6 

Dashti SG, Buchanan DD, Jayasekara H, Ouakrim DA, Clendenning M, Rosty C, Winship IM, Macrae FA, Giles GG, Parry S, et al: Alcohol consumption and the risk of colorectal cancer for mismatch repair gene mutation carriers. Cancer Epidemiol Biomarkers Prev. 26:366–375. 2017. View Article : Google Scholar : PubMed/NCBI

7 

Botteri E, Iodice S, Bagnardi V, Raimondi S, Lowenfels AB and Maisonneuve P: Smoking and colorectal cancer: A meta-analysis. JAMA. 300:2765–2778. 2008. View Article : Google Scholar : PubMed/NCBI

8 

Limsui D, Vierkant RA, Tillmans LS, Wang AH, Weisenberger DJ, Laird PW, Lynch CF, Anderson KE, French AJ, Haile RW, et al: Cigarette smoking and colorectal cancer risk by molecularly defined subtypes. J Natl Cancer Inst. 102:1012–1022. 2010. View Article : Google Scholar : PubMed/NCBI

9 

Ordóñez-Mena JM, Walter V, Schöttker B, Jenab M, O'Doherty MG, Kee F, Bueno-de-Mesquita B, Peeters PH, Stricker BH, Ruiter R, et al: Impact of prediagnostic smoking and smoking cessation on colorectal cancer prognosis: A meta-analysis of individual patient data from cohorts within the CHANCES consortium. Ann Oncol. 29:472–483. 2018. View Article : Google Scholar : PubMed/NCBI

10 

Thrift AP, Gong J, Peters U, Chang-Claude J, Rudolph A, Slattery ML, Chan AT, Locke AE, Kahali B, Justice AE, et al: Mendelian randomization study of body mass index and colorectal cancer risk. Cancer Epidemiol Biomarkers Prev. 24:1024–1031. 2015. View Article : Google Scholar : PubMed/NCBI

11 

Gharahkhani P, Ong JS, An J, Law MH, Whiteman DC, Neale RE and MacGregor S: Effect of increased body mass index on risk of diagnosis or death from cancer. Br J Cancer. 120:565–570. 2019. View Article : Google Scholar : PubMed/NCBI

12 

Dekker E, Tanis PJ, Vleugels JLA, Kasi PM and Wallace MB: Colorectal cancer. Lancet. 394:1467–1480. 2019. View Article : Google Scholar : PubMed/NCBI

13 

Lu Y, Kweon SS, Tanikawa C, Jia WH, Xiang YB, Cai Q, Zeng C, Schmit SL, Shin A, Matsuo K, et al: Large-scale genome-wide association study of east asians identifies loci associated with risk for colorectal cancer. Gastroenterology. 156:1455–1466. 2019. View Article : Google Scholar : PubMed/NCBI

14 

Zeng C, Matsuda K, Jia WH, Chang J, Kweon SS, Xiang YB, Shin A, Jee SH, Kim DH, Zhang B, et al: Identification of susceptibility loci and genes for colorectal cancer risk. Gastroenterology. 150:1633–1645. 2016. View Article : Google Scholar : PubMed/NCBI

15 

Al-Tassan NA, Whiffin N, Hosking FJ, Palles C, Farrington SM, Dobbins SE, Harris R, Gorman M, Tenesa A, Meyer BF, et al: A new GWAS and meta-analysis with 1000Genomes imputation identifies novel risk variants for colorectal cancer. Sci Rep. 5:104422015. View Article : Google Scholar : PubMed/NCBI

16 

Zhang B, Jia WH, Matsuda K, Kweon SS, Matsuo K, Xiang YB, Shin A, Jee SH, Kim DH, Cai Q, et al: Large-scale genetic study in East Asians identifies six new loci associated with colorectal cancer risk. Nat Genet. 46:533–542. 2014. View Article : Google Scholar : PubMed/NCBI

17 

Takahashi Y, Sugimachi K, Yamamoto K, Niida A, Shimamura T, Sato T, Watanabe M, Tanaka J, Kudo S, Sugihara K, et al: Japanese genome-wide association study identifies a significant colorectal cancer susceptibility locus at chromosome 10p14. Cancer Sci. 108:2239–2247. 2017. View Article : Google Scholar : PubMed/NCBI

18 

Schmit SL, Edlund CK, Schumacher FR, Gong J, Harrison TA, Huyghe JR, Qu C, Melas M, Van Den Berg DJ, Wang H, et al: Novel common genetic susceptibility loci for colorectal cancer. J Natl Cancer Inst. 111:146–157. 2019. View Article : Google Scholar : PubMed/NCBI

19 

Biller LH, Syngal S and Yurgelun MB: Recent advances in lynch syndrome. Fam Cancer. 18:211–219. 2019. View Article : Google Scholar : PubMed/NCBI

20 

Yurgelun MB and Hampel H: Recent advances in lynch syndrome: Diagnosis, treatment, and cancer prevention. Am Soc Clin Oncol Educ Book. 38:101–109. 2018. View Article : Google Scholar : PubMed/NCBI

21 

Møller P, Seppälä T, Bernstein I, Holinski-Feder E, Sala P, Evans DG, Lindblom A, Macrae F, Blanco I, Sijmons R, et al: Incidence of and survival after subsequent cancers in carriers of pathogenic MMR variants with previous cancer: A report from the prospective lynch syndrome database. Gut. 66:1657–1664. 2017. View Article : Google Scholar : PubMed/NCBI

22 

Møller P, Seppälä TT, Bernstein I, Holinski-Feder E, Sala P, Evans DG, Lindblom A, Macrae F, Blanco I, Sijmons RH, et al: Cancer risk and survival in carriers by gene and gender up to 75 years of age: A report from the prospective lynch syndrome database. Gut. 67:1306–1316. 2018. View Article : Google Scholar : PubMed/NCBI

23 

Soares BL, Brant AC, Gomes R, Pastor T, Schneider NB, Ribeiro-Dos-Santos Â, de Assumpção PP, Achatz MI, Ashton-Prolla P and Moreira MA: Screening for germline mutations in mismatch repair genes in patients with lynch syndrome by next generation sequencing. Fam Cancer. 17:387–394. 2018. View Article : Google Scholar : PubMed/NCBI

24 

Cox VL, Bamashmos AA, Foo WC, Gupta S, Yedururi S, Garg N and Kang HC: Lynch syndrome: Genomics update and imaging review. Radiographics. 38:483–499. 2018. View Article : Google Scholar : PubMed/NCBI

25 

Le S, Ansari U, Mumtaz A, Malik K, Patel P, Doyle A and Khachemoune A: Lynch syndrome and muir-torre syndrome: An update and review on the genetics, epidemiology, and management of two related disorders. Dermatol Online J. 23:130302017.

26 

Peltomäki P: Update on lynch syndrome genomics. Fam Cancer. 15:385–393. 2016. View Article : Google Scholar : PubMed/NCBI

27 

Duraturo F, Liccardo R, Cavallo A, De Rosa M, Grosso M and Izzo P: Association of low-risk MSH3 and MSH2 variant alleles with Lynch syndrome: Probability of synergistic effects. Int J Cancer. 129:1643–1650. 2011. View Article : Google Scholar : PubMed/NCBI

28 

Kuiper RP, Vissers LELM, Venkatachalam R, Bodmer D, Hoenselaar E, Goossens M, Haufe A, Kamping E, Niessen RC, Hogervorst FB, et al: Recurrence and variability of germline EPCAM deletions in Lynch syndrome. Hum Mutat. 32:407–414. 2011. View Article : Google Scholar : PubMed/NCBI

29 

Shah SN, Hile SE and Eckert KA: Defective mismatch repair, microsatellite mutation bias, and variability in clinical cancer phenotypes. Cancer Res. 70:431–435. 2010. View Article : Google Scholar : PubMed/NCBI

30 

Martin-Morales L, Rofes P, Diaz-Rubio E, Llovet P, Lorca V, Bando I, Perez-Segura P, de la Hoya M, Garre P, Garcia-Barberan V and Caldes T: Novel genetic mutations detected by multigene panel are associated with hereditary colorectal cancer predisposition. PLoS One. 13:e02038852018. View Article : Google Scholar : PubMed/NCBI

31 

Budis J, Gazdarica J, Radvanszky J, Harsanyova M, Gazdaricova I, Strieskova L, Frno R, Duris F, Minarik G, Sekelska M, et al: Non-invasive prenatal testing as a valuable source of population specific allelic frequencies. J Biotechnol. 299:72–78. 2019. View Article : Google Scholar : PubMed/NCBI

32 

Pös O, Budis J, Kubiritova Z, Kucharik M, Duris F, Radvanszky J and Szemes T: Identification of structural variation from NGS-Based non-invasive prenatal testing. Int J Mol Sci. 20:44032019. View Article : Google Scholar : PubMed/NCBI

33 

Liu S, Huang S, Chen F, Zhao L, Yuan Y, Francis SS, Fang L, Li Z, Lin L, Liu R, et al: Genomic analyses from non-invasive prenatal testing reveal genetic associations, patterns of viral infections, and Chinese population history. Cell. 175:347–359. 2018. View Article : Google Scholar : PubMed/NCBI

34 

Tran NH, Vo TB, Nguyen VT, Tran NT, Trinh THN, Pham HAT, Dao THT, Nguyen NM, Van YLT, Tran VU, et al: Genetic profiling of Vietnamese population from large-scale genomic analysis of non-invasive prenatal testing data. Sci Rep. 10:191422020. View Article : Google Scholar : PubMed/NCBI

35 

Pös O, Budiš J and Szemes T: Recent trends in prenatal genetic screening and testing. F1000Res. 8:F10002019. View Article : Google Scholar : PubMed/NCBI

36 

Minarik G, Repiska G, Hyblova M, Nagyova E, Soltys K, Budis J, Duris F, Sysak R, Bujalkova MG, Vlkova-Izrael B, et al: Utilization of benchtop next generation sequencing platforms ion torrent PGM and MiSeq in noninvasive prenatal testing for chromosome 21 trisomy and testing of impact of in silico and physical size selection on its analytical performance. PLoS One. 10:e01448112015. View Article : Google Scholar : PubMed/NCBI

37 

Beyene J and Pare G: Statistical genetics with application to population-based study design: A primer for clinicians. Eur Heart J. 35:495–500. 2014. View Article : Google Scholar : PubMed/NCBI

38 

Zhu L, Huang Y, Fang X, Liu C, Deng W, Zhong C, Xu J, Xu D and Yuan Y: A novel and reliable method to detect microsatellite instability in colorectal cancer by next-generation sequencing. J Mol Diagn. 20:225–231. 2018. View Article : Google Scholar : PubMed/NCBI

39 

Yurgelun MB, Allen B, Kaldate RR, Bowles KR, Judkins T, Kaushik P, Roa BB, Wenstrup RJ, Hartman AR and Syngal S: Identification of a variety of mutations in cancer predisposition genes in patients with suspected lynch syndrome. Gastroenterology. 149:604–613. 2015. View Article : Google Scholar : PubMed/NCBI

40 

Valle L, de Voer RM, Goldberg Y, Sjursen W, Försti A, Ruiz-Ponte C, Caldés T, Garré P, Olsen MF, Nordling M, et al: Update on genetic predisposition to colorectal cancer and polyposis. Mol Aspects Med. 69:10–26. 2019. View Article : Google Scholar : PubMed/NCBI

41 

Budiš J, Kucharík M, Ďuriš F, Gazdarica J, Zrubcová M, Ficek A, Szemes T, Brejová B and Radvanszky J: Dante: Genotyping of known complex and expanded short tandem repeats. Bioinformatics. 35:1310–1317. 2019. View Article : Google Scholar : PubMed/NCBI

42 

Jiao S, Peters U, Berndt S, Brenner H, Butterbach K, Caan BJ, Carlson CS, Chan AT, Chang-Claude J, Chanock S, et al: Estimating the heritability of colorectal cancer. Hum Mol Genet. 23:3898–3905. 2014. View Article : Google Scholar : PubMed/NCBI

43 

Law PJ, Timofeeva M, Fernandez-Rozadilla C, Broderick P, Studd J, Fernandez-Tajes J, Farrington S, Svinti V, Palles C, Orlando G, et al: Association analyses identify 31 new risk loci for colorectal cancer susceptibility. Nat Commun. 10:21542019. View Article : Google Scholar : PubMed/NCBI

44 

Zhang K, Civan J, Mukherjee S, Patel F and Yang H: Genetic variations in colorectal cancer risk and clinical outcome. World J Gastroenterol. 20:4167–4177. 2014. View Article : Google Scholar : PubMed/NCBI

45 

Hofer P, Hagmann M, Brezina S, Dolejsi E, Mach K, Leeb G, Baierl A, Buch S, Sutterlüty-Fall H, Karner-Hanusch J, et al: Bayesian and frequentist analysis of an Austrian genome-wide association study of colorectal cancer and advanced adenomas. Oncotarget. 8:98623–98634. 2017. View Article : Google Scholar : PubMed/NCBI

46 

Wang H, Schmit SL, Haiman CA, Keku TO, Kato I, Palmer JR, van den Berg D, Wilkins LR, Burnett T, Conti DV, et al: Novel colon cancer susceptibility variants identified from a genome-wide association study in African Americans. Int J Cancer. 140:2728–2733. 2017. View Article : Google Scholar : PubMed/NCBI

47 

Closa A, Cordero D, Sanz-Pamplona R, Solé X, Crous-Bou M, Paré-Brunet L, Berenguer A, Guino E, Lopez-Doriga A, Guardiola J, et al: Identification of candidate susceptibility genes for colorectal cancer through eQTL analysis. Carcinogenesis. 35:2039–2046. 2014. View Article : Google Scholar : PubMed/NCBI

48 

Dunlop MG, Dobbins SE, Farrington SM, Jones AM, Palles C, Whiffin N, Tenesa A, Spain S, Broderick P, Ooi LY, et al: Common variation near CDKN1A, POLD3 and SHROOM2 influences colorectal cancer risk. Nat Genet. 44:770–776. 2012. View Article : Google Scholar : PubMed/NCBI

49 

Wang H, Haiman CA, Burnett T, Fortini BK, Kolonel LN, Henderson BE, Signorello LB, Blot WJ, Keku TO, Berndt SI, et al: Fine-mapping of genome-wide association study-identified risk loci for colorectal cancer in African Americans. Hum Mol Genet. 22:5048–5055. 2013. View Article : Google Scholar : PubMed/NCBI

50 

Hong SN, Park C, Kim JI, Kim DH, Kim HC, Chang DK, Rhee PL, Kim JJ, Rhee JC, Son HJ and Kim YH: Colorectal cancer-susceptibility single-nucleotide polymorphisms in Korean population. J Gastroenterol Hepatol. 30:849–857. 2015. View Article : Google Scholar : PubMed/NCBI

51 

Moazeni-Roodi A, Ghavami S, Ansari H and Hashemi M: Association between the flap endonuclease 1 gene polymorphisms and cancer susceptibility: An updated meta-analysis. J Cell Biochem. 120:13583–13597. 2019. View Article : Google Scholar : PubMed/NCBI

52 

Chou AK, Shen MY, Chen FY, Hsiao CL, Shih LC, Chang WS, Tsai CW, Ying TH, Wu MH, Huang CY and Bau DT: The association of flap endonuclease 1 genotypes with the susceptibility of endometriosis. Cancer Genomics Proteomics. 14:455–460. 2017.PubMed/NCBI

53 

Kubiritova Z, Gyuraszova M, Nagyova E, Hyblova M, Harsanyova M, Budis J, Hekel R, Gazdarica J, Duris F, Kadasi L, et al: On the critical evaluation and confirmation of germline sequence variants identified using massively parallel sequencing. J Biotechnol. 298:64–75. 2019. View Article : Google Scholar : PubMed/NCBI

54 

Trost B, Engchuan W, Nguyen CM, Thiruvahindrapuram B, Dolzhenko E, Backstrom I, Mirceta M, Mojarad BA, Yin Y, Dov A, et al: Genome-wide detection of tandem DNA repeats that are expanded in autism. Nature. 586:80–86. 2020. View Article : Google Scholar : PubMed/NCBI

55 

Mahmood S, Sivoňová M, Matáková T, Dobrota D, Wsólová L, Dzian A, et al: Association of EGF and p53 gene polymorphisms and colorectal cancer risk in the Slovak population. Cent Eur J Med. 9:405–416. 2014.

56 

Škereňová M, Halašová E, Matáková T, Jesenská L, Jurečeková J, Šarlinová M, Čierny D and Dobrota D: Low variability and stable frequency of common haplotypes of the tp53 gene region in colorectal cancer patients in a Slovak population. Anticancer Res. 37:1901–1907. 2017. View Article : Google Scholar : PubMed/NCBI

57 

Kašubová I, Kalman M, Jašek K, Burjanivová T, Malicherová B, Vaňochová A, Meršaková S, Lasabová Z and Plank L: Stratification of patients with colorectal cancer without the recorded family history. Oncol Lett. 17:3649–3656. 2019.PubMed/NCBI

58 

Jia WH, Zhang B, Matsuo K, Shin A, Xiang YB, Jee SH, Kim DH, Ren Z, Cai Q, Long J, et al: Genome-wide association analyses in East Asians identify new susceptibility loci for colorectal cancer. Nat Genet. 45:191–196. 2013. View Article : Google Scholar : PubMed/NCBI

59 

Peters U, Jiao S, Schumacher FR, Hutter CM, Aragaki AK, Baron JA, Berndt SI, Bézieau S, Brenner H, Butterbach K, et al: Identification of genetic susceptibility loci for colorectal tumors in a genome-wide meta-analysis. Gastroenterology. 144:799–807. 2013. View Article : Google Scholar : PubMed/NCBI

60 

Whiffin N, Hosking FJ, Farrington SM, Palles C, Dobbins SE, Zgaga L, Lloyd A, Kinnersley B, Gorman M, Tenes A, et al: Identification of susceptibility loci for colorectal cancer in a genome-wide meta-analysis. Hum Mol Genet. 23:4729–4737. 2014. View Article : Google Scholar : PubMed/NCBI

61 

Houlston RS, Cheadle J, Dobbins SE, Tenesa A, Jones AM, Howarth K, Spain SL, Broderick P, Domingo E, Farrington S, et al: Meta-analysis of three genome-wide association studies identifies susceptibility loci for colorectal cancer at 1q41, 3q26.2, 12q13.13 and 20q13.33. Nat Genet. 42:973–977. 2010. View Article : Google Scholar : PubMed/NCBI

62 

Orlando G, Law PJ, Palin K, Tuupanen S, Gylfe A, Hänninen UA, Cajuso T, Tanskanen T, Kondelin J, Kaasinen E, et al: Variation at 2q35 (PNKD and TMBIM1) influences colorectal cancer risk and identifies a pleiotropic effect with inflammatory bowel disease. Hum Mol Genet. 25:2349–2359. 2016. View Article : Google Scholar : PubMed/NCBI

63 

Schmit SL, Schumacher FR, Edlund CK, Conti DV, Raskin L, Lejbkowicz F, Pinchev M, Rennert HS, Jenkins MA, Hopper JL, et al: A novel colorectal cancer risk locus at 4q32.2 identified from an international genome-wide association study. Carcinogenesis. 35:2512–2519. 2014. View Article : Google Scholar : PubMed/NCBI

64 

Peters U, Hutter CM, Hsu L, Schumacher FR, Conti DV, Carlson CS, Edlund CK, Haile RW, Gallinger S, Zanke BW, et al: Meta-analysis of new genome-wide association studies of colorectal cancer risk. Hum Genet. 131:217–234. 2012. View Article : Google Scholar : PubMed/NCBI

65 

Tomlinson IPM, Webb E, Carvajal-Carmona L, Broderick P, Howarth K, Pittman AM, Spain S, Lubbe S, Walther A, Sullivan K, et al: A genome-wide association study identifies colorectal cancer susceptibility loci on chromosomes 10p14 and 8q23.3. Nat Genet. 40:623–630. 2008. View Article : Google Scholar : PubMed/NCBI

66 

Tenesa A, Farrington SM, Prendergast JGD, Porteous ME, Walker M, Haq N, Barnetson RA, Theodoratou E, Cetnarskyj R, Cartwright N, et al: Genome-wide association scan identifies a colorectal cancer susceptibility locus on 11q23 and replicates risk loci at 8q24 and 18q21. Nat Genet. 40:631–637. 2008. View Article : Google Scholar : PubMed/NCBI

67 

Tomlinson IPM, Carvajal-Carmona LG, Dobbins SE, Tenesa A, Jones AM, Howarth K, Palles C, Broderick P, Jaeger EEM, Farrington S, et al: Multiple common susceptibility variants near BMP pathway loci GREM1, BMP4, and BMP2 explain part of the missing heritability of colorectal cancer. PLoS Genet. 7:e10021052011. View Article : Google Scholar : PubMed/NCBI

68 

COGENT Study, ; Houlston RS, Webb E, Broderick P, Pittman AM, Di Bernardo MC, Lubbe S, Chandler I, Vijayakrishnan J, Sullivan K, et al: Meta-analysis of genome-wide association data identifies four new susceptibility loci for colorectal cancer. Nat Genet. 40:1426–1435. 2008. View Article : Google Scholar : PubMed/NCBI

69 

Broderick P, Carvajal-Carmona L, Pittman AM, Webb E, Howarth K, Rowan A, Lubbe S, Spain S, Sullivan K, Fielding S, et al: A genome-wide association study shows that common alleles of SMAD7 influence colorectal cancer risk. Nat Genet. 39:1315–1317. 2007. View Article : Google Scholar : PubMed/NCBI

70 

Schumacher FR, Schmit SL, Jiao S, Edlund CK, Wang H, Zhang B, Hsu L, Huang SC, Fischer CP, Harju JF, et al: Genome-wide association study of colorectal cancer identifies six new susceptibility loci. Nat Commun. 6:71382015. View Article : Google Scholar : PubMed/NCBI

71 

Haiman CA, Le Marchand L, Yamamato J, Stram DO, Sheng X, Kolonel LN, Wu AH, Reich D and Henderson BE: A common genetic risk factor for colorectal and prostate cancer. Nat Genet. 39:954–956. 2007. View Article : Google Scholar : PubMed/NCBI

72 

Tomlinson I, Webb E, Carvajal-Carmona L, Broderick P, Kemp Z, Spain S, Penegar S, Chandler I, Gorman M, Wood W, et al: A genome-wide association scan of tag SNPs identifies a susceptibility variant for colorectal cancer at 8q24.21. Nat Genet. 39:984–988. 2007. View Article : Google Scholar : PubMed/NCBI

73 

Hutter CM, Slattery ML, Duggan DJ, Muehling J, Curtin K, Hsu L, Beresford SA, Rajkovic A, Sarto GE, Marshall JR, et al: Characterization of the association between 8q24 and colon cancer: Gene-environment exploration and meta-analysis. BMC Cancer. 10:6702010. View Article : Google Scholar : PubMed/NCBI

74 

Cui R, Okada Y, Jang SG, Ku JL, Park JG, Kamatani Y, Hosono N, Tsunoda T, Kumar V, Tanikawa C, et al: Common variant in 6q26-q27 is associated with distal colon cancer in an Asian population. Gut. 60:799–805. 2011. View Article : Google Scholar : PubMed/NCBI

75 

Wang H, Burnett T, Kono S, Haiman CA, Iwasaki M, Wilkens LR, Loo LW, Van Den Berg D, Kolonel LN, Henderson BE, et al: Trans-ethnic genome-wide association study of colorectal cancer identifies a new susceptibility locus in VTI1A. Nat Commun. 5:46132014. View Article : Google Scholar : PubMed/NCBI

Related Articles

Journal Cover

November-2021
Volume 22 Issue 5

Print ISSN: 1792-1074
Online ISSN:1792-1082

Sign up for eToc alerts

Recommend to Library

Copy and paste a formatted citation
x
Spandidos Publications style
Forgacova N, Gazdarica J, Budis J, Radvanszky J and Szemes T: Repurposing non‑invasive prenatal testing data: Population study of single nucleotide variants associated with colorectal cancer and Lynch syndrome. Oncol Lett 22: 779, 2021
APA
Forgacova, N., Gazdarica, J., Budis, J., Radvanszky, J., & Szemes, T. (2021). Repurposing non‑invasive prenatal testing data: Population study of single nucleotide variants associated with colorectal cancer and Lynch syndrome. Oncology Letters, 22, 779. https://doi.org/10.3892/ol.2021.13040
MLA
Forgacova, N., Gazdarica, J., Budis, J., Radvanszky, J., Szemes, T."Repurposing non‑invasive prenatal testing data: Population study of single nucleotide variants associated with colorectal cancer and Lynch syndrome". Oncology Letters 22.5 (2021): 779.
Chicago
Forgacova, N., Gazdarica, J., Budis, J., Radvanszky, J., Szemes, T."Repurposing non‑invasive prenatal testing data: Population study of single nucleotide variants associated with colorectal cancer and Lynch syndrome". Oncology Letters 22, no. 5 (2021): 779. https://doi.org/10.3892/ol.2021.13040