Introduction

Lynch syndrome or hereditary non-polyposis colorectal cancer, the most common form of inherited colorectal cancer (CRC), results from germline mutations within the genes of the mismatch repair (MMR) system, MSH2, MLH1, MSH6 and PMS2 (for a review, see Lynch et al1). In MMR mutation carriers, the main tumour risks are colorectal and endometrial cancers. The cumulative risk at 70 years of CRC has been estimated to be 47–78% in males2, 3, 4 and 30–57% in females,2, 3, 4 and the risk of endometrial cancer to be 25–61%.2, 3 Like in other Mendelian forms of cancer characterized by an incomplete penetrance, one of the main challenges is to identify modifier genetic factors that can modulate the mutation penetrance. Characterization of validated modifier genetic factors should have important clinical consequences in the future since it should then be possible to adapt the follow-up of MMR gene mutation carriers, in terms of nature and timing of investigations, according to modifier alleles.

The modifier genes that have been reported so far in Lynch syndrome have been characterized according to two strategies: several studies have been focused on genes whose implication in CRC development had previously been suggested. One of the first modifier genetic factors identified in Lynch syndrome corresponds to a CA repeat polymorphism present within the IGF1 promoter, 1 kb upstream from the transcriptional initiation site.5, 6 In a study performed in 121 MMR mutation carriers from 59 families, mainly of Caucasian origin, Zecevic et al5 reported that a CA repeat ≤17 was significantly associated with a higher CRC risk and an earlier age of tumour onset. This result was confirmed by an independent study performed in 443 Australian and Polish MMR mutation carriers originated from 269 distinct families.6 Another class of modifier genes reported in Lynch syndrome corresponds to genes encoding xenobiotic metabolizing enzymes involved in environmental carcinogen metabolism. In 129 subjects of South African origin and harbouring the same MLH1 missense mutation, males harbouring the null genotype for the GSTT1 and GSTM1 genes developed cancer earlier than the males harbouring the other genotypes.7 This effect of the GSTT1 and GSTM1 null genotypes was not confirmed in a second study including 257 MMR mutation carriers from 130 families.8 Nevertheless, the authors reported that subjects heterozygous for the CYP1A1 rs1048943 SNP (c.1384A>G; p.Ile462Val) developed CRC earlier than individuals with the homozygous wild-type genotype and that subjects heterozygous for this polymorphism and an additional SNP rs4646903 (Msp1; g.6235T>C) had an increased CRC risk.8

The second strategy, which has recently allowed the detection of modifier genes in Lynch patients, originated from the numerous genome-wide association studies reporting SNPs associated with CRC risk in the general population. Starting from the hypothesis that SNPs acting as risk factors for CRC in the general population might act as risk modifiers in patients harbouring a highly penetrant mutation, Wijnen et al9 recently reported a significant association of CRC risk with rs16892766 (8q23.3) and rs3802842 (11q23.1) in 675 Dutch MMR mutation carriers from 127 families.

In this study, we investigated the impact of these different genetic factors on CRC risk in a large series of French MMR mutation carriers.

Patients and methods

Patients

The study included 748 unselected patients (Table 1) derived from 359 families with Lynch syndrome and recruited from Rouen (n=494) and Lille (n=254) university hospitals in France. All these individuals were confirmed carriers of a deleterious mutation of MSH2, MLH1 or MSH6 and were Caucasian. Among these 748 mutation carriers, 329 (44%) had been diagnosed with CRC prior to inclusion in the study, with a mean age of CRC onset of 43 years (range: 18–82 years). Additionally, 51 mutation carriers (6.8%) had been diagnosed with another tumour belonging to the Lynch syndrome spectrum (endometrial carcinoma, tumour of the urinary tract, ovarian carcinoma, cancer of the stomach and of the small intestine) and 368 (49.2%) had developed no tumour at the time of inclusion.

Table 1 Clinical characteristics of MMR mutation carriers

Genotyping

The IGF1 CA repeat was PCR-amplified using dye-labelled primers and the length of the PCR products was determined after migration on an Applied Biosystems model 3130 Genetic Analyser (PE Applied Biosystems, Foster City, CA, USA) and analysed using the GeneMapper Analysis Software version 4.0 (Applied Biosystems). To calibrate the results, we sequenced, after cloning into a plasmid, the most frequent allele (n=19) and used the corresponding genomic DNA as a reference DNA. The rs16892766 (8q23.3), rs3802842 (11q23.1) SNPs and the rs1048943 (c.1384A>G; p.Ile462Val) and rs4646903 SNPs, both located within the CYP1A1 gene (15q24.1), were genotyped using SNaPshot multiplex assays based on primer extension with dye-labelled dideoxynucleotides (ABI PRISM SNaPshot Multiplex kit, Applied Biosystems). We also analysed, by SNaPshot multiplex analysis, the Novel 1 SNP located on 18q21, which had been shown to be associated with CRC risk in the general population, by altering SMAD7 expression.10 Labelled products were separated using a 25-min run on an ABI Prism 3100 DNA sequencer and data were analysed using the GeneMapper Analysis Software version 4.0 (Applied Biosystems). Primer sequences used for genotyping and conditions of the SNaPshot multiplex assays are available on request.

Statistical analyses

The genotype frequencies were tested for Hardy–Weinberg equilibrium. In all following analyses of time to CRC, first CRC onset was the event of interest, age (in years) was used at the time scale and retrospective follow-up started at birth for all subjects in the study and ended at either first CRC onset for incident CRC cases or age at last follow-up otherwise. These analyses were performed for all subjects and separately for males and females and were repeated for CRC cases only. The association of sex or mutated MMR gene with CRC risk was assessed using the log-rank test. The association of each polymorphism (ie, five SNPs and IGF1 CA repeat length) with CRC risk was assessed using the log-rank test and then the Cox proportional hazard model with stratification on gender (if applicable), mutated MMR gene (MSH2, MLH1, MSH6) and date of birth (<1940, 1940–1949, 1950–1959, 1960–1969, 1970–1979 and ≥1980) and with control for familial clustering. In the Cox model analyses of each of the five SNPs, an overall comparison of the three genotypes was performed and a hazard ratio (HR) and corresponding 95% confidence interval were estimated for homozygous and heterozygous subjects, relative to subjects homozygous for the wild-type allele. Moreover, alternative Cox model analyses of each of the five SNPs considered the number of alleles with the SNP (0, 1, 2) and fitted a trend producing an alternative test of association with CRC risk and an estimate of the CRC HR per mutated allele. Statistical analyses were performed using the SAS software version 9.1 (SAS Institute, Cary, NC, USA).

Results

Among the 748 MMR mutation carriers included in this study, 329 carriers (44%) developed CRC. As shown in Table 1, a significant difference was observed in the age at first CRC occurrence between males and females, with males developing CRC earlier, between MSH2 and MSH6 mutation carriers, between MLH1 and MSH6 mutation carriers, with MSH2 and MLH1 carriers developing CRC earlier, but not between MSH2 and MLH1 mutation carriers (log-rank test, P=0.26). The five SNPs analysed, located on 8q23.3, 11q23.1, 18q21 and within CYP1A1, were in Hardy–Weinberg equilibrium. The distribution of the IGF1 CA repeat alleles was the following: 12 repeats (n=8), 13–15 repeats (n=0), 16 repeats (n=2), 17 repeats (n=28), 18 repeats (n=85), 19 repeats (n=976), 20 repeats (n=262), 21 repeats (n=114) and 22 repeats (n=21) and alleles were dichotomized in ≤17 and ≥18 repeats for the statistical analyses, like in the previous studies.5, 6

In univariate analysis, age at first CRC occurrence did not vary significantly for the five SNPs considered or according to number of IGF1 CA repeats, whether for the whole sample or separately for males and females (Table 2). As indicated in Table 3, the multivariate Cox analysis performed in the total sample of mutation carriers revealed that, for rs16892766 (8q23.3), the homozygosity for the minor C allele present only in three mutation carriers was significantly associated with a decreased risk relative to subjects with AA alleles (HR=0.267, P=0.0271). For the CYPA1 SNPs, we observed a significant difference only in males for one of the SNPs analysed (rs4646903): the CC genotype, detected in four male mutation carriers, was associated with an increased risk relative to subjects with TT alleles (HR=3.496, P=0.0033). Considering that penetrance of MSH6 mutation was found in this study to be lower than that of MSH2 and MLH1 mutations, we restricted the log-rank and Cox analyses to MSH2 and MLH1 carriers and this restriction did not modify the results (data not shown).

Table 2 Allelic frequency of genetic variants and assessment of the associated colorectal risk by the log-rank test
Table 3 Evaluation in 748 MMR mutation carriers of the CRC risk associated with 8q23.3, 11q23.1, 18q21, CYP1A1 and IGF1 variants

As colonoscopy had not been performed in all the unaffected mutation carriers at the time of this study, we then focused the statistical analyses on the 329 MMR mutation carriers with CRC, in order to avoid a bias in the phenotypic evaluation. As indicated in Table 4, the unrestricted multivariate Cox analysis showed that the only remaining significant association was a decreased risk associated with the 8q23.3 CC genotype, but it should be emphasized that, because this genotype was only carried by two CRC cases, no HR estimate could be produced.

Table 4 Evaluation in 329 MMR mutation carriers with CRC of the risk associated with 8q23.3, 11q23.1, 18q21, CYP1A1 and IGF1 variants

Discussion

Like in other cases of Mendelian predisposition to cancer with an incomplete penetrance, the characterization of modifier loci for CRC in Lynch syndrome might allow identification of a subset of mutation carriers who could benefit from a reinforced tumour detection program, such as annual chromocolonoscopy from 20 years of age. This prospect prompted us to evaluate in a large series of 748 MMR carriers the risk conferred by the allelic variants that have recently been reported as risk modifiers for CRC in MMR mutation carriers. Within this series, 44% of the MMR mutation carriers had developed CRC at the time of the analysis; the median age of CRC onset was 43 years among them and we observed a significant difference of CRC penetrance between males and females and between MSH2 and MLH1 mutation versus MSH6 mutation carriers. These results are in agreement with the published studies, which have estimated that the cumulative CRC risk at 70 years is higher in males than in females2, 3, 4 and is lower in MSH6 mutation carriers as compared with MSH2 or MLH1 mutation carriers.2, 3, 11, 12 These findings are in favour of the representativeness of our series with respect to Lynch syndrome.

The distribution of IGF1 CA repeat length was similar to that reported by Zecevic et al4 in 121 MMR mutation carriers, the frequency of the major allele (19 repeats) being 65 and 58% in the two studies, respectively, although we found a lower proportion of alleles with length ≤17 (2.5 versus 5.1%). In the study published by Zecevic et al5 on 121 MMR mutation carriers, the presence of a CA repeat ≤17, detected in 12 subjects, was found to be associated with a higher CRC risk (HR 2.36) and an earlier median age of CRC onset (44 versus 56.5 years), but our study, performed on a larger series, did not allow us to reproduce these results.

For the rs1048943 (c.1384A>G; p.Ile462Val) and rs4646903 SNPs, located within the CYP1A1 gene, we found respective allelic frequencies for the minor allele of 0.03 and 0.12, which correspond to the HapMap values estimated in the European population. Pande et al8 had reported a higher MAF (0.074) for rs1048943 in 257 MMR mutation carriers, but it should be underlined that these carriers were from different ethnicities and included Hispanic and Asian subjects, in whom the allelic frequency of the G allele is known to be higher. In this series of 257 mutation carriers, CRC free survival statistically differed according to rs1048943 genotype, with a 1.78-fold increase of CRC hazard associated with the AG genotype, a result that we did not reproduce in our series of 748 MMR mutation carriers from Caucasian origin. For the second CYPA1 SNP (rs4646903), Pande et al8 had reported that the TC genotype was associated with an increased hazard for earlier CRC and we detected no significant difference in our study for this genotype (Table 3). In contrast, we found that the CC genotype was associated in males with a 3.5-fold increased HR, but this observation should be interpreted with utmost caution, considering the limited number of male mutation carriers presenting this genotype (n=4). For rs16892766 (8q23.3) and rs3802842 (11q23.1) SNPs, we observed respective allele frequencies of 7 and 28%: these frequencies had been, respectively, estimated in the European population at 7%13 and 29%,14 and in 675 Dutch MMR mutation carriers at 10 and 24.9%.9 We did not detect in our series a modification of CRC risk in mutation carriers harbouring the 11q23.1 rs3802842 CC genotype, as previously reported in females. Among the 675 Dutch mutation carriers, the 8q23.3 rs16892766 CC genotype, detected in 9 subjects, was found to be associated with a 2.16-fold CRC risk increase.9 In contrast, we found in our study (Table 3) that this genotype was significantly associated with a decreased CRC risk (HR=0.267), and this significant association was the only one remaining when we restricted the analysis to affected carriers, but the small number of subjects harbouring this genotype (n=3) and presenting a CRC (n=2) must be highlighted.

During the submission of this study, Talseth-Palmer et al15 reported that in MLH1 carriers, but not in MSH2 carriers, the 11q23.1 CC and 8q23.3 AC genotypes were associated with an increased risk, but this significant association detected in 373 Australian mutation carriers was not found in 311 Polish mutation carriers analysed in the same study. We did not replicate either of these associations in the French sample. Indeed, the only association that we detected for these two SNPs in the 267 French MLH1 mutation carriers was a decrease of CRC risk with the 8q23.3 CC genotype (P<0.0001), but as this genotype was carried by only one CRC case, it was not possible to estimate the corresponding HR and due caution should be applied to interpret this finding.

To increase the power of statistical analysis of the 8q23.3 and 11q23.1 SNPs’ impact on CRC risk in MMR mutation carriers, we performed a meta-analysis based on the study published by Wijnen et al9 and our results, corresponding to 1423 MMR mutation carriers (Figure 1). In univariate analysis, age at first CRC occurrence varied significantly in male carriers according to the 8q23.3 genotype (AA, AC or CC; P=0.0265 using Fisher's standard method to combine P-values from two log-rank tests), with increased CRC risk associated with the AC genotype and decreased risk associated with the CC genotype. In multivariate Cox analysis, as shown on the Forest plot for the total sample (Figure 1), CRC risk did not vary significantly according to the 8q23.3 genotype. In males, the multivariate Cox analysis showed that the 8q23.3 CC genotype was associated with a decrease in CRC risk (combined fixed-effect weighted sum HR=0.61 relative to AA, P=0.0135). It remains, however, that there were very few males and even fewer male CRC cases carrying the CC genotype in both studies and thus the results of this meta-analysis should be interpreted with caution. Univariate and multivariate Cox analyses did not reveal any significant association between CRC risk and the 11q23.1 genotype (Figure 1).

Figure 1
figure 1

Forest plots of the meta-analysis of the CRC risk associated with 8q23.3 and 11q23.1 variants in 1423 MMR mutation carriers. CI, confidence interval; HR, hazard rate; LCL, lower 95% confidence limit; UCL, upper 95% confidence limit. For each panel, the last two columns display the number of subjects with colorectal cancer and the indicated genotype, followed by the overall number of subjects with this genotype.

In conclusion, the evaluation in a series of 748 MMR mutation carriers representative of Lynch syndrome of the IGF1 CA repeat, CYP1A1, 8q23.3 and 11q23.1 SNPs, previously reported as modifier factors for CRC risk in MMR mutation carriers, did not allow us to reproduce the previously published results. The only significant association that we detected both in the whole sample and in the unbiased affected mutation carrier sample was between the 8q23.3 CC genotype and a decreased CRC risk in males, whereas previously published studies had reported an increased CRC risk associated either with this genotype9 or with the heterozygous AC genotype.15 This discrepancy might be explained by the genetic heterogeneity between the different analysed populations. Considering that, in MMR mutation carriers, cancer development requires somatic mutation of several target genes containing repeated sequences, it is likely that numerous polymorphisms affect either the DNA repair or target genes and/or that environmental factors having a role in colorectal carcinogenesis will modify the penetrance of MMR mutations. It is also possible that a key factor modulating the penetrance is the allelic heterogeneity of the MMR mutations. Finally, another explanation of this discrepancy between our results and previous studies is that the previously reported associations were spurious associations due to too small samples. Therefore, we conclude that, from the medical point of view, genotyping of these polymorphisms is not useful to evaluate CRC risk in MMR mutation carriers and to optimize their clinical follow-up.