Uncovering the protein translocon at the chloroplast inner envelope membrane

The Complete Chloroplast Genome of Euphrasia regelii, Pseudogenization of ndh Genes and the Phylogenetic Relationships Within Orobanchaceae

Euphrasia (Orobanchaceae) is a genus which is widely distributed in temperate regions of the southern and northern hemisphere. The taxonomy of Euphrasia is still controversial due to the similarity of morphological characters and a lack of genomic resources. Here, we present the first complete chloroplast (cp) genome of this taxonomically challenging genus. The cp genome of Euphrasia regelii consists of 153,026 bp, including a large single-copy region (83,893 bp), a small single-copy region (15,801 bp) and two inverted repeats (26,666 bp). There are 105 unique genes, including 71 protein-coding genes, 30 tRNA and 4 rRNA genes. Although the structure and gene order is comparable to the one in other angiosperm cp genomes, genes encoding the NAD(P)H dehydrogenase complex are widely pseudogenized due to mutations resulting in frameshifts, and stop codon positions. We detected 36 dispersed repeats, 7 tandem repeats and 65 simple sequence repeat loci in the E. regelii plastome. Comparative analyses indicated that the cp genome of E. regelii is more conserved compared to other hemiparasitic taxa in the Pedicularideae and Buchnereae. No structural rearrangements or loss of genes were detected. Our analyses suggested that three genes (clpP, ycf2 and rps14) were under positive selection and other genes under purifying selection. Phylogenetic analysis of monophyletic Orobanchaceae based on 45 plastomes indicated a close relationship between E. regelii and Neobartsia inaequalis. In addition, autotrophic lineages occupied the earliest diverging branches in our phylogeny, suggesting that autotrophy is the ancestral trait in this parasitic family.

Complete plastome sequences of Picea asperata and P. crassifolia and comparative analyses with P. abies and P. morrisonicola

Picea asperata and P. crassifolia have sympatric ranges and are closely related, but the differences between these species at the plastome level are unknown. To better understand the patterns of variation among Picea plastomes, the complete plastomes of P. asperata and P. crassifolia were sequenced. Then, the plastomes were compared with the complete plastomes of P. abies and P. morrisonicola, which are closely and distantly related to the focal species, respectively. We also used these sequences to construct phylogenetic trees to determine the relationships among and between the four species as well as additional taxa from Pinaceae and other gymnosperms. Analysis of our sequencing data allowed us to identify 438 single nucleotide polymorphism (SNPs) point mutation events, 95 indel events, four inversion events, and seven highly variable regions, including six gene spacer regions (psbJ-petA, trnT-psaM, trnS-trnD, trnL-rps4, psaC-ccsA, and rps7-trnL) and one gene (ycf1). The highly variable regions are appropriate targets for future use in the phylogenetic reconstructions of closely related, sympatric species of Picea as well as Pinaceae in general.

Reduced expression of the proteolytically inactive FtsH members has impacts on the Darwinian fitness of Arabidopsis thaliana

FtsH (filamentation-temperature-sensitive protein H) proteases are a family of membrane-bound enzymes present in eubacteria, animals, and plants. Besides the 12 genes encoding proteolytically active members of the FtsH family in the genome of Arabidopsis, there are five genes coding for members that are assumed to be proteolytically inactive due to mutations in the protease domain; these are termed FtsHi (i for inactive). Despite their lack of proteolytic activity, these FtsHi members seem to be important for chloroplast and plant development as four out of five homozygous knockout-mutants of FtsHis are embryo-lethal. Here, we analysed the Darwinian fitness of weak homozygous (ftshi1,3,4) and heterozygous (ftshi/FTSHi2,4,5) mutants. We compared the growth and development of these mutants to their respective wild-type Arabidopsis plants under controlled laboratory conditions and in the field, and we also evaluated the photosynthetic efficiency by pulse-amplitude modulation fluorescence. Homologous genotypes were subjected to various stress conditions in a greenhouse and gene co-expression as well as phylogenetic analyses were performed. Analysis of the gene-expression network of the five FTSHi genes indicated common clusters with genes encoding FtsH12, OTP51, and methylase. Phylogenetic analyses pointed to a common evolution (and common disappearance in grasses and gymnosperms) of FtsH12 and multiple presumably proteolytically inactive FtsHi enzymes. Our data show that the FtsHi enzymes are highly important during the seedling stage and for Darwinian fitness analyses in semi-natural conditions.

Chloroplast competition is controlled by lipid biosynthesis in evening primroses

In most eukaryotes, organellar genomes are transmitted preferentially by the mother, but molecular mechanisms and evolutionary forces underlying this fundamental biological principle are far from understood. It is believed that biparental inheritance promotes competition between the cytoplasmic organelles and allows the spread of so-called selfish cytoplasmic elements. Those can be, for example, fast-replicating or aggressive chloroplasts (plastids) that are incompatible with the hybrid nuclear genome and therefore maladaptive. Here we show that the ability of plastids to compete against each other is a metabolic phenotype determined by extremely rapidly evolving genes in the plastid genome of the evening primrose Oenothera Repeats in the regulatory region of accD (the plastid-encoded subunit of the acetyl-CoA carboxylase, which catalyzes the first and rate-limiting step of lipid biosynthesis), as well as in ycf2 (a giant reading frame of still unknown function), are responsible for the differences in competitive behavior of plastid genotypes. Polymorphisms in these genes influence lipid synthesis and most likely profiles of the plastid envelope membrane. These in turn determine plastid division and/or turnover rates and hence competitiveness. This work uncovers cytoplasmic drive loci controlling the outcome of biparental chloroplast transmission. Here, they define the mode of chloroplast inheritance, as plastid competitiveness can result in uniparental inheritance (through elimination of the "weak" plastid) or biparental inheritance (when two similarly "strong" plastids are transmitted).

Consequences of impaired 1-MDa TIC complex assembly for the abundance and composition of chloroplast high-molecular mass protein complexes

We report a systematic analysis of chloroplast high-molecular mass protein complexes using a combination of native gel electrophoresis and absolute protein quantification by MS^E. With this experimental setup, we characterized the effect of the tic56-3 mutation in the 1-MDa inner envelope translocase (TIC) on the assembly of the chloroplast proteome. We show that the tic56-3 mutation results in a reduction of the 1-MDa TIC complex to approximately 10% of wildtype levels. Hierarchical clustering confirmed the association of malate dehydrogenase (MDH) with an envelope-associated FtsH/FtsHi complex and suggested the association of a glycine-rich protein with the 1-MDa TIC complex. Depletion of this complex leads to a reduction of chloroplast ATPase to approx. 75% of wildtype levels, while the abundance of the FtsH/FtsHi complex is increased to approx. 140% of wildtype. The accumulation of the major photosynthetic complexes is not affected by the mutation, suggesting that tic56-3 plants can sustain a functional photosynthetic machinery despite a significant reduction of the 1-MDa TIC complex. Together our analysis expands recent efforts to catalogue the native molecular masses of chloroplast proteins and provides information on the consequences of impaired accumulation of the 1-MDa TIC translocase for chloroplast proteome assembly.

Structural considerations of folded protein import through the chloroplast TOC/TIC translocons

Protein import into chloroplasts is carried out by the protein translocons at the outer and inner envelope membranes (TOC and TIC). Detailed structures for these translocons are lacking, with only a low-resolution TOC complex structure available. Recently, we showed that the TOC/TIC translocons can import folded proteins, a rather unique feat for a coupled double membrane system. We also determined the maximum functional TOC/TIC pore size to be 30-35 Å. Here, we discuss how such large pores could form and compare the structural dynamics of the pore-forming Toc75 subunit to its bacterial/mitochondrial Omp85 family homologs. We put forward structural models that can be empirically tested and also briefly review the pore dynamics of other protein translocons with known structures.

Plastome-Wide Rearrangements and Gene Losses in Carnivorous Droseraceae

The plastid genomes of four related carnivorous plants (Drosera regia, Drosera erythrorhiza, Aldrovanda vesiculosa, and Dionaea muscipula) were sequenced to examine changes potentially induced by the transition to carnivory. The plastid genomes of the Droseraceae show multiple rearrangements, gene losses, and large expansions or contractions of the inverted repeat. All the ndh genes are lost or nonfunctional, as well as in some of the species, clpP1, ycf1, ycf2 and some tRNA genes. Uniquely, among land plants, the trnK gene has no intron. Carnivory in the Droseraceae coincides with changes in plastid gene content similar to those induced by parasitism and mycoheterotrophy, suggesting parallel changes in chloroplast function due to the similar switch from autotrophy to (mixo-) heterotrophy. A molecular phylogeny of the taxa based on all shared plastid genes indicates that the "snap-traps" of Aldrovanda and Dionaea have a common origin.

Comparative Analysis of the Complete Chloroplast Genome Sequences of Three Closely Related East-Asian Wild Roses (Rosa sect. Synstylae; Rosaceae)

Species belonging to Rosa section Synstylae (Rosaceae) are mainly distributed in East Asia, and represent recently diverged lineages within the genus. Over decades, inferring phylogenetic relationships within section Synstylae have been exceptional challenges, due to short branch lengths and low support values. Of approximately 36 species in the section Synstylae, Rosa multiflora, Rosa luciae and Rosa maximowicziana are widely distributed in the Sino-Japanese floristic region. In this study, we assembled chloroplast genomes of these three species to compare the genomic features within section Synstylae, and to compare with other infrageneric groups. We found that three Rosa sect. Synstylae species had lost infA genes with pseudogenization, and they were almost identical to each other. Two protein-coding gene regions (ndhF and ycf1) and five non-coding regions (5'matK-trnK, psbI-trnS-trnG, rps16-trnG, rpoB-trnC, and rps4-trnT) were identified as being highly informative markers. Within three section Synstylae chloroplast genomes, 85 simple sequence repeat (SSR) motifs were detected, of which at least 13 motifs were identified to be effective markers. The phylogenetic relationships of R. multiflora, R. luciae and R. maximowicziana could not be resolved, even with chloroplast genome-wide data. This study reveals the chloroplast genomic data of Rosa sect. Synstylae, and it provides valuable markers for DNA barcoding and phylogenetic analyses for further studies.

A Novel Prokaryote-Type ECF/ABC Transporter Module in Chloroplast Metal Homeostasis

During evolution, chloroplasts, which originated by endosymbiosis of a prokaryotic ancestor of today's cyanobacteria with a eukaryotic host cell, were established as the site for photosynthesis. Therefore, chloroplast organelles are loaded with transition metals including iron, copper, and manganese, which are essential for photosynthetic electron transport due to their redox capacity. Although transport, storage, and cofactor-assembly of metal ions in chloroplasts are tightly controlled and crucial throughout plant growth and development, knowledge on the molecular nature of chloroplast metal-transport proteins is still fragmentary. Here, we characterized the soluble, ATP-binding ABC-transporter subunits ABCI10 and ABCI11 in Arabidopsis thaliana, which show similarities to components of prokaryotic, multisubunit ABC transporters. Both ABCI10 and ABCI11 proteins appear to be strongly attached to chloroplast-intrinsic membranes, most likely inner envelopes for ABCI10 and possibly plastoglobuli for ABCI11. Loss of ABCI10 and ABCI11 gene products in Arabidopsis leads to extremely dwarfed, albino plants showing impaired chloroplast biogenesis and deregulated metal homeostasis. Further, we identified the membrane-intrinsic protein ABCI12 as potential interaction partner for ABCI10 in the inner envelope. Our results suggest that ABCI12 inserts into the chloroplast inner envelope membrane most likely with five predicted α-helical transmembrane domains and represents the membrane-intrinsic subunit of a prokaryotic-type, energy-coupling factor (ECF) ABC-transporter complex. In bacteria, these multisubunit ECF importers are widely distributed for the uptake of nickel and cobalt metal ions as well as for import of vitamins and several other metabolites. Therefore, we propose that ABCI10 (as the ATPase A-subunit) and ABCI12 (as the membrane-intrinsic, energy-coupling T-subunit) are part of a novel, chloroplast envelope-localized, AAT energy-coupling module of a prokaryotic-type ECF transporter, most likely involved in metal ion uptake.

Plastid Genomes of Five Species of Riverweeds (Podostemaceae): Structural Organization and Comparative Analysis in Malpighiales

With the advent of next-generation sequencing technologies, whole-plastome data can be obtained as a byproduct of low-coverage sequencing of the plant genomic DNA. This provides an opportunity to study plastid evolution across groups, as well as testing phylogenetic relationships among taxa. Within the order Malpighiales (∼16,000 spp.), the Podostemaceae (∼300 spp.) stand out for their unique habit, living attached to rocks in fast-flowing aquatic habitats, and displaying highly modified morphologies that confound our understanding of their classification, biology, and evolution. In this study, we used genome skimming data to assemble the full plastid genome of 5 species within Podostemaceae. We analyzed our data in a comparative framework within Malpighiales to determine the structure, gene content, and rearrangements in the plastomes of the family. The Podostemaceae have one of the smallest plastid genomes reported so far for the Malpighiales, possibly due to variation in length of inverted repeat (IR) regions, gene loss, and intergenic region variation. We also detected a major inversion in the large single-copy region unique to the family. The uncommon loss or pseudogenization of ycf1 and ycf2 in angiosperms and in land plants in general is also found to be characteristic of Podostemaceae, but the compensatory mechanisms and implications of this and of the pseudogenization of accD, rpl22, and clpP and loss of rps16 remain to be explained in this group. In addition, we estimated a phylogenetic tree among selected species in Malpighiales. Our findings indicate that the Podostemaceae are a distinct lineage with long branches that suggest faster rates of evolution in the plastome of the group, compared with other taxa in the order. This study lays the foundations for future phylogenomic studies in the family.

Targeted Control of Chloroplast Quality to Improve Plant Acclimation: From Protein Import to Degradation

The chloroplast is an important energy-producing organelle acting as an environmental sensor for the plant cell. The normal turnover of the entire damaged chloroplast and its specific components is required for efficient photosynthesis and other metabolic reactions under stress conditions. Nuclear-encoded proteins must be imported into the chloroplast through different membrane transport complexes, and the orderly protein import plays an important role in plant adaptive regulation. Under adverse environmental conditions, the damaged chloroplast or its specific components need to be degraded efficiently to ensure normal cell function. In this review, we discuss the molecular mechanism of protein import and degradation in the chloroplast. Specifically, quality control of chloroplast from protein import to degradation and associated regulatory pathways are discussed to better understand how plants adapt to environmental stress by fine-tuning chloroplast homeostasis, which will benefit breeding approaches to improve crop yield.

A Novel Prokaryote-Type ECF/ABC Transporter Module in Chloroplast Metal Homeostasis

During evolution, chloroplasts, which originated by endosymbiosis of a prokaryotic ancestor of today's cyanobacteria with a eukaryotic host cell, were established as the site for photosynthesis. Therefore, chloroplast organelles are loaded with transition metals including iron, copper, and manganese, which are essential for photosynthetic electron transport due to their redox capacity. Although transport, storage, and cofactor-assembly of metal ions in chloroplasts are tightly controlled and crucial throughout plant growth and development, knowledge on the molecular nature of chloroplast metal-transport proteins is still fragmentary. Here, we characterized the soluble, ATP-binding ABC-transporter subunits ABCI10 and ABCI11 in Arabidopsis thaliana, which show similarities to components of prokaryotic, multisubunit ABC transporters. Both ABCI10 and ABCI11 proteins appear to be strongly attached to chloroplast-intrinsic membranes, most likely inner envelopes for ABCI10 and possibly plastoglobuli for ABCI11. Loss of ABCI10 and ABCI11 gene products in Arabidopsis leads to extremely dwarfed, albino plants showing impaired chloroplast biogenesis and deregulated metal homeostasis. Further, we identified the membrane-intrinsic protein ABCI12 as potential interaction partner for ABCI10 in the inner envelope. Our results suggest that ABCI12 inserts into the chloroplast inner envelope membrane most likely with five predicted α-helical transmembrane domains and represents the membrane-intrinsic subunit of a prokaryotic-type, energy-coupling factor (ECF) ABC-transporter complex. In bacteria, these multisubunit ECF importers are widely distributed for the uptake of nickel and cobalt metal ions as well as for import of vitamins and several other metabolites. Therefore, we propose that ABCI10 (as the ATPase A-subunit) and ABCI12 (as the membrane-intrinsic, energy-coupling T-subunit) are part of a novel, chloroplast envelope-localized, AAT energy-coupling module of a prokaryotic-type ECF transporter, most likely involved in metal ion uptake.

The Chloroplast Envelope Protease FTSH11 - Interaction With CPN60 and Identification of Potential Substrates

FTSH proteases are membrane-bound, ATP-dependent metalloproteases found in bacteria, mitochondria and chloroplasts. The product of one of the 12 genes encoding FTSH proteases in Arabidopsis, FTSH11, has been previously shown to be essential for acquired thermotolerance. However, the substrates of this protease, as well as the mechanism linking it to thermotolerance are largely unknown. To get insight into these, the FTSH11 knockout mutant was complemented with proteolytically active or inactive variants of this protease, tagged with HA-tag, under the control of the native promoter. Using these plants in thermotolerance assay demonstrated that the proteolytic activity, and not only the ATPase one, is essential for conferring thermotolerance. Immunoblot analyses of leaf extracts, isolated organelles and sub-fractionated chloroplast membranes localized FTSH11 mostly to chloroplast envelopes. Affinity purification followed by mass spectrometry analysis revealed interaction between FTSH11 and different components of the CPN60 chaperonin. In affinity enrichment assays, CPN60s as well as a number of envelope, stroma and thylakoid proteins were found associated with proteolytically inactive FTSH11. Comparative proteomic analysis of WT and knockout plants, grown at 20°C or exposed to 30°C for 6 h, revealed a plethora of upregulated chloroplast proteins in the knockout, some of them might be candidate substrates. Among these stood out TIC40, which was stabilized in the knockout line after recovery from heat stress, and three proteins that were found trapped in the affinity enrichment assay: the nucleotide antiporter PAPST2, the fatty acid binding protein FAP1 and the chaperone HSP70. The consistent behavior of these four proteins in different assays suggest that they are potential FTSH11 substrates.

Evolutionary Analysis of Plastid Genomes of Seven Lonicera L. Species: Implications for Sequence Divergence and Phylogenetic Relationships

Plant plastomes play crucial roles in species evolution and phylogenetic reconstruction studies due to being maternally inherited and due to the moderate evolutionary rate of genomes. However, patterns of sequence divergence and molecular evolution of the plastid genomes in the horticulturally- and economically-important Lonicera L. species are poorly understood. In this study, we collected the complete plastomes of seven Lonicera species and determined the various repeat sequence variations and protein sequence evolution by comparative genomic analysis. A total of 498 repeats were identified in plastid genomes, which included tandem (130), dispersed (277), and palindromic (91) types of repeat variations. Simple sequence repeat (SSR) elements analysis indicated the enriched SSRs in seven genomes to be mononucleotides, followed by tetra-nucleotides, dinucleotides, tri-nucleotides, hex-nucleotides, and penta-nucleotides. We identified 18 divergence hotspot regions (rps15, rps16, rps18, rpl23, psaJ, infA, ycf1, trnN-GUU-ndhF, rpoC2-rpoC1, rbcL-psaI, trnI-CAU-ycf2, psbZ-trnG-UCC, trnK-UUU-rps16, infA-rps8, rpl14-rpl16, trnV-GAC-rrn16, trnL-UAA intron, and rps12-clpP) that could be used as the potential molecular genetic markers for the further study of population genetics and phylogenetic evolution of Lonicera species. We found that a large number of repeat sequences were distributed in the divergence hotspots of plastid genomes. Interestingly, 16 genes were determined under positive selection, which included four genes for the subunits of ribosome proteins (rps7, rpl2, rpl16, and rpl22), three genes for the subunits of photosystem proteins (psaJ, psbC, and ycf4), three NADH oxidoreductase genes (ndhB, ndhH, and ndhK), two subunits of ATP genes (atpA and atpB), and four other genes (infA, rbcL, ycf1, and ycf2). Phylogenetic analysis based on the whole plastome demonstrated that the seven Lonicera species form a highly-supported monophyletic clade. The availability of these plastid genomes provides important genetic information for further species identification and biological research on Lonicera.

En route into chloroplasts: preproteins' way home

Chloroplasts are the characteristic endosymbiotic organelles of plant cells which during the course of evolution lost most of their genetic information to the nucleus. Thus, they critically depend on the host cell for allocation of nearly their complete protein supply. This includes gene expression, translation, protein targeting, and transport-all of which need to be tightly regulated and perfectly coordinated to accommodate the cells' needs. To this end, multiple signaling pathways have been implemented that interchange information between the different cellular compartments. One of the most complex and energy consuming processes is the translocation of chloroplast-destined proteins into their target organelle. It is a concerted effort from chaperones, receptor proteins, channels, and regulatory elements to ensure correct targeting, efficient transport, and subsequent folding. Although we have discovered and learned a lot about protein import into chloroplasts in the last decades, there are still many open questions and debates about the roles of individual proteins as well as the mechanistic details. In this review, I will summarize and discuss the published data with a focus on the translocation complex in the chloroplast inner envelope membrane.

Evolution of protein transport to the chloroplast envelope membranes

Chloroplasts are descendants of an ancient endosymbiotic cyanobacterium that lived inside a eukaryotic cell. They inherited the prokaryotic double membrane envelope from cyanobacteria. This envelope contains prokaryotic protein sorting machineries including a Sec translocase and relatives of the central component of the bacterial outer membrane β-barrel assembly module. As the endosymbiont was integrated with the rest of the cell, the synthesis of most of its proteins shifted from the stroma to the host cytosol. This included nearly all the envelope proteins identified so far. Consequently, the overall biogenesis of the chloroplast envelope must be distinct from cyanobacteria. Envelope proteins initially approach their functional locations from the exterior rather than the interior. In many cases, they have been shown to use components of the general import pathway that also serves the stroma and thylakoids. If the ancient prokaryotic protein sorting machineries are still used for chloroplast envelope proteins, their activities must have been modified or combined with the general import pathway. In this review, we analyze the current knowledge pertaining to chloroplast envelope biogenesis and compare this to bacteria.

Complete Chloroplast Genome Sequence and Phylogenetic Analysis of Aster tataricus

We sequenced and analyzed the complete chloroplast genome of Aster tataricus (family Asteraceae), a Chinese herb used medicinally to relieve coughs and reduce sputum. The A. tataricus chloroplast genome was 152,992 bp in size, and harbored a pair of inverted repeat regions (IRa and IRb, each 24,850 bp) divided into a large single-copy (LSC, 84,698 bp) and a small single-copy (SSC, 18,250 bp) region. Our annotation revealed that the A. tataricus chloroplast genome contained 115 genes, including 81 protein-coding genes, 4 ribosomal RNA genes, and 30 transfer RNA genes. In addition, 70 simple sequence repeats (SSRs) were detected in the A. tataricus chloroplast genome, including mononucleotides (36), dinucleotides (1), trinucleotides (23), tetranucleotides (1), pentanucleotides (8), and hexanucleotides (1). Comparative chloroplast genome analysis of three Aster species indicated that a higher similarity was preserved in the IR regions than in the LSC and SSC regions, and that the differences in the degree of preservation were slighter between A. tataricus and A. altaicus than between A. tataricus and A. spathulifolius. Phylogenetic analysis revealed that A. tataricus was more closely related to A. altaicus than to A. spathulifolius. Our findings offer valuable information for future research on Aster species identification and selective breeding.

Complete Chloroplast Genome Sequence and Phylogenetic Analysis of Quercus acutissima

Quercus acutissima, an important endemic and ecological plant of the Quercus genus, is widely distributed throughout China. However, there have been few studies on its chloroplast genome. In this study, the complete chloroplast (cp) genome of Q. acutissima was sequenced, analyzed, and compared to four species in the Fagaceae family. The size of the Q. acutissima chloroplast genome is 161,124 bp, including one large single copy (LSC) region of 90,423 bp and one small single copy (SSC) region of 19,068 bp, separated by two inverted repeat (IR) regions of 51,632 bp. The GC content of the whole genome is 36.08%, while those of LSC, SSC, and IR are 34.62%, 30.84%, and 42.78%, respectively. The Q. acutissima chloroplast genome encodes 136 genes, including 88 protein-coding genes, four ribosomal RNA genes, and 40 transfer RNA genes. In the repeat structure analysis, 31 forward and 22 inverted long repeats and 65 simple-sequence repeat loci were detected in the Q. acutissima cp genome. The existence of abundant simple-sequence repeat loci in the genome suggests the potential for future population genetic work. The genome comparison revealed that the LSC region is more divergent than the SSC and IR regions, and there is higher divergence in noncoding regions than in coding regions. The phylogenetic relationships of 25 species inferred that members of the Quercus genus do not form a clade and that Q. acutissima is closely related to Q. variabilis. This study identified the unique characteristics of the Q. acutissima cp genome, which will provide a theoretical basis for species identification and biological research.

Complete chloroplast genome sequence of Fagopyrum dibotrys: genome features, comparative analysis and phylogenetic relationships

Fagopyrum dibotrys, belongs to Polygonaceae family, is one of national key conserved wild plants of China with important medicinal and economic values. Here, the complete chloroplast (cp) genome sequence of F. dibotrys is reported. The cp genome size is 159,919 bp with a typical quadripartite structure and consisting of a pair of inverted repeat regions (30,738 bp) separated by large single copy region (85,134 bp) and small single copy region (13,309 bp). Sequencing analyses indicated that the cp genome encodes 131 genes, including 80 protein-coding genes, 28 tRNA genes and 4 rRNA genes. The genome structure, gene order and codon usage are typical of angiosperm cp genomes. We also identified 48 simple sequence repeats (SSR) loci, fewer of them are distributed in the protein-coding sequences compared to the noncoding regions. Comparison of F. dibotrys cp genome to other Polygonaceae cp genomes indicated the inverted repeats (IRs) and coding regions were more conserved than single copy and noncoding regions, and several variation hotspots were detected. Coding gene sequence divergence analyses indicated that five genes (ndhK, petL rpoC2, ycf1, ycf2) were subject to positive selection. Phylogenetic analysis among 42 species based on cp genomes and 50 protein-coding genes indicated a close relationship between F. dibotrys and F. tataricum. In summary, the complete cp genome sequence of F. dibotrys reported in this study will provide useful plastid genomic resources for population genetics and pave the way for resolving phylogenetic relationships of order Caryophyllales.

RNA-seq highlights parallel and contrasting patterns in the evolution of the nuclear genome of fully mycoheterotrophic plants

BACKGROUND

While photosynthesis is the most notable trait of plants, several lineages of plants (so-called full heterotrophs) have adapted to obtain organic compounds from other sources. The switch to heterotrophy leads to profound changes at the morphological, physiological and genomic levels.

RESULTS

Here, we characterize the transcriptomes of three species representing two lineages of mycoheterotrophic plants: orchids (Epipogium aphyllum and Epipogium roseum) and Ericaceae (Hypopitys monotropa). Comparative analysis is used to highlight the parallelism between distantly related fully heterotrophic plants. In both lineages, we observed genome-wide elimination of nuclear genes that encode proteins related to photosynthesis, while systems associated with protein import to plastids as well as plastid transcription and translation remain active. Genes encoding components of plastid ribosomes that have been lost from the plastid genomes have not been transferred to the nuclear genomes; instead, some of the encoded proteins have been substituted by homologs. The nuclear genes of both Epipogium species accumulated nucleotide substitutions twice as rapidly as their photosynthetic relatives; in contrast, no increase in the substitution rate was observed in H. monotropa.

CONCLUSIONS

Full heterotrophy leads to profound changes in nuclear gene content. The observed increase in the rate of nucleotide substitutions is lineage specific, rather than a universal phenomenon among non-photosynthetic plants.

Plastome organization, genome-based phylogeny and evolution of plastid genes in Podophylloideae (Berberidaceae)

Species of Podophylloideae (Berberidaceae, Ranunculales) are of great pharmacogenetic importance and represent the classic biogeographic disjunction between eastern Asia (EA; 10 ssp.) and eastern North America (ENA; 2 ssp.). However, previous molecular studies of this group suffered from low phylogenetic resolution and/or insufficient marker variability. This study is the first to report whole-plastome sequence data for all 12 species of Podophylloideae (14 individuals) and a close relative, Achlys triphylla. These 15 plastomes proved highly similar in overall size (156,240-157,370 bp), structure, gene order and content, also when compared to other Ranunculales, but also revealed some structural variations caused by the expansion or contraction of the inverted repeats (IRs) into or out of adjacent single-copy regions. Our phylogenomic analysis, based on 63 plastome-derived protein-coding genes (CDS), supported the monophyly of Podophylloideae and its two major genera (EA: Dysosma, EA/ENA: Diphylleia), with Podophyllum peltatum L. (ENA) being more closely related to Diphylleia than to the group's earliest diverging species, Sinopodophyllum hexandrum (EA). Furthermore, within this subfamily/dataset, matK was identified as the fastest evolving gene, which proved to be under positive selection especially in more recently derived, lower-elevation lineages of Dysosma, possibly reflecting an adaptive response to novel environmental (i.e. subtropical compared to higher-elevation/alpine) conditions. Finally, several highly variable noncoding regions were identified in the plastomes of Podophylloideae and Ranunculales. These highly variable loci should be the best choices for future phylogenetic, phylogeographic, and population-level genetic studies. Overall, our results demonstrate the power of plastid phylogenomics to improve phylogenetic resolution, and contribute to a better understanding of plastid gene evolution in Podophylloideae.

The Direct Involvement of Dark-Induced Tic55 Protein in Chlorophyll Catabolism and Its Indirect Role in the MYB108-NAC Signaling Pathway during Leaf Senescence in Arabidopsis thaliana

The chloroplast relies on proteins encoded in the nucleus, synthesized in the cytosol and subsequently transported into chloroplast through the protein complexes Toc and Tic (Translocon at the outer/inner membrane of chloroplasts). A Tic complex member, Tic55, contains a redox-related motif essential for protein import into chloroplasts in peas. However, Tic55 is not crucial for protein import in Arabidopsis. Here, a tic55-II-knockout mutant of Arabidopsis thaliana was characterized for Tic55 localization, its relationship with other translocon proteins, and its association with plant leaf senescence when compared to the wild type. Individually darkened leaves (IDLs) obtained through dark-induced leaf senescence were used to demonstrate chlorophyll breakdown and its relationship with plant senescence in the tic55-II-knockout mutant. The IDLs of the tic55-II-knockout mutant contained higher chlorophyll concentrations than those of the wild type. Our microarray analysis of IDLs during leaf senescence identified seven senescence-associated genes (SAGs) that were downregulated in the tic55-II-knockout mutant: ASP3, APG7, DIN2, DIN11, SAG12, SAG13, and YLS9. Real-time quantitative PCR confirmed the reliability of microarray analysis by showing the same expression patterns with those of the microarray data. Thus, Tic55 functions in dark-induced aging in A. thaliana by indirectly regulating downstream SAGs expression. In addition, the expression of four NAC genes, including ANAC003, ANAC010, ANAC042, and ANAC075 of IDL treated tic55-II-knockout mutant appeared to be downregulated. Yeast one hybrid assay revealed that only ANAC003 promoter region can be bound by MYB108, suggesting that a MYB-NAC regulatory network is involved in dark-stressed senescence.

New Perspectives on Chloroplast Protein Import

Virtually all chloroplasts in extant photosynthetic eukaryotes derive from a single endosymbiotic event that probably occurred more than a billion years ago between a host eukaryotic cell and a cyanobacterium-like ancestor. Many endosymbiont genes were subsequently transferred to the host nuclear genome, concomitant with the establishment of a system for protein transport through the chloroplast double-membrane envelope. Presently, 2,000-3,000 different nucleus-encoded chloroplast proteins must be imported into the chloroplast following their synthesis in the cytosol. The TOC (translocon at the outer envelope membrane of chloroplasts) and TIC (translocon at the inner envelope membrane of chloroplasts) complexes are protein translocation machineries at the outer and inner envelope membranes, respectively, that facilitate this chloroplast protein import with the aid of a TIC-associated ATP-driven import motor. All the essential components of this protein import system seemed to have been identified through biochemical analyses and subsequent genetic studies that initiated in the late 1990s. However, in 2013, the Nakai group reported a novel inner envelope membrane TIC complex, for which a novel ATP-driven import motor associated with this TIC complex is likely to exist. In this mini review, I will summarize these recent discoveries together with new, or reanalyzed, data presented by other groups in recent years. Whereas the precise concurrent view of chloroplast protein import is still a matter of some debate, it is anticipated that the entire TOC/TIC/ATP motor system, including any novel components, will be conclusively established in the next decade. Such findings may lead to an extensively revised view of the evolution and molecular mechanisms of chloroplast protein import.

The Carboxy Terminus of YCF1 Contains a Motif Conserved throughout >500 Myr of Streptophyte Evolution

Plastids evolved from cyanobacteria by endosymbiosis. During the course of evolution, the coding capacity of plastid genomes shrinks due to gene loss or transfer to the nucleus. In the green lineage, however, there were apparent gene gains including that of ycf1. Although its function is still debated, YCF1 has proven to be a useful marker for plastid evolution. YCF1 sequence and predicted structural features unite the plastid genomes of land plants with those of their closest algal relatives, the higher streptophyte algae; YCF1 appears to have undergone pronounced changes during the course of streptophyte algal evolution. Using new data, we show that YCF1 underwent divergent evolution in the common ancestor of higher streptophyte algae and Klebsormidiophycae. This divergence resulted in the origin of an extreme, klebsormidiophycean-specific YCF1 and the higher streptophyte Ste-YCF1. Most importantly, our analysis uncovers a conserved carboxy-terminal sequence stretch within YCF1 that is unique to higher streptophytes and hints at an important, yet unexplored function.

Development of high transferability cpSSR markers for individual identification and genetic investigation in Cupressaceae species

Given the low substitution rate in plastomes, the polymorphic and codominant nature of chloroplast SSRs (cpSSRs) makes them ideal markers, complementing their nuclear counterpart. In Cupressaceae, cpSSRs are mostly paternally inherited, thus, they are useful in mating systems and pollen flow studies. Using e-PCR, 92 SSR loci were identified across six Cupressaceae plastomes, and primers were designed for 26 loci with potential interspecific transferability. The 26 developed cpSSRs were polymorphic in four genera, Platycladus, Sabina, Juniperus, and Cupressus and are suitable for Cupressaceae molecular genetic studies and utilization. We genotyped 192 Platycladus orientalis samples from a core breeding population using 10 of the developed cpSSRs and 10 nuclear SSRs, and these individuals were identified with high confidence. The developed cpSSRs can be used in (1) a marker-assisted breeding scheme, specifically when paternity identification is required, (2) population genetics investigations, and (3) biogeography of Cupressaceae and unraveling the genetic relationships between related species.

Genomics-Informed Insights into Endosymbiotic Organelle Evolution in Photosynthetic Eukaryotes

The conversion of free-living cyanobacteria to photosynthetic organelles of eukaryotic cells through endosymbiosis transformed the biosphere and eventually provided the basis for life on land. Despite the presumable advantage conferred by the acquisition of photoautotrophy through endosymbiosis, only two independent cases of primary endosymbiosis have been documented: one that gave rise to the Archaeplastida, and the other to photosynthetic species of the thecate, filose amoeba Paulinella. Here, we review recent genomics-informed insights into the primary endosymbiotic origins of cyanobacteria-derived organelles. Furthermore, we discuss the preconditions for the evolution of nitrogen-fixing organelles. Recent genomic data on previously undersampled cyanobacterial and protist taxa provide new clues to the origins of the host cell and endosymbiont, and proteomic approaches allow insights into the rearrangement of the endosymbiont proteome during organellogenesis. We conclude that in addition to endosymbiotic gene transfers, horizontal gene acquisitions from a broad variety of prokaryotic taxa were crucial to organelle evolution.

Evolution and Design Principles of the Diverse Chloroplast Transit Peptides

Chloroplasts are present in organisms belonging to the kingdom Plantae. These organelles are thought to have originated from photosynthetic cyanobacteria through endosymbiosis. During endosymbiosis, most cyanobacterial genes were transferred to the host nucleus. Therefore, most chloroplast proteins became encoded in the nuclear genome and must return to the chloroplast after translation. The N-terminal cleavable transit peptide (TP) is necessary and sufficient for the import of nucleus-encoded interior chloroplast proteins. Over the past decade, extensive research on the TP has revealed many important characteristic features of TPs. These studies have also shed light on the question of how the many diverse TPs could have evolved to target specific proteins to the chloroplast. In this review, we summarize the characteristic features of TPs. We also highlight recent advances in our understanding of TP evolution and provide future perspectives about this important research area.

Development of Chloroplast Genomic Resources in Chinese Yam (Dioscorea polystachya)

Chinese yam has been used both as a food and in traditional herbal medicine. Developing more effective genetic markers in this species is necessary to assess its genetic diversity and perform cultivar identification. In this study, new chloroplast genomic resources were developed using whole chloroplast genomes from six genotypes originating from different geographical locations. The Dioscorea polystachya chloroplast genome is a circular molecule consisting of two single-copy regions separated by a pair of inverted repeats. Comparative analyses of six D. polystachya chloroplast genomes revealed 141 single nucleotide polymorphisms (SNPs). Seventy simple sequence repeats (SSRs) were found in the six genotypes, including 24 polymorphic SSRs. Forty-three common indels and five small inversions were detected. Phylogenetic analysis based on the complete chloroplast genome provided the best resolution among the genotypes. Our evaluation of chloroplast genome resources among these genotypes led us to consider the complete chloroplast genome sequence of D. polystachya as a source of reliable and valuable molecular markers for revealing biogeographical structure and the extent of genetic variation in wild populations and for identifying different cultivars.

The complete chloroplast genome sequence of Epipremnum aureum and its comparative analysis among eight Araceae species

Epipremnum aureum is an important foliage plant in the Araceae family. In this study, we have sequenced the complete chloroplast genome of E. aureum by using Illumina Hiseq sequencing platforms. This genome is a double-stranded circular DNA sequence of 164,831 bp that contains 35.8% GC. The two inverted repeats (IRa and IRb; 26,606 bp) are spaced by a small single-copy region (22,868 bp) and a large single-copy region (88,751 bp). The chloroplast genome has 131 (113 unique) functional genes, including 86 (79 unique) protein-coding genes, 37 (30 unique) tRNA genes, and eight (four unique) rRNA genes. Tandem repeats comprise the majority of the 43 long repetitive sequences. In addition, 111 simple sequence repeats are present, with mononucleotides being the most common type and di- and tetranucleotides being infrequent events. Positive selection pressure on rps12 in the E. aureum chloroplast has been demonstrated via synonymous and nonsynonymous substitution rates and selection pressure sites analyses. Ycf15 and infA are pseudogenes in this species. We constructed a Maximum Likelihood phylogenetic tree based on the complete chloroplast genomes of 38 species from 13 families. Those results strongly indicated that E. aureum is positioned as the sister of Colocasia esculenta within the Araceae family. This work may provide information for further study of the molecular phylogenetic relationships within Araceae, as well as molecular markers and breeding novel varieties by chloroplast genetic-transformation of E. aureum in particular.

Molecular Evolution of Chloroplast Genomes of Orchid Species: Insights into Phylogenetic Relationship and Adaptive Evolution

Orchidaceae is the 3rd largest family of angiosperms, an evolved young branch of monocotyledons. This family contains a number of economically-important horticulture and flowering plants. However, the limited availability of genomic information largely hindered the study of molecular evolution and phylogeny of Orchidaceae. In this study, we determined the evolutionary characteristics of whole chloroplast (cp) genomes and the phylogenetic relationships of the family Orchidaceae. We firstly characterized the cp genomes of four orchid species: Cremastra appendiculata, Calanthe davidii, Epipactis mairei, and Platanthera japonica. The size of the chloroplast genome ranged from 153,629 bp (C. davidi) to 160,427 bp (E. mairei). The gene order, GC content, and gene compositions are similar to those of other previously-reported angiosperms. We identified that the genes of ndhC, ndhI, and ndhK were lost in C. appendiculata, in that the ndh I gene was lost in P. japonica and E. mairei. In addition, the four types of repeats (forward, palindromic, reverse, and complement repeats) were examined in orchid species. E. mairei had the highest number of repeats (81), while C. davidii had the lowest number (57). The total number of Simple Sequence Repeats is at least 50 in C. davidii, and, at most, 78 in P. japonica. Interestingly, we identified 16 genes with positive selection sites (the psbH, petD, petL, rpl22, rpl32, rpoC1, rpoC2, rps12, rps15, rps16, accD, ccsA, rbcL, ycf1, ycf2, and ycf4 genes), which might play an important role in the orchid species’ adaptation to diverse environments. Additionally, 11 mutational hotspot regions were determined, including five non-coding regions (ndhB intron, ccsA-ndhD, rpl33-rps18, ndhE-ndhG, and ndhF-rpl32) and six coding regions (rps16, ndhC, rpl32, ndhI, ndhK, and ndhF). The phylogenetic analysis based on whole cp genomes showed that C. appendiculata was closely related to C. striata var. vreelandii, while C. davidii and C. triplicate formed a small monophyletic evolutionary clade with a high bootstrap support. In addition, five subfamilies of Orchidaceae, Apostasioideae, Cypripedioideae, Epidendroideae, Orchidoideae, and Vanilloideae, formed a nested evolutionary relationship in the phylogenetic tree. These results provide important insights into the adaptive evolution and phylogeny of Orchidaceae.

The Linum usitatissimum L. plastome reveals atypical structural evolution, new editing sites, and the phylogenetic position of Linaceae within Malpighiales

The plastome of Linum usitatissimum was completely sequenced allowing analyses of evolution of genome structure, RNA editing sites, molecular markers, and indicating the position of Linaceae within Malpighiales. Flax (Linum usitatissimum L.) is an economically important crop used as food, feed, and industrial feedstock. It belongs to the Linaceae family, which is noted by high morphological and ecological diversity. Here, we reported the complete sequence of flax plastome, the first species within Linaceae family to have the plastome sequenced, assembled and characterized in detail. The plastome of flax is a circular DNA molecule of 156,721 bp with a typical quadripartite structure including two IRs of 31,990 bp separating the LSC of 81,767 bp and the SSC of 10,974 bp. It shows two expansion events from IRB to LSC and from IRB to SSC, and a contraction event in the IRA-LSC junction, which changed significantly the size and the gene content of LSC, SSC and IRs. We identified 109 unique genes and 2 pseudogenes (rpl23 and ndhF). The plastome lost the conserved introns of clpP gene and the complete sequence of rps16 gene. The clpP, ycf1, and ycf2 genes show high nucleotide and aminoacid divergence, but they still possibly retain the functionality. Moreover, we also identified 176 SSRs, 20 tandem repeats, and 39 dispersed repeats. We predicted in 18 genes a total of 53 RNA editing sites of which 32 were not found before in other species. The phylogenetic inference based on 63 plastid protein-coding genes of 38 taxa supports three major clades within Malpighiales order. One of these clades has flax (Linaceae) sister to Chrysobalanaceae family, differing from earlier studies that included Linaceae into the euphorbioid clade.

Complete Chloroplast Genome Sequence and Phylogenetic Analysis of Paeonia ostii

Paeonia ostii, a common oil-tree peony, is important ornamentally and medicinally. However, there are few studies on the chloroplast genome of Paeonia ostii. We sequenced and analyzed the complete chloroplast genome of P. ostii. The size of the P. ostii chloroplast genome is 152,153 bp, including a large single-copy region (85,373 bp), a small single-copy region (17,054 bp), and a pair of inverted repeats regions (24,863 bp). The P. ostii chloroplast genome encodes 111 genes, including 77 protein-coding genes, four ribosomal RNA genes, and 30 transfer RNA genes. The genome contains forward repeats (22), palindromic repeats (28), and tandem repeats (24). The presence of rich simple-sequence repeat loci in the genome provides opportunities for future population genetics work for breeding new varieties. A phylogenetic analysis showed that P. ostii is more closely related to Paeonia delavayi and Paeonialudlowii than to Paeoniaobovata and Paeoniaveitchii. The results of this study provide an assembly of the whole chloroplast genome of P. ostii, which may be useful for future breeding and further biological discoveries. It will provide a theoretical basis for the improvement of peony yield and the determination of phylogenetic status.

A novel chloroplast gene reported for flagellate plants

PREMISE OF THE STUDY

Gene space in plant plastid genomes is well characterized and annotated, yet we discovered an unrecognized open reading frame (ORF) in the fern lineage that is conserved across flagellate plants.

METHODS

We initially detected a putative uncharacterized ORF by the existence of a highly conserved region between rps16 and matK in a series of matK alignments of leptosporangiate ferns. We mined available plastid genomes for this ORF, which we now refer to as ycf94, to infer evolutionary selection pressures and assist in functional prediction. To further examine the transcription of ycf94, we assembled the plastid genome and sequenced the transcriptome of the leptosporangiate fern Adiantum shastense Huiet & A.R. Sm.

KEY RESULTS

The ycf94 predicted protein has a distinct transmembrane domain but with no sequence homology to other proteins with known function. The nonsynonymous/synonymous substitution rate ratio of ycf94 is on par with other fern plastid protein-encoding genes, and additional homologs can be found in a few lycophyte, moss, hornwort, and liverwort plastid genomes. Homologs of ycf94 were not found in seed plants. In addition, we report a high level of RNA editing for ycf94 transcripts-a hallmark of protein-coding genes in fern plastomes.

CONCLUSIONS

The degree of sequence conservation, together with the presence of a distinct transmembrane domain and RNA-editing sites, suggests that ycf94 is a protein-coding gene of functional significance in ferns and, potentially, bryophytes and lycophytes. However, the origin and exact function of this gene require further investigation.

The integration of chloroplast protein targeting with plant developmental and stress responses

The plastids, including chloroplasts, are a group of interrelated organelles that confer photoautotrophic growth and the unique metabolic capabilities that are characteristic of plant systems. Plastid biogenesis relies on the expression, import, and assembly of thousands of nuclear encoded preproteins. Plastid proteomes undergo rapid remodeling in response to developmental and environmental signals to generate functionally distinct plastid types in specific cells and tissues. In this review, we will highlight the central role of the plastid protein import system in regulating and coordinating the import of functionally related sets of preproteins that are required for plastid-type transitions and maintenance.

High-throughput m6A-seq reveals RNA m6A methylation patterns in the chloroplast and mitochondria transcriptomes of Arabidopsis thaliana

This study is the first to comprehensively characterize m6A patterns in the Arabidopsis chloroplast and mitochondria transcriptomes based on our open accessible data deposited in NCBI's Gene Expression Omnibus with GEO Series accession number of GSE72706. Over 86% of the transcripts were methylated by m6A in the two organelles. Over 550 and 350 m6A sites were mapped, with ~5.6 to ~5.8 and ~4.6 to ~4.9 m6A sites per transcript, to the chloroplast and mitochondria genome, respectively. The overall m6A methylation extent in the two organelles was greatly higher than that in the nucleus. The m6A motif sequences in the transcriptome of two organelles were similar to the nuclear motifs, suggesting that selection of the m6A motifs for RNA methylation was conserved between the nucleus and organelle transcriptomes. The m6A patterns of rRNAs and tRNAs in the organelle were similar to those in the nucleus. However, the m6A patterns in coding RNAs were distinct between the nucleus and the organelle, suggesting that that regulation of the m6A methylation patterns may be different between the nuclei and the organelles. The extensively methylated transcripts in the two organelles were mainly associated with rRNA, ribosomal proteins, photosystem reaction proteins, tRNA, NADH dehydrogenase and redox. On average, 64% and 79% of the transcripts in the two organelles showed differential m6A methylation across three organs of the leaves, flowers and roots. The m6A methylation extent in the chloroplast was higher than that in the mitochondria. This study provides deep insights into the m6A methylation topology and differentiation in the plant organelle transcriptomes.

Interspecific Plastome Recombination Reflects Ancient Reticulate Evolution in Picea (Pinaceae)

Plastid sequences are a cornerstone in plant systematic studies and key aspects of their evolution, such as uniparental inheritance and absent recombination, are often treated as axioms. While exceptions to these assumptions can profoundly influence evolutionary inference, detecting them can require extensive sampling, abundant sequence data, and detailed testing. Using advancements in high-throughput sequencing, we analyzed the whole plastomes of 65 accessions of Picea, a genus of ∼35 coniferous forest tree species, to test for deviations from canonical plastome evolution. Using complementary hypothesis and data-driven tests, we found evidence for chimeric plastomes generated by interspecific hybridization and recombination in the clade comprising Norway spruce (P. abies) and 10 other species. Support for interspecific recombination remained after controlling for sequence saturation, positive selection, and potential alignment artifacts. These results reconcile previous conflicting plastid-based phylogenies and strengthen the mounting evidence of reticulate evolution in Picea. Given the relatively high frequency of hybridization and biparental plastid inheritance in plants, we suggest interspecific plastome recombination may be more widespread than currently appreciated and could underlie reported cases of discordant plastid phylogenies.

The evolutionary fate of the chloroplast and nuclear rps16 genes as revealed through the sequencing and comparative analyses of four novel legume chloroplast genomes from Lupinus

The Fabaceae family is considered as a model system for understanding chloroplast genome evolution due to the presence of extensive structural rearrangements, gene losses and localized hypermutable regions. Here, we provide sequences of four chloroplast genomes from the Lupinus genus, belonging to the underinvestigated Genistoid clade. Notably, we found in Lupinus species the functional loss of the essential rps16 gene, which was most likely replaced by the nuclear rps16 gene that encodes chloroplast and mitochondrion targeted RPS16 proteins. To study the evolutionary fate of the rps16 gene, we explored all available plant chloroplast, mitochondrial and nuclear genomes. Whereas no plant mitochondrial genomes carry an rps16 gene, many plants still have a functional nuclear and chloroplast rps16 gene. Ka/Ks ratios revealed that both chloroplast and nuclear rps16 copies were under purifying selection. However, due to the dual targeting of the nuclear rps16 gene product and the absence of a mitochondrial copy, the chloroplast gene may be lost. We also performed comparative analyses of lupine plastomes (SNPs, indels and repeat elements), identified the most variable regions and examined their phylogenetic utility. The markers identified here will help to reveal the evolutionary history of lupines, Genistoids and closely related clades.

Predicting Amyloidogenic Proteins in the Proteomes of Plants

Amyloids are protein fibrils with characteristic spatial structure. Though amyloids were long perceived to be pathogens that cause dozens of incurable pathologies in humans and mammals, it is currently clear that amyloids also represent a functionally important form of protein structure implicated in a variety of biological processes in organisms ranging from archaea and bacteria to fungi and animals. Despite their social significance, plants remain the most poorly studied group of organisms in the field of amyloid biology. To date, amyloid properties have only been demonstrated in vitro or in heterologous systems for a small number of plant proteins. Here, for the first time, we performed a comprehensive analysis of the distribution of potentially amyloidogenic proteins in the proteomes of approximately 70 species of land plants using the Waltz and SARP (Sequence Analysis based on the Ranking of Probabilities) bioinformatic algorithms. We analyzed more than 2.9 million protein sequences and found that potentially amyloidogenic proteins are abundant in plant proteomes. We found that such proteins are overrepresented among membrane as well as DNA- and RNA-binding proteins of plants. Moreover, seed storage and defense proteins of most plant species are rich in amyloidogenic regions. Taken together, our data demonstrate the diversity of potentially amyloidogenic proteins in plant proteomes and suggest biological processes where formation of amyloids might be functionally important.

On the brink: the highly reduced plastomes of nonphotosynthetic Ericaceae

Ericaceae (the heather family) is a large and diverse group of plants that forms elaborate symbiotic relationships with mycorrhizal fungi, and includes several nonphotosynthetic lineages. Using an extensive sample of fully mycoheterotrophic (MH) species, we explored inter- and intraspecific variation as well as selective constraints acting on the plastomes of these unusual plants. The plastomes of seven MH genera were analysed in a phylogenetic context with two geographically disparate individuals sequenced for Allotropa, Monotropa, and Pityopus. The plastomes of nonphotosynthetic Ericaceae are highly reduced in size (c. 33-41 kbp) and content, having lost all photosynthesis-related genes, and are reduced to encoding housekeeping genes as well as a protease subunit (clpP)-like and acetyl-CoA carboxylase subunit D (accD)-like open reading frames. Despite an increase in the rate of their nucleotide substitutions, the remaining protein-coding genes are typically under purifying selection in full MHs. We also identified ribosomal proteins under relaxed or neutral selection. These plastomes also exhibit striking structural rearrangements. Intraspecific variation within MH Ericaceae ranges from a few differences (Allotropa) to extensive population divergences (Monotropa, Hypopitys), which indicates that cryptic speciation may be occurring in several lineages. The pattern of gene loss within fully MH Ericaceae plastomes suggests an advanced state of degradation.

Stable megadalton TOC-TIC supercomplexes as major mediators of protein import into chloroplasts

Preproteins are believed to be imported into chloroplasts through membrane contact sites where the translocon complexes of the outer (TOC) and inner (TIC) envelope membranes are assembled together. However, a single TOC-TIC supercomplex containing preproteins undergoing active import has not yet been directly observed. We optimized the blue native polyacrylamide gel electrophoresis (PAGE) (BN-PAGE) system to detect and resolve megadalton (MD)-sized complexes. Using this optimized system, the outer-membrane channel Toc75 from pea chloroplasts was found in at least two complexes: the 880-kD TOC complex and a previously undetected 1-MD complex. Two-dimensional BN-PAGE immunoblots further showed that Toc75, Toc159, Toc34, Tic20, Tic56 and Tic110 were all located in the 880-kD to 1.3-MD region. During active preprotein import, preproteins were transported mostly through the 1-MD complex and a smaller amount of preproteins was also detected in a complex of 1.25 MD. Antibody-shift assays showed that the 1-MD complex is a TOC-TIC supercomplex containing at least Toc75, Toc159, Toc34 and Tic110. Results from crosslinking and import with Arabidopsis chloroplasts suggest that the 1.25-MD complex is also a supercomplex. Our data provide direct evidence supporting that chloroplast preproteins are imported through TOC-TIC supercomplexes, and also provide the first size estimation of these supercomplexes. Furthermore, unlike in mitochondria where translocon supercomplexes are only transiently assembled during preprotein import, in chloroplasts at least some of the supercomplexes are preassembled stable structures.

Complete chloroplast genome sequences contribute to plant species delimitation: A case study of the Anemopaegma species complex

PREMISE OF THE STUDY

Bignoniaceae is an important component of neotropical forests and a model for evolutionary and biogeographical studies. A previous combination of molecular markers and morphological traits improved the phylogeny of the group. Here we demonstrate the value of next-generation sequencing (NGS) to assemble the chloroplast genome of eight Anemopaegma species and solve taxonomic problems.

METHODS

Three NGS platforms were used to sequence total DNA of Anemopaegma species. After genome assembly and annotation, we compared chloroplast genomes within Anemopaegma, with other Lamiales species, and the evolutionary rates of protein-coding genes using Tanaecium tetragonolobum as the outgroup. Phylogenetic analyses of Anemopaegma with different data sets were performed.

KEY RESULTS

Chloroplast genomes of Anemopaegma species ranged from 167,413 bp in A. foetidum to 168,987 bp in A. acutifolium ("typical" form). They exhibited a characteristic quadripartite structure with a large single-copy region (75,070-75,761 bp), a small single-copy region (12,766-12,817 bp) and a pair of inverted repeat regions (IRs) (39,480-40,481) encoding an identical set of 112 genes. An inversion of a fragment with ca. 8 kb, located in the IRs and containing the genes trnI-AAU, ycf2, and trnL-CAA, was observed in these chloroplast genomes when compared with those of other Lamiales.

CONCLUSIONS

Anemopaegma species have the largest genomes within the Lamiales possibly due to the large amount of repetitive sequences and IR expansion. Variation was higher in coding regions than in noncoding regions, and some genes were identified as markers for differentiation between species. The use of the entire chloroplast genome gave better phylogenetic resolution of the taxonomic groups. We found that two forms of A. acutifolium result from different maternal lineages.

Chloroplast genomes of Arabidopsis halleri ssp. gemmifera and Arabidopsis lyrata ssp. petraea: Structures and comparative analysis

We investigated the complete chloroplast (cp) genomes of non-model Arabidopsis halleri ssp. gemmifera and Arabidopsis lyrata ssp. petraea using Illumina paired-end sequencing to understand their genetic organization and structure. Detailed bioinformatics analysis revealed genome sizes of both subspecies ranging between 154.4~154.5 kbp, with a large single-copy region (84,197~84,158 bp), a small single-copy region (17,738~17,813 bp) and pair of inverted repeats (IRa/IRb; 26,264~26,259 bp). Both cp genomes encode 130 genes, including 85 protein-coding genes, eight ribosomal RNA genes and 37 transfer RNA genes. Whole cp genome comparison of A. halleri ssp. gemmifera and A. lyrata ssp. petraea, along with ten other Arabidopsis species, showed an overall high degree of sequence similarity, with divergence among some intergenic spacers. The location and distribution of repeat sequences were determined, and sequence divergences of shared genes were calculated among related species. Comparative phylogenetic analysis of the entire genomic data set and 70 shared genes between both cp genomes confirmed the previous phylogeny and generated phylogenetic trees with the same topologies. The sister species of A. halleri ssp. gemmifera is A. umezawana, whereas the closest relative of A. lyrata spp. petraea is A. arenicola.

Progress and challenges of engineering a biophysical CO2-concentrating mechanism into higher plants

Growth and productivity in important crop plants is limited by the inefficiencies of the C3 photosynthetic pathway. Introducing CO2-concentrating mechanisms (CCMs) into C3 plants could overcome these limitations and lead to increased yields. Many unicellular microautotrophs, such as cyanobacteria and green algae, possess highly efficient biophysical CCMs that increase CO2 concentrations around the primary carboxylase enzyme, Rubisco, to enhance CO2 assimilation rates. Algal and cyanobacterial CCMs utilize distinct molecular components, but share several functional commonalities. Here we outline the recent progress and current challenges of engineering biophysical CCMs into C3 plants. We review the predicted requirements for a functional biophysical CCM based on current knowledge of cyanobacterial and algal CCMs, the molecular engineering tools and research pipelines required to translate our theoretical knowledge into practice, and the current challenges to achieving these goals.

The novel chloroplast outer membrane kinase KOC1 is a required component of the plastid protein import machinery

The biogenesis and maintenance of cell organelles such as mitochondria and chloroplasts require the import of many proteins from the cytosol, a process that is controlled by phosphorylation. In the case of chloroplasts, the import of hundreds of different proteins depends on translocons at the outer and inner chloroplast membrane (TOC and TIC, respectively) complexes. The essential protein TOC159 functions thereby as an import receptor. It has an N-terminal acidic (A-) domain that extends into the cytosol, controls receptor specificity, and is highly phosphorylated in vivo However, kinases that phosphorylate the TOC159 A-domain to enable protein import have remained elusive. Here, using co-purification with TOC159 from Arabidopsis, we discovered a novel component of the chloroplast import machinery, the regulatory kinase at the outer chloroplast membrane 1 (KOC1). We found that KOC1 is an integral membrane protein facing the cytosol and stably associates with TOC. Moreover, KOC1 phosphorylated the A-domain of TOC159 in vitro, and in mutant koc1 chloroplasts, preprotein import efficiency was diminished. koc1 Arabidopsis seedlings had reduced survival rates after transfer from the dark to the light in which protein import into plastids is required to rapidly complete chloroplast biogenesis. In summary, our data indicate that KOC1 is a functional component of the TOC machinery that phosphorylates import receptors, supports preprotein import, and contributes to efficient chloroplast biogenesis.

Responses of the picoprasinophyte Micromonas commoda to light and ultraviolet stress

Micromonas is a unicellular marine green alga that thrives from tropical to polar ecosystems. We investigated the growth and cellular characteristics of acclimated mid-exponential phase Micromonas commoda RCC299 over multiple light levels and over the diel cycle (14:10 hour light:dark). We also exposed the light:dark acclimated M. commoda to experimental shifts from moderate to high light (HL), and to HL plus ultraviolet radiation (HL+UV), 4.5 hours into the light period. Cellular responses of this prasinophyte were quantified by flow cytometry and changes in gene expression by qPCR and RNA-seq. While proxies for chlorophyll a content and cell size exhibited similar diel variations in HL and controls, with progressive increases during day and decreases at night, both parameters sharply decreased after the HL+UV shift. Two distinct transcriptional responses were observed among chloroplast genes in the light shift experiments: i) expression of transcription and translation-related genes decreased over the time course, and this transition occurred earlier in treatments than controls; ii) expression of several photosystem I and II genes increased in HL relative to controls, as did the growth rate within the same diel period. However, expression of these genes decreased in HL+UV, likely as a photoprotective mechanism. RNA-seq also revealed two genes in the chloroplast genome, ycf2-like and ycf1-like, that had not previously been reported. The latter encodes the second largest chloroplast protein in Micromonas and has weak homology to plant Ycf1, an essential component of the plant protein translocon. Analysis of several nuclear genes showed that the expression of LHCSR2, which is involved in non-photochemical quenching, and five light-harvesting-like genes, increased 30 to >50-fold in HL+UV, but was largely unchanged in HL and controls. Under HL alone, a gene encoding a novel nitrite reductase fusion protein (NIRFU) increased, possibly reflecting enhanced N-assimilation under the 625 μmol photons m^-2 s^-1 supplied in the HL treatment. NIRFU’s domain structure suggests it may have more efficient electron transfer than plant NIR proteins. Our analyses indicate that Micromonas can readily respond to abrupt environmental changes, such that strong photoinhibition was provoked by combined exposure to HL and UV, but a ca. 6-fold increase in light was stimulatory.

Protein import-independent functions of Tic56, a component of the 1-MDa translocase at the inner chloroplast envelope membrane

Tic56 is an essential subunit of a 1-MDa protein complex at the inner chloroplast envelope membrane that comprises Tic100, Tic214 and the protein conducting channel protein Tic20-I. The complex was characterized as the "general protein import translocase" because mutants in either of its subunits have a severe growth phenotype and fail to assemble a photosynthetic machinery. In a recent publication we show that the albino phenotype of tic56-1 mutants results at least in part from a defect in ribosome assembly and a deficiency in plastid translation. We furthermore could not detect any impairment of protein import activity with plastids from tic56-3 mutants, despite a lack of full-length Tic56 and a decreased abundance of other 1-MDa complex subunits. These findings suggest that the 1-MDa complex consists of subunits that have functions other than protein import.

Surveying the Oligomeric State of Arabidopsis thaliana Chloroplasts

Blue native-PAGE (BN-PAGE) resolves protein complexes in their native state. When combined with immunoblotting, it can be used to identify the presence of high molecular weight complexes at high resolution for any protein, given a suitable antibody. To identify proteins in high molecular weight complexes on a large scale and to bypass the requirement for specific antibodies, we applied a tandem mass spectrometry (MS/MS) approach to BN-PAGE-resolved chloroplasts. Fractionation of the gel into six bands allowed identification and label-free quantification of 1000 chloroplast proteins with native molecular weight separation. Significantly, this approach achieves a depth of identification comparable with traditional shotgun proteomic analyses of chloroplasts, indicating much of the known chloroplast proteome is amenable to MS/MS identification under our fractionation scheme. By limiting the number of fractionation bands to six, we facilitate scaled-up comparative analyses, as we demonstrate with the reticulata chloroplast mutant displaying a reticulated leaf phenotype. Our comparative proteomics approach identified a candidate interacting protein of RETICULATA as well as effects on lipid remodeling proteins, amino acid metabolic enzymes, and plastid division machinery. We additionally highlight selected proteins from each sub-compartment of the chloroplast that provide novel insight on known or hypothesized protein complexes to further illustrate the utility of this approach. Our results demonstrate the high sensitivity and reproducibility of this technique, which is anticipated to be widely adaptable to other sub-cellular compartments.