Exome sequencing in a Swedish family with PMS2 mutation with varying penetrance of colorectal cancer: investigating the presence of genetic risk modifiers in colorectal cancer risk : European Journal of Cancer Prevention

Secondary Logo

Journal Logo

Research article

Exome sequencing in a Swedish family with PMS2 mutation with varying penetrance of colorectal cancer: investigating the presence of genetic risk modifiers in colorectal cancer risk

Bryant, Patricka,b; Walton Bernstedt, Sophiec,d; Thutkawkorapin, Jessadaa,e; Backman, Ann-Sofiec,f; Lindblom, Annikaa; Lagerstedt-Robinson, Kristinaa,g

Author Information
European Journal of Cancer Prevention ():10.1097/CEJ.0000000000000769, December 06, 2022. | DOI: 10.1097/CEJ.0000000000000769

Abstract

Introduction

Lynch syndrome is caused by germline mutations in the mismatch repair (MMR) genes MLH1, MSH2, MSH6, PMS2 and EPCAM. It is characterised by early onset, familial accumulation of cancer, predominantly colorectal cancer (CRC) and endometrial cancer (EC), but also a wide range of tumours such as gastric, small bowel, ovarian, pancreatic and uroepithelial cancer. Initially, the PMS2 gene was not considered to cause Lynch syndrome due to low penetrance of cancer compared to other MMR genes (Liu et al., 2001). Biallelic PMS2 mutations have been associated with early onset cancer in sites such as the brain and the colorectum (De Rosa et al., 2000). It was, therefore, suggested that mutations in the PMS2 gene are inherited in an autosomal recessive fashion (de Vos et al., 2005) and that biallelic germline inactivation is required for cancer predisposition (Nakagawa et al., 2004). It has later been shown (Hendriks et al., 2006) that heterozygous mutations in PMS2 cause Lynch syndrome in families where the analysis for MLH1, MSH2 and MSH6 have been negative for Lynch syndrome. However, the penetrance of cancer in PMS2 carriers is still not fully elucidated. A recent register paper (Møller et al., 2017) found that PMS2 mutation carriers were at a small increased risk for CRC and EC but not at significant risk for any other Lynch syndrome-associated cancer. Heterozygous PMS2 mutation carriers were at increased risk for CRC (cumulative risk to age 80 years of 13% for males and 12% for females) and endometrial cancer (13%), compared with the general population (6.6, 4.7 and 2.4%, respectively) (Ten Broeke et al., 2018). CRC risk has been displayed to vary widely between members of families having PMS2 mutations, demonstrated by the wide age range of initial CRC diagnosis (mean 52 years; range 26–86 years) as well as the difference in the mean age of CRC (P < 0.001) between probands with PMS2 mutation (mean 47 years; range 26–68 years) and other family members with a PMS2 mutation (mean, 58 years; range, 31–86 years) (Ten Broeke et al., 2015). This heterogeneity in risk is consistent with the lifetime CRC risk distribution for MLH1 carriers [34% (25–50%, 95% confidence interval {CI}) for males and 36% (25–51%, 95% CI) for females] and MSH2 carriers [47% (36–60%, 95% CI) for males and 37% (27–50%, 95% CI) for females] (Dowty et al., 2013), suggesting the presence of risk modifiers. These modifiers can be both genetic variations and lifestyle factors. To date, two main classes of cancer susceptibility variants have been detected: rare moderate-penetrance variants with risk allele frequencies <2% and odds ratios (OR) >2.0 and common low-penetrance alleles with risk allele frequencies >5% and OR <1.5 (Sud et al., 2017). However, it is considered likely that the spectrum of penetrance and frequency of risk alleles for many cancers occurs on a continuum. This would infer the possibility for any potential genetic modifiers to vary in both frequency and effect size (Sud et al., 2017). Lynch syndrome has been estimated to account for 2.2% of all cases of CRC (Hampel et al, 2005). In Sweden, CRC is the third most common cancer with approximately 6000 cases annually (Socialstyrelsen, 2018). It is estimated that the narrow sense heritability, the component of phenotypic variance that is accounted for by additive genetic variance, in the Swedish population is 13 and 12% for colon and rectal cancer, respectively. The estimated narrow sense heritability explained by common (>1%) genome-wide association studies (GWAS) CRC SNPs (single nucleotide polymorphism) is 0.65%, whereas all common GWAS SNPs are estimated to explain 7.42% (Jiao et al, 2014).

The Swedish family studied here carries a PMS2 missense mutation [c.2113G>A, p.(Glu705Lys)] (Miyaki et al, 1997; Lagerstedt-Robinson et al, 2007). This mutation has been found in other families with Lynch syndrome with similar clinical phenotype and does not segregate perfectly in families (Lagerstedt-Robinson et al, 2016). The studied family displays high penetrance for CRC in PMS2 carriers, as well as healthy carriers and noncarriers with late-onset CRC. Because of this obvious variation in penetrance, it is suitable for investigating the possibility of genetic risk modifiers. Whole exome sequencing was performed on five family members to identify potential genetic modifiers.

Objectives

To investigate the presence of genetic risk modifiers in a Swedish family with a PMS2 missense mutation [c.2113G>A, p.(Glu705Lys)] (Miyaki et al, 1997; Lagerstedt-Robinson et al, 2007) using whole exome sequencing data.

Methods

Ethics approval

The Regional Ethics Examination Board of Stockholm approved the study (ethical permits 2008/125-31.2, 2019-0233, 2020-00953).

Family 119

Whole exome sequencing was performed in the individuals I:V, II:I, III:II, II:V and II:II of family 119 (Fig. 1) on genomic DNA isolated from blood samples. As shown in Fig. 1, II:II does not carry the PMS2 mutation and got cancer at age 72 years, whereas the rest of the family members all carry the PMS2 mutation. II:IV got cancer at age 79 years. I:V got cancer at age 61 years and his two sons II:I and III:II at age 53 years.

F1
Fig. 1:
Pedigree of family 119².²Family member II:IV got cancer at age 79 years. The other mutation carriers had earlier onset cancer. A scenario where the remaining four family members share a genetic modifier that II:IV lacks was hypothesised (scenario 1). A second scenario where all sequenced family members share a genetic modifier that the sons to II:V (III:I and III:II) are heterozygous for was also investigated (scenario 2).

Exome sequencing of colorectal cancer samples

Genomic DNA was quantified using a Qubit Fluorometer (Life Technologies, Carlsbad, CA, USA). Sequencing libraries were prepared according to the TruSeq DNA Sample Preparation Kit EUC 15005180 or EUC 15026489 (Illumina, San Diego, CA, USA). Briefly, 1-1.5 µg of genomic DNA was fragmented using the Covaris 400 bp protocol (Covaris Inc, Woburn, MA, USA). After fragmentation, all samples were subjected to end-repair, A-tailing and adaptor ligation of Illumina Multiplexing PE adaptors. An additional gel-based size selection step was performed and the adapter-ligated fragments were subsequently enriched by PCR followed by purification using Agencourt AMPure Beads (Beckman Coulter, Pasadena, CA, USA). Exome capture was performed by pre-pooling equimolar amounts and performing enrichment in 5- or 6-plex reactions according to the TruSeq Exome Enrichment Kit Protocol (EUC 15013230). Library size was checked on a Bioanalyzer High Sensitivity DNA chip (Agilent Technologies, Santa Clara, CA, USA) whereas concentration was calculated by quantitative PCR. The pooled DNA libraries were clustered on a cBot instrument (Illumina) using the TruSeq PE Cluster Kit v3. Paired-end sequencing was performed for 100 cycles using a HiSeq 2000 instrument (Illumina) with TruSeq SBS Chemistry v3, according to the manufacturer’s protocol. Base calling was performed with real time analysis (RTA) (1.12.4.2 or 1.13.48) and the resulting binary base call (BCL) files were filtered, de-multiplexed and converted to FASTQ format using CASAVA 1.7 or 1.8 (Illumina). The sequencing was performed at an average coverage of 100×.

Bioinformatics workflow

Sequencing reads were aligned to the reference genome GRCh37 using BWA (Li and Durbin, 2010). Aligned reads were sorted and PCR-duplicated reads were removed using Picard (http://broadinstitute.github.io/picard/). The calculation of mapping and enrichment statistics were performed with Picard and genome analysis toolkit (GATK) (McKenna et al., 2010). Variants were called using GATK by following the best practice procedure implemented at the Broad Institute (DePristo et al., 2011). Variant Quality Score Recalibration from GATK was used for quality control of the variants. Variant annotation was performed by ANNOVAR (Wang et al., 2010). The annotated information includes RefSeq gene annotation (O’Leary et al., 2016) and dbSNP rs number (Sherry et al., 2001). Background allele frequencies are from ExAC/gnomAD (Karczewski et al., 2020), 200 Danes (Li et al., 2010), SweGen (Ameur et al., 2017) and 1000 Genomes Project allele frequencies (Auton et al., 2015). In silico predictors used for predicting pathogenic effects include sorting intolerant from tolerant (SIFT) (Kumar et al., 2009), PolyPhen2 (Adzhubei et al., 2010), Phylop (Cooper et al., 2005), likelihood ratio test (LRT) (Chun and Fay, 2009), Mutation Taster (Schwarz et al., 2010), Mutation Assessor (Reva et al., 2011), FATHMM (Shihab et al., 2015), GERP++ (Davydov et al., 2010), RadialSVM (Li et al., 2014), LR (Li et al., 2014), and CADD (Kircher et al., 2014).

Variant selection

Variants obtained from sequencing were first filtered using the minor allele frequencies (MAFs) ≤0.2 in population databases (SweGen, gnomAD). Only the resulting nonsynonymous coding exonic variants were investigated (see Supplementary Table A1 and A2). Because II:IV got cancer at age 79 years, whereas the other carriers had earlier onset, a scenario where the remaining four family members share a genetic modifier that II:IV lacks was hypothesised (scenario 1). A second scenario where all sequenced family members share a genetic modifier that the sons to II:V (III:I and III:II) are heterozygous for was also investigated (scenario 2). The homozygous alternative was not investigated because inbreeding is unlikely and the probability thereby low of the sons being of such zygosity. The criteria for scenario 1 and 2 were then selected for generating two different result sets. The genes present in these result sets were investigated by literature search using the databases Pubmed (www.ncbi.nlm.nih.gov/pubmed) and Genecards (www.genecards.org). The genes previously reported to be associated with CRC development, other cancer development or genes influencing replication, transcription (e.g. transcription factor binding), translation, apoptosis, T-cell signalling, immune modulation, growth factor regulation, tumour suppressor ability or act as metastasis drivers were selected and grouped. All other genes containing variants were disregarded, a subset was of unknown importance, whereas the remainder were considered likely not to be influencing (for example olfactory receptors, such as OR10J1, see Supplementary Table A3, appendix section, Supplemental digital content 1, https://links.lww.com/EJCP/A366).

The presence of the sequence variants in the selected genes was evaluated in 66 previously investigated CRC families. One previously investigated individual was randomly selected from each family and the presence of each specific sequence variant in any zygosity state was reported. Two quotas were also calculated by dividing the frequency in the CRC families with the largest reported frequency in the mentioned databases and with the SweGen frequency. The variants found to have a quota above 1.10 and thereby be more common in the CRC families and displaying a frequency of at least 0.5 % in SWEGEN were selected. After these filtering steps, the variants in the first scenario were evaluated using association analysis in 724 cases and 711 controls (Table 1).

Table 1 - Results for scenario 1¹
Variant avSNP150 MAF
SWEGEN Gene #cases/
controls Allele count Chi2 P val OR
Chr3:184020542G>T rs11545169 0.163 PSMD2, c.939G>T, p.(Glu313Asp) - - -
Chr5:139905844G>A rs1051309 0.028 ANKHD1-EIF4EBP3, c.4756G>A, p.(Gly1589Ser) 705/673 Ctrl
aG:1309/A:37
Case G:1372/A:40 0.99937 1.03
Chr5:149389968G>A rs142114383 0.021 HMGXB3, c.607G>A, p.Glu203Lys - - -
Chr5:156376703A>G rs6873053 0.9265 TIMD4,
c.719T>C, p.(Val240Ala) 704/681 Ctrl G:1249/A:113
Case G:1280/A:128 0.90781 1.11
Chr6:149783095T>C rs61997220 0.0585 ZC3H12D,
c.317A>G, p.(Lys106Arg) 701/672 0.88615 1.14
Chr7:27135096G>T rs45571645 0.0265 HOXA1,
c.436C>A, p.(His146Asn) 698/670 0.99929 1.03
Chr7: 27196116C>T rs78410337 0.014 HOXA7, c.49G>T, p.Gly17Trp - - -
Chr7:150557578G>A rs11771771 0.117 AOC1,
c.1903G>A, p.(Val635Met) 686/664 0.12024 0.73
Chr8:17731611T>A rs35431851 0.0205 FGL1,
c.419A>T, p.(Tyr140Phe) 712/692 0.96972 0.88
Chr8:95272605G>C rs2170363 0.1275 GEM,
c.127C>G, p.(Arg43Gly) 707/644 0.75167 0.89
Chr9:33264393_33264410del rs574287414 0.005 BAG1,c.263_280del, p.(Arg88_Thr93del) 528/504 0.98782 1.22
Chr9:140110555G>A rs62587579 0.043 NDOR1,
c.1591G>A, p.(Val531Ile) 709/697 0.94293 1.13
Chr11:102272839C>A rs11547915 0.0585 TMEM123, c.256G>T, p.Val86Phe - - -
¹After these filtering steps, the variants in the first scenario were evaluated using association analysis in 724 cases and 711 controls (see Table 1).

Association analysis

TaqMan was performed for 12 of the 13 selected candidate sequence variants. One of the variants (BAG1, rs574287414), consisting of a deletion of 18 bp was analysed using fragment length analysis. The 724 cases and 711 controls were divided into two 384 plates each. TaqMan was performed as previously described. The amount of DNA per well was between 100 and 250 ng (Supplementary Table A4, appendix section, Supplemental digital content 1, https://links.lww.com/EJCP/A366). Fragment analysis was performed using 250 ng DNA per reaction (Supplementary Table A5, appendix section, Supplemental digital content 1, https://links.lww.com/EJCP/A366). PCR-products were run on an ABI 3130 genetic analyzer and 5 µl Milli-Q purified water was added to each well and the plate was centrifuged for 20 s. Sizing of the PCR fragments was performed by KI Gene, Center for Molecular Medicine, Karolinska Institute, Solna, Sweden.

Results

In total five individuals within the family were selected for whole exome sequencing; two brothers with high age and cancer with (II:IV) and without (II:II) the PMS2 mutation in addition to three individuals with early onset CRC and the PMS2 mutation (II:V, III:I and III:II). It was hypothesised that II:V, II:I and III:I who has early onset CRC and the PMS2 mutation should have the modifier, as well as II:II who has late onset but lacks the PMS2 mutation. II:IV who has late onset as well as the PMS2 mutation should lack the modifier. The researchers identified 13 candidate variants which fulfil these criteria in genes HMGXB3, PSMD2, TMEM123, ANKHD1-EIF4EBP3, TIMD4, ZC3H12D, HOXA1, HOXA7, AOC1, FGL1, GEM, BAG1 and NDOR1 (Table 1). The following variants could not be evaluated due to technical difficulties: rs142114383 (HMGXB3), rs11545169 (PSMD2), rs11547915 (TMEM123) and rs78410337 (HOXA7). The SNP rs78410337 had two underlying SNPs in the probe design area, why two different probes were used for evaluation. No significant effects were found for any controlled variant. Odds ratios above 1 were observed for rs6873053 (TIMD4), rs62587579 (NDOR1), rs61997220 (ZC3H12D) and rs574287414 (BAG1) (Table 1).

Discussion

The Swedish family investigated here, family 119, harbours many CRC cases as well as a PMS2 missense mutation [c.2113G>A, p.(Glu705Lys)] (Fig. 1). Because the PMS2 mutation does not segregate perfectly with the CRC within the family, a modifying genetic effect is plausible. This makes the family suitable for studying the incomplete penetrance of PMS2.

The selected family has a promising structure for studying the incomplete penetrance of PMS2 and the possibility of a genetic modifier because it contains both healthy individuals and individuals with the disease of various ages, both with and without PMS2 mutations present. If such a modifier as the one hypothesised here does exist, it should likely be found. However, no significant effects were found in this study.

The ORs of importance are 1.11, 1.13, 1.14 and 1.22 for rs68730353, rs62587579, rs61997220 and rs574287414, respectively (Table 1). We note that these ORs are low. rs68730353 is in gene TIMD4, which has been found to be cancer-related (Yano et al., 2017). The gene coding for NDOR1 (rs62587579) has been found to be upregulated in colonic adenoma cells (Paine et al., 2000) and ZC3H12D (rs61997220) has been reported to be a possible tumour suppressor gene (Minagawa et al, 2007). The 18 bp deletion in BAG1 (rs574287414) displays the highest OR of 1.22. BAG1 (BAG cochaperone 1) is a protein-coding gene and binds to the oncogene BCL2. BCL2 is an anti-apoptotic membrane protein and BAG1 enhances the antiapoptotic effects of BCL2 (NCBI’s Gene, 2022). Numerous protein isoforms are encoded by BAG1 mRNA. Nuclear BAG1 expression has been linked to a poor prognosis in CRC and may be used as a predictive factor for distant metastasis by means of immunohistochemical staining (Kikuchi et al., 2002). Survival tended to be shorter for patients with BAG1 positive nuclei in tumours than those without nuclear BAG-1-staining (P = 0.011). By using gamma radiation or a vitamin D analogue to induce apoptosis, it has been demonstrated that apoptosis was preceded by a decrease in nuclear and an increase in cytoplasmic BAG1 expression (Barnes et al., 2005). The expression of the nuclear BAG-1L isoform enhanced cellular survival after induced apoptosis. Cytoplasmic BAG-1S fused with a nuclear localisation signal which seemed to protect cells against induced apoptosis. Knockdown of BAG-1 in a CRC cell line has been demonstrated to inhibit NF-κB (nuclear factor-κB) transcriptional activity (Clemo et al., 2008). It has therefore been proposed that BAG-1 promotes colorectal tumour cell survival through increased NF-κB activity. A link between BAG1 and TGFB1 has also been shown; induction of BAG-1L caused suppression of TGFB1 mRNA in colorectal tumour cells (Skeen et al., 2013). Even though the 18 bp deletion is likely debilitating to the protein structure, regulatory relationships are often complicated, and one cannot exclude that the effect may be driving in combination with certain genetic backgrounds (Huang et al., 2014); it has not been explored how the BAG1 gene interacts with the PMS2 missense mutation in this family.

The four variants that could not be evaluated due to technical difficulties, rs142114383, rs11545169, rs11547915 and rs78410337, are related to carcinogenic processes. rs11545169 is in gene PSMD2 and may be related to the tumour necrosis factor signalling pathway (Tsurumi et al., 1996). rs142114383 in HMGXB3, which suppression is found to result in tumorigenesis in the colorectal epithelium (Sun et al., 2017). TMEM123 (rs11547915) is proposed to function as a cell surface receptor that mediates cell death (Huang et al., 2013) and rs78410337 in gene HOXA7 is related to acute myeloid leukaemia (Luo et al., 2018).

A possible explanation for why the selected variants do not display significant effects is that they do not impact the penetrance of CRC. Likely the relationship between CRC and the incomplete penetrance are much more complicated than hypothesised here. Even if effects were detected, there is a limit to their ability in explaining CRC heritability. Because the variants studied were considered in combination only with one PMS2 variant, this investigation is unlikely to explain much more of the complex heritability of CRC in the Swedish family studied here. Furthermore, only coding variants were considered, whereas most of the variation is situated within noncoding regions.

Limitations

Limitations include the small sample size and the rarity of the selected variants in the Swedish population.

Conclusion

In this study, whole exome sequencing was performed in a Swedish family carrying a PMS2 missense mutation [c.2113G>T;A, p.(Glu705Lys)]. Thirteen genetic sequence variants were further selected and analysed in a case-control study (724 cases and 711 controls). The variants with the highest ORs were in the genes TIMD4, NDOR1, ZC3H12D and BAG1. The most interesting sequence found was an 18 bp deletion rs68730353 in gene BAG1 and may play a role in colorectal tumorigenesis. It is not known how the BAG1 gene interacts with the PMS2 mutation in the Swedish family. The heritability of CRC is likely more complex and needs further work needs to be done to elucidate this concept. Further theoretical development and a better understanding of molecular biology are likely needed to reach the next major breakthrough in the heritability of CRC.

Acknowledgements

The authors would like to thank Anna-Lena Kastman for contributing with help regarding TaqMan analysis.

This work was mainly funded by grants to A.L. from the Swedish Research Council, the Swedish Cancer Society and the Cancer Research Funds of Radiumhemmet. All participants gave written consent to participate in the studies. All participants gave written consent to participate in the studies.

Conflicts of interest

There are no conflicts of interest.

References

Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, et al. (2010). A method and server for predicting damaging missense mutations. Nat Methods 7:248–249.
Ameur A, Dahlberg J, Olason P, Vezzi F, Karlsson R, Martin M, et al. (2017). SweGen: a whole-genome data resource of genetic variability in a cross-section of the Swedish population. Eur J Hum Genet 25:1253–1260.
Auton A, Brooks LD, Durbin RM, Garrison EP, Kang HM, et al.; The 1000 Genomes Project Consortium (2015). A global reference for human genetic variation. Nature 526:68–74.
Barnes JD, Arhel NJ, Lee SS, Sharp A, Al-Okail M, Packham G, et al. (2005). Nuclear BAG-1 expression inhibits apoptosis in colorectal adenoma-derived epithelial cells. Apoptosis 10:301–311.
Chun S, Fay JC (2009). Identification of deleterious mutations within three human genomes. Genome Res 19:1553–1561.
Clemo NK, Collard TJ, Southern SL, Edwards KD, Moorghen M, Packham G, et al. (2008). BAG-1 is up-regulated in colorectal tumour progression and promotes colorectal tumour cell survival through increased NF-kappaB activity. Carcinogenesis 29:849–857.
Cooper GM, Stone EA, Asimenos G, Green ED, Batzoglou S, Sidow A; NISC Comparative Sequencing Program. (2005). Distribution and intensity of constraint in mammalian genomic sequence. Genome Res 15:901–913.
Davydov EV, Goode DL, Sirota M, Cooper GM, Sidow A, Batzoglou S (2010). Identifying a high fraction of the human genome to be under selective constraint using GERP++. PLoS Comput Biol 6:e1001025.
De Rosa M, Fasano C, Panariello L, Scarano MI, Belli G, Iannelli A, et al. (2000). Evidence for a recessive inheritance of Turcot’s syndrome caused by compound heterozygous mutations within the PMS2 gene. Oncogene 19:1719–1723.
de Vos M, Hayward B, Bonthron DT, Sheridan E (2005). Phenotype associated with recessively inherited mutations in DNA mismatch repair (MMR) genes. Biochem Soc Trans 33:718–720.
DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, et al. (2011). A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet 43:491–498.
Dowty JG, Win AK, Buchanan DD, Lindor NM, Macrae FA, Clendenning M, et al. (2013). Cancer risks for MLH1 and MSH2 mutation carriers. Hum Mutat 34:490–497.
Hampel H, Frankel WL, Martin E, Arnold M, Khanduja K, Kuebler P, et al. (2005). Screening for the Lynch syndrome (hereditary nonpolyposis colorectal cancer). N Engl J Med 352:1851–1860.
Hendriks YM, Jagmohan-Changur S, van der Klift HM, Morreau H, van Puijenbroek M, Tops C, et al. (2006). Heterozygous mutations in PMS2 cause hereditary nonpolyposis colorectal carcinoma (Lynch syndrome). Gastroenterology 130:312–322.
Huang B, Zhou H, Lang X, Liu Z, Xiong F, Wang S (2014). Expression of BAG-1 is closely related to cell differentiation and TNM stage in esophageal cancer and its downregulation inhibits the proliferation and invasion of human esophageal carcinoma cells. Oncol Rep 32:1441–1446.
Huang Q, Zhang R, Zou Ly, Cao X, Chu X (2013). Cell death pathways in astrocytes with a modified model of oxygen-glucose deprivation. PLoS One 8:e61345.
Jiao S, Peters U, Berndt S, Brenner H, Butterbach K, Caan BJ, et al. (2014). Estimating the heritability of colorectal cancer. Hum Mol Genet 23:3898–3905.
Karczewski KJ, Francioli LC, Tiao G, Cummings BB, Alföldi J, Wang Q, et al.; Genome Aggregation Database Consortium (2020). The mutational constraint spectrum quantified from variation in 141,456 humans. Nature 581:434–443.
Kikuchi R, Noguchi T, Takeno S, Funada Y, Moriyama H, Uchida Y (2002). Nuclear BAG-1 expression reflects malignant potential in colorectal carcinomas. Br J Cancer 87:1136–1139.
Kircher M, Witten DM, Jain P, O’Roak BJ, Cooper GM, Shendure J (2014). A general framework for estimating the relative pathogenicity of human genetic variants. Nat Genet 46:310–315.
Kumar P, Henikoff S, Ng PC (2009). Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat Protoc 4:1073–1081.
Lagerstedt Robinson K, Liu T, Vandrovcova J, Halvarsson B, Clendenning M, Frebourg T, et al. (2007). Lynch syndrome (hereditary nonpolyposis colorectal cancer) diagnostics. J Natl Cancer Inst 99:291–299.
Lagerstedt-Robinson K, Rohlin A, Aravidis C, Melin B, Nordling M, Stenmark-Askmalm M, et al. (2016). Mismatch repair gene mutation spectrum in the Swedish Lynch syndrome population. Oncol Rep 36:2823–2835.
Li H, Durbin R (2010). Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics 26:589–595.
Li Q, Liu X, Gibbs RA, Boerwinkle E, Polychronakos C, Qu HQ (2014). Gene-specific function prediction for non-synonymous mutations in monogenic diabetes genes. PLoS One 9:e104452.
Li Y, Vinckenbosch N, Tian G, Huerta-Sanchez E, Jiang T, Jiang H, et al. (2010). Resequencing of 200 human exomes identifies an excess of low-frequency non-synonymous coding variants. Nat Genet 42:969–972.
Liu T, Yan H, Kuismanen S, Percesepe A, Bisgaard ML, Pedroni M, et al. (2001). The role of hPMS1 and hPMS2 in predisposing to colorectal cancer. Cancer Res 61:7798–7802.
Luo H, Wang F, Zha J, Li H, Yan B, Du Q, et al. (2018). CTCF boundary remodels chromatin domain and drives aberrant HOX gene transcription in acute myeloid leukemia. Blood 132:837–848.
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. (2010). The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res 20:1297–1303.
Minagawa K, Yamamoto K, Nishikawa S, Ito M, Sada A, Yakushijin K, et al. (2007). Deregulation of a possible tumour suppressor gene, ZC3H12D, by translocation of [email protected] in transformed follicular lymphoma with t(2;6)(p12;q25). Br J Haematol 139:161–163.
Miyaki M, Nishio J, Konishi M, Kikuchi-Yanoshita R, Tanaka K, Muraoka M, et al. (1997). Drastic genetic instability of tumors and normal tissues in Turcot syndrome. Oncogene 15:2877–2881.
Møller P, Seppälä T, Bernstein I, Holinski-Feder E, Sala P, Evans DG, et al. Mallorca Group (http://mallorca-group.eu) (2017). Cancer incidence and survival in Lynch syndrome patients receiving colonoscopic and gynaecological surveillance: first report from the prospective Lynch syndrome database. Gut 66:464–472.
Nakagawa H, Lockman JC, Frankel WL, Hampel H, Steenblock K, Burgart LJ, et al. (2004). Mismatch repair gene PMS2: disease-causing germline mutations are frequent in patients whose tumors stain negative for PMS2 protein, but paralogous genes obscure mutation detection and interpretation. Cancer Res 64:4721–4727.
NCBI. NCBI’s Gene [database on the Internet]. Available from: https://www.ncbi.nlm.nih.gov/gene/573. [Accessed 17 May 2022].
O’Leary NA, Wright MW, Brister JR, Ciufo S, Haddad D, McVeigh R, et al. (2016). Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation. Nucleic Acids Res 44:D733–D745.
Paine MJ, Garner AP, Powell D, Sibbald J, Sales M, Pratt N, et al. (2000). Cloning and characterization of a novel human dual flavin reductase. J Biol Chem 275:1471–1478.
Reva B, Antipin Y, Sander C (2011). Predicting the functional impact of protein mutations: application to cancer genomics. Nucleic Acids Res 39:e118.
Schwarz JM, Rödelsperger C, Schuelke M, Seelow D (2010). MutationTaster evaluates disease-causing potential of sequence alterations. Nat Methods 7:575–576.
Sherry ST, Ward MH, Kholodov M, Baker J, Phan L, Smigielski EM, Sirotkin K (2001). dbSNP: the NCBI database of genetic variation. Nucleic Acids Res 29:308–311.
Shihab HA, Rogers MF, Gough J, Mort M, Cooper DN, Day IN, et al. (2015). An integrative approach to predicting the functional effects of non-coding and coding sequence variation. Bioinformatics 31:1536–1543.
Skeen VR, Collard TJ, Southern SL, Greenhough A, Hague A, Townsend PA, et al. (2013). BAG-1 suppresses expression of the key regulatory cytokine transforming growth factor β (TGF-β1) in colorectal tumour cells. Oncogene 32:4490–4499.
Socialstyrelsen. Statistikdatabas för cancer [database on the Internet]. Available from: http://www.socialstyrelsen.se/statistik/statistikdatabas/cancer. [Accessed 19 October 2018].
Sud A, Kinnersley B, Houlston RS (2017). Genome-wide association studies of cancer: current insights and future perspectives. Nat Rev Cancer 17:692–704.
Sun LN, Xing C, Zhi Z, Liu Y, Chen LY, Shen T, et al. (2017). Dicer suppresses cytoskeleton remodeling and tumorigenesis of colorectal epithelium by miR-324-5p mediated suppression of HMGXB3 and WASF-2. Oncotarget 8:55776–55789.
Ten Broeke SW, van der Klift HM, Tops CMJ, Aretz S, Bernstein I, Buchanan DD, et al. (2018). Cancer risks for PMS2-associated Lynch syndrome. J Clin Oncol 36:2961–2968.
Tsurumi C, Shimizu Y, Saeki M, Kato S, Demartino GN, Slaughter CA, et al. (1996). cDNA cloning and functional analysis of the p97 subunit of the 26S proteasome, a polypeptide identical to the type-1 tumor-necrosis-factor-receptor-associated protein-2/55.11. Eur J Biochem 239:912–921.
Ten Broeke SW, Brohet RM, Tops CM, van der Klift HM, Velthuizen ME, Bernstein I, et al. (2015). Lynch syndrome caused by germline PMS2 mutations: delineating the cancer risk. J Clin Oncol 33:319–325.
Wang K, Li M, Hakonarson H (2010). ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res 38:e164.
Yano H, Motoshima T, Ma C, Pan C, Yamada S, Nakayama T, et al. (2017). The significance of TIMD4 expression in clear cell renal cell carcinoma. Med Mol Morphol 50:220–226.
Keywords:

colorectal cancer; gene exome sequencing; genetic risk modifier; Lynch syndrome; PMS2

Supplemental Digital Content

Copyright © 2022 The Author(s). Published by Wolters Kluwer Health, Inc.