Phylogenetics of early branching eudicots:Comparing phylogenetic signal across plastid introns,spacers,and genes

<正>Recent phylogenetic analyses revealed a grade with Ranunculales,Sabiales,Proteales,Trochodendrales, and Buxales as first branching eudicots,with the respective positions of Proteales and Sabiales still lacking statistical confidence.As previous analyses of conserved plastid genes remain inconclusive,we aimed to use and evaluate a representative set of plastid introns(group I:trnL;group II:petD,rpll6,trnK) and intergenic spacers(trnL-F,petB-petD, atpB-rbcL,rps3-rpl16) in comparison to the rapidly evolving matK and slowly evolving atpB and rbcL genes. Overall patterns of micro structural mutations converged across genomic regions,underscoring the existence of a general mutational pattern throughout the plastid genome.Phylogenetic signal differed strongly between functionally and structurally different genomic regions and was highest in matK,followed by spacers,then group II and group I introns.The more conserved atpB and rbcL coding regions showed distinctly lower phylogenetic information content.Parsimony,maximum likelihood,and Bayesian phylogenetic analyses based on the combined dataset of non-coding and rapidly evolving regions(14 000 aligned characters) converged to a backbone topology of eudicots with Ranunculales branching first,a Proteales-Sabiales clade second,followed by Trochodendrales and Buxales. Gunnerales generally appeared as sister to all remaining core eudicots with maximum support.Our results show that a small number of intron and spacer sequences allow similar insights into phylogenetic relationships of eudicots compared to datasets of many combined genes.The non-coding proportion of the plastid genome thus can be considered an important information source for plastid phylogenomics.     

Highly variable chloroplast markers for evaluating plant phylogeny at low taxonomic levels and for DNA barcoding

Background

At present, plant molecular systematics and DNA barcoding techniques rely heavily on the use of chloroplast gene sequences. Because of the relatively low evolutionary rates of chloroplast genes, there are very few choices suitable for molecular studies on angiosperms at low taxonomic levels, and for DNA barcoding of species.

Methodology/Principal Findings

We scanned the entire chloroplast genomes of 12 genera to search for highly variable regions. The sequence data of 9 genera were from GenBank and 3 genera were of our own. We identified nearly 5% of the most variable loci from all variable loci in the chloroplast genomes of each genus, and then selected 23 loci that were present in at least three genera. The 23 loci included 4 coding regions, 2 introns, and 17 intergenic spacers. Of the 23 loci, the most variable (in order from highest variability to lowest) were intergenic regions ycf1-a, trnK, rpl32-trnL, and trnH-psbA, followed by trnS^UGA-trnG^UCC, petA-psbJ, rps16-trnQ, ndhC-trnV, ycf1-b, ndhF, rpoB-trnC, psbE-petL, and rbcL-accD. Three loci, trnS^UGA-trnG^UCC, trnT-psbD, and trnW-psaJ, showed very high nucleotide diversity per site (π values) across three genera. Other loci may have strong potential for resolving phylogenetic and species identification problems at the species level. The loci accD-psaI, rbcL-accD, rpl32-trnL, rps16-trnQ, and ycf1 are absent from some genera. To amplify and sequence the highly variable loci identified in this study, we designed primers from their conserved flanking regions. We tested the applicability of the primers to amplify target sequences in eight species representing basal angiosperms, monocots, eudicots, rosids, and asterids, and confirmed that the primers amplified the desired sequences of these species.

Significance/Conclusions

Chloroplast genome sequences contain regions that are highly variable. Such regions are the first consideration when screening the suitable loci to resolve closely related species or genera in phylogenetic analyses, and for DNA barcoding.     

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.     

Comparative analysis of inverted repeats of polypod fern (Polypodiales) plastomes reveals two hypervariable regions

Background

Ferns are large and underexplored group of vascular plants (~ 11 thousands species). The genomic data available by now include low coverage nuclear genomes sequences and partial sequences of mitochondrial genomes for six species and several plastid genomes.

Results

We characterized plastid genomes of three species of Dryopteris, which is one of the largest fern genera, using sequencing of chloroplast DNA enriched samples and performed comparative analysis with available plastomes of Polypodiales, the most species-rich group of ferns. We also sequenced the plastome of Adianthum hispidulum (Pteridaceae). Unexpectedly, we found high variability in the IR region, including duplication of rrn16 in D. blanfordii, complete loss of trnI-GAU in D. filix-mas, its pseudogenization due to the loss of an exon in D. blanfordii. Analysis of previously reported plastomes of Polypodiales demonstrated that Woodwardia unigemmata and Lepisorus clathratus have unusual insertions in the IR region. The sequence of these inserted regions has high similarity to several LSC fragments of ferns outside of Polypodiales and to spacer between tRNA-CGA and tRNA-TTT genes of mitochondrial genome of Asplenium nidus. We suggest that this reflects the ancient DNA transfer from mitochondrial to plastid genome occurred in a common ancestor of ferns. We determined the marked conservation of gene content and relative evolution rate of genes and intergenic spacers in the IRs of Polypodiales. Faster evolution of the four intergenic regions had been demonstrated (trnA- orf42, rrn16-rps12, rps7-psbA and ycf2-trnN).

Conclusions

IRs of Polypodiales plastomes are dynamic, driven by such events as gene loss, duplication and putative lateral transfer from mitochondria.

     

Tobacco plastid ribosomal protein S18 is essential for cell survival

Plastid genomes contain a conserved set of genes most of which are involved in either photosynthesis or gene expression. Among the ribosomal protein genes present in higher plant plastid genomes, rps18 is special in that it is absent from the plastid genomes of several non-green unicellular organisms, including Euglena longa and Toxoplasma gondii. Here we have tested whether the ribosomal protein S18 is required for translation by deleting the rps18 gene from the tobacco plastid genome. We report that, while deletion of the rps18 gene was readily obtained, no homoplasmic Δrps18 plants or leaf sectors could be isolated. Instead, segregation into homoplasmy led to severe defects in leaf development suggesting that the knockout of rps18 is lethal and the S18 protein is required for cell survival. Our data demonstrate that S18 is indispensable for plastid ribosome function in tobacco and support an essential role for plastid translation in plant development. Moreover, we demonstrate the occurrence of flip-flop recombination on short inverted repeat sequences which generates different isoforms of the transformed plastid genome that differ in the orientation a 70 kb segment in the large single-copy region. However, infrequent occurrence of flip-flop recombination and random segregation of plastid genomes result in the predominant presence of only one of the isoforms in many tissue samples. Implications for the interpretation of chloroplast transformation experiments and vector design are discussed.     

Selection on Accessible Chromatin Regions in Capsella grandiflora

Accurate estimates of genome-wide rates and fitness effects of new mutations are essential for an improved understanding of molecular evolutionary processes. Although eukaryotic genomes generally contain a large noncoding fraction, functional noncoding regions and fitness effects of mutations in such regions are still incompletely characterized. A promising approach to characterize functional noncoding regions relies on identifying accessible chromatin regions (ACRs) tightly associated with regulatory DNA. Here, we applied this approach to identify and estimate selection on ACRs in Capsella grandiflora, a crucifer species ideal for population genomic quantification of selection due to its favorable population demography. We describe a population-wide ACR distribution based on ATAC-seq data for leaf samples of 16 individuals from a natural population. We use population genomic methods to estimate fitness effects and proportions of positively selected fixations (α) in ACRs and find that intergenic ACRs harbor a considerable fraction of weakly deleterious new mutations, as well as a significantly higher proportion of strongly deleterious mutations than comparable inaccessible intergenic regions. ACRs are enriched for expression quantitative trait loci (eQTL) and depleted of transposable element insertions, as expected if intergenic ACRs are under selection because they harbor regulatory regions. By integrating empirical identification of intergenic ACRs with analyses of eQTL and population genomic analyses of selection, we demonstrate that intergenic regulatory regions are an important source of nearly neutral mutations. These results improve our understanding of selection on noncoding regions and the role of nearly neutral mutations for evolutionary processes in outcrossing Brassicaceae species.     

Complete Plastid Genome Sequencing of Trochodendraceae Reveals a Significant Expansion of the Inverted Repeat and Suggests a Paleogene Divergence between the Two Extant Species

The early-diverging eudicot order Trochodendrales contains only two monospecific genera, Tetracentron and Trochodendron. Although an extensive fossil record indicates that the clade is perhaps 100 million years old and was widespread throughout the Northern Hemisphere during the Paleogene and Neogene, the two extant genera are both narrowly distributed in eastern Asia. Recent phylogenetic analyses strongly support a clade of Trochodendrales, Buxales, and Gunneridae (core eudicots), but complete plastome analyses do not resolve the relationships among these groups with strong support. However, plastid phylogenomic analyses have not included data for Tetracentron. To better resolve basal eudicot relationships and to clarify when the two extant genera of Trochodendrales diverged, we sequenced the complete plastid genome of Tetracentron sinense using Illumina technology. The Tetracentron and Trochodendron plastomes possess the typical gene content and arrangement that characterize most angiosperm plastid genomes, but both genomes have the same unusual ∼4 kb expansion of the inverted repeat region to include five genes (rpl22, rps3, rpl16, rpl14, and rps8) that are normally found in the large single-copy region. Maximum likelihood analyses of an 83-gene, 88 taxon angiosperm data set yield an identical tree topology as previous plastid-based trees, and moderately support the sister relationship between Buxaceae and Gunneridae. Molecular dating analyses suggest that Tetracentron and Trochodendron diverged between 44-30 million years ago, which is congruent with the fossil record of Trochodendrales and with previous estimates of the divergence time of these two taxa. We also characterize 154 simple sequence repeat loci from the Tetracentron sinense and Trochodendron aralioides plastomes that will be useful in future studies of population genetic structure for these relict species, both of which are of conservation concern.     

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.     

Synthetic Lethality in the Tobacco Plastid Ribosome and Its Rescue at Elevated Growth Temperatures

Consistent with their origin from cyanobacteria, plastids (chloroplasts) perform protein biosynthesis on bacterial-type 70S ribosomes. The plastid genomes of seed plants contain a conserved set of ribosomal protein genes. Three of these have proven to be nonessential for translation and, thus, for cellular viability: rps15, rpl33, and rpl36. To help define the minimum ribosome, here, we examined whether more than one of these nonessential plastid ribosomal proteins can be removed from the 70S ribosome. To that end, we constructed all possible double knockouts for the S15, L33, and L36 ribosomal proteins by stable transformation of the tobacco (Nicotiana tabacum) plastid genome. We find that, although S15 and L33 function in different ribosomal particles (30S and 50S, respectively), their combined deletion from the plastid genome results in synthetic lethality under autotrophic conditions. Interestingly, the lethality can be overcome by growth under elevated temperatures due to an improved efficiency of plastid ribosome biogenesis. Our results reveal functional interactions between protein and RNA components of the 70S ribosome and uncover the interdependence of the biogenesis of the two ribosomal subunits. In addition, our findings suggest that defining a minimal set of plastid genes may prove more complex than generally believed.     

Evidence for replication slippage in the evolution of Oenothera chloroplast DNA.

Evolutionary relationships of four plastid genomes (plastomes) from different Oenothera species have been assessed by sequence comparisons of two intergenic regions that separate the ribosomal protein genes rpl16, rpl14, and rps8. Sequence changes include base substitutions, the occurrence of a 29-base tandem duplication, and variation in the length of two poly-A stretches. Additions/deletions in chloroplast DNA may not be useful for evolutionary comparisons more distant than these, particularly if the sequences undergo divergence after the initial event, but the length mutations reported here allow a finer resolution of the phylogeny of the closely related Oenothera plastomes than would have been possible if only base substitutions had been considered. Comparisons with the orthogous sequence from tobacco chloroplast DNA indicate the direction of change at most of the sites. The results suggest that plastomes I and II are closely related to each other, as are plastomes III and IV. Replication slippage is proposed as a mechanism to explain the length mutations.     

Long-run bacteria-phage coexistence dynamics under natural habitat conditions in an environmental biotechnology system

Bacterial viruses are widespread and abundant across natural and engineered habitats. They influence ecosystem functioning through interactions with their hosts. Laboratory studies of phage–host pairs have advanced our understanding of phenotypic and genetic diversification in bacteria and phages. However, the dynamics of phage–host interactions have been seldom recorded in complex natural environments. We conducted an observational metagenomic study of the dynamics of interaction between Gordonia and their phages using a three-year data series of samples collected from a full-scale wastewater treatment plant. The aim was to obtain a comprehensive picture of the coevolution dynamics in naturally evolving populations at relatively high time resolution. Coevolution was followed by monitoring changes over time in the CRISPR loci of Gordonia metagenome-assembled genome, and reciprocal changes in the viral genome. Genome-wide analysis indicated low strain variability of Gordonia, and almost clonal conservation of the trailer end of the CRISPR loci. Incorporation of newer spacers gave rise to multiple coexisting bacterial populations. The host population carrying a shorter CRISPR locus that contain only ancestral spacers, which has not acquired newer spacers against the coexisting phages, accounted for more than half of the total host abundance in the majority of samples. Phages genome co-evolved by introducing directional changes, with no preference for mutations within the protospacer and PAM regions. Metagenomic reconstruction of time-resolved variants of host and viral genomes revealed how the complexity at the population level has important consequences for bacteria-phage coexistence.Subject terms: Microbial ecology, Metagenomics, Bacteriophages     

The Chloroplast Genome Sequence of Mungbean (Vigna radiata) Determined by High-throughput Pyrosequencing: Structural Organization and Phylogenetic Relationships

Mungbean is an economically important crop which is grown principally for its protein-rich dry seeds. However, genomic research of mungbean has lagged behind other species in the Fabaceae family. Here, we reported the complete chloroplast (cp) genome sequence of mungbean obtained by the 454 pyrosequencing technology. The mungbean cp genome is 151 271 bp in length which includes a pair of inverted repeats (IRs) of 26 474 bp separated by a small single-copy region of 17 427 bp and a large single-copy region of 80 896 bp. The genome contains 108 unique genes and 19 of these genes are duplicated in the IR. Of these, 75 are predicted protein-coding genes, 4 ribosomal RNA genes and 29 tRNA genes. Relative to other plant cp genomes, we observed two distinct rearrangements: a 50-kb inversion between accD/rps16 and rbcL/trnK-UUU, and a 78-kb rearrangement between trnH/rpl14 and rps19/rps8. We detected sequence length polymorphism in the cp homopolymeric regions at the intra- and inter-specific levels in the Vigna species. Phylogenetic analysis demonstrated a close relationship between Vigna and Phaseolus in the phaseolinae subtribe and provided a strong support for a monophyletic group of the eurosid I.