Skip to main content

Homoplastic microinversions and the avian tree of life

Abstract

Background

Microinversions are cytologically undetectable inversions of DNA sequences that accumulate slowly in genomes. Like many other rare genomic changes (RGCs), microinversions are thought to be virtually homoplasy-free evolutionary characters, suggesting that they may be very useful for difficult phylogenetic problems such as the avian tree of life. However, few detailed surveys of these genomic rearrangements have been conducted, making it difficult to assess this hypothesis or understand the impact of microinversions upon genome evolution.

Results

We surveyed non-coding sequence data from a recent avian phylogenetic study and found substantially more microinversions than expected based upon prior information about vertebrate inversion rates, although this is likely due to underestimation of these rates in previous studies. Most microinversions were lineage-specific or united well-accepted groups. However, some homoplastic microinversions were evident among the informative characters. Hemiplasy, which reflects differences between gene trees and the species tree, did not explain the observed homoplasy. Two specific loci were microinversion hotspots, with high numbers of inversions that included both the homoplastic as well as some overlapping microinversions. Neither stem-loop structures nor detectable sequence motifs were associated with microinversions in the hotspots.

Conclusions

Microinversions can provide valuable phylogenetic information, although power analysis indicates that large amounts of sequence data will be necessary to identify enough inversions (and similar RGCs) to resolve short branches in the tree of life. Moreover, microinversions are not perfect characters and should be interpreted with caution, just as with any other character type. Independent of their use for phylogenetic analyses, microinversions are important because they have the potential to complicate alignment of non-coding sequences. Despite their low rate of accumulation, they have clearly contributed to genome evolution, suggesting that active identification of microinversions will prove useful in future phylogenomic studies.

Background

Reconstructing the evolutionary relationships among organisms and changes in their genomes are major goals of phylogenomics [1–3]. The characteristics of genomes that have been used to reconstruct evolutionary history reflect the multitude of changes that arise due to distinct mutational mechanisms and accumulate at a variety of rates (Figure 1). The most slowly accumulating changes, collectively designated rare genomic changes (RGCs), reflect a heterogeneous set of mutational processes. RGCs include transposable element insertions (e.g., Kriegs et al. [4]), gene order changes [5], and additional less-studied phenomena [6–8]. Microinversions [6] are one of these relatively poorly-studied types of RGCs.

Figure 1
figure 1

Approximate rates of accumulation for different genomic changes over evolutionary time. Details of the literature survey used to estimate these rates are provided in Additional file 2. The estimate of the avian microinversion rate reflects the results of this paper. Estimates of evolutionary rates for nucleotide substitutions and indels in birds appear lower than those for mammals, consistent with some previous publications [59], but it is important to note that substantial rate variation occurs within each group (e.g., [27, 60]). As described in the text, it may be better to interpret prior estimates of the mammalian microinversion rate as the rate at which relatively long microinversions accumulate.

Despite this heterogeneity, RGCs are thought to exhibit less homoplasy (evolutionary convergence and reversals) than nucleotide substitutions [9]. Indeed, some RGCs have been viewed as "perfect" homoplasy-free (or virtually homoplasy-free) characters. Establishing that specific types of RGCs, like microinversions, are perfect characters is important for two reasons. First, it would provide information about the mutational and evolutionary processes that underlie their accumulation, illuminating processes that contribute to genome evolution. Second, perfect RGCs could provide a practical means to assemble the tree of life because phylogenetic reconstruction is straightforward when homoplasy is absent [6].

Even perfect RGCs can appear homoplastic when found in genomic regions with an evolutionary history incongruent with the species tree [5, 10]. The appearance of homoplasy due to incomplete lineage sorting, called hemiplasy [11], typically occurs in trees with short internal branches [12, 13]. However, rapid radiations with short internal branches ("bushes" or "biological big bangs") may be relatively common events in the tree of life [14, 15]. This suggests that analyses of RGC data should consider hemiplasy explicitly.

Microinversions are defined as cytologically undetectable inversions [6], although in practice the size range considered depends on the type of data examined and method used for detection. Feuk et al. [16] classified inversions ranging in size from 23 base pairs (bp) to 62 megabases (Mb) as microinversions, whereas Ma et al. [1] considered all inversions greater than 50 kilobases (kb) to be "large" inversions rather than microinversions. The lower limit also varies, going down to 4 bp [17]. Not surprisingly, studies using whole genomes (e.g., [1, 16]) have identified larger inversions, while phylogenetic studies (often restricted to a single locus or region of an organellar genome) have typically revealed much smaller microinversions (e.g., [17–21]). Nonetheless, the size spectra reported for genome-scale and phylogenetic studies overlap, suggesting that both types of studies include at least some inversions that result from similar biological phenomena. Using the term "microinversion" to refer to inversions that are long enough to include one or more complete genes seems inappropriate, suggesting that it should be reserved for shorter inversions. However, this criterion may be difficult to apply in practice, since the length of genes exhibits substantial variation among organisms and within genomes. The majority of genes are <50 kb in length in most vertebrate lineages, suggesting that the Ma et al. [1] size criterion may be appropriate and simple to use. Therefore, we recommend using 50 kb as the maximum size for microinversions in most vertebrate genomes, although we also note that the most appropriate size criterion is likely to depend upon the focal organism.

The hypothesis that microinversions and other RGCs are perfect characters reflects both their large state space (number of potential character states) and slow rate of accumulation over evolutionary time, making independent changes to the same state unlikely. The state space for different RGCs will depend upon the details of each type of genomic change, but it seems likely that the state space for microinversions is large; they can be of a variety of lengths and have any specific nucleotide for endpoints, making it unlikely that independent microinversions will appear identical. Previous studies have also suggested that microinversions accumulate at a very low rate (Figure 1), although this observation may be biased by the size spectrum of the inversions that were identified and considered to be microinversions. Ma et al. [1] reported that smaller microinversions (they identified inversions as short as 31 bp) occur more frequently than larger ones. However, the rate of accumulation for inversions that are even shorter than those identified by Ma et al. [1] remains unclear and these differences among previous studies make direct comparisons challenging. Nonetheless, it seems certain that microinversions accumulate at least several orders of magnitude more slowly than nucleotide substitutions. Thus, the hypothesis that microinversions are perfect characters that will be very useful for assembling the tree of life remains reasonable.

The mechanism(s) responsible for microinversion accumulation remain poorly characterized, making empirical tests of the "perfect character hypothesis" for these relatively poorly studied RGCs critical. Indeed, homoplastic microinversions have been identified in angiosperm chloroplast genomes [17, 19], in contrast to expectation based upon the perfect character hypothesis. Most chloroplast microinversions appear to be associated with palindromic sequences that have the potential to form stem-loop structures in transcripts [17, 19] and these palindromes may facilitate inversion. Indeed, Catalano et al. [21] reported that microinversions are correlated with higher stability of the hairpins that have the potential to form at these stem-loop regions, in agreement with the hypothesis that hairpin formation facilitates inversion. Since many chloroplast stem-loop structures have regulatory functions (e.g., Stern et al. [22]) they are typically conserved, creating the potential for recurrent inversions at specific sites. Regulatory stem-loops are present in vertebrate introns (e.g., Hugo et al. [23]) and at least one vertebrate microinversion noted in a vertebrate phylogenetic study was associated with an inverted repeat [18]. However, conserved stem-loops appear to be uncommon in vertebrate introns whereas chloroplast stem-loops are relatively common [22, 24]. This difference is consistent with the observation that few animal microinversions appear homoplastic [6, 25]. Indeed, all microinversions observed in those studies were either homoplasy-free or conflicted with short branches. Thus, the small number of animal microinversions that appear to conflict with the species tree based upon other data may result from hemiplasy rather than homoplasy. Thus, microinversions in animal nuclear genomes remain candidates for "ideal RGCs", able to recover branches in gene trees accurately.

Microinversions can be difficult to identify, making the study of these interesting and phylogenetically useful genomic changes challenging. In fact, ~80% of the inversions identified in the Feuk et al. [16] comparison of the human and chimpanzee genomes were later suggested to be contig assembly artifacts [6]. This problem can be solved by restricting the term microinversion to the shortest part of the inversion spectrum, limiting the maximum size of the microinversions to less than the length of an individual sequencing read (i.e., focusing on inversions that are <400 bp for Sanger sequencing). Comparing closely related taxa also has the potential to facilitate microinversion identification. Indeed, most microinversions identified in a comparison of four mammalian genomes were found in the two most closely related taxa [1]. Here we use these strategies to identify microinversions in non-coding regions associated with 17 loci from 169 birds. We examined variation among loci in the microinversion rate (hereafter abbreviated λMI), identified phylogenetically informative and homoplastic microinversions, and found evidence that the number of microinversions has been underestimated in previous large-scale studies.

Methods

Sequencing, Alignment and Microinversion Identification

We primarily used published data [26–28], although some novel CLTCL1 sequences were generated using the primers and PCR conditions from Kimball et al. [29] (for details, see Additional file 1). For this study, we focused on shorter sequences with extensive taxon sampling (Table 1) instead of complete genomic sequences [26–28]. Sequences were aligned manually, sometimes starting from an alignment produced in an automated manner (i.e., using Clustal [30] or MAFFT [31]). Alignments were refined iteratively with input from at least two different individuals. During this process alignments were examined carefully; this resulted in the identification of a number of microinversions "by eye" (Additional file 2, Table S2).

Table 1 Estimates of the microinversion rate (λMI) for different loci.

Microinversions were also identified by a computational method that combined the multiple sequence alignments with the results of complementary strand alignments for all pairs of sequences (Additional file 2, Figure S1). The pairwise complementary strand alignments were generated using bl2seq [32] and YASS [33] and mapped onto the multiple sequence alignments using a program written by ELB. This program saved a table that included the first and last positions of each pairwise complementary strand alignment in the multiple sequence alignment and highlighted the overlapping pairwise complementary strand alignments (an example is presented in Additional File 3 along with a description of the algorithm in pseudocode). Microinversions are expected to result in complementary strand alignments that either overlap or are located near each other in the sequence alignment. The presence or absence of microinversions at each position identified as a significant complementary strand hit involving sequences that were overlapping or located near each other in the multiple sequence alignment was then validated by visual inspection. Microinversion endpoints were assigned based upon the length of the complementary strand alignments, although there were some cases where inversion endpoints were difficult to identify (e.g., Figure 2). Validating microinversions shorter than 5 bp was difficult, so that was the minimum size considered.

Figure 2
figure 2

Example of a microinversion. (a) A conserved region in TPM1 intron 6 with a 24 bp microinversion (outlined in white) in Trogon personatus. (b) Inverting the Trogon sequence (indicated in lower-case) results in a sequence identical to Pharomachrus auriceps, its sister taxon in the tree.

The DNA mfold server (http://mfold.bioinfo.rpi.edu/cgi-bin/dna-form1.cgi; [34]) was used to search for stem-loop structures, and the MEME server (http://meme.sdsc.edu/meme4_4_0/intro.html) was used to search for sequence motifs that might be associated with inversions.

Patterns and Rates of Microinversion Evolution

Microinversions were coded as binary characters, and PAUP* 4.0b10 [35] was used to calculate numbers of inversion events using maximum-parsimony (MP) and the Hackett et al. [27] topology. λMI was expressed as microinversions Mb-1 MY-1 to facilitate comparison to other studies [6]. The null hypothesis of equal genome-wide microinversion rates was tested as described by Han et al. [36]. Briefly, a global Poisson model (which assumes equal genome-wide rates) was used as the null hypothesis, and the fit of that null model was compared to that of the more general negative binomial (NB) model (which permits variation in λMI) using a likelihood ratio test (LRT). See Additional file 2 for details.

Phylogenetic Analyses

Phylogenetic analyses of the CLTC alignment, conducted to provide an estimate of the CLTC gene tree, used RAxML 7.0.4 [37]. Microinversions and sites with gaps and/or missing data in more than 50% of taxa were excluded before conducting the RAxML search. See Additional file 2 for details.

Results and Discussion

Many Avian Microinversions were Identified

Manual and automated searches revealed that non-coding regions associated with 11 of the 17 loci we examined contained microinversions (e.g., Figure 2) ranging from 5 bp to 38 bp (Additional file 2, Table S2). Their median length was 22 bp. A number of the microinversions identified here were much shorter than those reported in genome-scale comparisons of mammals [1, 16], where the smallest microinversions were 23 bp and 31 bp, respectively. Although it is possible that birds and mammals have distinct microinversion size spectra, it seems more likely that the large-scale surveys of mammalian data failed to identify the shortest microinversions.

If λMI was similar in birds and mammals, fewer than four microinversions would be expected given the amount of sequence data examined; instead, microinversions were identified at 49 positions (Table 1). Ma et al. [1] reported that short inversions are more common than long inversions. If this pattern continues as microinversions become even shorter than those they identified, the larger number of microinversions that we observed could reflect our identification of smaller inversions rather than any inherent difference between mammalian and avian genomes. The denser taxon sampling in our study, relative to whole genome studies in mammals, is also likely to have improved microinversion identification. Taken as a whole, our results suggest that previous studies that used mammalian data [1, 6] underestimated λMI.

The identification of microinversions can be difficult because point mutations and insertion-deletion events (indels) continue to accumulate after inversions. This has the potential to make ancient microinversions particularly difficult, or impossible, to identify. Denser taxon sampling can help by increasing the number of sequences closely related to those with the microinversion and by providing multiple versions of the inverted sequence (Additional file 2, Figure S1). Although the taxon sampling for this study was denser than previous surveys that used mammalian data, computational searches for microinversions were difficult. Many complementary strand alignments were not validated as actual inversions; the false positives reflected palindromes and other phenomena. bl2seq performed better than YASS, producing fewer false positives while still identifying all of the microinversions also found by YASS. However, even after employing two computational approaches, some microinversions were only identified "by eye" (Additional file 2, Table S2), suggesting that further improvements to the methods used to identify microinversions are required.

Most microinversions were assigned to terminal branches in the Hackett et al. [27] phylogeny (Figure 3) when the MP criterion was used. This raises the question of whether an acquisition bias caused us to miss a number of ancient microinversions that occurred closer to the base of the tree. However, the structure of the avian tree of life is dominated by a rapid radiation at the base of Neoaves, the most speciose avian supergroup (identified in Figure 3), leading to a tree dominated by terminal branches. Indeed, 70.8% of the overall treelength in the Hackett et al. ML tree [27] comprises terminal branches. The number of microinversions observed on terminal branches was not significantly different from expectation given the proportion of the tree that reflected internal and terminal branches (χ2 = 3.0; P = 0.08). Thus, acquisition bias did not have a major impact upon our ability to identify ancient inversions.

Figure 3
figure 3

Microinversions indicated on the Hackett et al. [27]phylogeny. Inversions in introns are indicated with tick marks (blue for no homoplasy, green for the homoplastic inversions in CLTC intron 6, and red for the homoplastic inversions in CLTC intron 7). The 3' UTR inversion the PCBD1, which was obtained from selected galliform (see Results and Discussion), is indicated with a blue diamond. This mapping of character state changes assumes a reversal to the ancestral state in Psittaciformes for the CLTC intron 7 microinversion (indicated by an X over the red tick mark). An inversion in CLTCL1 where Palaeognathae and Neognathae differ is shown along the root branch. Orders united by microinversions are indicated using names above the branch uniting them and brackets to the right. The order Galliformes is emphasized because 3' UTRs were sequenced from additional taxa in that order (see text). This phylogeny is presented as a cladogram because many internal branches are very short and this presentation makes it easier to locate the inversion events. For branch length information refer to Figure 3 in Hackett et al. [27] and the chronogram presented for this publication (Additional file 2, Figure S3).

Avian Microinversion Rates Vary Among Loci

Estimates of λMI differ among loci (Table 1). The Poisson model of microinversion accumulation (the null hypothesis) was rejected in favour of the NB model (which includes rate variation) using the LRT (2δlnL = 27.55; P < 10-6). Excluding the highest-rate loci (CLTC and CLTCL1) eliminated our ability to reject the Poisson model (2δlnL = 2.29; P = 0.13) and reduced the λMI estimate to 0.25 microinversions Mb-1 MY-1 (the value presented in Figure 1; 95% confidence interval of 0.17 - 0.36). This suggests a "hotspot" model in which CLTC and CLTCL1 are inversion-prone. However, even the lower estimate of λMI for "non-hotspot" loci greatly exceeded previous estimates of λMI, consistent with our hypothesis that the identification of microinversions, especially the shortest inversions, has been improved relative to prior studies.

Surprisingly, both hotspot loci encode clathrin heavy chains, which are proteins critical for endocytosis [38], suggesting that the high microinversion rates could reflect their functional similarities. However, these clathrin heavy chain paralogs arose by duplication early in vertebrate evolution [39], and the homologous introns in CLTC and CLTCL1 do not exhibit detectable sequence similarity. Although specific intronic motifs can be overrepresented in functionally related genes [40], motifs common to the CLTC and CLTCL1 introns were not identified (data not shown). This suggests that it will be necessary to identify additional hotspot loci to understand the basis for inversion hotspots.

Microinversions were absent in some loci (Table 1), but it is unclear whether this reflects stochastic variation or the existence of "coldspots". 3' UTRs are coldspot candidates because they exhibit a lower rate of sequence evolution than introns [29, 41] and they are known to include regulatory elements [42]. Many of these regulatory sequences are non-palindromic [43, 44] and are unlikely to remain functional after inversion. Two to three microinversions were expected in our 3' UTR data (assuming equal rates for non-hotspot loci), but none were identified. We examined 3' UTRs from five additional loci (ALDOB, CRYAA, EEF2, HMGN2, and PCBD1), four of which have intronic microinversions (Table 1), by examining 23 members of the avian order Galliformes [41]. A 36 bp microinversion is present in the Rollulus roulroul PCBD1 3' UTR, indicating that these regions are not absolutely refractory to microinversions. Thus, future surveys should include 3' UTRs to improve λMI estimates for those regions and establish whether they exhibit among-locus rate variation similar to introns.

Homoplastic and Overlapping Microinversions Exist

Two microinversions in CLTC appeared homoplastic because the inverted forms were present in divergent lineages (e.g., Additional File 2, Figure S2). These homoplastic microinversions required at least three (CLTC intron 6) or four (CLTC intron 7) changes on the Hackett et al. [27] phylogeny using the MP criterion to explain the observed distribution of character states (Figure 3). Errors in the phylogeny are unlikely to explain this observation, since the relevant branches are well supported (compare Figure 3 to Figure 2 of Hackett et al. [27]; also see Additional File 2, Figure S2). Moreover, when these microinversions were mapped on other recent estimates of avian phylogeny using the MP criterion they require similar levels of homoplasy. These other estimates of phylogeny are based upon nuclear [26, 45], mitochondrial [46–48], and morphological data [49, 50], as well as expert opinion (e.g., Figure 27.10 in Cracraft et al. [51] and Figure 5 in Mayr [52]).

Hemiplasy is unlikely to explain the observed homoplastic microinversions for two reasons. First, hemiplasy would require maintenance of polymorphic inversions over multiple, long internal branches (estimates of branch lengths are presented as a chronogram in Additional File 2, Figure S3). Second, the estimate of the CLTC gene tree was not consistent with the microinversion distribution (Additional file 2, Figure S4), even in the single case in which branch lengths are short enough that hemiplasy is plausible. Thus, the CLTC inversions reflect genuine homoplasy, not hemiplasy, a novel finding for microinversions in animal nuclear genomes.

In addition to the homoplastic microinversions in CLTC, we also found several overlapping microinversions (Additional file 2, Table S2). All of these overlapping microinversions reflected independent inversions in distinct lineages. We identified two overlapping microinversions in CLTC and one in CLTCL1; the two overlapping microinversions in CLTC (INV-14 and INV-15; see Additional file 2, Table S2) also overlapped with one of the homoplastic microinversions in CLTC (INV-13). Thus, there were at least 12 inversion events in four specific regions of the two hotspot loci. There were also two additional overlapping inversions in low-rate loci (EEF2 and IRF2). Neither the homoplastic nor the overlapping microinversions were associated with stem-loop motifs (e.g., Additional file 2, Figure S4) or any other motifs that could be identified using MEME. These homoplastic and overlapping microinversions indicate that the actual state space for microinversions is likely to be smaller than their potential state space.

Are Microinversions useful for Phylogenetics?

Although the existence of homoplastic microinversions demonstrates that they are not perfect characters, they still have the potential to be useful phylogenetic markers. The retention index of microinversions (RIMI = 0.949) given the Hackett et al. [27] tree is substantially higher than the retention index for nucleotide changes (RIintron = 0.52, RIcoding exon = 0.54, RIUTR = 0.58). Such low amount of homoplasy suggests that an appropriate analytical approach (that accommodates homoplasy and hemiplasy) should yield an accurate species tree given a sufficient number of inversions.

Branches at the base of Neoaves are very short and this radiation is a classic example of a "bush" phylogeny [27]. In fact, the base of Neoaves has even been suggested to be a "hard" polytomy [53]. Hard polytomies reflect genuine multiple speciation events, so they cannot be represented as bifurcating trees. Even if Neoaves is a "soft" polytomy, many branches are likely to be <1 MY in length (Additional File 2, Figure S3; also see [26, 45]). The low estimates of λMI imply that microinversions will seldom occur along these short branches. How much sequence data would be necessary to resolve internodes of this length using microinversions? Power analysis assuming 1 MY branch lengths using the rate estimate that excludes the hotspot loci [54] indicates ~1.2 Mbp of non-coding sequence per taxon is needed to find at least one informative inversion and ~12 Mbp per taxon to identify an inversion on a specific branch (Additional file 2, Table S3). This estimate is orders of magnitude larger than the amount needed for of conventional analyses of sequence data (cf. Chojnowski et al. [26]). Moreover, it is desirable to identify multiple informative inversions along internodes given the potential for hemiplasy and homoplasy, suggesting that the use of microinversions as the sole source of information to estimate a phylogeny similar to the avian tree of life would require even more data (Additional file 2, Table S3).

Microinversions and Multiple Sequence Alignment

The identification of microinversions is also important to ensure correct sequence alignment. Otherwise estimates of the amount of evolutionary change will be distorted, potentially resulting in incorrect phylogenetic estimation [19]. Algorithms for sequence alignment that include the possibility of inversions have been proposed [55–57], and they have the potential advantage of incorporating explicit penalties for inversion events. However, the optimal inversion penalty to limit false positives may be difficult to determine and the available algorithms are limited to the identification of non-overlapping microinversions. Overlapping microinversions were found at four loci that we examined, suggesting that the inability to identify overlapping inversions may represent a major limitation. Overlapping and homoplastic microinversions can be divided into three basic categories (Additional file 2, Figure S6), and the strategy we employed should be able to detect two of these categories efficiently. The third category (type III in Additional file 2, Figure S6, which corresponds to the case of multiple homoplastic or overlapping inversion events on a single branch) is expected to be rare. It may be possible to overcome this problem in a multiple sequence alignment framework using a divide-and-conquer approach by selecting subsets of taxa for which overlapping microinversions are less likely to be present. This would necessitate a subsequent assembly of the alignments. Moreover, such an approach might eliminate the benefits of dense taxon sampling. Despite these limitations, fully automated approaches could be less labour intensive than our approach. However, it is unclear whether microinversion identification can be fully automated since our results suggest that short microinversions may always require manual validation. Taken as a whole, these issues further emphasize the need to continue to improve algorithms for the detection and alignment of these interesting genomic changes.

Conclusions

These analyses demonstrate that the identification of microinversions is important, despite the relatively low rate of accumulation of these genomic changes. This study revealed that microinversions accumulate more rapidly in avian genomes than expected based upon prior analyses of mammalian genomes, although this difference is likely to reflect the failure to identify very short inversions in the large-scale comparisons of mammalian data. If this failure to identify short microinversion does explain the differences among this and previous studies, the estimates of λMI presented here, which are similar to the rate of accumulation of the most common type of avian TE insertion (Figure 1), may be more typical of vertebrate genomes. This likelihood that typical vertebrate λMI values may be higher than suggested by previous studies emphasizes the importance of understanding the impact of microinversions upon genome evolution. We also documented the existence of microinversion hotspots, suggesting that some regions of the genome are especially prone to these mutations. The identification of additional hotspots may provide information about the mechanistic basis of these mutations. Indeed, we were able to exclude one proposed mechanism, the existence of conserved stem-loops, based upon an examination of the inversion hotspots identified here. Despite our observation that microinversions can exhibit homoplasy, they are still relatively reliable RGCs and as such may define gene tree bipartitions more accurately than conventional sequence data (see Nishihara et al. [58]). In the future, analytical methods that integrate microinversions with sequence data and information about other RGCs (and incorporate the potential for both hemiplasy and homoplasy) will facilitate robust resolution of difficult nodes in the tree of life and provide additional insights into the mechanism(s) responsible for their accumulation over evolutionary time.

References

  1. Ma J, Zhang L, Suh BB, Raney BJ, Burhans RC, Kent WJ, Blanchette M, Haussler D, Miller W: Reconstructing contiguous regions of an ancestral genome. Genome Res. 2006, 16: 1557-1565. 10.1101/gr.5383506.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Rascol VL, Pontarotti P, Levasseur A: Ancestral animal genomes reconstruction. Curr Opin Immunol. 2007, 19: 542-546. 10.1016/j.coi.2007.06.009.

    Article  CAS  PubMed  Google Scholar 

  3. Levasseur A, Pontarotti P, Poch O, Thompson JD: Strategies for reliable exploitation of evolutionary concepts in high throughput biology. Evol Bioinform Online. 2008, 121-137. 4

  4. Kriegs JO, Churakov G, Kiefmann M, Jordan U, Brosius J, Schmitz J: Retroposed elements as archives for the evolutionary history of placental mammals. PLoS Biol. 2006, 4: e91-10.1371/journal.pbio.0040091.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Boore JL: The use of genome-level characters for phylogenetic reconstruction. Trends Ecol Evol. 2006, 21: 439-446. 10.1016/j.tree.2006.05.009.

    Article  PubMed  Google Scholar 

  6. Chaisson MJ, Raphael BJ, Pevzner PA: Microinversions in mammalian evolution. Proc Natl Acad Sci USA. 2006, 103: 19824-19829. 10.1073/pnas.0603984103.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Krauss V, Thümmler C, Georgi F, Lehmann J, Stadler PF, Eisenhardt C: Near intron positions are reliable phylogenetic markers: An application to holometabolous insects. Mol Biol Evol. 2008, 25: 821-830. 10.1093/molbev/msn013.

    Article  CAS  PubMed  Google Scholar 

  8. Rogozin IB, Thomson K, Csürös M, Carmel L, Koonin EV: Homoplasy in genome-wide analysis of rare amino acid replacements: The molecular-evolutionary basis for Vavilov's law of homologous series. Biol Direct. 2008, 3: 7-10.1186/1745-6150-3-7.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Rokas A, Holland PWH: Rare genomic changes as a tool for phylogenetics. Trends Ecol Evol. 2000, 15: 454-459. 10.1016/S0169-5347(00)01967-4.

    Article  PubMed  Google Scholar 

  10. Hillis DM: SINEs of the perfect character. Proc Natl Acad Sci USA. 1999, 96: 9979-9981. 10.1073/pnas.96.18.9979.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Avise JC, Robinson TJ: Hemiplasy: a new term in the lexicon of phylogenetics. Syst Biol. 2008, 57: 503-507. 10.1080/10635150802164587.

    Article  PubMed  Google Scholar 

  12. Pamilo P, Nei M: Relationships between gene trees and species trees. Mol Biol Evol. 1988, 5: 568-583.

    CAS  PubMed  Google Scholar 

  13. Moore WS: Inferring phylogenies from mtDNA variation: Mitochondrial-gene trees versus nuclear-gene trees. Evolution. 1995, 49: 718-726. 10.2307/2410325.

    Article  Google Scholar 

  14. Rokas A, Carroll SB: Bushes in the Tree of Life. PLoS Biol. 2006, 4: e352-10.1371/journal.pbio.0040352.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Koonin EV: The Biological Big Bang model for the major transitions in evolution. Biol Direct. 2007, 2: 21-10.1186/1745-6150-2-21.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Feuk L, MacDonald JR, Tang T, Carson AR, Li M, Rao G, Khaja R, Scherer SW: Discovery of human inversion polymorphisms by comparative analysis of human and chimpanzee DNA sequence assemblies. PLoS Genet. 2005, 1: e56-10.1371/journal.pgen.0010056.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Kelchner SA, Wendel JF: Hairpins create minute inversions in non-coding regions of chloroplast DNA. Curr Genet. 1996, 30: 259-262. 10.1007/s002940050130.

    Article  CAS  PubMed  Google Scholar 

  18. Harshman J, Huddleston CJ, Bollback JP, Parsons TJ, Braun MJ: True and false gharials: A nuclear gene phylogeny of crocodylia. Syst Biol. 2003, 52: 386-402. 10.1080/10635150390197028.

    Article  PubMed  Google Scholar 

  19. Kim KJ, Lee HL: Widespread occurrence of small inversions in the chloroplast genomes of land plants. Mol Cells. 2005, 19: 104-113.

    CAS  PubMed  Google Scholar 

  20. Kimball RT, Braun EL: A multigene phylogeny of Galliformes supports a single origin of erectile ability in non-feathered facial traits. J Avian Biol. 2008, 39: 438-445.

    Google Scholar 

  21. Catalano SA, Saidman BO, Vilardi JC: Evolution of small inversions in chloroplast genome: a case study from a recurrent inversion in angiosperms. Cladistics. 2009, 25: 93-104. 10.1111/j.1096-0031.2008.00236.x.

    Article  Google Scholar 

  22. Stern DB, Jones H, Gruissem W: Function of plastid mRNA 3' inverted repeats. RNA stabilization and gene-specific protein binding. J Biol Chem. 1989, 264: 18742-18750.

    CAS  PubMed  Google Scholar 

  23. Hugo H, Cures A, Suraweera N, Drabsch Y, Purcell D, Mantamadiotis T, Phillips W, Dobrovic A, Zupi G, Gonda TJ, Iacopetta B, Ramsay RG: Mutations in the MYB intron I regulatory sequence increase transcription in colon cancers. Genes Chromosomes Cancer. 2006, 45: 1143-1154. 10.1002/gcc.20378.

    Article  CAS  PubMed  Google Scholar 

  24. Rott R, Liveanu V, Drager RG, Stern DB, Schuster G: The sequence and structure of the 3'-untranslated regions of chloroplast transcripts are important determinants of mRNA accumulation and stability. Plant Mol Biol. 1998, 36: 307-314. 10.1023/A:1005943701253.

    Article  CAS  PubMed  Google Scholar 

  25. Macdonald SJ, Long AD: Fine scale structural variants distinguish the genomes of Drosophila melanogaster and D. pseudoobscura. Genome Biol. 2006, 7: R67-10.1186/gb-2006-7-7-r67.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Chojnowski JL, Kimball RT, Braun EL: Introns outperform exons in analyses of basal avian phylogeny using clathrin heavy chain genes. Gene. 2008, 410: 89-96. 10.1016/j.gene.2007.11.016.

    Article  CAS  PubMed  Google Scholar 

  27. Hackett SJ, Kimball RT, Reddy S, Bowie RCK, Braun EL, Braun MJ, Chojnowski JL, Cox WA, Han KL, Harshman J, Huddleston C, Marks BD, Miglia KJ, Moore WS, Sheldon FH, Steadman DW, Witt CC, Yuri T: A phylogenomic study of birds reveals their evolutionary history. Science. 2008, 320: 1763-1768. 10.1126/science.1157704.

    Article  CAS  PubMed  Google Scholar 

  28. Harshman J, Braun EL, Braun MJ, Huddleston CJ, Bowie RCK, Chojnowski JL, Hackett SJ, Han KL, Kimball RT, Marks BD, Miglia KJ, Moore WS, Reddy S, Sheldon FH, Steppan SJ, Witt CC, Yuri T: Phylogenomic evidence for multiple losses of flight in ratite birds. Proc Natl Acad Sci USA. 2008, 105: 13462-13467. 10.1073/pnas.0803242105.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Kimball RT, Braun EL, Barker FK, Bowie RCK, Braun MJ, Chojnowski JL, Hackett SJ, Han KL, Harshman J, Heimer-Torres V, Holznagel W, Huddleston CJ, Marks BD, Miglia KJ, Moore WS, Reddy S, Sheldon FH, Smith JV, Witt CC, Yuri T: A well-tested set of primers to amplify regions spread across the avian genome. Mol Phylogenet Evol. 2009, 50: 654-660. 10.1016/j.ympev.2008.11.018.

    Article  CAS  PubMed  Google Scholar 

  30. Chenna R, Sugawara H, Koike T, Lopez R, Gibson TJ, Higgins DG, Thompson JD: Multiple sequence alignment with the Clustal series of programs. Nucleic Acids Res. 2003, 31: 3497-3500. 10.1093/nar/gkg500.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Katoh M, Kuma M: MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002, 30: 3059-3066. 10.1093/nar/gkf436.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Tatusova TA, Madden TL: BLAST 2 Sequences, a new tool for comparing protein and nucleotide sequences. FEMS Microbiol Lett. 1999, 174: 247-250. 10.1111/j.1574-6968.1999.tb13575.x.

    Article  CAS  PubMed  Google Scholar 

  33. Noé L, Kucherov G: YASS: enhancing the sensitivity of DNA similarity search. Nucleic Acids Res. 2005, 33: W540-W543. 10.1093/nar/gki478.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31: 3406-3415. 10.1093/nar/gkg595.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Swofford DL: PAUP*: Phylogenetic analysis using parsimony (*and other methods), Version 4. 2003, Sunderland, MA: Sinauer Associates

    Google Scholar 

  36. Han KL, Braun EL, Kimball RT, Reddy S, Bowie RCK, Braun MJ, Chojnowski JL, Hackett SJ, Harshman J, Huddleston CJ, Marks BD, Miglia KJ, Moore WS, Sheldon FH, Steadman DW, Witt CC, Yuri T: Are transposable element insertions homoplasy free? An examination using the avian tree of life. Syst Biol. 2011, 60: 375-386. 10.1093/sysbio/syq100.

    Article  PubMed  Google Scholar 

  37. Stamatakis A: Maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006, 22: 2688-2690. 10.1093/bioinformatics/btl446.

    Article  CAS  PubMed  Google Scholar 

  38. Kirchhausen T: Clathrin. Annu Rev Biochem. 2000, 69: 699-727. 10.1146/annurev.biochem.69.1.699.

    Article  CAS  PubMed  Google Scholar 

  39. Wakeham DE, Abi-Rached L, Towler MC, Wilbur JD, Parham P, Brodsky FM: Clathrin heavy and light chain isoforms originated by independent mechanisms of gene duplication during chordate evolution. Proc Natl Acad Sci USA. 2005, 102: 7209-7214. 10.1073/pnas.0502058102.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Tsirigos A, Rigoutsos I: Human and mouse introns are linked to the same processes and functions through each genome's most frequent non-conserved motifs. Nucleic Acids Res. 2008, 36: 3484-3493. 10.1093/nar/gkn155.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Bonilla AJ, Braun EL, Kimball RT: Comparative molecular evolution and phylogenetic utility of 3'-UTRs and introns in Galliformes. Mol Phylogenet Evol. 2010, 56: 536-542. 10.1016/j.ympev.2010.04.006.

    Article  CAS  PubMed  Google Scholar 

  42. Warren WC, Clayton DF, Ellegren H, Arnold AP, Hillier LW, Kunstner A, Searle S, White S, Vilella AJ, Fairley S, Heger A, Kong L, Ponting CP, Jarvis ED, Mello CV, Minx P, Lovell P, Velho TAF, Ferris M, Balakrishnan CN, Sinha S, Blatti C, London SE, Li Y, Lin YC, George J, Sweedler J, Southey B, Gunaratne P, Watson M, et al: The genome of a songbird. Nature. 2010, 464: 757-762. 10.1038/nature08819.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Khabar KSA: The AU-rich transcriptome: More than interferons and cytokines, and its role in disease. J Interferon Cytokine Res. 2005, 25: 1-10. 10.1089/jir.2005.25.1.

    Article  CAS  PubMed  Google Scholar 

  44. Chen JM, Férec C, Cooper DN: A systematic analysis of disease-associated variants in the 3' regulatory regions of human protein-coding genes II: The importance of mRNA secondary structure in assessing the functionality of 3' UTR variants. Hum Genet. 2006, 120: 301-333. 10.1007/s00439-006-0218-x.

    Article  CAS  PubMed  Google Scholar 

  45. Ericson PGP, Anderson CL, Britton T, Elzanowki A, Johansson US, Källersjö M, Ohlson JI, Parsons TJ, Zuccon D, Mayr G: Diversification of Neoaves: Integration of molecular sequence data and fossils. Biol Lett. 2006, 2: 543-547. 10.1098/rsbl.2006.0523.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Gibb GC, Kardailsky O, Kimball RT, Braun EL, Penny D: Mitochondrial genomes and avian phylogeny: Complex characters and resolvability without explosive radiations. Mol Biol Evol. 2007, 24: 269-280.

    Article  CAS  PubMed  Google Scholar 

  47. Brown JW, Payne RB, Mindell DP: Nuclear DNA does not reconcile 'rocks' and 'clocks' in Neoaves: a comment on Ericson et al. Biol Lett. 2007, 3: 257-259. 10.1098/rsbl.2006.0611.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Pratt RC, Gibb GC, Morgan-Richards M, Phillips MJ, Hendy MD, Penny D: Toward resolving deep Neoaves phylogeny: Data, signal enhancement, and priors. Mol Biol Evol. 2009, 26: 313-326. 10.1093/molbev/msn248.

    Article  CAS  PubMed  Google Scholar 

  49. Mayr G, Clarke J: The deep divergences of neornithine birds: A phylogenetic analysis of morphological characters. Cladistics. 2003, 19: 527-553. 10.1111/j.1096-0031.2003.tb00387.x.

    Article  Google Scholar 

  50. Livezey BC, Zusi RL: Higher-order phylogeny of modern birds (Theropoda, Aves: Neornithes) based on comparative anatomy. II. Analysis and discussion. Zool J Linn Soc. 2007, 149: 1-95.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Cracraft J, Barker FK, Braun M, Harshman J, Dyke GJ, Feinstein J, Stanley S, Cibois A, Schikler P, Beresford P, García-Moreno J, Yuri T, Mindell DP: Phylogenetic relationships among modern birds (Neornithes): Toward an avian tree of life. Assembling the tree of life. Edited by: Cracraft J, Donoghue MJ. 2004, New York: Oxford University Press, 468-489.

    Google Scholar 

  52. Mayr G: Metaves, Mirandornithes, Strisores and other novelties - a critical review of the higher-level phylogeny of neornithine birds. J Zool Syst Evol Res. 2011, 49: 58-76.

    Article  Google Scholar 

  53. Poe S, Chubb AL: Birds in a bush: five genes indicate explosive evolution of avian orders. Evolution. 2004, 58: 404-415.

    Article  PubMed  Google Scholar 

  54. Braun EL, Kimball RT: Polytomies, the power of phylogenetic inference, and the stochastic nature of molecular evolution: A comment on Walsh et al. (1999). Evolution. 2001, 55: 1261-1263.

    Article  CAS  PubMed  Google Scholar 

  55. Schöniger M, Waterman MS: A local algorithm for DNA sequence alignment with inversions. Bull Math Biol. 1992, 54: 521-536.

    Article  PubMed  Google Scholar 

  56. Vellozo AF, Alves CER, do Lago AP: Alignment with non-overlapping inversions in O(n3)-time. Lect Notes Comput Sc. 2006, 4175: 186-196. 10.1007/11851561_18.

    Article  Google Scholar 

  57. Ledergerber C, Dessimoz C: Alignments with non-overlapping moves, inversions and tandem duplications in O(n4) time. J Comb Optim. 2008, 16: 263-278. 10.1007/s10878-007-9132-y.

    Article  Google Scholar 

  58. Nishihara H, Maruyama S, Okada N: Retroposon analysis and recent geological data suggest near-simultaneous divergence of the three superorders of mammals. Proc Natl Acad Sci USA. 2009, 106: 5235-5240. 10.1073/pnas.0809297106.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Mindell DP, Knight A, Baer C, Huddleston CJ: Slow rates of molecular evolution in birds and the metabolic rate and body temperature hypotheses. Mol Biol Evol. 1996, 13: 422-426.

    Article  CAS  Google Scholar 

  60. Cooper GM, Brudno M, NISC Comparative Sequencing Program, Green ED, Batzoglou S, Sidow A: Quantitative estimates of sequence divergence for comparative analyses of mammalian genomes. Genome Res. 2003, 13: 813-820. 10.1101/gr.1064503.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We thank Clare Rittschof and members of the Kimball-Braun lab group for helpful comments on this manuscript and we are grateful to the museums and collectors (Additional file 1) for loaning samples. This work was supported by the U.S. National Science Foundation Assembling the Tree of Life program, grants DEB-0228682 (to RTK, ELB, and David W Steadman), DEB-0228675 (to SJH), DEB-0228688 (to FHS), and DEB-0228617 (to WSM). Publication of this article was funded in part by the University of Florida Open-Access Publishing Fund.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Edward L Braun.

Additional information

Authors' contributions

ELB designed the study, wrote many of the computer programs, conducted analyses, and wrote the manuscript. RTK helped design the study, validated computational microinversion searches, and helped draft the manuscript. K-LH manually identified the homoplastic microinversions and some overlapping microinversions. NRI-V contributed to analyses and wrote a novel computer program. AJB, JLC, JVS, RCKB, MJB, JH, SR, CCW, and TY helped to draft and revise the manuscript. All authors (except NRI-V) contributed to the sequence data analysed, helped construct and edit alignments, and manually searched the alignments for microinversions. All authors read and approved the final manuscript.

Electronic supplementary material

12862_2010_1746_MOESM1_ESM.XLS

Additional file 1: Taxon list. List of the taxa used for this analysis and the accession numbers for the novel CLTCL1 sequences collected for this study, in Microsoft Excel format. (XLS 98 KB)

12862_2010_1746_MOESM2_ESM.PDF

Additional file 2: Supplementary information. Six figures, three tables, and supplementary methods (including the details of the literature survey used to estimate the rates of various types of genomic changes and the power analysis described in the main text), in pdf format. (PDF 6 MB)

12862_2010_1746_MOESM3_ESM.XLS

Additional file 3: Details of a microinversion search. An example of a microinversion search (of TPM1 intron 6) is presented along with a description of the search algorithm using pseudocode, in Microsoft Excel format. (XLS 72 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to 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.

Reprints and permissions

About this article

Cite this article

Braun, E.L., Kimball, R.T., Han, KL. et al. Homoplastic microinversions and the avian tree of life. BMC Evol Biol 11, 141 (2011). https://doi.org/10.1186/1471-2148-11-141

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2148-11-141

Keywords