Impacts of pretranscriptional DNA methylation, transcriptional transcription factor, and posttranscriptional microRNA regulations on protein evolutionary rate

MicroRNA in the Diagnosis and Treatment of Doxorubicin-Induced Cardiotoxicity

Doxorubicin (DOX), a broad-spectrum chemotherapy drug, is widely applied to the treatment of cancer; however, DOX-induced cardiotoxicity (DIC) limits its clinical therapeutic utility. However, it is difficult to monitor and detect DIC at an early stage using conventional detection methods. Thus, sensitive, accurate, and specific methods of diagnosis and treatment are important in clinical practice. MicroRNAs (miRNAs) belong to non-coding RNAs (ncRNAs) and are stable and easy to detect. Moreover, miRNAs are expected to become biomarkers and therapeutic targets for DIC; thus, there are currently many studies focusing on the role of miRNAs in DIC. In this review, we list the prominent studies on the diagnosis and treatment of miRNAs in DIC, explore the feasibility and difficulties of using miRNAs as diagnostic biomarkers and therapeutic targets, and provide recommendations for future research.

Genome-Wide DNA Methylation Signatures Predict the Early Asymptomatic Doxorubicin-Induced Cardiotoxicity in Breast Cancer

Chemotherapy with doxorubicin (DOX) may cause unpredictable cardiotoxicity. This study aimed to determine whether the methylation signature of peripheral blood mononuclear cells (PBMCs) prior to and after the first cycle of DOX-based chemotherapy could predict the risk of cardiotoxicity in breast cancer patients. Cardiotoxicity was defined as a decrease in left ventricular ejection fraction (LVEF) by >10%. DNA methylation of PBMCs from 9 patients with abnormal LVEF and 10 patients with normal LVEF were examined using Infinium HumanMethylation450 BeadChip. We have identified 14,883 differentially methylated CpGs at baseline and 18,718 CpGs after the first cycle of chemotherapy, which significantly correlated with LVEF status. Significant differentially methylated regions (DMRs) were found in the promoter and the gene body of SLFN12, IRF6 and RNF39 in patients with abnormal LVEF. The pathway analysis found enrichment for regulation of transcription, mRNA splicing, pathways in cancer and ErbB2/4 signaling. The preliminary results from this study showed that the DNA methylation profile of PBMCs may predict the risk of DOX-induced cardiotoxicity prior to chemotherapy. Further studies with larger cohorts of patients are needed to confirm these findings.

Histone H3K27 methylation perturbs transcriptional robustness and underpins dispensability of highly conserved genes in fungi

Epigenetic modifications are key regulators of gene expression and underpin genome integrity. Yet, how epigenetic changes affect the evolution and transcriptional robustness of genes remains largely unknown. Here, we show how the repressive histone mark H3K27me3 underpins the trajectory of highly conserved genes in fungi. We first performed transcriptomic profiling on closely related species of the plant pathogen Fusarium graminearum species complex. We determined transcriptional responsiveness of genes across environmental conditions to determine expression robustness. To infer evolutionary conservation, we used a framework of 23 species across the Fusarium genus including three species covered with histone methylation data. Gene expression variation is negatively correlated with gene conservation confirming that highly conserved genes show higher expression robustness. In contrast, genes marked by H3K27me3 do not show such associations. Furthermore, highly conserved genes marked by H3K27me3 encode smaller proteins, exhibit weaker codon usage bias, higher levels of hydrophobicity, show lower intrinsically disordered regions, and are enriched for functions related to regulation and membrane transport. The evolutionary age of conserved genes with H3K27me3 histone marks falls typically within the origins of the Fusarium genus. We show that highly conserved genes marked by H3K27me3 are more likely to be dispensable for survival during host infection. Lastly, we show that conserved genes exposed to repressive H3K27me3 marks across distantly related Fusarium fungi are associated with transcriptional perturbation at the microevolutionary scale. In conclusion, we show how repressive histone marks are entangled in the evolutionary fate of highly conserved genes across evolutionary timescales.

Molecular Profiling of DNA Methylation and Alternative Splicing of Genes in Skeletal Muscle of Obese Rabbits

DNA methylation and the alternative splicing of precursor messenger RNAs (pre-mRNAs) are two important genetic modification mechanisms. However, both are currently uncharacterized in the muscle metabolism of rabbits. Thus, we constructed the Tianfu black rabbit obesity model (obese rabbits fed with a 10% high-fat diet and control rabbits from 35 days to 70 days) and collected the skeletal muscle samples from the two groups for Genome methylation sequencing and RNA sequencing. DNA methylation data showed that the promoter regions of 599 genes and gene body region of 2522 genes had significantly differential methylation rates between the two groups, of which 288 genes had differential methylation rates in promoter and gene body regions. Analysis of alternative splicing showed 555 genes involved in exon skipping (ES) patterns, and 15 genes existed in differential methylation regions. Network analysis showed that 20 hub genes were associated with ubiquitinated protein degradation, muscle development pathways, and skeletal muscle energy metabolism. Our findings suggest that the two types of genetic modification have potential regulatory effects on skeletal muscle development and provide a basis for further mechanistic studies in the rabbit.

A FITM1-Related Methylation Signature Predicts the Prognosis of Patients With Non-Viral Hepatocellular Carcinoma

Although great progress has been made in treatment against hepatitis virus infection, the prognosis of hepatocellular carcinoma (HCC) remains unsatisfied. Therefore, there is an unmet need to explore biomarkers or prognostic models for monitoring non-viral hepatocellular carcinoma. Accumulating evidence indicates that DNA methylation participates in carcinogenesis of malignancies. In the present study, we analyzed 101 non-viral HCC patients from TCGA database to figure out methylation-driven genes (MDGs) that might get involved in non-viral HCC pathogenesis using MethyMix algorithm. Then we picked out 8 key genes out of 137 MDGs that could affect the overall survival (OS) of both methylation and expression level. Using PCA, Uni-variate, Multi-variate, and LASSO cox regression analyses, we confirmed the potential prognostic value of these eight epigenetic genes. Ultimately, combined with immunohistochemistry (IHC), ROC, OS, and GSEA analyses, fat storage-inducing transmembrane protein1 (FITM1) was identified as a novel tumor suppressor gene in non-viral HCC and an applicable FITM1-methylation-based signature was built in a training set and validated in a testing set. Briefly, our work provides several potential biomarkers, especially FITM1, as well as a new method for disease surveillance and treatment strategy development.

A-to-I RNA editing contributes to the persistence of predicted damaging mutations in populations

Adenosine-to-inosine (A-to-I) RNA editing is a very common co-/posttranscriptional modification that can lead to A-to-G changes at the RNA level and compensate for G-to-A genomic changes to a certain extent. It has been shown that each healthy individual can carry dozens of missense variants predicted to be severely deleterious. Why strongly detrimental variants are preserved in a population and not eliminated by negative natural selection remains mostly unclear. Here, we ask if RNA editing correlates with the burden of deleterious A/G polymorphisms in a population. Integrating genome and transcriptome sequencing data from 447 human lymphoblastoid cell lines, we show that nonsynonymous editing activities (prevalence/level) are negatively correlated with the deleteriousness of A-to-G genomic changes and positively correlated with that of G-to-A genomic changes within the population. We find a significantly negative correlation between nonsynonymous editing activities and allele frequency of A within the population. This negative editing-allele frequency correlation is particularly strong when editing sites are located in highly important genes/loci. Examinations of deleterious missense variants from the 1000 Genomes Project further show a significantly higher proportion of rare missense mutations for G-to-A changes than for other types of changes. The proportion for G-to-A changes increases with increasing deleterious effects of the changes. Moreover, the deleteriousness of G-to-A changes is significantly positively correlated with the percentage of editing enzyme binding motifs at the variants. Overall, we show that nonsynonymous editing is associated with the increased burden of G-to-A missense mutations in healthy individuals, expanding RNA editing in pathogenomics studies.

Demethylation and microRNA differential expression regulate plasma-induced improvement of chicken sperm quality

The sperm quality is a vital economical requisite of poultry production. Our previous study found non-thermal dielectric barrier discharge plasma exposure on fertilized eggs could increase the chicken growth and the male reproduction. However, it is unclear how plasma treatment regulates the reproductive capacity in male chickens. In this study, we used the optimal plasma treatment condition (2.81 W for 2 min) which has been applied on 3.5-day-incubated fertilized eggs in the previous work and investigated the reproductive performance in male chickens aged at 20 and 40 weeks. The results showed that plasma exposure increased sperm count, motility, fertility rate, and fertilization period of male chickens. The sperm quality-promoting effect of plasma treatment was regulated by the significant improvements of adenosine triphosphate production and testosterone level, and by the modulation of reactive oxygen species balance and adenosine monophosphate-activated protein kinase and mammalian target of rapamycin pathway in the spermatozoa. Additionally, the plasma effect suggested that DNA demethylation and microRNA differential expression (a total number of 39 microRNAs were up-regulated whereas 53 microRNAs down-regulated in the testis) regulated the increases of adenosine triphosphate synthesis and testosterone level for promoting the chicken sperm quality. This finding might be beneficial to elevate the fertilization rate and embryo quality for the next generation in poultry breeding.

Allele-specific DNA methylation of disease susceptibility genes in Japanese patients with inflammatory bowel disease

BACKGROUND

Inflammatory bowel disease (IBD) has an unknown etiology; however, accumulating evidence suggests that IBD is a multifactorial disease influenced by a combination of genetic and environmental factors. The influence of genetic variants on DNA methylation in cis and cis effects on expression have been demonstrated. We hypothesized that IBD susceptibility single-nucleotide polymorphisms (SNPs) regulate susceptibility gene expressions in cis by regulating DNA methylation around SNPs. For this, we determined cis-regulated allele-specific DNA methylation (ASM) around IBD susceptibility genes in CD4+ effector/memory T cells (Tem) in lamina propria mononuclear cells (LPMCs) in patients with IBD and examined the association between the ASM SNP genotype and neighboring susceptibility gene expressions.

METHODS

CD4+ effector/memory T cells (Tem) were isolated from LPMCs in 15 Japanese IBD patients (ten Crohn's disease [CD] and five ulcerative colitis [UC] patients). ASM analysis was performed by methylation-sensitive SNP array analysis. We defined ASM as a changing average relative allele score ([Formula: see text]) >0.1 after digestion by methylation-sensitive restriction enzymes. Among SNPs showing [Formula: see text] >0.1, we extracted the probes located on tag-SNPs of 200 IBD susceptibility loci and around IBD susceptibility genes as candidate ASM SNPs. To validate ASM, bisulfite-pyrosequencing was performed. Transcriptome analysis was examined in 11 IBD patients (seven CD and four UC patients). The relation between rs36221701 genotype and neighboring gene expressions were analyzed.

RESULTS

We extracted six candidate ASM SNPs around IBD susceptibility genes. The top of [Formula: see text] (0.23) was rs1130368 located on HLA-DQB1. ASM around rs36221701 ([Formula: see text] = 0.14) located near SMAD3 was validated using bisulfite pyrosequencing. The SMAD3 expression was significantly associated with the rs36221701 genotype (p = 0.016).

CONCLUSIONS

We confirmed the existence of cis-regulated ASM around IBD susceptibility genes and the association between ASM SNP (rs36221701) genotype and SMAD3 expression, a susceptibility gene for IBD. These results give us supporting evidence that DNA methylation mediates genetic effects on disease susceptibility.

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

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

Complex Patterns of Association between Pleiotropy and Transcription Factor Evolution

Pleiotropy has been claimed to constrain gene evolution but specific mechanisms and extent of these constraints have been difficult to demonstrate. The expansion of molecular data makes it possible to investigate these pleiotropic effects. Few classes of genes have been characterized as intensely as human transcription factors (TFs). We therefore analyzed the evolutionary rates of full TF proteins, along with their DNA binding domains and protein-protein interacting domains (PID) in light of the degree of pleiotropy, measured by the number of TF-TF interactions, or the number of DNA-binding targets. Data were extracted from the ENCODE Chip-Seq dataset, the String v 9.2 database, and the NHGRI GWAS catalog. Evolutionary rates of proteins and domains were calculated using the PAML CodeML package. Our analysis shows that the numbers of TF-TF interactions and DNA binding targets associated with constrained gene evolution; however, the constraint caused by the number of DNA binding targets was restricted to the DNA binding domains, whereas the number of TF-TF interactions constrained the full protein and did so more strongly. Additionally, we found a positive correlation between the number of protein-PIDs and the evolutionary rates of the protein-PIDs. These findings show that not only does pleiotropy associate with constrained protein evolution but the constraint differs by domain function. Finally, we show that GWAS associated TF genes are more highly pleiotropic : The GWAS data illustrates that mutations in highly pleiotropic genes are more likely to be associated with disease phenotypes.

Transcriptome sequencing reveals genome-wide variation in molecular evolutionary rate among ferns

BACKGROUND

Transcriptomics in non-model plant systems has recently reached a point where the examination of nuclear genome-wide patterns in understudied groups is an achievable reality. This progress is especially notable in evolutionary studies of ferns, for which molecular resources to date have been derived primarily from the plastid genome. Here, we utilize transcriptome data in the first genome-wide comparative study of molecular evolutionary rate in ferns. We focus on the ecologically diverse family Pteridaceae, which comprises about 10 % of fern diversity and includes the enigmatic vittarioid ferns-an epiphytic, tropical lineage known for dramatically reduced morphologies and radically elongated phylogenetic branch lengths. Using expressed sequence data for 2091 loci, we perform pairwise comparisons of molecular evolutionary rate among 12 species spanning the three largest clades in the family and ask whether previously documented heterogeneity in plastid substitution rates is reflected in their nuclear genomes. We then inquire whether variation in evolutionary rate is being shaped by genes belonging to specific functional categories and test for differential patterns of selection.

RESULTS

We find significant, genome-wide differences in evolutionary rate for vittarioid ferns relative to all other lineages within the Pteridaceae, but we recover few significant correlations between faster/slower vittarioid loci and known functional gene categories. We demonstrate that the faster rates characteristic of the vittarioid ferns are likely not driven by positive selection, nor are they unique to any particular type of nucleotide substitution.

CONCLUSIONS

Our results reinforce recently reviewed mechanisms hypothesized to shape molecular evolutionary rates in vittarioid ferns and provide novel insight into substitution rate variation both within and among fern nuclear genomes.

Extensive ceRNA-ceRNA interaction networks mediated by miRNAs regulate development in multiple rhesus tissues

Crosstalk between RNAs mediated by shared microRNAs (miRNAs) represents a novel layer of gene regulation, which plays important roles in development. In this study, we analyzed time series expression data for coding genes and long non-coding RNAs (lncRNAs) to identify thousands of interactions among competitive endogenous RNAs (ceRNAs) in four rhesus tissues. The ceRNAs exhibited dynamic expression and regulatory patterns during each tissue development process, which suggests that ceRNAs might work synergistically during different developmental stages or tissues to control specific functions. In addition, lncRNAs exhibit higher specificity as ceRNAs than coding-genes and their functions were predicted based on their competitive coding-gene partners to discover their important developmental roles. In addition to the specificity of tissue development, functional analyses demonstrated that the combined effects of multiple ceRNAs can have major impacts on general developmental and metabolic processes in multiple tissues, especially transcription-related functions where competitive interactions. Moreover, ceRNA interactions could sequentially and/or synergistically mediate the crosstalk among different signaling pathways during brain development. Analyzing ceRNA interactions during the development of multiple tissues will provideinsights in the regulation of normal development and the dysregulation of key mechanisms during pathogenesis.

Genome-wide differential expression of genes and small RNAs in testis of two different porcine breeds and at two different ages

Some documented evidences proved small RNAs (sRNA) and targeted genes are involved in mammalian testicular development and spermatogenesis. However, the detailed molecular regulation mechanisms of them remain largely unknown so far. In this study, we obtained a total of 10,716 mRNAs, 67 miRNAs and 16,953 piRNAs which were differentially expressed between LC and LW pig breeds or between the two sexual maturity stages. Of which, we identified 16 miRNAs and 28 targeted genes possibly related to spermatogenesis; 14 miRNA and 18 targeted genes probably associated with cell adhesion related testis development. We also annotated 579 piRNAs which could potentially regulate cell death, nucleosome organization and other basic biology process, which implied that those piRNAs might be involved in sexual maturation difference. The integrated network analysis results suggested that some differentially expressed genes were involved in spermatogenesis through the ECM-receptor interaction, focal adhesion, Wnt and PI3K-Akt signaling pathways, some particular miRNAs have the negative regulation roles and some special piRNAs have the positive and negative regulation roles in testicular development. Our data provide novel insights into the molecular expression and regulation similarities and diversities of spermatogenesis and testicular development in different pig breeds at different stages of sexual maturity.

Effects of DNA Methylation and Chromatin State on Rates of Molecular Evolution in Insects

Epigenetic information is widely appreciated for its role in gene regulation in eukaryotic organisms. However, epigenetic information can also influence genome evolution. Here, we investigate the effects of epigenetic information on gene sequence evolution in two disparate insects: the fly Drosophila melanogaster, which lacks substantial DNA methylation, and the ant Camponotus floridanus, which possesses a functional DNA methylation system. We found that DNA methylation was positively correlated with the synonymous substitution rate in C. floridanus, suggesting a key effect of DNA methylation on patterns of gene evolution. However, our data suggest the link between DNA methylation and elevated rates of synonymous substitution was explained, in large part, by the targeting of DNA methylation to genes with signatures of transcriptionally active chromatin, rather than the mutational effect of DNA methylation itself. This phenomenon may be explained by an elevated mutation rate for genes residing in transcriptionally active chromatin, or by increased structural constraints on genes in inactive chromatin. This result highlights the importance of chromatin structure as the primary epigenetic driver of genome evolution in insects. Overall, our study demonstrates how different epigenetic systems contribute to variation in the rates of coding sequence evolution.

DNA Methylation in Basal Metazoans: Insights from Ctenophores

Epigenetic modifications control gene expression without altering the primary DNA sequence. However, little is known about DNA methylation in invertebrates and its evolution. Here, we characterize two types of genomic DNA methylation in ctenophores, 5-methyl cytosine (5-mC) and the unconventional form of methylation 6-methyl adenine (6-mA). Using both bisulfite sequencing and an ELISA-based colorimetric assay, we experimentally confirmed the presence of 5-mC DNA methylation in ctenophores. In contrast to other invertebrates studied, Mnemiopsis leidyi has lower levels of genome-wide 5-mC methylation, but higher levels of 5-mC methylation in promoters when compared with gene bodies. Phylogenetic analysis showed that ctenophores have distinct forms of DNA methyltransferase 1 (DNMT1); the zf-CXXC domain type, which localized DNMT1 to CpG sites, and is a metazoan specific innovation. We also show that ctenophores encode the full repertoire of putative enzymes for 6-mA DNA methylation, and these genes are expressed in the aboral organ of Mnemiopsis. Using an ELISA-based colorimetric assay, we experimentally confirmed the presence of 6-mA methylation in the genomes of three different species of ctenophores, M. leidyi, Beroe abyssicola, and Pleurobrachia bachei. The functional role of this novel epigenomic mark is currently unknown. In summary, despite their compact genomes, there is a wide variety of epigenomic mechanisms employed by basal metazoans that provide novel insights into the evolutionary origins of biological novelties.

DNA Methylation within Transcribed Regions

DNA methylation within transcribed genes is commonly found in diverse animals and plants. Here, we provide an overview of recent advances and the remaining mystery regarding intragenic DNA methylation.

Determinants of the rate of protein sequence evolution

The rate and mechanism of protein sequence evolution have been central questions in evolutionary biology since the 1960s. Although the rate of protein sequence evolution depends primarily on the level of functional constraint, exactly what determines functional constraint has remained unclear. The increasing availability of genomic data has enabled much needed empirical examinations on the nature of functional constraint. These studies found that the evolutionary rate of a protein is predominantly influenced by its expression level rather than functional importance. A combination of theoretical and empirical analyses has identified multiple mechanisms behind these observations and demonstrated a prominent role in protein evolution of selection against errors in molecular and cellular processes.

Apparent microRNA-Target-specific Histone Modification in Mammalian Spermatogenesis

BACKGROUND

Epigenetics is an important mRNA expression regulator. However, how distinct epigenetic factors, such as microRNAs (miRNAs) and promoter methylation, cooperatively regulate mRNA expression is rarely discussed. Recently, apparent miRNA regulation of promoter methylation was identified by bioinformatic analysis; however, it has not yet been experimentally confirmed. If miRNA regulation of other epigenetic factors were identified, it would reveal another layer of epigenetic regulation. In this paper, histone modifications (H3K4me1, H3K4me3, H3K27me3, H3K27ac, H3K9ac, and H2AZ) during mammalian spermatogenesis were studied and the apparent miRNA-target-specific histone modification was investigated by bioinformatic analyses of publicly available datasets.

RESULTS

We identified several miRNAs’ target genes that are significantly associated with histone modification during mammalian spermatogenesis. MiRNAs that target genes associated with the most significant histone modifications are expressed before or during spermatogenesis; thus the results were convincing.

CONCLUSIONS

In this paper, we identified apparent miRNA regulation of histone modifications using a bioinformatics approach. The biological mechanisms of this effect should be further experimentally investigated.

CPuORF correlates with miRNA responsive elements on protein evolutionary rates

miRNA is increasingly being recognized as a key regulator of metabolism in animals. A wealth of evidence has suggested that miRNA mainly binds 3' UTR of mRNA and modulates the cell activities via repressing the mRNA translation. However, as the translation initiates at 5' UTR, cis elements like upstream open reading frame (uORF) resided in 5' UTR may also affect the translation efficiency or elongation. In this study, we performed a systematic analysis of miRNA responsive elements (MREs) and uORF of the same transcript in three model organisms (human, mouse, and Drosophila). Intriguingly, we found that the 3' UTR length grew with the complexity of species (human>mouse>Drosophila), in sharp contrast with the invariability of 5' UTR. Additionally, MRE number correlated well with the 3' UTR length, while uORF number showed a weak correlation with the 5' UTR length. Further, we found that human genes with conserved peptide upstream open reading frame (CPuORF) tend to have more MREs and lower evolutionary rates, which provides new insights into the correlation between UTR properties and translational control in animals.

The evolution of the coding exome of the Arabidopsis species--the influences of DNA methylation, relative exon position, and exon length

Background

The evolution of the coding exome is a major driving force of functional divergence both between species and between protein isoforms. Exons at different positions in the transcript or in different transcript isoforms may (1) mutate at different rates due to variations in DNA methylation level; and (2) serve distinct biological roles, and thus be differentially targeted by natural selection. Furthermore, intrinsic exonic features, such as exon length, may also affect the evolution of individual exons. Importantly, the evolutionary effects of these intrinsic/extrinsic features may differ significantly between animals and plants. Such inter-lineage differences, however, have not been systematically examined.

Results

Here we examine how DNA methylation at CpG dinucleotides (CpG methylation), in the context of intrinsic exonic features (exon length and relative exon position in the transcript), influences the evolution of coding exons of Arabidopsis thaliana. We observed fairly different evolutionary patterns in A. thaliana as compared with those reported for animals. Firstly, the mutagenic effect of CpG methylation is the strongest for internal exons and the weakest for first exons despite the stringent selective constraints on the former group. Secondly, the mutagenic effect of CpG methylation increases significantly with length in first exons but not in the other two exon groups. Thirdly, CpG methylation level is correlated with evolutionary rates (d_S, d_N, and the d_N/d_S ratio) with markedly different patterns among the three exon groups. The correlations are generally positive, negative, and mixed for first, last, and internal exons, respectively. Fourthly, exon length is a CpG methylation-independent indicator of evolutionary rates, particularly for d_N and the d_N/d_S ratio in last and internal exons. Finally, the evolutionary patterns of coding exons with regard to CpG methylation differ significantly between Arabidopsis species and mammals.

Conclusions

Our results suggest that intrinsic features, including relative exonic position in the transcript and exon length, play an important role in the evolution of A. thaliana coding exons. Furthermore, CpG methylation is correlated with exonic evolutionary rates differentially between A. thaliana and animals, and may have served different biological roles in the two lineages.