Phylogenetic and evolutionary implications of complete chloroplast genome sequences of four early-diverging angiosperms: Buxus (Buxaceae), Chloranthus (Chloranthaceae), Dioscorea (Dioscoreaceae), and Illicium (Schisandraceae)

We have determined the complete chloroplast genome sequences of four early-diverging lineages of angiosperms, Buxus (Buxaceae), Chloranthus (Chloranthaceae), Dioscorea (Dioscoreaceae), and Illicium (Schisandraceae), to examine the organization and evolution of plastid genomes and to estimate phylogenetic relationships among angiosperms. For the most part, the organization of these plastid genomes is quite similar to the ancestral angiosperm plastid genome with a few notable exceptions. Dioscorea has lost one protein-coding gene, rps16; this gene loss has also happened independently in four other land plant lineages, liverworts, conifers, Populus, and legumes. There has also been a small expansion of the inverted repeat (IR) in Dioscorea that has duplicated trnH-GUG. This event has also occurred multiple times in angiosperms, including in monocots, and in the two basal angiosperms Nuphar and Drimys. The Illicium chloroplast genome is unusual by having a 10 kb contraction of the IR. The four taxa sequenced represent key groups in resolving phylogenetic relationships among angiosperms. Illicium is one of the basal angiosperms in the Austrobaileyales, Chloranthus (Chloranthales) remains unplaced in angiosperm classifications, and Buxus and Dioscorea are early-diverging eudicots and monocots, respectively. We have used sequences for 61 shared protein-coding genes from these four genomes and combined them with sequences from 35 other genomes to estimate phylogenetic relationships using parsimony, likelihood, and Bayesian methods. There is strong congruence among the trees generated by the three methods, and most nodes have high levels of support. The results indicate that Amborella alone is sister to the remaining angiosperms; the Nymphaeales represent the next-diverging clade followed by Illicium; Chloranthus is sister to the magnoliids and together this group is sister to a large clade that includes eudicots and monocots; and Dioscorea represents an early-diverging lineage of monocots just internal to Acorus.     

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 chicken or the egg? Plastome evolution and an independent loss of the inverted repeat in papilionoid legumes

The plastid genome (plastome), while surprisingly constant in gene order and content across most photosynthetic angiosperms, exhibits variability in several unrelated lineages. During the diversification history of the legume family Fabaceae, plastomes have undergone many rearrangements, including inversions, expansion, contraction and loss of the typical inverted repeat (IR), gene loss and repeat accumulation in both shared and independent events. While legume plastomes have been the subject of study for some time, most work has focused on agricultural species in the IR-lacking clade (IRLC) and the plant model Medicago truncatula. The subfamily Papilionoideae, which contains virtually all of the agricultural legume species, also comprises most of the plastome variation detected thus far in the family. In this study three non-papilioniods were included among 34 newly sequenced legume plastomes, along with 33 publicly available sequences, to assess plastome structural evolution in the subfamily. In an effort to examine plastome variation across the subfamily, approximately 20% of the sampling represents the IRLC with the remainder selected to represent the early-branching papilionoid clades. A number of IR-related and repeat-mediated changes were identified and examined in a phylogenetic context. Recombination between direct repeats associated with ycf2 resulted in intraindividual plastome heteroplasmy. Although loss of the IR has not been reported in legumes outside of the IRLC, one genistoid taxon was found to completely lack the typical plastome IR. The role of the IR and non-IR repeats in the progression of plastome change is discussed.     

Extensive Reorganization of the Plastid Genome of <Emphasis Type="Italic">Trifolium subterraneum</Emphasis> (Fabaceae) Is Associated with Numerous Repeated Sequences and Novel DNA Insertions

The plastid genome of Trifolium subterraneum is 144,763 bp, about 20 kb longer than those of closely related legumes, which also lost one copy of the large inverted repeat (IR). The genome has undergone extensive genomic reconfiguration, including the loss of six genes (accD, infA, rpl22, rps16, rps18, and ycf1) and two introns (clpP and rps12) and numerous gene order changes, attributable to 14–18 inversions. All endpoints of rearranged gene clusters are flanked by repeated sequences, tRNAs, or pseudogenes. One unusual feature of the Trifolium subterraneum genome is the large number of dispersed repeats, which comprise 19.5% (ca. 28 kb) of the genome (versus about 4% for other angiosperms) and account for part of the increase in genome size. Nine genes (psbT, rbcL, clpP, rps3, rpl23, atpB, psbN, trnI-cau, and ycf3) have also been duplicated either partially or completely. rpl23 is the most highly duplicated gene, with portions of this gene duplicated six times. Comparisons of the Trifolium plastid genome with the Plant Repeat Database and searches for flanking inverted repeats suggest that the high incidence of dispersed repeats and rearrangements is not likely the result of transposition. Trifolium has 19.5 kb of unique DNA distributed among 160 fragments ranging in size from 30 to 494 bp, greatly surpassing the other five sequenced legume plastid genomes in novel DNA content. At least some of this unique DNA may represent horizontal transfer from bacterial genomes. These unusual features provide direction for the development of more complex models of plastid genome evolution. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Plastid genome sequences of legumes reveal parallel inversions and multiple losses of rps16 in papilionoids

To date, publicly available plastid genomes of legumes have for the most part been limited to the subfamily Papilionoideae. Here we report 13 new plastid genomes of legumes spanning all three subfamilies. The genomes representing Caesalpinioideae and Mimosoideae are highly conserved in gene content and gene order, similar to the ancestral angiosperm genome organization. Genomes within the Papilionoideae, however, have reduced sizes due to deletions in nine intergenic spacers primarily in the large single copy region. Our study also indicates that rps16 has been independently lost at least five times in legumes, with additional gene and intron losses scattered among the papilionoids. Additionally, genera from two distinct lineages within the papilionoids, Lupinus and Robinia, have a parallel inversion of 36 and 39 kb, respectively. This parallel inversion is novel as it appears to be caused by a 29 bp repeat within two trnS genes. This repeat is present in all available legume plastid genomes indicating that there is the potential for this inversion to be present in more species. This case of a homoplasious inversion is also evidence that some inversion events may not be reliable phylogenetic markers.     

Phylogenomics and biogeography of Wisteria: Implications on plastome evolution among inverted repeat-lacking clade (IRLC) legumes

The genus Wisteria (Fabaceae) is disjunctly distributed in eastern Asian and eastern North American temperate deciduous forests, and it is widely cultivated around the world as spectacular garden plants. It is a member of inverted repeat-lacking clade (IRLC). The IRLC Species are characterized by the loss of an IR region in their plastomes, which has long been of great interest. In this research, we report whole plastome sequences from all four Wisteria species and a Wisteriopsis japonica, combining these with existing data to explore phylogenetic relationships and biogeography of Wisteria, as well as plastome evolution of IRLC species. Phylogenetic analyses recognized a clade containing Glycyrrhiza–Wisteriopsis–Wisteria as sister to the remaining genera of IRLC. North American Wisteria frutescens and the three Asian species formed reciprocal clades, and Wisteria brachybotrys was sister to Wisteria floribunda and Wisteria sinensis. Wisteria may have originated in Japan near the boundary of the Oligocene and Miocene. The disappearance of Bering Land Bridge in the late Miocene might lead to the Eastern Asian–Eastern North American disjunction of Wisteria. Allopatric speciation of Wisteria between the Japanese archipelago and the Asian continent in the Quaternary increased the species richness of eastern Asia in comparison with eastern North America. Synonymous substitution rates (d_S) of protein-coding genes in the IRLC species were around 2-fold (SC genes) or 11-fold (IR genes) higher than those of non-IRLC species. For both SC and IR genes, herbaceous legumes have around 3-fold higher d_S than woody ones. Both loss of one IR region and herbaceous habit elevated substitution rates of the plastomes.     

The complete nucleotide sequence of the cassava (Manihot esculenta) chloroplast genome and the evolution of atpF in Malpighiales: RNA editing and multiple losses of a group II intron

The complete sequence of the chloroplast genome of cassava (Manihot esculenta, Euphorbiaceae) has been determined. The genome is 161,453 bp in length and includes a pair of inverted repeats (IR) of 26,954 bp. The genome includes 128 genes; 96 are single copy and 16 are duplicated in the IR. There are four rRNA genes and 30 distinct tRNAs, seven of which are duplicated in the IR. The infA gene is absent; expansion of IRb has duplicated 62 amino acids at the 3′ end of rps19 and a number of coding regions have large insertions or deletions, including insertions within the 23S rRNA gene. There are 17 intron-containing genes in cassava, 15 of which have a single intron while two (clpP, ycf3) have two introns. The usually conserved atpF group II intron is absent and this is the first report of its loss from land plant chloroplast genomes. The phylogenetic distribution of the atpF intron loss was determined by a PCR survey of 251 taxa representing 34 families of Malpighiales and 16 taxa from closely related rosids. The atpF intron is not only missing in cassava but also from closely related Euphorbiaceae and other Malpighiales, suggesting that there have been at least seven independent losses. In cassava and all other sequenced Malphigiales, atpF gene sequences showed a strong association between C-to-T substitutions at nucleotide position 92 and the loss of the intron, suggesting that recombination between an edited mRNA and the atpF gene may be a possible mechanism for the intron loss.     

Gene relocations within chloroplast genomes of Jasminum and Menodora (Oleaceae) are due to multiple, overlapping inversions

The chloroplast (cp) DNA sequence of Jasminum nudiflorum (Oleaceae-Jasmineae) is completed and compared with the large single-copy region sequences from 6 related species. The cp genomes of the tribe Jasmineae (Jasminum and Menodora) show several distinctive rearrangements, including inversions, gene duplications, insertions, inverted repeat expansions, and gene and intron losses. The ycf4-psaI region in Jasminum section Primulina was relocated as a result of 2 overlapping inversions of 21,169 and 18,414 bp. The 1st, larger inversion is shared by all members of the Jasmineae indicating that it occurred in the common ancestor of the tribe. Similar rearrangements were also identified in the cp genome of Menodora. In this case, 2 fragments including ycf4 and rps4-trnS-ycf3 genes were moved by 2 additional inversions of 14 and 59 kb that are unique to Menodora. Other rearrangements in the Oleaceae are confined to certain regions of the Jasminum and Menodora cp genomes, including the presence of highly repeated sequences and duplications of coding and noncoding sequences that are inserted into clpP and between rbcL and psaI. These insertions are correlated with the loss of 2 introns in clpP and a serial loss of segments of accD. The loss of the accD gene and clpP introns in both the monocot family Poaceae and the eudicot family Oleaceae are clearly independent evolutionary events. However, their genome organization is surprisingly similar despite the distant relationship of these 2 angiosperm families.     

A phylogeny of legumes (Leguminosae) based on analysis of the plastid matK gene resolves many well-supported subclades within the family

Phylogenetic analysis of 330 plastid matK gene sequences, representing 235 genera from 37 of 39 tribes, and four outgroup taxa from eurosids I supports many well-resolved subclades within the Leguminosae. These results are generally consistent with those derived from other plastid sequence data (rbcL and trnL), but show greater resolution and clade support overall. In particular, the monophyly of subfamily Papilionoideae and at least seven major subclades are well-supported by bootstrap and Bayesian credibility values. These subclades are informally recognized as the Cladrastis clade, genistoid sensu lato, dalbergioid sensu lato, mirbelioid, millettioid, and robinioid clades, and the inverted-repeat-lacking clade (IRLC). The genistoid clade is expanded to include genera such as Poecilanthe, Cyclolobium, Bowdichia, and Diplotropis and thus contains the vast majority of papilionoids known to produce quinolizidine alkaloids. The dalbergioid clade is expanded to include the tribe Amorpheae. The mirbelioids include the tribes Bossiaeeae and Mirbelieae, with Hypocalypteae as its sister group. The millettioids comprise two major subclades that roughly correspond to the tribes Millettieae and Phaseoleae and represent the only major papilionoid clade marked by a macromorphological apomorphy, pseudoracemose inflorescences. The robinioids are expanded to include Sesbania and members of the tribe Loteae. The IRLC, the most species-rich subclade, is sister to the robinioids. Analysis of the matK data consistently resolves but modestly supports a clade comprising papilionoid taxa that accumulate canavanine in the seeds. This suggests a single origin for the biosynthesis of this most commonly produced of the nonprotein amino acids in legumes.     

Lineage-specific group II intron gains and losses of the mitochondrial rps3 gene in gymnosperms

According to PCR assays and sequencing, we now report the shared presence of two rps3 introns, namely the rps3i74 and the rps3i249, in the mitochondria of all the classes representing the surviving lineages of gymnosperms, and unveil several lineages experiencing intron loss.Interestingly, the rps3 intron gains and losses within the four groups of gymnosperms let us sort out the Pinaceae and the non-Pinaceae into intron (+)- and intron (?)-lineages, respectively. Worthy of mention is also the finding that only Gnetum within the Gnetales harbours both the rps3 introns.This intron distribution pattern is consistent with the hypothesis that the two rps3 introns were likely present in the common ancestor of the seed plants and, then, independently lost in the non-Pinaceae during gymnosperm evolution.The derived secondary structural model of the novel group IIA intron improves our understanding of the significance and origin of the extraordinary length polymorphisms observed among rps3i249 orthologs.Despite the remarkable structural plasticity to adopt and reject introns, the rps3 mRNAs undergo accurate processing by splicing and extensive editing in gymnosperm mitochondria.This study provides additional insights into the evolutionarily high dynamics of mitochondrial introns which may come and go in closely related plant species.The turnover of the mitochondrial rps3 group II introns seen among lineages of seed plants further suggests that these introns might be an additional signature to discriminate between particularly cryptical taxonomic groups for which there is a need of a further evaluation of their evolutionary affiliation.     

Nonessential plastid-encoded ribosomal proteins in tobacco: a developmental role for plastid translation and implications for reductive genome evolution

Plastid genomes of higher plants contain a conserved set of ribosomal protein genes. Although plastid translational activity is essential for cell survival in tobacco (Nicotiana tabacum), individual plastid ribosomal proteins can be nonessential. Candidates for nonessential plastid ribosomal proteins are ribosomal proteins identified as nonessential in bacteria and those whose genes were lost from the highly reduced plastid genomes of nonphotosynthetic plastid-bearing lineages (parasitic plants, apicomplexan protozoa). Here we report the reverse genetic analysis of seven plastid-encoded ribosomal proteins that meet these criteria. We have introduced knockout alleles for the corresponding genes into the tobacco plastid genome. Five of the targeted genes (ribosomal protein of the large subunit22 [rpl22], rpl23, rpl32, ribosomal protein of the small subunit3 [rps3], and rps16) were shown to be essential even under heterotrophic conditions, despite their loss in at least some parasitic plastid-bearing lineages. This suggests that nonphotosynthetic plastids show elevated rates of gene transfer to the nuclear genome. Knockout of two ribosomal protein genes, rps15 and rpl36, yielded homoplasmic transplastomic mutants, thus indicating nonessentiality. Whereas Δrps15 plants showed only a mild phenotype, Δrpl36 plants were severely impaired in photosynthesis and growth and, moreover, displayed greatly altered leaf morphology. This finding provides strong genetic evidence that chloroplast translational activity influences leaf development, presumably via a retrograde signaling pathway.     

Plastid phylogenomics and molecular evolution of Thismiaceae (Dioscoreales)

Premise

Species in Thismiaceae can no longer photosynthesize and instead obtain carbon from soil fungi. Here we infer Thismiaceae phylogeny using plastid genome data and characterize the molecular evolution of this genome.

Methods

We assembled five Thismiaceae plastid genomes from genome skimming data, adding to previously published data for phylogenomic inference. We investigated plastid-genome structural changes, considering locally colinear blocks (LCBs). We also characterized possible shifts in selection pressure in retained genes by considering changes in the ratio of nonsynonymous to synonymous changes (ω).

Results

Thismiaceae experienced two major pulses of gene loss around the early diversification of the family, with subsequent scattered gene losses across descendent lineages. In addition to massive size reduction, Thismiaceae plastid genomes experienced occasional inversions, and there were likely two independent losses of the plastid inverted repeat (IR) region. Retained plastid genes remain under generally strong purifying selection (ω << 1), with significant and sporadic weakening or strengthening in several instances. The bifunctional trnE-UUC gene of Thismia huangii may retain a secondary role in heme biosynthesis, despite a probable loss of functionality in protein translation. Several cis-spliced group IIA introns have been retained, despite the loss of the plastid intron maturase, matK.

Conclusions

We infer that most gene losses in Thismiaceae occurred early and rapidly, following the initial loss of photosynthesis in its stem lineage. As a species-rich, fully mycoheterotrophic lineage, Thismiaceae provide a model system for uncovering the unique and divergent ways in which plastid genomes evolve in heterotrophic plants.     

Evidence for Transitional Stages in the Evolution of Euglenid Group II Introns and Twintrons in the Monomorphina aenigmatica Plastid Genome

Background

Photosynthetic euglenids acquired their plastid by secondary endosymbiosis of a prasinophyte-like green alga. But unlike its prasinophyte counterparts, the plastid genome of the euglenid Euglena gracilis is riddled with introns that interrupt almost every protein-encoding gene. The atypical group II introns and twintrons (introns-within-introns) found in the E. gracilis plastid have been hypothesized to have been acquired late in the evolution of euglenids, implying that massive numbers of introns may be lacking in other taxa. This late emergence was recently corroborated by the plastid genome sequences of the two basal euglenids, Eutreptiella gymnastica and Eutreptia viridis, which were found to contain fewer introns.

Methodology/Principal Findings

To gain further insights into the proliferation of introns in euglenid plastids, we have characterized the complete plastid genome sequence of Monomorphina aenigmatica, a freshwater species occupying an intermediate phylogenetic position between early and late branching euglenids. The M. aenigmatica UTEX 1284 plastid genome (74,746 bp, 70.6% A+T, 87 genes) contains 53 intron insertion sites, of which 41 were found to be shared with other euglenids including 12 of the 15 twintron insertion sites reported in E. gracilis.

Conclusions

The pattern of insertion sites suggests an ongoing but uneven process of intron gain in the lineage, with perhaps a minimum of two bursts of rapid intron proliferation. We also identified several sites that represent intermediates in the process of twintron evolution, where the external intron is in place, but not the internal one, offering a glimpse into how these convoluted molecular contraptions originate.     

Reticulate phylogenetics and phytogeographical structure of Heliosperma (Sileneae, Caryophyllaceae) inferred from chloroplast and nuclear DNA sequences

The Balkan Peninsula is known to be one of the most diverse and species-rich parts of Europe, but its biota has gained much less attention in phylogenetic and evolutionary studies compared to other southern European mountain systems. We used nuclear ribosomal internal transcribed spacer (ITS) sequences and intron sequences of the chloroplast gene rps16 to examine phylogenetic and biogeographical patterns within the genus Heliosperma (Sileneae, Caryophyllaceae). The ITS and rps16 intron sequences both support monophyly of Heliosperma, but the data are not conclusive with regard to its exact origin. Three strongly supported clades are found in both data sets, corresponding to Heliosperma alpestre, Heliosperma macranthum and the Heliosperma pusillum clade, including all other taxa. The interrelationships among these three differ between the nuclear and the plastid data sets. Hierarchical relationships within the H. pusillum clade are poorly resolved by the ITS data, but the rps16 intron sequences form two well-supported clades which are geographically, rather than taxonomically, correlated. A similar geographical structure is found in the ITS data, when analyzed with the NeighbourNet method. The apparent rate of change within Heliosperma is slightly higher for rps16 as compared to ITS. In contrast, in the Sileneae outgroup, ITS substitution rates are more than twice as high as those for rps16, a situation more in agreement with what has been found in other rate comparisons of noncoding cpDNA and ITS. Unlike most other Sileneae ITS sequences, the H. pusillum group sequences display extensive polymorphism. A possible explanation to these patterns is extensive hybridization and gene flow within Heliosperma, which together with concerted evolution may have eradicated the ancient divergence suggested by the rps16 data. The morphological differentiation into high elevation, mainly widely distributed taxa, and low elevation narrow endemics is not correlated with the molecular data, and is possibly a result of ecological differentiation.