Comparative analysis of Illumina Mouse Methylation BeadChip and reduced-representation bisulfite sequencing for routine DNA methylation analysis

Researching the murine epigenome in disease models has been hampered by the lack of appropriate and cost-effective DNA methylation arrays. Here we perform a comprehensive, comparative analysis between the Mouse Methylation BeadChip (MMB) and reduced-representation bisulfite sequencing (RRBS) in two murine models of colorectal carcinogenesis. We evaluate the coverage, variability, and ability to identify differential DNA methylation of RRBS and MMB. We show that MMB is an effective tool for profiling the murine methylome that performs comparably with RRBS, identifying similar differentially methylated pathways. Although choice of technology is experiment dependent and will be predicated on the underlying biology being probed, these analyses provide insights into the relative strengths and weaknesses of each approach.

Global Methylome Scores Correlate with Histological Subtypes of Colorectal Carcinoma and Show Different Associations with Common Clinical and Molecular Features

BACKGROUND

The typical methylation patterns associated with cancer are hypermethylation at gene promoters and global genome hypomethylation. Aberrant CpG island hypermethylation at promoter regions and global genome hypomethylation have not been associated with histological colorectal carcinomas (CRC) subsets. Using Illumina's 450 k Infinium Human Methylation beadchip, the methylome of 82 CRCs were analyzed, comprising different histological subtypes: 40 serrated adenocarcinomas (SAC), 32 conventional carcinomas (CC) and 10 CRCs showing histological and molecular features of microsatellite instability (hmMSI-H), and, additionally, 35 normal adjacent mucosae. Scores reflecting the overall methylation at 250 bp, 1 kb and 2 kb from the transcription starting site (TSS) were studied.

RESULTS

SAC has an intermediate methylation pattern between CC and hmMSI-H for the three genome locations. In addition, the shift from promoter hypermethylation to genomic hypomethylation occurs at a small sequence between 250 bp and 1 Kb from the gene TSS, and an asymmetric distribution of methylation was observed between both sides of the CpG islands (N vs. S shores).

CONCLUSION

These findings show that different histological subtypes of CRC have a particular global methylation pattern depending on sequence distance to TSS and highlight the so far underestimated importance of CpGs aberrantly hypomethylated in the clinical phenotype of CRCs.

Aspirin reduces the incidence of metastasis in a pre-clinical study of Braf mutant serrated colorectal neoplasia

BACKGROUND

Aspirin reduces the incidence of conventional adenomas driven by APC mutation and thus colorectal cancer. The effect of aspirin on the ~20% of colorectal cancers arising via BRAF mutation is yet to be established.

METHODS

Braf^V637E/+;Villin-Cre^ERT2/+ mice were allocated to a control (n = 86) or aspirin-supplemented (n = 83) diet. After 14 months the incidence of murine serrated lesions, carcinoma and distant metastases were measured by histological examination. RNA was extracted from carcinomas from each cohort and subjected to sequencing to identify differentially expressed genes and molecular pathways.

RESULTS

Aspirin did not reduce the incidence of murine serrated lesions or carcinoma when compared to control, however, did significantly reduce lesion size (P = 0.0042). Among the mice with carcinoma there was a significant reduction in the incidence of distant metastasis with aspirin treatment (RR 0.69, 95% CI 0.48-0.90, P = 0.0134). Key pathways underlying metastasis of carcinoma cells include NOTCH, FGFR and PI3K signalling, were significantly downregulated in carcinomas sampled from mice on an aspirin-supplemented diet.

CONCLUSIONS

Aspirin reduces the incidence of metastatic Braf mutant carcinoma, although this is not due to a reduction in primary disease. The reduction in metastasis could be attributed to a delay or prevention of molecular changes within the primary site driving metastatic growth.

APC Mutation Marks an Aggressive Subtype of BRAF Mutant Colorectal Cancers

BACKGROUND

WNT activation is a hallmark of colorectal cancer. BRAF mutation is present in 15% of colorectal cancers, and the role of mutations in WNT signaling regulators in this context is unclear. Here, we evaluate the mutational landscape of WNT signaling regulators in BRAF mutant cancers.

METHODS

we performed exome-sequencing on 24 BRAF mutant colorectal cancers and analyzed these data in combination with 175 publicly available BRAF mutant colorectal cancer exomes. We assessed the somatic mutational landscape of WNT signaling regulators, and performed hotspot and driver mutation analyses to identify potential drivers of WNT signaling. The effects of Apc and Braf mutation were modelled, in vivo, using the Apcmin/+ and BrafV637/Villin-CreERT2/+ mouse, respectively.

RESULTS

RNF43 was the most frequently mutated WNT signaling regulator (41%). Mutations in the beta-catenin destruction complex occurred in 48% of cancers. Hotspot analyses identified potential cancer driver genes in the WNT signaling cascade, including MEN1, GNG12 and WNT16. Truncating APC mutation was identified in 20.8% of cancers. Truncating APC mutation was associated with early age at diagnosis (p < 2 × 10-5), advanced stage (p < 0.01), and poor survival (p = 0.026). Apcmin/+/BrafV637 animals had more numerous and larger SI and colonic lesions (p < 0.0001 and p < 0.05, respectively), and a markedly reduced survival (median survival: 3.2 months, p = 8.8 × 10-21), compared to animals with Apc or Braf mutation alone.

CONCLUSIONS

the WNT signaling axis is frequently mutated in BRAF mutant colorectal cancers. WNT16 and MEN1 may be novel drivers of aberrant WNT signaling in colorectal cancer. Co-mutation of BRAF and APC generates an extremely aggressive neoplastic phenotype that is associated with poor patient outcome.

Alterations in signaling pathways that accompany spontaneous transition to malignancy in a mouse model of BRAF mutant microsatellite stable colorectal cancer

The serrated neoplasia pathway gives rise to a distinct subgroup of colorectal cancers distinguished by the presence of mutant BRAF^V600E and the CpG Island Methylator Phenotype (CIMP). BRAF mutant CRC are commonly associated with microsatellite instability, which have an excellent clinical outcome. However, a proportion of BRAF mutant CRC retain microsatellite stability and have a dismal prognosis. The molecular drivers responsible for the development of this cancer subgroup are unknown. To address this, we established a murine model of BRAF^V600E mutant microsatellite stable CRC and comprehensively investigated the exome and transcriptome to identify molecular alterations in signaling pathways that drive malignancy. Exome sequencing of murine serrated lesions (mSL) and carcinomas identified frequent hot spot mutations within the gene encoding β-catenin (Ctnnb1). Immunohistochemical staining of β-catenin indicated that these mutations led to an increase in the presence of aberrant nuclear β-catenin that resulted in gene expression changes in targets of β-catenin transcription. Gene expression profiling identified a significant enrichment for transforming growth factor-β (TGF-β) signaling that was present in mSL and carcinomas. Early activation of TGF-β suggests that this pathway may be an early cue directing mSL to microsatellite stable carcinoma. These findings in the mouse model support the importance of alterations in WNT and TGF-β signaling during the transition of human sessile serrated lesions to malignancy.

DNA methylation changes that precede onset of dysplasia in advanced sessile serrated adenomas

Background

Sessile serrated adenomas (SSAs) are common polyps which give rise to 20–30% of colorectal cancer (CRC). SSAs display clinicopathologic features which present challenges in surveillance, including overrepresentation in young patients, proclivity for the proximal colon and rarity of histologic dysplasia (referred to then as SSAs with dysplasia, SSADs). Once dysplasia develops, there is rapid progression to CRC, even at a small size. There is therefore a clinical need to separate the “advanced” SSAs at high risk of progression to SSAD and cancer from ordinary SSAs. Since SSAs are known to accumulate methylation over time prior to the development of dysplasia, SSAD backgrounds (the remnant SSA present within an SSAD) likely harbour additional methylation events compared with ordinary SSAs. We therefore performed MethyLight and comprehensive methylation array (Illumina MethylationEPIC) on 40 SSAD backgrounds and 40 matched ordinary SSAs, and compared the methylation results with CRC methylation, CRC expression and immunohistochemical data.

Results

SSAD backgrounds demonstrated significant hypermethylation of CpG islands compared with ordinary SSAs, and the proportion of hypermethylated probes decreased progressively in the shore, shelf and open sea regions. Hypomethylation occurred in concert with hypermethylation, which showed a reverse pattern, increasing progressively away from the island regions. These methylation changes were also identified in BRAF-mutant hypermethylated CRCs. When compared with CRC expression data, SV2B, MLH1/EPM2AIP1, C16orf62, RCOR3, BAIAP3, OGDHL, HDHD3 and ATP1B2 demonstrated both promoter hypermethylation and decreased expression. Although SSAD backgrounds were histologically indistinguishable from ordinary SSAs, MLH1 methylation was detectable via MethyLight in 62.9% of SSAD backgrounds, and focal immunohistochemical MLH1 loss was seen in 52.5% of SSAD backgrounds.

Conclusions

Significant hyper- and hypomethylation events occur during SSA progression well before the development of histologically identifiable changes. Methylation is a heterogeneous process within individual SSAs, as typified by MLH1, where both MLH1 methylation and focal immunohistochemical MLH1 loss can be seen in the absence of dysplasia. This heterogeneity is likely a generalised phenomenon and should be taken into account in future methylation-based studies and the development of clinical methylation panels.

Electronic supplementary material

The online version of this article (10.1186/s13148-019-0691-4) contains supplementary material, which is available to authorized users.

RNF43 is mutated less frequently in Lynch Syndrome compared with sporadic microsatellite unstable colorectal cancers

The WNT signaling pathway is commonly altered during colorectal cancer development. The E3 ubiquitin ligase, RNF43, negatively regulates the WNT signal through increased ubiquitination and subsequent degradation of the Frizzled receptor. RNF43 has recently been reported to harbor frequent truncating frameshift mutations in sporadic microsatellite unstable (MSI) colorectal cancers. This study assesses the relative frequency of RNF43 mutations in hereditary colorectal cancers arising in the setting of Lynch syndrome. The entire coding region of RNF43 was Sanger sequenced in 24 colorectal cancers from 23 patients who either (i) carried a germline mutation in one of the DNA mismatch repair genes (MLH1, MSH6, MSH2, PMS2), or (ii) showed immunohistochemical loss of expression of one or more of the DNA mismatch repair proteins, was BRAF wild type at V600E, were under 60 years of age at diagnosis, and demonstrated no promoter region methylation for MLH1 in tumor DNA. A validation cohort of 44 colorectal cancers from mismatch repair germline mutation carriers from the Australasian Colorectal Cancer Family Registry (ACCFR) were sequenced for the most common truncating mutation hotspots (X117 and X659). RNF43 mutations were found in 9 of 24 (37.5%) Lynch syndrome colorectal cancers. The majority of mutations were frameshift deletions in the G659 G7 repeat tract (29%); 2 cancers (2/24, 8%) from the one patient harbored frameshift mutations at codon R117 (C6 repeat tract) within exon 3. In the ACCFR validation cohort, RNF43 hotspot mutations were identified in 19/44 (43.2%) of samples, which was not significantly different to the initial series. The proportion of mutant RNF43 in Lynch syndrome related colorectal cancers is significantly lower than the previously reported mutation rate found in sporadic MSI colorectal cancers. These findings identify further genetic differences between sporadic and hereditary colorectal cancers. This may be because Lynch Syndrome cancers commonly arise in colorectal adenomas already bearing the APC mutation, whereas sporadic microsatellite unstable colorectal cancers arise from serrated polyps typically lacking APC mutation, decreasing the selection pressure on other WNT signaling related loci in Lynch syndrome.

MLH1-93 G/a polymorphism is associated with MLH1 promoter methylation and protein loss in dysplastic sessile serrated adenomas with BRAFV600E mutation

BACKGROUND

Sessile serrated adenomas with BRAF mutation progress rapidly to cancer following the development of dysplasia (SSAD). Approximately 75% of SSADs methylate the mismatch repair gene MLH1, develop mismatch repair deficiency and the resultant cancers have a good prognosis. The remaining SSADs and BRAF mutant traditional serrated adenomas (TSA) develop into microsatellite stable cancers with a poor prognosis. The reason for this dichotomy is unknown. In this study, we assessed the genotypic frequency of the MLH1-93 polymorphism rs1800734 in SSADs and TSAs to determine if the uncommon variant A allele predisposes to MLH1 promoter hypermethylation.

METHODS

We performed genotyping for the MLH1-93 polymorphism, quantitative methylation specific PCR, and MLH1 immunohistochemistry on 124 SSAD, 128 TSA, 203 BRAF mutant CRCs and 147 control subjects with normal colonoscopy.

RESULTS

The minor A allele was significantly associated with a dose dependent increase in methylation at the MLH1 promoter in SSADs (p = 0.022). The AA genotype was only observed in SSADs with MLH1 loss. The A allele was also overrepresented in BRAF mutant cancers with MLH1 loss. Only one of the TSAs showed loss of MLH1 and the overall genotype distribution in TSAs did not differ from controls.

CONCLUSIONS

The MLH1-93 AA genotype is significantly associated with promoter hypermethylation and MLH1 loss in the context of SSADs. BRAF mutant microsatellite stable colorectal cancers with the AA genotype most likely arise in TSAs since the A allele does not predispose to methylation in this context.