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.

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.

A supergene-linked estrogen receptor drives alternative phenotypes in a polymorphic songbird.

Behavioral evolution relies on genetic changes, yet few behaviors can be traced to specific genetic sequences in vertebrates. Here we provide experimental evidence showing that differentiation of a single gene has contributed to the evolution of divergent behavioral phenotypes in the white-throated sparrow, a common backyard songbird. In this species, a series of chromosomal inversions has formed a supergene that segregates with an aggressive phenotype. The supergene has captured ESR1, the gene that encodes estrogen receptor α (ERα); as a result, this gene is accumulating changes that now distinguish the supergene allele from the standard allele. Our results show that in birds of the more aggressive phenotype, ERα knockdown caused a phenotypic change to that of the less aggressive phenotype. We next showed that in a free-living population, aggression is predicted by allelic imbalance favoring the supergene allele. Finally, we identified cis-regulatory features, both genetic and epigenetic, that explain the allelic imbalance. This work provides a rare illustration of how genotypic divergence has led to behavioral phenotypic divergence in a vertebrate.

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.

Myxosporea (Myxozoa, Cnidaria) Lack DNA Cytosine Methylation.

DNA cytosine methylation is central to many biological processes, including regulation of gene expression, cellular differentiation, and development. This DNA modification is conserved across animals, having been found in representatives of sponges, ctenophores, cnidarians, and bilaterians, and with very few known instances of secondary loss in animals. Myxozoans are a group of microscopic, obligate endoparasitic cnidarians that have lost many genes over the course of their evolution from free-living ancestors. Here, we investigated the evolution of the key enzymes involved in DNA cytosine methylation in 29 cnidarians and found that these enzymes were lost in an ancestor of Myxosporea (the most speciose class of Myxozoa). Additionally, using whole-genome bisulfite sequencing, we confirmed that the genomes of two distant species of myxosporeans, Ceratonova shasta and Henneguya salminicola, completely lack DNA cytosine methylation. Our results add a notable and novel taxonomic group, the Myxosporea, to the very short list of animal taxa lacking DNA cytosine methylation, further illuminating the complex evolutionary history of this epigenetic regulatory mechanism.

Regional epigenetic differentiation of the Z Chromosome between sexes in a female heterogametic system.

In male heterogametic systems, the X Chromosome is epigenetically differentiated between males and females, to facilitate dosage compensation. For example, the X Chromosome in female mammals is largely inactivated. Relative to well-studied male heterogametic systems, the extent of epigenetic differentiation between male and female Z Chromosomes in female heterogametic species, which often lack complete dosage compensation, is poorly understood. Here, we examined the chromosomal DNA methylation landscapes of male and female Z Chromosomes in two distantly related avian species, namely chicken and white-throated sparrow. We show that, in contrast to the pattern in mammals, male and female Z Chromosomes in these species exhibit highly similar patterns of DNA methylation, which is consistent with weak or absent dosage compensation. We further demonstrate that the epigenetic differences between male and female chicken Z Chromosomes are localized to a few regions, including a previously identified male hypermethylated region 1 (MHM1; CGNC: 80601). We discovered a novel region with elevated male-to-female methylation ratios on the chicken Z Chromosome (male hypermethylated region 2 [MHM2]; CGNC: 80602). The MHM1 and MHM2, despite little sequence similarity between them, bear similar molecular features that are likely associated with their functions. We present evidence consistent with female hypomethylation of MHMs and up-regulation of nearby genes. Therefore, despite little methylation differentiation between sexes, extremely localized DNA methylation differences between male and female chicken Z Chromosomes have evolved and affect expression of nearby regions. Our findings offer new insights into epigenetic regulation of gene expression between sexes in female heterogametic systems.

Distinct epigenomic and transcriptomic modifications associated with Wolbachia-mediated asexuality.

Wolbachia are maternally transmitted intracellular bacteria that induce a range of pathogenic and fitness-altering effects on insect and nematode hosts. In parasitoid wasps of the genus Trichogramma, Wolbachia infection induces asexual production of females, thus increasing transmission of Wolbachia. It has been hypothesized that Wolbachia infection accompanies a modification of the host epigenome. However, to date, data on genome-wide epigenomic changes associated with Wolbachia are limited, and are often confounded by background genetic differences. Here, we took sexually reproducing Trichogramma free of Wolbachia and introgressed their genome into a Wolbachia-infected cytoplasm, converting them to Wolbachia-mediated asexuality. Wolbachia was then cured from replicates of these introgressed lines, allowing us to examine the genome-wide effects of wasps newly converted to asexual reproduction while controlling for genetic background. We thus identified gene expression and DNA methylation changes associated with Wolbachia-infection. We found no overlaps between differentially expressed genes and differentially methylated genes, indicating that Wolbachia-infection associated DNA methylation change does not directly modulate levels of gene expression. Furthermore, genes affected by these mechanisms exhibit distinct evolutionary histories. Genes differentially methylated due to the infection tended to be evolutionarily conserved. In contrast, differentially expressed genes were significantly more likely to be unique to the Trichogramma lineage, suggesting host-specific transcriptomic responses to infection. Nevertheless, we identified several novel aspects of Wolbachia-associated DNA methylation changes. Differentially methylated genes included those involved in oocyte development and chromosome segregation. Interestingly, Wolbachia-infection was associated with higher levels of DNA methylation. Additionally, Wolbachia infection reduced overall variability in gene expression, even after accounting for the effect of DNA methylation. We also identified specific cases where alternative exon usage was associated with DNA methylation changes due to Wolbachia infection. These results begin to reveal distinct genes and molecular pathways subject to Wolbachia induced epigenetic modification and/or host responses to Wolbachia-infection.

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.

Accelerated evolution of oligodendrocytes in the human brain.

Recent discussions of human brain evolution have largely focused on increased neuron numbers and changes in their connectivity and expression. However, it is increasingly appreciated that oligodendrocytes play important roles in cognitive function and disease. Whether both cell types follow similar or distinctive evolutionary trajectories is not known. We examined the transcriptomes of neurons and oligodendrocytes in the frontal cortex of humans, chimpanzees, and rhesus macaques. We identified human-specific trajectories of gene expression in neurons and oligodendrocytes and show that both cell types exhibit human-specific up-regulation. Moreover, oligodendrocytes have undergone more pronounced accelerated gene expression evolution in the human lineage compared to neurons. We highlighted human-specific coexpression networks with specific functions. Our data suggest that oligodendrocyte human-specific networks are enriched for alternative splicing and transcriptional regulation. Oligodendrocyte networks are also enriched for variants associated with schizophrenia and other neuropsychiatric disorders. Such enrichments were not found in neuronal networks. These results offer a glimpse into the molecular mechanisms of oligodendrocytes during evolution and how such mechanisms are associated with neuropsychiatric disorders.

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.

Detecting differential DNA methylation from sequencing of bisulfite converted DNA of diverse species.

DNA methylation is one of the most extensively studied epigenetic modifications of genomic DNA. In recent years, sequencing of bisulfite-converted DNA, particularly via next-generation sequencing technologies, has become a widely popular method to study DNA methylation. This method can be readily applied to a variety of species, dramatically expanding the scope of DNA methylation studies beyond the traditionally studied human and mouse systems. In parallel to the increasing wealth of genomic methylation profiles, many statistical tools have been developed to detect differentially methylated loci (DMLs) or differentially methylated regions (DMRs) between biological conditions. We discuss and summarize several key properties of currently available tools to detect DMLs and DMRs from sequencing of bisulfite-converted DNA. However, the majority of the statistical tools developed for DML/DMR analyses have been validated using only mammalian data sets, and less priority has been placed on the analyses of invertebrate or plant DNA methylation data. We demonstrate that genomic methylation profiles of non-mammalian species are often highly distinct from those of mammalian species using examples of honey bees and humans. We then discuss how such differences in data properties may affect statistical analyses. Based on these differences, we provide three specific recommendations to improve the power and accuracy of DML and DMR analyses of invertebrate data when using currently available statistical tools. These considerations should facilitate systematic and robust analyses of DNA methylation from diverse species, thus advancing our understanding of DNA methylation.

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.

Rapid regulatory evolution of a nonrecombining autosome linked to divergent behavioral phenotypes.

In the white-throated sparrow (Zonotrichia albicollis), the second chromosome bears a striking resemblance to sex chromosomes. First, within each breeding pair of birds, one bird is homozygous for the standard arrangement of the chromosome (ZAL2/ZAL2) and its mate is heterozygous for a different version (ZAL2/ZAL2m). Second, recombination is profoundly suppressed between the two versions, leading to genetic differentiation between them. Third, the ZAL2m version is linked with phenotypic traits, such as bright plumage, high aggression, and low parental behavior, which are usually associated with males. These similarities to sex chromosomes suggest that the evolutionary mechanisms that shape sex chromosomes, in particular genetic degeneration of the heterogametic version due to the suppression of recombination, are likely important in this system as well. Here, we investigated patterns of protein sequence evolution and gene expression evolution between the ZAL2 and ZAL2m chromosomes by whole-genome sequencing and transcriptome analyses. Patterns of protein evolution exhibited only weak signals of genetic degeneration, and few genes harbored signatures of positive selection. We found substantial evidence of transcriptome evolution, such as significant expression divergence between ZAL2 and ZAL2m alleles and signatures of dosage compensation for highly expressed genes. These results suggest that, early in the evolution of heteromorphic chromosomes, gene expression divergence and dosage compensation can prevail before large-scale genetic degeneration. Our results show further that suppression of recombination between heteromorphic chromosomes can lead to the evolution of alternative (sex-like) behavioral phenotypes before substantial genetic degeneration.

Selection on the regulation of sympathetic nervous activity in humans and chimpanzees.

Adrenergic α2C receptor (ADRA2C) is an inhibitory modulator of the sympathetic nervous system. Knockout mice for this gene show physiological and behavioural alterations that are associated with the fight-or-flight response. There is evidence of positive selection on the regulation of this gene during chicken domestication. Here, we find that the neuronal expression of ADRA2C is lower in human and chimpanzee than in other primates. On the basis of three-dimensional chromatin structure, we identified a cis-regulatory region whose DNA sequences have been significantly accelerated in human and chimpanzee. Active histone modification marks this region in rhesus macaque but not in human and chimpanzee; instead, repressive marks are enriched in various human brain samples. This region contains two neuron-restrictive silencer factor (NRSF) binding motifs, each of which harbours a polymorphism. Our genotyping and analysis of population genome data indicate that at both polymorphic sites, the derived allele has reached fixation in humans and chimpanzees but not in bonobos, whereas only the ancestral allele is present among macaques. Our CRISPR/Cas9 genome editing and reporter assays show that both derived nucleotides repress ADRA2C, most likely by increasing NRSF binding. In addition, we detected signatures of recent positive selection for lower neuronal ADRA2C expression in humans. Our findings indicate that there has been selective pressure for enhanced sympathetic nervous activity in the evolution of humans and chimpanzees.

Inside the supergene of the bird with four sexes.

The white-throated sparrow (Zonotrichia albicollis) offers unique opportunities to understand the adaptive value of supergenes, particularly their role in alternative phenotypes. In this species, alternative plumage morphs segregate with a nonrecombining segment of chromosome 2, which has been called a 'supergene'. The species mates disassortatively with respect to the supergene; that is, each breeding pair consists of one individual with it and one without it. This species has therefore been called the "bird with four sexes". The supergene segregates with a behavioral phenotype; birds with it are more aggressive and less parental than birds without it. Here, we review our efforts to identify the genes inside the supergene that are responsible for the behavioral polymorphism. The gene ESR1, which encodes estrogen receptor α, differs between the morphs and predicts both territorial and parental behavior. Variation in the regulatory regions of ESR1 causes an imbalance in expression of the two alleles, and the degree to which this imbalance favors the supergene allele predicts territorial singing. In heterozygotes, knockdown of ESR1 causes a phenotypic switch, from more aggressive to less aggressive. We recently showed that another gene important for social behavior, vasoactive intestinal peptide (VIP), is differentially expressed between the morphs and predicts territorial singing. We hypothesize that ESR1 and VIP contribute to behavior in a coordinated way and could represent co-adapted alleles. Because the supergene contains more than 1000 individual genes, this species provides rich possibilities for discovering alleles that work together to mediate life-history trade-offs and maximize the fitness of alternative complex phenotypes.

Functional conservation of sequence determinants at rapidly evolving regulatory regions across mammals.

Recent advances in epigenomics have made it possible to map genome-wide regulatory regions using empirical methods. Subsequent comparative epigenomic studies have revealed that regulatory regions diverge rapidly between genome of different species, and that the divergence is more pronounced in enhancers than in promoters. To understand genomic changes underlying these patterns, we investigated if we can identify specific sequence fragments that are over-enriched in regulatory regions, thus potentially contributing to regulatory functions of such regions. Here we report numerous sequence fragments that are statistically over-enriched in enhancers and promoters of different mammals (which we refer to as 'sequence determinants'). Interestingly, the degree of statistical enrichment, which presumably is associated with the degree of regulatory impacts of the specific sequence determinant, was significantly higher for promoter sequence determinants than enhancer sequence determinants. We further used a machine learning method to construct prediction models using sequence determinants. Remarkably, prediction models constructed from one species could be used to predict regulatory regions of other species with high accuracy. This observation indicates that even though the precise locations of regulatory regions diverge rapidly during evolution, the functional potential of sequence determinants underlying regulatory sequences may be conserved between species.

Body-hypomethylated human genes harbor extensive intragenic transcriptional activity and are prone to cancer-associated dysregulation.

Genomic DNA methylation maps (methylomes) encode genetic and environmental effects as stable chemical modifications of DNA. Variations in DNA methylation, especially in regulatory regions such as promoters and enhancers, are known to affect numerous downstream processes. In contrast, most transcription units (gene bodies) in the human genome are thought to be heavily methylated. However, epigenetic reprogramming in cancer often involves gene body hypomethylation with consequences on gene expression. In this study, we focus on the relatively unexplored phenomenon that some gene bodies are devoid of DNA methylation under normal conditions. Utilizing nucleotide-resolution methylomes of diverse samples, we show that nearly 2000 human genes are commonly hypomethylated. Remarkably, these genes occupy highly specialized genomic, epigenomic, evolutionary and functional niches in our genomes. For example, hypomethylated genes tend to be short yet encode significantly more transcripts than expected based upon their lengths, include many genes involved in nucleosome and chromatin formation, and are extensively and significantly enriched for histone-tail modifications and transcription factor binding with particular relevance for cis-regulation. Furthermore, they are significantly more prone to cancer-associated hypomethylation and mutation. Consequently, gene body hypomethylation represents an additional layer of epigenetic regulatory complexity, with implications on cancer-associated epigenetic reprogramming.

Comparative genomics of the miniature wasp and pest control agent Trichogramma pretiosum.

BACKGROUND: Trichogrammatids are minute parasitoid wasps that develop within other insect eggs. They are less than half a millimeter long, smaller than some protozoans. The Trichogrammatidae are one of the earliest branching families of Chalcidoidea: a diverse superfamily of approximately half a million species of parasitoid wasps, proposed to have evolved from a miniaturized ancestor. Trichogramma are frequently used in agriculture, released as biological control agents against major moth and butterfly pests. Additionally, Trichogramma are well known for their symbiotic bacteria that induce asexual reproduction in infected females. Knowledge of the genome sequence of Trichogramma is a major step towards further understanding its biology and potential applications in pest control. RESULTS: We report the 195-Mb genome sequence of Trichogramma pretiosum and uncover signatures of miniaturization and adaptation in Trichogramma and related parasitoids. Comparative analyses reveal relatively rapid evolution of proteins involved in ribosome biogenesis and function, transcriptional regulation, and ploidy regulation. Chalcids also show loss or especially rapid evolution of 285 gene clusters conserved in other Hymenoptera, including many that are involved in signal transduction and embryonic development. Comparisons between sexual and asexual lineages of Trichogramma pretiosum reveal that there is no strong evidence for genome degradation (e.g., gene loss) in the asexual lineage, although it does contain a lower repeat content than the sexual lineage. Trichogramma shows particularly rapid genome evolution compared to other hymenopterans. We speculate these changes reflect adaptations to miniaturization, and to life as a specialized egg parasitoid. CONCLUSIONS: The genomes of Trichogramma and related parasitoids are a valuable resource for future studies of these diverse and economically important insects, including explorations of parasitoid biology, symbiosis, asexuality, biological control, and the evolution of miniaturization. Understanding the molecular determinants of parasitism can also inform mass rearing of Trichogramma and other parasitoids for biological control.

Whole-genome bisulfite sequencing maps from multiple human tissues reveal novel CpG islands associated with tissue-specific regulation.

CpG islands (CGIs) are one of the most widely studied regulatory features of the human genome, with critical roles in development and disease. Despite such significance and the original epigenetic definition, currently used CGI sets are typically predicted from DNA sequence characteristics. Although CGIs are deeply implicated in practical analyses of DNA methylation, recent studies have shown that such computational annotations suffer from inaccuracies. Here we used whole-genome bisulfite sequencing from 10 diverse human tissues to identify a comprehensive, experimentally obtained, single-base resolution CGI catalog. In addition to the unparalleled annotation precision, our method is free from potential bias due to arbitrary sequence features or probe affinity differences. In addition to clarifying substantial false positives in the widely used University of California Santa Cruz (UCSC) annotations, our study identifies numerous novel epigenetic loci. In particular, we reveal significant impact of transposable elements on the epigenetic regulatory landscape of the human genome and demonstrate ubiquitous presence of transcription initiation at CGIs, including alternative promoters in gene bodies and non-coding RNAs in intergenic regions. Moreover, coordinated DNA methylation and chromatin modifications mark tissue-specific enhancers at novel CGIs. Enrichment of specific transcription factor binding from ChIP-seq supports mechanistic roles of CGIs on the regulation of tissue-specific transcription. The new CGI catalog provides a comprehensive and integrated list of genomic hotspots of epigenetic regulation.

Evolutionary Transition of Promoter and Gene Body DNA Methylation across Invertebrate-Vertebrate Boundary.

Genomes of invertebrates and vertebrates exhibit highly divergent patterns of DNA methylation. Invertebrate genomes tend to be sparsely methylated, and DNA methylation is mostly targeted to a subset of transcription units (gene bodies). In a drastic contrast, vertebrate genomes are generally globally and heavily methylated, punctuated by the limited local hypo-methylation of putative regulatory regions such as promoters. These genomic differences also translate into functional differences in DNA methylation and gene regulation. Although promoter DNA methylation is an important regulatory component of vertebrate gene expression, its role in invertebrate gene regulation has been little explored. Instead, gene body DNA methylation is associated with expression of invertebrate genes. However, the evolutionary steps leading to the differentiation of invertebrate and vertebrate genomic DNA methylation remain unresolved. Here we analyzed experimentally determined DNA methylation maps of several species across the invertebrate-vertebrate boundary, to elucidate how vertebrate gene methylation has evolved. We show that, in contrast to the prevailing idea, a substantial number of promoters in an invertebrate basal chordate Ciona intestinalis are methylated. Moreover, gene expression data indicate significant, epigenomic context-dependent associations between promoter methylation and expression in C. intestinalis. However, there is no evidence that promoter methylation in invertebrate chordate has been evolutionarily maintained across the invertebrate-vertebrate boundary. Rather, body-methylated invertebrate genes preferentially obtain hypo-methylated promoters among vertebrates. Conversely, promoter methylation is preferentially found in lineage- and tissue-specific vertebrate genes. These results provide important insights into the evolutionary origin of epigenetic regulation of vertebrate gene expression.