DNA methylation rebalances gene dosage after mammalian gene duplications

6mA DNA Methylation on Genes in Plants Is Associated with Gene Complexity, Expression and Duplication

N6-methyladenine (6mA) DNA methylation has emerged as an important epigenetic modification in eukaryotes. Nevertheless, the evolution of the 6mA methylation of homologous genes after species and after gene duplications remains unclear in plants. To understand the evolution of 6mA methylation, we detected the genome-wide 6mA methylation patterns of four lotus plants (Nelumbo nucifera) from different geographic origins by nanopore sequencing and compared them to patterns in Arabidopsis and rice. Within lotus, the genomic distributions of 6mA sites are different from the widely studied 5mC methylation sites. Consistently, in lotus, Arabidopsis and rice, 6mA sites are enriched around transcriptional start sites, positively correlated with gene expression levels, and preferentially retained in highly and broadly expressed orthologs with longer gene lengths and more exons. Among different duplicate genes, 6mA methylation is significantly more enriched and conserved in whole-genome duplicates than in local duplicates. Overall, our study reveals the convergent patterns of 6mA methylation evolution based on both lineage and duplicate gene divergence, which underpin their potential role in gene regulatory evolution in plants.

DNA methylation signatures of duplicate gene evolution in angiosperms

Gene duplication is a source of evolutionary novelty. DNA methylation may play a role in the evolution of duplicate genes (paralogs) through its association with gene expression. While this relationship has been examined to varying extents in a few individual species, the generalizability of these results at either a broad phylogenetic scale with species of differing duplication histories or across a population remains unknown. We applied a comparative epigenomics approach to 43 angiosperm species across the phylogeny and a population of 928 Arabidopsis (Arabidopsis thaliana) accessions, examining the association of DNA methylation with paralog evolution. Genic DNA methylation was differentially associated with duplication type, the age of duplication, sequence evolution, and gene expression. Whole genome duplicates were typically enriched for CG-only gene-body methylated or unmethylated genes, while single-gene duplications were typically enriched for non-CG methylated or unmethylated genes. Non-CG methylation, in particular, was characteristic of more recent single-gene duplicates. Core angiosperm gene families differentiated into those which preferentially retain paralogs and 'duplication-resistant' families, which convergently reverted to singletons following duplication. Duplication-resistant families that still have paralogous copies were, uncharacteristically for core angiosperm genes, enriched for non-CG methylation. Non-CG methylated paralogs had higher rates of sequence evolution, higher frequency of presence-absence variation, and more limited expression. This suggests that silencing by non-CG methylation may be important to maintaining dosage following duplication and be a precursor to fractionation. Our results indicate that genic methylation marks differing evolutionary trajectories and fates between paralogous genes and have a role in maintaining dosage following duplication.

Evolutionary fates of gene-body methylation and its divergent association with gene expression in pigeonpea

Pigeonpea (Cajanus cajan L. Huth) is an agronomically important legume cultivated worldwide. In this study, we extensively analyzed gene-body methylation (GbM) patterns in pigeonpea. We found a bimodal distribution of CG and CHG methylation patterns. GbM features- slow evolution rate and increased length remained conserved. Genes with moderate CG body methylation showed highest expression where as highly-methylated genes showed lowest expression. Transposable element (TE)-related genes were methylated in multiple contexts and hence classified as C-methylated genes. A low expression among C-methylated genes was associated with transposons insertion in gene-body and upstream regulatory regions. The CG methylation patterns were found to be conserved in orthologs compared with non-CG methylation. By comparing methylation patterns between differentially methylated regions (DMRs) of the three genotypes, we found that variably methylated marks are less likely to target evolutionary conserved sequences. Finally, our analysis showed enrichment of nitrogen-related genes in GbM orthologs of legumes, which could be promising candidates for generating epialleles for crop improvement.

New insights into the evolution of CAF1 family and utilization of TaCAF1Ia1 specificity to reveal the origin of the maternal progenitor for common wheat

INTRODUCTION

Until now, the most likely direct maternal progenitor (AABB) for common wheat (AABBDD) has yet to be identified. Here, we try to solve this particular problem with the specificity of a novel gene family in wheat and by using large population of rare germplasm resources.

OBJECTIVES

Dissect the novelty of TaCAF1Ia subfamily in wheat. Exploit the conservative and specific characteristics of TaCAF1Ia1 to reveal the origin of the maternal progenitor for common wheat.

METHODS

Phylogenetic and collinear analysis of TaCAF1 genes were performed to identify the evolutionary specificity of TaCAF1Ia subfamily. The large-scale expression patterns and interaction patterns analysis of CCR4-NOT complex were used to clarify the expressed and structural specificity of TaCAF1Ia subfamily in wheat. The population resequencing and phylogeny analysis of the TaCAF1Ia1 were utilized for the traceability analysis to understand gene-pool exchanges during the transferring and subsequent development from tetraploid to hexaploidy wheat.

RESULTS

TaCAF1Ia is a novel non-typical CAF1 subfamily without DEDD (Asp-Glu-Asp-Asp) domain, whose members were extensively duplicated in wheat genome. The replication events had started and constantly evolved from ancestor species. Specifically, it was found that a key member CAF1Ia1 was highly specialized and only existed in the subB genome and S genome. Unlike CAF1s reported in other plants, TaCAF1Ia genes may be new factors for anther development. These atypical TaCAF1s could also form CCR4-NOT complex in wheat but with new interaction sites. Utilizing the particular but conserved characteristics of the TaCAF1Ia1 gene, the comparative analysis of haplotypes composition for TaCAF1Ia1 were identified among wheat populations with different ploidy levels. Based on this, the dual-lineages origin model of maternal progenitor for common wheat and potential three-lineages domestication model for cultivated tetraploid wheat were proposed.

CONCLUSION

This study brings fresh insights for revealing the origin of wheat and the function of CAF1 in plants.

Evolutionary implications of the RNA N6-methyladenosine methylome in plants

Epigenetic modifications play important roles in genome evolution and innovation. However, most analyses have focused on the evolutionary role of DNA modifications, and little is understood about the influence of post-transcriptional RNA modifications on genome evolution. To explore the evolutionary significance of RNA modifications, we generated transcriptome-wide profiles of N6-methyladenosine (m6A), the most prevalent internal modification of mRNA, for 13 representative plant species spanning over half a billion years of evolution. These data reveal the evolutionary conservation and divergence of m6A methylomes in plants, uncover the preference of m6A modifications on ancient orthologous genes, and demonstrate less m6A divergence between orthologous gene pairs with earlier evolutionary origins. Further investigation revealed that the evolutionary divergence of m6A modifications is related to sequence variation between homologs from whole genome duplication and gene family expansion from local genome duplication. Unexpectedly, a significant negative correlation was found between the retention ratio of m6A modifications and the number of family members. Moreover, the divergence of m6A modifications is accompanied by variation in the expression level and translation efficiency of duplicated genes from whole and local genome duplication. Our work reveals new insights into evolutionary patterns of m6A methylomes in plant species and their implications, and provides a resource of plant m6A profiles for further studies of m6A regulation and function in an evolutionary context.

Ampliconic Genes on the Great Ape Y Chromosomes: Rapid Evolution of Copy Number but Conservation of Expression Levels

Multicopy ampliconic gene families on the Y chromosome play an important role in spermatogenesis. Thus, studying their genetic variation in endangered great ape species is critical. We estimated the sizes (copy number) of nine Y ampliconic gene families in population samples of chimpanzee, bonobo, and orangutan with droplet digital polymerase chain reaction, combined these estimates with published data for human and gorilla, and produced genome-wide testis gene expression data for great apes. Analyzing this comprehensive data set within an evolutionary framework, we, first, found high inter- and intraspecific variation in gene family size, with larger families exhibiting higher variation as compared with smaller families, a pattern consistent with random genetic drift. Second, for four gene families, we observed significant interspecific size differences, sometimes even between sister species—chimpanzee and bonobo. Third, despite substantial variation in copy number, Y ampliconic gene families’ expression levels did not differ significantly among species, suggesting dosage regulation. Fourth, for three gene families, size was positively correlated with gene expression levels across species, suggesting that, given sufficient evolutionary time, copy number influences gene expression. Our results indicate high variability in size but conservation in gene expression levels in Y ampliconic gene families, significantly advancing our understanding of Y-chromosome evolution in great apes.

BLAST from the Past: Impacts of Evolving Approaches on Studies of Evolution by Gene Duplication

In 1970, Susumu Ohno hypothesized that gene duplication was a major reservoir of adaptive innovation. However, it was not until over two decades later that DNA sequencing studies uncovered the ubiquity of gene duplication across all domains of life, highlighting its global importance in the evolution of phenotypic complexity and species diversification. Today, it seems that there are no limits to the study of evolution by gene duplication, as it has rapidly coevolved with numerous experimental and computational advances in genomics. In this perspective, we examine word stem usage in PubMed abstracts to infer how evolving discoveries and technologies have shaped the landscape of studying evolution by gene duplication, leading to a more refined understanding of its role in the emergence of novel phenotypes.

The impact of epigenetic information on genome evolution

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

Selection following Gene Duplication Shapes Recent Genome Evolution in the Pea Aphid Acyrthosiphon pisum

Ecology of insects is as wide as their diversity, which reflects their high capacity of adaptation in most of the environments of our planet. Aphids, with over 4,000 species, have developed a series of adaptations including a high phenotypic plasticity and the ability to feed on the phloem sap of plants, which is enriched in sugars derived from photosynthesis. Recent analyses of aphid genomes have indicated a high level of shared ancestral gene duplications that might represent a basis for genetic innovation and broad adaptations. In addition, there are a large number of recent, species-specific gene duplications whose role in adaptation remains poorly understood. Here, we tested whether duplicates specific to the pea aphid Acyrthosiphon pisum are related to genomic innovation by combining comparative genomics, transcriptomics, and chromatin accessibility analyses. Consistent with large levels of neofunctionalization, we found that most of the recent pairs of gene duplicates evolved asymmetrically, showing divergent patterns of positive selection and gene expression. Genes under selection involved a plethora of biological functions, suggesting that neofunctionalization and tissue specificity, among other evolutionary mechanisms, have orchestrated the evolution of recent paralogs in the pea aphid and may have facilitated host-symbiont cooperation. Our comprehensive phylogenomics analysis allowed us to tackle the history of duplicated genes to pave the road toward understanding the role of gene duplication in ecological adaptation.

Gene Duplication in the Honeybee: Patterns of DNA Methylation, Gene Expression, and Genomic Environment

Gene duplication serves a critical role in evolutionary adaptation by providing genetic raw material to the genome. The evolution of duplicated genes may be influenced by epigenetic processes such as DNA methylation, which affects gene function in some taxa. However, the manner in which DNA methylation affects duplicated genes is not well understood. We studied duplicated genes in the honeybee Apis mellifera, an insect with a highly sophisticated social structure, to investigate whether DNA methylation was associated with gene duplication and genic evolution. We found that levels of gene body methylation were significantly lower in duplicate genes than in single-copy genes, implicating a possible role of DNA methylation in postduplication gene maintenance. Additionally, we discovered associations of gene body methylation with the location, length, and time since divergence of paralogous genes. We also found that divergence in DNA methylation was associated with divergence in gene expression in paralogs, although the relationship was not completely consistent with a direct link between DNA methylation and gene expression. Overall, our results provide further insight into genic methylation and how its association with duplicate genes might facilitate evolutionary processes and adaptation.

Reduced Translational Efficiency of Eukaryotic Genes after Duplication Events

Control of gene expression has been found to be predominantly determined at the level of protein translation. However, to date, reduced expression from duplicated genes in eukaryotes for dosage maintenance has only been linked to transcriptional control involving epigenetic mechanisms. Here, we hypothesize that dosage maintenance following gene duplication also involves regulation at the protein level. To test this hypothesis, we compared transcriptome and proteome data of yeast models, Saccharomyces cerevisiae and Schizosaccharomyces pombe, and worm models, Caenorhabditis elegans and Caenorhabditis briggsae, to investigate lineage-specifically duplicated genes. Duplicated genes in both eukaryotic models exhibited a reduced protein-to-mRNA abundance ratio. Moreover, dosage sensitive genes, represented by genes encoding protein complex subunits, reduced their protein-to-mRNA abundance ratios more significantly than the other genes after duplication events. An analysis of ribosome profiling (Ribo-Seq) data further showed that reduced translational efficiency was more prominent for dosage sensitive genes than for the other genes. Meanwhile, no difference in protein degradation rate was associated with duplication events. Translationally repressed duplicated genes were also more likely to be inhibited at the level of transcription. Taken together, these results suggest that translation-mediated dosage control is partially contributed by natural selection and it enhances transcriptional control in maintaining gene dosage after gene duplication events during eukaryotic genome evolution.

Copy number variations and young duplicate genes have high methylation levels in sticklebacks

Gene duplication is an important driver of genomic diversity that can promote adaptive evolution. However, like most mutations, a newly duplicated gene is often deleterious and removed from the genome by drift or natural selection. The early molecular changes that occur soon after duplication therefore may influence the long-term survival of gene duplicates, but relatively little empirical data exist on the events near the onset of duplication before mutations have time to accumulate. In this study, we contrast gene expression and DNA methylation levels of duplicate genes in the threespine stickleback, Gasterosteus aculeatus, including recently emerged duplications that segregate as copy number variations (CNVs). We find that younger duplicate genes have higher levels of promoter methylation than older genes, and that gene CNVs have higher promoter methylation than non-CNVs. These results suggest preferential duplication of highly methylated genes or rapid methylation changes soon after duplication. We also find a negative association between methylation and expression, providing a putative role for methylation in suppressing transcription that compensates for increases in gene copy numbers and promoting paralog retention. We propose that methylation contributes to the longevity of young duplicate genes, extending the window of opportunity for functional divergence via mutation.

Transcriptomic analysis reveals a controlling mechanism for NLRP3 and IL-17A in dextran sulfate sodium (DSS)-induced colitis

The incidence of inflammatory bowel disease (IBD) has markedly increased. Our research findings during the past showed that medicinal plant extracts and the derived phytochemical components from Wedelia chinensis (WC) can have strong anti-colitis activities. Here, we further identified the key component phytochemicals from active fractions of different WC preparations (WCHA) that are responsible for the protective effect of WCHA in colitis mice. Of the 3 major compounds (wedelolactone, luteolin and apigenin) in this fraction, luteolin had the highest anti-inflammatory effect in vivo. Using a next-generation sequencing (NGS) (e.g., RNA-seq) system to analyze the transcriptome of colorectal cells/tissues in mice with dextran sulfate sodium (DSS)-induced colitis with/without phytochemicals treatment, luteolin was found to strongly suppress the DSS-activated IL-17 pathway in colon tissue. In addition, co-treatment with wedelolactone and luteolin had a synergistic effect on the expression level of some IL-17 pathway-related genes. Interestingly, our NGS analyses also indicated that luteolin and wedelolactone can specifically suppress the expression of NLRP3 and NLRP1. Using a 3-dimensional cell co-culture system, we further demonstrated that luteolin could efficiently suppress NLRP3 expression via disruption of IL-17A signaling in inflamed colon tissue, which also indicates the pharmacological potential of luteolin and wedelolactone in treating IBD.

Genome-wide epigenetic variation among ash trees differing in susceptibility to a fungal disease

Background

European ash trees (Fraxinus excelsior) are currently threatened by ash dieback (ADB) caused by the fungus Hymenoscyphus fraxineus but a small percentage of the population possesses natural low susceptibility. The genome of a European ash tree has recently been sequenced. Here, we present whole genome DNA methylation data for two F. excelsior genotypes with high susceptibility to ADB, and two genotypes with low susceptibility, each clonally replicated. We also include two genotypes of Manchurian ash (F. mandshurica), an ash species which has co-evolved with H. fraxineus and also has low susceptibility to ADB.

Results

In F. excelsior, we find an average methylation level of 76.2% in the CG context, 52.0% in the CHG context, and 13.9% in the CHH context; similar levels to those of tomato. We find higher methylation in transposable elements as opposed to non-mobile elements, and high densities of Non-Differentially Methylation Positions (N-DMPs) in genes with housekeeping functions. Of genes putatively duplicated in whole genome duplication (WGD) events, an average of 25.9% are differentially methylated in at least one cytosine context, potentially indicative of unequal silencing. Variability in methylation patterns exists among clonal replicates, and this is only slightly less than the variability found between different genotypes. Of twenty genes previously found to have expression levels associated with ADB susceptibility, we find only two of these have differential methylation between high and low susceptibility F. excelsior trees. In addition, we identify 1683 significant Differentially Methylated Regions (DMRs) (q-value< 0.001) between the high and low susceptibility genotypes of F. excelsior trees, of which 665 remain significant when F. mandshurica samples are added to the low susceptibility group.

Conclusions

We find a higher frequency of differentially methylated WGD-derived gene duplicates in ash than other plant species previously studied. We also identify a set of genes with differential methylation between genotypes and species with high versus low susceptibility to ADB. This provides valuable foundational data for future work on the role that epigenetics may play in gene dosage compensation and susceptibility to ADB in ash.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4874-8) contains supplementary material, which is available to authorized users.

DNA methylation changes and evolution of RNA-based duplication in Sus scrofa: based on a two-step strategy

AIM

This study aims to couple DNA methylation changes and evolution of retrogenes.

MATERIALS & METHODS

A new two-step strategy was developed to screen retrogenes. Further, reduced representation bisulfite sequencing and RNA-seq data of eight tissues were used to analyze retrogenes.

RESULTS

A total of 964 retrocopies were identified and new retrocopies were available for the synthesis of glycans and lipids corresponding to pig phenotypic traits. Retrogenes were consistently hypermethylated. Hypomethylation of parental genes presented more susceptibility to retroposition. Promoter DNA methylation of retrogenes was negatively correlated with evolutionary time and played important roles in regulating retrogene tissue-specific expression pattern.

CONCLUSION

A two-step procedure is effective and necessary for identifying retrogenes. DNA methylation drives origination, survival, evolution and expression of retrogenes.

Recruitment of histone modifications to assist mRNA dosage maintenance after degeneration of cytosine DNA methylation during animal evolution

Following gene duplication, mRNA expression of the duplicated gene is reduced to maintain mRNA dosage. In mammals, this process is achieved with increased cytosine DNA methylation of the promoters of duplicated genes to suppress transcriptional initiation. However, not all animal species possess a full apparatus for cytosine DNA methylation. For such species, such as the roundworm (Caenorhabditis elegans, "worm" hereafter) or fruit fly (Drosophila melanogaster, "fly" hereafter), it is unclear how reduced expression of duplicated genes has been achieved evolutionarily. Here, we hypothesize that in the absence of a classical cytosine DNA methylation pathway, histone modifications play an increasing role in maintaining mRNA dosage following gene duplication. We initially verified that reduced gene expression of duplicated genes had occurred in the worm, fly, and mouse (Mus musculus). Next, several histone marks, with the capacity to control mRNA abundance in the models studied, were examined. In the worm and fly, but not in the mouse, multiple histone modifications were found to assist mRNA dosage maintenance following gene duplication events and the possible involvement of adenine DNA methylation in this process was excluded. Furthermore, the histone marks and acting regions that mediated the reduction in duplicated gene expression were found to be largely organism specific. Thus, it appears that many of the histone marks that maintain mRNA dosage were independently recruited during the evolution of worms and flies to compensate for the loss of cytosine DNA methylation machinery from their genomes.

Underrepresentation of active histone modification marks in evolutionarily young genes

It is known that evolutionarily new genes can rapidly evolve essential roles in fundamental biological processes. Nevertheless, the underlying molecular mechanism of how they acquire their novel transcriptional pattern is less characterized except for the role of cis-regulatory evolution. Epigenetic modification offers an alternative possibility. Here, we examined how histone modifications have changed among different gene age groups in Drosophila melanogaster by integrative analyses of an updated new gene dataset and published epigenomic data. We found a robust pattern across various datasets where both the coverage and intensity of active histone modifications, histone 3 lysine 4 trimethylation and lysine 36 trimethylation, increased with evolutionary age. Such a temporal correlation is negative and much weaker for the repressive histone mark, lysine 9 trimethylation, which is expected given its major association with heterochromatin. By further comparison with neighboring old genes, the depletion of active marks of new genes could be only partially explained by the local epigenetic context. All these data are consistent with the observation that older genes bear relatively higher expression levels and suggest that the evolution of histone modifications could be implicated in transcriptional evolution after gene birth.

Facilitation of environmental adaptation and evolution by epigenetic phenotype variation: insights from clonal, invasive, polyploid, and domesticated animals

There is increasing evidence, particularly from plants, that epigenetic mechanisms can contribute to environmental adaptation and evolution. The present article provides an overview on this topic for animals and highlights the special suitability of clonal, invasive, hybrid, polyploid, and domesticated species for environmental and evolutionary epigenetics. Laboratory and field studies with asexually reproducing animals have shown that epigenetically diverse phenotypes can be produced from the same genome either by developmental stochasticity or environmental induction. The analysis of invasions revealed that epigenetic phenotype variation may help to overcome genetic barriers typically associated with invasions such as bottlenecks and inbreeding. Research with hybrids and polyploids established that epigenetic mechanisms are involved in consolidation of speciation by contributing to reproductive isolation and restructuring of the genome in the neo-species. Epigenetic mechanisms may even have the potential to trigger speciation but evidence is still meager. The comparison of domesticated animals and their wild ancestors demonstrated heritability and selectability of phenotype modulating DNA methylation patterns. Hypotheses, model predictions, and empirical results are presented to explain how epigenetic phenotype variation could facilitate adaptation and speciation. Clonal laboratory lineages, monoclonal invaders, and adaptive radiations of different evolutionary age seem particularly suitable to empirically test the proposed ideas. A respective research agenda is presented.

Global DNA methylation variations after short-term heat shock treatment in cultured microspores of Brassica napus cv. Topas

Heat stress can induce the cultured microspores into embryogenesis. In this study, whole genome bisulphite sequencing was employed to study global DNA methylation variations after short-term heat shock (STHS) treatments in cultured microspores of Brassica napus cv. Topas. Our results indicated that treatment on cultured Topas microspores at 32 °C for 6 h triggered DNA hypomethylation, particularly in the CG and CHG contexts. And the total number of T32 (Topas 32 °C for 6 h) vs. T0 (Topas 0 h) differentially methylated region-related genes (DRGs) was approximately two-fold higher than that of T18 (Topas 18 °C for 6 h) vs. T0 DRGs, which suggested that 32 °C might be a more intense external stimulus than 18 °C resulting in more changes in the DNA methylation status of cultured microspores. Additionally, 32 °C treatment for 6 h led to increased CHG differential methylations of transposons (DMTs), which were mainly constituted by overlaps between the hypomethylated differentially methylated regions (hypo-DMRs) and transposon elements (TEs). Further analysis demonstrated that the DRGs and their paralogs exhibited differential methylated/demethylated patterns. To summarize, the present study is the first methylome analysis of cultured microspores in response to STHS and may provide valuable information on the roles of DNA methylation in heat response.

Divergent DNA Methylation Provides Insights into the Evolution of Duplicate Genes in Zebrafish

The evolutionary mechanism, fate and function of duplicate genes in various taxa have been widely studied; however, the mechanism underlying the maintenance and divergence of duplicate genes in Danio rerio remains largely unexplored. Whether and how the divergence of DNA methylation between duplicate pairs is associated with gene expression and evolutionary time are poorly understood. In this study, by analyzing bisulfite sequencing (BS-seq) and RNA-seq datasets from public data, we demonstrated that DNA methylation played a critical role in duplicate gene evolution in zebrafish. Initially, we found promoter methylation of duplicate genes generally decreased with evolutionary time as measured by synonymous substitution rate between paralogous duplicates (Ks). Importantly, promoter methylation of duplicate genes was negatively correlated with gene expression. Interestingly, for 665 duplicate gene pairs, one gene was consistently promoter methylated, while the other was unmethylated across nine different datasets we studied. Moreover, one motif enriched in promoter methylated duplicate genes tended to be bound by the transcription repression factor FOXD3, whereas a motif enriched in the promoter unmethylated sequences interacted with the transcription activator Sp1, indicating a complex interaction between the genomic environment and epigenome. Besides, body-methylated genes showed longer length than body-unmethylated genes. Overall, our results suggest that DNA methylation is highly important in the differential expression and evolution of duplicate genes in zebrafish.

Evolution of Epigenetic Regulation in Vertebrate Genomes

Empirical models of sequence evolution have spurred progress in the field of evolutionary genetics for decades. We are now realizing the importance and complexity of the eukaryotic epigenome. While epigenome analysis has been applied to genomes from single-cell eukaryotes to human, comparative analyses are still relatively few and computational algorithms to quantify epigenome evolution remain scarce. Accordingly, a quantitative model of epigenome evolution remains to be established. We review here the comparative epigenomics literature and synthesize its overarching themes. We also suggest one mechanism, transcription factor binding site (TFBS) turnover, which relates sequence evolution to epigenetic conservation or divergence. Lastly, we propose a framework for how the field can move forward to build a coherent quantitative model of epigenome evolution.

Global analysis of human duplicated genes reveals the relative importance of whole-genome duplicates originated in the early vertebrate evolution

Background

Gene duplication is a genetic mutation that creates functionally redundant gene copies that are initially relieved from selective pressures and may adapt themselves to new functions with time. The levels of gene duplication may vary from small-scale duplication (SSD) to whole genome duplication (WGD). Studies with yeast revealed ample differences between these duplicates: Yeast WGD pairs were functionally more similar, less divergent in subcellular localization and contained a lesser proportion of essential genes. In this study, we explored the differences in evolutionary genomic properties of human SSD and WGD genes, with the identifiable human duplicates coming from the two rounds of whole genome duplication occurred early in vertebrate evolution.

Results

We observed that these two groups of duplicates were also dissimilar in terms of their evolutionary and genomic properties. But interestingly, this is not like the same observed in yeast. The human WGDs were found to be functionally less similar, diverge more in subcellular level and contain a higher proportion of essential genes than the SSDs, all of which are opposite from yeast. Additionally, we explored that human WGDs were more divergent in their gene expression profile, have higher multifunctionality and are more often associated with disease, and are evolutionarily more conserved than human SSDs.

Conclusions

Our study suggests that human WGD duplicates are more divergent and entails the adaptation of WGDs to novel and important functions that consequently lead to their evolutionary conservation in the course of evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2392-0) contains supplementary material, which is available to authorized users.

Rapid divergence and diversification of mammalian duplicate gene functions

BACKGROUND

Gene duplication provides raw material for the evolution of functional innovation. We recently developed a phylogenetic method that classifies evolutionary processes driving the retention of duplicate genes by quantifying divergence between their spatial gene expression profiles and that of their single-copy orthologous gene in a closely related sister species.

RESULTS

Here, we apply our classification method to pairs of duplicate genes in eight mammalian genomes, using data from 11 tissues to construct spatial gene expression profiles. We find that young mammalian duplicates are often functionally conserved, and that expression divergence rapidly increases over evolutionary time. Moreover, expression divergence results in increased tissue specificity, with an overrepresentation of expression in male kidney, underrepresentation of expression in female liver, and strong underrepresentation of expression in testis. Thus, duplicate genes acquire a diversity of new tissue-specific functions outside of the testis, possibly contributing to the origin of a multitude of complex phenotypes during mammalian evolution.

CONCLUSIONS

Our findings reveal that mammalian duplicate genes are initially functionally conserved, and then undergo rapid functional divergence over evolutionary time, acquiring diverse tissue-specific biological roles. These observations are in stark contrast to the much faster expression divergence and acquisition of broad housekeeping roles we previously observed in Drosophila duplicate genes. Due to the smaller effective population sizes of mammals relative to Drosophila, these analyses implicate natural selection in the functional evolution of duplicate genes.

Closing the genotype-phenotype gap: emerging technologies for evolutionary genetics in ecological model vertebrate systems

The analysis of genetic and epigenetic mechanisms of the genotype-phenotypic connection has, so far, only been possible in a handful of genetic model systems. Recent technological advances, including next-generation sequencing methods such as RNA-seq, ChIP-seq and RAD-seq, and genome-editing approaches including CRISPR-Cas, now permit to address these fundamental questions of biology also in organisms that have been studied in their natural habitats. We provide an overview of the benefits and drawbacks of these novel techniques and experimental approaches that can now be applied to ecological and evolutionary vertebrate models such as sticklebacks and cichlid fish. We can anticipate that these new methods will increase the understanding of the genetic and epigenetic factors influencing adaptations and phenotypic variation in ecological settings. These new arrows in the methodological quiver of ecologist will drastically increase the understanding of the genetic basis of adaptive traits - leading to a further closing of the genotype-phenotype gap.

Divergence of gene body DNA methylation and evolution of plant duplicate genes

It has been shown that gene body DNA methylation is associated with gene expression. However, whether and how deviation of gene body DNA methylation between duplicate genes can influence their divergence remains largely unexplored. Here, we aim to elucidate the potential role of gene body DNA methylation in the fate of duplicate genes. We identified paralogous gene pairs from Arabidopsis and rice (Oryza sativa ssp. japonica) genomes and reprocessed their single-base resolution methylome data. We show that methylation in paralogous genes nonlinearly correlates with several gene properties including exon number/gene length, expression level and mutation rate. Further, we demonstrated that divergence of methylation level and pattern in paralogs indeed positively correlate with their sequence and expression divergences. This result held even after controlling for other confounding factors known to influence the divergence of paralogs. We observed that methylation level divergence might be more relevant to the expression divergence of paralogs than methylation pattern divergence. Finally, we explored the mechanisms that might give rise to the divergence of gene body methylation in paralogs. We found that exonic methylation divergence more closely correlates with expression divergence than intronic methylation divergence. We show that genomic environments (e.g., flanked by transposable elements and repetitive sequences) of paralogs generated by various duplication mechanisms are associated with the methylation divergence of paralogs. Overall, our results suggest that the changes in gene body DNA methylation could provide another avenue for duplicate genes to develop differential expression patterns and undergo different evolutionary fates in plant genomes.

Distinct Copy Number, Coding Sequence, and Locus Methylation Patterns Underlie Rhg1-Mediated Soybean Resistance to Soybean Cyst Nematode

Copy number variation of kilobase-scale genomic DNA segments, beyond presence/absence polymorphisms, can be an important driver of adaptive traits. Resistance to Heterodera glycines (Rhg1) is a widely utilized quantitative trait locus that makes the strongest known contribution to resistance against soybean cyst nematode (SCN), Heterodera glycines, the most damaging pathogen of soybean (Glycine max). Rhg1 was recently discovered to be a complex locus at which resistance-conferring haplotypes carry up to 10 tandem repeat copies of a 31-kb DNA segment, and three disparate genes present on each repeat contribute to SCN resistance. Here, we use whole-genome sequencing, fiber-FISH (fluorescence in situ hybridization), and other methods to discover the genetic variation at Rhg1 across 41 diverse soybean accessions. Based on copy number variation, transcript abundance, nucleic acid polymorphisms, and differentially methylated DNA regions, we find that SCN resistance is associated with multicopy Rhg1 haplotypes that form two distinct groups. The tested high-copy-number Rhg1 accessions, including plant introduction (PI) 88788, contain a flexible number of copies (seven to 10) of the 31-kb Rhg1 repeat. The identified low-copy-number Rhg1 group, including PI 548402 (Peking) and PI 437654, contains three copies of the Rhg1 repeat and a newly identified allele of Glyma18g02590 (a predicted α-SNAP [α-soluble N-ethylmaleimide-sensitive factor attachment protein]). There is strong evidence for a shared origin of the two resistance-conferring multicopy Rhg1 groups and subsequent independent evolution. Differentially methylated DNA regions also were identified within Rhg1 that correlate with SCN resistance. These data provide insights into copy number variation of multigene segments, using as the example a disease resistance trait of high economic importance.

DNA methylation and evolution of duplicate genes

The evolutionary mechanisms underlying duplicate gene maintenance and divergence remain highly debated. Epigenetic modifications, such as DNA methylation, may contribute to duplicate gene evolution by facilitating tissue-specific regulation. However, the role of epigenetic divergence on duplicate gene evolution remains little understood. Here we show, using comprehensive data across 10 diverse human tissues, that DNA methylation plays critical roles in several aspects of duplicate gene evolution. We first demonstrate that duplicate genes are initially heavily methylated, before gradually losing DNA methylation as they age. Within each pair, DNA methylation divergence between duplicate partners increases with evolutionary age. Importantly, tissue-specific DNA methylation of duplicates correlates with tissue-specific expression, implicating DNA methylation as a causative factor for functional divergence of duplicate genes. These patterns are apparent in promoters but not in gene bodies, in accord with the complex relationship between gene-body DNA methylation and transcription. Remarkably, many duplicate gene pairs exhibit consistent division of DNA methylation across multiple, divergent tissues: For the majority (73%) of duplicate gene pairs, one partner is always hypermethylated compared with the other. This is indicative of a common underlying determinant of DNA methylation. The division of DNA methylation is also consistent with their chromatin accessibility profiles. Moreover, at least two sequence motifs known to interact with the Sp1 transcription factor mark promoters of more hypomethylated duplicate partners. These results demonstrate critical roles of DNA methylation, as well as complex interaction between genome and epigenome, on duplicate gene evolution.

Accumulation of CTCF-binding sites drives expression divergence between tandemly duplicated genes in humans

Background

During eukaryotic genome evolution, tandem gene duplication is the most frequent event giving rise to clustered gene families. However, how expression divergence between tandemly duplicated genes has emerged and maintained remain unclear. In particular, it is unknown if epigenetic regulators have been involved in the process.

Results

We demonstrate that CCCTC-binding factor (CTCF), the master epigenetic regulator and the only known insulator protein in humans, has played a predominant role in generating divergence in both expression profiles and expression levels between adjacent paralogs in the human genome. This phenomenon was not observed for non-paralogous adjacent genes. After tandem duplication events, CTCF-binding sites gradually accumulate between paralogs. This trend was more prominent for genes involved in particular functions.

Conclusions

The accumulation of CTCF-binding sites drives expression divergence of tandemly duplicated genes. This process is likely targeted by natural selection. Our study reveals the importance of CTCF to the evolution of animal diversity and complexity.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-S1-S8) contains supplementary material, which is available to authorized users.

Flagellated algae protein evolution suggests the prevalence of lineage-specific rules governing evolutionary rates of eukaryotic proteins

Understanding the general rules governing the rate of protein evolution is fundamental to evolutionary biology. However, attempts to address this issue in yeasts and mammals have revealed considerable differences in the relative importance of determinants for protein evolutionary rates. This phenomenon was previously explained by the fact that yeasts and mammals are different in many cellular and genomic properties. Flagellated algae species have several cellular and genomic characteristics that are intermediate between yeasts and mammals. Using partial correlation analyses on the evolution of 6,921 orthologous proteins from Chlamydomonas reinhardtii and Volvox carteri, we examined factors influencing evolutionary rates of proteins in flagellated algae. Previous studies have shown that mRNA abundance and gene compactness are strong determinants for protein evolutionary rates in yeasts and mammals, respectively. We show that both factors also influence algae protein evolution with mRNA abundance having a larger impact than gene compactness on the rates of algae protein evolution. More importantly, among all the factors examined, coding sequence (CDS) length has the strongest (positive) correlation with protein evolutionary rates. This correlation between CDS length and the rates of protein evolution is not due to alignment-related issues or domain density. These results suggest no simple and universal rules governing protein evolutionary rates across different eukaryotic lineages. Instead, gene properties influence the rate of protein evolution in a lineage-specific manner.

Functional characterization of motif sequences under purifying selection

Diverse life forms are driven by the evolution of gene regulatory programs including changes in regulator proteins and cis-regulatory elements. Alterations of cis-regulatory elements are likely to dominate the evolution of the gene regulatory networks, as they are subjected to smaller selective constraints compared with proteins and hence may evolve quickly to adapt the environment. Prior studies on cis-regulatory element evolution focus primarily on sequence substitutions of known transcription factor-binding motifs. However, evolutionary models for the dynamics of motif occurrence are relatively rare, and comprehensive characterization of the evolution of all possible motif sequences has not been pursued. In the present study, we propose an algorithm to estimate the strength of purifying selection of a motif sequence based on an evolutionary model capturing the birth and death of motif occurrences on promoters. We term this measure as the 'evolutionary retention coefficient', as it is related yet distinct from the canonical definition of selection coefficient in population genetics. Using this algorithm, we estimate and report the evolutionary retention coefficients of all possible 10-nucleotide sequences from the aligned promoter sequences of 27 748. orthologous gene families in 34 mammalian species. Intriguingly, the evolutionary retention coefficients of motifs are intimately associated with their functional relevance. Top-ranking motifs (sorted by evolutionary retention coefficients) are significantly enriched with transcription factor-binding sequences according to the curated knowledge from the TRANSFAC database and the ChIP-seq data generated from the ENCODE Consortium. Moreover, genes harbouring high-scoring motifs on their promoters retain significantly coherent expression profiles, and those genes are over-represented in the functional classes involved in gene regulation. The validation results reveal the dependencies between natural selection and functions of cis-regulatory elements and shed light on the evolution of gene regulatory networks.

Escape from Adaptive Conflict follows from weak functional trade-offs and mutational robustness

A fundamental question in molecular evolution is how proteins can adapt to new functions while being conserved for an existing function at the same time. Several theoretical models have been put forward to explain this apparent paradox. The most popular models include neofunctionalization, subfunctionalization (SUBF) by degenerative mutations, and dosage models. All of these models focus on adaptation after gene duplication. A newly proposed model named "Escape from Adaptive Conflict" (EAC) includes adaptive processes before and after gene duplication that lead to multifunctional proteins, and divergence (SUBF). Support for the importance of multifunctionality for the evolution of new protein functions comes from two experimental observations. First, many enzymes have highly evolvable promiscuous side activities. Second, different structural states of the same protein can be associated with different functions. How these observations may be related to the EAC model, under which conditions EAC is possible, and how the different models relate to each other is still unclear. Here, we present a theoretical framework that uses biophysical principles to infer the roles of functional promiscuity, gene dosage, gene duplication, point mutations, and selection pressures in the evolution of proteins. We find that selection pressures can determine whether neofunctionalization or SUBF is the more likely evolutionary process. Multifunctional proteins, arising during EAC evolution, allow rapid adaptation independent of gene duplication. This becomes a crucial advantage when gene duplications are rare. Finally, we propose that an increase in mutational robustness, not necessarily functional optimization, can be the sole driving force behind SUBF.

Transcriptomes of mouse olfactory epithelium reveal sexual differences in odorant detection

To sense numerous odorants and chemicals, animals have evolved a large number of olfactory receptor genes (Olfrs) in their genome. In particular, the house mouse has ∼1,100 genes in the Olfr gene family. This makes the mouse a good model organism to study Olfr genes and olfaction-related genes. To date, whether male and female mice possess the same ability in detecting environmental odorants is still unknown. Using the next generation sequencing technology (paired-end mRNA-seq), we detected 1,088 expressed Olfr genes in both male and female olfactory epithelium. We found that not only Olfr genes but also odorant-binding protein (Obp) genes have evolved rapidly in the mouse lineage. Interestingly, Olfr genes tend to express at a higher level in males than in females, whereas the Obp genes clustered on the X chromosome show the opposite trend. These observations may imply a more efficient odorant-transporting system in females, whereas a more active Olfr gene expressing system in males. In addition, we detected the expression of two genes encoding major urinary proteins, which have been proposed to bind and transport pheromones or act as pheromones in mouse urine. This observation suggests a role of main olfactory system (MOS) in pheromone detection, contrary to the view that only accessory olfactory system (AOS) is involved in pheromone detection. This study suggests the sexual differences in detecting environmental odorants in MOS and demonstrates that mRNA-seq provides a powerful tool for detecting genes with low expression levels and with high sequence similarities.

DNA methylation in amphioxus: from ancestral functions to new roles in vertebrates

In vertebrates, DNA methylation is an epigenetic mechanism that modulates gene transcription, and plays crucial roles during development, cell fate maintenance, germ cell pluripotency and inheritable genome imprinting. DNA methylation might also play a role as a genome defense mechanism against the mutational activity derived from transposon mobility. In contrast to the heavily methylated genomes in vertebrates, most genomes in invertebrates are poorly or just moderately methylated, and the function of DNA methylation remains unclear. Here, we review the DNA methylation system in the cephalochordate amphioxus, which belongs to the most basally divergent group of our own phylum, the chordates. First, surveys of the amphioxus genome database reveal the presence of the DNA methylation machinery, DNA methyltransferases and methyl-CpG-binding domain proteins. Second, comparative genomics and analyses of conserved synteny between amphioxus and vertebrates provide robust evidence that the DNA methylation machinery of amphioxus represents the ancestral toolkit of chordates, and that its expansion in vertebrates was originated by the two rounds of whole-genome duplication that occurred in stem vertebrates. Third, in silico analysis of CpGo/e ratios throughout the amphioxus genome suggests a bimodal distribution of DNA methylation, consistent with a mosaic pattern comprising domains of methylated DNA interspersed with domains of unmethylated DNA, similar to the situation described in ascidians, but radically different to the globally methylated vertebrate genomes. Finally, we discuss potential roles of the DNA methylation system in amphioxus in the context of chordate genome evolution and the origin of vertebrates.