Manipulating base quality scores enables variant calling from bisulfite sequencing alignments using conventional bayesian approaches

DNA Methylation Carries Signatures of Sublethal Effects Under Thermal Stress in Loggerhead Sea Turtles

To date, studies of the impacts of climate warming on individuals and populations have mostly focused on mortality and thermal tolerance. In contrast, much less is known about the consequences of sublethal effects, which are more challenging to detect, particularly in wild species with cryptic life histories. This necessitates the development of molecular tools to identify their signatures. In a split-clutch field experiment, we relocated clutches of wild, nesting loggerhead sea turtles (Caretta caretta) to an in situ hatchery. Eggs were then split into two sub-clutches and incubated under shallow or deep conditions, with those in the shallow treatment experiencing significantly higher temperatures in otherwise natural conditions. Although no difference in hatching success was observed between treatments, hatchlings from the shallow, warmer treatment had different length-mass relationships and were weaker at locomotion tests than their siblings incubated in the deep, cooler treatment. To characterise the molecular signatures of these thermal effects, we performed whole genome bisulfite sequencing on blood samples collected upon emergence. We identified 287 differentially methylated sites between hatchlings from different treatments, including on genes with neurodevelopmental, cytoskeletal, and lipid metabolism functions. Taken together, our results show that higher incubation temperatures induce sublethal effects in hatchlings, which are reflected in their DNA methylation status at identified sites. These sites could be used as biomarkers of thermal stress, especially if they are retained across life stages. Overall, this study suggests that global warming reduces hatchling fitness, which has implications for dispersal capacity and ultimately a population's adaptive potential. Conservation efforts for these endangered species and similar climate-threatened taxa will therefore benefit from strategies for monitoring and mitigating exposure to temperatures that induce sublethal effects.

Beyond the base pairs: comparative genome-wide DNA methylation profiling across sequencing technologies

Deoxyribonucleic acid (DNA) methylation plays a key role in gene regulation and is critical for development and human disease. Techniques such as whole-genome bisulfite sequencing (WGBS) and reduced representation bisulfite sequencing (RRBS) allow DNA methylation analysis at the genome scale, with Illumina NovaSeq 6000 and MGI Tech DNBSEQ-T7 being popular due to their efficiency and affordability. However, detailed comparative studies of their performance are not available. In this study, we constructed 60 WGBS and RRBS libraries for two platforms using different types of clinical samples and generated approximately 2.8 terabases of sequencing data. We systematically compared quality control metrics, genomic coverage, CpG methylation levels, intra- and interplatform correlations, and performance in detecting differentially methylated positions. Our results revealed that the DNBSEQ platform exhibited better raw read quality, although base quality recalibration indicated potential overestimation of base quality. The DNBSEQ platform also showed lower sequencing depth and less coverage uniformity in GC-rich regions than did the NovaSeq platform and tended to enrich methylated regions. Overall, both platforms demonstrated robust intra- and interplatform reproducibility for RRBS and WGBS, with NovaSeq performing better for WGBS, highlighting the importance of considering these factors when selecting a platform for bisulfite sequencing.

Methylation marks in blood DNA reveal breast cancer risk in patients fulfilling hereditary disease criteria

Less than 15-20% of patients who meet the criteria for hereditary breast and ovarian cancer (HBOC) carry pathogenic coding genetic mutations, implying that other molecular mechanisms may contribute to the increased risk of this condition. DNA methylation in peripheral blood has been suggested as a potential epigenetic marker for the risk of breast cancer (BC). We aimed to discover methylation marks in peripheral blood associated with BC in 231 pre-treatment BC patients meeting HBOC criteria, testing negative for coding pathogenic variants, and 156 healthy controls, through methylation analysis by targeted bisulfite sequencing on 18 tumor suppressor gene promoters (330 CpG sites). We found i) hypermethylation in EPCAM (17 CpG sites; p = 0.017) and RAD51C (27 CpG sites; p = 0.048); ii) hypermethylation in 36 CpG-specific sites (FDR q < 0.05) in the BC patients; iii) four specific CpG sites were associated with a higher risk of BC (FDR q < 0.01, Bonferroni p < 0.001): cg89786999-FANCI (OR = 1.65; 95% CI:1.2-2.2), cg23652916-PALB2 (OR = 2.83; 95% CI:1.7-4.7), cg47630224-MSH2 (OR = 4.17; 95% CI:2.1-8.5), and cg47596828-EPCAM (OR = 1.84; 95% CI:1.5-2.3). Validation of cg47630224-MSH2 methylation in one Australian cohort showed an association with 3-fold increased BC risk (AUC: 0.929; 95% CI: 0.904-0.955). Our findings suggest that four DNA methylation CpG sites may be associated with a higher risk of BC, potentially serving as biomarkers in patients without detectable coding mutations.

Mixed Patterns of Intergenerational DNA Methylation Inheritance in Acropora

For sessile organisms at high risk from climate change, phenotypic plasticity can be critical to rapid acclimation. Epigenetic markers like DNA methylation are hypothesized as mediators of plasticity; methylation is associated with the regulation of gene expression, can change in response to ecological cues, and is a proposed basis for the inheritance of acquired traits. Within reef-building corals, gene-body methylation (gbM) can change in response to ecological stressors. If coral DNA methylation is transmissible across generations, this could potentially facilitate rapid acclimation to environmental change. We investigated methylation heritability in Acropora, a stony reef-building coral. Two Acropora millepora and two Acropora selago adults were crossed, producing eight offspring crosses (four hybrid, two of each species). We used whole-genome bisulfite sequencing to identify methylated loci and allele-specific alignments to quantify per-locus inheritance. If methylation is heritable, differential methylation (DM) between the parents should equal DM between paired offspring alleles at a given locus. We found a mixture of heritable and nonheritable loci, with heritable portions ranging from 44% to 90% among crosses. gBM was more heritable than intergenic methylation, and most loci had a consistent degree of heritability between crosses (i.e. the deviation between parental and offspring DM were of similar magnitude and direction). Our results provide evidence that coral methylation can be inherited but that heritability is heterogenous throughout the genome. Future investigations into this heterogeneity and its phenotypic implications will be important to understanding the potential capability of intergenerational environmental acclimation in reef building corals.

DNA methylation in the wild: epigenetic transgenerational inheritance can mediate adaptation in clones of wild strawberry (Fragaria vesca)

Due to the accelerating climate change, it is crucial to understand how plants adapt to rapid environmental changes. Such adaptation may be mediated by epigenetic mechanisms like DNA methylation, which could heritably alter phenotypes without changing the DNA sequence, especially across clonal generations. However, we are still missing robust evidence of the adaptive potential of DNA methylation in wild clonal populations. Here, we studied genetic, epigenetic and transcriptomic variation of Fragaria vesca, a predominantly clonally reproducing herb. We examined samples from 21 natural populations across three climatically distinct geographic regions, as well as clones of the same individuals grown in a common garden. We found that epigenetic variation was partly associated with climate of origin, particularly in non-CG contexts. Importantly, a large proportion of this variation was heritable across clonal generations. Additionally, a subset of these epigenetic changes affected the expression of genes mainly involved in plant growth and responses to pathogen and abiotic stress. These findings highlight the potential influence of epigenetic changes on phenotypic traits. Our findings indicate that variation in DNA methylation, which can be environmentally inducible and heritable, may enable clonal plant populations to adjust to their environmental conditions even in the absence of genetic adaptation.

Comprehensive DNA Methylation Analysis Indicates That Pancreatic Intraepithelial Neoplasia Lesions Are Acinar-Derived and Epigenetically Primed for Carcinogenesis

Pancreatic ductal adenocarcinoma (PDAC) is believed to arise from the accumulation of a series of somatic mutations and is also frequently associated with pancreatic intraepithelial neoplasia (PanIN) lesions. However, there is still debate as to whether the cell type-of-origin of PanINs and PDACs in humans is acinar or ductal. As cell type identity is maintained epigenetically, DNA methylation changes during pancreatic neoplasia can provide a compelling perspective to examine this question. Here, we performed laser-capture microdissection on surgically resected specimens from 18 patients to isolate, with high purity, DNA for whole-genome bisulfite sequencing from four relevant cell types: acini, nonneoplastic ducts, PanIN lesions, and PDAC lesions. Differentially methylated regions (DMR) were identified using two complementary analytical approaches: bsseq, which identifies any DMRs but is particularly useful for large block-like DMRs, and informME, which profiles the potential energy landscape across the genome and is particularly useful for identifying differential methylation entropy. Both global methylation profiles and block DMRs clearly implicated an acinar origin for PanINs. At the gene level, PanIN lesions exhibited an intermediate acinar-ductal phenotype resembling acinar-to-ductal metaplasia. In 97.6% of PanIN-specific DMRs, PanIN lesions had an intermediate methylation level between normal and PDAC, which suggests from an information theory perspective that PanIN lesions are epigenetically primed to progress to PDAC. Thus, epigenomic analysis complements histopathology to define molecular progression toward PDAC. The shared epigenetic lineage between PanIN and PDAC lesions could provide an opportunity for prevention by targeting aberrantly methylated progression-related genes.

SIGNIFICANCE

Analysis of DNA methylation landscapes provides insights into the cell-of-origin of PanIN lesions, clarifies the role of PanIN lesions as metaplastic precursors to human PDAC, and suggests potential targets for chemoprevention.

On the causes of gene-body methylation variation in Arabidopsis thaliana

Gene-body methylation (gbM) refers to sparse CG methylation of coding regions, which is especially prominent in evolutionarily conserved house-keeping genes. It is found in both plants and animals, but is directly and stably (epigenetically) inherited over multiple generations in the former. Studies in Arabidopsis thaliana have demonstrated that plants originating from different parts of the world exhibit genome-wide differences in gbM, which could reflect direct selection on gbM, but which could also reflect an epigenetic memory of ancestral genetic and/or environmental factors. Here we look for evidence of such factors in F2 plants resulting from a cross between a southern Swedish line with low gbM and a northern Swedish line with high gbM, grown at two different temperatures. Using bisulfite-sequencing data with nucleotide-level resolution on hundreds of individuals, we confirm that CG sites are either methylated (nearly 100% methylation across sampled cells) or unmethylated (approximately 0% methylation across sampled cells), and show that the higher level of gbM in the northern line is due to more sites being methylated. Furthermore, methylation variants almost always show Mendelian segregation, consistent with their being directly and stably inherited through meiosis. To explore how the differences between the parental lines could have arisen, we focused on somatic deviations from the inherited state, distinguishing between gains (relative to the inherited 0% methylation) and losses (relative to the inherited 100% methylation) at each site in the F2 generation. We demonstrate that deviations predominantly affect sites that differ between the parental lines, consistent with these sites being more mutable. Gains and losses behave very differently in terms of the genomic distribution, and are influenced by the local chromatin state. We find clear evidence for different trans-acting genetic polymorphism affecting gains and losses, with those affecting gains showing strong environmental interactions (G×E). Direct effects of the environment were minimal. In conclusion, we show that genetic and environmental factors can change gbM at a cellular level, and hypothesize that these factors can also lead to transgenerational differences between individuals via the inclusion of such changes in the zygote. If true, this could explain genographic pattern of gbM with selection, and would cast doubt on estimates of epimutation rates from inbred lines in constant environments.