Profiling DNA methylation differences between inbred mouse strains on the Illumina Human Infinium MethylationEPIC microarray

Integrative ensemble modelling of cetuximab sensitivity in colorectal cancer patient-derived xenografts

Patient-derived xenografts (PDXs) are tumour fragments engrafted into mice for preclinical studies. PDXs offer clear advantages over simpler in vitro cancer models - such as cancer cell lines (CCLs) and organoids - in terms of structural complexity, heterogeneity, and stromal interactions. Here, we characterise 231 colorectal cancer PDXs at the genomic, transcriptomic, and epigenetic levels, along with their response to cetuximab, an EGFR inhibitor used clinically for metastatic colorectal cancer. After evaluating the PDXs' quality, stability, and molecular concordance with publicly available patient cohorts, we present results from training, interpreting, and validating the integrative ensemble classifier CeSta. This model takes in input the PDXs' multi-omic characterisation and predicts their sensitivity to cetuximab treatment, achieving an area under the receiver operating characteristics curve > 0.88. Our study demonstrates that large PDX collections can be leveraged to train accurate, interpretable drug sensitivity models that: (1) better capture patient-derived therapeutic biomarkers compared to models trained on CCL data, (2) can be robustly validated across independent PDX cohorts, and (3) could contribute to the development of future therapeutic biomarkers.

Pleiotropic influence of DNA methylation QTLs on physiological and ageing traits

DNA methylation is influenced by genetic and non-genetic factors. Here, we chart quantitative trait loci (QTLs) that modulate levels of methylation at highly conserved CpGs using liver methylome data from mouse strains belonging to the BXD family. A regulatory hotspot on chromosome 5 had the highest density of trans-acting methylation QTLs (trans-meQTLs) associated with multiple distant CpGs. We refer to this locus as meQTL.5a. Trans-modulated CpGs showed age-dependent changes and were enriched in developmental genes, including several members of the MODY pathway (maturity onset diabetes of the young). The joint modulation by genotype and ageing resulted in a more 'aged methylome' for BXD strains that inherited the DBA/2J parental allele at meQTL.5a. Further, several gene expression traits, body weight, and lipid levels mapped to meQTL.5a, and there was a modest linkage with lifespan. DNA binding motif and protein-protein interaction enrichment analyses identified the hepatic nuclear factor, Hnf1a (MODY3 gene in humans), as a strong candidate. The pleiotropic effects of meQTL.5a could contribute to variations in body size and metabolic traits, and influence CpG methylation and epigenetic ageing that could have an impact on lifespan.

Lack of Methylation Changes in GJB2 and RB1 Non-coding Regions of Cochlear Implant Patients with Sensorineural Hearing Loss

Objective

Recent advances in epigenetic studies continue to reveal novel mechanisms of gene regulation and control, however little is known on the role of epigenetics in sensorineural hearing loss (SNHL) in humans. We aimed to investigate the methylation patterns of two regions, one in RB1 and another in GJB2 in Filipino patients with SNHL compared to hearing control individuals.

Methods

We investigated an RB1 promoter region that was previously identified as differentially methylated in children with SNHL and lead exposure. Additionally, we investigated a sequence in an enhancer-like region within GJB2 that contains four CpGs in close proximity. Bisulfite conversion was performed on salivary DNA samples from 15 children with SNHL and 45 unrelated ethnically-matched individuals. We then performed methylation-specific real-time PCR analysis (qMSP) using TaqMan® probes to determine percentage methylation of the two regions.

Results

Using qMSP, both our cases and controls had zero methylation at the targeted GJB2 and RB1 regions.

Conclusion

Our study showed no changes in methylation at the selected CpG regions in RB1 and GJB2 in the two comparison groups with or without SNHL. This may be due to a lack of environmental exposures to these target regions. Other epigenetic marks may be present around these regions as well as those of other HL-associated genes.

Speg interactions that regulate the stability of excitation-contraction coupling protein complexes in triads and dyads

Here we show that striated muscle preferentially expressed protein kinase α (Spegα) maintains cardiac function in hearts with Spegβ deficiency. Speg is required for stability of excitation-contraction coupling (ECC) complexes and interacts with esterase D (Esd), Cardiomyopathy-Associated Protein 5 (Cmya5), and Fibronectin Type III and SPRY Domain Containing 2 (Fsd2) in cardiac and skeletal muscle. Mice with a sequence encoding a V5/HA tag inserted into the first exon of the Speg gene (HA-Speg mice) display a >90% decrease in Spegβ but Spegα is expressed at ~50% of normal levels. Mice deficient in both Spegα and Speg β (Speg KO mice) develop a severe dilated cardiomyopathy and muscle weakness and atrophy, but HA-Speg mice display mild muscle weakness with no cardiac involvement. Spegα in HA-Speg mice suppresses Ca2+ leak, proteolytic cleavage of Jph2, and disruption of transverse tubules. Despite it's low levels, HA-Spegβ immunoprecipitation identified Esd, Cmya5 and Fsd2 as Spegβ binding partners that localize to triads and dyads to stabilize ECC complexes. This study suggests that Spegα and Spegβ display functional redundancy, identifies Esd, Cmya5 and Fsd2 as components of both cardiac dyads and skeletal muscle triads and lays the groundwork for the identification of new therapeutic targets for centronuclear myopathy.

Variation in histone configurations correlates with gene expression across nine inbred strains of mice

The Diversity Outbred (DO) mice and their inbred founders are widely used models of human disease. However, although the genetic diversity of these mice has been well documented, their epigenetic diversity has not. Epigenetic modifications, such as histone modifications and DNA methylation, are important regulators of gene expression and, as such, are a critical mechanistic link between genotype and phenotype. Therefore, creating a map of epigenetic modifications in the DO mice and their founders is an important step toward understanding mechanisms of gene regulation and the link to disease in this widely used resource. To this end, we performed a strain survey of epigenetic modifications in hepatocytes of the DO founders. We surveyed four histone modifications (H3K4me1, H3K4me3, H3K27me3, and H3K27ac), as well as DNA methylation. We used ChromHMM to identify 14 chromatin states, each of which represents a distinct combination of the four histone modifications. We found that the epigenetic landscape is highly variable across the DO founders and is associated with variation in gene expression across strains. We found that epigenetic state imputed into a population of DO mice recapitulated the association with gene expression seen in the founders, suggesting that both histone modifications and DNA methylation are highly heritable mechanisms of gene expression regulation. We illustrate how DO gene expression can be aligned with inbred epigenetic states to identify putative cis-regulatory regions. Finally, we provide a data resource that documents strain-specific variation in the chromatin state and DNA methylation in hepatocytes across nine widely used strains of laboratory mice.

Comparative epigenome analysis using Infinium DNA methylation BeadChips

The arrival of the Infinium DNA methylation BeadChips for mice and other nonhuman mammalian species has outpaced the development of the informatics that supports their use for epigenetics study in model organisms. Here, we present informatics infrastructure and methods to allow easy DNA methylation analysis on multiple species, including domesticated animals and inbred laboratory mice (in SeSAMe version 1.16.0+). First, we developed a data-driven analysis pipeline covering species inference, genome-specific data preprocessing and regression modeling. We targeted genomes of 310 species and 37 inbred mouse strains and showed that genome-specific preprocessing prevents artifacts and yields more accurate measurements than generic pipelines. Second, we uncovered the dynamics of the epigenome evolution in different genomic territories and tissue types through comparative analysis. We identified a catalog of inbred mouse strain-specific methylation differences, some of which are linked to the strains' immune, metabolic and neurological phenotypes. By streamlining DNA methylation array analysis for undesigned genomes, our methods extend epigenome research to broad species contexts.

Differential Methylation Profile in Fragile X Syndrome-Prone Offspring Mice after in Utero Exposure to Lactobacillus Reuteri

Environmental factors such as diet, gut microbiota, and infections have proven to have a significant role in epigenetic modifications. It is known that epigenetic modifications may cause behavioral and neuronal changes observed in neurodevelopmental disabilities, including fragile X syndrome (FXS) and autism (ASD). Probiotics are live microorganisms that provide health benefits when consumed, and in some cases are shown to decrease the chance of developing neurological disorders. Here, we examined the epigenetic outcomes in offspring mice after feeding of a probiotic organism, Lactobacillus reuteri (L. reuteri), to pregnant mother animals. In this study, we tested a cohort of Western diet-fed descendant mice exhibiting a high frequency of behavioral features and lower FMRP protein expression similar to what is observed in FXS in humans (described in a companion manuscript in this same GENES special topic issue). By investigating 17,735 CpG sites spanning the whole mouse genome, we characterized the epigenetic profile in two cohorts of mice descended from mothers treated and non-treated with L. reuteri to determine the effect of prenatal probiotic exposure on the prevention of FXS-like symptoms. We found several genes involved in different neurological pathways being differentially methylated (p ≤ 0.05) between the cohorts. Among the key functions, synaptogenesis, neurogenesis, synaptic modulation, synaptic transmission, reelin signaling pathway, promotion of specification and maturation of neurons, and long-term potentiation were observed. The results of this study are relevant as they could lead to a better understanding of the pathways involved in these disorders, to novel therapeutics approaches, and to the identification of potential biomarkers for early detection of these conditions.

DNA methylation dynamics and dysregulation delineated by high-throughput profiling in the mouse

We have developed a mouse DNA methylation array that contains 296,070 probes representing the diversity of mouse DNA methylation biology. We present a mouse methylation atlas as a rich reference resource of 1,239 DNA samples encompassing distinct tissues, strains, ages, sexes, and pathologies. We describe applications for comparative epigenomics, genomic imprinting, epigenetic inhibitors, patient-derived xenograft assessment, backcross tracing, and epigenetic clocks. We dissect DNA methylation processes associated with differentiation, aging, and tumorigenesis. Notably, we find that tissue-specific methylation signatures localize to binding sites for transcription factors controlling the corresponding tissue development. Age-associated hypermethylation is enriched at regions of Polycomb repression, while hypomethylation is enhanced at regions bound by cohesin complex members. Apc ^Min/+ polyp-associated hypermethylation affects enhancers regulating intestinal differentiation, while hypomethylation targets AP-1 binding sites. This Infinium Mouse Methylation BeadChip (version MM285) is widely accessible to the research community and will accelerate high-sample-throughput studies in this important model organism.

Epigenetic Clocks for Mice Based on Age-Associated Regions That are Conserved Between Mouse Strains and Human

Aging of mice can be tracked by DNA methylation changes at specific sites in the genome. In this study, we used the recently released Infinium Mouse Methylation BeadChip to compare such epigenetic modifications in C57BL/6 (B6) and DBA/2J (DBA) mice. We observed marked differences in age-associated DNA methylation in these commonly used inbred mouse strains, indicating that epigenetic clocks for one strain cannot be simply applied to other strains without further verification. In B6 mice age-associated hypomethylation prevailed with focused hypermethylation at CpG islands, whereas in DBA mice CpG islands revealed rather hypomethylation upon aging. Interestingly, the CpGs with highest age-correlation were still overlapping in B6 and DBA mice and included the genes Hsf4, Prima1, Aspa, and Wnt3a. Notably, Hsf4 and Prima1 were also top candidates in previous studies based on whole genome deep sequencing approaches. Furthermore, Hsf4, Aspa, and Wnt3a revealed highly significant age-associated DNA methylation in the homologous regions in human. Subsequently, we used pyrosequencing of the four relevant regions to establish a targeted epigenetic clock that provided very high correlation with chronological age in independent cohorts of B6 (R2 = 0.98) and DBA (R2 = 0.91). Taken together, the methylome differs extensively between B6 and DBA mice, while prominent age-associated changes are conserved among these strains and even in humans. Our new targeted epigenetic clock with 4 CpGs provides a versatile tool for other researchers analyzing aging in mice.

Analysis and annotation of DNA methylation in two nonhuman primate species using the Infinium Human Methylation 450K and EPIC BeadChips

Aim: Nonhuman primates are essential for research on many human diseases. The Infinium Human Methylation450/EPIC BeadChips are popular tools for the study of the methylation state across the human genome at affordable cost. Methods: We performed a precise evaluation and re-annotation of the BeadChip probes for the analysis of genome-wide DNA methylation patterns in rhesus macaques and African green monkeys through in silico analyses combined with functional validation by pyrosequencing. Results: Up to 165,847 of the 450K and 261,545 probes of the EPIC BeadChip can be reliably used. The annotation files are provided in a format compatible with a variety of standard bioinformatic pipelines. Conclusion: Our study will facilitate high-throughput DNA methylation analyses in Macaca mulatta and Chlorocebus sabaeus.

Genome-wide DNA methylation profiling in human breast tissue by Illumina TruSeq methyl capture EPIC sequencing and infinium methylationEPIC beadchip microarray

A newly-developed platform, the Illumina TruSeq Methyl Capture EPIC library prep (TruSeq EPIC), builds on the content of the Infinium MethylationEPIC Beadchip Microarray (EPIC-array) and leverages the power of next-generation sequencing for targeted bisulphite sequencing. We empirically examined the performance of TruSeq EPIC and EPIC-array in assessing genome-wide DNA methylation in breast tissue samples. TruSeq EPIC provided data with a much higher density in the regions when compared to EPIC-array (~2.74 million CpGs with at least 10X coverage vs ~752 K CpGs, respectively). Approximately 398 K CpGs were common and measured across the two platforms in every sample. Overall, there was high concordance in methylation levels between the two platforms (Pearson correlation r = 0.98, P < 0.0001). However, we observed that TruSeq EPIC measurements provided a wider dynamic range and likely a higher quantitative sensitivity for CpGs that were either hypo- or hyper-methylated (β close to 0 or 1, respectively). In addition, when comparing different breast tissue types TruSeq EPIC identified more differentially methylated CpGs than EPIC-array, not only out of additional sites interrogated by TruSeq EPIC alone, but also out of common sites interrogated by both platforms. Our results suggest that both platforms show high reproducibility and reliability in genome-wide DNA methylation profiling, while TruSeq EPIC had a significant improvement over EPIC-array regarding genomic resolution and coverage. The wider dynamic range and likely higher precision of the estimates by the TruSeq EPIC may lead to the identification of novel differentially methylated markers that are associated with disease risk.

Body weight and high-fat diet are associated with epigenetic aging in female members of the BXD murine family

DNA methylation (DNAm) is shaped by genetic and environmental factors and modulated by aging. Here, we examine interrelations between epigenetic aging, body weight (BW), and life span in 12 isogenic strains from the BXD family of mice that exhibit over twofold variation in longevity. Genome-wide DNAm was assayed in 70 liver specimens from predominantly female cases, 6-25 months old, that were maintained on normal chow or high-fat diet (HFD). We defined subsets of CpG regions associated with age, BW at young adulthood, and strain-by-diet-dependent life span. These age-associated differentially methylated CpG regions (age-DMRs) featured distinct genomic characteristics, with DNAm gains over time occurring in sites such as promoters and exons that have high CpG density and low average methylation. CpG regions associated with BW were enriched in introns, tended to have lower methylation in mice with higher BW, and were inversely correlated with gene expression (i.e., higher mRNA levels in mice with higher BW). CpG regions associated with life span were linked to genes involved in life span modulation, including the telomerase reverse transcriptase gene, Tert, which had both lower methylation and higher expression in long-lived strains. An epigenetic clock defined from age-DMRs revealed accelerated aging in mice belonging to strains with shorter life spans. Both higher BW and the HFD were associated with accelerated epigenetic aging. Our results highlight the age-accelerating effect of heavier BW. Furthermore, we demonstrate that the measure of epigenetic aging derived from age-DMRs can predict genotype and diet-induced differences in life span among female BXD members.

DNA methylation aging clocks: challenges and recommendations

Epigenetic clocks comprise a set of CpG sites whose DNA methylation levels measure subject age. These clocks are acknowledged as a highly accurate molecular correlate of chronological age in humans and other vertebrates. Also, extensive research is aimed at their potential to quantify biological aging rates and test longevity or rejuvenating interventions. Here, we discuss key challenges to understand clock mechanisms and biomarker utility. This requires dissecting the drivers and regulators of age-related changes in single-cell, tissue- and disease-specific models, as well as exploring other epigenomic marks, longitudinal and diverse population studies, and non-human models. We also highlight important ethical issues in forensic age determination and predicting the trajectory of biological aging in an individual.

Epigenome-wide association studies in asthma: A systematic review

OBJECTIVE

Asthma is a common chronic respiratory airway disease influenced by environmental factors and possibly their interaction with the human genome causing epigenetic changes. Epigenome-wide association studies (EWAS) have mainly investigated DNA methylation and its association with disease or traits, exposure factors or gene expression. This systematic review aimed to identify all EWAS assessing differentially methylated sites associated with asthma in humans.

DESIGN

Structured systematic literature search following PRISMA guidelines, Newcastle-Ottawa Scale (NOS) for cohort studies was used for bias assessment.

DATA SOURCES

We searched PubMed and Embase databases from 2005 to 2019.

ELIGIBILITY CRITERIA

Epigenome-wide association studies testing association between differential methylation and asthma in humans.

RESULTS

Overall, we identified 16 EWAS studies complying with our search criteria. Twelve studies were conducted on children, and 10 were conducted on sample sizes <150 subjects. Four hundred and nineteen CpGs were reported in children studies after correction for multiple testing. In the adult studies, thousands of differentially methylated sites were identified. Differential methylation in inflammatory-related genes correlated with higher levels of gene expressions of inflammatory modulators in asthma. Differentially methylated genes associated with asthma included SMAD3, SERPINC1, PROK1, IL13, RUNX3 and TIGIT. Forty-one CpGs were replicated at least once in blood samples, and 28 CpGs were replicated in nasal samples.

CONCLUSION

Although many differentially methylated CpGs in genes known to be involved in asthma have been identified in EWAS to date, we conclude that further studies of larger sample sizes and analyses of differential methylation between different phenotypes are needed in order to comprehensively evaluate the role of epigenetic factors in the pathophysiology and heterogeneity of asthma, and the potential clinical utility to predict or classify patients with asthma.

Methylation at cg05575921 of a smoking-related gene (AHRR) in non-smoking Taiwanese adults residing in areas with different PM2.5 concentrations

BACKGROUND

DNA methylation is associated with cancer, metabolic, neurological, and autoimmune disorders. Hypomethylation of aryl hydrocarbon receptor repressor (AHRR) especially at cg05575921 is associated with smoking and lung cancer. Studies on the association between AHRR methylation at cg05575921 and sources of polycyclic aromatic hydrocarbon (PAH) other than smoking are limited. The aim of our study was to assess the pattern of blood DNA methylation at cg05575921 in non-smoking Taiwanese adults living in areas with different PM_2.5 levels.

METHODS

Data on blood DNA methylation, smoking, and residence were retrieved from the Taiwan Biobank dataset (2008-2015). Current and former smokers, as well as individuals with incomplete information were excluded from the current study. The final analysis included 708 participants (279 men and 429 women) aged 30-70 years. PM_2.5 levels have been shown to increase as one moves from the northern through central towards southern Taiwan. Based on this trend, the study areas were categorized into northern, north-central, central, and southern regions.

RESULTS

Living in PM_2.5 areas was associated with lower methylation levels: compared with the northern area (reference area), living in north-central, central, and southern areas was associated with lower methylation levels at cg05575921. However, only methylation levels in those living in central and southern areas were significant (β = - 0.01003, P = 0.009 and β = - 0.01480, P < 0.001, respectively. Even though methylation levels in those living in the north-central area were not statistically significant, the test for linear trend was significant (P < 0.001). When PM_2.5 was included in the regression model, a unit increase in PM_2.5 was associated with 0.00115 (P < 0.001) lower cg05575921 methylation levels.

CONCLUSION

Living in PM_2.5 areas was inversely associated with blood AHRR methylation levels at cg05575921. The methylation levels were lowest in participants residing in southern followed by central and north-central areas. Moreover, when PM_2.5 was included in the regression model, it was inversely associated with methylation levels at cg05575921. Blood methylation at cg05575921 (AHRR) in non-smokers might indicate different exposures to PM_2.5 and lung cancer which is a PM_2.5-related disease.

SOX2 promoter hypermethylation in non-smoking Taiwanese adults residing in air pollution areas

BACKGROUND

Both SOX2 promoter methylation and air pollution have been associated with lung cancer risk. However, little has been done to assess SOX2 promoter methylation in individuals living in air pollution areas. The aim of this study was to investigate SOX2 promoter methylation in non-smoking Taiwanese adults living in areas with different levels of air pollution especially particulate matter with diameter < 2.5 μm (PM2.5).

METHODS

A total of 1142 individuals aged 30-70 years were recruited. Data on SOX2 methylation, residence, age, and exposure to second-hand smoke (SHS) among others were extracted from the Taiwan Biobank dataset (2008-2015). After excluding former and current smokers, alongside those with incomplete information, a total of 461 non-smokers comprising 176 men and 285 women were included in the study. Participants' residences were grouped under northern and central/southern areas because air pollution (PM2.5) is lower in northern compared to central and southern areas.

RESULTS

The methylation levels in men (0.16310 ± 0.01230) and women (0.15740 ± 0.01240) were significantly different (P < .0001). In both sexes, the SOX2 promoter region was shown to be significantly hypermethylated in central and southern areas compared with the northern areas. The regression coefficient (β) was 0.00331 (P = 0.0257) in men and 0.00514 (P < .0001) in women.

CONCLUSION

SOX2 was significantly hypermethylated in both men and women residing in central and southern areas. The consistency in the results for both sexes shows that SOX2 promoter methylation could serve as a potential biomarker for industrial air pollution exposure. Moreover, it might reflect predisposition to cancer. Hence, healthy non-smokers at precancerous stages who have not been clinically diagnosed could be identified.