Email updates

Keep up to date with the latest news and content from BMC Systems Biology and BioMed Central.

This article is part of the supplement: Proceedings of the 23rd International Conference on Genome Informatics (GIW 2012)

Open Access Proceedings

Natural selection drives rapid evolution of mouse embryonic heart enhancers

Ben-Yang Liao* and Meng-Pin Weng

Author Affiliations

Division of Biostatistics & Bioinformatics, Institute of Population Health Sciences, National Health Research Institutes, Zhunan Town, Miaoli County 350, Taiwan, ROC

For all author emails, please log on.

BMC Systems Biology 2012, 6(Suppl 2):S1  doi:10.1186/1752-0509-6-S2-S1

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1752-0509/6/S2/S1


Published:12 December 2012

© 2012 Liao and Weng; licensee BioMed Central Ltd.

This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Mouse E11.5 embryonic heart enhancers were found to exhibit exceptionally weak sequence conservation during vertebrate evolution compared to enhancers of other developing organs. However, it is unknown whether this phenomenon is due to elevated mutation rates, or is a consequence of natural selection.

Results

In this study, based on the aligned orthologous genomic sequences of mouse and other closely related mammals, the substitution rates of fourfold degenerate sites or intron sequences in neighboring genes were used as neutral references to normalize substitution rates of mouse enhancers. Subsequent comparisons indicated that heart enhancers' evolutionary rates were increased by natural selection. Correspondingly, the results of Fisher's exact tests to examine the differential enrichment of substitutions between enhancers and neutral sequences suggest that both relaxed purifying selection and positive selection caused the rapid evolution of heart enhancers. Analyses on recombination rates and substitution patterns indicated that GC-biased gene conversion does not contribute to evolutionary rate variations among enhancers. In general, pleiotropic enhancers and enhancers in proximity to weakly expressed genes, tend to evolve slowly. Although heart enhancers are less pleiotropic and are adjacent to highly expressed genes, these biases do not account for the rapid evolution observed.

Conclusions

In combination, the results of the present study suggest that factors associated with functions or characteristics of the tissue may exert direct and profound effects on the intensity and direction of the natural selection applied to regulatory DNAs, such as enhancers.

Background

Evolution of gene expression may underlie the morphological diversity of animal species [1-3]. Therefore, it is important to understand the rules that govern the evolution of regulatory DNA sequences. Enhancers are a major category of regulatory DNAs, and are short genomic segments that bind to regulatory proteins and increase gene transcription [4,5]. Correspondingly, enhancers have been shown to regulate many of the processes that affect organismal development [6,7]. To locate these important regulatory non-coding DNAs in the mouse (Mus musculus) genome, several thousand mouse enhancers that function in the embryonic heart, forebrain, midbrain, or limbs were recently identified using a p300 ChIP-Seq approach [8-10]. Subsequent analyses showed that the evolutionary signatures of the enhancers associated with different tissue types were variable: Embryonic heart enhancers exhibited an exceptionally low sequence conservation compared to enhancers from other tissues [9]. Since the heart is a highly conserved organ critical to the survival of all vertebrate organisms, DNAs associated with the regulatory mechanisms that control heart development were expected to be conserved across a long phylogenetic distance [11]. Although the underlying cause for the rapid sequence evolution observed in embryonic heart enhancers remains unknown, this unexpected phenomenon has challenged the use of sequence conservation as a criterion for predicting functionally important non-coding DNAs [12,13].

There are three potential explanations for the fast evolution rate observed for heart enhancers. First, the mutation rate of heart enhancers are elevated, similar to that of fast evolving DNAs that reside in genomic regions prone to replication errors [14-18]. Second, heart enhancers are subject to a more frequent, or intense, positive selection, similar to that of genes encoding adaptive proteins where new mutations that develop tend to be fixed more frequently [19,20]. Third, heart enhancers have evolved under a more relaxed purifying selection, similar to that of genes with non-essential functions that have high evolution rates relative to essential genes [21,22]. To examine if differential frequencies/intensities in mutation rate or natural selection have led to an increased evolution of heart enhancers, based on mouse-rat or mouse-human genomic alignments, substitution rates of fourfold degenerate sites (d4), or intron sites (di), of neighboring genes were used as neutral references to normalize the substitution rate of each enhancer. Substitution rates of mouse enhancers from different tissues were then compared. Fisher's exact tests were also used to examine the differential enrichment of substitutions between enhancers and the corresponding neutral sequences, and these results were used to determine the strength of the purifying selection, or the frequency/intensity of the positive selection, experienced by the enhancers. We also examined the potential role for GC-biased gene conversion events in affecting the interpretation of processes that accelerate sequence substitutions of heart enhancers in mammals.

Studies on protein evolution have shown that unanticipated confounding factors could influence the identification and interpretation of previously reported protein evolutionary rate determinants. For instance, gene essentiality was asserted to influence bacterial protein evolutionary rates [23]. However, a later study found that this phenomenon was simply due to that highly expressed genes evolved slowly and essential genes were highly expressed [24]. Similarly, extracellularity was suggested to be correlated with the rate of protein evolution in yeasts [25]. Nevertheless, this was because extracellular proteins tend to be non-essential to the organism and non-essential genes were less constrained evolutionarily [26]. To determine whether the rapid evolution of heart enhancers is due to factors directly associated with functions or characteristics of the embryonic tissues or organs where a DNA performs functions (termed as "tissue factors" hereafter) [27,28], the influence of enhancer properties and neighboring genes to the evolution of these heart enhancers were examined. We found that the evolution rates of heart enhancers were correlated with pleiotropy (see Methods), as well as the expression levels of neighboring genes. Embryonic heart enhancers generally were less pleiotropic and were in proximity of highly expressed neighboring genes. After controlling for these biases, heart enhancers still exhibited highest rates of evolution. The results suggested that tissue factors may exert direct effects on the intensity and direction of the natural selection on enhancers. Further studies to elucidate how tissue factors specifically cause variations in selective forces on enhancers are needed.

Methods

Enhancers of developing mouse embryonic tissues

Based on the Ensembl mouse genome assembly v59 (NCBI m37), the co-activator protein, p300, was found to bind 2759, 2786, 3839, and 3597 enhancers in forebrain, midbrain, limb, and heart mouse embryonic tissues (Supplementary material in [9]). When enhancers with overlapping coordinates were consolidated into a single enhancer region, 11,332 non-overlapping enhancer regions were determined. These regions were further classified into "specific enhancers" (i.e., regions consolidated from enhancers of a single embryonic tissue) (n = 10,030) and "pleiotropic enhancers", the remaining set of enhancers (n = 1,302). It should be noted that, since only four tissues from the E11.5 stage were examined in p300 ChIP-seq experiments [9], the "specific enhancers" defined in the present study may function in other tissues or developmental stages during mouse embryogenesis. Furthermore, the term "pleiotropic" was used to describe enhancers that are likely to significantly affect the developmental processes of at least two of the examined embryonic organs in [9] when they are absent from the genome. The recombination rate of the mouse genome segment in which an enhancer was located was obtained from the Supplementary table of [29].

Properties of the mouse genes

The coordinates of the mouse genes analyzed were retrieved from the BioMart (http://www.biomart.org/ webcite). For mouse genes with a null phenotype, the essentiality of each gene was defined based on phenotypic annotations of Mouse Genome Informatics 4.21 (http://www.informatics.jax.org/ webcite), according to [30]. When premature death or infertility were the knockout phenotype, these genes were considered essential. All other genes with at least one documented null phenotype were considered non-essential. The gene expression profiles for the mouse genes analyzed were obtained from microarray data collected using 61 mouse tissues, including 53 adult tissues, 3 cell lines, and 5 early embryos from E6.5 to E10.5 [31]. Expression levels were calculated by averaging microarray-based expression signals from all tissues according to the previous studies [22,26].

Tests of natural selection on enhancer sequences

To study the patterns and determinants of enhancer evolution rates, aligned mouse-rat (Rattus norvegicus), mouse-human (Homo sapiens), and mouse-rat-human orthologous genomic regions were retrieved from Ensembl (http://www.ensembl.org/ webcite) [32], or the UCSC Genome Browser (http://genome.ucsc.edu webcite) [33] (genome versions: mouse: mm9; rat: rn4; human: hg18). Sequence divergence was subsequently calculated, and the numbers of substitutions were computed. Based on the genomic alignments retrieved from the UCSC Genome Browser (http://genome.ucsc.edu webcite) [33], the sequence divergence of enhancer regions (D) from mouse versus other mammal sequences was calculated according to baseml of PAML [34] using a GTR model with gamma distribution for site heterogeneity (model = 7; ncatG = 5). Fourfold degenerate sites and intron sequences of the nearest gene with an ortholog in the rat or human genome were used as neutral references to an enhancer. Since mammalian genes are often alternatively spliced and produce multiple transcripts [35], the coding sequence and intron sequence of a gene were defined using the exon-intron structure of the longest isoform annotated in Ensembl v59, with alignments of coding sequences and intron sequences retrieved from the UCSC Genome Browser. Fourfold degenerate sites of mouse genes were defined according to W-H Li [36]. When intron sequences were used as the neutral reference, similar to R Haygood, O Fedrigo, B Hanson, K-D Yokoyama and GA Wray [37], the first intron was removed due to the possibility it could contain regulatory motifs [38]. The latter criterion resulted in that intronless genes and single-intron genes have only the neutral reference of d4, and not di. Similarly, sequence divergences between mouse and another mammal under neutrality, d4 and di, were calculated by baseml of PAML [34] using a GTR model with gamma distribution for site heterogeneity (model = 7; ncatG = 5). As indicated above, the calculation of d4 or di required the presence of an Ensembl annotation of a mouse gene to its one-to-one ortholog to ensure the reliability of neutral estimation. The normalized enhancer sequence divergences (D/d4 or D/di) were then used as indexes for the relative strength of selection. Lower D/d4 or D/di values indicated a weaker purifying selection, or a stronger positive selection, had occurred for the enhancer sequence examined [1].

According to AP Rooney and J Zhang [39], Fisher's exact test was used to examine the natural selection on enhancers based on the 2x2 contingency table. The row and column categories of the table are "the enhancer vs. the neutral reference used" and "the number of substituted sites vs. the number of sites with no substitution", respectively. Both fourfold degenerate sites, or intron sequences of the neighboring gene (defined above), were used as neutral references to examine the differential enrichment of substitutions between enhancers and neutral references. Compared to the neutral reference, significantly more substitutions in an enhancer implied that positive selection had occurred. Alternatively, fewer substitutions in an enhancer suggested that purifying selection had occurred. Otherwise, the model of neutrality cannot be rejected statistically, and the enhancer would be considered to have evolved free from natural selection.

Results and discussion

Natural selection leads to rapid evolution of embryonic heart enhancers

In a previous study, an analysis of conservation depth using DNA alignments from a wide range of vertebrates identified the rapid evolution of embryonic heart enhancers [9]. In the present study, sequence divergence, D (see Methods), was used to examine the evolution rate of enhancers, while D/d4 or D/di (see Methods) ratios were used to examine the presence and direction of natural selection on non-coding sequences after the divergence of mouse-rat (Figure 1), or mouse-human (Additional file 1), genomes. Heart enhancers were found to have the highest D (indicating the weakest sequence conservation), while forebrain enhancers had the lowest D (indicating strongest sequence conservation) (mouse-rat: Figure 1A; mouse-human: Additional file 1). These results are consistent with the order of enhancer types determined by the previous analysis of conservation depth for these tissues (Figure 1b of [9]). If the differences in D determined for enhancers of forebrain, midbrain, limb and heart tissues is only due to variations in local mutation rates, these differences should disappear when D is normalized to neutral substitution rates (e.g., approximately by d4 or di). Based on the analysis performed, the orders of the enhancers of four tissue types in D/d4 (mouse-rat: Figure 1B; mouse-human: Additional file 1) or D/di (mouse-rat: Figure 1C; mouse-human: Additional file 1) followed the order shown in D, while such an order was neither observed for d4 (mouse-rat: Figure 1D; mouse-human: Additional file 1) nor di (mouse-rat: Figure 1E; mouse-human: Additional file 1). An estimation of mouse-rat D, D/d4, and D/di with the inclusion of human orthologous sequence as the outgroup in the model of baseml yielded virtually identical results (Additional file 2). Taken together, these results clearly indicate that natural selection, and not mutation, is responsible for the rapid sequence divergence (e.g., large D) associated with heart enhancers.

thumbnailFigure 1. Evolutionary rates of enhancers. Rates of evolution calculated for mouse enhancers of embryonic forebrain (FB), midbrain (MB), limb (LM), and heart (HT), included D (A), normalized evolutionary rates D/d4 (B) or D/di (C), and neutral substitution rates d4 (D) or di (E). The values of upper quartile, median, and lower quartile are indicated in each box, whereas the bars outside the box indicate semi-quartile ranges. D, D/d4, D/di, d4 and di were computed based on mouse-rat alignments. Pairwise comparisons showing significant differences in D, D/d4, D/di, d4, or di are connected with gray lines (Mann-Whitney U test).

Additional file 1. Supplementary figure S1. Evolutionary rates of enhancers. Rates of evolution calculated for mouse enhancers of embryonic forebrain (FB), midbrain (MB), limb (LM), and heart (HT), included D (A), normalized evolutionary rates D/d4 (B) or D/di (C), and neutral substitution rates d4 (D) or di (E). The values of upper quartile, median, and lower quartile are indicated in each box, whereas the bars outside the box indicate semi-quartile ranges. D, D/d4, D/di, d4 and di were computed based on mouse-human alignments. Pairwise comparisons showing significant differences in D, D/d4, D/di, d4, or di are connected with gray lines (Mann-Whitney U test).

Format: PDF Size: 248KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 2. Supplementary figure S2. Evolutionary rates of enhancers. Rates of evolution calculated for mouse enhancers of embryonic forebrain (FB), midbrain (MB), limb (LM), and heart (HT), included D (A), normalized evolutionary rates D/d4 (B) or D/di (C), and neutral substitution rates d4 (D) or di (E). The values of upper quartile, median, and lower quartile are indicated in each box, whereas the bars outside the box indicate semi-quartile ranges. D, D/d4, D/di, d4 and di were computed based on mouse-rat-human multiple alignments. Pairwise comparisons showing significant differences in D, D/d4, D/di, d4, or di are connected with gray lines (Mann-Whitney U test).

Format: PDF Size: 259KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

The higher D of heart enhancers may be due to a more relaxed purifying selection, or to a more frequent and stronger positive selection, that occurred during evolution. To distinguish between these possibilities, Fisher's exact test [39] was used on each enhancer to examine differences in the enrichment of substitutions between enhancers and the neutral reference [i.e., fourfold degenerate sites present in the neighboring gene (Table 1 and Additional file 3) or intron sites in the neighboring gene (sites of the first intron were excluded, see Methods) (Table 2 and Additional file 4)]. If the enhancer evolved without the influence of natural selection, no significant difference in the proportion of nucleotides that underwent substitution after mouse-rat or mouse-human divergence between the enhancer and the neutral reference would be observed. Otherwise, enrichment or depletion of substitutions in the enhancer compared to the neutral reference would suggest that the enhancer had undergone positive selection or purifying selection, respectively. To minimize the effect of enhancer pleiotropy (see Methods, and below), only "specific enhancers" (see Methods) with available neutral references were used in the analysis. The sample size of the enhancers tested for Table 2 (or Additional file 4) was smaller than that of Table 1 (or Additional file 4) since enhancers with a nearby gene that did not have introns or single-exon genes were excluded due to the lack of neutral intron sites. Analyses of both mouse-rat alignments (Table 1, 2) and mouse-human alignments (Additional file 3, 4) were performed. Several consistent observations were made. First, in comparison with enhancers of the forebrain, midbrain, and limb, a smaller proportion of heart enhancers were subject to natural selection (see item "Under selection/Total" in Table 1, 2 and Additional file 3, 4). These results indicate that heart enhancers have undergone a more relaxed purifying selection. Second, a larger proportion of heart enhancers had higher substitution rates than the neutral reference, compared to the enhancers of other tissues (see items "Positively selected/Under selection" and "Positively selected/Total" in Table 1, 2 and Additional file 3, 4). These results indicate that mouse embryonic heart enhancers experienced more frequent or more intense, positive selections during their evolution.

Table 1. The proportions of enhancers that underwent selection based on mouse-rat alignments (neutral reference: fourfold degenerate sites)

Additional file 3. Supplementary table S1. The proportions of enhancers that underwent selection based on mouse-human alignments (neutral reference: fourfold degenerate sites)

Format: PDF Size: 7KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Table 2. The proportions of enhancers that underwent selection based on mouse-rat alignments (neutral reference: intron sites)

Additional file 4. Supplementary table S2. The proportions of enhancers that underwent selection based on mouse-human alignments (neutral reference: intron sites)

Format: PDF Size: 7KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

It is important to note that estimates of the proportion of enhancers under selection were subject to the neutral references used. For example, when intron sites were served as the neutral reference, a higher proportion of enhancers were found to undergo natural selection (item "Under selection/Total"). In addition, "Positively selected/Under selection" and "Positively selected/Total" values were found to slightly differ. There are several possible explanations for these observations. First, the available numbers and samples of enhancers used to generate Table 1 (or Additional file 3) and Table 2 (or Additional file 4) varied. Second, an enhancer's neighboring gene generally had a larger number of intron sites than fourfold degenerate sites for Fisher's exact tests. With an increase in neutral reference sites, a statistically significant result to reject the neutral model was more likely to be obtained. Although our estimations of "Under selection/Total", "Positively selected/Under selection", and "Positively selected/Total" values may be affected by changes in the neutral reference used, these changes were small, especially for results based on mouse-human alignments (Additional file 3 and 4). In addition, it should be noted that the present study focused on the relative values for "Under selection/Total", "Positively selected/Under selection", and "Positively selected/Total", not the absolute values, to compare enhancers from different tissues. Estimations of these three indexes using a similar approach and neutral references should be subjected to the same biases. The consistent patterns observed for data in Tables 1 and 2 and Additional files 3 and 4 unambiguously indicate that a higher tendency of being targeted by relaxed purifying selection or positive selection was associated with the more rapid evolution of mouse embryonic heart enhancers.

GC-biased gene conversion does not explain the rapid evolution of heart enhancers

GC-biased gene conversion can mimic positive selection by generating similar genomic patterns, thereby confounding inferences regarding the type of selection that has occurred (reviewed in [40]). To investigate the possibility that a GC-biased gene conversion pattern could be misinterpreted as a positive selection of heart enhancers, we examined whether heart enhancers are located on genomic regions with higher recombination rates (recombination involves the formation of heteroduplex DNA that triggers gene conversion). Using the most recent recombination map of the mouse genome [29], although genomic regions containing heart enhancers were found to have relatively high recombination rates, these recombination rates were not higher than those of the genomic regions containing limb enhancers (Figure 2A) (the D, D/d4, and D/di for heart enhancers were higher than those of limb enhancers, see Figure 1, A-C).

thumbnailFigure 2. Compared to other enhancers, heart enhancers are not subject to more frequent GC-biased gene conversions. Heart enhancers (A) are not preferentially localized to genomic regions with higher recombination rates and (B) do not contain more A/T to G/C substitutions in the mouse lineages. The values of upper quartile, median, and lower quartile are indicated in each box, whereas the bars outside the box indicate semi-quartile ranges. Pairwise comparisons showing significant differences are connected with gray lines (*0.01<P≦0.05 and ***P≦0.001 by Mann-Whitney U test). FB, forebrain; MB, midbrain; LB, limb; HT, heart.

Furthermore, based on alignments of mouse-rat-human orthologous DNAs, the number of nucleotides that were ancestrally A or T before mouse-rat divergence (human and rat are both either A or T of the aligned site) and had a substitution event in the mouse lineage was calculated and defined as "Nsubstituted A/T" for each enhancer. The number of mouse lineage-specific substitutions from A/T to G/C (NA/T→G/C) was also calculated for each enhancer. If heart enhancers underwent more GC-biased gene conversion events, they would be predicted to have a higher NA/T→G/C to Nsubstituted A/T ratio. However, heart enhancers were associated with a significantly lower ratio of NA/T→G/C to Nsubstituted AT (Figure 2B) compared to the enhancers of other tissues.

In combination, the results shown in Figure 2 indicate that heart enhancers did not experience more GC-biased gene conversion events over their evolution. Therefore, mechanisms associated with GC-biased gene conversion did not explain the rapid evolution of heart enhancers and the selective pressure acting on them.

Determinants for the evolution rate of enhancer sequences

To understand the source of the selective forces acting on heart enhancers, the factors potentially associated with the evolution rate of enhancer sequences, and the properties of heart enhancers, were investigated. Two types of enhancer characteristics were examined: enhancer pleiotropy and the properties of neighboring genes. As mentioned in Methods, enhancer regions that only functioned in a single tissue were defined as specific enhancer regions, while all others were defined as pleiotropic enhancer regions. Initially, specific enhancer regions were found to have a significantly higher D (P < e-300, Mann-Whitney U test, Figure 3A) compared to pleiotropic enhancer regions. These results indicated that specific enhancer regions had a faster rate of sequence evolution. Furthermore, higher D/d4 (P = 6.7e-16, U test, Figure 3D) ratios of specific enhancer regions suggests that the faster evolution rate observed is due to a fixation bias, rather than a mutation bias. To avoid redundancy, analyses based on D/di were not conducted. Previous studies reported that pleiotropic [41], or housekeeping protein coding genes [42], are subject to more stringent selective constraints. Similarly, the evolution of pleiotropic enhancer regions could be constrained by having a role in the development of a wide range of biological systems. As a result, mutations in these enhancer regions would have greater deleterious effects. These hypotheses are consistent with the patterns observed in Figure 3A and 3D.

thumbnailFigure 3. Correlations between the evolutionary rates of mouse enhancers and enhancer properties: the tissue-specificity of the enhancers (A, D), essentiality of the adjacent gene (B, E), and the expression level of the adjacent gene (C, F). Evolutionary rates are represented by D (A-C) or D/d4 (D-F). Spearman's rank correlation coefficient is (C) 0.040 (P <10-3) and (F) 0.046 (P <10-4) for the unbinned data. The values of upper quartile, median, and lower quartile are indicated in each box, whereas the bars outside the box indicate semi-quartile ranges. Pairwise comparisons showing significant differences in D or D/d4 are connected with gray lines (*0.01<P≦0.05 and ***P≦0.001 by Mann-Whitney U test).

Enhancers tend to regulate neighboring genes that are within several kilobases [9]. Therefore, we also investigated the association between the evolution rate of enhancer regions and the properties of adjacent genes (e.g., the closest coding gene in the mouse genome). The median distance for adjacent genes was found to be 49.7 kb. It is important to note that the adjacent gene of an enhancer may not be the same gene used to compute d4 or di, since the calculation of d4 and di require the presence of a one-to-one ortholog in another genome. For enhancer regions with an adjacent gene that is considered to be essential, a significantly smaller D (P = 3.83e-6, U test) was observed compared to enhancer regions adjacent to a non-essential gene (Figure 3B). These results suggest that faster rates of evolution occurred in the latter case. It has previously been shown that genomic regions with high gene density have low mutation rates [43]. Moreover, highly expressed genes, which tend to be essential [24,44], are located in gene-dense genomic regions [45]. Thus, we expected that genomic regions enriched with essential genes would have a lowered mutation rate which potentially explains Figure 3B. As expected, a smaller D for enhancer regions with an adjacent gene that is essential was observed, while no significant difference in D/d4 was observed for enhancers with an adjacent essential gene versus a non-essential gene (P = 0.27, U test, Figure 3E).

In addition, enhancer regions with a highly expressed adjacent gene tended to evolve faster than enhancer regions adjacent to a weakly expressed gene (Spearman's rank correlation coefficient σ of D vs. the expression level of the adjacent gene = 0.040, P = 1.9e-4; Figure 3C, for U test). A slightly more significant positive correlation was found between D/d4 and the expression level of an adjacent gene (Spearman's σ = 0.044, P = 9.8e-5; Figure 3F, for U test), suggesting that embryonic enhancers near weakly expressed genes are subject to a more biased selective force. The results remained unchanged when gene expression levels were defined by the maximum expression signal detected for the 61 mouse tissues analyzed (Additional file 5). However, at the present stage, it remains unclear whether expression levels of an adjacent gene directly, or indirectly, affect D/d4. Determining the cause of such an effect needs future studies.

Additional file 5. Supplementary figure S3. Enhancers with a highly expressed adjacent gene evolve more rapidly. Evolutionary rates are represented by D (A) or D/d4 (B). Gene expression level was defined by the maximum expression signal of the 61 mouse tissues. Spearman's rank correlation coefficient is (C) 0.044 (P <10-4) and (F) 0.042 (P <10-3) for the unbinned data. The values of upper quartile, median, and lower quartile are indicated in each box, whereas the bars outside the box indicate semi-quartile ranges. Pairwise comparisons showing significant differences in D or D/d4 are connected with gray lines (*0.01<P≦0.05 and ***P≦0.001 by Mann-Whitney U test).

Format: PDF Size: 211KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Local genomic and regulatory characteristics do not explain the rapid evolution of heart enhancers

As described in the Introduction, studies of protein evolution have shown that unanticipated confounding factors can influence the identification and the interpretation of determinants for coding DNA evolutionary rates. Therefore, variations in the evolution rates of different embryonic tissues can result from factors indirectly related to the tissue factor. In the present study, enhancers from the forebrain, midbrain, limb, and heart did not only differ in the tissue types from which they were identified, but also in several of their general properties that influenced D as well as D/d4. For embryonic heart enhancers, they typically were not pleiotropic (Figure 4A), and they tended to be physically close to highly expressed genes (Figure 4B). Since both tissue-specific enhancers and enhancers near highly expressed genes have been associated with fast evolution rates due to natural selection (Figure 3D and 3F), these biases could also contribute to the rapid evolution of embryonic heart enhancers.

thumbnailFigure 4. Properties of heart enhancers. Heart enhancers include more specific enhancers (A) and enhancers adjacent to highly expressed genes (B). (A) The numbers in boxes represent the number of enhancers in a group. (B) The values of upper quartile, median, and lower quartile are indicated in each box, whereas the bars outside the box indicate semi-quartile ranges. Pairwise comparisons with significantly different expression levels of the adjacent gene are connected with gray lines (*0.01<P≦0.05 and ***P≦0.001 by (A) × 2 test or (B) Mann-Whitney U test). FB, forebrain; MB, midbrain; LB, limb; HT, heart.

To determine the role of enhancer pleiotropy and proximity to highly expressed genes, D (Figure 4A) and D/d4 (Figure 4B) values were compared for enhancers of the embryonic forebrain, midbrain, limbs, and heart, after controlling for the aforementioned two biases. To control for tissue specificity, only specific enhancers were examined. To control for the expression level of adjacent genes, enhancers were divided into three groups depending on the expression level of the adjacent gene, < 200, 200-400, and ≥ 400. D and D/d4 were then compared for the enhancers of the four tissues within a given group. After such controls, embryonic heart enhancers still had the highest D (Figure 5A) and D/d4 (Figure 5B) values. These results suggest that the relatively low pleiotropy of embryonic heart enhancers, and their physical proximity to highly expressed genes, only partially accounts for their high rate of evolution and biases in selection.

thumbnailFigure 5. Pleiotropy and expression levels of adjacent genes do not explain enhancer evolutionary rate variations. Comparisons of D (A) and D/d4 (B) for mouse forebrain (FB), midbrain (MB), limb (LM), and heart (HT) enhancers after controlling for the enhancer pleiotropy and the expression level of adjacent genes. Only specific enhancers are included in the analysis. The values of upper quartile, median, and lower quartile are indicated in each box, whereas the bars outside the box indicate semi-quartile ranges. Pairwise comparisons showing significant difference in D or D/d4 are connected with gray lines (***P≦0.001 by Mann-Whitney U test).

Conclusions

The rapid evolution rates of mouse embryonic heart enhancers can potentially be explained by elevated rates of mutation and/or changes in the direction or intensity of natural selection. In the present study, neutral sites of the neighboring gene were used as references to examine mutation bias, as well as selective bias, on the evolution of enhancer sequences. Subsequent analyses demonstrated that the rapid evolution of mouse E11.5 heart enhancers cannot be explained by potential confounding factors examined in the present study, and therefore, may be directly associated with embryonic tissue factors.

The identification and understanding of regulatory DNAs in the mammalian genome is a major goal in the post-genomic era. This study employed a simple, yet effective, method to explore the role of natural selection on the evolution of non-coding sequences. Although the present study indicates that the rapid evolution of heart enhancers is likely to have been the result of a selective force associated with the type of embryonic tissue involved, the underlying causes for the relaxed selective constraint, or the higher frequency/intensity of positive selection, remain to be discovered. One possible underlying cause for a more relaxed purifying selection on heart enhancer is that during mouse embryogenesis, heart development begins and finishes earlier than brain and limb developments. If all enhancer sequences identified in [9] are functional regulatory DNA sequences at E11.5, these enhancers are regulating late-stage heart organogenesis but early-stage organogenesis of the brain or limbs. Mutations resulting in abnormal organ development at a later stage of organogenesis could bring a milder defect than those resulting in defects occur at an earlier stage, and therefore are more tolerable to the organism.

Our finding that heart enhancers tend to be positively selected is particularly intriguingly. In the future, it would be interesting to determine whether heart enhancers have been repeatedly shaped by positive selection because they are also involved in other developmental processes associated with morphological or physiological diversity. With the increasing abundance of functional genomic data and phenotypic data for mouse genes, such investigations will be feasible when a more complete set of spatial and temporal p300 ChIP-Seq data for mouse embryos and other species becomes available.

List of abbreviations

d4: substitution rate of fourfold degenerate sites; di: substitution rate of intron sites; D: sequence divergence of enhancer regions; FB: forebrain; MB: midbrain; LB: limb; HT: heart.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

B.-Y.L. designed research; B.-Y.L. and M.-P.W. performed research; M.-P.W. analyzed data; B.-Y.L. wrote the paper. All authors read and approved the final manuscript.

Acknowledgements and funding

This work was supported by intramural funding from the National Health Research Institutes, Taiwan, and research grant (NSC 101-2311-B-400-001-MY3) and (NSC 99-2311-B-400-003-MY2) from the National Science Council, Taiwan, to B.-Y.L.

This article has been published as part of BMC Systems Biology Volume 6 Supplement 2, 2012: Proceedings of the 23rd International Conference on Genome Informatics (GIW 2012). The full contents of the supplement are available online at http://www.biomedcentral.com/bmcsystbiol/supplements/6/S2.

References

  1. Liao B-Y, Weng M-P, Zhang J: Contrasting genetic paths to morphological and physiological evolution.

    Proc Natl Acad Sci USA 2010, 107(16):7353-7358. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  2. Wray GA: The evolutionary significance of cis-regulatory mutations.

    Nat Rev Genet 2007, 8(3):206-216. PubMed Abstract | Publisher Full Text OpenURL

  3. Carroll SB: Evo-devo and an expanding evolutionary synthesis: a genetic theory of morphological evolution.

    Cell 2008, 134(1):25-36. PubMed Abstract | Publisher Full Text OpenURL

  4. Ptashne M: Gene regulation by proteins acting nearby and at a distance.

    Nature 1986, 322(6081):697-701. PubMed Abstract | Publisher Full Text OpenURL

  5. Dillon N, Sabbattini P: Functional gene expression domains: defining the functional unit of eukaryotic gene regulation.

    BioEssays 2000, 22(7):657-665. PubMed Abstract | Publisher Full Text OpenURL

  6. Pennacchio LA, Ahituv N, Moses AM, Prabhakar S, Nobrega MA, Shoukry M, Minovitsky S, Dubchak I, Holt A, Lewis KD, et al.: In vivo enhancer analysis of human conserved non-coding sequences.

    Nature 2006, 444(7118):499-502. PubMed Abstract | Publisher Full Text OpenURL

  7. Borok MJ, Tran DA, Ho MCW, Drewell RA: Dissecting the regulatory switches of development: lessons from enhancer evolution in Drosophila.

    Development 2010, 137(1):5-13. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Robertson G, Hirst M, Bainbridge M, Bilenky M, Zhao Y, Zeng T, Euskirchen G, Bernier B, Varhol R, Delaney A, et al.: Genome-wide profiles of STAT1 DNA association using chromatin immunoprecipitation and massively parallel sequencing.

    Nat Methods 2007, 4(8):651-657. PubMed Abstract | Publisher Full Text OpenURL

  9. Blow MJ, McCulley DJ, Li Z, Zhang T, Akiyama JA, Holt A, Plajzer-Frick I, Shoukry M, Wright C, Chen F, et al.: ChIP-Seq identification of weakly conserved heart enhancers.

    Nat Genet 2010, 42(9):806-810. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Visel A, Blow MJ, Li Z, Zhang T, Akiyama JA, Holt A, Plajzer-Frick I, Shoukry M, Wright C, Chen F, et al.: ChIP-seq accurately predicts tissue-specific activity of enhancers.

    Nature 2009, 457(7231):854-858. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  11. Hardison RC: Variable evolutionary signatures at the heart of enhancers.

    Nat Genet 2010, 42(9):734-735. PubMed Abstract | Publisher Full Text OpenURL

  12. Frazer KA, Sheehan JB, Stokowski RP, Chen X, Hosseini R, Cheng JF, Fodor SP, Cox DR, Patil N: Evolutionarily conserved sequences on human chromosome 21.

    Genome Res 2001, 11(10):1651-1659. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Siepel A, Bejerano G, Pedersen JS, Hinrichs AS, Hou M, Rosenbloom K, Clawson H, Spieth J, Hillier LW, Richards S, et al.: Evolutionarily conserved elements in vertebrate, insect, worm, and yeast genomes.

    Genome Res 2005, 15(8):1034-1050. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Jeffreys AJ, Royle NJ, Wilson V, Wong Z: Spontaneous mutation rates to new length alleles at tandem-repetitive hypervariable loci in human DNA.

    Nature 1988, 332(6161):278-281. PubMed Abstract | Publisher Full Text OpenURL

  15. Duret L, Arndt PF: The impact of recombination on nucleotide substitutions in the human genome.

    PLoS Genet 2008, 4(5):e1000071. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Hodgkinson A, Ladoukakis E, Eyre-Walker A: Cryptic variation in the human mutation rate.

    PLoS Biol 2009, 7(2):e1000027. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Tian D, Wang Q, Zhang P, Araki H, Yang S, Kreitman M, Nagylaki T, Hudson R, Bergelson J, Chen J-Q: Single-nucleotide mutation rate increases close to insertions/deletions in eukaryotes.

    Nature 2008, 455(7209):105-108. PubMed Abstract | Publisher Full Text OpenURL

  18. Bird AP: DNA methylation and the frequency of CpG in animal DNA.

    Nucleic Acids Res 1980, 8(7):1499-1504. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Vallender EJ, Lahn BT: Positive selection on the human genome.

    Hum Mol Genet 2004, 13(Spec No 2):R245-254. PubMed Abstract | Publisher Full Text OpenURL

  20. Zhang J: Evolution of the human ASPM gene, a major determinant of brain size.

    Genetics 2003, 165(4):2063-2070. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  21. Zhang J, He X: Significant impact of protein dispensability on the instantaneous rate of protein evolution.

    Mol Biol Evol 2005, 22(4):1147-1155. PubMed Abstract | Publisher Full Text OpenURL

  22. Liao B-Y, Scott NM, Zhang J: Impacts of gene essentiality, expression pattern, and gene compactness on the evolutionary rate of mammalian proteins.

    Mol Biol Evol 2006, 23(11):2072-2080. PubMed Abstract | Publisher Full Text OpenURL

  23. Jordan IK, Rogozin IB, Wolf YI, Koonin EV: Essential genes are more evolutionarily conserved than are nonessential genes in bacteria.

    Genome Res 2002, 12(6):962-968. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Rocha EP, Danchin A: An analysis of determinants of amino acids substitution rates in bacterial proteins.

    Mol Biol Evol 2004, 21(1):108-116. PubMed Abstract | Publisher Full Text OpenURL

  25. Julenius K, Pedersen AG: Protein evolution is faster outside the cell.

    Mol Biol Evol 2006, 23(11):2039-2048. PubMed Abstract | Publisher Full Text OpenURL

  26. Liao B-Y, Weng M-P, Zhang J: Impact of extracellularity on the evolutionary rate of mammalian proteins.

    Genome Biol Evol 2010, 2010:39-43. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Gu X, Su Z: Tissue-driven hypothesis of genomic evolution and sequence-expression correlations.

    Proc Natl Acad Sci USA 2007, 104(8):2779-2784. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Tuller T, Kupiec M, Ruppin E: Evolutionary rate and gene expression across different brain regions.

    Genome Biol 2008, 9(9):R142. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  29. Jensen-Seaman MI, Furey TS, Payseur BA, Lu Y, Roskin KM, Chen CF, Thomas MA, Haussler D, Jacob HJ: Comparative recombination rates in the rat, mouse, and human genomes.

    Genome Res 2004, 14(4):528-538. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Liao B-Y, Zhang J: Mouse duplicate genes are as essential as singletons.

    Trends Genet 2007, 23(8):378-381. PubMed Abstract | Publisher Full Text OpenURL

  31. Su AI, Wiltshire T, Batalov S, Lapp H, Ching KA, Block D, Zhang J, Soden R, Hayakawa M, Kreiman G, et al.: A gene atlas of the mouse and human protein-encoding transcriptomes.

    Proc Natl Acad Sci USA 2004, 101(16):6062-6067. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Hubbard TJ, Aken BL, Beal K, Ballester B, Caccamo M, Chen Y, Clarke L, Coates G, Cunningham F, Cutts T, et al.: Ensembl 2007.

    Nucleic Acids Res 2007, 35(Database issue):D610-617. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  33. Jukes TH, Cantor CR: Evolution of Protein Molecules. New York: Academic Press; 1969.

  34. Yang Z: PAML 4: phylogenetic analysis by maximum likelihood.

    Mol Biol Evol 2007, 24(8):1586-1591. PubMed Abstract | Publisher Full Text OpenURL

  35. Sultan M, Schulz MH, Richard H, Magen A, Klingenhoff A, Scherf M, Seifert M, Borodina T, Soldatov A, Parkhomchuk D, et al.: A global view of gene activity and alternative splicing by deep sequencing of the human transcriptome.

    Science 2008, 321(5891):956-960. PubMed Abstract | Publisher Full Text OpenURL

  36. Li W-H: Molecular Evolution. Sunderland, MA: Sinauer Associates; 1997.

  37. Haygood R, Fedrigo O, Hanson B, Yokoyama K-D, Wray GA: Promoter regions of many neural- and nutrition-related genes have experienced positive selection during human evolution.

    Nat Genet 2007, 39(9):1140-1144. PubMed Abstract | Publisher Full Text OpenURL

  38. Blanchette M, Bataille AR, Chen X, Poitras C, Laganiere J, Lefebvre C, Deblois G, Giguere V, Ferretti V, Bergeron D, et al.: Genome-wide computational prediction of transcriptional regulatory modules reveals new insights into human gene expression.

    Genome Res 2006, 16(5):656-668. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  39. Rooney AP, Zhang J: Rapid evolution of a primate sperm protein: relaxation of functional constraint or positive Darwinian selection?

    Mol Biol Evol 1999, 16(5):706-710. PubMed Abstract | Publisher Full Text OpenURL

  40. Duret L, Galtier N: Biased gene conversion and the evolution of mammalian genomic landscapes.

    Annu Rev Genomics Hum Genet 2009, 10:285-311. PubMed Abstract | Publisher Full Text OpenURL

  41. He X, Zhang J: Toward a molecular understanding of pleiotropy.

    Genetics 2006, 173(4):1885-1891. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  42. Zhang L, Li WH: Mammalian housekeeping genes evolve more slowly than tissue-specific genes.

    Mol Biol Evol 2004, 21(2):236-239. PubMed Abstract | Publisher Full Text OpenURL

  43. Payseur BA, Nachman MW: Gene density and human nucleotide polymorphism.

    Mol Biol Evol 2002, 19(3):336-340. PubMed Abstract | Publisher Full Text OpenURL

  44. Xia Y, Franzosa EA, Gerstein MB: Integrated assessment of genomic correlates of protein evolutionary rate.

    PLoS Comput Biol 2009, 5(6):e1000413. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  45. Versteeg R, van Schaik BD, van Batenburg MF, Roos M, Monajemi R, Caron H, Bussemaker HJ, van Kampen AH: The human transcriptome map reveals extremes in gene density, intron length, GC content, and repeat pattern for domains of highly and weakly expressed genes.

    Genome Res 2003, 13(9):1998-2004. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL