A genome-wide study of DNA methylation patterns and gene expression levels in multiple human and chimpanzee tissues

Inferring DNA methylation in non-skeletal tissues of ancient specimens

Genome-wide premortem DNA methylation patterns can be computationally reconstructed from high-coverage DNA sequences of ancient samples. Because DNA methylation is more conserved across species than across tissues, and ancient DNA is typically extracted from bones and teeth, previous works utilizing ancient DNA methylation maps focused on studying evolutionary changes in the skeletal system. Here we suggest that DNA methylation patterns in one tissue may, under certain conditions, be informative on DNA methylation patterns in other tissues of the same individual. Using the fact that tissue-specific DNA methylation builds up during embryonic development, we identified the conditions that allow for such cross-tissue inference and devised an algorithm that carries it out. We trained the algorithm on methylation data from extant species and reached high precisions of up to 0.92 for validation datasets. We then used the algorithm on archaic humans, and identified more than 1,850 positions for which we were able to observe differential DNA methylation in prefrontal cortex neurons. These positions are linked to hundreds of genes, many of which are involved in neural functions such as structural and developmental processes. Six positions are located in the neuroblastoma breaking point family (NBPF) gene family, which probably played a role in human brain evolution. The algorithm we present here allows for the examination of epigenetic changes in tissues and cell types that are absent from the palaeontological record, and therefore provides new ways to study the evolutionary impacts of epigenetic changes.

A role for DNA methylation in bumblebee morphogenesis hints at female-specific developmental erasure

Epigenetic mechanisms, such as DNA methylation, are crucial factors in animal development. In some mammals, almost all DNA methylation is erased during embryo development and re-established in a sex- and cell-specific manner. This erasure and re-establishment is thought to primarily be a vertebrate-specific trait. Insects are particularly interesting in terms of development as many species often undergo remarkable morphological changes en route to maturity, that is, morphogenesis. However, little is known about the role of epigenetic mechanisms in this process across species. We have used whole-genome bisulfite sequencing to track genome-wide DNA methylation changes through the development of an economically and environmentally important pollinator species, the bumblebee Bombus terrestris (Hymenoptera:Apidae Linnaeus). We find overall levels of DNA methylation vary throughout development, and we find developmentally relevant differentially methylated genes throughout. Intriguingly, we have identified a depletion of DNA methylation in ovaries/eggs and an enrichment of highly methylated genes in sperm. We suggest this could represent a sex-specific DNA methylation erasure event. To our knowledge, this is the first suggestion of possible developmental DNA methylation erasure in an insect species. This study lays the required groundwork for functional experimental work to determine if there is a causal nature to the DNA methylation differences identified. Additionally, the application of single-cell methylation sequencing to this system will enable more accurate identification of if or when DNA methylation is erased during development.

Cross-Species Comparative DNA Methylation Reveals Novel Insights into Complex Trait Genetics among Cattle, Sheep, and Goats

The cross-species characterization of evolutionary changes in the functional genome can facilitate the translation of genetic findings across species and the interpretation of the evolutionary basis underlying complex phenotypes. Yet, this has not been fully explored between cattle, sheep, goats, and other mammals. Here, we systematically characterized the evolutionary dynamics of DNA methylation and gene expression in 3 somatic tissues (i.e. brain, liver, and skeletal muscle) and sperm across 7 mammalian species, including 3 ruminant livestock species (cattle, sheep, and goats), humans, pigs, mice, and dogs, by generating and integrating 160 DNA methylation and transcriptomic data sets. We demonstrate dynamic changes of DNA hypomethylated regions and hypermethylated regions in tissue-type manner across cattle, sheep, and goats. Specifically, based on the phylo-epigenetic model of DNA methylome, we identified a total of 25,074 hypomethylated region extension events specific to cattle, which participated in rewiring tissue-specific regulatory network. Furthermore, by integrating genome-wide association studies of 50 cattle traits, we provided novel insights into the genetic and evolutionary basis of complex phenotypes in cattle. Overall, our study provides a valuable resource for exploring the evolutionary dynamics of the functional genome and highlights the importance of cross-species characterization of multiomics data sets for the evolutionary interpretation of complex phenotypes in cattle livestock.

The relationship between regulatory changes in cis and trans and the evolution of gene expression in humans and chimpanzees

BACKGROUND

Comparative gene expression studies in apes are fundamentally limited by the challenges associated with sampling across different tissues. Here, we used single-cell RNA sequencing of embryoid bodies to collect transcriptomic data from over 70 cell types in three humans and three chimpanzees.

RESULTS

We find hundreds of genes whose regulation is conserved across cell types, as well as genes whose regulation likely evolves under directional selection in one or a handful of cell types. Using embryoid bodies from a human-chimpanzee fused cell line, we also infer the proportion of inter-species regulatory differences due to changes in cis and trans elements between the species. Using the cis/trans inference and an analysis of transcription factor binding sites, we identify dozens of transcription factors whose inter-species differences in expression are affecting expression differences between humans and chimpanzees in hundreds of target genes.

CONCLUSIONS

Here, we present the most comprehensive dataset of comparative gene expression from humans and chimpanzees to date, including a catalog of regulatory mechanisms associated with inter-species differences.

Brain DNA methylomic analysis of frontotemporal lobar degeneration reveals OTUD4 in shared dysregulated signatures across pathological subtypes

Frontotemporal lobar degeneration (FTLD) is an umbrella term describing the neuropathology of a clinically, genetically and pathologically heterogeneous group of diseases, including frontotemporal dementia (FTD) and progressive supranuclear palsy (PSP). Among the major FTLD pathological subgroups, FTLD with TDP-43 positive inclusions (FTLD-TDP) and FTLD with tau-positive inclusions (FTLD-tau) are the most common, representing about 90% of the cases. Although alterations in DNA methylation have been consistently associated with neurodegenerative diseases, including Alzheimer's disease and Parkinson's disease, little is known for FTLD and its heterogeneous subgroups and subtypes. The main goal of this study was to investigate DNA methylation variation in FTLD-TDP and FTLD-tau. We used frontal cortex genome-wide DNA methylation profiles from three FTLD cohorts (142 FTLD cases and 92 controls), generated using the Illumina 450K or EPIC microarrays. We performed epigenome-wide association studies (EWAS) for each cohort followed by meta-analysis to identify shared differentially methylated loci across FTLD subgroups/subtypes. In addition, we used weighted gene correlation network analysis to identify co-methylation signatures associated with FTLD and other disease-related traits. Wherever possible, we also incorporated relevant gene/protein expression data. After accounting for a conservative Bonferroni multiple testing correction, the EWAS meta-analysis revealed two differentially methylated loci in FTLD, one annotated to OTUD4 (5'UTR-shore) and the other to NFATC1 (gene body-island). Of these loci, OTUD4 showed consistent upregulation of mRNA and protein expression in FTLD. In addition, in the three independent co-methylation networks, OTUD4-containing modules were enriched for EWAS meta-analysis top loci and were strongly associated with the FTLD status. These co-methylation modules were enriched for genes implicated in the ubiquitin system, RNA/stress granule formation and glutamatergic synaptic signalling. Altogether, our findings identified novel FTLD-associated loci, and support a role for DNA methylation as a mechanism involved in the dysregulation of biological processes relevant to FTLD, highlighting novel potential avenues for therapeutic development.

The contribution of DNA methylation to the (dys)function of oligodendroglia in neurodegeneration

Neurodegenerative diseases encompass a heterogeneous group of conditions characterised by the progressive degeneration of the structure and function of the central or peripheral nervous systems. The pathogenic mechanisms underlying these diseases are not fully understood. However, a central feature consists of regional aggregation of proteins in the brain, such as the accumulation of β-amyloid plaques in Alzheimer's disease (AD), inclusions of hyperphosphorylated microtubule-binding tau in AD and other tauopathies, or inclusions containing α-synuclein in Parkinson's disease (PD), dementia with Lewy bodies (DLB) and multiple system atrophy (MSA). Various pathogenic mechanisms are thought to contribute to disease, and an increasing number of studies implicate dysfunction of oligodendrocytes (the myelin producing cells of the central nervous system) and myelin loss. Aberrant DNA methylation, the most widely studied epigenetic modification, has been associated with many neurodegenerative diseases, including AD, PD, DLB and MSA, and recent findings highlight aberrant DNA methylation in oligodendrocyte/myelin-related genes. Here we briefly review the evidence showing that changes to oligodendrocytes and myelin are key in neurodegeneration, and explore the relevance of DNA methylation in oligodendrocyte (dys)function. As DNA methylation is reversible, elucidating its involvement in pathogenic mechanisms of neurodegenerative diseases and in dysfunction of specific cell-types such as oligodendrocytes may bring opportunities for therapeutic interventions for these diseases.

Comparative three-dimensional genome architectures of adipose tissues provide insight into human-specific regulation of metabolic homeostasis

Elucidating the regulatory mechanisms of human adipose tissues (ATs) evolution is essential for understanding human-specific metabolic regulation, but the functional importance and evolutionary dynamics of three-dimensional (3D) genome organizations of ATs are not well defined. Here, we compared the 3D genome architectures of anatomically distinct ATs from humans and six representative mammalian models. We recognized evolutionarily conserved and human-specific chromatin conformation in ATs at multiple scales, including compartmentalization, topologically associating domain (TAD), and promoter-enhancer interactions (PEI), which have not been described previously. We found PEI are much more evolutionarily dynamic with respect to compartmentalization and topologically associating domain. Compared to conserved PEIs, human-specific PEIs are enriched for human-specific sequence, and the binding motifs of their potential mediators (transcription factors) are less conserved. Our data also demonstrated that genes involved in the evolutionary dynamics of chromatin organization have weaker transcriptional conservation than those associated with conserved chromatin organization. Furthermore, the genes involved in energy metabolism and the maintenance of metabolic homeostasis are enriched in human-specific chromatin organization, while housekeeping genes, health-related genes, and genetic variations are enriched in evolutionarily conserved compared to human-specific chromatin organization. Finally, we showed extensively divergent human-specific 3D genome organizations among one subcutaneous and three visceral ATs. Together, these findings provide a global overview of 3D genome architecture dynamics between ATs from human and mammalian models and new insights into understanding the regulatory evolution of human ATs.

Comparative analysis reveals epigenomic evolution related to species traits and genomic imprinting in mammals

DNA methylation is an epigenetic modification that plays a crucial role in various regulatory processes, including gene expression regulation, transposable element repression, and genomic imprinting. However, most studies on DNA methylation have been conducted in humans and other model species, whereas the dynamics of DNA methylation across mammals remain poorly explored, limiting our understanding of epigenomic evolution in mammals and the evolutionary impacts of conserved and lineage-specific DNA methylation. Here, we generated and gathered comparative epigenomic data from 13 mammalian species, including two marsupial species, to demonstrate that DNA methylation plays critical roles in several aspects of gene evolution and species trait evolution. We found that the species-specific DNA methylation of promoters and noncoding elements correlates with species-specific traits such as body patterning, indicating that DNA methylation might help establish or maintain interspecies differences in gene regulation that shape phenotypes. For a broader view, we investigated the evolutionary histories of 88 known imprinting control regions across mammals to identify their evolutionary origins. By analyzing the features of known and newly identified potential imprints in all studied mammals, we found that genomic imprinting may function in embryonic development through the binding of specific transcription factors. Our findings show that DNA methylation and the complex interaction between the genome and epigenome have a significant impact on mammalian evolution, suggesting that evolutionary epigenomics should be incorporated to develop a unified evolutionary theory.

Associations Between Levels of Peripheral NCAPH2 Promoter Methylation and Different Stages of Alzheimer's Disease: A Cross-Sectional Study

BACKGROUND

Several studies have examined NCAPH2 methylation in amnestic mild cognitive impairment (aMCI) and Alzheimer's disease (AD), but little is known of NCAPH2 methylation in subjective cognitive decline (SCD).

OBJECTIVE

To examine whether methylation of peripheral NCAPH2 are differentially changed at various phases of AD, and whether it could serve as a diagnostic biomarker for SCD.

METHODS

A total of 40 AD patients, 52 aMCI patients, 148 SCD patients, and 193 cognitively normal controls (NCs) were recruited in the current case-control study. Besides, 54 cognitively normal individuals have received amyloid positron emission tomography (amyloid PET) scans. Using bisulfite pyrosequencing method, we measured blood DNA methylation in the NCAPH2 gene promoter.

RESULTS

The main outcomes were: 1) For SCD, there was no significant difference between SCD and NC regarding NCAPH2 methylation; 2) For aMCI, NCAPH2 methylation at CpG2 were significantly lower in aMCI compared with NC and SCD in the entire population and male subgroup; 3) For AD, NCAPH2 methylation at CpG1 were significantly lower in AD compared with NC among females; 4) A relationship with apolipoprotein E (APOE) ɛ4 status was shown. Receiver operating characteristic (ROC) analysis by combining NCAPH2 methylation, age, education, and APOEɛ4 status could distinguish between patients with aMCI (area under the curve (AUC): 0.742) and AD (AUC: 0.873) from NCs.

CONCLUSION

NCAPH2 methylation levels were altered at the aMCI and AD stage and may be convenient and cost-effective biomarkers of AD and aMCI.

Exploiting Signal Joint T Cell Receptor Excision Circle to Investigate the Impact of COVID-19 and Autoimmune Diseases on Age Prediction and Immunosenescence

Signal joint T cell receptor excision circles (sjTRECs) are a promising marker for age estimation and immunosenescence in different ethnic groups. Several limitations are expected to overshadow their use as accurate markers for age prediction. The current study was conducted to determine the influence of immunologic disorders, such as autoimmune diseases and COVID-19, on the accuracy of sjTRECs as molecular markers for age estimation and immunosenescence among living Egyptians. Peripheral blood sjTRECs level was measured by qPCR in 90 autoimmune patients, 58 COVID-19 patients, and 85 healthy controls. The mean dCt values were significantly (p = 0.0002) different between the three groups, with the highest values in healthy subjects, followed by autoimmune and COVID-19 patients. A significant negative correlation was identified between the sjTRECs levels and ages in all studied cases. There were significant positive correlations between chronological age and predicted age for healthy individuals, autoimmune, and COVID-19 patients with mean absolute deviations (MAD) of 9.40, 11.04, and 9.71, respectively. The two patients' groups exhibited early immunosenescence, which was more noticeable among the young adults with COVID-19 and autoimmune patients of age range (18-49 years). Autoimmunity may represent a critical factor impacting the accuracy of sjTRECs quantitation for age prediction.

Genome-Wide DNA Methylation Profile Indicates Potential Epigenetic Regulation of Aging in the Rhesus Macaque Thymus

We analyzed whole-genome bisulfite sequencing (WGBS) and RNA sequencing data of two young (1 year old) and two adult (9 years old) rhesus macaques (Macaca mulatta) to characterize the genomic DNA methylation profile of the thymus and explore the molecular mechanism of age-related changes in the thymus. Combining the two-omics data, we identified correlations between DNA methylation and gene expression and found that DNA methylation played an essential role in the functional changes of the aging thymus, especially in immunity and coagulation. The hypomethylation levels of C3 and C5AR2 and the hypermethylation level of C7 may lead to the high expressions of these genes in adult rhesus macaque thymuses, thus activating the classical complement pathway and the alternative pathway and enhancing their innate immune function. Adult thymuses had an enhanced coagulation pathway, which may have resulted from the hypomethylation and upregulated expressions of seven coagulation-promoting factor genes (F13A1, CLEC4D, CLEC4E, FCN3, PDGFRA, FGF2 and FGF7) and the hypomethylation and low expression of CPB2 to inhibit the degradation of blood clots. Furthermore, the functional decline in differentiation, activation and maturation of T cells in adult thymuses was also closely related to the changes in methylation levels and gene expression levels of T cell development genes (CD3G, GAD2, ADAMDEC1 and LCK) and the thymogenic hormone gene TMPO. A comparison of the age-related methylated genes among four mammal species revealed that most of the epigenetic clocks were species-specific. Furthermore, based on the genomic landscape of allele-specific DNA methylation, we identified several age-related clustered sequence-dependent allele-specific DNA methylated (cS-ASM) genes. Overall, these DNA methylation patterns may also help to assist with understanding the mechanisms of the aging thymus with the epigenome.

Sequence divergence and retrotransposon insertion underlie interspecific epigenetic differences in primates

Changes in the epigenome can affect the phenotype without the presence of changes in the genomic sequence. Given the high identity of the human and chimpanzee genome sequences, a substantial portion of their phenotypic divergence likely arises from epigenomic differences between the two species. In this study, the transcriptome and epigenome were determined for induced pluripotent stem cells (iPSCs) generated from human and chimpanzee individuals. The transcriptome and epigenomes for trimethylated histone H3 at lysine-4 (H3K4me3) and lysine-27 (H3K27me3) showed high levels of similarity between the two species. However, there were some differences in histone modifications. Although such regions, in general, did not show significant enrichment of interspecies nucleotide variations, gains in binding motifs for pluripotency-related transcription factors, especially POU5F1 and SOX2, were frequently found in species-specific H3K4me3 regions. We also revealed that species-specific insertions of retrotransposons, including the LTR5_Hs subfamily in human and a newly identified LTR5_Pt subfamily in chimpanzee, created species-specific H3K4me3 regions associated with increased expression of nearby genes. Human iPSCs have more species-specific H3K27me3 regions, resulting in more abundant bivalent domains. Only a limited number of these species-specific H3K4me3 and H3K27me3 regions overlap with species-biased enhancers in cranial neural crest cells, suggesting that differences in the epigenetic state of developmental enhancers appear late in development. Therefore, iPSCs serve as a suitable starting material for studying evolutionary changes in epigenome dynamics during development.

Epigenetic aging of the prefrontal cortex and cerebellum in humans and chimpanzees

Epigenetic age has emerged as an important biomarker of biological aging. It has revealed that some tissues age faster than others, which is vital to understanding the complex phenomenon of aging and developing effective interventions. Previous studies have demonstrated that humans exhibit heterogeneity in pace of epigenetic aging among brain structures that are consistent with differences in structural and microanatomical deterioration. Here, we add comparative data on epigenetic brain aging for chimpanzees, humans' closest relatives. Such comparisons can further our understanding of which aspects of human aging are evolutionarily conserved or specific to our species, especially given that humans are distinguished by a long lifespan, large brain, and, potentially, more severe neurodegeneration with age. Specifically, we investigated epigenetic aging of the dorsolateral prefrontal cortex and cerebellum, of humans and chimpanzees by generating genome-wide CpG methylation data and applying established epigenetic clock algorithms to produce estimates of biological age for these tissues. We found that both species exhibit relatively slow epigenetic aging in the brain relative to blood. Between brain structures, humans show a faster rate of epigenetic aging in the dorsolateral prefrontal cortex compared to the cerebellum, which is consistent with previous findings. Chimpanzees, in contrast, show comparable rates of epigenetic aging in the two brain structures. Greater epigenetic change in the human dorsolateral prefrontal cortex compared to the cerebellum may reflect both the protracted development of this structure in humans and its greater age-related vulnerability to neurodegenerative pathology.

Long-range promoter-enhancer contacts are conserved during evolution and contribute to gene expression robustness

Gene expression is regulated through complex molecular interactions, involving cis-acting elements that can be situated far away from their target genes. Data on long-range contacts between promoters and regulatory elements are rapidly accumulating. However, it remains unclear how these regulatory relationships evolve and how they contribute to the establishment of robust gene expression profiles. Here, we address these questions by comparing genome-wide maps of promoter-centered chromatin contacts in mouse and human. We show that there is significant evolutionary conservation of cis-regulatory landscapes, indicating that selective pressures act to preserve not only regulatory element sequences but also their chromatin contacts with target genes. The extent of evolutionary conservation is remarkable for long-range promoter–enhancer contacts, illustrating how the structure of regulatory landscapes constrains large-scale genome evolution. We show that the evolution of cis-regulatory landscapes, measured in terms of distal element sequences, synteny, or contacts with target genes, is significantly associated with gene expression evolution.

Skeletal muscle methylome and transcriptome integration reveals profound sex differences related to muscle function and substrate metabolism

Nearly all human complex traits and diseases exhibit some degree of sex differences, with epigenetics being one of the main contributing factors. Various tissues display sex differences in DNA methylation; however, this has not yet been explored in skeletal muscle, despite skeletal muscle being among the tissues with the most transcriptomic sex differences. For the first time, we investigated the effect of sex on autosomal DNA methylation in human skeletal muscle across three independent cohorts (Gene SMART, FUSION, and GSE38291) using a meta-analysis approach, totalling 369 human muscle samples (222 males and 147 females), and integrated this with known sex-biased transcriptomics. We found 10,240 differentially methylated regions (DMRs) at FDR < 0.005, 94% of which were hypomethylated in males, and gene set enrichment analysis revealed that differentially methylated genes were involved in muscle contraction and substrate metabolism. We then investigated biological factors underlying DNA methylation sex differences and found that circulating hormones were not associated with differential methylation at sex-biased DNA methylation loci; however, these sex-specific loci were enriched for binding sites of hormone-related transcription factors (with top TFs including androgen (AR), estrogen (ESR1), and glucocorticoid (NR3C1) receptors). Fibre type proportions were associated with differential methylation across the genome, as well as across 16% of sex-biased DNA methylation loci (FDR < 0.005). Integration of DNA methylomic results with transcriptomic data from the GTEx database and the FUSION cohort revealed 326 autosomal genes that display sex differences at both the epigenome and transcriptome levels. Importantly, transcriptional sex-biased genes were overrepresented among epigenetic sex-biased genes (p value = 4.6e−13), suggesting differential DNA methylation and gene expression between male and female muscle are functionally linked. Finally, we validated expression of three genes with large effect sizes (FOXO3A, ALDH1A1, and GGT7) in the Gene SMART cohort with qPCR. GGT7, involved in antioxidant metabolism, displays male-biased expression as well as lower methylation in males across the three cohorts. In conclusion, we uncovered 8420 genes that exhibit DNA methylation differences between males and females in human skeletal muscle that may modulate mechanisms controlling muscle metabolism and health.

MOBP and HIP1 in multiple system atrophy: New α-synuclein partners in glial cytoplasmic inclusions implicated in the disease pathogenesis

AIMS

Multiple system atrophy (MSA) is a fatal neurodegenerative disease. Similar to Parkinson's disease (PD), MSA is an α-synucleinopathy, and its pathological hallmark consists of glial cytoplasmic inclusions (GCIs) containing α-synuclein (SNCA) in oligodendrocytes. We previously identified consistent changes in myelin-associated oligodendrocyte basic protein (MOBP) and huntingtin interacting protein 1 (HIP1) DNA methylation status in MSA. We hypothesized that if differential DNA methylation at these loci is mechanistically relevant for MSA, it should have downstream consequences on gene regulation.

METHODS

We investigated the relationship between MOBP and HIP1 DNA methylation and mRNA levels in cerebellar white matter from MSA and healthy controls. Additionally, we analysed protein expression using western blotting, immunohistochemistry and proximity ligation assays.

RESULTS

We found decreased MOBP mRNA levels significantly correlated with increased DNA methylation in MSA. For HIP1, we found a distinct relationship between DNA methylation and gene expression levels in MSA compared to healthy controls, suggesting this locus may be subjected to epigenetic remodelling in MSA. Although soluble protein levels for MOBP and HIP1 in cerebellar white matter were not significantly different between MSA cases and controls, we found striking differences between MSA and other neurodegenerative diseases, including PD and Huntington's disease. We also found that MOBP and HIP1 are mislocalized into the GCIs in MSA, where they appear to interact with SNCA.

CONCLUSIONS

This study supports a role for DNA methylation in downregulation of MOBP mRNA in MSA. Most importantly, the identification of MOBP and HIP1 as new constituents of GCIs emphasizes the relevance of these two loci to the pathogenesis of MSA.

Neurodegenerative movement disorders: An epigenetics perspective and promise for the future

Neurodegenerative movement disorders (NMDs) are age‐dependent disorders that are characterised by the degeneration and loss of neurons, typically accompanied by pathological accumulation of different protein aggregates in the brain, which lead to motor symptoms. NMDs include Parkinson's disease, multiple system atrophy, progressive supranuclear palsy, and Huntington's disease, among others. Epigenetic modifications are responsible for functional gene regulation during development, adult life and ageing and have progressively been implicated in complex diseases such as cancer and more recently in neurodegenerative diseases, such as NMDs. DNA methylation is by far the most widely studied epigenetic modification and consists of the reversible addition of a methyl group to the DNA without changing the DNA sequence. Although this research field is still in its infancy in relation to NMDs, an increasing number of studies point towards a role for DNA methylation in disease processes. This review addresses recent advances in epigenetic and epigenomic research in NMDs, with a focus on human brain DNA methylation studies. We discuss the current understanding of the DNA methylation changes underlying these disorders, the potential for use of these DNA modifications in peripheral tissues as biomarkers in early disease detection, classification and progression as well as a promising role in future disease management and therapy.

Epigenomic profiling of primate lymphoblastoid cell lines reveals the evolutionary patterns of epigenetic activities in gene regulatory architectures

Changes in the epigenetic regulation of gene expression have a central role in evolution. Here, we extensively profiled a panel of human, chimpanzee, gorilla, orangutan, and macaque lymphoblastoid cell lines (LCLs), using ChIP-seq for five histone marks, ATAC-seq and RNA-seq, further complemented with whole genome sequencing (WGS) and whole genome bisulfite sequencing (WGBS). We annotated regulatory elements (RE) and integrated chromatin contact maps to define gene regulatory architectures, creating the largest catalog of RE in primates to date. We report that epigenetic conservation and its correlation with sequence conservation in primates depends on the activity state of the regulatory element. Our gene regulatory architectures reveal the coordination of different types of components and highlight the role of promoters and intragenic enhancers (gE) in the regulation of gene expression. We observe that most regulatory changes occur in weakly active gE. Remarkably, novel human-specific gE with weak activities are enriched in human-specific nucleotide changes. These elements appear in genes with signals of positive selection and human acceleration, tissue-specific expression, and particular functional enrichments, suggesting that the regulatory evolution of these genes may have contributed to human adaptation.

Characterizing Genetic Regulatory Elements in Ovine Tissues

The Ovine Functional Annotation of Animal Genomes (FAANG) project, part of the broader livestock species FAANG initiative, aims to identify and characterize gene regulatory elements in domestic sheep. Regulatory element annotation is essential for identifying genetic variants that affect health and production traits in this important agricultural species, as greater than 90% of variants underlying genetic effects are estimated to lie outside of transcribed regions. Histone modifications that distinguish active or repressed chromatin states, CTCF binding, and DNA methylation were used to characterize regulatory elements in liver, spleen, and cerebellum tissues from four yearling sheep. Chromatin immunoprecipitation with sequencing (ChIP-seq) was performed for H3K4me3, H3K27ac, H3K4me1, H3K27me3, and CTCF. Nine chromatin states including active promoters, active enhancers, poised enhancers, repressed enhancers, and insulators were characterized in each tissue using ChromHMM. Whole-genome bisulfite sequencing (WGBS) was performed to determine the complement of whole-genome DNA methylation with the ChIP-seq data. Hypermethylated and hypomethylated regions were identified across tissues, and these locations were compared with chromatin states to better distinguish and validate regulatory elements in these tissues. Interestingly, chromatin states with the poised enhancer mark H3K4me1 in the spleen and cerebellum and CTCF in the liver displayed the greatest number of hypermethylated sites. Not surprisingly, active enhancers in the liver and spleen, and promoters in the cerebellum, displayed the greatest number of hypomethylated sites. Overall, chromatin states defined by histone marks and CTCF occupied approximately 22% of the genome in all three tissues. Furthermore, the liver and spleen displayed in common the greatest percent of active promoter (65%) and active enhancer (81%) states, and the liver and cerebellum displayed in common the greatest percent of poised enhancer (53%), repressed enhancer (68%), hypermethylated sites (75%), and hypomethylated sites (73%). In addition, both known and de novo CTCF-binding motifs were identified in all three tissues, with the highest number of unique motifs identified in the cerebellum. In summary, this study has identified the regulatory regions of genes in three tissues that play key roles in defining health and economically important traits and has set the precedent for the characterization of regulatory elements in ovine tissues using the Rambouillet reference genome.

Comparative analysis reveals distinctive epigenetic features of the human cerebellum

Identifying the molecular underpinnings of the neural specializations that underlie human cognitive and behavioral traits has long been of considerable interest. Much research on human-specific changes in gene expression and epigenetic marks has focused on the prefrontal cortex, a brain structure distinguished by its role in executive functions. The cerebellum shows expansion in great apes and is gaining increasing attention for its role in motor skills and cognitive processing, including language. However, relatively few molecular studies of the cerebellum in a comparative evolutionary context have been conducted. Here, we identify human-specific methylation in the lateral cerebellum relative to the dorsolateral prefrontal cortex, in a comparative study with chimpanzees (Pan troglodytes) and rhesus macaques (Macaca mulatta). Specifically, we profiled genome-wide methylation levels in the three species for each of the two brain structures and identified human-specific differentially methylated genomic regions unique to each structure. We further identified which differentially methylated regions (DMRs) overlap likely regulatory elements and determined whether associated genes show corresponding species differences in gene expression. We found greater human-specific methylation in the cerebellum than the dorsolateral prefrontal cortex, with differentially methylated regions overlapping genes involved in several conditions or processes relevant to human neurobiology, including synaptic plasticity, lipid metabolism, neuroinflammation and neurodegeneration, and neurodevelopment, including developmental disorders. Moreover, our results show some overlap with those of previous studies focused on the neocortex, indicating that such results may be common to multiple brain structures. These findings further our understanding of the cerebellum in human brain evolution.

Humans are distinguished from other species by several aspects of cognition. While much comparative evolutionary neuroscience has focused on the neocortex, increasing recognition of the cerebellum’s role in cognition and motor processing has inspired considerable new research. Comparative molecular studies, however, generally continue to focus on the neocortex. We sought to characterize potential genetic regulatory traits distinguishing the human cerebellum by undertaking genome-wide epigenetic profiling of the lateral cerebellum, and compared this to the prefrontal cortex of humans, chimpanzees, and rhesus macaque monkeys. We found that humans showed greater differential CpG methylation–an epigenetic modification of DNA that can reflect past or present gene expression–in the cerebellum than the prefrontal cortex, highlighting the importance of this structure in human brain evolution. Humans also specifically show methylation differences at genes involved in neurodevelopment, neuroinflammation, synaptic plasticity, and lipid metabolism. These differences are relevant for understanding processes specific to humans, such as extensive plasticity, as well as pronounced and prevalent neurodegenerative conditions associated with aging.

Role of promoters in regulating alternative splicing

Alternative splicing (AS) plays a critical role in enhancing proteome complexity in higher eukaryotes. Almost all the multi intron-containing genes undergo AS in humans. Splicing mainly occurs co-transcriptionally, where RNA polymerase II (RNA pol II) plays a crucial role in coordinating transcription and pre-mRNA splicing. Aberrant AS leads to non-functional proteins causative in various pathophysiological conditions such as cancers, neurodegenerative diseases, and muscular dystrophies. Transcription and pre-mRNA splicing are deeply interconnected and can influence each other's functions. Several studies evinced that specific promoters employed by RNA pol II dictate the RNA processing decisions. Promoter-specific recruitment of certain transcriptional factors or transcriptional coactivators influences splicing, and the extent to which these factors affect splicing has not been discussed in detail. Here, in this review, various DNA-binding proteins and their influence on promoter-specific AS are extensively discussed. Besides, this review highlights how the promoter-specific epigenetic changes might regulate AS.

Divergence in alternative polyadenylation contributes to gene regulatory differences between humans and chimpanzees

While comparative functional genomic studies have shown that inter-species differences in gene expression can be explained by corresponding inter-species differences in genetic and epigenetic regulatory mechanisms, co-transcriptional mechanisms, such as alternative polyadenylation (APA), have received little attention. We characterized APA in lymphoblastoid cell lines from six humans and six chimpanzees by identifying and estimating the usage for 44,432 polyadenylation sites (PAS) in 9518 genes. Although APA is largely conserved, 1705 genes showed significantly different PAS usage (FDR 0.05) between species. Genes with divergent APA also tend to be differentially expressed, are enriched among genes showing differences in protein translation, and can explain a subset of observed inter-species protein expression differences that do not differ at the transcript level. Finally, we found that genes with a dominant PAS, which is used more often than other PAS, are particularly enriched for differentially expressed genes.

Epigenomic differences in the human and chimpanzee genomes are associated with structural variation

Structural variation (SV), including insertions and deletions (indels), is a primary mechanism of genome evolution. However, the mechanism by which SV contributes to epigenome evolution is poorly understood. In this study, we characterized the association between lineage-specific indels and epigenome differences between human and chimpanzee to investigate how SVs might have shaped the epigenetic landscape. By intersecting medium-to-large human-chimpanzee indels (20 bp-50 kb) with putative promoters and enhancers in cranial neural crest cells (CNCCs) and repressed regions in induced pluripotent cells (iPSCs), we found that 12% of indels overlap putative regulatory and repressed regions (RRRs), and 15% of these indels are associated with lineage-biased RRRs. Indel-associated putative enhancer and repressive regions are approximately 1.3 times and approximately three times as likely to be lineage-biased, respectively, as those not associated with indels. We found a twofold enrichment of medium-sized indels (20-50 bp) in CpG island (CGI)-containing promoters than expected by chance. Lastly, from human-specific transposable element insertions, we identified putative regulatory elements, including NR2F1-bound putative CNCC enhancers derived from SVAs and putative iPSC promoters derived from LTR5s. Our results show that different types of indels are associated with specific epigenomic diversity between human and chimpanzee.

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.

Integrated Analysis of Tissue-Specific Promoter Methylation and Gene Expression Profile in Complex Diseases

This study investigated whether the promoter region of DNA methylation positively or negatively regulates tissue-specific genes (TSGs) and if it correlates with disease pathophysiology. We assessed tissue specificity metrics in five human tissues, using sequencing-based approaches, including 52 whole genome bisulfite sequencing (WGBS), 52 RNA-seq, and 144 chromatin immunoprecipitation sequencing (ChIP-seq) data. A correlation analysis was performed between the gene expression and DNA methylation levels of the TSG promoter region. The TSG enrichment analyses were conducted in the gene-disease association network (DisGeNET). The epigenomic association analyses of CpGs in enriched TSG promoters were performed using 1986 Infinium MethylationEPIC array data. A correlation analysis showed significant associations between the promoter methylation and 449 TSGs' expression. A disease enrichment analysis showed that diabetes- and obesity-related diseases were high-ranked. In an epigenomic association analysis based on obesity, 62 CpGs showed statistical significance. Among them, three obesity-related CpGs were newly identified and replicated with statistical significance in independent data. In particular, a CpG (cg17075888 of PDK4), considered as potential therapeutic targets, were associated with complex diseases, including obesity and type 2 diabetes. The methylation changes in a substantial number of the TSG promoters showed a significant association with metabolic diseases. Collectively, our findings provided strong evidence of the relationship between tissue-specific patterns of epigenetic changes and metabolic diseases.

Harnessing epigenetics to study human evolution

Recent advances in ancient DNA extraction and high-throughput sequencing technologies enabled the high-quality sequencing of archaic genomes, including the Neanderthal and the Denisovan. While comparisons with modern humans revealed both archaic-specific and human-specific sequence changes, in the absence of gene expression information, understanding the functional implications of such genetic variations remains a major challenge. To study gene regulation in archaic humans, epigenetic research comes to our aid. DNA methylation, which is highly correlated with transcription, can be directly measured in modern samples, as well as reconstructed in ancient samples. This puts DNA methylation as a natural basis for comparative epigenetics between modern humans, archaic humans and nonhuman primates.

TypeTE: a tool to genotype mobile element insertions from whole genome resequencing data

Alu retrotransposons account for more than 10% of the human genome, and insertions of these elements create structural variants segregating in human populations. Such polymorphic Alus are powerful markers to understand population structure, and they represent variants that can greatly impact genome function, including gene expression. Accurate genotyping of Alus and other mobile elements has been challenging. Indeed, we found that Alu genotypes previously called for the 1000 Genomes Project are sometimes erroneous, which poses significant problems for phasing these insertions with other variants that comprise the haplotype. To ameliorate this issue, we introduce a new pipeline - TypeTE - which genotypes Alu insertions from whole-genome sequencing data. Starting from a list of polymorphic Alus, TypeTE identifies the hallmarks (poly-A tail and target site duplication) and orientation of Alu insertions using local re-assembly to reconstruct presence and absence alleles. Genotype likelihoods are then computed after re-mapping sequencing reads to the reconstructed alleles. Using a high-quality set of PCR-based genotyping of >200 loci, we show that TypeTE improves genotype accuracy from 83% to 92% in the 1000 Genomes dataset. TypeTE can be readily adapted to other retrotransposon families and brings a valuable toolbox addition for population genomics.

Metabolic changes in human brain evolution

Because the human brain is considerably larger than those of other primates, it is not surprising that its energy requirements would far exceed that of any of the species within the order. Recently, the development of stem cell technologies and single-cell transcriptomics provides novel ways to address the question of what specific genomic changes underlie the human brain's unique phenotype. In this review, we consider what is currently known about human brain metabolism using a variety of methods from brain imaging and stereology to transcriptomics. Next, we examine novel opportunities that stem cell technologies and single-cell transcriptomics provide to further our knowledge of human brain energetics. These new experimental approaches provide the ability to elucidate the functional effects of changes in genetic sequence and expression levels that potentially had a profound impact on the evolution of the human brain.

Intraspecific and interspecific investigations of skeletal DNA methylation and femur morphology in primates

OBJECTIVES

Epigenetic mechanisms influence the development and maintenance of complex phenotypes and may also contribute to the evolution of species-specific phenotypes. With respect to skeletal traits, little is known about the gene regulation underlying these hard tissues or how tissue-specific patterns are associated with bone morphology or vary among species. To begin exploring these topics, this study evaluates one epigenetic mechanism, DNA methylation, in skeletal tissues from five nonhuman primate species which display anatomical and locomotor differences representative of their phylogenetic groups.

MATERIALS AND METHODS

First, we test whether intraspecific variation in skeletal DNA methylation is associated with intraspecific variation in femur morphology. Second, we identify interspecific differences in DNA methylation and assess whether these lineage-specific patterns may have contributed to species-specific morphologies. Specifically, we use the Illumina Infinium MethylationEPIC BeadChip to identify DNA methylation patterns in femur trabecular bone from baboons (n = 28), macaques (n = 10), vervets (n = 10), chimpanzees (n = 4), and marmosets (n = 6).

RESULTS

Significant differentially methylated positions (DMPs) were associated with a subset of morphological variants, but these likely have small biological effects and may be confounded by other variables associated with morphological variation. Conversely, several species-specific DMPs were identified, and these are found in genes enriched for functions associated with complex skeletal traits.

DISCUSSION

Overall, these findings reveal that while intraspecific epigenetic variation is not readily associated with skeletal morphology differences, some interspecific epigenetic differences in skeletal tissues exist and may contribute to evolutionarily distinct phenotypes. This work forms a foundation for future explorations of gene regulation and skeletal trait evolution in primates.

Dynamic Signatures of the Epigenome: Friend or Foe?

Highly dynamic epigenetic signaling is influenced mainly by (micro)environmental stimuli and genetic factors. The exact mechanisms affecting particular epigenomic patterns differ dependently on the context. In the current review, we focus on the causes and effects of the dynamic signatures of the human epigenome as evaluated with the high-throughput profiling data and single-gene approaches. We will discuss three different aspects of phenotypic outcomes occurring as a consequence of epigenetics interplaying with genotype and environment. The first issue is related to the cases of environmental impacts on epigenetic profile, and its adverse and advantageous effects related to human health and evolutionary adaptation. The next topic will present a model of the interwoven co-evolution of genetic and epigenetic patterns exemplified with transposable elements (TEs) and their epigenetic repressors Krüppel-associated box zinc finger proteins (KRAB–ZNFs). The third aspect concentrates on the mitosis-based microevolution that takes place during carcinogenesis, leading to clonal diversity and expansion of tumor cells. The whole picture of epigenome plasticity and its role in distinct biological processes is still incomplete. However, accumulating data define epigenomic dynamics as an essential co-factor driving adaptation at the cellular and inter-species levels with a benefit or disadvantage to the host.

Comparative Molecular Life History of Spontaneous Canine and Human Gliomas

Sporadic gliomas in companion dogs provide a window on the interaction between tumorigenic mechanisms and host environment. We compared the molecular profiles of canine gliomas with those of human pediatric and adult gliomas to characterize evolutionarily conserved mammalian mutational processes in gliomagenesis. Employing whole-genome, exome, transcriptome, and methylation sequencing of 83 canine gliomas, we found alterations shared between canine and human gliomas such as the receptor tyrosine kinases, TP53 and cell-cycle pathways, and IDH1 R132. Canine gliomas showed high similarity with human pediatric gliomas per robust aneuploidy, mutational rates, relative timing of mutations, and DNA-methylation patterns. Our cross-species comparative genomic analysis provides unique insights into glioma etiology and the chronology of glioma-causing somatic alterations.

A comparison of gene expression and DNA methylation patterns across tissues and species

Previously published comparative functional genomic data sets from primates using frozen tissue samples, including many data sets from our own group, were often collected and analyzed using nonoptimal study designs and analysis approaches. In addition, when samples from multiple tissues were studied in a comparative framework, individuals and tissues were confounded. We designed a multitissue comparative study of gene expression and DNA methylation in primates that minimizes confounding effects by using a balanced design with respect to species, tissues, and individuals. We also developed a comparative analysis pipeline that minimizes biases attributable to sequence divergence. Thus, we present the most comprehensive catalog of similarities and differences in gene expression and DNA methylation levels between livers, kidneys, hearts, and lungs, in humans, chimpanzees, and rhesus macaques. We estimate that overall, interspecies and inter-tissue differences in gene expression levels can only modestly be accounted for by corresponding differences in promoter DNA methylation. However, the expression pattern of genes with conserved inter-tissue expression differences can be explained by corresponding interspecies methylation changes more often. Finally, we show that genes whose tissue-specific regulatory patterns are consistent with the action of natural selection are highly connected in both gene regulatory and protein–protein interaction networks.

IMAGE: high-powered detection of genetic effects on DNA methylation using integrated methylation QTL mapping and allele-specific analysis

Identifying genetic variants that are associated with methylation variation-an analysis commonly referred to as methylation quantitative trait locus (mQTL) mapping-is important for understanding the epigenetic mechanisms underlying genotype-trait associations. Here, we develop a statistical method, IMAGE, for mQTL mapping in sequencing-based methylation studies. IMAGE properly accounts for the count nature of bisulfite sequencing data and incorporates allele-specific methylation patterns from heterozygous individuals to enable more powerful mQTL discovery. We compare IMAGE with existing approaches through extensive simulation. We also apply IMAGE to analyze two bisulfite sequencing studies, in which IMAGE identifies more mQTL than existing approaches.

Evaluating Chromatin Accessibility Differences Across Multiple Primate Species Using a Joint Modeling Approach

Changes in transcriptional regulation are thought to be a major contributor to the evolution of phenotypic traits, but the contribution of changes in chromatin accessibility to the evolution of gene expression remains almost entirely unknown. To address this important gap in knowledge, we developed a new method to identify DNase I Hypersensitive (DHS) sites with differential chromatin accessibility between species using a joint modeling approach. Our method overcomes several limitations inherent to conventional threshold-based pairwise comparisons that become increasingly apparent as the number of species analyzed rises. Our approach employs a single quantitative test which is more sensitive than existing pairwise methods. To illustrate, we applied our joint approach to DHS sites in fibroblast cells from five primates (human, chimpanzee, gorilla, orangutan, and rhesus macaque). We identified 89,744 DHS sites, of which 41% are identified as differential between species using the joint model compared with 33% using the conventional pairwise approach. The joint model provides a principled approach to distinguishing single from multiple chromatin accessibility changes among species. We found that nondifferential DHS sites are enriched for nucleotide conservation. Differential DHS sites with decreased chromatin accessibility relative to rhesus macaque occur more commonly near transcription start sites (TSS), while those with increased chromatin accessibility occur more commonly distal to TSS. Further, differential DHS sites near TSS are less cell type-specific than more distal regulatory elements. Taken together, these results point to distinct classes of DHS sites, each with distinct characteristics of selection, genomic location, and cell type specificity.

Reorganization of 3D genome structure may contribute to gene regulatory evolution in primates

A growing body of evidence supports the notion that variation in gene regulation plays a crucial role in both speciation and adaptation. However, a comprehensive functional understanding of the mechanisms underlying regulatory evolution remains elusive. In primates, one of the crucial missing pieces of information towards a better understanding of regulatory evolution is a comparative annotation of interactions between distal regulatory elements and promoters. Chromatin conformation capture technologies have enabled genome-wide quantifications of such distal 3D interactions. However, relatively little comparative research in primates has been done using such technologies. To address this gap, we used Hi-C to characterize 3D chromatin interactions in induced pluripotent stem cells (iPSCs) from humans and chimpanzees. We also used RNA-seq to collect gene expression data from the same lines. We generally observed that lower-order, pairwise 3D genomic interactions are conserved in humans and chimpanzees, but higher order genomic structures, such as topologically associating domains (TADs), are not as conserved. Inter-species differences in 3D genomic interactions are often associated with gene expression differences between the species. To provide additional functional context to our observations, we considered previously published chromatin data from human stem cells. We found that inter-species differences in 3D genomic interactions, which are also associated with gene expression differences between the species, are enriched for both active and repressive marks. Overall, our data demonstrate that, as expected, an understanding of 3D genome reorganization is key to explaining regulatory evolution.

The way in which a genome folds affects the regulation of gene expression. This is often due to loops in the three-dimensional structure that bring linearly distant genes and regulatory elements into close proximity. Most studies examining three-dimensional structure genome-wide are limited to a single species. In this study, we compared three-dimensional structure in the genomes of induced pluripotent stem cells from humans and chimpanzees. We collected gene expression data from the same samples, which allowed us to assess the contribution of three-dimensional chromatin conformation to gene regulatory evolution in primates. Our results demonstrate that gene expression differences between the species may often be mediated by differences in three-dimensional genomic interactions. Our data also suggest that large-scale chromatin structures (i.e. topologically associating domains, TADs) are not well conserved in their placement across species. We hope the analytical paradigms we present here could serve as a basis for future comparative studies of three-dimensional genome organization, elucidating the putative functional regulatory loci driving speciation.

Assessment of DNA Methylation Patterns in the Bone and Cartilage of a Nonhuman Primate Model of Osteoarthritis

OBJECTIVE

Osteoarthritis (OA) affects humans and several other animals. Thus, the mechanisms underlying this disorder, such as specific skeletal tissue DNA methylation patterns, may be evolutionary conserved. However, associations between methylation and OA have not been readily studied in nonhuman animals. Baboons serve as important models of disease and develop OA at rates similar to those in humans. Therefore, this study investigated the associations between methylation and OA in baboons to advance the evolutionary understanding of OA.

DESIGN

Trabecular bone and cartilage was collected from the medial condyles of adult female baboon femora, 5 with and 5 without knee OA. The Infinium HumanMethylation450 BeadChip (450K array) was used to identify DNA methylation patterns in these tissues.

RESULTS

Approximately 44% of the 450K array probes reliably align to the baboon genome, contain a CpG site of interest, and maintain a wide distribution throughout the genome. Of the 2 filtering methods tested, both identified significantly differentially methylated positions (DMPs) between healthy and OA individuals in cartilage tissues, and some of these patterns overlap with those previously identified in humans. Conversely, no DMPs were found between tissue types or between disease states in bone tissues.

CONCLUSIONS

Overall, the 450K array can be used to measure genome-wide DNA methylation in baboon tissues and identify significant associations with complex traits. The results of this study indicate that some DNA methylation patterns associated with OA are evolutionarily conserved, while others are not. This warrants further investigation in a larger and more phylogenetically diverse sample set.

Oxytocin alleviates cellular senescence through oxytocin receptor-mediated extracellular signal-regulated kinase/Nrf2 signalling

BACKGROUND

Oxytocin (OT) is a neuropeptide hormone that has many beneficial biological effects, including protection against age-related disorders. However, less is known about its role in intrinsic skin ageing, which is accelerated by an increase in senescent cell fraction in skin tissue.

OBJECTIVES

To investigate the novel function and the underlying mechanism of OT in preventing cellular senescence in normal human dermal fibroblasts (NHDFs) isolated from the skin of female donors of different ages.

METHODS

NHDFs from young and old donors were exposed to conditioned medium from senescent or control NHDFs in the presence or absence of 10 nmol L^-1 OT for 3 days, and were continuously subcultured for 12 days. Subsequently, various age-associated signs of senescence including decreased proliferation rate, elevated p16 and p21 levels, and positivity for senescence-associated β-galactosidase expression were examined.

RESULTS

We found that OT suppressed senescence-associated secretory phenotype-induced senescence in NHDFs, and its effect depended on the age of the donor's NHDFs. The inhibitory effects of OT required signalling by OT receptor-mediated extracellular signal-regulated kinase/Nrf2 (nuclear factor erythroid 2-related factor 2). The age-dependent antisenescence effects of OT are closely related to hypermethylation of the OT receptor gene (OXTR).

CONCLUSIONS

Our findings bring to light the role of OT in the prevention of skin ageing, which might allow development of new clinical strategies. What's already known about this topic? Senescent keratinocytes and fibroblasts accumulate with age in the skin and contribute to the loss of skin function and integrity during ageing. Senescent cells secrete senescence-associated secretory phenotype (SASP), which includes the release of proinflammatory cytokines such as interleukin (IL)-6 and IL-1, chemokines, extracellular matrix-remodelling proteases and growth factors. The neuropeptide oxytocin (OT) and its receptor (OXTR) have protective effects against various age-related disorders. What does this study add? OT suppressed SASP-induced cellular senescence in normal human dermal fibroblasts (NHDFs), depending on the age of the NHDFs' donor. The inhibitory effects of OT on cellular senescence required OXTR-mediated phosphorylation of extracellular signal-regulated kinase, which enhanced nuclear localization of Nrf2, a vital factor in the antioxidant defence system. The age-specific antisenescent effects of OT were closely related to hypermethylation of OXTR. What is the translational message? Our results suggest that OT and OXTR agonists could be clinically promising agents for the improvement of age-associated skin ageing, especially in women.

Evolution of DNA Methylation in Papio Baboons

Changes in gene regulation have long been thought to play an important role in primate evolution. However, although a number of studies have compared genome-wide gene expression patterns across primate species, fewer have investigated the gene regulatory mechanisms that underlie such patterns, or the relative contribution of drift versus selection. Here, we profiled genome-scale DNA methylation levels in blood samples from five of the six extant species of the baboon genus Papio (4-14 individuals per species). This radiation presents the opportunity to investigate DNA methylation divergence at both shallow and deeper timescales (0.380-1.4 My). In contrast to studies in human populations, but similar to studies in great apes, DNA methylation profiles clearly mirror genetic and geographic structure. Divergence in DNA methylation proceeds fastest in unannotated regions of the genome and slowest in regions of the genome that are likely more constrained at the sequence level (e.g., gene exons). Both heuristic approaches and Ornstein-Uhlenbeck models suggest that DNA methylation levels at a small set of sites have been affected by positive selection, and that this class is enriched in functionally relevant contexts, including promoters, enhancers, and CpG islands. Our results thus indicate that the rate and distribution of DNA methylation changes across the genome largely mirror genetic structure. However, at some CpG sites, DNA methylation levels themselves may have been a target of positive selection, pointing to loci that could be important in connecting sequence variation to fitness-related traits.

Multiple Innovations in Genetic and Epigenetic Mechanisms Cooperate to Underpin Human Brain Evolution

Our knowledge of how the human brain differs from those of other species in terms of evolutionary adaptations and functionality is limited. Comparative genomics reveal valuable insight, especially the expansion of human-specific noncoding regulatory and repeat-containing regions. Recent studies add to our knowledge of evolving brain function by investigating cellular mechanisms such as protein emergence, extensive sequence editing, retrotransposon activity, dynamic epigenetic modifications, and multiple noncoding RNA functions. These findings present an opportunity to combine newly discovered genetic and epigenetic mechanisms with more established concepts into a more comprehensive picture to better understand the uniquely evolved human brain.

Exploring the genetic basis of human population differences in DNA methylation and their causal impact on immune gene regulation

BACKGROUND

DNA methylation is influenced by both environmental and genetic factors and is increasingly thought to affect variation in complex traits and diseases. Yet, the extent of ancestry-related differences in DNA methylation, their genetic determinants, and their respective causal impact on immune gene regulation remain elusive.

RESULTS

We report extensive population differences in DNA methylation between 156 individuals of African and European descent, detected in primary monocytes that are used as a model of a major innate immunity cell type. Most of these differences (~ 70%) are driven by DNA sequence variants nearby CpG sites, which account for ~ 60% of the variance in DNA methylation. We also identify several master regulators of DNA methylation variation in trans, including a regulatory hub nearby the transcription factor-encoding CTCF gene, which contributes markedly to ancestry-related differences in DNA methylation. Furthermore, we establish that variation in DNA methylation is associated with varying gene expression levels following mostly, but not exclusively, a canonical model of negative associations, particularly in enhancer regions. Specifically, we find that DNA methylation highly correlates with transcriptional activity of 811 and 230 genes, at the basal state and upon immune stimulation, respectively. Finally, using a Bayesian approach, we estimate causal mediation effects of DNA methylation on gene expression in ~ 20% of the studied cases, indicating that DNA methylation can play an active role in immune gene regulation.

CONCLUSION

Using a system-level approach, our study reveals substantial ancestry-related differences in DNA methylation and provides evidence for their causal impact on immune gene regulation.

Modeling three-dimensional genomic organization in evolution and pathogenesis

The regulation of gene expression is mediated via the complex three-dimensional (3D) conformation of the genetic material and its interactions with various intracellular factors. Various experimental and computational approaches have been developed in recent years for understating the relation between the 3D conformation of the genome and the phenotypes of cells in normal condition and diseases. In this review, we will discuss novel approaches for analyzing and modeling the 3D genomic conformation, focusing on deciphering disease-causing mutations that affect gene expression. We conclude that as this is a very challenging mission, an important direction should involve the comparative analysis of various 3D models from various organisms or cells.

Towards a map of cis-regulatory sequences in the human genome

Accumulating evidence indicates that transcription factor (TF) binding sites, or cis-regulatory elements (CREs), and their clusters termed cis-regulatory modules (CRMs) play a more important role than do gene-coding sequences in specifying complex traits in humans, including the susceptibility to common complex diseases. To fully characterize their roles in deriving the complex traits/diseases, it is necessary to annotate all CREs and CRMs encoded in the human genome. However, the current annotations of CREs and CRMs in the human genome are still very limited and mostly coarse-grained, as they often lack the detailed information of CREs in CRMs. Here, we integrated 620 TF ChIP-seq datasets produced by the ENCODE project for 168 TFs in 79 different cell/tissue types and predicted an unprecedentedly completely map of CREs in CRMs in the human genome at single nucleotide resolution. The map includes 305 912 CRMs containing a total of 1 178 913 CREs belonging to 736 unique TF binding motifs. The predicted CREs and CRMs tend to be subject to either purifying selection or positive selection, thus are likely to be functional. Based on the results, we also examined the status of available ChIP-seq datasets for predicting the entire regulatory genome of humans.

TOMM40 and APOE Gene Expression and Cognitive Decline in Japanese Alzheimer's Disease Subjects

BACKGROUND

TOMM40 is located on chromosome 19, is in linkage disequilibrium with apolipoprotein E (APOE), andis reported in several genome-wide association studies to be associated with Alzheimer's disease (AD).

OBJECTIVE

Assess APOE and TOM40 and mitochondrial genes as blood biomarkers for AD.

METHODS

We examined TOMM40, PTEN-induced putative kinase 1 (PINK1), Parkin RBR E3 ubiquitin protein ligase (PARK2), and APOE mRNA expression in relation to the methylation rates of CpG sites in the upstream region of TOMM40exon 1 in peripheral leukocytes and TOMM40523 polyT genotypes in 60 AD and age- and sex-matched control subjects.

RESULTS

TOMM40 mRNA expression was significantly lower in AD subjects (0.87±0.18 versus 1.0±0.23, p = 0.005), and PINK1 mRNA expression was higher in AD subjects (1.5±0.61 versus 1.0±0.52, p < 0.001). TOMM40 mRNA expression was significantly correlated with the Mini-Mental State Examination total score (r = 0.290, p = 0.027). There was no expressional change in peripheral APOE mRNA in either AD or control subjects (p = 0.32). Methylation rates in the upstream region of TOMM40exon 1 were not different between AD and control subjects (average rate: 1.37±0.99 versus 1.39±1.20, p = 0.885), and TOMM40523 polyT genotypes were also not different between AD and control subjects (p = 0.67).

CONCLUSION

TOMM40 mRNA expression was lower in AD subjects and was correlated with cognitive decline. Significant changes in both TOMM40 and PINK1 mRNA may be related to mitochondrial dysfunction.

Of Men and Mice: Modeling the Fragile X Syndrome

The Fragile X Syndrome (FXS) is one of the most common forms of inherited intellectual disability in all human societies. Caused by the transcriptional silencing of a single gene, the fragile x mental retardation gene FMR1, FXS is characterized by a variety of symptoms, which range from mental disabilities to autism and epilepsy. More than 20 years ago, a first animal model was described, the Fmr1 knock-out mouse. Several other models have been developed since then, including conditional knock-out mice, knock-out rats, a zebrafish and a drosophila model. Using these model systems, various targets for potential pharmaceutical treatments have been identified and many treatments have been shown to be efficient in preclinical studies. However, all attempts to turn these findings into a therapy for patients have failed thus far. In this review, I will discuss underlying difficulties and address potential alternatives for our future research.

Ageing-associated DNA methylation dynamics are a molecular readout of lifespan variation among mammalian species

BACKGROUND

Mammalian species exhibit a wide range of lifespans. To date, a robust and dynamic molecular readout of these lifespan differences has not yet been identified. Recent studies have established the existence of ageing-associated differentially methylated positions (aDMPs) in human and mouse. These are CpG sites at which DNA methylation dynamics show significant correlations with age. We hypothesise that aDMPs are pan-mammalian and are a dynamic molecular readout of lifespan variation among different mammalian species.

RESULTS

A large-scale integrated analysis of aDMPs in six different mammals reveals a strong negative relationship between rate of change of methylation levels at aDMPs and lifespan. This relationship also holds when comparing two different dog breeds with known differences in lifespans. In an ageing cohort of aneuploid mice carrying a complete copy of human chromosome 21, aDMPs accumulate far more rapidly than is seen in human tissues, revealing that DNA methylation at aDMP sites is largely shaped by the nuclear trans-environment and represents a robust molecular readout of the ageing cellular milieu.

CONCLUSIONS

Overall, we define the first dynamic molecular readout of lifespan differences among mammalian species and propose that aDMPs will be an invaluable molecular tool for future evolutionary and mechanistic studies aimed at understanding the biological factors that determine lifespan in mammals.

Evolutionary expansion of DNA hypomethylation in the mammalian germline genome

DNA methylation in the germline is among the most important factors influencing the evolution of mammalian genomes. Yet little is known about its evolutionary rate or the fraction of the methylome that has undergone change. We compared whole-genome, single-CpG DNA methylation profiles in sperm of seven species-human, chimpanzee, gorilla, rhesus macaque, mouse, rat, and dog-to investigate epigenomic evolution. We developed a phylo-epigenetic model for DNA methylation that accommodates the correlation of states at neighboring sites and allows for inference of ancestral states. Applying this model to the sperm methylomes, we uncovered an overall evolutionary expansion of the hypomethylated fraction of the genome, driven both by the birth of new hypomethylated regions and by extensive widening of hypomethylated intervals in ancestral species. This expansion shows strong lineage-specific aspects, most notably that hypomethylated intervals around transcription start sites have evolved to be considerably wider in primates and dog than in rodents, whereas rodents show evidence of a greater trend toward birth of new hypomethylated regions. Lineage-specific hypomethylated regions are enriched near sets of genes with common developmental functions and significant overlap across lineages. Rodent-specific and primate-specific hypomethylated regions are enriched for binding sites of similar transcription factors, suggesting that the plasticity accommodated by certain regulatory factors is conserved, despite substantial change in the specific sites of regulation. Overall our results reveal substantial global epigenomic change in mammalian sperm methylomes and point to a divergence in trans-epigenetic mechanisms that govern the organization of epigenetic states at gene promoters.

Elevated mRNA Expression and Low Methylation of SNCA in Japanese Alzheimer's Disease Subjects

Despite the continuing debate about the amyloid hypothesis in Alzheimer's disease (AD), the precise pathogenesis is still unclear. Mixed pathology is common and multiple different protein aggregates are seen in human postmortem brains. Aggregates consisting of the alpha-synuclein protein encoded by the Synuclein Alpha gene (SCNA) are common in both dementia with Lewy bodies and AD. We examined SNCA mRNA expression and methylation rates of the CpG island at intron 1 of SNCA in peripheral leukocytes in 50 AD and age- and sex-matched control subjects to verify whether alpha-synuclein pathology affects the AD pathogenesis. SNCA mRNA expression in AD subjects was significantly higher than that in control subjects (1.62±0.73 versus 0.98±0.50, p < 0.001). We found significant differences between AD and control subjects at seven CpG sites (average rate; 8.8±2.7 versus 9.5±2.5, respectively: p = 0.027). The methylation rates tended to be lower in AD subjects at all CpG sites. We conclude that mRNA expression and methylation of SNCA intron 1 are altered in AD, which may be caused by Lewy body pathology in AD.

Complexity and conservation of regulatory landscapes underlie evolutionary resilience of mammalian gene expression

To gain insight into how mammalian gene expression is controlled by rapidly evolving regulatory elements, we jointly analysed promoter and enhancer activity with downstream transcription levels in liver samples from 15 species. Genes associated with complex regulatory landscapes generally exhibit high expression levels that remain evolutionarily stable. While the number of regulatory elements is the key driver of transcriptional output and resilience, regulatory conservation matters: elements active across mammals most effectively stabilize gene expression. In contrast, recently evolved enhancers typically contribute weakly, consistent with their high evolutionary plasticity. These effects are observed across the entire mammalian clade and are robust to potential confounders, such as the gene expression level. Using liver as a representative somatic tissue, our results illuminate how the evolutionary stability of gene expression is profoundly entwined with both the number and conservation of surrounding promoters and enhancers.

DNA Methylation Divergence and Tissue Specialization in the Developing Mouse Placenta

The placental epigenome plays a vital role in regulating mammalian growth and development. Aberrations in placental DNA methylation are linked to several disease states, including intrauterine growth restriction and preeclampsia. Studying the evolution and development of the placental epigenome is critical to understanding the origin and progression of such diseases. Although high-resolution studies have found substantial variation between placental methylomes of different species, the nature of methylome variation has yet to be characterized within any individual species. We conducted a study of placental DNA methylation at high resolution in multiple strains and closely related species of house mice (Mus musculus musculus, Mus m. domesticus, and M. spretus), across developmental timepoints (embryonic days 15-18), and between two distinct layers (labyrinthine transport and junctional endocrine). We observed substantial genome-wide methylation heterogeneity in mouse placenta compared with other differentiated tissues. Species-specific methylation profiles were concentrated in retrotransposon subfamilies, specifically RLTR10 and RLTR20 subfamilies. Regulatory regions such as gene promoters and CpG islands displayed cross-species conservation, but showed strong differences between layers and developmental timepoints. Partially methylated domains exist in the mouse placenta and widen during development. Taken together, our results characterize the mouse placental methylome as a highly heterogeneous and deregulated landscape globally, intermixed with actively regulated promoter and retrotransposon sequences.