Evolutionary constraints on the plastid tRNA set decoding methionine and isoleucine

A comparative analysis of plastome evolution in autotrophic Piperales

PREMISE

Many plastomes of autotrophic Piperales have been reported to date, describing a variety of differences. Most studies focused only on a few species or a single genus, and extensive, comparative analyses have not been done. Here, we reviewed publicly available plastome reconstructions for autotrophic Piperales, reanalyzed publicly available raw data, and provided new sequence data for all previously missing genera. Comparative plastome genomics of >100 autotrophic Piperales were performed.

METHODS

We performed de novo assemblies to reconstruct the plastomes of newly generated sequence data. We used Sanger sequencing and read mapping to verify the assemblies and to bridge assembly gaps. Furthermore, we reconstructed the phylogenetic relationships as a foundation for comparative plastome genomics.

RESULTS

We identified a plethora of assembly and annotation issues in published plastome data, which, if unattended, will lead to an artificial increase of diversity. We were able to detect patterns of missing and incorrect feature annotation and determined that the inverted repeat (IR) boundaries were the major source for erroneous assembly. Accounting for the aforementioned issues, we discovered relatively stable junctions of the IRs and the small single-copy region (SSC), whereas the majority of plastome variations among Piperales stems from fluctuations of the boundaries of the IR and the large single-copy (LSC) region.

CONCLUSIONS

This study of all available plastomes of autotrophic Piperales, expanded by new data for previously missing genera, highlights the IR-LSC junctions as a potential marker for discrimination of various taxonomic levels. Our data indicates a pseudogene-like status for cemA and ycf15 in various Piperales. Based on a review of published data, we conclude that incorrect IR-SSC boundary identification is the major source for erroneous plastome assembly. We propose a gold standard for assembly and annotation of high-quality plastomes based on de novo assembly methods and appropriate references for gene annotation.

Plastome variations reveal the distinct evolutionary scenarios of plastomes in the subfamily Cereoideae (Cactaceae)

BACKGROUND

The cactus family (Cactaceae) has been reported to have evolved a minimal photosynthetic plastome size, with the loss of inverted-repeat (IR) regions and NDH gene suites. However, there are very limited genomic data on the family, especially Cereoideae, the largest subfamily of cacti.

RESULTS

In the present study, we assembled and annotated 35 plastomes, 33 of which were representatives of Cereoideae, alongside 2 previously published plastomes. We analyzed the organelle genomes of 35 genera in the subfamily. These plastomes have variations rarely observed in those of other angiosperms, including size differences (with ~ 30 kb between the shortest and longest), dramatic dynamic changes in IR boundaries, frequent plastome inversions, and rearrangements. These results suggested that cacti have the most complex plastome evolution among angiosperms.

CONCLUSION

These results provide unique insight into the dynamic evolutionary history of Cereoideae plastomes and refine current knowledge of the relationships within the subfamily.

The plastome of Melocactus glaucescens Buining & Brederoo reveals unique evolutionary features and loss of essential tRNA genes

The plastome of Melocactus glaucescens shows unique rearrangements, IR expansion, and unprecedented gene losses in Cactaceae. Our data indicate tRNA import from the cytosol to the plastids in this species. Cactaceae represents one of the richest families in keystone species of arid and semiarid biomes. This family shows various specific features comprehending morphology, anatomy, and metabolism, which allow them to grow under unfavorable environmental conditions. The subfamily Cactoideae contains the most divergence of species, which are highly variable in growth habit and morphology. This subfamily includes the endangered species Melocactus glaucescens (tribe Cereeae), which is a cactus endemic to the biome Caatinga in Brazil. Aiming to analyze the plastid evolution and develop molecular markers, we sequenced and analyzed in detail the plastome of M. glaucescens. Our analyses revealed that the M. glaucescens plastome is the most divergent among the species of the family Cactaceae sequenced so far. We characterized here unique rearrangements, expanded IRs containing an unusual set of genes, and several gene losses. Some genes related to the ndh complex were lost during the plastome evolution, while others have lost their functionality. Additionally, the loss of three tRNA genes (trnA-UGC, trnV-UAC, and trnV-GAC) suggests tRNA import from the cytosol to the plastids in M. glaucescens. Moreover, we identified high gene divergence, several putative positive signatures, and possible unique RNA-editing sites. Furthermore, we mapped 169 SSRs in the plastome of M. glaucescens, which are helpful to access the genetic diversity of natural populations and conservation strategies. Finally, our data provide new insights into the evolution of plastids in Cactaceae, which is an outstanding lineage adapted to extreme environmental conditions and a notorious example of the atypical evolution of plastomes.

Overview of tRNA Modifications in Chloroplasts

The chloroplast is a promising platform for biotechnological innovation due to its compact translation machinery. Nucleotide modifications within a minimal set of tRNAs modulate codon–anticodon interactions that are crucial for translation efficiency. However, a comprehensive assessment of these modifications does not presently exist in chloroplasts. Here, we synthesize all available information concerning tRNA modifications in the chloroplast and assign translation efficiency for each modified anticodon–codon pair. In addition, we perform a bioinformatics analysis that links enzymes to tRNA modifications and aminoacylation in the chloroplast of Chlamydomonas reinhardtii. This work provides the first comprehensive analysis of codon and anticodon interactions of chloroplasts and its implication for translation efficiency.

Evolution and structural variations in chloroplast tRNAs in gymnosperms

Background

Chloroplast transfer RNAs (tRNAs) can participate in various vital processes. Gymnosperms have important ecological and economic value, and they are the dominant species in forest ecosystems in the Northern Hemisphere. However, the evolution and structural changes in chloroplast tRNAs in gymnosperms remain largely unclear.

Results

In this study, we determined the nucleotide evolution, phylogenetic relationships, and structural variations in 1779 chloroplast tRNAs in gymnosperms. The numbers and types of tRNA genes present in the chloroplast genomes of different gymnosperms did not differ greatly, where the average number of tRNAs was 33 and the frequencies of occurrence for various types of tRNAs were generally consistent. Nearly half of the anticodons were absent. Molecular sequence variation analysis identified the conserved secondary structures of tRNAs. About a quarter of the tRNA genes were found to contain precoded 3′ CCA tails. A few tRNAs have undergone novel structural changes that are closely related to their minimum free energy, and these structural changes affect the stability of the tRNAs. Phylogenetic analysis showed that tRNAs have evolved from multiple common ancestors. The transition rate was higher than the transversion rate in gymnosperm chloroplast tRNAs. More loss events than duplication events have occurred in gymnosperm chloroplast tRNAs during their evolutionary process.

Conclusions

These findings provide novel insights into the molecular evolution and biological characteristics of chloroplast tRNAs in gymnosperms.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08058-3.

Chloroplast Genome Evolution in Four Montane Zingiberaceae Taxa in China

Chloroplasts are critical to plant survival and adaptive evolution. The comparison of chloroplast genomes could provide insight into the adaptive evolution of closely related species. To identify potential adaptive evolution in the chloroplast genomes of four montane Zingiberaceae taxa (Cautleya, Roscoea, Rhynchanthus, and Pommereschea) that inhabit distinct habitats in the mountains of Yunnan, China, the nucleotide sequences of 13 complete chloroplast genomes, including five newly sequenced species, were characterized and compared. The five newly sequenced chloroplast genomes (162,878-163,831 bp) possessed typical quadripartite structures, which included a large single copy (LSC) region, a small single copy (SSC) region, and a pair of inverted repeat regions (IRa and IRb), and even though the structure was highly conserved among the 13 taxa, one of the rps19 genes was absent in Cautleya, possibly due to expansion of the LSC region. Positive selection of rpoA and ycf2 suggests that these montane species have experienced adaptive evolution to habitats with different sunlight intensities and that adaptation related to the chloroplast genome has played an important role in the evolution of Zingiberaceae taxa.

Combining tRNA sequencing methods to characterize plant tRNA expression and post-transcriptional modification

Differences in tRNA expression have been implicated in a remarkable number of biological processes. There is growing evidence that tRNA genes can play dramatically different roles depending on both expression and post-transcriptional modification, yet sequencing tRNAs to measure abundance and detect modifications remains challenging. Their secondary structure and extensive post-transcriptional modifications interfere with RNA-seq library preparation methods and have limited the utility of high-throughput sequencing technologies. Here, we combine two modifications to standard RNA-seq methods by treating with the demethylating enzyme AlkB and ligating with tRNA-specific adapters in order to sequence tRNAs from four species of flowering plants, a group that has been shown to have some of the most extensive rates of post-transcriptional tRNA modifications. This protocol has the advantage of detecting full-length tRNAs and sequence variants that can be used to infer many post-transcriptional modifications. We used the resulting data to produce a modification index of almost all unique reference tRNAs in Arabidopsis thaliana, which exhibited many anciently conserved similarities with humans but also positions that appear to be 'hot spots' for modifications in angiosperm tRNAs. We also found evidence based on northern blot analysis and droplet digital PCR that, even after demethylation treatment, tRNA-seq can produce highly biased estimates of absolute expression levels most likely due to biased reverse transcription. Nevertheless, the generation of full-length tRNA sequences with modification data is still promising for assessing differences in relative tRNA expression across treatments, tissues or subcellular fractions and help elucidate the functional roles of tRNA modifications.

The Origin and Evolution of Plastid Genome Downsizing in Southern Hemispheric Cypresses (Cupressaceae)

Plastome downsizing is rare in photosynthetic seed plants. However, a large-scale study of five cupressophyte families (conifers II) indicated that the plastomes of some Cupressaceous genera are notably reduced and compact. Here, we enriched taxon sampling in Cupressaceae by adding plastomes of ten previously unreported genera to determine the origin, evolution, and consequences of plastome reduction in this family. We discovered that plastome downsizing is specific to Callitroideae (a Southern Hemispheric subfamily). Their plastomes are the smallest, encode the fewest plastid genes, and contain the fewest GC-end codons among Cupressaceae. We show that repeated tRNA losses and shrinkage of intergenic spacers together contributed to the plastome downsizing in Callitroideae. Moreover, our absolute nucleotide substitution rate analyses suggest relaxed functional constraints in translation-related plastid genes (clpP, infA, rpl, and rps), but not in photosynthesis- or transcription-related ones, of Callitris (the most diverse genus in Callitroideae). We hypothesize that the small and low-GC plastomes of Callitroideae emerged ca. 112-75 million years ago as an adaptation to increased competition with angiosperms on the Gondwana supercontinent. Our findings highlight Callitroideae as another case of plastome downsizing in photosynthetic seed plant lineages.

Plastome of the mycoheterotrophic eudicot Exacum paucisquama (Gentianaceae) exhibits extensive gene loss and a highly expanded inverted repeat region

Mycoheterotrophic plants are highly specialized species able to acquire organic carbon from symbiotic fungi, with relaxed dependence on photosynthesis for carbon fixation. The relaxation of the functional constraint of photosynthesis and thereby the relaxed selective pressure on functional photosynthetic genes usually lead to substantial gene loss and a highly degraded plastid genome in heterotrophs. In this study, we sequenced and analyzed the plastome of the eudicot Exacum paucisquama, providing the first plastid genome of a mycoheterotroph in the family Gentianaceae to date. The E. paucisquama plastome was 44,028 bp in length, which is much smaller than the plastomes of autotrophic eudicots. Although the E. paucisquama plastome had a quadripartite structure, a distinct boundary shift was observed in comparison with the plastomes of other eudicots. We detected extensive gene loss and only 21 putative functional genes (15 protein-coding genes, four rRNA genes and two tRNA genes). Our results provide valuable information for comparative evolutionary analyses of plastomes of heterotrophic species belonging to different phylogenetic groups.

The Functions of Chloroplast Glutamyl-tRNA in Translation and Tetrapyrrole Biosynthesis

The chloroplast glutamyl-tRNA (tRNA^Glu) is unique in that it has two entirely different functions. In addition to acting in translation, it serves as the substrate of glutamyl-tRNA reductase (GluTR), the enzyme catalyzing the committed step in the tetrapyrrole biosynthetic pathway. How the tRNA^Glu pool is distributed between the two pathways and whether tRNA^Glu allocation limits tetrapyrrole biosynthesis and/or protein biosynthesis remains poorly understood. We generated a series of transplastomic tobacco (Nicotiana tabacum) plants to alter tRNA^Glu expression levels and introduced a point mutation into the plastid trnE gene, which has been reported to uncouple protein biosynthesis from tetrapyrrole biosynthesis in chloroplasts of the protist Euglena gracilis We show that, rather than comparable uncoupling of the two pathways, the trnE mutation is lethal in tobacco because it inhibits tRNA processing, thus preventing translation of Glu codons. Ectopic expression of the mutated trnE gene uncovered an unexpected inhibition of glutamyl-tRNA reductase by immature tRNA^Glu We further demonstrate that whereas overexpression of tRNA^Glu does not affect tetrapyrrole biosynthesis, reduction of GluTR activity through inhibition by tRNA^Glu precursors causes tetrapyrrole synthesis to become limiting in early plant development when active photosystem biogenesis provokes a high demand for de novo chlorophyll biosynthesis. Taken together, our findings provide insight into the roles of tRNA^Glu at the intersection of protein biosynthesis and tetrapyrrole biosynthesis.

Molecular evolution of chloroplast genomes in Monsteroideae (Araceae)

This study provides broad insight into the chloroplast genomes of the subfamily Monsteroideae. The identified polymorphic regions may be suitable for designing unique and robust molecular markers for phylogenetic inference. Monsteroideae is the third largest subfamily (comprises 369 species) and one of the early diverging lineages of the monocot plant family Araceae. The phylogeny of this important subfamily is not well resolved at the species level due to scarcity of genomic resources and suitable molecular markers. Here, we report annotated chloroplast genome sequences of four Monsteroideae species: Spathiphyllum patulinervum, Stenospermation multiovulatum, Monstera adansonii, and Rhaphidophora amplissima. The quadripartite chloroplast genomes (size range 163,335-164,751 bp) consist of a pair of inverted repeats (25,270-25,931 bp), separating a small single copy region (21,448-22,346 bp) from a large single copy region (89,714-91,841 bp). The genomes contain 114 unique genes, including four rRNA genes, 80 protein-coding genes, and 30 tRNA genes. Gene features, amino acid frequencies, codon usage, GC contents, oligonucleotide repeats, and inverted repeats dynamics exhibit similarities among the four genomes. Higher rate of synonymous substitutions was observed as compared to non-synonymous substitutions in 76 protein-coding genes. Positive selection was observed in seven protein-coding genes, including psbK, ndhK, ndhD, rbcL, accD, rps8, and ycf2. Our included species of Araceae showed the monophyly in Monsteroideae and other subfamilies. We report 30 suitable polymorphic regions. The polymorphic regions identified here might be suitable for designing unique and robust markers for inferring the phylogeny and phylogeography among closely related species within the genus Spathiphyllum and among distantly related species within the subfamily Monsteroideae. The chloroplast genomes presented here are a valuable contribution towards understanding the molecular evolutionary dynamics in the family Araceae.

tRNA functional signatures classify plastids as late-branching cyanobacteria

BACKGROUND

Eukaryotes acquired the trait of oxygenic photosynthesis through endosymbiosis of the cyanobacterial progenitor of plastid organelles. Despite recent advances in the phylogenomics of Cyanobacteria, the phylogenetic root of plastids remains controversial. Although a single origin of plastids by endosymbiosis is broadly supported, recent phylogenomic studies are contradictory on whether plastids branch early or late within Cyanobacteria. One underlying cause may be poor fit of evolutionary models to complex phylogenomic data.

RESULTS

Using Posterior Predictive Analysis, we show that recently applied evolutionary models poorly fit three phylogenomic datasets curated from cyanobacteria and plastid genomes because of heterogeneities in both substitution processes across sites and of compositions across lineages. To circumvent these sources of bias, we developed CYANO-MLP, a machine learning algorithm that consistently and accurately phylogenetically classifies ("phyloclassifies") cyanobacterial genomes to their clade of origin based on bioinformatically predicted function-informative features in tRNA gene complements. Classification of cyanobacterial genomes with CYANO-MLP is accurate and robust to deletion of clades, unbalanced sampling, and compositional heterogeneity in input tRNA data. CYANO-MLP consistently classifies plastid genomes into a late-branching cyanobacterial sub-clade containing single-cell, starch-producing, nitrogen-fixing ecotypes, consistent with metabolic and gene transfer data.

CONCLUSIONS

Phylogenomic data of cyanobacteria and plastids exhibit both site-process heterogeneities and compositional heterogeneities across lineages. These aspects of the data require careful modeling to avoid bias in phylogenomic estimation. Furthermore, we show that amino acid recoding strategies may be insufficient to mitigate bias from compositional heterogeneities. However, the combination of our novel tRNA-specific strategy with machine learning in CYANO-MLP appears robust to these sources of bias with high accuracy in phyloclassification of cyanobacterial genomes. CYANO-MLP consistently classifies plastids as late-branching Cyanobacteria, consistent with independent evidence from signature-based approaches and some previous phylogenetic studies.

De Novo Assembly Discovered Novel Structures in Genome of Plastids and Revealed Divergent Inverted Repeats in Mammillaria (Cactaceae, Caryophyllales)

The complete sequence of chloroplast genome (cpDNA) has been documented for single large columnar species of Cactaceae, lacking inverted repeats (IRs). We sequenced cpDNA for seven species of the short-globose cacti of Mammillaria and de novo assembly revealed three novel structures in land plants. These structures have a large single copy (LSC) that is 2.5 to 10 times larger than the small single copy (SSC), and two IRs that contain strong differences in length and gene composition. Structure 1 is distinguished by short IRs of <1 kb composed by rpl23-trnI-CAU-ycf2; with a total length of 110,189 bp and 113 genes. In structure 2, each IR is approximately 7.2 kb and is composed of 11 genes and one Intergenic Spacer-(psbK-trnQ)-trnQ-UUG-rps16-trnK-UUU-matK-trnK-UUU-psbA-trnH-GUG-rpl2-rpl23-trnI-CAU-ycf2; with a total size of 116,175 bp and 120 genes. Structure 3 has divergent IRs of approximately 14.1 kb, where IRA is composed of 20 genes: psbA-trnH-GUG-rpl23-trnI-CAU-ycf2-ndhB-rps7-rps12-trnV-GAC-rrn16-ycf68-trnI-GAU-trnA-AGC-rrn23-rrn4.5-rrn5-trnR-ACG-trnN-GUU-ndhF-rpl32; and IRB is identical to the IRA, but lacks rpl23. This structure has 131 genes and, by pseudogenization, it is shown to have the shortest cpDNA, of just 107,343 bp. Our findings show that Mammillaria bears an unusual structural diversity of cpDNA, which supports the elucidation of the evolutionary processes involved in cacti lineages.

Mitochondrial genome evolution in parasitic plants

Background

Parasitic plants rely on their host to cover their nutritional requirements either for their entire life or a smaller part of it. Depending on the level of parasitism, a proportional reduction on the plastid genome has been found. However, knowledge on gene loss and evolution of the mitogenome of parasitic plants is only available for four hemiparasitic Viscum species (Viscaceae), which lack many of the mitochondrial genes, while the remaining genes exhibit very fast molecular evolution rates. In this study, we include another genus, Phoradendron, from the Viscaceae, as well as 10 other hemiparasitic or holoparasitic taxa from across the phylogeny of the angiosperms to investigate how fast molecular evolution works on their mitogenomes, and the extent of gene loss.

Results

Our observations from Viscum were replicated in Phoradendron liga, whereas the remaining parasitic plants in the study have a complete set of the core mitochondrial genes and exhibit moderate or only slightly raised substitution rates compared to most autotrophic taxa, without any statistically significant difference between the different groups (autotrophs, hemiparasites and holoparasites). Additionally, further evidence is provided for the placement of Balanophoraceae within the order Santalales, while the exact placement of Cynomoriaceae still remains elusive.

Conclusions

We examine the mitochondrial gene content of 11 hemiparasitic and holoparasitic plants and confirm previous observations in Viscaceae. We show that the remaining parasitic plants do not have significantly higher substitution rates than autotrophic plants in their mitochondrial genes. We provide further evidence for the placement of Balanophoraceae in the Santalales.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1401-8) contains supplementary material, which is available to authorized users.

Plastid Genomes from Diverse Glaucophyte Genera Reveal a Largely Conserved Gene Content and Limited Architectural Diversity

Plastid genome (ptDNA) data of Glaucophyta have been limited for many years to the genus Cyanophora. Here, we sequenced the ptDNAs of Gloeochaete wittrockiana, Cyanoptyche gloeocystis, Glaucocystis incrassata, and Glaucocystis sp. BBH. The reported sequences are the first genome-scale plastid data available for these three poorly studied glaucophyte genera. Although the Glaucophyta plastids appear morphologically “ancestral,” they actually bear derived genomes not radically different from those of red algae or viridiplants. The glaucophyte plastid coding capacity is highly conserved (112 genes shared) and the architecture of the plastid chromosomes is relatively simple. Phylogenomic analyses recovered Glaucophyta as the earliest diverging Archaeplastida lineage, but the position of viridiplants as the first branching group was not rejected by the approximately unbiased test. Pairwise distances estimated from 19 different plastid genes revealed that the highest sequence divergence between glaucophyte genera is frequently higher than distances between species of different classes within red algae or viridiplants. Gene synteny and sequence similarity in the ptDNAs of the two Glaucocystis species analyzed is conserved. However, the ptDNA of Gla. incrassata contains a 7.9-kb insertion not detected in Glaucocystis sp. BBH. The insertion contains ten open reading frames that include four coding regions similar to bacterial serine recombinases (two open reading frames), DNA primases, and peptidoglycan aminohydrolases. These three enzymes, often encoded in bacterial plasmids and bacteriophage genomes, are known to participate in the mobilization and replication of DNA mobile elements. It is therefore plausible that the insertion in Gla. incrassata ptDNA is derived from a DNA mobile element.

AtRsgA from Arabidopsis thaliana is important for maturation of the small subunit of the chloroplast ribosome

Plastid ribosomes are very similar in structure and function to the ribosomes of their bacterial ancestors. Since ribosome biogenesis is not thermodynamically favorable under biological conditions it requires the activity of many assembly factors. Here we have characterized a homolog of bacterial RsgA in Arabidopsis thaliana and show that it can complement the bacterial homolog. Functional characterization of a strong mutant in Arabidopsis revealed that the protein is essential for plant viability, while a weak mutant produced dwarf, chlorotic plants that incorporated immature pre-16S ribosomal RNA into translating ribosomes. Physiological analysis of the mutant plants revealed smaller, but more numerous, chloroplasts in the mesophyll cells, reduction of chlorophyll a and b, depletion of proplastids from the rib meristem and decreased photosynthetic electron transport rate and efficiency. Comparative RNA sequencing and proteomic analysis of the weak mutant and wild-type plants revealed that various biotic stress-related, transcriptional regulation and post-transcriptional modification pathways were repressed in the mutant. Intriguingly, while nuclear- and chloroplast-encoded photosynthesis-related proteins were less abundant in the mutant, the corresponding transcripts were increased, suggesting an elaborate compensatory mechanism, potentially via differentially active retrograde signaling pathways. To conclude, this study reveals a chloroplast ribosome assembly factor and outlines the transcriptomic and proteomic responses of the compensatory mechanism activated during decreased chloroplast function.

Reconfiguration of the plastid genome in Lamprocapnos spectabilis: IR boundary shifting, inversion, and intraspecific variation

We generated a complete plastid genome (plastome) sequence for Lamprocapnos spectabilis, providing the first complete plastome from the subfamily Fumarioideae (Papaveraceae). The Lamprocapnos plastome shows large differences in size, structure, gene content, and substitution rates compared with two sequenced Papaveraceae plastomes. We propose a model that explains the major rearrangements observed, involving at least six inverted repeat (IR) boundary shifts and five inversions, generating a number of gene duplications and relocations, as well as a two-fold expansion of the IR and miniaturized small single-copy region. A reduction in the substitution rates for genes transferred from the single-copy regions to the IR was observed. Accelerated substitution rates of plastid accD and clpP were detected in the Lamprocapnos plastome. The accelerated substitution rate for the accD gene was correlated with a large insertion of amino acid repeat (AAR) motifs in the middle region, but the forces driving the higher substitution rate of the clpP gene are unclear. We found a variable number of AARs in Lamprocapnos accD and ycf1 genes within individuals, and the repeats were associated with coiled-coil regions. In addition, comparative analysis of three Papaveraceae plastomes revealed loss of rps15 in Papaver, and functional replacement to the nucleus was identified.

Ectopic Transplastomic Expression of a Synthetic MatK Gene Leads to Cotyledon-Specific Leaf Variegation

Chloroplasts (and other plastids) harbor their own genetic material, with a bacterial-like gene-expression systems. Chloroplast RNA metabolism is complex and is predominantly mediated by nuclear-encoded RNA-binding proteins. In addition to these nuclear factors, the chloroplast-encoded intron maturase MatK has been suggested to perform as a splicing factor for a subset of chloroplast introns. MatK is essential for plant cell survival in tobacco, and thus null mutants have not yet been isolated. We therefore attempted to over-express MatK from a neutral site in the chloroplast, placing it under the control of a theophylline-inducible riboswitch. This ectopic insertion of MatK lead to a variegated cotyledons phenotype. The addition of the inducer theophylline exacerbated the phenotype in a concentration-dependent manner. The extent of variegation was further modulated by light, sucrose and spectinomycin, suggesting that the function of MatK is intertwined with photosynthesis and plastid translation. Inhibiting translation in the transplastomic lines has a profound effect on the accumulation of several chloroplast mRNAs, including the accumulation of an RNA antisense to rpl33, a gene coding for an essential chloroplast ribosomal protein. Our study further supports the idea that MatK expression needs to be tightly regulated to prevent detrimental effects and establishes another link between leaf variegation and chloroplast translation.

GeSeq - versatile and accurate annotation of organelle genomes

We have developed the web application GeSeq (https://chlorobox.mpimp-golm.mpg.de/geseq.html) for the rapid and accurate annotation of organellar genome sequences, in particular chloroplast genomes. In contrast to existing tools, GeSeq combines batch processing with a fully customizable reference sequence selection of organellar genome records from NCBI and/or references uploaded by the user. For the annotation of chloroplast genomes, the application additionally provides an integrated database of manually curated reference sequences. GeSeq identifies genes or other feature-encoding regions by BLAT-based homology searches and additionally, by profile HMM searches for protein and rRNA coding genes and two de novo predictors for tRNA genes. These unique features enable the user to conveniently compare the annotations of different state-of-the-art methods, thus supporting high-quality annotations. The main output of GeSeq is a GenBank file that usually requires only little curation and is instantly visualized by OGDRAW. GeSeq also offers a variety of optional additional outputs that facilitate downstream analyzes, for example comparative genomic or phylogenetic studies.

Birth of Four Chimeric Plastid Gene Clusters in Japanese Umbrella Pine

Many genes in the plastid genomes (plastomes) of plants are organized as gene clusters, in which genes are co-transcribed, resembling bacterial operons. These plastid operons are highly conserved, even among conifers, whose plastomes are highly rearranged relative to other seed plants. We have determined the complete plastome sequence of Sciadopitys verticillata (Japanese umbrella pine), the sole member of Sciadopityaceae. The Sciadopitys plastome is characterized by extensive inversions, pseudogenization of four tRNA genes after tandem duplications, and a unique pair of 370-bp inverted repeats involved in the formation of isomeric plastomes. We showed that plastomic inversions in Sciadopitys have led to shuffling of the remote conserved operons, resulting in the birth of four chimeric gene clusters. Our data also demonstrated that the relocated genes can be co-transcribed in these chimeric gene clusters. The plastome of Sciadopitys advances our current understanding of how the conifer plastomes have evolved toward increased diversity and complexity.

Detecting and Characterizing the Highly Divergent Plastid Genome of the Nonphotosynthetic Parasitic Plant Hydnora visseri (Hydnoraceae)

Plastid genomes of photosynthetic flowering plants are usually highly conserved in both structure and gene content. However, the plastomes of parasitic and mycoheterotrophic plants may be released from selective constraint due to the reduction or loss of photosynthetic ability. Here we present the greatly reduced and highly divergent, yet functional, plastome of the nonphotosynthetic holoparasite Hydnora visseri (Hydnoraceae, Piperales). The plastome is 27 kb in length, with 24 genes encoding ribosomal proteins, ribosomal RNAs, tRNAs, and a few nonbioenergetic genes, but no genes related to photosynthesis. The inverted repeat and the small single copy region are only approximately 1.5 kb, and intergenic regions have been drastically reduced. Despite extreme reduction, gene order and orientation are highly similar to the plastome of Piper cenocladum, a related photosynthetic plant in Piperales. Gene sequences in Hydnora are highly divergent and several complementary approaches using the highest possible sensitivity were required for identification and annotation of this plastome. Active transcription is detected for all of the protein-coding genes in the plastid genome, and one of two introns is appropriately spliced out of rps12 transcripts. The whole-genome shotgun read depth is 1,400× coverage for the plastome, whereas the mitochondrial genome is covered at 40× and the nuclear genome at 2×. Despite the extreme reduction of the genome and high sequence divergence, the presence of syntenic, long transcriptionally active open-reading frames with distant similarity to other plastid genomes and a high plastome stoichiometry relative to the mitochondrial and nuclear genomes suggests that the plastome remains functional in H. visseri. A four-stage model of gene reduction, including the potential for complete plastome loss, is proposed to account for the range of plastid genomes in nonphotosynthetic plants.

A Member of the Arabidopsis Mitochondrial Transcription Termination Factor Family Is Required for Maturation of Chloroplast Transfer RNAIle(GAU)

Plastid gene expression is crucial for organelle function, but the factors that control it are still largely unclear. Members of the so-called mitochondrial transcription termination factor (mTERF) family are found in metazoans and plants and regulate organellar gene expression at different levels. Arabidopsis (Arabidopsis thaliana) mTERF6 is localized in chloroplasts and mitochondria, and its knockout perturbs plastid development and results in seedling lethality. In the leaky mterf6-1 mutant, a defect in photosynthesis is associated with reduced levels of photosystem subunits, although corresponding messenger RNA levels are unaffected, whereas translational capacity and maturation of chloroplast ribosomal RNAs (rRNAs) are perturbed in mterf6-1 mutants. Bacterial one-hybrid screening, electrophoretic mobility shift assays, and coimmunoprecipitation experiments reveal a specific interaction between mTERF6 and an RNA sequence in the chloroplast isoleucine transfer RNA gene (trnI.2) located in the rRNA operon. In vitro, recombinant mTERF6 bound to its plastid DNA target site can terminate transcription. At present, it is unclear whether disturbed rRNA maturation is a primary or secondary defect. However, it is clear that mTERF6 is required for the maturation of trnI.2. This points to an additional function of mTERFs.

The Highly Reduced Plastome of Mycoheterotrophic Sciaphila (Triuridaceae) Is Colinear with Its Green Relatives and Is under Strong Purifying Selection

The enigmatic monocot family Triuridaceae provides a potentially useful model system for studying the effects of an ancient loss of photosynthesis on the plant plastid genome, as all of its members are mycoheterotrophic and achlorophyllous. However, few studies have placed the family in a comparative context, and its phylogenetic placement is only partly resolved. It was also unclear whether any taxa in this family have retained a plastid genome. Here, we used genome survey sequencing to retrieve plastid genome data for Sciaphila densiflora (Triuridaceae) and ten autotrophic relatives in the orders Dioscoreales and Pandanales. We recovered a highly reduced plastome for Sciaphila that is nearly colinear with Carludovica palmata, a photosynthetic relative that belongs to its sister group in Pandanales, Cyclanthaceae-Pandanaceae. This phylogenetic placement is well supported and robust to a broad range of analytical assumptions in maximum-likelihood inference, and is congruent with recent findings based on nuclear and mitochondrial evidence. The 28 genes retained in the S. densiflora plastid genome are involved in translation and other nonphotosynthetic functions, and we demonstrate that nearly all of the 18 protein-coding genes are under strong purifying selection. Our study confirms the utility of whole plastid genome data in phylogenetic studies of highly modified heterotrophic plants, even when they have substantially elevated rates of substitution.

The Highly Reduced Plastome of Mycoheterotrophic Sciaphila (Triuridaceae) Is Colinear with Its Green Relatives and Is under Strong Purifying Selection

The enigmatic monocot family Triuridaceae provides a potentially useful model system for studying the effects of an ancient loss of photosynthesis on the plant plastid genome, as all of its members are mycoheterotrophic and achlorophyllous. However, few studies have placed the family in a comparative context, and its phylogenetic placement is only partly resolved. It was also unclear whether any taxa in this family have retained a plastid genome. Here, we used genome survey sequencing to retrieve plastid genome data for Sciaphila densiflora (Triuridaceae) and ten autotrophic relatives in the orders Dioscoreales and Pandanales. We recovered a highly reduced plastome for Sciaphila that is nearly colinear with Carludovica palmata, a photosynthetic relative that belongs to its sister group in Pandanales, Cyclanthaceae-Pandanaceae. This phylogenetic placement is well supported and robust to a broad range of analytical assumptions in maximum-likelihood inference, and is congruent with recent findings based on nuclear and mitochondrial evidence. The 28 genes retained in the S. densiflora plastid genome are involved in translation and other nonphotosynthetic functions, and we demonstrate that nearly all of the 18 protein-coding genes are under strong purifying selection. Our study confirms the utility of whole plastid genome data in phylogenetic studies of highly modified heterotrophic plants, even when they have substantially elevated rates of substitution.

Engineering plastid genomes: methods, tools, and applications in basic research and biotechnology

The small bacterial-type genome of the plastid (chloroplast) can be engineered by genetic transformation, generating cells and plants with transgenic plastid genomes, also referred to as transplastomic plants. The transformation process relies on homologous recombination, thereby facilitating the site-specific alteration of endogenous plastid genes as well as the precisely targeted insertion of foreign genes into the plastid DNA. The technology has been used extensively to analyze chloroplast gene functions and study plastid gene expression at all levels in vivo. Over the years, a large toolbox has been assembled that is now nearly comparable to the techniques available for plant nuclear transformation and that has enabled new applications of transplastomic technology in basic and applied research. This review describes the state of the art in engineering the plastid genomes of algae and land plants (Embryophyta). It provides an overview of the existing tools for plastid genome engineering, discusses current technological limitations, and highlights selected applications that demonstrate the immense potential of chloroplast transformation in several key areas of plant biotechnology.

Identification of enzymes for adenosine-to-inosine editing and discovery of cytidine-to-uridine editing in nucleus-encoded transfer RNAs of Arabidopsis

In all organisms, transfer RNAs (tRNAs) contain numerous modified nucleotides. For many base modifications in tRNAs, the functional significance is not well understood, and the enzymes performing the modification reactions are unknown. Here, we have studied members of a family of putative nucleotide deaminases in the model plant Arabidopsis (Arabidopsis thaliana). We show that two Arabidopsis genes encoding homologs of yeast (Saccharomyces cerevisiae) tRNA adenosine deaminases catalyze adenosine-to-inosine editing in position 34 of several cytosolic tRNA species. The encoded proteins (AtTAD2 and AtTAD3, for tRNA-specific adenosine deaminase) localize to the nucleus and interact with each other in planta in bimolecular fluorescence complementation and coimmunoprecipitation assays. Both AtTAD2 and AtTAD3 are encoded by essential genes whose knockout is lethal and leads to arrested embryo development at the globular stage. Knockdown mutants for AtTAD2 and AtTAD3 display reduced growth and inefficient editing from adenosine to inosine in six nucleus-encoded tRNA species. Moreover, upon comparison of DNA and complementary DNA sequences, we discovered cytidine-to-uridine RNA editing in position 32 of two nucleus-encoded serine tRNAs, tRNA-serine(AGA) and tRNA-serine(GCT). This adds a unique type of RNA editing to the modifications occurring in nuclear genome-encoded RNAs in plants.

The translational apparatus of plastids and its role in plant development

Chloroplasts (plastids) possess a genome and their own machinery to express it. Translation in plastids occurs on bacterial-type 70S ribosomes utilizing a set of tRNAs that is entirely encoded in the plastid genome. In recent years, the components of the chloroplast translational apparatus have been intensely studied by proteomic approaches and by reverse genetics in the model systems tobacco (plastid-encoded components) and Arabidopsis (nucleus-encoded components). This work has provided important new insights into the structure, function, and biogenesis of chloroplast ribosomes, and also has shed fresh light on the molecular mechanisms of the translation process in plastids. In addition, mutants affected in plastid translation have yielded strong genetic evidence for chloroplast genes and gene products influencing plant development at various levels, presumably via retrograde signaling pathway(s). In this review, we describe recent progress with the functional analysis of components of the chloroplast translational machinery and discuss the currently available evidence that supports a significant impact of plastid translational activity on plant anatomy and morphology.

Synthetic biology in plastids

Plastids (chloroplasts) harbor a small gene-dense genome that is amenable to genetic manipulation by transformation. During 1 billion years of evolution from the cyanobacterial endosymbiont to present-day chloroplasts, the plastid genome has undergone a dramatic size reduction, mainly as a result of gene losses and the large-scale transfer of genes to the nuclear genome. Thus the plastid genome can be regarded as a naturally evolved miniature genome, the gradual size reduction and compaction of which has provided a blueprint for the design of minimum genomes. Furthermore, because of the largely prokaryotic genome structure and gene expression machinery, the high transgene expression levels attainable in transgenic chloroplasts and the very low production costs in plant systems, the chloroplast lends itself to synthetic biology applications that are directed towards the efficient synthesis of green chemicals, biopharmaceuticals and other metabolites of commercial interest. This review describes recent progress with the engineering of plastid genomes with large constructs of foreign or synthetic DNA, and highlights the potential of the chloroplast as a model system in bottom-up and top-down synthetic biology approaches.

Analysis of the complete plastid genome of the unicellular red alga Porphyridium purpureum

We determined the complete nucleotide sequence of the plastid genome of the unicellular marine red alga Porphyridium purpureum strain NIES 2140, belonging to the unsequenced class Porphyridiophyceae. The genome is a circular DNA composed of 217,694 bp with the GC content of 30.3%. Twenty-nine of the 224 protein-coding genes contain one or multiple intron(s). A group I intron was found in the rpl28 gene, whereas the other introns were group II introns. The P. purpureum plastid genome has one non-coding RNA (ncRNA) gene, 29 tRNA genes and two nonidentical ribosomal RNA operons. One rRNA operon has a tRNA(Ala)(UGC) gene between the rrs and the rrl genes, whereas another has a tRNA(Ile)(GAU) gene. Phylogenetic analyses suggest that the plastids of Heterokontophyta, Cryptophyta and Haptophyta originated from the subphylum Rhodophytina. The order of the genes in the ribosomal protein cluster of the P. purpureum plastid genome differs from that of other Rhodophyta and Chromalveolata. These results suggest that a large-scale rearrangement occurred in the plastid genome of P. purpureum after its separation from other Rhodophyta.

A plastid without a genome: evidence from the nonphotosynthetic green algal genus Polytomella

Polytomella spp. are free-living, nonphotosynthetic green algae closely related to the model organism Chlamydomonas reinhardtii. Although colorless, Polytomella spp. have a plastid, but it is still unknown whether they harbor a plastid genome. We took a next generation sequencing approach, along with transcriptome sequencing, to search for a plastid genome and an associated gene expression system in Polytomella spp. Illumina sequencing of total DNA from four Polytomella spp. did not produce any recognizable plastid-derived reads but did generate a large number of mitochondrial DNA sequences. Transcriptomic analysis of Polytomella parva uncovered hundreds of putative nuclear-encoded, plastid-targeted proteins, which support the presence of plastid-based metabolic functions, similar to those observed in the plastids of other nonphotosynthetic algae. Conspicuously absent, however, were any plastid-targeted proteins involved in the expression, replication, or repair of plastid DNA. Based on these findings and earlier findings, we argue that the Polytomella genus represents the first well-supported example, to our knowledge, of a primary plastid-bearing lineage without a plastid genome.

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.

The tRNA 3'-end processing enzyme tRNase Z2 contributes to chloroplast biogenesis in rice

tRNase Z (TRZ) is a ubiquitous endonuclease that removes the 3'-trailer from precursor tRNAs during maturation. In yeast and animals, TRZ regulates the cell cycle via its (t)RNA processing activity; however, its physiological function in higher plants has not been well characterized. This study describes the identification of a rice (Oryza sativa) TRZ2 mutant; plants homozygous for the osatrz2 mutation were albinos with deficient chlorophyll content. A microscopic analysis of the mutant plants revealed that the transition of proplastids to chloroplasts was arrested at an early stage, and the number and size of the plastids in callus cells was substantially decreased. A genetic complementation test and an RNA interference analysis confirmed that disruption of OsaTRZ2 was responsible for the mutant phenotype. OsaTRZ2 is expressed in all rice tissues, but is preferentially expressed in leaves, sheathes, and calli. OsaTRZ2 was subcellularly localized in chloroplasts, and displayed tRNA 3'-end processing activity in both in vitro and in vivo assays. In the osatrz2 mutants, transcription of plastid-encoded and nucleus-encoded RNA polymerases was severely reduced and moderately increased, respectively. These results suggest that the tRNA 3' processing activity of OsaTRZ2 contributes to chloroplast biogenesis.

Reverse genetics in complex multigene operons by co-transformation of the plastid genome and its application to the open reading frame previously designated psbN

Reverse genetics approaches have contributed enormously to the elucidation of gene functions in plastid genomes and the determination of structure-function relationships in chloroplast multiprotein complexes. Gene knock-outs are usually performed by disrupting the reading frame of interest with a selectable marker cassette. Site-directed mutagenesis is done by placing the marker into the adjacent intergenic spacer and relying on co-integration of the desired mutation by homologous recombination. These strategies are not applicable to genes residing in large multigene operons or other gene-dense genomic regions, because insertion of the marker cassette into an operon-internal gene or into the nearest intergenic spacer is likely to interfere with expression of adjacent genes in the operon or disrupt cis-elements for the expression of neighboring genes and operons. Here we have explored the possibility of using a co-transformation strategy to mutate a small gene of unknown function (psbN) that is embedded in a complex multigene operon. Although inactivation of psbN resulted in strong impairment of photosynthesis, homoplasmic knock-out lines were readily recovered by co-transformation with a selectable marker integrating >38 kb away from the targeted psbN. Our results suggest co-transformation as a suitable strategy for the functional analysis of plastid genes and operons, which allows the recovery of unselected homoplasmic mutants even if the introduced mutations entail a significant selective disadvantage. Moreover, our data provide evidence for involvement of the psbN gene product in the biogenesis of both photosystem I and photosystem II. We therefore propose to rename the gene product 'photosystem biogenesis factor 1' and the gene pbf1.

Importance of adenosine-to-inosine editing adjacent to the anticodon in an Arabidopsis alanine tRNA under environmental stress

In all organisms, transfer RNAs (tRNAs) undergo extensive post-transcriptional modifications. Although base modifications in the anticodon are known to alter decoding specificity or improve decoding accuracy, much less is known about the functional relevance of modifications in other positions of tRNAs. Here, we report the identification of an A-to-I tRNA editing enzyme that modifies the tRNA-Ala(AGC) in the model plant Arabidopsis thaliana. The enzyme is homologous to Tad1p, a yeast tRNA-specific adenosine deaminase, and it selectively deaminates the adenosine in the position 3'-adjacent to the anticodon (A37) to inosine. We show that the AtTAD1 protein is exclusively localized in the nucleus. The tad1 loss-of-function mutants isolated in Arabidopsis show normal accumulation of the tRNA-Ala(AGC), suggesting that the loss of the I37 modification does not affect tRNA stability. The tad1 knockout mutants display no discernible phenotype under standard growth conditions, but produce less biomass under environmental stress conditions. Our results provide the first evidence in support of a physiological relevance of the A37-to-I modification in eukaryotes.

The contributions of wobbling and superwobbling to the reading of the genetic code

Reduced bacterial genomes and most genomes of cell organelles (chloroplasts and mitochondria) do not encode the full set of 32 tRNA species required to read all triplets of the genetic code according to the conventional wobble rules. Superwobbling, in which a single tRNA species that contains a uridine in the wobble position of the anticodon reads an entire four-fold degenerate codon box, has been suggested as a possible mechanism for how tRNA sets can be reduced. However, the general feasibility of superwobbling and its efficiency in the various codon boxes have remained unknown. Here we report a complete experimental assessment of the decoding rules in a typical prokaryotic genetic system, the plastid genome. By constructing a large set of transplastomic knock-out mutants for pairs of isoaccepting tRNA species, we show that superwobbling occurs in all codon boxes where it is theoretically possible. Phenotypic characterization of the transplastomic mutant plants revealed that the efficiency of superwobbling varies in a codon box-dependent manner, but--contrary to previous suggestions--it is independent of the number of hydrogen bonds engaged in codon-anticodon interaction. Finally, our data provide experimental evidence of the minimum tRNA set comprising 25 tRNA species, a number lower than previously suggested. Our results demonstrate that all triplets with pyrimidines in third codon position are dually decoded: by a tRNA species utilizing standard base pairing or wobbling and by a second tRNA species employing superwobbling. This has important implications for the interpretation of the genetic code and will aid the construction of synthetic genomes with a minimum-size translational apparatus.