Multiple acquisitions via horizontal transfer of a group I intron in the mitochondrial cox1 gene during evolution of the Araceae family.

A group I intron has recently been shown to have invaded mitochondrial cox1 genes by horizontal transfer many times during the broad course of angiosperm evolution. To investigate the frequency of acquisition of this intron within a more closely related group of plants, we determined its distribution and inferred its evolutionary history among 14 genera of the monocot family Araceae. Southern blot hybridizations showed that 6 of the 14 genera contain this intron in their cox1 genes. Nucleotide sequencing showed that these six introns are highly similar in sequence (97.7%-99.4% identity) and identical in length (966 nt). Phylogenetic evidence from parsimony reconstructions of intron distribution and phylogenetic analyses of intron sequences is consistent with a largely vertical history of intron transmission in the family; the simplest scenarios posit but one intron gain and two losses. Despite this, however, striking differences in lengths of exonic co-conversion tracts, coupled with the absence of co-conversion in intron-lacking taxa, indicate that the six intron-containing Araceae probably acquired their introns by at least three and quite possibly five separate horizontal transfers. The highly similar nature of these independently acquired introns implies a closely related set of donor organisms.     

Nuclear-encoded rDNA group I introns: origin and phylogenetic relationships of insertion site lineages in the green algae

Group I introns are widespread in eukaryotic organelles and nuclear- encoded ribosomal DNAs (rDNAs). The green algae are particularly rich in rDNA group I introns. To better understand the origins and phylogenetic relationships of green algal nuclear-encoded small subunit rDNA group I introns, a secondary structure-based alignment was constructed with available intron sequences and 11 new subgroup ICI and three new subgroup IB3 intron sequences determined from members of the Trebouxiophyceae (common phycobiont components of lichen) and the Ulvophyceae. Phylogenetic analyses using a weighted maximum-parsimony method showed that most group I introns form distinct lineages defined by insertion sites within the SSU rDNA. The comparison of topologies defining the phylogenetic relationships of 12 members of the 1512 group I intron insertion site lineage (position relative to the E. coli SSU rDNA coding region) with that of the host cells (i.e., SSU rDNAs) that contain these introns provided insights into the possible origin, stability, loss, and lateral transfer of ICI group I introns. The phylogenetic data were consistent with a viral origin of the 1512 group I intron in the green algae. This intron appears to have originated, minimally, within the SSU rDNA of the common ancestor of the trebouxiophytes and has subsequently been vertically inherited within this algal lineage with loss of the intron in some taxa. The phylogenetic analyses also suggested that the 1512 intron was laterally transferred among later-diverging trebouxiophytes; these algal taxa may have coexisted in a developing lichen thallus, thus facilitating cell- to-cell contact and the lateral transfer. Comparison of available group I intron sequences from the nuclear-encoded SSU rDNA of phycobiont and mycobiont components of lichens demonstrated that these sequences have independent origins and are not the result of lateral transfer from one component to the other.      

Evolution of rbcL group IA introns and intron open reading frames within the colonial Volvocales (Chlorophyceae)

Mobile group I introns sometimes contain an open reading frame (ORF) possibly encoding a site-specific DNA endonuclease. However, previous phylogenetic studies have not clearly deduced the evolutionary roles of the group I intron ORFs. In this paper, we examined the phylogeny of group IA2 introns inserted in the position identical to that of the chloroplast-encoded rbcL coding region (rbcL-462 introns) and their ORFs from 13 strains of five genera (Volvox, Pleodorina, Volvulina, Astrephomene, and Gonium) of the colonial Volvocales (Chlorophyceae) and a related unicellular green alga, Vitreochlamys. The rbcL-462 introns contained an intact or degenerate ORF of various sizes except for the Gonium multicoccum rbcL-462 intron. Partial amino acid sequences of some rbcL-462 intron ORFs exhibited possible homology to the endo/excinuclease amino acid terminal domain. The distribution of the rbcL-462 introns is sporadic in the phylogenetic trees of the colonial Volvocales based on the five chloroplast exon sequences (6021 bp). Phylogenetic analyses of the conserved intron sequences resolved that the G. multicoccum rbcL-462 intron had a phylogenetic position separate from those of other colonial volvocalean rbcL-462 introns, indicating the recent horizontal transmission of the intron in the G. multicoccum lineage. However, the combined data set from conserved intron sequences and ORFs from most of the rbcL-462 introns resolved robust phylogenetic relationships of the introns that were consistent with those of the host organisms. Therefore, most of the extant rbcL-462 introns may have been vertically inherited from the common ancestor of their host organisms, whereas such introns may have been lost in other lineages during evolution of the colonial Volvocales. In addition, apparently higher synonymous substitutions than nonsynonymous substitutions in the rbcL-462 intron ORFs indicated that the ORFs might evolve under functional constraint, which could result in homing of the rbcL-462 intron in cases of spontaneous intron loss. On the other hand, the presence of intact to largely degenerate ORFs of the rbcL-462 introns within the three isolates of Gonium viridistellatum and the rare occurrence of the ORF-lacking rbcL-462 intron suggested that the ORFs might degenerate to result in the spontaneous intron loss during a very short evolutionary time following the loss of the ORF function. Thus, the sporadic distribution of the rbcL-462 introns within the colonial Volvocales can be largely explained by an equilibrium between maintenance of the introns by the intron ORF and spontaneous loss of introns when the introns do not have a functional ORF.     

Wagner and Dollo: a stochastic duet by composing two parsimonious solos

New contributions toward generalizing evolutionary models expand greatly our ability to analyze complex evolutionary characters and advance phylogeny reconstruction. In this article, we extend the binary stochastic Dollo model to allow for multi-state characters. In doing so, we align previously incompatible Wagner and Dollo parsimony principles under a common probabilistic framework by embedding arbitrary continuous-time Markov chains into the binary stochastic Dollo model. This approach enables us to analyze character traits that exhibit both Dollo and Wagner characteristics throughout their evolutionary histories. Utilizing Bayesian inference, we apply our novel model to analyze intron conservation patterns and the evolution of alternatively spliced exons. The generalized framework we develop demonstrates potential in distinguishing between phylogenetic hypotheses and providing robust estimates of evolutionary rates. Moreover, for the two applications analyzed here, our framework is the first to provide an adequate stochastic process for the data. We discuss possible extensions to the framework from both theoretical and applied perspectives.     

Evidence of the presence of tRNAfMet group I intron in the marine cyanobacterium Synechococcus elongatus

Self-splicing group I introns in tRNA anticodon loops have been found in diverse groups of bacteria. In this work, we identified tRNAfMet group I introns in six strains of marine Synechococcus elongatus. Introns with sizes around 280 bp were consistently obtained in all the strains tested. In a phylogenetic analysis using the nucleotide sequence determined in this study with other cyanobacterial tRNAfMet and tRNALeu intron sequences, the Synechococcus sequence was grouped together with the sequences from other unicellular cyanobacterial strains. Interestingly, the phylogenetic tree inferred from the intronic sequences clearly separates the different tRNA introns, suggesting that each family has its own evolutionary history.     

Exclusive conservation of mitochondrial group II intron nad4i548 among liverworts and its use for phylogenetic studies in this ancient plant clade

Liverworts occupy a pivotal position in land plant (embryophyte) phylogeny as the presumed earliest-branching major clade, sister to all other land plants, including the mosses, hornworts, lycophytes, monilophytes and seed plants. Molecular support for this earliest dichotomy in land plant phylogeny comes from strikingly different occurrences of introns in mitochondrial genes distinguishing liverworts from all other embryophytes. Exceptionally, however, the nad5 gene--the mitochondrial locus hitherto used most widely to elucidate early land plant phylogeny--carries a group I type intron that is shared between liverworts and mosses. We here explored whether a group II intron, the other major type of organellar intron, would similarly be conserved in position across the entire diversity of extant liverworts and could be of use for phylogenetic analyses in this supposedly most ancient embryophyte clade. To this end, we investigated the nad4 gene as a candidate locus possibly featuring different introns in liverworts as opposed to the non-liverwort embryophyte (NLE) lineage. We indeed found group II intron nad4i548 universally conserved in a wide phylogenetic sampling of 55 liverwort taxa, confirming clade specificity and surprising evolutionary stability of plant mitochondrial introns. As expected, intron nad4i548g2 carries phylogenetic information in its variable sequences, which confirms and extends previous cladistic insights on liverwort evolution. We integrate the new nad4 data with those of the previously established mitochondrial nad5 and the chloroplast rbcL and rps4 genes and present a phylogeny based on the fused datasets. Notably, the phylogenetic analyses suggest a reconsideration of previous phylogenetic and taxonomic assignments for the genera Calycularia and Mylia and resolve a sister group relationship of Ptilidiales and Porellales.     

Dynamic Insertion–Deletion of Introns in Deuterostome EF-1α Genes

To test the validity of intron–exon structure as a phylogenetic marker, the intron–exon structure of EF-1α genes was investigated for starfish, acornworms, ascidians, larvaceans, and amphioxus and compared with that of vertebrates. Of the 11 distinct intron insertion sites found within the coding regions of the deuterostome EF-1α genes, 7 are shared by several taxa, while the remainder are unique to certain taxa. Examination of the shared introns of the deuterostome EF-1α gene revealed that independent intron loss or intron insertion must have occurred in separate lineages of the deuterostome taxa. Maximum parsimony analysis of the intron–exon data matrix recovered five parsimonious trees (consistency index = 0.867). From this result, we concluded that the intron–exon structure of deuterostome EF-1α has evolved more dynamically than previously thought, rendering it unsuitable as a phylogenetic marker. We also reconstructed an evolutionary history of intron insertion–deletion events on the deuterostome phylogeny, based on several molecular phylogenetic studies. These analyses revealed that the deuterostome EF-1α gene has lost individual introns more frequently than all introns simultaneously.     

Origin and evolution of group I introns in cyanobacterial tRNA genes.

Many tRNA(Leu)UAA genes from plastids contain a group I intron. An intron is also inserted in the same gene at the same position in cyanobacteria, the bacterial progenitors of plastids, suggesting an ancient bacterial origin for this intron. A group I intron has also been found in the tRNA(fMet) gene of some cyanobacteria but not in plastids, suggesting a more recent origin for this intron. In this study, we investigate the phylogenetic distributions of the two introns among cyanobacteria, from the earliest branching to the more derived species. The phylogenetic distribution of the tRNA(Leu)UAA intron follows the clustering of rRNA sequences, being either absent or present in clades of closely related species, with only one exception in the Pseudanabaena group. Our data support the notion that the tRNA(Leu)UAA intron was inherited by cyanobacteria and plastids through a common ancestor. Conversely, the tRNA(fMet) intron has a sporadic distribution, implying that many gains and losses occurred during cyanobacterial evolution. Interestingly, a phylogenetic tree inferred from intronic sequences clearly separates the different tRNA introns, suggesting that each family has its own evolutionary history.     

A phylogeny of the damselfly genus calopteryx (Odonata) using mitochondrial 16S rDNA markers

We seek to reconstruct the phylogenetic relationships of the damselfly genus Calopteryx, for which extensive behavioral and morphological knowledge already exists. To date, analyses of the evolutionary pathways of different life history traits have been hampered by the absence of a robust phylogeny based on morphological data. In this study, we concentrate on establishing phylogenetic information from parts of the 16S rDNA gene, which we sequenced for nine Calopteryx species and five outgroup species. The mt 16S rDNA data set did not show signs of saturated variation for ingroup taxa, and phylogenetic reconstructions were insensitive to variation of outgroup taxa. Parsimony, neighbor-joining, and maximum-likelihood reconstructions agreed on parts of the tree. A consensus tree summarizes the significant results and indicates problematic nodes. The 16S rDNA sequences support monophyly of the genera Mnais, Matrona, and Calopteryx. However, the genus Calopteryx may not be monophyletic, since Matrona basilaris and Calopteryx atrata are sister taxa under every parameter setting. The North American and European taxa each appear as monophyletic clades, while the Asian Calopteryx atrata and Calopteryx cornelia are not monophyletic. Our data implies a different paleobiogeographic history of the Eurasian and North American species, with extant Eurasian species complexes shaped by glacial periods, in contrast to extant North American species groups.     

Spliceosomal introns as tools for genomic and evolutionary analysis

Over the past 5 years, the availability of dozens of whole genomic sequences from a wide variety of eukaryotic lineages has revealed a very large amount of information about the dynamics of intron loss and gain through eukaryotic history, as well as the evolution of intron sequences. Implicit in these advances is a great deal of information about the structure and evolution of surrounding sequences. Here, we review the wealth of ways in which structures of spliceosomal introns as well as their conservation and change through evolution may be harnessed for evolutionary and genomic analysis. First, we discuss uses of intron length distributions and positions in sequence assembly and annotation, and for improving alignment of homologous regions. Second, we review uses of introns in evolutionary studies, including the utility of introns as indicators of rates of sequence evolution, for inferences about molecular evolution, as signatures of orthology and paralogy, and for estimating rates of nucleotide substitution. We conclude with a discussion of phylogenetic methods utilizing intron sequences and positions.     

Evolutionary implications of intron–exon distribution and the properties and sequences of the RPL10A gene in eukaryotes

The RPL10A gene encodes the RPL10 protein, required for joining 40S and 60S subunits into a functional 80S ribosome. This highly conserved gene, ubiquitous across all eukaryotic super-groups, is characterized by a variable number of spliceosomal introns, present in most organisms. These properties facilitate the recognition of orthologs among distant taxa and thus comparative studies of sequences as well as the distribution and properties of introns in taxonomically distant groups of eukaryotes. The present study examined the multiple ways in which RPL10A conservation vs. sequence changes in the gene over the course of evolution, including in exons, introns, and the encoded proteins, can be exploited for evolutionary analysis at different taxonomic levels. At least 25 different positions harboring introns within the RPL10A gene were determined in different taxa, including animals, plants, fungi, and alveolates. Generally, intron positions were found to be well conserved even across different kingdoms. However, certain introns seemed to be restricted to specific groups of organisms. Analyses of several properties of introns, including insertion site, phase, and length, along with exon and intron GC content and exon–intron boundaries, suggested biases within different groups of organisms. The use of a standard primer pair to analyze a portion of the intron-containing RPL10A gene in 12 genera of green algae within Chlorophyta is presented as a case study for evolutionary analyses of introns at intermediate and low taxonomic levels. Our study shows that phylogenetic reconstructions at different depths can be achieved using RPL10A nucleotide sequences from both exons and introns as well as the amino acid sequences of the encoded protein.     

Introns outperform exons in analyses of basal avian phylogeny using clathrin heavy chain genes

Neoaves is the most diverse major avian clade, containing ~95% of avian species, and it underwent an ancient but rapid diversification that has made resolution of relationships at the base of the clade difficult. In fact, Neoaves has been suggested to be a "hard" polytomy that cannot be resolved with any amount of data. However, this conclusion was based on slowly evolving coding sequences and ribosomal RNAs and some recent studies using more rapidly evolving intron sequences have suggested some resolution at the base of Neoaves. To further examine the utility of introns and exons for phylogenetics, we sequenced parts of two unlinked clathrin heavy chain genes (CLTC and CLTCL1). Comparisons of phylogenetic trees based upon individual partitions (i.e. introns and exons), the combined dataset, and published phylogenies using Robinson-Foulds distances (a metric of topological differences) revealed more similarity than expected by chance, suggesting there is structure at the base of Neoaves. We found that introns provided more informative sites, were subject to less homoplasy, and provided better support for well-accepted clades, suggesting that intron evolution is better suited to determining closely-spaced branching events like the base of Neoaves. Furthermore, phylogenetic power analyses indicated that existing molecular datasets for birds are unlikely to provide sufficient phylogenetic information to resolve relationships at the base of Neoaves, especially when comprised of exon or other slowly evolving regions. Although relationships among the orders in Neoaves cannot be definitively established using available data, the base of Neoaves does not appear to represent a hard polytomy. Our analyses suggest that large intron datasets have the best potential to resolve relationships among avian orders and indicate that the utility of intron data for other phylogenetic questions should be examined.     

Phylogenetic trees based on gene content

Comparing gene content between species can be a useful approach for reconstructing phylogenetic trees. In this paper, we derive a maximum-likelihood estimation of evolutionary distance between species under a simple model of gene genesis and gene loss. Using simulated data on a biological tree with 107 taxa (and on a number of randomly generated trees), we compare the accuracy of tree reconstruction using this ML distance measure to an earlier ad hoc distance. We then compare these distance-based approaches to a character-based tree reconstruction method (Dollo parsimony) which seems well suited to the analysis of gene content data. To simplify simulations, we give a formal proof of the well-known 'fact' that the Dollo parsimony score is independent of the choice of root. Our results show a consistent trend, with the character-based method and ML distance measure outperforming the earlier ad hoc distance method. AVAILABILITY: http://www.ab.informatik.uni-tuebingen.de/software/genecontent/welcome_en.html     

Evolutionary convergence on highly-conserved 3' intron structures in intron-poor eukaryotes and insights into the ancestral eukaryotic genome

The presence of spliceosomal introns in eukaryotes raises a range of questions about genomic evolution. Along with the fundamental mysteries of introns' initial proliferation and persistence, the evolutionary forces acting on intron sequences remain largely mysterious. Intron number varies across species from a few introns per genome to several introns per gene, and the elements of intron sequences directly implicated in splicing vary from degenerate to strict consensus motifs. We report a 50-species comparative genomic study of intron sequences across most eukaryotic groups. We find two broad and striking patterns. First, we find that some highly intron-poor lineages have undergone evolutionary convergence to strong 3' consensus intron structures. This finding holds for both branch point sequence and distance between the branch point and the 3' splice site. Interestingly, this difference appears to exist within the genomes of green alga of the genus Ostreococcus, which exhibit highly constrained intron sequences through most of the intron-poor genome, but not in one much more intron-dense genomic region. Second, we find evidence that ancestral genomes contained highly variable branch point sequences, similar to more complex modern intron-rich eukaryotic lineages. In addition, ancestral structures are likely to have included polyT tails similar to those in metazoans and plants, which we found in a variety of protist lineages. Intriguingly, intron structure evolution appears to be quite different across lineages experiencing different types of genome reduction: whereas lineages with very few introns tend towards highly regular intronic sequences, lineages with very short introns tend towards highly degenerate sequences. Together, these results attest to the complex nature of ancestral eukaryotic splicing, the qualitatively different evolutionary forces acting on intron structures across modern lineages, and the impressive evolutionary malleability of eukaryotic gene structures.     

Molecular evolution and phylogenetic utility of the petD group II intron: a case study in basal angiosperms

Sequences of spacers and group I introns in plant chloroplast genomes have recently been shown to be very effective in phylogenetic reconstruction at higher taxonomic levels and not only for inferring relationships among species. Group II introns, being more frequent in those genomes than group I introns, may be further promising markers. Because group II introns are structurally constrained, we assumed that sequences of a group II intron should be alignable across seed plants. We designed universal amplification primers for the petD intron and sequenced this intron in a representative selection of 47 angiosperms and three gymnosperms. Our sampling of taxa is the most representative of major seed plant lineages to date for group II introns. Through differential analysis of structural partitions, we studied patterns of molecular evolution and their contribution to phylogenetic signal. Nonpairing stretches (loops, bulges, and interhelical nucleotides) were considerably more variable in both substitutions and indels than in helical elements. Differences among the domains are basically a function of their structural composition. After the exclusion of four mutational hotspots accounting for less than 18% of sequence length, which are located in loops of domains I and IV, all sequences could be aligned unambiguously across seed plants. Microstructural changes predominantly occurred in loop regions and are mostly simple sequence repeats. An indel matrix comprising 241 characters revealed microstructural changes to be of lower homoplasy than are substitutions. In showing Amborella first branching and providing support for a magnoliid clade through a synapomorphic indel, the petD data set proved effective in testing between alternative hypotheses on the basal nodes of the angiosperm tree. Within angiosperms, group II introns offer phylogenetic signal that is intermediate in information content between that of spacers and group I introns on the one hand and coding sequences on the other.     

Evolution and biogeography of Alectryon (Sapindaceae).

Phylogenetic analysis of nucleotide sequences from four plastid loci (matK, partial trnK-matK introns, rps16 intron) and one nuclear locus (the internal transcribed spacer of rDNA; ITS-1) was conducted for 14 species of Alectryon and five related genera in Sapindaceae. Both matK and rps16 intron provide few informative characters within Alectryon, whereas ITS-1 provides the largest number of parsimony-informative characters and has the greatest sequence divergence between taxa. Support for branches in cladograms produced in PAUP increased markedly upon inclusion of ITS-1 data to matK and rps16 intron data. Analyses of each region alone or combined produced congruent results, suggesting that the regions are complementary. Phylogenetic analysis indicates that there are two main lineages within Alectryon, with A. subcinereus sister to the remaining sampled Alectryon taxa. Two morphological characters, presence/absence of petals and aril patterning, are congruent with the molecular phylogeny. One robustly supported clade is characterized by smooth arils and petals, in contrast to the taxa in the other major clade which have patterned arils and an absence of petals. These analyses also support a number of revised subgeneric groupings for Alectryon. The decision to submerge Heterodendrum in Alectryon is supported, although taxa belonging to Heterodendrum do not form a clade. The majority of the Australian Alectryon appear to belong to the tropical monsoonal/arid flora with species from both lineages being found in representative vine thickets across northern Australia. It appears that the seasonally dry rainforest communities comprise a number of elements that do not share common evolutionary histories within this genus.     

The phylogenetic analysis of variable-length sequence data: elongation factor-1alpha introns in European populations of the parasitoid wasp genus Pauesia (Hymenoptera: Braconidae: Aphidiinae)

Elongation factor-1alpha (EF-1alpha) is a highly conserved nuclear coding gene that can be used to investigate recent divergences due to the presence of rapidly evolving introns. However, a universal feature of intron sequences is that even closely related species exhibit insertion and deletion events, which cause variation in the lengths of the sequences. Indels are frequently rich in evolutionary information, but most investigators ignore sites that fall within these variable regions, largely because the analytical tools and theory are not well developed. We examined this problem in the taxonomically problematic parasitoid wasp genus Pauesia (Hymenoptera: Braconidae: Aphidiinae) using congruence as a criterion for assessing a range of methods for aligning such variable-length EF-1alpha intron sequences. These methods included distance- and parsimony-based multiple-alignment programs (CLUSTAL W and MALIGN), direct optimization (POY), and two "by eye" alignment strategies. Furthermore, with one method (CLUSTAL W) we explored in detail the robustness of results to changes in the gap cost parameters. Phenetic-based alignments ("by eye" and CLUSTAL W) appeared, under our criterion, to perform as well as more readily defensible, but computationally more demanding, methods. In general, all of our alignment and tree-building strategies recovered the same basic topological structure, which means that an underlying phylogenetic signal remained regardless of the strategy chosen. However, several relationships between clades were sensitive both to alignment and to tree-building protocol. Further alignments, considering only sequences belonging to the same group, allowed us to infer a range of phylogenetic relationships that were highly robust to tree-building protocol. By comparing these topologies with those obtained by varying the CLUSTAL parameters, we generated the distribution area of congruence and taxonomic compatibility. Finally, we present the first robust estimate of the European Pauesia phylogeny by using two EF-1alpha introns and 38 taxa (plus 3 outgroups). This estimate conflicts markedly with the traditional subgeneric classification. We recommend that this classification be abandoned, and we propose a series of monophyletic species groups.     

Stability of characters and construction of phylogenetic trees.

Parsimony methods infer phylogenetic trees by minimizing number of character changes required to explain observed character states. From the perspective of applicability of parsimony methods, it is important to assess whether the characters used to infer phylogeny are likely to provide a correct tree. We introduce a graph theoretical characterization that helps to assess whether given set of characters is appropriate to use with parsimony methods. Given a set of characters and a set of taxa, we construct a network called character overlap graph. We show that the character overlap graph for characters that are appropriate to use in parsimony methods is characterized by significant under-representation of subnetworks known as holes, and provide a validation for this observation. This characterization explains success in constructing evolutionary trees using parsimony method for some characters (e.g., protein domains) and lack of such success for other characters (e.g., introns). In the latter case, the understanding of obstacles to applying parsimony methods in a direct way has lead us to a new approach for detecting inconsistent and/or noisy data. Namely, we introduce the concept of stable characters which is similar but less restrictive than the well known concept of pairwise compatible characters. Application of this approach to introns produces the evolutionary tree consistent with the Coelomata hypothesis.     

Group II introns as phylogenetic tools: structure, function, and evolutionary constraints

Group II introns comprise the majority of noncoding DNA in many plant chloroplast genomes and include the commonly sequenced regions trnK/matK, the rps16 intron, and the rpl16 intron. As demand increases for nucleotide characters at lower taxonomic levels, chloroplast introns may come to provide the bulk of plastome sequence data for assessment of evolutionary relationships in infrageneric, intergeneric, and interfamilial studies. Group II introns have many attractive properties for the molecular systematist: they are confined to organellar genomes in eukaryotes and the majority are single-copy; they share a well-defined and empirically tested secondary and tertiary structure; and many are easily amplified due to highly conserved sequence in flanking exons. However, structure-linked mutation patterns in group II intron sequences are more complex than generally supposed and have important implications for aligning nucleotides, assessing mutational biases in the data, and selecting appropriate models of character evolution for phylogenetic analysis. This paper presents a summary of group II intron function and structure, reviews the link between that structure and specific mutational constraints in group II intron sequences, and discusses strategies for accommodating the resulting complex mutational patterns in subsequent phylogenetic analyses.