Search

Daudi Jjingo 1 , Andrew B. Conley 1 , Soojin V. Yi 1 , Victoria V. Lunyak 2 and I. King Jordan 1,3 1 School of Biology, Georgia Institute of Technology, Atlanta GA, 30332 USA 2 Buck Institute for Research on Aging, Novato CA, 94945 USA 3 PanAmerican Bioinformatics Institute, Santa Marta, Magdalena, Colombia Received: March 29, 2012, Accepted: April 28, 2012, Published: May 9, 2012 Keywords: genome-wide methylation, epigenetic mark, intragenic transcription, methylating enzyme complexes Correspondence: I. King Jordan, // // Abstract DNA methylation of promoter sequences is a repressive epigenetic mark that down-regulates gene expression. However, DNA methylation is more prevalent within gene-bodies than seen for promoters, and gene-body methylation has been observed to be positively correlated with gene expression levels. This paradox remains unexplained, and accordingly the role of DNA methylation in gene-bodies is poorly understood. We addressed the presence and role of human gene-body DNA methylation using a meta-analysis of human genome-wide methylation, expression and chromatin data sets. Methylation is associated with transcribed regions as genic sequences have higher levels of methylation than intergenic or promoter sequences. We also find that the relationship between gene-body DNA methylation and expression levels is non-monotonic and bell-shaped. Mid-level expressed genes have the highest levels of gene-body methylation, whereas the most lowly and highly expressed sets of genes both have low levels of methylation. While gene-body methylation can be seen to efficiently repress the initiation of intragenic transcription, the vast majority of methylated sites within genes are not associated with intragenic promoters. In fact, highly expressed genes initiate the most intragenic transcription, which is inconsistent with the previously held notion that gene-body methylation serves to repress spurious intragenic transcription to allow for efficient transcriptional elongation. These observations lead us to propose a model to explain the presence of human gene-body methylation. This model holds that the repression of intragenic transcription by gene-body methylation is largely epiphenomenal, and suggests that gene-body methylation levels are predominantly shaped via the accessibility of the DNA to methylating enzyme complexes.

Theories predict that the evolutionary rates of X-linked regions can differ from those of autosomal regions. The male-biased mutation theory predicts a slower rate of neutral substitution on the X chromosome (slow-X evolution), as the X spends less time in male germlines, where more mutations originate per generation than in female germlines. The fast-X theory, however, predicts a faster rate of adaptive substitution on the X chromosome when newly arising beneficial mutations are, on average, partially recessive (fast-X evolution), as the X enjoys a greater efficacy of positive selection. The slow- and fast-X processes are expected to interact as the degree of male-biased mutation can in turn influence the relative rate of adaptive evolution on the X. Here, we investigate lineage-specific variation in, and the interaction of, slow- and fast-X processes using genomic data from four primates. We find consistent evidence for slow-X evolution in all lineages. In contrast, evidence for fast-X evolution exists in only a subset of lineages. In particular, the marmoset lineage, which shows the strongest evidence of fast-X, exhibits the lowest male mutation bias. We discuss the possible interaction between slow- and fast-X evolution and other factors that influence the degrees of slow- and fast-X evolution.

In complex organisms, different tissues express different genes, which ultimately shape the function and phenotype of each tissue. An important goal of modern biology is to understand how some genes are turned on and off in specific tissues and how the numbers of different gene expression products are determined. These aspects are named 'expression breadth' (or 'tissue specificity') and 'expression level', respectively. Here, we show that we can predict substantial amount of variation in levels and breadths of gene expression using genomic information of each gene. Interestingly, many genomic traits are correlated with both aspects of gene expression in similar directions, suggesting shared molecular pathways. However, to elucidate distinctive molecular mechanisms governing gene expression levels and breadths, we need to identify the relative significance of each genomic trait on these two aspects of gene expression. To this end, we developed a novel multivariate multiple regression method. Using this new method, we show that gene compactness (in particular, the mean size of exons), codon usage bias and non-synonymous rates have a stronger influence on expression levels compared with their effects on expression breadths. In contrast, the propensity of promoter DNA methylation is a stronger indicator of expression breadths than of expression levels. Interestingly, intron DNA methylation exhibits an opposite pattern to the promoter DNA methylation in the human genome, suggesting that DNA methylation may play multiple roles depending upon its genomic targets. Furthermore, synonymous rates have stronger associations with expression breadths than with expression levels in the human genome. These findings provide clues toward distinctive molecular mechanisms regulating different aspects of gene expression.

imb_1092 553..566 DNA methylation plays an important role in gene regulation in animals. However, the evolution and function of DNA methylation has only recently emerged as the subject of widespread study in insects. In this review we profile the known distribution of DNA methylation systems across insect taxa and synthesize functional inferences from studies of DNA methylation in insects and vertebrates. Unlike vertebrate genomes, which tend to be globally methylated, DNA methylation is primarily targeted to genes in insects. Nevertheless, mounting evidence suggests that a specialized role exists for genic methylation in the regulation of transcription, and possibly mRNA splicing, in both insects and mammals. Investigations in several insect taxa further reveal that DNA methylation is preferentially targeted to ubiquitously expressed genes and may play a key role in the regulation of phenotypic plasticity. We suggest that insects are particularly amenable to advancing our understanding of the biological functions of DNA methylation, because insects are evolutionarily diverse, display several lineage-specific losses of DNA methylation and possess tractable patterns of DNA methylation in moderately sized genomes.

The functional and evolutionary significance of DNA methylation in insect genomes remains to be resolved. Nasonia is well situated for comparative analyses of DNA methylation and genome evolution, since the genomes of a moderately distant outgroup species as well as closely related sibling species are available. Using direct sequencing of bisulfite-converted DNA, we uncovered a substantial level of DNA methylation in 17 of 18 Nasonia vitripennis genes and a strong correlation between methylation level and CpG depletion. Notably, in the sex-determining locus transformer, the exon that is alternatively spliced between the sexes is heavily methylated in both males and females, whereas other exons are only sparsely methylated. Orthologous genes of the honeybee and Nasonia show highly similar relative levels of CpG depletion, despite ~190 My divergence. Densely and sparsely methylated genes in these species also exhibit similar functional enrichments. We found that the degree of CpG depletion is negatively correlated with substitution rates between closely related Nasonia species for synonymous, nonsynonymous, and intron sites. This suggests that mutation rates increase with decreasing levels of germ line methylation. Thus, DNA methylation is prevalent in the Nasonia genome, may participate in regulatory processes such as sex determination and alternative splicing, and is correlated with several aspects of genome and sequence evolution.

We investigated the biological significance of path lengths in 12 protein-protein interaction (PPI) networks. We put forward three predictions, based on the idea that biological complexity influences path lengths. First, at the network level, path lengths are generally longer in PPIs than in random networks. Second, this pattern is more pronounced in more complex organisms. Third, within a PPI network, path lengths of individual proteins are biologically significant. We found that in 11 of the 12 species, average path lengths in PPI networks are significantly longer than those in randomly rewired networks. The PPI network of the malaria parasite Plasmodium falciparum, however, does not exhibit deviation from rewired networks. Furthermore, eukaryotic PPIs exhibit significantly greater deviation from randomly rewired networks than prokaryotic PPIs. Thus our study highlights the potentially meaningful variation in path lengths of PPI networks. Moreover, node eccentricity, defined as the longest path from a protein to others, is significantly correlated with the levels of gene expression and dispensability in the yeast PPI network. We conclude that biological complexity influences both global and local properties of path lengths in PPI networks. Investigating variation of path lengths may provide new tools to analyze the evolution of functional modules in biological systems.

Ants have evolved very complex societies and are key ecosystem members. Some ants, such as the fire ant Solenopsis invicta , are also major pests. Here, we present a draft genome of S. invicta , assembled from Roche 454 and Illumina sequencing reads obtained from a focal haploid male and his brothers. We used comparative genomic methods to obtain insight into the unique features of the S. invicta genome. For example, we found that this genome harbors four adjacent copies of vitellogenin. A phylogenetic analysis revealed that an ancestral vitellogenin gene first underwent a duplication that was followed by possibly independent duplications of each of the daughter vitellogenins. The vitellogenin genes have undergone subfunctionalization with queen- and worker-specific expression, possibly reflecting differential selection acting on the queen and worker castes. Additionally, we identified more than 400 putative olfactory receptors of which at least 297 are intact. This represents the largest repertoire reported so far in insects. S. invicta also harbors an expansion of a specific family of lipid-processing genes, two putative orthologs to the transformer/feminizer sex differentiation gene, a functional DNA methylation system, and a single putative telomerase ortholog. EST data indicate that this S. invicta telomerase ortholog has at least four spliceforms that differ in their use of two sets of mutually exclusive exons. Some of these and other unique aspects of the fire ant genome are likely linked to the complex social behavior of this species.

Phenotypic plasticity allows organisms to produce alternative phenotypes under different conditions and represents one of the most important ways by which organisms adaptively respond to the environment. However, the relationship between phenotypic plasticity and molecular evolution remains poorly understood. We addressed this issue by investigating the evolution of genes associated with phenotypically plastic castes, sexes, and developmental stages of the fire ant Solenopsis invicta . We first determined if genes associated with phenotypic plasticity in S. invicta evolved at a rapid rate, as predicted under theoretical models. We found that genes differentially expressed between S. invicta castes, sexes, and developmental stages all exhibited elevated rates of evolution compared with ubiquitously expressed genes. We next investigated the evolutionary history of genes associated with the production of castes. Surprisingly, we found that orthologs of caste-biased genes in S. invicta and the social bee Apis mellifera evolved rapidly in lineages without castes. Thus, in contrast to some theoretical predictions, our results suggest that rapid rates of molecular evolution may not arise primarily as a consequence of phenotypic plasticity. Instead, genes evolving under relaxed purifying selection may more readily adopt new forms of biased expression during the evolution of alternate phenotypes. These results suggest that relaxed selective constraint on protein-coding genes is an important and underappreciated element in the evolutionary origin of phenotypic plasticity.

MicroRNAs (miRNAs) are ancient, short noncoding RNA molecules that regulate the transcriptome through post-transcriptional mechanisms. miRNA riboregulation is involved in a diverse range of biological processes, and misregulation is implicated in disease. It is generally thought that miRNAs function to canalize cellular outputs, for instance as "fail-safe" repressors of gene misexpression. Genomic surveys in humans have revealed reduced genetic polymorphism and the signature of negative selection for both miRNAs themselves and the target sequences to which they are predicted to bind. We investigated the evolution of miRNAs and their binding sites across cichlid fishes from Lake Malawi (East Africa), where hundreds of diverse species have evolved in the last million years. Using low-coverage genome sequence data, we identified 100 cichlid miRNA genes with mature regions that are highly conserved in other animal species. We computationally predicted target sites on the 3'-untranslated regions (3'-UTRs) of cichlid genes to which miRNAs may bind and found that these sites possessed elevated single nucleotide polymorphism (SNP) densities. Furthermore, polymorphic sites in predicted miRNA targets showed higher minor allele frequencies on average and greater genetic differentiation between Malawi lineages when compared with a neutral expectation and nontarget 3'-UTR SNPs. Our data suggest that divergent selection on miRNA riboregulation may have contributed to the diversification of cichlid species and may similarly play a role in rapid phenotypic evolution of other natural systems.

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., Author summary Adrenergic α2C receptor (ADRA2C) is a regulator of the fight-or-flight response. It has been shown in mice that repression of this gene can result in relevant physiological and behavioral alterations. A strong selection signature in the genomes of domesticated chickens has been reported for this gene, suggesting that less aggression toward humans has been positively selected during chicken domestication. In this work, we analyze the genomes, transcriptomes, and epigenomes of a large number of humans and chimpanzees along with those of other primates to propose that repression of this gene has been positively selected in the evolution of humans and chimpanzees.

The co-existence of several ploidy types in natural populations makes the cyprinid loach Misgurnus anguillicaudatus an exciting model system to study the genetic and phenotypic consequences of ploidy variations. A first step in such effort is to identify the specific ploidy of an individual. Currently popular methods of karyotyping via cytological preparation or flow cytometry require a large amount of tissue (such as blood) samples, which can be damaging or fatal to the fishes. Here, we developed novel microsatellite markers (SSR markers) from M. anguillicaudatus and show that they can effectively discriminate ploidy using samples collected in a minimally invasive way. Specifically, we generated whole genome transcriptomes from multiple M. anguillicaudatus using the Illumina paired-end sequencing. Approximately 150 million raw reads were assembled into 76,544 non-redundant unigenes. A total of 8,194 potential SSR markers were identified. We selected 98 pairs with more than five tandem repeats for further assays. Out of 45 putative EST-SSR markers that successfully amplified and harbored polymorphism in diploids, 11 markers displayed high variability in tetraploids. We further demonstrate that a set of five EST-SSR markers selected from these are sufficient to distinguish ploidy levels, by first validating them on 69 reference specimens with known ploidy levels and then subsequently using fresh-collected 96 ploidy-unknown specimens. The results from EST-SSR markers are highly concordant with those from independent flow cytometry analysis. The novel EST-SSR markers developed here should facilitate genetic studies of polyploidy in the emerging model system M. anguillicaudatus.

Eusocial insects exhibit unparalleled levels of cooperation and dominate terrestrial ecosystems. The success of eusocial insects stems from the presence of specialized castes that undertake distinct tasks. We investigated whether the evolutionary transition to societies with discrete castes was associated with changes in protein evolution. We predicted that proteins with caste-biased gene expression would evolve rapidly due to reduced antagonistic pleiotropy. We found that queen-biased proteins of the honeybee Apis mellifera did indeed evolve rapidly, as predicted. However, worker-biased proteins exhibited slower evolutionary rates than queen-biased or nonbiased proteins. We suggest that distinct selective pressures operating on caste-biased genes, rather than a general reduction in pleiotropy, explain the observed differences in evolutionary rates. Our study highlights, for the first time, the interaction between highly social behavior and dynamics of protein evolution.

A growing body of evidence suggests that DNA methylation is functionally divergent among different taxa. The recently discovered functional methylation system in the honeybee Apis mellifera presents an attractive invertebrate model system to study evolution and function of DNA methylation. In the honeybee, DNA methylation is mostly targeted toward transcription units (gene bodies) of a subset of genes. Here, we report an intriguing covariation of length and epigenetic status of honeybee genes. Hypermethylated and hypomethylated genes in honeybee are dramatically different in their lengths for both exons and introns. By analyzing orthologs in Drosophila melanogaster, Acyrthosiphon pisum, and Ciona intestinalis, we show genes that were short and long in the past are now preferentially situated in hyper- and hypomethylated classes respectively, in the honeybee. Moreover, we demonstrate that a subset of high-CpG genes are conspicuously longer than expected under the evolutionary relationship alone and that they are enriched in specific functional categories. We suggest that gene length evolution in the honeybee is partially driven by evolutionary forces related to regulation of gene expression, which in turn is associated with DNA methylation. However, lineage-specific patterns of gene length evolution suggest that there may exist additional forces underlying the observed interaction between DNA methylation and gene lengths in the honeybee.

There are more than 125 species of extant New World monkeys (Primates: Platyrrhini) found in approximately 15 genera. The phylogenetic relationships of these neotropical primates have been extensively studied from a molecular perspective. While these studies have been successful at inferring many of the relationships within the platyrrhines, key questions remain. The current study provides a framework for using non-genic, non-coding markers in comparative primate phylogenomic studies in species whose genomes are not yet scheduled for complete sequencing. A random genomic shotgun library was generated from the nocturnal Owl monkey Aotus lemurinus. Eleven unlinked, non-coding, non-genic, non-repetitive, nuclear DNA markers derived from this library were sequenced in at least one representative species of every platyrrhine genus. The combined sequence from these markers yielded a 7.7 kb multiple sequence alignment of 22 taxa. We analyzed these markers independently and combined with a 10 kb dataset consisting of "traditional," previously published markers located within or directly adjacent to genes. Parsimony, maximum likelihood, and Bayesian analysis converged on a single topology for the platyrrhine generic relationships. Notably, we confidently inferred that Pitheciidae is the sister taxon to the other two platyrrhine families (Cebidae, Atelidae). This relationship is supported by high values of branch support as well as topology tests. Additionally, Aotus formed a sister taxon to a clade comprising Cebus and Saimiri. With a fully resolved platyrrhine phylogeny in place it is now possible to design and test hypotheses regarding the evolution and diversification of platyrrhine phenotypes and life histories.

Two patterns from large-scale DNA sequence data have been put forward as evidence that speciation between humans and chimpanzees was complex, involving hybridization and strong selection. First, divergence between humans and chimpanzees varies considerably across the autosomes. Second, divergence between humans and chimpanzees (but not gorillas) is markedly lower on the X chromosome. Here, we describe how simple speciation and neutral molecular evolution explain both patterns. In particular, the wide range in autosomal divergence is consistent with stochastic variation in coalescence times in the ancestral population; and the lower human-chimpanzee divergence on the X chromosome is consistent with species differences in the strength of male-biased mutation caused by differences in mating system. We also highlight two further patterns of divergence that are problematic for the complex speciation model. Our conclusions raise doubts about complex speciation between humans and chimpanzees.

Inactivation of mismatch repair (MMR) has been shown to increase the accumulation of spontaneous mutations and frequency of recombination for diverse pathogenic bacteria. Currently, little is known regarding the role of mutator phenotypes for the diversification of natural populations of opportunistic human pathogens in marine environments. In this study, a higher frequency of mutators was detected among V. parahaemolyticus strains obtained from environmental sources compared with clinical sources. Inactivation of the MMR gene mutS caused increased antibiotic resistance and phase variation resulting in translucent colony morphologies. Increased nucleotide diversity in mutS and rpoB alleles from mutator compared with wild-type strains indicated a significant contribution of the mutator phenotype to the evolution of select genes. The results of this study indicate that the inactivation of MMR in V. parahaemolyticus leads to increased genetic and phenotypic diversity. This study is the first to report a higher frequency of natural mutators among Vibrio environmental strains and to provide evidence that inactivation of MMR increases the diversity of V. parahaemolyticus.

On the mechanistic level, epigenetic modifications include two major components: DNA methylation and histone tail modifications. Explicit evolutionary models of how epigenetic modifications evolve are currently lacking. Nevertheless, there are several examples of epigenetics having a significant evolutionary impact. A classic example is the impact of DNA methylation on mutational processes and genome sequence composition. More recently, evolution of mammalian recombination hotspots and transposable element methylation and gene expression in plants provide two solid examples where evolutionary theories incorporating epigenetic effects can be tested. The era of comparative epigenomics has already arrived.

The molecular clock hypothesis is a central concept in molecular evolution and has inspired much research into why evolutionary rates vary between and within genomes. In the age of modern comparative genomics, understanding the neutral genomic molecular clock occupies a critical place. It has been demonstrated that molecular clocks run differently between closely related species, and generation time is an important determinant of lineage specific molecular clocks. Moreover, it has been repeatedly shown that regional molecular clocks vary even within a genome, which should be taken into account when measuring evolutionary constraint of specific genomic regions. With the availability of a large amount of genomic sequence data, new insights into the patterns and causes of variation in molecular clocks are emerging. In particular, factors such as nucleotide composition, molecular origins of mutations, weak selection and recombination rates are important determinants of neutral genomic molecular clocks.

DNA methylation is one of the most important phenomena in epigenetics. Because the DNA methylation process regulates gene expression without altering the corresponding DNA sequences, it can explain many biological processes such as onset of cancer, aging, and differentiation of species. Therefore, detecting differentially methylated region (DMR) between groups has been considered one of important analyses in epigenetics. For the datasets generated by the methylated DNA immunoprecipitation (MeDIP) technique, DMR analysis has employed the previous methods developed for the microarray-based gene expression dataset. However, for the datasets generated by more recent bisulfite-sequencing (BS-seq) technologies, totally different approaches need to be used. Unfortunately, only a few methods are available for DMR analysis among which the Fisher's exact test is most widely used. We first propose a new DMR test for BS-seq data. Our approach employs Cochran-Mantel-Hantzel (CMH) test and extend it by incorporating correlation structures between adjacent CpG cytosine sites. Then we apply it to a real dataset of two subspecies of honeybee for detecting DMR regions.

The genome of the white-throated sparrow (Zonotrichia albicollis) contains an inversion polymorphism on chromosome 2 that is linked to predictable variation in a suite of phenotypic traits including plumage color, aggression, and parental behavior. Differences in gene expression between the two color morphs, which represent the two common inversion genotypes (ZAL2/ZAL2 and ZAL2/ZAL2m), are therefore of potential interest toward understanding the molecular underpinnings of these phenotypes. To identify genes that are differentially expressed between the two morphs and correlated with behavior, we quantified both behavior and brain gene expression in a population of free-living white-throated sparrows. We quantified behavioral responses to simulated territorial intrusions (STIs) early during the breeding season. In the same birds, we then performed a transcriptome-wide analysis of gene expression in two regions, the medial amygdala and hypothalamus. Both regions are part of a ‘social behavior network’, which is rich in steroid hormone receptors and previously linked with territorial behavior. Using network analyses, we identified modules of genes that were correlated with both morph and STI-induced singing behavior. The majority of these genes were located within the inversion, demonstrating the profound effect the inversion has on the expression of genes captured by the rearrangement. Gene pathway analyses revealed that in the medial amygdala, the most prominent pathways were those related to steroid hormone receptor activity. Within these pathways, the only gene encoding such a receptor was ESR1 (estrogen receptor alpha). Our results thus suggest that ESR1 and related genes are important for behavioral differences between the morphs.

Epigenetic changes can occur due to extracellular environmental conditions. Consequently, epigenetic mechanisms can play an intermediate role to translate environmental signals to intracellular changes. Such a role might be particularly important in plants, which often show strong local adaptation and have the potential for heritable epigenetic states. However, little is currently known about the role of epigenetic variation in the ecological mechanisms of adaptation. Here, we used multivariate redundancy analyses to examine genomewide associations between DNA methylation polymorphisms and climate variation in two independent panels of Arabidopsis accessions, including 122 Eurasian accessions as well as in a regional panel of 148 accessions in Sweden. At the single-nucleotide methylation level, climate and space (geographic spatial structure) explain small yet significant amount of variation in both panels. On the other hand, when viewed in a context of genomic clusters of methylated and unmethylated cytosines, climate and space variables explain much greater amounts of variation in DNA methylation than those explained by variation at the single-nucleotide level. We found that the single-nucleotide methylation polymorphisms with the strongest associations with climate were enriched in transposable elements and in potentially RNA-directed methylation contexts. When viewed in the context of genomic clusters, variation of DNA methylation at different sequence contexts exhibit distinctive segregation along different axes of variation in the redundancy analyses. Genomewide methylation showed much stronger associations with climate within the regional panel (Sweden) compared to the global (Eurasia). Together, these findings indicate that genetic and epigenetic variation across the genome may play a role in response to climate conditions and local adaptation.

Understanding the fundamental dynamics of epigenome variation during normal aging is critical for elucidating key epigenetic alterations that affect development, cell differentiation and diseases. Advances in the field of aging and DNA methylation strongly support the aging epigenetic drift model. Although this model aligns with previous studies, the role of other epigenetic marks, such as histone modification, as well as the impact of sampling specific CpGs, must be evaluated. Ultimately, it is crucial to investigate how all CpGs in the human genome change their methylation with aging in their specific genomic and epigenomic contexts. Here, we analyze whole genome bisulfite sequencing DNA methylation maps of brain frontal cortex from individuals of diverse ages. Comparisons with blood data reveal tissue-specific patterns of epigenetic drift. By integrating chromatin state information, divergent degrees and directions of aging-associated methylation in different genomic regions are revealed. Whole genome bisulfite sequencing data also open a new door to investigate whether adjacent CpG sites exhibit coordinated DNA methylation changes with aging. We identified significant 'aging-segments', which are clusters of nearby CpGs that respond to aging by similar DNA methylation changes. These segments not only capture previously identified aging-CpGs but also include specific functional categories of genes with implications on epigenetic regulation of aging. For example, genes associated with development are highly enriched in positive aging segments, which are gradually hyper-methylated with aging. On the other hand, regions that are gradually hypo-methylated with aging ('negative aging segments') in the brain harbor genes involved in metabolism and protein ubiquitination. Given the importance of protein ubiquitination in proteome homeostasis of aging brains and neurodegenerative disorders, our finding suggests the significance of epigenetic regulation of this posttranslational modification pathway in the aging brain. Utilizing aging segments rather than individual CpGs will provide more comprehensive genomic and epigenomic contexts to understand the intricate associations between genomic neighborhoods and developmental and aging processes. These results complement the aging epigenetic drift model and provide new insights.

Populations of honey bees are declining throughout the world, with US beekeepers losing 30% of their colonies each winter. Though multiple factors are driving these colony losses, it is increasingly clear that viruses play a major role. However, information about the molecular mechanisms mediating antiviral immunity in honey bees is surprisingly limited. Here, we examined the transcriptional and epigenetic (DNA methylation) responses to viral infection in honey bee workers. One-day old worker honey bees were fed solutions containing Israeli Acute Paralysis Virus (IAPV), a virus which causes muscle paralysis and death and has previously been associated with colony loss. Uninfected control and infected, symptomatic bees were collected within 20–24 hours after infection. Worker fat bodies, the primary tissue involved in metabolism, detoxification and immune responses, were collected for analysis. We performed transcriptome- and bisulfite-sequencing of the worker fat bodies to identify genome-wide gene expression and DNA methylation patterns associated with viral infection. There were 753 differentially expressed genes (FDR, Author Summary Honey bees are a critical pollinator of a wide variety of agricultural crops, but beekeepers experience heavy annual losses of honey bee colonies. Several factors are associated with colony losses, in particular infection with Israeli Acute Paralysis Virus (IAPV). Despite the importance of viruses to honey bee health, our understanding of the molecular mechanisms mediating host-pathogen interactions is limited. Here, we characterized the epigenomic and transcriptomic responses of honey bees to short term (

The role of sequence divergence in functional divergence of duplicate genes is a topic of great interest. In this study, we compare the numbers of amino acid substitutions in each sequence since two yeast duplicates diverged, using a preduplication ancestral outgroup. Using this strategy, we explored the relationship between sequence divergence and functional divergence between duplicate partners. We show that the degree of relative functional asymmetry between duplicate proteins is proportional to the relative sequence divergence between them. Furthermore, of the two duplicates, the copy closer to their ancestral sequence (fewer number of amino acid substitutions) interacts with more proteins and affects fitness more severely when deleted. Therefore, asymmetric sequence divergence between duplicates is correlated with asymmetric functional divergence and may underlie the duplicate's role in genetic robustness against mutations. Among the functional traits considered, protein abundance appears to have the strongest correlation with the nonsynonymous divergence between duplicates. Taken together with the results from whole-genome analyses, our results indicate that within-species duplicates are subject to the same evolutionary force that acts on interspecific sequence and functional divergence. In particular, we detect signs of purifying selection on the more slowly evolving duplicate.

A recent theory suggesting that genome size and complexity can increase as a passive consequence of small effective population size has generated much controversy. In this article, we demonstrate that freshwater fish species, which have smaller effective population sizes than marine fish species, have larger genomes. We show that genome size is negatively correlated with genetic variability, independent of phylogeny, body size and generation time. Genome duplication is also observed predominantly in freshwater fish. These results suggest that the raw materials of complexity originate under conditions of reduced selection efficiency.

We have investigated patterns of within-species polymorphism and between-species divergence for synonymous and nonsynonymous variants at a set of autosomal and X-linked loci of Drosophila miranda. D. pseudoobscura and D. affinis were used for the between-species comparisons. The results suggest the action of purifying selection on nonsynonymous, polymorphic variants. Among synonymous polymorphisms, there is a significant excess of synonymous mutations from preferred to unpreferred codons and of GC to AT mutations. There was no excess of GC to AT mutations among polymorphisms at noncoding sites. This suggests that selection is acting to maintain the use of preferred codons. Indirect evidence suggests that biased gene conversion in favor of GC base pairs may also be operating. The joint intensity of selection and biased gene conversion, in terms of the product of effective population size and the sum of the selection and conversion coefficients, was estimated to be ∼0.65.

Recent studies indicated that recombination is strongly mutagenic. In particular, data from the mouse pseudoautosomal boundary (PAB) suggested that locally intensive recombination increased the nucleotide substitution rate by more than 100-fold and greatly increased the GC content. Here we study the rates of nucleotide substitution in eight introns of the human and great ape XG gene, which spans the boundary between the pseudoautosomal region 1 (PAR1) and the X-specific region. Contrary to what is expected under the above hypothesis, our sequence data from humans and great apes reveal that the PAR1 introns of XG have actually evolved slightly slower than X-specific introns. Only when a New World monkey was compared with hominoids were the rates slightly increased in the PAR1 introns. In terms of base composition, although the intergenic regions of the human PAR1 show a significant increase of G and C nucleotides, the base composition of the surveyed PAR1 introns is similar to that of the X-specific introns. Direct and indirect evidence indicates that the recombination rate is, indeed, much higher in PAR1 introns than in X-specific introns, and that the present PAB has persisted since the common ancestor of hominoids. Therefore, the mutagenic effect of recombination is far weaker than previously proposed, at least in hominoid PABs.

Two longstanding issues on the molecular clock hypothesis are studied in this article. First, is there a global molecular clock in mammals? Although many authors have observed unequal rates of nucleotide substitution among mammalian lineages, some authors have proposed a global clock for all eutherians, i.e., a single global rate of 2.2 x 10(-9) substitutions per nucleotide site per year. We reexamine this issue using noncoding, nonrepetitive DNA from Old World monkeys (OWMs), chimpanzee, and human. First, using the minimal date of 6 MYA for the human-chimpanzee divergence and more than 2.5 million base pairs of genomic sequences from human and chimpanzee, we estimate a maximal rate of 0.99 x 10(-9) for noncoding, nonrepetitive genomic regions for these two species. This estimate is less than half of the proposed global rate and much smaller than the commonly used rate (3.5 x 10(-9)) for eutherians. Further, using a minimal date of 23 MYA for the human-OWM divergence, we estimate a maximal rate of 1.5 x 10(-9) for both introns and fourfold degenerate sites in humans and OWMs. In addition, with the New World monkey (NWM) lineage as an outgroup, we estimate that the rate of substitution in introns is 30% higher in the OWM lineage than in the human lineage. Clearly, there is no global molecular clock in eutherians. Second, although many studies have indicated considerable variation in the mutation rate among regions of the mammalian genome, a recent study proposed a uniform rate. Using new and existing intron sequence data from higher primates, we find significant rate variation among genomic regions and a positive correlation between the rate of substitution and the GC content, refuting the claim of a uniform rate.

Epigenetic mechanisms traditionally have been studied in the domains of development and disease, but they may also play important roles in ecological and evolutionary processes. In this article, we revisit historical as well as recent studies that indicate significant impacts of epigenetic processes on evolution. Our main focus is DNA methylation, which is a prevalent chemical modification of genomic DNA. First, it has been long known that DNA methylation acts as a major mutational facilitator in animal genomes and influences nucleotide compositions of genomes. More recently, genome-wide analyses have demonstrated that the current levels of DNA methylation can be predicted from the evolutionary signatures of DNA methylation, indicating that these two processes are intimately correlated. Indeed, the recent explosive growth in the knowledge of genomic DNA methylation in wide-ranging taxa has revealed that patterns of DNA methylation are surprisingly conserved across deep phylogenies. Interestingly, comparative analyses of humans and closely related primate species show that genomic regions that do show evolutionary divergence of DNA methylation are enriched for developmental and tissue specializations. A key question is how epigenetic patterns transmit between generations and impact evolutionary dynamics. On the one hand, some studies report direct transmissions of epigenetic features to the next generation. On the other hand, it is becoming clear that genomic sequence variants exist that encode and presumably regulate distinctive epigenetic patterns. For instance, numerous single-nucleotide polymorphisms that affect DNA-methylation patterns have been discovered in human populations. These studies begin to unveil a dynamic interplay between genomic and epigenomic factors across long and short evolutionary timescales.

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.

In Drosophila miranda, a chromosome fusion between the Y chromosome and the autosome corresponding to Muller’s element C has created a new sex chromosome system. The chromosome attached to the ancestral Y chromosome is transmitted paternally and hence is not exposed to crossing over. This chromosome, conventionally called the neo-Y, and the homologous neo-X chromosome display many properties of evolving sex chromosomes. We report here the transposition of the exuperantia1 (exu1) locus from a neo-sex chromosome to the ancestral X chromosome of D. miranda. Exu1 is known to have several critical developmental functions, including a male-specific role in spermatogenesis. The ancestral location of exu1 is conserved in the sibling species of D. miranda, as well as in a more distantly related species. The transposition of exu1 can be interpreted as an adaptive fixation, driven by a selective advantage conferred by its effect on dosage compensation. This explanation is supported by the pattern of within-species sequence variation at exu1 and the nearby exu2 locus. The implications of this phenomenon for genome evolution are discussed.

CpG islands (CGIs) are commonly used as genomic markers to study the patterns and regulatory consequences of DNA methylation. Interestingly, recent studies reveal a substantial diversity among CGIs: long and short CGIs, for example, exhibit contrasting patterns of gene expression complexity and nucleosome occupancy. Evolutionary origins of CGIs are also highly heterogeneous. In order to systematically evaluate potential diversities among CGIs and ultimately to illuminate the link between diversity of CGIs and their epigenetic variation, we analyzed the nucleotide-resolution DNA methylation maps (methylomes) of multiple cellular origins. We discover novel ‘clusters’ of CGIs according to their patterns of DNA methylation; the stably hypomethylated CGI cluster (cluster I), sperm-hypomethylated CGI cluster (cluster II), and variably methylated CGI cluster (cluster III). These epigenomic CGI clusters are strikingly distinct at multiple biological features including genomic, evolutionary, and functional characteristics. At the genomic level, the stably hypomethylated CGI cluster tends to be longer and harbors many more CpG dinucleotides than those in other clusters. They are also frequently associated with promoters, while CGI clusters II and III mostly reside in intragenic or intergenic regions and exhibit highly tissue-specific DNA methylation. Functional ontology terms and transcriptional profiles co-vary with CGI clusters, indicating that the regulatory functions of CGIs are tightly linked to their heterogeneity. Finally, CGIs associated with distinctive biological processes, such as diseases, aging, and imprinting, occur disproportionately across CGI clusters. These new findings provide an effective means to combine existing knowledge on CGIs into a genomic context while bringing new insights that elucidate the significance of DNA methylation across different biological conditions and demography.

Background Whole genome sequencing of bisulfite converted DNA (‘methylC-seq’) method provides comprehensive information of DNA methylation. An important application of these whole genome methylation maps is classifying each position as a methylated versus non-methylated nucleotide. A widely used current method for this purpose, the so-called binomial method, is intuitive and straightforward, but lacks power when the sequence coverage and the genome-wide methylation level are low. These problems present a particular challenge when analyzing sparsely methylated genomes, such as those of many invertebrates and plants. Results We demonstrate that the number of sequence reads per position from methylC-seq data displays a large variance and can be modeled as a shifted negative binomial distribution. We also show that DNA methylation levels of adjacent CpG sites are correlated, and this similarity in local DNA methylation levels extends several kilobases. Taking these observations into account, we propose a new method based on Bayesian classification to infer DNA methylation status while considering the neighborhood DNA methylation levels of a specific site. We show that our approach has higher sensitivity and better classification performance than the binomial method via multiple analyses, including computational simulations, Area Under Curve (AUC) analyses, and improved consistencies across biological replicates. This method is especially advantageous in the analyses of sparsely methylated genomes with low coverage. Conclusions Our method improves the existing binomial method for binary methylation calls by utilizing a posterior odds framework and incorporating local methylation information. This method should be widely applicable to the analyses of methylC-seq data from diverse sparsely methylated genomes. Bis-Class and example data are provided at a dedicated website (http://bibs.snu.ac.kr/software/Bisclass). Electronic supplementary material The online version of this article (doi:10.1186/1471-2164-15-608) contains supplementary material, which is available to authorized users.

Kin selection is a fundamentally important process that affects the evolution of social behaviours. The genomics revolution now provides the opportunity to test kin selection theory using genomic data. In this commentary, we discuss previous studies that explored the link between kin selection and patterns of variation within the genome. We then present a new theory aimed at understanding the evolution of genes involved in the development of social insects. Specifically, we investigate caste-antagonistic pleiotropy, which occurs when the phenotypes of distinct castes are optimized by different genotypes at a single locus. We find that caste-antagonistic pleiotropy leads to narrow regions where polymorphism can be maintained. Furthermore, multiple mating by queens reduces the region in which worker-favoured alleles fix, which suggests that multiple mating impedes worker caste evolution. We conclude by discussing ways to test these and other facets of kin selection using newly emerging genomic data.

The "developmental hourglass'' describes a pattern of increasing morphological divergence towards earlier and later embryonic development, separated by a period of significant conservation across distant species (the "phylotypic stage''). Recent studies have found evidence in support of the hourglass effect at the genomic level. For instance, the phylotypic stage expresses the oldest and most conserved transcriptomes. However, the regulatory mechanism that causes the hourglass pattern remains an open question. Here, we use an evolutionary model of regulatory gene interactions during development to identify the conditions under which the hourglass effect can emerge in a general setting. The model focuses on the hierarchical gene regulatory network that controls the developmental process, and on the evolution of a population under random perturbations in the structure of that network. The model predicts, under fairly general assumptions, the emergence of an hourglass pattern in the structure of a temporal representation of the underlying gene regulatory network. The evolutionary age of the corresponding genes also follows an hourglass pattern, with the oldest genes concentrated at the hourglass waist. The key behind the hourglass effect is that developmental regulators should have an increasingly specific function as development progresses. Analysis of developmental gene expression profiles from Drosophila melanogaster and Arabidopsis thaliana provide consistent results with our theoretical predictions.

Background Fig pollinating wasps form obligate symbioses with their fig hosts. This mutualism arose approximately 75 million years ago. Unlike many other intimate symbioses, which involve vertical transmission of symbionts to host offspring, female fig wasps fly great distances to transfer horizontally between hosts. In contrast, male wasps are wingless and cannot disperse. Symbionts that keep intimate contact with their hosts often show genome reduction, but it is not clear if the wide dispersal of female fig wasps will counteract this general tendency. We sequenced the genome of the fig wasp Ceratosolen solmsi to address this question. Results The genome size of the fig wasp C. solmsi is typical of insects, but has undergone dramatic reductions of gene families involved in environmental sensing and detoxification. The streamlined chemosensory ability reflects the overwhelming importance of females finding trees of their only host species, Ficus hispida, during their fleeting adult lives. Despite long-distance dispersal, little need exists for detoxification or environmental protection because fig wasps spend nearly all of their lives inside a largely benign host. Analyses of transcriptomes in females and males at four key life stages reveal that the extreme anatomical sexual dimorphism of fig wasps may result from a strong bias in sex-differential gene expression. Conclusions Our comparison of the C. solmsi genome with other insects provides new insights into the evolution of obligate mutualism. The draft genome of the fig wasp, and transcriptomic comparisons between both sexes at four different life stages, provide insights into the molecular basis for the extreme anatomical sexual dimorphism of this species.

The development and evolution of organisms is heavily influenced by their environment. Thus, understanding the historical biogeography of taxa can provide insights into their evolutionary history, adaptations and trade-offs realized throughout time. In the present study we have taken a phylogenomic approach to infer New World monkey phylogeny, upon which we have reconstructed the biogeographic history of extant platyrrhines. In order to generate sufficient phylogenetic signal within the New World monkey clade, we carried out a large-scale phylogenetic analysis of approximately 40 kb of non-genic genomic DNA sequence in a 36 species subset of extant New World monkeys. Maximum parsimony, maximum likelihood and Bayesian inference analysis all converged on a single optimal tree topology. Divergence dating and biogeographic analysis reconstruct the timing and geographic location of divergence events. The ancestral area reconstruction describes the geographic locations of the last common ancestor of extant platyrrhines and provides insight into key biogeographic events occurring during platyrrhine diversification. Through these analyses we conclude that the diversification of the platyrrhines took place concurrently with the establishment and diversification of the Amazon rainforest. This suggests that an expanding rainforest environment rather than geographic isolation drove platyrrhine diversification.

Epigenetic information is an important mediator of the relationship between genotype and phenotype in eukaryotic organisms. One of the most important and widely conserved forms of epigenetic information is the methylation of genes. However, the function of intragenic DNA methylation remains poorly understood. The goal of this study was to gain greater understanding of the nature of intragenic methylation by determining its role in the multilayered epigenetic landscape of insects. We first investigated the evolutionary lability of DNA methylation by examining whether methylation patterns were conserved in the fire ant and honey bee. We found that DNA methylation was targeted to largely overlapping sets of orthologs in both species. Next, we compared intragenic DNA methylation levels in the fire ant and honey bee to comprehensive epigenetic and gene-regulatory data from Drosophila melanogaster orthologs. We observed striking evidence of a conserved association between DNA methylation in fire ants and honey bees, and several active histone modifications, constitutive gene expression, and "broad" promoter architecture in D. melanogaster. Overall, our study illustrates that DNA methylation is a single component of a conserved, integrated, multilayered epigenetic and regulatory landscape in insect genomes.

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.

Half a century ago, when the field of molecular evolution did not even exist, Morris Goodman analyzed profiles of immunological interactions between species and reached the following two remarkable conclusions: first, protein evolution slowed down in the human lineage compared to other primate lineages; second, this slowdown was more pronounced for proteins whose functions were likely to be neutral. It took several decades of research to fully grasp these ideas and document the pattern of hominoid rate slowdown. Along the way, studies of hominoid rate slowdown led to major progresses in understanding determinants of neutral molecular evolution, which in turn is used to calibrate rates of adaptive evolution. Furthermore, the growing knowledge on the origin of mutations provides a basis for understanding differential evolutionary rates between sex chromosomes and autosomes, which has deep implications for inferring human evolutionary histories, and other aspects of molecular evolution. Primate genomics in particular stand to provide critical information in these pursuits, due to the abundance of genomic data, relatively rich documentation of life history traits, and several model systems, including our own species.

DNA methylation of transcription units (gene bodies) occurs in the genomes of many animal and plant species. Phylogenetic persistence of gene body methylation implies biological significance; yet, the functional roles of gene body methylation remain elusive. In this study, we analyzed methylation levels of orthologs from four distantly related invertebrate species, including the honeybee, silkworm, sea squirt, and sea anemone. We demonstrate that in all four species, gene bodies distinctively cluster to two groups, which correspond to high and low methylation levels. This pattern resembles that of sequence composition arising from the mutagenetic effect of DNA methylation. In spite of this effect, our results show that protein sequences of genes targeted by high levels of methylation are conserved relative to genes lacking methylation. Our investigation identified many genes that either gained or lost methylation during the course of invertebrate evolution. Most of these genes appear to have lost methylation in the insect lineages we investigated, particularly in the honeybee. We found that genes that are methylated in all four invertebrate taxa are enriched for housekeeping functions related to transcription and translation, whereas the loss of DNA methylation occurred in genes whose functions include cellular signaling and reproductive processes. Overall, our study helps to illuminate the functional significance of gene body methylation and its impacts on genome evolution in diverse invertebrate taxa.

6 páginas, 5 figuras -- PAGS nros. 487-492, Grbic, Miodrag et.al., The spider mite Tetranychus urticae is a cosmopolitan agricultural pest with an extensive host plant range and an extreme record of pesticide resistance. Here we present the completely sequenced and annotated spider mite genome, representing the first complete chelicerate genome. At 90 megabases T. urticae has the smallest sequenced arthropod genome. Compared with other arthropods, the spider mite genome shows unique changes in the hormonal environment and organization of the Hox complex, and also reveals evolutionary innovation of silk production. We find strong signatures of polyphagy and detoxification in gene families associated with feeding on different hosts and in new gene families acquired by lateral gene transfer. Deep transcriptome analysis of mites feeding on different plants shows how this pest responds to a changing host environment. The T. urticae genome thus offers new insights into arthropod evolution and plant–herbivore interactions, and provides unique opportunities for developing novel plant protection strategies, M.G. and V.G. acknowledge support from NSERC Strategic Grant STPGP 322206-05, Marie Curie Incoming International Fellowship, OECD Co-operative Research Programme: Biological resource management for Sustainable Agricultural Systems JA00053351, and Ontario Research Fund–Global Leadership in Genomics and Life Sciences GL2-01-035. The genome and transcriptome sequencing projects were funded by the Government of Canada through Genome Canada and the Ontario Genomics Institute (OGI-046), JGI Community Sequencing Program grant 777506 to M.G., a University of Utah SEED grant (to R.M.C.), and National Science Foundation (NSF) grant 0820985 (to R.M.C., Principal Investigator L. Sieburth); work conducted by the US Department of Energy Joint Genome Institute is supported by the Office of Science of the US Department of Energy under contract No. DE-AC02-05CH11231

DNA methylation is a fundamental epigenetic mark known to have wide-ranging effects on gene regulation in a variety of animal taxa. Comparative genomic analyses can help elucidate the function of DNA methylation by identifying conserved features of methylated genes and other genomic regions. In this study, we used computational approaches to distinguish genes marked by heavy methylation from those marked by little or no methylation in the pea aphid, Acyrthosiphon pisum. We investigated if these two classes had distinct evolutionary histories and functional roles by conducting comparative analysis with the honeybee, Apis (Ap.) mellifera. We found that highly methylated orthologs in A. pisum and Ap. mellifera exhibited greater conservation of methylation status, suggesting that highly methylated genes in ancestral species may remain highly methylated over time. We also found that methylated genes tended to show different rates of evolution than unmethylated genes. In addition, genes targeted by methylation were enriched for particular biological processes that differed from those in relatively unmethylated genes. Finally, methylated genes were preferentially ubiquitously expressed among alternate phenotypes in both species, whereas genes lacking signatures of methylation were preferentially associated with condition-specific gene regulation expression. Overall, our analyses support a conserved role for DNA methylation in insects with comparable methylation systems.

DNA methylation is implicated in modulating gene transcription in a wide array of taxa. However, a complete understanding of the specific functions of DNA methylation in diverse organisms requires study of the evolutionary patterns of genome methylation. Unfortunately, investigating DNA methylation through empirical means remains challenging due to the transient nature of DNA methylation, the paucity of model systems in which to study genome methylation, and the costs associated with experimental methodology. Here we review how computational methods can be used to complement experimental approaches to further our understanding of DNA methylation in animals. For instance, comparative analyses of the molecular machinery involved in DNA methylation have been informative in revealing the dynamics of genome methylation in many organisms. In addition, analyses of specific genomic signatures of DNA methylation have furthered our understanding of the patterns of methylation within species. Finally, insight into the role of DNA methylation has resulted from computational methods used to identify specific sets of methylated genes. We suggest that an understanding of the evolution of genomic DNA methylation can be most readily achieved by integrating computational and empirical methods.

The recent, unexpected discovery of a functional DNA methylation system in the genome of the social bee Apis mellifera underscores the potential importance of DNA methylation in invertebrates. The extent of genomic DNA methylation and its role in A. mellifera remain unknown, however. Here we show that genes in A. mellifera can be divided into 2 distinct classes, one with low-CpG dinucleotide content and the other with high-CpG dinucleotide content. This dichotomy is explained by the gradual depletion of CpG dinucleotides, a well-known consequence of DNA methylation. The loss of CpG dinucleotides associated with DNA methylation also may explain the unusual mutational patterns seen in A. mellifera that lead to AT-rich regions of the genome. A detailed investigation of this dichotomy implicates DNA methylation in A. mellifera development. High-CpG genes, which are predicted to be hypomethylated in germlines, are enriched with functions associated with developmental processes, whereas low-CpG genes, predicted to be hypermethylated in germlines, are enriched with functions associated with basic biological processes. Furthermore, genes more highly expressed in one caste than another are overrepresented among high-CpG genes. Our results highlight the potential significance of epigenetic modifications, such as DNA methylation, in developmental processes in social insects. In particular, the pervasiveness of DNA methylation in the genome of A. mellifera provides fertile ground for future studies of phenotypic plasticity and genomic imprinting.

We analysed over 8 million base pairs of bacterial artificial chromosome-based sequence alignments of four Old World monkeys and the human genome. Our findings are as follows. (i) Genomic divergences among several Old World monkeys mirror those between well-studied hominoids. (ii) The X-chromosome evolves slower than autosomes, in accord with ‘male-driven evolution’. However, the degree of male mutation bias is lower in Old World monkeys than in hominoids. (iii) Evolutionary rates vary significantly between lineages. The baboon branch shows a particularly slow molecular evolution. Thus, lineage-specific evolutionary rate variation is a common theme of primate genome evolution. (iv) In contrast to the overall pattern, mutations originating from DNA methylation exhibit little variation between lineages. Our study illustrates the potential of primates as a model system to investigate genome evolution, in particular to elucidate molecular mechanisms of substitution rate variation.

Functional RNAs (fRNAs) are being recognized as an important regulatory component in biological processes. Interestingly, recent computational studies suggest that the number and biological significance of functional RNAs within coding regions (coding fRNAs) may have been underestimated. We hypothesized that such coding fRNAs will impose additional constraint on sequence evolution because the DNA primary sequence has to simultaneously code for functional RNA secondary structures on the messenger RNA in addition to the amino acid codons for the protein sequence. To test this prediction, we first utilized computational methods to predict conserved fRNA secondary structures within multiple species alignments of Saccharomyces sensu strico genomes. We predict that as much as 5% of the genes in the yeast genome contain at least one functional RNA secondary structure within their protein-coding region. We then analyzed the impact of coding fRNAs on the evolutionary rate of protein-coding genes because a decrease in evolutionary rate implies constraint due to biological functionality. We found that our predicted coding fRNAs have a significant influence on evolutionary rates (especially at synonymous sites), independent of other functional measures. Thus, coding fRNA may play a role on sequence evolution. Given that coding regions of humans and flies contain many more predicted coding fRNAs than yeast, the impact of coding fRNAs on sequence evolution may be substantial in genomes of higher eukaryotes.

It is often stated that patterns of nonsynonymous rate variation among mammalian lineages are more irregular than expected or overdispersed under the neutral model, whereas synonymous sites conform to the neutral model. Here we reexamined genome-wide patterns of the variance to mean ratio, or index of dispersion (R), of substitutions in proteins from human, mouse, and dog. Contrary to the prevailing notion, we found that the mean index of dispersion for nonsynonymous sites of mammalian proteins is not significantly different from 1. We propose that earlier analyses were biased because the data included disproportionately more protein hormones, which tend to be more dispersed than genes in other functional categories. Synonymous sites exhibit greater degree of dispersion than nonsynonymous sites, although similar to earlier estimates and potentially due to errors associated with correction for multiple hits. Overall, our analysis identifies strong genome-wide generation-time effect and natural selection as important determinants of among-lineage variation of protein evolutionary rates. Furthermore, patterns of lineage-specific selective constraint are consistent with the nearly neutral model of molecular evolution.

When comparing protein-coding DNA sequences between species, it is useful to distinguish between synonymous (silent) substitutions and nonsynonymous (amino-acid changing) substitutions. Because the two types of change are under different intensities of selective constraints, comparing their rates provides insights into the mechanisms of evolution of the gene under study. Keywords: coding regions; nucleotide substitution; degenerate code; functional constraints; positive selection