DNA methylation signatures of early-life adversity are exposure-dependent in wild baboons

DNA Methylation Dynamics Reflect Sex and Status Differences in Mortality Rates in a Polygynous Bat

Males of polygynous mammals often do not live as long as females and, in some cases, exhibit evidence of earlier senescence. Patterns of DNA methylation (DNAm) have recently been used to predict chronological age in mammals. Whether DNAm also changes as a consequence of survival and senescence is largely untested in wild animals. In this study, we estimate mortality rates using recaptures of 2700 greater spear-nosed bats, Phyllostomus hastatus, over 34 years and DNAm profiled for over 300 adult bats. In this species, one male typically controls mating access to a group of unrelated females. Bayesian analysis reveals that mortality risk in males is 1.8 times that of females, and comparison of age-associated differences in DNAm indicates that DNAm changes 1.4 times faster in males than females. Therefore, even though the age of either sex is predicted by a common set of sites, the methylome of males is more dynamic than that of females. Sites associated with sex differences in the rate of DNAm change are sensitive to androgens and enriched on the X chromosome. Sites that exhibit hypermethylation are enriched in promoters of genes involved in the regulation of metabolic processes. Unexpectedly, subordinate males have higher mortality rates than reproductively dominant males and exhibit faster DNAm change than dominants at dozens of sites. Our results reveal that differences in mortality associated with sex and social status are reflected by changes in DNA methylation, providing novel insights into mechanisms of aging and mortality in this and likely other wild animal populations.

Addressing missing context in regulatory variation across primate evolution

In primates, loci associated with adaptive trait variation often fall in non-coding regions. Understanding the mechanisms linking these regulatory variants to fitness-relevant phenotypes remains challenging, but can be addressed using functional genomic data. However, such data are rarely generated at scale in non-human primates. When they are, only select tissues, cell types, developmental stages, and cellular environments are typically considered, despite appreciation that adaptive variants often exhibit context-dependent effects. In this review, we 1) discuss why context-dependent regulatory loci might be especially evolutionarily relevant in primates, 2) explore challenges and emerging solutions for mapping such context-dependent variation, and 3) discuss the scientific questions these data could address. We argue that filling this gap will provide critical insights into evolutionary processes, human disease, and regulatory adaptation.

The interplay between epigenomic and transcriptomic variation during ecotype divergence in stickleback

BACKGROUND

Populations colonizing contrasting environments are likely to undergo adaptive divergence and evolve ecotypes with locally adapted phenotypes. While diverse molecular mechanisms underlying ecotype divergence have been identified, less is known about their interplay and degree of divergence.

RESULTS

Here we integrated epigenomic and transcriptomic data to explore the interactions among gene expression, alternative splicing, DNA methylation, and microRNA expression to gauge the extent to which patterns of divergence at the four molecular levels are aligned in a case of postglacial divergence between marine and freshwater ecotypes of nine-spined sticklebacks (Pungitius pungitius). Despite significant genome-wide associations between epigenomic and transcriptomic variation, we found largely non-parallel patterns of ecotype divergence across epigenomic and transcriptomic levels, with predominantly nonoverlapping (ranging from 43.40 to 87.98%) sets of differentially expressed, spliced and methylated genes, and candidate genes targeted by differentially expressed miRNA between the ecotypes. Furthermore, we found significant variation in the extent of ecotype divergence across different molecular mechanisms, with differential methylation and differential splicing showing the highest and lowest extent of divergence between ecotypes, respectively. Finally, we found a significant enrichment of genes associated with ecotype divergence in differential methylation.

CONCLUSIONS

Our results suggest a nuanced relationship between epigenomic and transcriptomic processes, with alignment at the genome-wide level masking relatively independent effects of different molecular mechanisms on ecotype divergence at the gene level.

Embryonic temperature influences transcriptomic and methylation profiles in the liver of juvenile largemouth bass

Understanding the impacts of environmental conditions at early life stages on phenotypes and physiological responses to thermal variability at later stages is crucial for elucidating adaptive strategies in fish species. This study investigated the lasting effects of embryonic temperature on the growth performance, transcriptomic profiles, and CpG methylation status of juvenile largemouth bass (Micropterus salmoides) under normal and heat stress (HS) conditions. Embryos were incubated at three temperatures (22 °C, 25 °C, and 28 °C), reared at a constant 25 °C for three months, and subjected to acute HS at 37 °C. Liver samples were collected before and after HS for mRNA sequencing and reduced representation bisulfite sequencing. Significant differences in body size, body weight, condition factor, and hepatosomatic index were observed among groups. Fish hatched at 28 °C displayed a significantly higher standard length and body weight and lower hepatosomatic index than those hatched at lower temperatures. PCA analysis and Venn diagrams based on transcriptomes revealed transcriptomic response to HS differed at 28 °C while 22 and 25 °C were similar. Heat shock protein genes followed a similar trend. Epigenetic analyses revealed distinct CpG methylation patterns across incubation groups, while DNA methylation barely contributes to transcriptional differences. Under HS, different incubation groups exhibit various DNA methylation alterations. The "Neuroactive ligand-receptor interaction" pathway appeared to play an important role in the response to HS, suggesting a potential involvement of epigenetic regulation. Additionally, the atrnl1 gene may be involved in a DNA methylation-mediated regulatory mechanism in response to HS.

Baboons, bacteria, and biological clocks address an age-old question

Studying the fecal microbiota of wild baboons helps provide new insight into the factors that influence biological aging.

Cost-effective solutions for high-throughput enzymatic DNA methylation sequencing

Characterizing DNA methylation patterns is important for addressing key questions in evolutionary biology, geroscience, and medical genomics. While costs are decreasing, whole-genome DNA methylation profiling remains prohibitively expensive for most population-scale studies, creating a need for cost-effective, reduced representation approaches (i.e., assays that rely on microarrays, enzyme digests, or sequence capture to target a subset of the genome). Most common whole genome and reduced representation techniques rely on bisulfite conversion, which can damage DNA resulting in DNA loss and sequencing biases. Enzymatic methyl sequencing (EM-seq) was recently proposed to overcome these issues, but thorough benchmarking of EM-seq combined with cost-effective, reduced representation strategies has not yet been performed. To do so, we optimized Targeted Methylation Sequencing protocol (TMS)-which profiles ∼4 million CpG sites-for miniaturization, flexibility, and multispecies use at a cost of ∼$80. First, we tested modifications to increase throughput and reduce cost, including increasing multiplexing, decreasing DNA input, and using enzymatic rather than mechanical fragmentation to prepare DNA. Second, we compared our optimized TMS protocol to commonly used techniques, specifically the Infinium MethylationEPIC BeadChip (n=55 paired samples) and whole genome bisulfite sequencing (n=6 paired samples). In both cases, we found strong agreement between technologies (R² = 0.97 and 0.99, respectively). Third, we tested the optimized TMS protocol in three non-human primate species (rhesus macaques, geladas, and capuchins). We captured a high percentage (mean=77.1%) of targeted CpG sites and produced methylation level estimates that agreed with those generated from reduced representation bisulfite sequencing (R² = 0.98). Finally, we applied our protocol to profile age-associated DNA methylation variation in two subsistence-level populations-the Tsimane of lowland Bolivia and the Orang Asli of Peninsular Malaysia-and found age-methylation patterns that were strikingly similar to those reported in high income cohorts, despite known differences in age-health relationships between lifestyle contexts. Altogether, our optimized TMS protocol will enable cost-effective, population-scale studies of genome-wide DNA methylation levels across human and non-human primate species.

Uncovering methylation-dependent genetic effects on regulatory element function in diverse genomes

A major goal in evolutionary biology and biomedicine is to understand the complex interactions between genetic variants, the epigenome, and gene expression. However, the causal relationships between these factors remain poorly understood. mSTARR-seq, a methylation-sensitive massively parallel reporter assay, is capable of identifying methylation-dependent regulatory activity at many thousands of genomic regions simultaneously, and allows for the testing of causal relationships between DNA methylation and gene expression on a region-by-region basis. Here, we developed a multiplexed mSTARR-seq protocol to assay naturally occurring human genetic variation from 25 individuals sampled from 10 localities in Europe and Africa. We identified 6,957 regulatory elements in either the unmethylated or methylated state, and this set was enriched for enhancer and promoter annotations, as expected. The expression of 58% of these regulatory elements was modulated by methylation, which was generally associated with decreased RNA expression. Within our set of regulatory elements, we used allele-specific expression analyses to identify 8,020 sites with genetic effects on gene regulation; further, we found that 42.3% of these genetic effects varied between methylated and unmethylated states. Sites exhibiting methylation-dependent genetic effects were enriched for GWAS and EWAS annotations, implicating them in human disease. Compared to datasets that assay DNA from a single European individual, our multiplexed assay uncovers dramatically more genetic effects and methylation-dependent genetic effects, highlighting the importance of including diverse individuals in assays which aim to understand gene regulatory processes.

DNA Methylation of PXDN Is Associated with Early-Life Adversity in Adult Mental Disorders

Early-life adversity (ELA) is characterized by exposure to traumatic events during early periods of life, particularly involving emotional, sexual and/or physical adversities during childhood. Mental disorders are strongly influenced by environmental and lifestyle-related risk factors including ELA. However, the molecular link between ELA and the risk of an adult mental disorder is still not fully understood. Evidence is emerging that long-lasting changes in the epigenetic processes regulating gene expression, such as DNA methylation, play an important role in the biological mechanisms linking ELA and mental disorders. Based on a recent study, we analyzed the DNA methylation of a specific CpG site within the gene PXDN-cg10888111-in blood in the context of ELA across a set of psychiatric disorders, namely Borderline Personality Disorder (BPD), Major Depressive Disorder (MDD) and Social Anxiety Disorder (SAD), and its potential contribution to their pathogenesis. We found significant hypermethylation in mentally ill patients with high levels of ELA compared to patients with low levels of ELA, whereas cg10888111 methylation in healthy control individuals was not affected by ELA. Further investigations revealed that this effect was driven by the MDD cohort. Providing a direct comparison of cg10888111 DNA methylation in blood in the context of ELA across three mental disorders, our results indicate the role of PXDN regulation in the response to ELA in the pathogenesis of mental disorders, especially MDD. Further studies will be needed to validate these results and decipher the corresponding biological network that is involved in the transmission of ELA to an adult mental disorder in general.