Regulation of gene expression in chloroplasts of higher plants

Identification of a highly efficient chloroplast-targeting peptide for plastid engineering

Plastids are pivotal target organelles for comprehensively enhancing photosynthetic and metabolic traits in plants via plastid engineering. Plastidial proteins predominantly originate in the nucleus and must traverse membrane-bound multiprotein translocons to access these organelles. This import process is meticulously regulated by chloroplast-targeting peptides (cTPs). Whereas many cTPs have been employed to guide recombinantly expressed functional proteins to chloroplasts, there is a critical need for more efficient cTPs. Here, we performed a comprehensive exploration and comparative assessment of an advanced suite of cTPs exhibiting superior targeting capabilities. We employed a multifaceted approach encompassing computational prediction, in planta expression, fluorescence tracking, and in vitro chloroplast import studies to identify and analyze 88 cTPs associated with Arabidopsis thaliana mutants with phenotypes linked to chloroplast function. These polypeptides exhibited distinct abilities to transport green fluorescent protein (GFP) to various compartments within leaf cells, particularly chloroplasts. A highly efficient cTP derived from Arabidopsis plastid ribosomal protein L35 (At2g24090) displayed remarkable effectiveness in chloroplast localization. This cTP facilitated the activities of chloroplast-targeted RNA-processing proteins and metabolic enzymes within plastids. This cTP could serve as an ideal transit peptide for precisely targeting biomolecules to plastids, leading to advancements in plastid engineering.

Comparative Phylogenetic Analysis of Ancient Korean Tea "Hadong Cheon-Nyeon Cha (Camellia sinensis var. sinensis)" Using Complete Chloroplast Genome Sequences

Wild teas are valuable genetic resources for studying evolution and breeding. Here, we report the complete chloroplast genome of the ancient Korean tea 'Hadong Cheon-nyeon Cha' (C. sinensis var. sinensis), which is known as the oldest tea tree in Korea. This study determined seven Camellia sinensis var. sinenesis, including Hadong Cheon-nyeon Cha (HCNC) chloroplast genome sequences, using Illumina sequencing technology via de novo assembly. The chloroplast genome sizes ranged from 157,019 to 157,114 bp and were organized into quadripartite regions with the typical chloroplast genomes. Further, differences in SNPs and InDels were detected across the seven chloroplast genomes through variance analysis. Principal component and phylogenetic analysis suggested that regional constraints, rather than functional constraints, strongly affected the sequence evolution of the cp genomes in this study. These genomic resources provide evolutionary insight into Korean tea plant cultivars and lay the foundation for a better understanding of the ancient Korean tea plant HCNC.

Comparative chloroplast genomes provided insights into the evolution and species identification on the Datureae plants

Generally, chloroplast genomes of angiosperms are always highly conserved but carry a certain number of variation among species. In this study, chloroplast genomes of 13 species from Datureae tribe that are of importance both in ornamental gardening and medicinal usage were studied. In addition, seven chloroplast genomes from Datureae together with two from Solanaceae species retrieved from the National Center for Biotechnology Information (NCBI) were integrated into this study. The chloroplast genomes ranged in size from 154,686 to 155,979 and from 155,497 to 155,919 bp for species of Datura and Brugmansia, respectively. As to Datura and Brugmansia, a total of 128 and 132 genes were identified, in which 83 and 87 protein coding genes were identified, respectively; Furthermore, 37 tRNA genes and 8 rRNA genes were both identified in Datura and Brugmansia. Repeats analysis indicated that the number and type varied among species for Simple sequence repeat (SSR), long repeats, and tandem repeats ranged in number from 53 to 59, 98 to 99, and 22 to 30, respectively. Phylogenetic analysis based on the plastid genomes supported the monophyletic relationship among Datura and Brugmansia and Trompettia, and a refined phylogenic relationships among each individual was resolved. In addition, a species-specific marker was designed based on variation spot that resulted from a comparative analysis of chloroplast genomes and verified as effective maker for identification of D. stramonium and D. stramonium var. inermis. Interestingly, we found that 31 genes were likely to be under positive selection, including genes encoding ATP protein subunits, photosystem protein subunit, ribosome protein subunits, NAD(P)H dehydrogenase complex subunits, and clpP, petB, rbcL, rpoCl, ycf4, and cemA genes. These genes may function as key roles in the adaption to diverse environment during evolution. The diversification of Datureae members was dated back to the late Oligocene periods. These chloroplast genomes are useful genetic resources for taxonomy, phylogeny, and evolution for Datureae.

Chloroplasts evolved an additional layer of translational regulation based on non-AUG start codons for proteins with different turnover rates

Chloroplasts have evolved from photosynthetic cyanobacteria-like progenitors through endosymbiosis. The chloroplasts of present-day land plants have their own transcription and translation systems that show several similarities with prokaryotic organisms. A remarkable feature of the chloroplast translation system is the use of non-AUG start codons in the protein synthesis of certain genes that are evolutionarily conserved from Algae to angiosperms. However, the biological significance of such use of non-AUG codons is not fully understood. The present study was undertaken to unravel the significance of non-AUG start codons in vivo using the chloroplast genetic engineering approach. For this purpose, stable transplastomic tobacco plants expressing a reporter gene i.e. uidA (GUS) under four different start codons (AUG/UUG/GUG/CUG) were generated and β-glucuronidase (GUS) expression was compared. To investigate further the role of promoter sequences proximal to the start codon, uidA was expressed under two different chloroplast gene promoters psbA and psbC that use AUG and a non-AUG (GUG) start codons, respectively, and also showed significant differences in the DNA sequence surrounding the start codon. Further, to delineate the role of RNA editing that creates AUG start codon by editing non-AUG codons, if any, which is another important feature of the chloroplast transcription and translation system, transcripts were sequenced. In addition, a proteomic approach was used to identify the translation initiation site(s) of GUS and the N-terminal amino acid encoded when expressed under different non-AUG start codons. The results showed that chloroplasts use non-AUG start codons in combination with the translation initiation site as an additional layer of gene regulation to over-express proteins that are required at high levels due to their high rates of turnover.

An Overview of Pentatricopeptide Repeat (PPR) Proteins in the Moss Physcomitrium patens and Their Role in Organellar Gene Expression

Pentatricopeptide repeat (PPR) proteins are one type of helical repeat protein that are widespread in eukaryotes. In particular, there are several hundred PPR members in flowering plants. The majority of PPR proteins are localized in the plastids and mitochondria, where they play a crucial role in various aspects of RNA metabolism at the post-transcriptional and translational steps during gene expression. Among the early land plants, the moss Physcomitrium (formerly Physcomitrella) patens has at least 107 PPR protein-encoding genes, but most of their functions remain unclear. To elucidate the functions of PPR proteins, a reverse-genetics approach has been applied to P. patens. To date, the molecular functions of 22 PPR proteins were identified as essential factors required for either mRNA processing and stabilization, RNA splicing, or RNA editing. This review examines the P. patens PPR gene family and their current functional characterization. Similarities and a diversity of functions of PPR proteins between P. patens and flowering plants and their roles in the post-transcriptional regulation of organellar gene expression are discussed.

Physcomitrium patens pentatricopeptide repeat protein PpPPR_32 is involved in the accumulation of psaC mRNA encoding the iron sulfur protein of photosystem I

Pentatricopeptide repeat (PPR) proteins are involved in RNA metabolism and also play a role in posttranscriptional regulation during plant organellar gene expression. Although a hundred of PPR proteins exist in the moss Physcomitrium patens, their functions are not fully understood. Here, we report the function of P-class PPR protein PpPPR_32 in P. patens. A transient expression assay using green fluorescent protein demonstrated that the N-terminal region of PpPPR_32 functions as a chloroplast-targeting transit peptide, indicating that PpPPR_32 is localized in chloroplasts. PpPPR_32 knockout (KO) mutants grew autotrophically but with reduced protonema growth and the poor formation of photosystem I (PSI) complexes. Quantitative real-time reverse transcription-polymerase chain reaction and RNA gel blot hybridization analyses revealed a significant reduction in the transcript level of the psaC gene encoding the iron sulfur protein of PSI but no alteration to the transcript levels of other PSI genes. This suggests that PpPPR_32 is specifically involved in the expression level of the psaC gene. Our results indicate that PpPPR_32 is essential for the accumulation of psaC transcript and PSI complexes.

An overview of RNA splicing and functioning of splicing factors in land plant chloroplasts

RNA splicing refers to a process by which introns of a pre-mRNA are excised and the exons at both ends are joined together. Chloroplast introns are inherently self-splicing ribozymes, but over time, they have lost self-splicing ability due to the degeneration of intronic elements. Thus, the splicing of chloroplast introns relies heavily on nuclear-encoded splicing factors, which belong to diverse protein families. Different splicing factors and their shared intron targets are supposed to form ribonucleoprotein particles (RNPs) to facilitate intron splicing. As characterized in a previous review, around 14 chloroplast intron splicing factors were identified until 2010. However, only a few genetic and biochemical evidence has shown that these splicing factors are required for the splicing of one or several introns. The roles of splicing factors are generally believed to facilitate intron folding; however, the precise role of each protein in RNA splicing remains ambiguous. This may be because the precise binding site of most of these splicing factors remains unexplored. In the last decade, several new splicing factors have been identified. Also, several splicing factors were found to bind to specific sequences within introns, which enhanced the understanding of splicing factors. Here, we summarize recent progress on the splicing factors in land plant chloroplasts and discuss their possible roles in chloroplast RNA splicing based on previous studies.

Chloroplast genomes of two Pueraria DC. species: sequencing, comparative analysis and molecular marker development

Puerariae lobatae radix (Ge‐Gen in Chinese) and Puerariae thomsonii radix (Fen‐Ge) are widely used as medicine and health products, particularly in Chinese medicine. Puerarin and daidzein are the primary bioactive compounds in Puerariae radix. These isoflavones have been used to treat cardiovascular and cerebrovascular diseases, hypertension, diabetes, and osteoporosis. The content of puerarin in Ge‐Gen is about six times higher than that in Fen‐Ge, so its use has a higher pharmacological effect. It is therefore of great importance to effectively distinguish between these two species. However, because their basal plants, P. lobata (Willd.) Ohwi and P. thomsonii Benth., possess an extremely similar appearance, and detecting the level of chemical constituents is just a rough distinction, it is necessary to develop more efficient identification approaches. Here the complete chloroplast genomes of P. lobata and P. thomsonii were deciphered, including sequencing, assembly, comparative analysis, and molecular marker development. The results showed that they are 153,393 and 153,442 bp in length, respectively; both contain 124 annotated genes, including eight encoding rRNA, 29 encoding tRNA, and 87 encoding proteins. Phylogenetic analysis showed that they form a clade, indicating that they originate from the same ancestor. After obtaining 10 intergenic/intronic regions with a genetic distance greater than 0.5 cm, primers were designed to amplify regions of high variability in P. lobata and P. thomsonii. Finally, a 60‐bp differential base fragment, located in the intron of rpl16, was developed as a molecular marker to efficiently distinguish between these two species.

Moss PPR-SMR protein PpPPR_64 influences the expression of a psaA-psaB-rps14 gene cluster and processing of the 23S-4.5S rRNA precursor in chloroplasts

Moss PPR-SMR protein PpPPR_64 is a pTAC2 homolog but is functionally distinct from pTAC2. PpPPR_64 is required for psaA gene expression and its function may have evolved in mosses. The pentatricopeptide repeat (PPR) proteins are key regulatory factors responsible for the control of plant organellar gene expression. A small subset of PPR proteins possess a C-terminal small MutS-related (SMR) domain and have diverse roles in plant organellar biogenesis. However, the function of PPR-SMR proteins is not fully understood. Here, we report the function of PPR-SMR protein PpPPR_64 in the moss Physcomitrium patens. Phylogenetic analysis indicated that PpPPR_64 belongs to the same clade as the Arabidopsis PPR-SMR protein pTAC2. PpPPR_64 knockout (KO) mutants grew autotrophically but with reduced protonemata growth and the poor formation of photosystems' antenna complexes. Quantitative reverse transcription-polymerase chain reaction and RNA gel blot hybridization analyses revealed a significant reduction in transcript levels of the psaA-psaB-rps14 gene cluster but no alteration to transcript levels of most photosynthesis- and non-photosynthesis-related genes. In addition, RNA processing of 23S-4.5S rRNA precursor was impaired in the PpPPR_64 KO mutants. This suggests that PpPPR_64 is specifically involved in the expression level of the psaA-psaB-rps14 gene and in processing of the 23S-4.5S rRNA precursor. Our results indicate that PpPPR_64 is functionally distinct from pTAC2 and is a novel PPR-SMR protein required for proper chloroplast biogenesis in P. patens.

Oryza sativa ObgC1 Acts as a Key Regulator of DNA Replication and Ribosome Biogenesis in Chloroplast Nucleoids

BACKGROUND

The Spo0B-associated GTP-binding protein (Obg) GTPase, has diverse and important functions in bacteria, including morphological development, DNA replication and ribosome maturation. Homologs of the Bacillus subtilis Obg have been also found in chloroplast of Oryza sativa, but their primary roles remain unknown.

RESULTS

We clarify that OsObgC1 is a functional homolog of AtObgC. The mutant obgc1-d1 exhibited hypersensitivity to the DNA replication inhibitor hydroxyurea. Quantitative PCR results showed that the ratio of chloroplast DNA to nuclear DNA in the mutants was higher than that of the wild-type plants. After DAPI staining, OsObgC1 mutants showed abnormal nucleoid architectures. The specific punctate staining pattern of OsObgC1-GFP signal suggests that this protein localizes to the chloroplast nucleoids. Furthermore, loss-of-function mutation in OsObgC1 led to a severe suppression of protein biosynthesis by affecting plastid rRNA processing. It was also demonstrated through rRNA profiling that plastid rRNA processing was decreased in obgc1-d mutants, which resulted in impaired ribosome biogenesis. The sucrose density gradient profiles revealed a defective chloroplast ribosome maturation of obgc1-d1 mutants.

CONCLUSION

Our findings here indicate that the OsObgC1 retains the evolutionarily biological conserved roles of prokaryotic Obg, which acts as a signaling hub that regulates DNA replication and ribosome biogenesis in chloroplast nucleoids.

Tissue and regional expression patterns of dicistronic tRNA-mRNA transcripts in grapevine (Vitis vinifera) and their evolutionary co-appearance with vasculature in land plants

Transfer RNAs (tRNA) are crucial adaptor molecules between messenger RNA (mRNA) and amino acids. Recent evidence in plants suggests that dicistronic tRNA-like structures also act as mobile signals for mRNA transcripts to move between distant tissues. Co-transcription is not a common feature in the plant nuclear genome and, in the few cases where polycistronic transcripts have been found, they include non-coding RNA species, such as small nucleolar RNAs and microRNAs. It is not known, however, the extent to which dicistronic transcripts of tRNA and mRNAs are expressed in field-grown plants, or the factors contributing to their expression. We analysed tRNA-mRNA dicistronic transcripts in the major horticultural crop grapevine (Vitis vinifera) using a novel pipeline developed to identify dicistronic transcripts from high-throughput RNA-sequencing data. We identified dicistronic tRNA-mRNA in leaf and berry samples from 22 commercial vineyards. Of the 124 tRNA genes that were expressed in both tissues, 18 tRNA were expressed forming part of 19 dicistronic tRNA-mRNAs. The presence and abundance of dicistronic molecules was tissue and geographic sub-region specific. In leaves, the expression patterns of dicistronic tRNA-mRNAs significantly correlated with tRNA expression, suggesting that their transcriptional regulation might be linked. We also found evidence of syntenic genomic arrangements of tRNAs and protein-coding genes between grapevine and Arabidopsis thaliana, and widespread prevalence of dicistronic tRNA-mRNA transcripts among vascular land plants but no evidence of these transcripts in non-vascular lineages. This suggests that the appearance of plant vasculature and tRNA-mRNA occurred concurrently during the evolution of land plants.

Plastomes in the holoparasitic family Balanophoraceae: Extremely high AT content, severe gene content reduction, and two independent genetic code changes

The transition to a heterotrophic lifestyle in angiosperms is characterized by convergent evolutionary changes. Plastid genome remodeling includes dramatic functional and physical reductions with the highest degrees observed in fully heterotrophic plants. Genes related to photosynthesis are generally absent or pseudogenized, while a few genes related to other metabolic processes that take place within the plastid are almost invariably maintained. The family Balanophoraceae consists of root holoparasites that present reduced plastid genomes with an extraordinarily elevated AT content and the single genetic code change ever documented in land plant plastomes (the stop codon TAG now codes for tryptophan). Here, we studied the plastomes of Lophophytum leandri and Ombrophytum subterraneum (Balanophoraceae) that showed the remarkable absence of the gene trnE, a highly biased nucleotide composition, and an independent genetic code change (the standard stop codon TGA codes for tryptophan). This is the second genetic code change identified in land plant plastomes. Analysis of the transcriptome of Lophophytum indicated that the entire C5 pathway typical of plants is conserved despite the lack of trnE in its plastome. A hypothetical model of plastome evolution in the Balanophoraceae is presented.

Two Novel PLS-Class Pentatricopeptide Repeat Proteins Are Involved in the Group II Intron Splicing of Mitochondrial Transcripts in the Moss Physcomitrella patens

Pentatricopeptide repeat (PPR) proteins are RNA-binding proteins that function in posttranscriptional regulation as gene-specific regulators of RNA metabolism in plant organelles. Plant PPR proteins are divided into four classes: P, PLS, E and DYW. The E- and DYW-class proteins are mainly implicated in RNA editing, whereas most of the P-class proteins predominantly participate in RNA cleavage, splicing and stabilization. In contrast, the functions of PLS-class proteins still remain obscure. Here, we report the function of PLS-class PpPPR_31 and PpPPR_9 in Physcomitrella patens. The knockout (KO) mutants of PpPPR_31 and PpPPR_9 exhibited slower protonema growth compared to the wild type. The PpPPR_31 KO mutants showed a considerable reduction in the splicing of nad5 intron 3 and atp9 intron 1. The PpPPR_9 KO mutants displayed severely reduced splicing of cox1 intron 3. An RNA electrophoresis mobility shift assay showed that the recombinant PpPPR_31 protein bound to the 5' region of nad5 exon 4 and the bulged A region in domain VI of atp9 group II intron 1 while the recombinant PpPPR_9 bound to the translated region of ORF622 in cox1 intron 3. These results suggest that a certain set of PLS-class PPR proteins may influence the splicing efficiency of mitochondrial group II introns.

Systematic sequencing of chloroplast transcript termini from Arabidopsis thaliana reveals >200 transcription initiation sites and the extensive imprints of RNA-binding proteins and secondary structures

Chloroplast transcription requires numerous quality control steps to generate the complex but selective mixture of accumulating RNAs. To gain insight into how this RNA diversity is achieved and regulated, we systematically mapped transcript ends by developing a protocol called Terminome-seq. Using Arabidopsis thaliana as a model, we catalogued >215 primary 5′ ends corresponding to transcription start sites (TSS), as well as 1628 processed 5′ ends and 1299 3′ ends. While most termini were found in intergenic regions, numerous abundant termini were also found within coding regions and introns, including several major TSS at unexpected locations. A consistent feature was the clustering of both 5′ and 3′ ends, contrasting with the prevailing description of discrete 5′ termini, suggesting an imprecision of the transcription and/or RNA processing machinery. Numerous termini correlated with the extremities of small RNA footprints or predicted stem-loop structures, in agreement with the model of passive RNA protection. Terminome-seq was also implemented for pnp1–1, a mutant lacking the processing enzyme polynucleotide phosphorylase. Nearly 2000 termini were altered in pnp1–1, revealing a dominant role in shaping the transcriptome. In summary, Terminome-seq permits precise delineation of the roles and regulation of the many factors involved in organellar transcriptome quality control.

Expanded inverted repeat region with large scale inversion in the first complete plastid genome sequence of Plantago ovata

Plantago ovata (Plantaginaceae) is an economically and medicinally important species, however, least is known about its genomics and evolution. Here, we report the first complete plastome genome of P. ovata and comparison with previously published genomes of related species from Plantaginaceae. The results revealed that P. ovata plastome size was 162,116 bp and that it had typical quadripartite structure containing a large single copy region of 82,084 bp and small single copy region of 5,272 bp. The genome has a markedly higher inverted repeat (IR) size of 37.4 kb, suggesting large-scale inversion of 13.8 kb within the expanded IR regions. In addition, the P. ovata plastome contains 149 different genes, including 43 tRNA, 8 rRNA, and 98 protein-coding genes. The analysis revealed 139 microsatellites, of which 71 were in the non-coding regions. Approximately 32 forward, 34 tandem, and 17 palindromic repeats were detected. The complete genome sequences, 72 shared genes, matK gene, and rbcL gene from related species generated the same phylogenetic signals, and phylogenetic analysis revealed that P. ovata formed a single clade with P. maritima and P. media. The divergence time estimation as employed in BEAST revealed that P. ovata diverged from P. maritima and P. media about 11.0 million years ago (Mya; 95% highest posterior density, 10.06-12.25 Mya). In conclusion, P. ovata had significant variation in the IR region, suggesting a more stable P. ovata plastome genome than that of other Plantaginaceae species.

PDM4, a Pentatricopeptide Repeat Protein, Affects Chloroplast Gene Expression and Chloroplast Development in Arabidopsis thaliana

Extensive studies have been carried out on chloroplast gene expression and chloroplast development; however, the regulatory mechanism is still largely unknown. Here, we characterized Pigment-Defective Mutant4 (PDM4), a P-type PPR protein localized in chloroplast. The pdm4 mutant showed seedling-lethal and albino phenotype under heterotrophic growth conditions. Transmission electron microscopic analysis revealed that thylakoid structure was totally disrupted in pdm4 mutant and eventually led to the breakdown of chloroplasts. The levels of several chloroplast- and nuclear-encoded proteins are strongly reduced in pdm4 mutant. Besides, transcript profile analysis detected that, in pdm4 mutant, the expression of plastid-encoded RNA polymerase-dependent genes was markedly affected, and deviant chloroplast rRNA processing was also observed. In addition, we found that PDM4 functions in the splicing of group II introns and may also be involved in the assembly of the 50S ribosomal particle. Our results demonstrate that PDM4 plays an important role in chloroplast gene expression and chloroplast development in Arabidopsis.

Complete chloroplast genomes of two Siraitia Merrill species: Comparative analysis, positive selection and novel molecular marker development

Siraitia grosvenorii fruit, known as Luo-Han-Guo, has been used as a traditional Chinese medicine for many years, and mogrosides are its primary active ingredients. Unfortunately, Siraitia siamensis, its wild relative, might be misused due to its indistinguishable appearance, not only threatening the reliability of the medication but also partly exacerbating wild resource scarcity. Therefore, high-resolution genetic markers must be developed to discriminate between these species. Here, the complete chloroplast genomes of S. grosvenorii and S. siamensis were assembled and analyzed for the first time; they were 158,757 and 159,190 bp in length, respectively, and possessed conserved quadripartite circular structures. Both contained 134 annotated genes, including 8 rRNA, 37 tRNA and 89 protein-coding genes. Twenty divergences (Pi > 0.03) were found in the intergenic regions. Nine protein-coding genes, accD, atpA, atpE, atpF, clpP, ndhF, psbH, rbcL, and rpoC2, underwent selection within Cucurbitaceae. Phylogenetic relationship analysis indicated that these two species originated from the same ancestor. Finally, four pairs of molecular markers were developed to distinguish the two species. The results of this study will be beneficial for taxonomic research, identification and conservation of Siraitia Merrill wild resources in the future.

Intraspecific Variation within the Utricularia amethystina Species Morphotypes Based on Chloroplast Genomes

Utricularia amethystina Salzm. ex A.St.-Hil. & Girard (Lentibulariaceae) is a highly polymorphic carnivorous plant taxonomically rearranged many times throughout history. Herein, the complete chloroplast genomes (cpDNA) of three U. amethystina morphotypes: purple-, white-, and yellow-flowered, were sequenced, compared, and putative markers for systematic, populations, and evolutionary studies were uncovered. In addition, RNA-Seq and RNA-editing analysis were employed for functional cpDNA evaluation. The cpDNA of three U. amethystina morphotypes exhibits typical quadripartite structure. Fine-grained sequence comparison revealed a high degree of intraspecific genetic variability in all morphotypes, including an exclusive inversion in the psbM and petN genes in U. amethystina yellow. Phylogenetic analyses indicate that U. amethystina morphotypes are monophyletic. Furthermore, in contrast to the terrestrial Utricularia reniformis cpDNA, the U. amethystina morphotypes retain all the plastid NAD(P)H-dehydrogenase (ndh) complex genes. This observation supports the hypothesis that the ndhs in terrestrial Utricularia were independently lost and regained, also suggesting that different habitats (aquatic and terrestrial) are not related to the absence of Utricularia ndhs gene repertoire as previously assumed. Moreover, RNA-Seq analyses recovered similar patterns, including nonsynonymous RNA-editing sites (e.g., rps14 and petB). Collectively, our results bring new insights into the chloroplast genome architecture and evolution of the photosynthesis machinery in the Lentibulariaceae.

Determination of the evolutionary pressure on Camellia oleifera on Hainan Island using the complete chloroplast genome sequence

Camellia oleifera is one of the four largest woody edible oil plants in the world with high ecological and medicinal values. Due to frequent interspecific hybridization, it was difficult to study its genetics and evolutionary history. This study used C. oleifera that was collected on Hainan Island to conduct our research. The unique island environment makes the quality of tea oil higher than that of other species grown in the mainland. Moreover, a long-term geographic isolation might affect gene structure. In order to better understand the molecular biology of this species, protect excellent germplasm resources, and promote the population genetics and phylogenetic studies of Camellia plants, high-throughput sequencing technology was used to obtain the chloroplast genome sequence of Hainan C. oleifera. The results showed that the whole chloroplast genome of C. oleifera in Hainan was 156,995 bp in length, with a typical quadripartite structure of a large single copy (LSC) region of 86,648 bp, a small single copy (SSC) region of 18,297 bp, and a pair of inverted repeats (IRs) of 26,025 bp. The whole genome encoded a total of 141 genes (115 different genes), including 88 protein-coding genes, 45 tRNA genes, and eight rRNA genes. Among these genes, nine genes contained one intron, two genes contained two introns, and four overlapping genes were also detected. The total GC content of Hainan C. oleifera's chloroplast genome was 37.29%. The chloroplast genome structure characteristics of Hainan C. oleifera were compared with mainland C. oleifera and those of the other eight closely related Theaceae species; it was found that the contractions and expansions of the IR/LSC and IR/SSC regions affected the length of chloroplast genome. The chloroplast genome sequences of these Theaceae species were highly similar. A comparative analysis indicated that the Theaceae species were conserved in structure and evolution. A total of 51 simple sequence repeat (SSR) loci were detected in the chloroplast genome of Hainan C. oleifera, and all Camellia plants did not have pentanucleotide repeats, which could be used as a good marker in phylogenetic studies. We also detected seven long repeats, the base composition of all repeats was biased toward A/T, which was consistent with the codon bias. It was found that Hainan C. oleifera had a similar evolutionary relationship with C. crapnelliana, through the use of codons and phylogenetic analysis. This study can provide an effective genomic resource for the evolutionary history of Theaceae family.

Prediction and large-scale analysis of primary operons in plastids reveals unique genetic features in the evolution of chloroplasts

While bacterial operons have been thoroughly studied, few analyses of chloroplast operons exist, limiting the ability to study fundamental elements of these structures and utilize them for synthetic biology. Here, we describe the creation of a plastome-specific operon database (link provided below) achieved by combining experimental tools and predictive modeling. Using a Reverse-Transcription-PCR based method and published data, we determined the transcription-state of 213 gene pairs from four plastomes of evolutionary distinct organisms. By analyzing sequence-based features computed for our dataset, we were able to highlight fundamental characteristics differentiating between operon pairs and non-operon pairs. These include an interesting tendency toward maintaining similar messenger RNA-folding profiles in operon gene pairs, a feature that failed to yield any informative separation in cyanobacteria, suggesting that it catches unique traits of operon gene expression, which have evolved post-endosymbiosis. Subsequently, we used this feature set to train a random-forest classifier for operon prediction. As our results demonstrate the ability of our predictor to obtain accurate (84%) and robust predictions on unlabeled datasets, we proceeded to building operon maps for 2018 sequenced plastids. Our database may now present new opportunities for promoting metabolic engineering and synthetic biology in chloroplasts.

Divergence of RNA editing among Arabidopsis species

RNA editing altered the RNA sequence by replacing the C nucleotide to U in the organellar genomes of plants. RNA editing status sometimes differed among distant species. The pattern of conservation and variation of RNA editing status made it possible to evaluate evolutionary mechanisms impacting functional aspects of RNA editing. In this study, divergence of RNA editing in the chloroplast genome among Arabidopsis species was analyzed to determine 9 losses and 1 gain in RNA editing. All changes in A. thaliana lineage resulted from changes to the chloroplast genome sequence, whereas changes in the A. lyrata / halleri lineage were possibly due to exclusive changes in the nuclear editing factors. One loss of RNA editing in A. lyrata was caused by a deficiency in the PPR gene OTP80. The changes in RNA editing occurred approximately every two million years and were not observed at functionally important sites. These results highlight the conserved nature of RNA editing status suggesting the importance of RNA editing during evolution.

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 chloroplast genomes of Paris Sect. Marmorata: insights into repeat regions and evolutionary implications

Background

Species of Paris Sect. Marmorata are valuable medicinal plants to synthesize steroidal saponins with effective pharmacological therapy. However, the wild resources of the species are threatened by plundering exploitation before the molecular genetics studies uncover the genomes and evolutionary significance. Thus, the availability of complete chloroplast genome sequences of Sect. Marmorata is necessary and crucial to the understanding the plastome evolution of this section and facilitating future population genetics studies. Here, we determined chloroplast genomes of Sect. Marmorata, and conducted the whole chloroplast genome comparison.

Results

This study presented detailed sequences and structural variations of chloroplast genomes of Sect. Marmorata. Over 40 large repeats and approximately 130 simple sequence repeats as well as a group of genomic hotspots were detected. Inverted repeat contraction of this section was inferred via comparing the chloroplast genomes with the one of P. verticillata. Additionally, almost all the plastid protein coding genes were found to prefer ending with A/U. Mutation bias and selection pressure predominately shaped the codon bias of most genes. And most of the genes underwent purifying selection, whereas photosynthetic genes experienced a relatively relaxed purifying selection.

Conclusions

Repeat sequences and hotspot regions can be scanned to detect the intraspecific and interspecific variability, and selected to infer the phylogenetic relationships of Sect. Marmorata and other species in subgenus Daiswa. Mutation and natural selection were the main forces to drive the codon bias pattern of most plastid protein coding genes. Therefore, this study enhances the understanding about evolution of Sect. Marmorata from the chloroplast genome, and provide genomic insights into genetic analyses of Sect. Marmorata.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5281-x) contains supplementary material, which is available to authorized users.

Molecular and Functional Diversity of RNA Editing in Plant Mitochondria

RNA editing is a fundamental biochemical process relating to the modification of nucleotides in messenger RNAs of functional genes in cells. RNA editing leads to re-establishment of conserved amino acid residues for functional proteins in nuclei, chloroplasts, and mitochondria. Identification of RNA editing factors that contributes to target site recognition increases our understanding of RNA editing mechanisms. Significant progress has been made in recent years in RNA editing studies for both animal and plant cells. RNA editing in nuclei and mitochondria of animal cells and in chloroplast of plant cells has been extensively documented and reviewed. RNA editing has been also extensively documented on plant mitochondria. However, functional diversity of RNA editing factors in plant mitochondria is not overviewed. Here, we review the biological significance of RNA editing, recent progress on the molecular mechanisms of RNA editing process, and function diversity of editing factors in plant mitochondrial research. We will focus on: (1) pentatricopeptide repeat proteins in Arabidopsis and in crop plants; (2) the progress of RNA editing process in plant mitochondria; (3) RNA editing-related RNA splicing; (4) RNA editing associated flower development; (5) RNA editing modulated male sterile; (6) RNA editing-regulated cell signaling; and (7) RNA editing involving abiotic stress. Advances described in this review will be valuable in expanding our understanding in RNA editing. The diverse functions of RNA editing in plant mitochondria will shed light on the investigation of molecular mechanisms that underlies plant development and abiotic stress tolerance.

Complete Chloroplast Genome Sequence of Broomcorn Millet (Panicum miliaceum L.) and Comparative Analysis with Other Panicoideae Species

Broomcorn millet (Panicum miliaceum L.) is one of the earliest domesticated cereals worldwide, holding significant agricultural, historical, and evolutionary importance. However, our genomic knowledge of it is rather limited at present, hampering further genetic and evolutionary studies. Here, we sequenced and assembled the chloroplast genome (cp) of broomcorn millet and compared it with five other Panicoideae species. Results showed that the cp genome of broomcorn millet was 139,826 bp in size, with a typical quadripartite structure. In total, 108 genes were annotated and 18 genes were duplicated in the IR (inverted region) region, which was similar to other Panicoideae species. Comparative analysis showed a rather conserved genome structure between them, with three common regions. Furthermore, RNA editing, codon usage, and expansion of the IR, as well as simple sequence repeat (SSR) elements, were systematically investigated and 13 potential DNA markers were developed for Panicoideae species identification. Finally, phylogenetic analysis implied that broomcorn millet was a sister species to Panicum virgatum within the tribe Paniceae, and supported a monophyly of the Panicoideae. This study has reported for the first time the genome organization, gene content, and structural features of the chloroplast genome of broomcorn millet, which provides valuable information for genetic and evolutionary studies in the genus Panicum and beyond.

A nuclear-encoded protein, mTERF6, mediates transcription termination of rpoA polycistron for plastid-encoded RNA polymerase-dependent chloroplast gene expression and chloroplast development

The expression of plastid genes is regulated by two types of DNA-dependent RNA polymerases, plastid-encoded RNA polymerase (PEP) and nuclear-encoded RNA polymerase (NEP). The plastid rpoA polycistron encodes a series of essential chloroplast ribosome subunits and a core subunit of PEP. Despite the functional importance, little is known about the regulation of rpoA polycistron. In this work, we show that mTERF6 directly associates with a 3′-end sequence of rpoA polycistron in vitro and in vivo, and that absence of mTERF6 promotes read-through transcription at this site, indicating that mTERF6 acts as a factor required for termination of plastid genes’ transcription in vivo. In addition, the transcriptions of some essential ribosome subunits encoded by rpoA polycistron and PEP-dependent plastid genes are reduced in the mterf6 knockout mutant. RpoA, a PEP core subunit, accumulates to about 50% that of the wild type in the mutant, where early chloroplast development is impaired. Overall, our functional analyses of mTERF6 provide evidence that it is more likely a factor required for transcription termination of rpoA polycistron, which is essential for chloroplast gene expression and chloroplast development.

First chloroplast genomics study of Phoenix dactylifera (var. Naghal and Khanezi): A comparative analysis

Date palm (Phoenix dactylifera L.) is one of the oldest fruit crops in the arid regions of the Middle East. However, little information is available regarding its plastid genomes. In this study, we sequenced the chloroplast (cp) genomes of two economically important but genomically unexplored date palm cultivars of Phoenix dactylifera (var. Naghal and Khanezi). The data assembly and genome annotation revealed a typical quadripartite structure similar to Arecaceae, and the genome sizes of Naghal and Khanezi were 158,210 bp and 158,211 bp, respectively. Structurally, both cp genomes were comprised of four regions: a pair of inverted repeats (27,273 bp for Khanezi and for Naghal 27,272 bp), a large single-copy region (86,090 bp and 86,092 bp) and a small single-copy region (17,575 bp and 17,574 bp). Both genomes had 138 representative genes, whereas 227 and 229 randomly distributed microsatellites were also observed in Khanezi and Naghal, respectively. Phylogenetic analysis based on the whole cp genomes and 68 shared genes showed identical phylogenetic trees of Khanezi and Naghal forming clades with Khalas and Aseel cultivars, respectively. The current study showed detailed comparative cp genome analysis, which could be essential for broader population genetics and molecular studies of these four date palm cultivars.

Generation and characterization of a collection of knock-down lines for the chloroplast Clp protease complex in tobacco

Protein degradation in chloroplasts is carried out by a set of proteases that eliminate misfolded, damaged, or superfluous proteins. The ATP-dependent caseinolytic protease (Clp) is the most complex protease in plastids and has been implicated mainly in stromal protein degradation. In contrast, FtsH, a thylakoid membrane-associated metalloprotease, is believed to participate mainly in the degradation of thylakoidal proteins. To determine the role of specific Clp and FtsH subunits in plant growth and development, RNAi lines targeting at least one subunit of each Clp ring and FtsH were generated in tobacco. In addition, mutation of the translation initiation codon was employed to down-regulate expression of the plastid-encoded ClpP1 subunit. These protease lines cover a broad range of reductions at the transcript and protein levels of the targeted genes. A wide spectrum of phenotypes was obtained, including pigment deficiency, alterations in leaf development, leaf variegations, and impaired photosynthesis. When knock-down lines for the different protease subunits were compared, both common and specific phenotypes were observed, suggesting distinct functions of at least some subunits. Our work provides a well-characterized collection of knock-down lines for plastid proteases in tobacco and reveals the importance of the Clp protease in physiology and plant development.

Epstein-Barr virus gene expression and latency pattern in gastric carcinomas: a systematic review

METHODS

A systematic review of literature was conducted to identify all published reports regarding the expression of Epstein-Barr Virus (EBV) proteins/transcripts and EBV latency patterns in EBV-associated gastric carcinomas (EBVaGC).

RESULTS

The literature search retrieved 247 papers, of which 25 papers matched the inclusion criteria. The analysis reveals that the most frequently expressed EBV latent proteins are EBNA1 (98.1%) and LMP2A (53.8%), while LMP1 and LMP2B are present in only 10% of cases. Lytic proteins, such as BARF0 and BARF1, and other lytic transcripts are present in almost half of cases.

CONCLUSION

EBVaGC seems to display a unique transcription/latency pattern that does not fit the 'standard' EBV latency patterns and therefore should be further studied to better understand EBVaGC carcinogenesis.

Plastome Rearrangements in the "Adenocalymma-Neojobertia" Clade (Bignonieae, Bignoniaceae) and Its Phylogenetic Implications

The chloroplast is one of the most important organelles of plants. This organelle has a circular DNA with approximately 130 genes. The use of plastid genomic data in phylogenetic and evolutionary studies became possible with high-throughput sequencing methods, which allowed us to rapidly obtain complete genomes at a reasonable cost. Here, we use high-throughput sequencing to study the "Adenocalymma-Neojobertia" clade (Bignonieae, Bignoniaceae). More specifically, we use Hi-Seq Illumina technology to sequence 10 complete plastid genomes. Plastomes were assembled using selected plastid reads and de novo approach with SPAdes. The 10 assembled genomes were analyzed in a phylogenetic context using five different partition schemes: (1) 91 protein-coding genes ("coding"); (2) 76 introns and spacers with alignment manually edited ("non-coding edited"); (3) 76 non-coding regions with poorly aligned regions removed using T-Coffee ("non-coding filtered"); (4) 91 coding regions plus 76 non-coding regions edited ("coding + non-coding edited"); and, (5) 91 protein-coding regions plus the 76 filtered non-coding regions ("coding + non-coding filtered"). Fragmented regions were aligned using Mafft. Phylogenetic analyses were conducted using Maximum Likelihood (ML) and Bayesian Criteria (BC). The analyses of the individual plastomes consistently recovered an expansion of the Inverted Repeated (IRs) regions and a compression of the Small Single Copy (SSC) region. Major genomic translocations were observed at the Large Single Copy (LSC) and IRs. ML phylogenetic analyses of the individual datasets led to the same topology, with the exception of the analysis of the "non-coding filtered" dataset. Overall, relationships were strongly supported, with the highest support values obtained through the analysis of the "coding + non-coding edited" dataset. Four regions at the LSC, SSC, and IR were selected for primer development. The "Adenocalymma-Neojobertia" clade shows an unusual pattern of plastid structure variation, including four major genomic translocations. These rearrangements challenge the current view of conserved plastid genome architecture in terms of gene order. It also complicates both genomic assemblies using reference genomes and sequence alignments using whole plastomes. Therefore, strategies that employ de novo assemblies and manual evaluation of sequence alignments are required to prevent assembly and alignment errors.

DELAYED GREENING 238, a Nuclear-Encoded Chloroplast Nucleoid Protein, Is Involved in the Regulation of Early Chloroplast Development and Plastid Gene Expression in Arabidopsis thaliana

Chloroplast development is an essential process for plant growth that is regulated by numerous proteins. Plastid-encoded plastid RNA polymerase (PEP) is a large complex that regulates plastid gene transcription and chloroplast development. However, many proteins in this complex remain to be identified. Here, through large-scale screening of Arabidopsis mutants by Chl fluorescence imaging, we identified a novel protein, DELAYED GREENING 238 (DG238), which is involved in regulating chloroplast development and plastid gene expression. Loss of DG238 retards plant growth, delays young leaf greening, affects chloroplast development and lowers photosynthetic efficiency. Moreover, blue-native PAGE (BN-PAGE) and Western blot analysis indicated that PSII and PSI protein levels are reduced in dg238 mutants. DG238 is mainly expressed in young tissues and is regulated by light signals. Subcellular localization analysis showed that DG238 is a nuclear-encoded chloroplast nucleoid protein. More interestingly, DG238 was co-expressed with FLN1, which encodes an essential subunit of the PEP complex. Bimolecular fluorescence complementation (BiFC) and co-immunoprecipitation (Co-IP) assays showed that DG238 can also interact with FLN1. Taken together, these results suggest that DG238 may function as a component of the PEP complex that is important for the early stage of chloroplast development and helps regulate PEP-dependent plastid gene expression.

Maize Cytolines Unmask Key Nuclear Genes That Are under the Control of Retrograde Signaling Pathways in Plants

The genomes of the two plant organelles encode for a relatively small number of proteins. Thus, nuclear genes encode the vast majority of their proteome. Organelle-to-nucleus communication takes place through retrograde signaling (RS) pathways. Signals relayed through RS pathways have an impact on nuclear gene expression but their target-genes remain elusive in a normal state of the cell (considering that only mutants and stress have been used so far). Here, we use maize cytolines as an alternative. The nucleus of a donor line was transferred into two other cytoplasmic environments through at least nine back-crosses, in a time-span of > 10 years. The transcriptomes of the resulting cytolines were sequenced and compared. There are 96 differentially regulated nuclear genes in two cytoplasm-donor lines when compared with their nucleus-donor. They are expressed throughout plant development, in various tissues and organs. One-third of the 96 proteins have a human homolog, stressing their potential role in mitochondrial RS. We also identified syntenic orthologous genes in four other grasses and homologous genes in Arabidopsis thaliana. These findings contribute to the paradigm we use to describe the RS in plants. The 96 nuclear genes identified here are not differentially regulated as a result of mutation, or any kind of stress. They are rather key players of the organelle-to-nucleus communication in a normal state of the cell.

Full transcription of the chloroplast genome in photosynthetic eukaryotes

Prokaryotes possess a simple genome transcription system that is different from that of eukaryotes. In chloroplasts (plastids), it is believed that the prokaryotic gene transcription features govern genome transcription. However, the polycistronic operon transcription model cannot account for all the chloroplast genome (plastome) transcription products at whole-genome level, especially regarding various RNA isoforms. By systematically analyzing transcriptomes of plastids of algae and higher plants, and cyanobacteria, we find that the entire plastome is transcribed in photosynthetic green plants, and that this pattern originated from prokaryotic cyanobacteria - ancestor of the chloroplast genomes that diverged about 1 billion years ago. We propose a multiple arrangement transcription model that multiple transcription initiations and terminations combine haphazardly to accomplish the genome transcription followed by subsequent RNA processing events, which explains the full chloroplast genome transcription phenomenon and numerous functional and/or aberrant pre-RNAs. Our findings indicate a complex prokaryotic genome regulation when processing primary transcripts.

Dysfunctional chloroplasts up-regulate the expression of mitochondrial genes in Arabidopsis seedlings

Chloroplasts and mitochondria play important roles in maintaining metabolic and energy homeostasis in the plant cell. The interactions between these two organelles, especially photosynthesis and respiration, have been intensively studied. Still, little is known about the regulation of mitochondrial gene expression by chloroplasts and vice versa. The gene expression machineries in chloroplasts and mitochondria rely heavily on the nuclear genome. Thus, the interactions between nucleus and these organelles, including anterograde and retrograde regulation, have been actively investigated in the last two decades. Norflurazon (NF) and lincomycin (Lin) are two commonly used inhibitors to study chloroplast-to-nucleus retrograde signaling in plants. We used NF and Lin to block the development and functions of chloroplasts and examined their effects on mitochondrial gene expression, RNA editing and splicing. The editing of most mitochondrial transcripts was not affected, but the editing extents of nad4-107, nad6-103, and ccmFc-1172 decreased slightly in NF- and Lin-treated seedlings. While the splicing of mitochondrial transcripts was not significantly affected, steady-state mRNA levels of several mitochondrial genes increased significantly in NF- and Lin-treated seedlings. Moreover, Lin seemed to have more profound effects than NF on the expression of mitochondrial genes, indicating that signals derived from these two inhibitors might be distinct. NF and Lin also significantly induced the expression of nuclear genes encoding subunits of mitochondrial electron transport chain complexes. Thus, dysfunctional chloroplasts may coordinately up-regulate the expression of nuclear and mitochondrial genes encoding subunits of respiratory complexes.

Expression of chloroperoxidase from Pseudomonas pyrrocinia in tobacco plastids for fungal resistance

The chloroperoxidase (cpo) gene from Pseudomonas pyrrocinia was transformed into the plastid genome (plastome) of Nicotiana tabacum var. Petit Havana and transplastomic lines were compared with a nuclear transformant for the same gene. Southern analysis confirmed integration in the plastome and western blotting confirmed the presence of the chloroperoxidase protein (CPO) in higher abundance in transplastomic plants than in cpo nuclear transformants. Northern analysis of primary plastome transformants for cpo showed 15-fold higher transcript abundance than in the nuclear transformant, yet this extent of enhancement was not observed in western blot, enzyme or bioassay, indicating a bottleneck at the post-transcriptional level. Representative plants from the two transplastomic lines showed resistance to fungal pathogens in vitro (Aspergillus flavus, Fusarium verticillioides, and Verticillium dahliae) and in planta (Alternaria alternata).

Analyses of charophyte chloroplast genomes help characterize the ancestral chloroplast genome of land plants

Despite the significance of the relationships between embryophytes and their charophyte algal ancestors in deciphering the origin and evolutionary success of land plants, few chloroplast genomes of the charophyte algae have been reconstructed to date. Here, we present new data for three chloroplast genomes of the freshwater charophytes Klebsormidium flaccidum (Klebsormidiophyceae), Mesotaenium endlicherianum (Zygnematophyceae), and Roya anglica (Zygnematophyceae). The chloroplast genome of Klebsormidium has a quadripartite organization with exceptionally large inverted repeat (IR) regions and, uniquely among streptophytes, has lost the rrn5 and rrn4.5 genes from the ribosomal RNA (rRNA) gene cluster operon. The chloroplast genome of Roya differs from other zygnematophycean chloroplasts, including the newly sequenced Mesotaenium, by having a quadripartite structure that is typical of other streptophytes. On the basis of the improbability of the novel gain of IR regions, we infer that the quadripartite structure has likely been lost independently in at least three zygnematophycean lineages, although the absence of the usual rRNA operonic synteny in the IR regions of Roya may indicate their de novo origin. Significantly, all zygnematophycean chloroplast genomes have undergone substantial genomic rearrangement, which may be the result of ancient retroelement activity evidenced by the presence of integrase-like and reverse transcriptase-like elements in the Roya chloroplast genome. Our results corroborate the close phylogenetic relationship between Zygnematophyceae and land plants and identify 89 protein-coding genes and 22 introns present in the chloroplast genome at the time of the evolutionary transition of plants to land, all of which can be found in the chloroplast genomes of extant charophytes.

The circadian regulation of photosynthesis

Correct circadian regulation increases plant productivity, and photosynthesis is circadian-regulated. Here, we discuss the regulatory basis for the circadian control of photosynthesis. We discuss candidate mechanisms underpinning circadian oscillations of light harvesting and consider how the circadian clock modulates CO2 fixation by Rubisco. We show that new techniques may provide a platform to better understand the signalling pathways that couple the circadian clock with the photosynthetic apparatus. Finally, we discuss how understanding circadian regulation in model systems is underpinning research into the impact of circadian regulation in crop species.

Plastid mRNA translation

Overall translational machinery in plastids is similar to that of E. coli. Initiation is the crucial step for translation and this step in plastids is somewhat different from that of E. coli. Unlike the Shine-Dalgarno sequence in E. coli, cis-elements for translation initiation are not well conserved in plastid mRNAs. Specific trans-acting factors are generally required for translation initiation and its regulation in plastids. During translation elongation, ribosomes pause sometimes on photosynthesis-related mRNAs due probably to proper insertion of nascent polypeptides into membrane complexes. Codon usage of plastid mRNAs is different from that of E. coli and mammalian cells. Plastid mRNAs do not have the so-called rare codons. Translation efficiencies of several synonymous codons are not always correlated with codon usage in plastid mRNAs.

Complete sequencing of five araliaceae chloroplast genomes and the phylogenetic implications

Background

The ginseng family (Araliaceae) includes a number of economically important plant species. Previously phylogenetic studies circumscribed three major clades within the core ginseng plant family, yet the internal relationships of each major group have been poorly resolved perhaps due to rapid radiation of these lineages. Recent studies have shown that phyogenomics based on chloroplast genomes provides a viable way to resolve complex relationships.

Methodology/Principal Findings

We report the complete nucleotide sequences of five Araliaceae chloroplast genomes using next-generation sequencing technology. The five chloroplast genomes are 156,333–156,459 bp in length including a pair of inverted repeats (25,551–26,108 bp) separated by the large single-copy (86,028–86,566 bp) and small single-copy (18,021–19,117 bp) regions. Each chloroplast genome contains the same 114 unique genes consisting of 30 transfer RNA genes, four ribosomal RNA genes, and 80 protein coding genes. Gene size, content, and order, AT content, and IR/SC boundary structure are similar among all Araliaceae chloroplast genomes. A total of 140 repeats were identified in the five chloroplast genomes with palindromic repeat as the most common type. Phylogenomic analyses using parsimony, likelihood, and Bayesian inference based on the complete chloroplast genomes strongly supported the monophyly of the Asian Palmate group and the Aralia-Panax group. Furthermore, the relationships among the sampled taxa within the Asian Palmate group were well resolved. Twenty-six DNA markers with the percentage of variable sites higher than 5% were identified, which may be useful for phylogenetic studies of Araliaceae.

Conclusion

The chloroplast genomes of Araliaceae are highly conserved in all aspects of genome features. The large-scale phylogenomic data based on the complete chloroplast DNA sequences is shown to be effective for the phylogenetic reconstruction of Araliaceae.

Strategies for metabolic pathway engineering with multiple transgenes

The engineering of metabolic pathways in plants often requires the concerted expression of more than one gene. While with traditional transgenic approaches, the expression of multiple transgenes has been challenging, recent progress has greatly expanded our repertoire of powerful techniques making this possible. New technological options include large-scale co-transformation of the nuclear genome, also referred to as combinatorial transformation, and transformation of the chloroplast genome with synthetic operon constructs. This review describes the state of the art in multigene genetic engineering of plants. It focuses on the methods currently available for the introduction of multiple transgenes into plants and the molecular mechanisms underlying successful transgene expression. Selected examples of metabolic pathway engineering are used to illustrate the attractions and limitations of each method and to highlight key factors that influence the experimenter's choice of the best strategy for multigene engineering.

Comparative chloroplast genomes of camellia species

BACKGROUND

Camellia, comprising more than 200 species, is a valuable economic commodity due to its enormously popular commercial products: tea leaves, flowers, and high-quality edible oils. It is the largest and most important genus in the family Theaceae. However, phylogenetic resolution of the species has proven to be difficult. Consequently, the interspecies relationships of the genus Camellia are still hotly debated. Phylogenomics is an attractive avenue that can be used to reconstruct the tree of life, especially at low taxonomic levels.

METHODOLOGY/PRINCIPAL FINDINGS

Seven complete chloroplast (cp) genomes were sequenced from six species representing different subdivisions of the genus Camellia using Illumina sequencing technology. Four junctions between the single-copy segments and the inverted repeats were confirmed and genome assemblies were validated by PCR-based product sequencing using 123 pairs of primers covering preliminary cp genome assemblies. The length of the Camellia cp genome was found to be about 157kb, which contained 123 unique genes and 23 were duplicated in the IR regions. We determined that the complete Camellia cp genome was relatively well conserved, but contained enough genetic differences to provide useful phylogenetic information. Phylogenetic relationships were analyzed using seven complete cp genomes of six Camellia species. We also identified rapidly evolving regions of the cp genome that have the potential to be used for further species identification and phylogenetic resolution.

CONCLUSIONS/SIGNIFICANCE

In this study, we wanted to determine if analyzing completely sequenced cp genomes could help settle these controversies of interspecies relationships in Camellia. The results demonstrate that cp genome data are beneficial in resolving species definition because they indicate that organelle-based "barcodes", can be established for a species and then used to unmask interspecies phylogenetic relationships. It reveals that phylogenomics based on cp genomes is an effective approach for achieving phylogenetic resolution between Camellia species.

Differential regulation of Arabidopsis plastid gene expression and RNA editing in non-photosynthetic tissues

RNA editing is one of the post-transcriptional processes that commonly occur in plant plastids and mitochondria. In Arabidopsis, 34 C-to-U RNA editing events, affecting transcripts of 18 plastid genes, have been identified. Here, we examined the editing and expression of these transcripts in different organs, and in green and non-green seedlings (etiolated, cia5-2, ispF and ispG albino mutants, lincomycin-, and norflurazon-treated). The editing efficiency of Arabidopsis plastid transcripts varies from site to site, and may be specifically regulated in different tissues. Steady state levels of plastid transcripts are low or undetectable in etiolated seedlings, but most editing sites are edited with efficiencies similar to those observed in green seedlings. By contrast, the editing of some sites is completely lost or significantly reduced in other non-green tissues; for instance, the editing of ndhB-149, ndhB-1255, and ndhD-2 is completely lost in roots and in lincomycin-treated seedlings. The editing of ndhD-2 is also completely lost in albino mutants and norflurazon-treated seedlings. However, matK-640 is completely edited, and accD-794, atpF-92, psbE-214, psbF-77, psbZ-50, and rps14-50 are completely or highly edited in both green and non-green tissues. In addition, the expression of nucleus-encoded RNA polymerase dependent transcripts is specifically induced by lincomycin, and the splicing of ndhB transcripts is significantly reduced in the albino mutants and inhibitor-treated seedlings. Our results indicate that plastid gene expression, and the splicing and editing of plastid transcripts are specifically and differentially regulated in various types of non-green tissues.

Complete chloroplast genome of the genus Cymbidium: lights into the species identification, phylogenetic implications and population genetic analyses

BACKGROUND

Cymbidium orchids, including some 50 species, are the famous flowers, and they possess high commercial value in the floricultural industry. Furthermore, the values of different orchids are great differences. However, species identification is very difficult. To a certain degree, chloroplast DNA sequence data are a versatile tool for species identification and phylogenetic implications in plants. Different chloroplast loci have been utilized for evaluating phylogenetic relationships at each classification level among plant species, including at the interspecies and intraspecies levels. However, there is no evidence that a short sequence can distinguish all plant species from each other in order to infer phylogenetic relationships. Molecular markers derived from the complete chloroplast genome can provide effective tools for species identification and phylogenetic resolution.

RESULTS

The complete nucleotide sequences of eight individuals from a total of five Cymbidium species' chloroplast (cp) genomes were determined using Illumina sequencing technology of the total DNA via a combination of de novo and reference-guided assembly. The length of the Cymbidium cp genome is about 155 kb. The cp genomes contain 123 unique genes, and the IR regions contain 24 duplicates. Although the genomes, including genome structure, gene order and orientation, are similar to those of other orchids, they are not evolutionarily conservative. The cp genome of Cymbidium evolved moderately with more than 3% sequence divergence, which could provide enough information for phylogeny. Rapidly evolving chloroplast genome regions were identified and 11 new divergence hotspot regions were disclosed for further phylogenetic study and species identification in Orchidaceae.

CONCLUSIONS

Phylogenomic analyses were conducted using 10 complete chloroplast genomes from seven orchid species. These data accurately identified the individuals and established the phylogenetic relationships between the species. The results reveal that phylogenomics based on organelle genome sequencing lights the species identification-organelle-scale "barcodes", and is also an effective approach for studying whole populations and phylogenetic characteristics of Cymbidium.

Contradiction between plastid gene transcription and function due to complex posttranscriptional splicing: an exemplary study of ycf15 function and evolution in angiosperms

Plant chloroplast genes are usually co-transcribed while its posttranscriptional splicing is fairly complex and remains largely unsolved. On basis of sequencing the three complete Camellia (Theaceae) chloroplast genomes for the first time, we comprehensively analyzed the evolutionary patterns of ycf15, a plastid gene quite paradoxical in terms of its function and evolution, along the inferred angiosperm phylogeny. Although many species in separate lineages including the three species reported here contained an intact ycf15 gene in their chloroplast genomes, the phylogenetic mixture of both intact and obviously disabled ycf15 genes imply that they are all non-functional. Both intracellular gene transfer (IGT) and horizontal gene transfer (HGT) failed to explain such distributional anomalies. While, transcriptome analyses revealed that ycf15 was transcribed as precursor polycistronic transcript which contained ycf2, ycf15 and antisense trnL-CAA. The transcriptome assembly was surprisingly found to cover near the complete Camellia chloroplast genome. Many non-coding regions including pseudogenes were mapped by multiple transcripts, indicating the generality of pseudogene transcriptions. Our results suggest that plastid DNA posttranscriptional splicing may involve complex cleavage of non-functional genes.

SIG1, a sigma factor for the chloroplast RNA polymerase, differently associates with multiple DNA regions in the chloroplast chromosomes in vivo

Chloroplasts have their own DNA and gene expression systems. Transcription in chloroplasts is regulated by two types of RNA polymerase, nuclear-encoded plastid RNA polymerase (NEP) and plastid-encoded plastid RNA polymerase (PEP), and multiple sigma factors for PEP. To study transcriptional regulation in chloroplasts, a molecular genetic approach has extensively been used. However, this method may include indirect effects, and it cannot be applied to the analysis of factors essential to survival. These limitations make understanding specific regulation by transcription factors difficult. Chromatin immunoprecipitation (ChIP) is a powerful and useful tool for obtaining information on transcription-factor binding sites; it can directly detect dynamic changes in their interaction patterns in vivo. To further understand transcriptional regulation in chloroplasts, we here established a ChIP-based method in Arabidopsis thaliana and analyzed the binding pattern of a chloroplast sigma factor, SIG1. We found that SIG1 specifically binds to newly identified target promoters as well as to a set of promoters of genes whose mRNA expression is dependent on OsSIG1 in rice and that this binding changed in response to high-light stress. These results suggested that the ChIP-based approach is very useful in understanding transcriptional regulation of chloroplast genes and can overcome several problems posed by conventional methods.

Functional characterization of ObgC in ribosome biogenesis during chloroplast development

The Spo0B-associated GTP-binding protein (Obg) GTPase, essential for bacterial viability, is also conserved in eukaryotes, but its primary role in eukaryotes remains unknown. Here, our functional characterization of Arabidopsis and rice obgc mutants strongly underlines the evolutionarily conserved role of eukaryotic Obgs in organellar ribosome biogenesis. The mutants exhibited a chlorotic phenotype, caused by retarded chloroplast development. A plastid DNA macroarray revealed a plastid-encoded RNA polymerase (PEP) deficiency in an obgc mutant, caused by incompleteness of the PEP complex, as its western blot exhibited reduced levels of RpoA protein, a component of PEP. Plastid rRNA profiling indicated that plastid rRNA processing is defective in obgc mutants, probably resulting in impaired ribosome biogenesis and, in turn, in reduced levels of RpoA protein. RNA co-immunoprecipitation revealed that ObgC specifically co-precipitates with 23S rRNA in vivo. These findings indicate that ObgC functions primarily in plastid ribosome biogenesis during chloroplast development. Furthermore, complementation analysis can provide new insights into the functional modes of three ObgC domains, including the Obg fold, G domain and OCT.

External light conditions and internal cell cycle phases coordinate accumulation of chloroplast and mitochondrial transcripts in the red alga Cyanidioschyzon merolae

The mitochondria and chloroplasts in plant cells are originated from bacterial endosymbioses, and they still replicate their own genome and divide in a similar manner as their ancestors did. It is thus likely that the organelle transcription is coordinated with its proliferation cycle. However, this possibility has not extensively been explored to date, because in most plant cells there are many mitochondria and chloroplasts that proliferate asynchronously. It is generally believed that the gene transfer from the organellar to nuclear genome has enabled nuclear control of the organelle functions during the evolution of eukaryotic plant cells. Nevertheless, no significant relationship has been reported between the organelle transcriptome and the host cell cycle even in Chlamydomonas reinhardtii. While the organelle proliferation cycle is not coordinated with the cell cycle in vascular plants, in the unicellular red alga Cyanidioschyzon merolae that contains only one mitochondrion, one chloroplast, and one nucleus per cell, each of the organelles is known to proliferate at a specific phase of the cell cycle. Here, we show that the expression of most of the organelle genes is highly coordinated with the cell cycle phases as well as with light regimes in clustering analyses. In addition, a strong correlation was observed between the gene expression profiles in the mitochondrion and chloroplast, resulting in the identification of a network of functionally related genes that are co-expressed during organelle proliferation.

Complete chloroplast genome sequence of a major invasive species, crofton weed (Ageratina adenophora)

Background

Crofton weed (Ageratina adenophora) is one of the most hazardous invasive plant species, which causes serious economic losses and environmental damages worldwide. However, the sequence resource and genome information of A. adenophora are rather limited, making phylogenetic identification and evolutionary studies very difficult. Here, we report the complete sequence of the A. adenophora chloroplast (cp) genome based on Illumina sequencing.

Methodology/Principal Findings

The A. adenophora cp genome is 150, 689 bp in length including a small single-copy (SSC) region of 18, 358 bp and a large single-copy (LSC) region of 84, 815 bp separated by a pair of inverted repeats (IRs) of 23, 755 bp. The genome contains 130 unique genes and 18 duplicated in the IR regions, with the gene content and organization similar to other Asteraceae cp genomes. Comparative analysis identified five DNA regions (ndhD-ccsA, psbI-trnS, ndhF-ycf1, ndhI-ndhG and atpA-trnR) containing parsimony-informative characters higher than 2%, which may be potential informative markers for barcoding and phylogenetic analysis. Repeat structure, codon usage and contraction of the IR were also investigated to reveal the pattern of evolution. Phylogenetic analysis demonstrated a sister relationship between A. adenophora and Guizotia abyssinica and supported a monophyly of the Asterales.

Conclusion

We have assembled and analyzed the chloroplast genome of A. adenophora in this study, which was the first sequenced plastome in the Eupatorieae tribe. The complete chloroplast genome information is useful for plant phylogenetic and evolutionary studies within this invasive species and also within the Asteraceae family.

A PPR-DYW protein is required for splicing of a group II intron of cox1 pre-mRNA in Physcomitrella patens

The pentatricopeptide repeat (PPR) protein family is involved in various steps of RNA metabolism in plastids and mitochondria. To investigate the function of a DYW sub-class PPR protein in the moss Physcomitrella patens, we constructed and characterized knockout mutants of the PpPPR_43 gene, which encodes a mitochondrial localized PPR protein with a C-terminal DYW domain. The disruptants showed poor growth of moss protonemata. To investigate whether mitochondrial transcripts were affected by disruption of PpPPR_43, we sequenced the cDNA to detect RNA editing events and performed RT-PCR analyses to measure steady-state mitochondrial transcript levels. Disruption of PpPPR_43 did not result in defective RNA editing, but a substantial reduction in the level of mature cox1 transcript was observed in the disruptants. RT-PCR analysis showed that the 3rd intron of cox1 pre-mRNA was not spliced out in the disruptants, but the 1st, 2nd and 4th introns were efficiently spliced out. This suggests that PpPPR_43 is an intron 3-specific splicing factor. The role of the C-terminal domains of PpPPR_43 in intron 3 splicing was analyzed by complementation experiments with truncated constructs lacking the DYW domain or both the E and DYW domains. Both truncated genes completely restored splicing in the PpPPR_43 knockout mutant. This indicates that the E and DYW domains of PpPPR_43 are not required for splicing, and can be deleted without loss of cox1 intron 3 splicing.

A functional component of the transcriptionally active chromosome complex, Arabidopsis pTAC14, interacts with pTAC12/HEMERA and regulates plastid gene expression

The SET domain-containing protein, pTAC14, was previously identified as a component of the transcriptionally active chromosome (TAC) complexes. Here, we investigated the function of pTAC14 in the regulation of plastid-encoded bacterial-type RNA polymerase (PEP) activity and chloroplast development. The knockout of pTAC14 led to the blockage of thylakoid formation in Arabidopsis (Arabidopsis thaliana), and ptac14 was seedling lethal. Sequence and transcriptional analysis showed that pTAC14 encodes a specific protein in plants that is located in the chloroplast associated with the thylakoid and that its expression depends on light. In addition, the transcript levels of all investigated PEP-dependent genes were clearly reduced in the ptac14-1 mutants, while the accumulation of nucleus-encoded phage-type RNA polymerase-dependent transcripts was increased, indicating an important role of pTAC14 in maintaining PEP activity. pTAC14 was found to interact with pTAC12/HEMERA, another component of TACs that is involved in phytochrome signaling. The data suggest that pTAC14 is essential for proper chloroplast development, most likely by affecting PEP activity and regulating PEP-dependent plastid gene transcription in Arabidopsis together with pTAC12.