Epigenetics Research in Evolutionary Biology: Perspectives on Timescales and Mechanisms.

Epigenetics research in evolutionary biology encompasses a variety of research areas, from regulation of gene expression to inheritance of environmentally mediated phenotypes. Such divergent research foci can occasionally render the umbrella term "epigenetics" ambiguous. Here I discuss several areas of contemporary epigenetics research in the context of evolutionary biology, aiming to provide balanced views across timescales and molecular mechanisms. The importance of epigenetics in development is now being assessed in many nonmodel species. These studies not only confirm the importance of epigenetic marks in developmental processes, but also highlight the significant diversity in epigenetic regulatory mechanisms across taxa. Further, these comparative epigenomic studies have begun to show promise toward enhancing our understanding of how regulatory programs evolve. A key property of epigenetic marks is that they can be inherited along mitotic cell lineages, and epigenetic differences that occur during early development can have lasting consequences on the organismal phenotypes. Thus, epigenetic marks may play roles in short-term (within an organism's lifetime or to the next generation) adaptation and phenotypic plasticity. However, the extent to which observed epigenetic variation occurs independently of genetic influences remains uncertain, due to the widespread impact of genetics on epigenetic variation and the limited availability of comprehensive (epi)genomic resources from most species. While epigenetic marks can be inherited independently of genetic sequences in some species, there is little evidence that such "transgenerational inheritance" is a general phenomenon. Rather, molecular mechanisms of epigenetic inheritance are highly variable between species.

Development of human amniotic epithelial cell-derived extracellular vesicles as cell-free therapy for dry eye disease.

PURPOSE: This study aimed to investigate the therapeutic potential of extracellular vesicles (EVs) derived from human amniotic epithelial cells (hAEC-EVs) for Dry Eye Disease (DED) treatment. METHODS: Highly purified EVs were isolated from the culture supernatants of hAECs, which obtained from term placenta and characterized. Proteomic contents were analyzed for assessing its biological function related to the therapeutic potentials for DED. Subsequently, we examined the therapeutic efficacy of hAEC-EVs on human corneal epithelial cells exposed to hyperosmotic stress and in an experimental DED mouse model induced by desiccation stress. RESULTS: Proteomic analysis of hAEC-EVs revealed proteins linked to cell proliferation and anti-inflammatory responses. We demonstrated efficient uptake of hAEC-EVs by ocular surface cells. Under DED conditions, EV treatment increased corneal epithelial cell proliferation and migration, and concurrently reducing inflammatory cytokines. In the DED mouse model, hAEC-EVs showed significant improvements in corneal staining score, tear secretion, corneal irregularity, and conjunctival goblet cell density. Additionally, hAEC-EVs exhibited a mitigating effect on ocular surface inflammation induced by desiccation. CONCLUSIONS: These findings suggest that hAEC-EVs hold potential as a cell-free therapy for corneal epithelial defects and ocular surface diseases, presenting a promising treatment option for DED.

Whole-genome DNA methylomes of Tree shrew brains reveal conserved and divergent roles of DNA methylation on sex chromosome regulation.

The tree shrew (Tupaia belangeri) is a promising emerging model organism in biomedical studies, notably due to their evolutionary proximity to primates. To enhance our understanding of how DNA methylation is implicated in regulation of gene expression and the X chromosome inactivation (XCI) in tree shrew brains, here we present their first genome-wide, single-base-resolution methylomes integrated with transcriptomes from prefrontal cortices. We discovered both divergent and conserved features of tree shrew DNA methylation compared to that of other mammals. DNA methylation levels of promoter and gene body regions are negatively correlated with gene expression, consistent with patterns in other mammalian brains studied. Comparing DNA methylation patterns of the female and male X chromosomes, we observed a clear and significant global reduction (hypomethylation) of DNA methylation across the entire X chromosome in females. Our data suggests that the female X hypomethylation does not directly contribute to the gene silencing of the inactivated X chromosome nor does it significantly drive sex-specific gene expression of tree shrews. However, we identified a putative regulatory region in the 5' end of the X inactive specific transcript (Xist) gene, a key gene for XCI, whose pattern of differential DNA methylation strongly relate to its differential expression between male and female tree shrews. We show that differential methylation of this region is conserved across different species. Moreover, we provide evidence suggesting that the observed difference between human and tree shrew X-linked promoter methylation is associated with the difference in genomic CpG contents. Our study offers novel information on genomic DNA methylation of tree shrews, as well as insights into the evolution of X chromosome regulation in mammals.

Human Brain Aging is Associated with Dysregulation of Cell-Type Epigenetic Identity.

Significant links between aging and DNA methylation are emerging from recent studies. On the one hand, DNA methylation undergoes changes with age, a process termed as epigenetic drift. On the other hand, DNA methylation serves as a readily accessible and accurate biomarker for aging. A key missing piece of information, however, is the molecular mechanisms underlying these processes, and how they are related, if any. Addressing the limitations of previous research due to the limited number of investigated CpGs and the heterogeneous nature of tissue samples, here we have examined DNA methylation of over 20 million CpGs across a broad age span in neurons and non-neuronal cells, primarily oligodendrocytes. We show that aging is a primary predictor of DNA methylation variation, surpassing the influence of factors such as sex and schizophrenia diagnosis, among others. On the genome-wide scale, epigenetic drift manifests as significant yet subtle trends that are influenced by the methylation level of individual CpGs. We reveal that CpGs that are highly differentiated between cell types are especially prone to age-associated DNA methylation alterations, leading to the divergence of epigenetic cell type identities as individuals age. On the other hand, CpGs that are included in commonly used epigenetic clocks tend to be those sites that are not highly cell type differentiated. Therefore, dysregulation of epigenetic cell-type identities and current DNA epigenetic clocks represent distinct features of age-associated DNA methylation alterations.

DNA methylation differences between the female and male X chromosomes in human brain.

The mechanisms of X chromosome inactivation suggest fundamental epigenetic differences between the female and male X chromosomes. However, DNA methylation studies often exclude the X chromosomes. In addition, many previous studies relied on techniques that examine non-randomly selected subsets of positions such as array-based methods, rather than assessing the whole X chromosome. Consequently, our understanding of X chromosome DNA methylation lags behind that of autosomes. Here we addressed this gap of knowledge by studying X chromosome DNA methylation using 89 whole genome bisulfite sequencing (WGBS) maps from neurons and oligodendrocytes. Using this unbiased and comprehensive data, we show that DNA methylation of the female X chromosomes is globally reduced (hypomethylated) across the entire chromosome compared to the male X chromosomes and autosomes. On the other hand, the majority of X-linked promoters were more highly methylated (hypermethylated) in females compared to males, consistent with the role of DNA methylation in X chromosome inactivation and dosage compensation. Remarkably, hypermethylation of female X promoters was limited to a group of previously lowly methylated promoters. The other group of highly methylated promoters were both hyper- and hypo-methylated in females with no obvious association with gene expression. Therefore, X chromosome inactivation by DNA methylation was exclusive to a subset of promoters with distinctive epigenetic feature. Apart from this group of promoters, differentially methylated regions in the female and male X chromosomes were dominated by female hypomethylation. Our study furthers the understanding of X-chromosome dosage regulation by DNA methylation on the chromosomal level as well as on individual gene level.

Potential Role of DNA Methylation as a Driver of Plastic Responses to the Environment Across Cells, Organisms, and Populations.

There is great interest in exploring epigenetic modifications as drivers of adaptive organismal responses to environmental change. Extending this hypothesis to populations, epigenetically driven plasticity could influence phenotypic changes across environments. The canonical model posits that epigenetic modifications alter gene regulation and subsequently impact phenotypes. We first discuss origins of epigenetic variation in nature, which may arise from genetic variation, spontaneous epimutations, epigenetic drift, or variation in epigenetic capacitors. We then review and synthesize literature addressing three facets of the aforementioned model: (i) causal effects of epigenetic modifications on phenotypic plasticity at the organismal level, (ii) divergence of epigenetic patterns in natural populations distributed across environmental gradients, and (iii) the relationship between environmentally induced epigenetic changes and gene expression at the molecular level. We focus on DNA methylation, the most extensively studied epigenetic modification. We find support for environmentally associated epigenetic structure in populations and selection on stable epigenetic variants, and that inhibition of epigenetic enzymes frequently bears causal effects on plasticity. However, there are pervasive confounding issues in the literature. Effects of chromatin-modifying enzymes on phenotype may be independent of epigenetic marks, alternatively resulting from functions and protein interactions extrinsic of epigenetics. Associations between environmentally induced changes in DNA methylation and expression are strong in plants and mammals but notably absent in invertebrates and nonmammalian vertebrates. Given these challenges, we describe emerging approaches to better investigate how epigenetic modifications affect gene regulation, phenotypic plasticity, and divergence among populations.

Molecular features driving cellular complexity of human brain evolution.

Human-specific genomic changes contribute to the unique functionalities of the human brain1-5. The cellular heterogeneity of the human brain6,7 and the complex regulation of gene expression highlight the need to characterize human-specific molecular features at cellular resolution. Here we analysed single-nucleus RNA-sequencing and single-nucleus assay for transposase-accessible chromatin with sequencing datasets for human, chimpanzee and rhesus macaque brain tissue from posterior cingulate cortex. We show a human-specific increase of oligodendrocyte progenitor cells and a decrease of mature oligodendrocytes across cortical tissues. Human-specific regulatory changes were accelerated in oligodendrocyte progenitor cells, and we highlight key biological pathways that may be associated with the proportional changes. We also identify human-specific regulatory changes in neuronal subtypes, which reveal human-specific upregulation of FOXP2 in only two of the neuronal subtypes. We additionally identify hundreds of new human accelerated genomic regions associated with human-specific chromatin accessibility changes. Our data also reveal that FOS::JUN and FOX motifs are enriched in the human-specifically accessible chromatin regions of excitatory neuronal subtypes. Together, our results reveal several new mechanisms underlying the evolutionary innovation of human brain at cell-type resolution.

Comprehensive analysis of TCGA data reveals correlation between DNA methylation and alternative splicing.

The effect of DNA methylation on the regulation of gene expression has been extensively discussed in the literature. However, the potential association between DNA methylation and alternative splicing is not understood well. In this study, we integrated multiple omics data types from The Cancer Genome Atlas (TCGA) and systematically examined the relationship between DNA methylation and alternative splicing. Using the methylation data and exon expression data, we identified many CpG sites significantly associated with exon expression in various types of cancers. We further observed that the direction and strength of significant CpG-exon correlation tended to be consistent across different cancer contexts, indicating that some CpG-exon correlation patterns reflect fundamental biological mechanisms that transcend tissue- and cancer- types. We also discovered that CpG sites correlated with exon expressions were more likely to be associated with patient survival outcomes compared to CpG sites that did not correlate with exon expressions. Furthermore, we found that CpG sites were more strongly correlated with exon expression than expression of isoforms harboring the corresponding exons. This observation suggests that a major effect of CpG methylation on alternative splicing may be related to the inclusion or exclusion of exons, which subsequently impacts the relative usage of various isoforms. Overall, our study revealed correlation patterns between DNA methylation and alternative splicing, which provides new insights into the role of methylation in the transcriptional process.

Comparative studies of genomic and epigenetic factors influencing transcriptional variation in two insect species.

Different genes show different levels of expression variability. For example, highly expressed genes tend to exhibit less expression variability. Genes whose promoters have TATA box and initiator motifs tend to have increased expression variability. On the other hand, DNA methylation of transcriptional units, or gene body DNA methylation, is associated with reduced gene expression variability in many species. Interestingly, some insect lineages, most notably Diptera including the canonical model insect Drosophila melanogaster, have lost DNA methylation. Therefore, it is of interest to determine whether genomic features similarly influence gene expression variability in lineages with and without DNA methylation. We analyzed recently generated large-scale data sets in D. melanogaster and honey bee (Apis mellifera) to investigate these questions. Our analysis shows that increased gene expression levels are consistently associated with reduced expression variability in both species, while the presence of TATA box is consistently associated with increased gene expression variability. In contrast, initiator motifs and gene lengths have weak effects limited to some data sets. Importantly, we show that a sequence characteristics indicative of gene body DNA methylation is strongly and negatively associate with gene expression variability in honey bees, while it shows no such association in D. melanogaster. These results suggest the evolutionary loss of DNA methylation in some insect lineages has reshaped the molecular mechanisms concerning the regulation of gene expression variability.

Novel genome sequence of Chinese cavefish (Triplophysa rosa) reveals pervasive relaxation of natural selection in cavefish genomes.

All cavefishes, living exclusively in caves across the globe, exhibit similar phenotypic traits, including the characteristic loss of eyes. To understand whether such phenotypic convergence shares similar genomic bases, here we investigated genome-wide evolutionary signatures of cavefish phenotypes by comparing whole-genome sequences of three pairs of cavefishes and their surface fish relatives. Notably, we newly sequenced and generated a whole-genome assembly of the Chinese cavefish Triplophysa rosa. Our comparative analyses revealed several shared features of cavefish genome evolution. Cavefishes had lower mutation rates than their surface fish relatives. In contrast, the ratio of nonsynonymous to synonymous substitutions (ω) was significantly elevated in cavefishes compared to in surface fishes, consistent with the relaxation of purifying selection. In addition, cavefish genomes had an increased mutational load, including mutations that alter protein hydrophobicity profiles, which were considered harmful. Interestingly, however, we found no overlap in positively selected genes among different cavefish lineages, indicating that the phenotypic convergence in cavefishes was not caused by positive selection of the same sets of genes. Analyses of previously identified candidate genes associated with cave phenotypes supported this conclusion. Genes belonging to the lipid metabolism functional ontology were under relaxed purifying selection in all cavefish genomes, which may be associated with the nutrient-poor habitat of cavefishes. Our work reveals previously uncharacterized patterns of cavefish genome evolution and provides comparative insights into the evolution of cave-associated phenotypic traits.

Dynamic molecular evolution of a supergene with suppressed recombination in white-throated sparrows.

In white-throated sparrows, two alternative morphs differing in plumage and behavior segregate with a large chromosomal rearrangement. As with sex chromosomes such as the mammalian Y, the rearranged version of chromosome two (ZAL2m) is in a near-constant state of heterozygosity, offering opportunities to investigate both degenerative and selective processes during the early evolutionary stages of 'supergenes.' Here, we generated, synthesized, and analyzed extensive genome-scale data to better understand the forces shaping the evolution of the ZAL2 and ZAL2m chromosomes in this species. We found that features of ZAL2m are consistent with substantially reduced recombination and low levels of degeneration. We also found evidence that selective sweeps took place both on ZAL2m and its standard counterpart, ZAL2, after the rearrangement event. Signatures of positive selection were associated with allelic bias in gene expression, suggesting that antagonistic selection has operated on gene regulation. Finally, we discovered a region exhibiting long-range haplotypes inside the rearrangement on ZAL2m. These haplotypes appear to have been maintained by balancing selection, retaining genetic diversity within the supergene. Together, our analyses illuminate mechanisms contributing to the evolution of a young chromosomal polymorphism, revealing complex selective processes acting concurrently with genetic degeneration to drive the evolution of supergenes.

Target ischemic stroke model creation method using photoacoustic microscopy with simultaneous vessel monitoring and dynamic photothrombosis induction.

The ischemic stroke animal model evaluates the efficacy of reperfusion and neuroprotective strategies for ischemic injuries. Various conventional methods have been reported to induce the ischemic models; however, controlling specific neurological deficits, mortality rates, and the extent of the infarction is difficult as the size of the affected region is not precisely controlled. In this paper, we report a single laser-based localized target ischemic stroke model development method by simultaneous vessel monitoring and photothrombosis induction using photoacoustic microscopy (PAM), which has minimized the infarct size at precise location with high reproducibility. The proposed method has significantly reduced the infarcted region by illuminating the precise localization. The reproducibility and validity of suggested method have been demonstrated through repeated experiments and histological analyses. These results demonstrate that our method can provide the ischemic stroke model closest to the clinical pathology for brain ischemia research from inducement, occurrence mechanisms to the recovery process.

Dnmt1a is essential for gene body methylation and the regulation of the zygotic genome in a wasp.

Gene body methylation (GBM) is an ancestral mode of DNA methylation whose role in development has been obscured by the more prominent roles of promoter and CpG island methylation. The wasp Nasonia vitripennis has little promoter and CpG island methylation, yet retains strong GBM, making it an excellent model for elucidating the roles of GBM. Here we show that N. vitripennis DNA methyltransferase 1a (Nv-Dnmt1a) knockdown leads to failures in cellularization and gastrulation of the embryo. Both of these disrupted events are hallmarks of the maternal-zygotic transition (MZT) in insects. Analysis of the embryonic transcriptome and methylome revealed strong reduction of GBM and widespread disruption of gene expression during embryogenesis after Nv-Dnmt1a knockdown. Strikingly, there was a strong correlation between loss of GBM and reduced gene expression in thousands of methylated loci, consistent with the hypothesis that GBM directly facilitates high levels of transcription. We propose that lower expression levels of methylated genes due to reduced GBM is the crucial direct effect of Nv-Dnmt1 knockdown. Subsequently, the disruption of methylated genes leads to downstream dysregulation of the MZT, culminating in developmental failure at gastrulation.

On the origin and evolution of SARS-CoV-2.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the ongoing global outbreak of a coronavirus disease (herein referred to as COVID-19). Other viruses in the same phylogenetic group have been responsible for previous regional outbreaks, including SARS and MERS. SARS-CoV-2 has a zoonotic origin, similar to the causative viruses of these previous outbreaks. The repetitive introduction of animal viruses into human populations resulting in disease outbreaks suggests that similar future epidemics are inevitable. Therefore, understanding the molecular origin and ongoing evolution of SARS-CoV-2 will provide critical insights for preparing for and preventing future outbreaks. A key feature of SARS-CoV-2 is its propensity for genetic recombination across host species boundaries. Consequently, the genome of SARS-CoV-2 harbors signatures of multiple recombination events, likely encompassing multiple species and broad geographic regions. Other regions of the SARS-CoV-2 genome show the impact of purifying selection. The spike (S) protein of SARS-CoV-2, which enables the virus to enter host cells, exhibits signatures of both purifying selection and ancestral recombination events, leading to an effective S protein capable of infecting human and many other mammalian cells. The global spread and explosive growth of the SARS-CoV-2 population (within human hosts) has contributed additional mutational variability into this genome, increasing opportunities for future recombination.

The impact of epigenetic information on genome evolution.

Epigenetic information affects gene function by interacting with chromatin, while not changing the DNA sequence itself. However, it has become apparent that the interactions between epigenetic information and chromatin can, in fact, indirectly lead to DNA mutations and ultimately influence genome evolution. This review evaluates the ways in which epigenetic information affects genome sequence and evolution. We discuss how DNA methylation has strong and pervasive effects on DNA sequence evolution in eukaryotic organisms. We also review how the physical interactions arising from the connections between histone proteins and DNA affect DNA mutation and repair. We then discuss how a variety of epigenetic mechanisms exert substantial effects on genome evolution by suppressing the movement of transposable elements. Finally, we examine how genome expansion through gene duplication is also partially controlled by epigenetic information. Overall, we conclude that epigenetic information has widespread indirect effects on DNA sequences in eukaryotes and represents a potent cause and constraint of genome evolution. This article is part of the theme issue 'How does epigenetics influence the course of evolution?'

Sharing and Teaching Electrocardiograms to Minimize Infarction (STEMI): reducing diagnostic time for acute coronary occlusion in the emergency department.

BACKGROUND: Limits to ST-Elevation Myocardial Infarction (STEMI) criteria may lead to prolonged diagnostic time for acute coronary occlusion. We aimed to reduce ECG-to-Activation (ETA) time through audit and feedback on STEMI-equivalents and subtle occlusions, without increasing Code STEMIs without culprit lesions. METHODS: This multi-centre, quality improvement initiative reviewed all Code STEMI patients from the emergency department (ED) over a one-year baseline and one-year intervention period. We measured ETA time, from the first ED ECG to the time a Code STEMI was activated. Our intervention strategy involved a grand rounds presentation and an internal website presenting weekly local challenging cases, along with literature on STEMI-equivalents and subtle occlusions. Our outcome measure was ETA time for culprit lesions, our process measure was website views/visits, and our balancing measure was the percentage of Code STEMIs without culprit lesions. RESULTS: There were 51 culprit lesions in the baseline period, and 64 in the intervention period. Median ETA declined from 28.0 min (95% confidence interval [CI] 15.0-45.0) to 8.0 min (95%CI 6.0-15.0). The website garnered 70.4 views/week and 27.7 visitors/week in a group of 80 physicians. There was no change in percentage of Code STEMIs without culprit lesions: 28.2% (95%CI 17.8-38.6) to 20.0% (95%CI 11.2-28.8%). Conclusions Our novel weekly web-based feedback to all emergency physicians was associated with a reduction in ETA time by 20 min, without increasing Code STEMIs without culprit lesions. Local ECG audit and feedback, guided by ETA as a quality metric for acute coronary occlusion, could be replicated in other settings to improve care.

Evolution of DNA methylation in the human brain.

DNA methylation is a critical regulatory mechanism implicated in development, learning, memory, and disease in the human brain. Here we have elucidated DNA methylation changes during recent human brain evolution. We demonstrate dynamic evolutionary trajectories of DNA methylation in cell-type and cytosine-context specific manner. Specifically, DNA methylation in non-CG context, namely CH methylation, has increased (hypermethylation) in neuronal gene bodies during human brain evolution, contributing to human-specific down-regulation of genes and co-expression modules. The effects of CH hypermethylation is particularly pronounced in early development and neuronal subtypes. In contrast, DNA methylation in CG context shows pronounced reduction (hypomethylation) in human brains, notably in cis-regulatory regions, leading to upregulation of downstream genes. We show that the majority of differential CG methylation between neurons and oligodendrocytes originated before the divergence of hominoids and catarrhine monkeys, and harbors strong signal for genetic risk for schizophrenia. Remarkably, a substantial portion of differential CG methylation between neurons and oligodendrocytes emerged in the human lineage since the divergence from the chimpanzee lineage and carries significant genetic risk for schizophrenia. Therefore, recent epigenetic evolution of human cortex has shaped the cellular regulatory landscape and contributed to the increased vulnerability to neuropsychiatric diseases.

Enhancer Pleiotropy, Gene Expression, and the Architecture of Human Enhancer-Gene Interactions.

Enhancers are often studied as noncoding regulatory elements that modulate the precise spatiotemporal expression of genes in a highly tissue-specific manner. This paradigm has been challenged by recent evidence of individual enhancers acting in multiple tissues or developmental contexts. However, the frequency of these enhancers with high degrees of "pleiotropy" out of all putative enhancers is not well understood. Consequently, it is unclear how the variation of enhancer pleiotropy corresponds to the variation in expression breadth of target genes. Here, we use multi-tissue chromatin maps from diverse human tissues to investigate the enhancer-gene interaction architecture while accounting for 1) the distribution of enhancer pleiotropy, 2) the variations of regulatory links from enhancers to target genes, and 3) the expression breadth of target genes. We show that most enhancers are tissue-specific and that highly pleiotropy enhancers account for <1% of all putative regulatory sequences in the human genome. Notably, several genomic features are indicative of increasing enhancer pleiotropy, including longer sequence length, greater number of links to genes, increasing abundance and diversity of encoded transcription factor motifs, and stronger evolutionary conservation. Intriguingly, the number of enhancers per gene remains remarkably consistent for all genes (∼14). However, enhancer pleiotropy does not directly translate to the expression breadth of target genes. We further present a series of Gaussian Mixture Models to represent this organization architecture. Consequently, we demonstrate that a modest trend of more pleiotropic enhancers targeting more broadly expressed genes can generate the observed diversity of expression breadths in the human genome.

Koala methylomes reveal divergent and conserved DNA methylation signatures of X chromosome regulation.

X chromosome inactivation (XCI) mediated by differential DNA methylation between sexes is an iconic example of epigenetic regulation. Although XCI is shared between eutherians and marsupials, the role of DNA methylation in marsupial XCI remains contested. Here, we examine genome-wide signatures of DNA methylation across fives tissues from a male and female koala (Phascolarctos cinereus), and present the first whole-genome, multi-tissue marsupial 'methylome atlas'. Using these novel data, we elucidate divergent versus common features of representative marsupial and eutherian DNA methylation. First, tissue-specific differential DNA methylation in koalas primarily occurs in gene bodies. Second, females show significant global reduction (hypomethylation) of X chromosome DNA methylation compared to males. We show that this pattern is also observed in eutherians. Third, on average, promoter DNA methylation shows little difference between male and female koala X chromosomes, a pattern distinct from that of eutherians. Fourth, the sex-specific DNA methylation landscape upstream of Rsx, the primary lncRNA associated with marsupial XCI, is consistent with the epigenetic regulation of female-specific (and presumably inactive X chromosome-specific) expression. Finally, we use the prominent female X chromosome hypomethylation and classify 98 previously unplaced scaffolds as X-linked, contributing an additional 14.6 Mb (21.5%) to genomic data annotated as the koala X chromosome. Our work demonstrates evolutionarily divergent pathways leading to functionally conserved patterns of XCI in two deep branches of mammals.

Genome-wide variation in DNA methylation linked to developmental stage and chromosomal suppression of recombination in white-throated sparrows.

Much of our knowledge on regulatory impacts of DNA methylation has come from laboratory-bred model organisms, which may not exhibit the full extent of variation found in wild populations. Here, we investigated naturally-occurring variation in DNA methylation in a wild avian species, the white-throated sparrow (Zonotrichia albicollis). This species offers exceptional opportunities for studying the link between genetic differentiation and phenotypic traits because of a nonrecombining chromosome pair linked to both plumage and behavioural phenotypes. Using novel single-nucleotide resolution methylation maps and gene expression data, we show that DNA methylation and the expression of DNA methyltransferases are significantly higher in adults than in nestlings. Genes for which DNA methylation varied between nestlings and adults were implicated in development and cell differentiation and were located throughout the genome. In contrast, differential methylation between plumage morphs was concentrated in the nonrecombining chromosome pair. Interestingly, a large number of CpGs on the nonrecombining chromosome, localized to transposable elements, have undergone dramatic loss of DNA methylation since the split of the ZAL2 and ZAL2m chromosomes. Changes in methylation predicted changes in gene expression for both chromosomes. In summary, we demonstrate changes in genome-wide DNA methylation that are associated with development and with specific functional categories of genes in white-throated sparrows. Moreover, we observe substantial DNA methylation reprogramming associated with the suppression of recombination, with implications for genome integrity and gene expression divergence. These results offer an unprecedented view of ongoing epigenetic reprogramming in a wild population.