Complete structure of the chloroplast genome of Arabidopsis thaliana

Flipping the switch on some of the slowest mutating genomes: Direct measurements of plant mitochondrial and plastid mutation rates in msh1 mutants

Plant mitochondrial and plastid genomes have exceptionally slow rates of sequence evolution, and recent work has identified an unusual member of the MutS gene family ("plant MSH1") as being instrumental in preventing point mutations in these genomes. However, the effects of disrupting MSH1-mediated DNA repair on "germline" mutation rates have not been quantified. Here, we used Arabidopsis thaliana mutation accumulation (MA) lines to measure mutation rates in msh1 mutants and matched wild type (WT) controls. We detected 124 single nucleotide variants (SNVs: 49 mitochondrial and 75 plastid) and 668 small insertions and deletions (indels: 258 mitochondrial and 410 plastid) in msh1 MA lines. In striking contrast, we did not find any organelle mutations in the WT MA lines, and reanalysis of data from a much larger WT MA experiment also failed to detect any variants. The observed number of SNVs in the msh1 MA lines corresponds to estimated mutation rates of 6.1×10-7 and 3.2 ×10-6 per bp per generation in mitochondrial and plastid genomes, respectively. These rates exceed those of species known to have very high mitochondrial mutation rates (e.g., nematodes and fruit flies) by an order of magnitude or more and are on par with estimated rates in humans despite the generation times of A. thaliana being nearly 100-fold shorter. Therefore, disruption of a single plant-specific genetic factor in A. thaliana is sufficient to erase or even reverse the enormous difference in organelle mutation rates between plants and animals.

Targeted C-to-T Base Editing in the Arabidopsis Plastid Genome

Arabidopsis thaliana, particularly the ecotype Columbia-0 (Col-0), has been extensively employed in the study of genetics of the nuclear genome. However, the difficulty of modifying the plastid genome of Col-0, the most widely used ecotype, has hindered investigation of the functional interactions between nuclear-encoded and plastid-encoded genes in this ecotype. Recently, we achieved targeted base editing, substituting a specific C:G pair with a T:A pair in the plastid genome of Col-0 through the application of genome-editing technology. This article introduces the method employed to accomplish this targeted base editing. The process involves four steps: (i) designing and constructing a binary vector encoding the genome-editing enzyme, (ii) introducing the binary vector into the nuclear genome of Col-0 via floral dipping, (iii) identifying base-edited plants, and (iv) verifying inheritance of the edited base(s) by the next generation. © 2025 The Author(s). Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Design and construction of a binary vector encoding ptpTALECD or ptpTALECD_v2 Basic Protocol 2: Agrobacterium-mediated introduction of a binary vector into the Arabidopsis nuclear genome Basic Protocol 3: Selection of plants harboring T-DNA in the nucleus and detection of base editing in the plastid genome.

Comprehensive comparative analysis and development of molecular markers for Lasianthus species based on complete chloroplast genome sequences

BACKGROUND

Lasianthus species are widely used in traditional Chinese folk medicine with high medicinal value. However, source materials and herbarium specimens are often misidentified due to morphological characteristics and commonly used DNA barcode fragments are not sufficient for accurately identifying Lasianthus species. To improve the molecular methods for distinguishing among Lasianthus species, we report the complete chloroplast (CP) genomes of Lasianthus attenuatus, Lasianthus henryi, Lasianthus hookeri, Lasianthus sikkimensis, obtained via high-throughput Illumina sequencing.

RESULTS

These showed CP genomes size of 160164-160246 bp and a typical quadripartite structure, including a large single-copy region (86675-86848 bp), a small single-copy region (17177-17326 bp), and a pair of inverted repeats (28089-28135 bp). As a whole, the gene order, GC content and IR/SC boundary structure were remarkably similar among of the four Lasianthus CP genomes, the partial gene length and IR, LSC and SSC regions length are still different. The average GC content of the CP genomes was 36.71-36.75%, and a total of 129 genes were detected, including 83 different protein-coding genes, 8 different rRNA genes and 38 different tRNA genes. Furthermore, we compared our 4 complete CP genomes data with publicly available CP genome data from six other Lasianthus species, and we initially screened eleven highly variable region fragments were initially screened. We then evaluated the identification efficiency of eleven highly variable region fragments and 5 regular barcode fragments. Ultimately, we found that the optimal combination fragment' ITS2 + psaI-ycf4' could authenticated the Lasianthus species well. Additionally, the results of genome comparison of Rubiaceae species showed that the coding region is more conservative than the non-coding region, and the ycf1 gene shows the most significant variation. Finally, 49 species of CP genome sequences belonging to 16 genera of the Rubiaceae family were used to construct phylogenetic trees.

CONCLUSIONS

Our research is the first to analyze the chloroplast genomes of four species of Lasianthus in detail and we ultimately determined that the combination fragment' ITS2 + psaI-ycf4' is the optimal barcode combination for identifying the genus of Lasianthus. Meanwhile, we gathered the available CP genome sequences from the Rubiaceae and used them to construct the most comprehensive phylogenetic tree for the Rubiaceae family. These investigations provide an important reference point for further studies in the species identification, genetic diversity, and phylogenetic analyses of Rubiaceae species.

Genome editing in angiosperm chloroplasts: targeted DNA double-strand break and base editing

Chloroplasts are organelles that are derived from a photosynthetic bacterium and have their own genome. Genome editing is a recently developing technology that allows for specific modifications of target sequences. The first successful application of genome editing in chloroplasts was reported in 2021, and since then, this research field has been expanding. Although the chloroplast genome of several dicot species can be stably modified by a conventional method, which involves inserting foreign DNAs into the chloroplast genome via homologous recombination, genome editing offers several advantages over this method. In this review, we introduce genome editing methods targeting the chloroplast genome and describe their advantages and limitations. So far, CRISPR/Cas systems are inapplicable for editing the chloroplast genome because guide RNAs, unlike proteins, cannot be efficiently delivered into chloroplasts. Therefore, protein-based enzymes are used to edit the chloroplast genome. These enzymes contain a chloroplast-transit peptide, the DNA-binding domain of transcription activator-like effector nuclease (TALEN), or a catalytic domain that induces DNA modifications. To date, genome editing methods can cause DNA double-strand break or introduce C:G-to-T:A and A:T-to-G:C base edits at or near the target sequence. These methods are expected to contribute to basic research on the chloroplast genome in many species and to be fundamental methods of plant breeding utilizing the chloroplast genome.

The yellow-cotyledon gene (ATYCO) is a crucial factor for thylakoid formation and photosynthesis regulation in Arabidopsis

Chloroplast development underpins plant growth, by facilitating not only photosynthesis but also other essential biochemical processes. Nonetheless, the regulatory mechanisms and functional components of chloroplast development remain largely uncharacterized due to their complexity. In our study, we identified a plastid-targeted gene, ATYCO/RP8/CDB1, as a critical factor in early chloroplast development in Arabidopsis thaliana. YCO knock-out mutant (yco) exhibited a seedling-lethal, albino phenotype, resulting from dysfunctional chloroplasts lacking thylakoid membranes. Conversely, YCO knock-down mutants produced a chlorophyll-deficient cotyledon and normal leaves when supplemented with sucrose. Transcription analysis also revealed that YCO deficiency could be partially compensated by sucrose supplementation, and that YCO played different roles in the cotyledons and the true leaves. In YCO knock-down mutants, the transcript levels of plastid-encoded RNA polymerase (PEP)-dependent genes and nuclear-encoded photosynthetic genes, as well as the accumulation of photosynthetic proteins, were significantly reduced in the cotyledons. Moreover, the chlorophyll-deficient phenotype in YCO knock-down line can be effectively suppressed by inhibition of PSI cyclic electron transport activity, implying an interaction between YCO and PSI cyclic electron transport. Taken together, our findings de underscore the vital role of YCO in early chloroplast development and photosynthesis.

ScRNA-seq reveals the spatiotemporal distribution of camptothecin pathway and transposon activity in Camptotheca acuminata shoot apexes and leaves

Camptotheca acuminata Decne., a significant natural source of the anticancer drug camptothecin (CPT), synthesizes CPT through the monoterpene indole alkaloid (MIA) pathway. In this study, we used single-cell RNA sequencing (scRNA-seq) to generate datasets encompassing over 60,000 cells from C. acuminata shoot apexes and leaves. After cell clustering and annotation, we identified five major cell types in shoot apexes and four in leaves. Analysis of MIA pathway gene expression revealed that most of them exhibited heightened expression in proliferating cells (PCs) and vascular cells (VCs). In contrast to MIA biosynthesis in Catharanthus roseus, CPT biosynthesis in C. acuminata did not exhibit multicellular compartmentalization. Some putative genes encoding enzymes and transcription factors (TFs) related to the biosynthesis of CPT and its derivatives were identified through co-expression analysis. These include 19 cytochrome P450 genes, 8 O-methyltransferase (OMT) genes, and 62 TFs. Additionally, these pathway genes exhibited dynamic expression patterns during VC and EC development. Furthermore, by integrating gene and transposable element (TE) expression data, we constructed novel single-cell transcriptome atlases for C. acuminata. This approach significantly facilitated the identification of rare cell types, including peripheral zone cells (PZs). Some TE families displayed cell type specific, tissue specific, or developmental stage-specific expression patterns, suggesting crucial roles for these TEs in cell differentiation and development. Overall, this study not only provides novel insights into CPT biosynthesis and spatial-temporal TE expression characteristics in C. acuminata, but also serves as a valuable resource for further comprehensive investigations into the development and physiology of this species.

Comparative analysis of 12 water lily plastid genomes reveals genomic divergence and evolutionary relationships in early flowering plants

The aquatic plant Nymphaea, a model genus of the early flowering plant lineage Nymphaeales and family Nymphaeaceae, has been extensively studied. However, the availability of chloroplast genome data for this genus is incomplete, and phylogenetic relationships within the order Nymphaeales remain controversial. In this study, 12 chloroplast genomes of Nymphaea were assembled and analyzed for the first time. These genomes were 158,290-160,042 bp in size and contained 113 non-repeat genes, including 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. We also report on codon usage, RNA editing sites, microsatellite structures, and new repetitive sequences in this genus. Comparative genomics revealed that expansion and contraction of IR regions can lead to changes in the gene numbers. Additionally, it was observed that the highly variable regions of the chloroplast genome were mainly located in intergenic regions. Furthermore, the phylogenetic tree showed the order Nymphaeales was divided into three families, and the genus Nymphaea can be divided into five (or three) subgenera, with the subgenus Nymphaea being the oldest. The divergence times of nymphaealean taxa were analyzed, with origins of the order Nymphaeales and family Nymphaeaceae being about 194 and 131 million years, respectively. The results of the phylogenetic analysis and estimated divergence times will be useful for future evolutionary studies of basal angiosperm lineages.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-024-00242-0.

Genetic characterization and phylogenetic analysis of the Nigella sativa (black seed) plastome

In this study, the complete plastome sequence of Nigella sativa (black seed), was analyzed for the first time. The plastome spans approximately 154,120 bp, comprising four sections: the Large Single-Copy (LSC) (85,538 bp), the Small Single-Copy (SSC) (17,984 bp), and two Inverted Repeat (IR) regions (25,299 bp). A comparative study of N. sativa's plastome with ten other species from various genera in the Ranunculaceae family reveals substantial structural variations. The contraction of the inverted repeat region in N. sativa influences the boundaries of single-copy regions, resulting in a shorter plastome size than other species. When comparing the plastome of N. sativa with those of its related species, significant divergence is observed, particularly except for N. damascena. Among these, the plastome of A. glaucifolium displays the highest average pairwise sequence divergence (0.2851) with N. sativa, followed by A. raddeana (0.2290) and A. coerulea (0.1222). Furthermore, the study identified 12 distinct hotspot regions characterized by elevated Pi values (> 0.1). These regions include trnH-GUG-psbA, matK-trnQ-UUG, psbK-trnR-UCU, atpF-atpI, rpoB-psbD, ycf3-ndhJ, ndhC-cemA, petA-psaJ, trnN-GUU-ndhF, trnV-GAC-rps12, ycf2-trnI-CAU, and ndhA-ycf1. Approximately, 24 tandem and 48 palindromic and forward repeats were detected in N. sativa plastome. The analysis revealed 32 microsatellites with the majority being mononucleotide repeats. In the N. sativa plastome, phenylalanine had the highest number of codons (1982 codons), while alanine was the least common amino acid with 260 codons. A phylogenetic tree, constructed using protein-coding genes, revealed a distinct monophyletic clade comprising N. sativa and N. damascene, closely aligned with the Cimicifugeae tribe and exhibiting robust support. This plastome provides valuable genetic information for precise species identification, phylogenetic resolution, and evolutionary studies of N. sativa.

Characterization of organelle DNA degradation mediated by DPD1 exonuclease in the rice genome-edited line

Mitochondria and plastids, originated as ancestral endosymbiotic bacteria, contain their own DNA sequences. These organelle DNAs (orgDNAs) are, despite the limited genetic information they contain, an indispensable part of the genetic systems but exist as multiple copies, making up a substantial amount of total cellular DNA. Given this abundance, orgDNA is known to undergo tissue-specific degradation in plants. Previous studies have shown that the exonuclease DPD1, conserved among seed plants, degrades orgDNAs during pollen maturation and leaf senescence in Arabidopsis. However, tissue-specific orgDNA degradation was shown to differ among species. To extend our knowledge, we characterized DPD1 in rice in this study. We created a genome-edited (GE) mutant in which OsDPD1 and OsDPD1-like were inactivated. Characterization of this GE plant demonstrated that DPD1 was involved in pollen orgDNA degradation, whereas it had no significant effect on orgDNA degradation during leaf senescence. Comparison of transcriptomes from wild-type and GE plants with different phosphate supply levels indicated that orgDNA had little impact on the phosphate starvation response, but instead had a global impact in plant growth. In fact, the GE plant showed lower fitness with reduced grain filling rate and grain weight in natural light conditions. Taken together, the presented data reinforce the important physiological roles of orgDNA degradation mediated by DPD1.

Chloroplast Ribosome Biogenesis Factors

The formation of chloroplasts can be traced back to an ancient event in which a eukaryotic host cell containing mitochondria ingested a cyanobacterium. Since then, chloroplasts have retained many characteristics of their bacterial ancestor, including their transcription and translation machinery. In this review, recent research on the maturation of rRNA and ribosome assembly in chloroplasts is explored, along with their crucial role in plant survival and their implications for plant acclimation to changing environments. A comparison is made between the ribosome composition and auxiliary factors of ancient and modern chloroplasts, providing insights into the evolution of ribosome assembly factors. Although the chloroplast contains ancient proteins with conserved functions in ribosome assembly, newly evolved factors have also emerged to help plants acclimate to changes in their environment and internal signals. Overall, this review offers a comprehensive analysis of the molecular mechanisms underlying chloroplast ribosome assembly and highlights the importance of this process in plant survival, acclimation and adaptation.

Nitrogen-fixing organelle in a marine alga

Symbiotic interactions were key to the evolution of chloroplast and mitochondria organelles, which mediate carbon and energy metabolism in eukaryotes. Biological nitrogen fixation, the reduction of abundant atmospheric nitrogen gas (N2) to biologically available ammonia, is a key metabolic process performed exclusively by prokaryotes. Candidatus Atelocyanobacterium thalassa, or UCYN-A, is a metabolically streamlined N2-fixing cyanobacterium previously reported to be an endosymbiont of a marine unicellular alga. Here we show that UCYN-A has been tightly integrated into algal cell architecture and organellar division and that it imports proteins encoded by the algal genome. These are characteristics of organelles and show that UCYN-A has evolved beyond endosymbiosis and functions as an early evolutionary stage N2-fixing organelle, or "nitroplast."

The newly assembled chloroplast genome of Aeluropus littoralis: molecular feature characterization and phylogenetic analysis with related species

Aeluropus littoralis, a halophyte grass, is widely distributed from the Mediterranean to the Indian subcontinent through the Mongolian Gobi. This model halophyte has garnered increasing attention owing to its use as forage and its high tolerance to environmental stressors. The chloroplast genomes of many plants have been extensively examined for molecular, phylogenetic and transplastomic applications. However, no published research on the A. littoralis chloroplast (cp) genome was discovered. Here, the entire chloroplast genome of A. littoralis was assembled implementing accurate long-read sequences. The entire chloroplast genome, with an estimated length of 135,532 bp (GC content: 38.2%), has a quadripartite architecture and includes a pair of inverted repeat (IR) regions, IRa and IRb (21,012 bp each), separated by a large and a small single-copy regions (80,823 and 12,685 bp, respectively). The features of A. littoralis consist of 133 genes that synthesize 87 peptides, 38 transfer RNAs, and 8 ribosomal RNAs. Of these genes, 86 were unique, whereas 19 were duplicated in IR regions. Additionally, a total of forty-six simple sequence repeats, categorized into 32-mono, four-di, two-tri, and eight-tetranucleotides, were discovered. Furthermore, ten sets of repeats greater than 20 bp were located primarily in the LSC region. Evolutionary analysis based on chloroplast sequence data revealed that A. littoralis with A. lagopoides and A. sinensis belong to the Aeluropodinae subtribe, which is a sister to the Eleusininae in the tribe Cynodonteae and the subfamily Chloridoideae. This subfamily belongs to the PACMAD clade, which contains the majority of the C4 photosynthetic plants in the Poaceae. The newly constructed A. littoralis cp genome offers valuable knowledge for DNA barcoding, phylogenetic, transplastomic research, and other biological studies.

The complete chloroplast genome sequence of an invasive plant, Tragopogon dubius Scopoli (asteraceae)

Tragopogon dubius Scopoli is native to Europe and western Asia and is considered an invasive plant in China. In this study, the complete chloroplast genome of T. dubius was obtained using high-throughput next-generation sequencing technology. The whole chloroplast genome was 153,017 bp long with a GC content of 38% and comprised 130 genes (86 protein-coding genes, 36 tRNA genes, and 8 rRNA genes). Phylogenetic analysis based on the concatenated chloroplast protein-coding sequences showed that T. dubius is most closely related to Tragopogon pratensis. This study provides valuable genetic data for further phylogenetic analysis and molecular identification of species in the genus Tragopogon.

Detection and editing of the updated Arabidopsis plastid- and mitochondrial-encoded proteomes through PeptideAtlas

Arabidopsis (Arabidopsis thaliana) ecotype Col-0 has plastid and mitochondrial genomes encoding over 100 proteins. Public databases (e.g. Araport11) have redundancy and discrepancies in gene identifiers for these organelle-encoded proteins. RNA editing results in changes to specific amino acid residues or creation of start and stop codons for many of these proteins, but the impact of RNA editing at the protein level is largely unexplored due to the complexities of detection. Here, we assembled the nonredundant set of identifiers, their correct protein sequences, and 452 predicted nonsynonymous editing sites of which 56 are edited at lower frequency. We then determined accumulation of edited and/or unedited proteoforms by searching ∼259 million raw tandem MS spectra from ProteomeXchange, which is part of PeptideAtlas (www.peptideatlas.org/builds/arabidopsis/). We identified all mitochondrial proteins and all except 3 plastid-encoded proteins (NdhG/Ndh6, PsbM, and Rps16), but no proteins predicted from the 4 ORFs were identified. We suggest that Rps16 and 3 of the ORFs are pseudogenes. Detection frequencies for each edit site and type of edit (e.g. S to L/F) were determined at the protein level, cross-referenced against the metadata (e.g. tissue), and evaluated for technical detection challenges. We detected 167 predicted edit sites at the proteome level. Minor frequency sites were edited at low frequency at the protein level except for cytochrome C biogenesis 382 at residue 124 (Ccb382-124). Major frequency sites (>50% editing of RNA) only accumulated in edited form (>98% to 100% edited) at the protein level, with the exception of Rpl5-22. We conclude that RNA editing for major editing sites is required for stable protein accumulation.

The complete plastome sequences of invasive weed Parthenium hysterophorus: genome organization, evolutionary significance, structural features, and comparative analysis

Parthenium hysterophorus, a globally widespread weed, poses a significant threat to agricultural ecosystems due to its invasive nature. We investigated the chloroplast genome of P. hysterophorus in this study. Our analysis revealed that the chloroplast genome of P. hysterophorus spans a length of 151,881 base pairs (bp). It exhibits typical quadripartite structure commonly found in chloroplast genomes, including inverted repeat regions (IR) of 25,085 bp, a small single copy (SSC) region of 18,052 bp, and a large single copy (LSC) region of 83,588 bp. A total of 129 unique genes were identified in P. hysterophorus chloroplast genomes, including 85 protein-coding genes, 36 tRNAs, and eight rRNAs genes. Comparative analysis of the P. hysterophorus plastome with those of related species from the tribe Heliantheae revealed both conserved structures and intriguing variations. While many structural elements were shared among the species, we identified a rearrangement in the large single-copy region of P. hysterophorus. Moreover, our study highlighted notable gene divergence in several specific genes, namely matK, ndhF, clpP, rps16, ndhA, rps3, and ndhD. Phylogenetic analysis based on the 72 shared genes placed P. hysterophorus in a distinct clade alongside another species, P. argentatum. Additionally, the estimated divergence time between the Parthenium genus and Helianthus (sunflowers) was approximately 15.1 million years ago (Mya). These findings provide valuable insights into the evolutionary history and genetic relationships of P. hysterophorus, shedding light on its divergence and adaptation over time.

CLAW: An automated Snakemake workflow for the assembly of chloroplast genomes from long-read data

Chloroplasts are photosynthetic organelles in algal and plant cells that contain their own genome. Chloroplast genomes are commonly used in evolutionary studies and taxonomic identification and are increasingly becoming a target for crop improvement studies. As DNA sequencing becomes more affordable, researchers are collecting vast swathes of high-quality whole-genome sequence data from laboratory and field settings alike. Whole tissue read libraries sequenced with the primary goal of understanding the nuclear genome will inadvertently contain many reads derived from the chloroplast genome. These whole-genome, whole-tissue read libraries can additionally be used to assemble chloroplast genomes with little to no extra cost. While several tools exist that make use of short-read second generation and third-generation long-read sequencing data for chloroplast genome assembly, these tools may have complex installation steps, inadequate error reporting, poor expandability, and/or lack scalability. Here, we present CLAW (Chloroplast Long-read Assembly Workflow), an easy to install, customise, and use Snakemake tool to assemble chloroplast genomes from chloroplast long-reads found in whole-genome read libraries (https://github.com/aaronphillips7493/CLAW). Using 19 publicly available reference chloroplast genome assemblies and long-read libraries from algal, monocot and eudicot species, we show that CLAW can rapidly produce chloroplast genome assemblies with high similarity to the reference assemblies. CLAW was designed such that users have complete control over parameterisation, allowing individuals to optimise CLAW to their specific use cases. We expect that CLAW will provide researchers (with varying levels of bioinformatics expertise) with an additional resource useful for contributing to the growing number of publicly available chloroplast genome assemblies.

PEP-ASSOCIATED PROTEIN 3 regulates rice tiller formation and grain yield by controlling chloroplast biogenesis

Plastid-encoded RNA polymerase (PEP) plays a pivotal role in chloroplast development by governing the transcription of chloroplast genes, and PEP-associated proteins (PAPs) modulate PEP transcriptional activity. Therefore, PAPs provide an intriguing target for those efforts to improve yield, by enhancing chloroplast development. In this study, we identified the rice (Oryza sativa) OsPAP3 gene and characterized its function in chloroplast development. OsPAP3 expression was light-dependent and leaf-specific, similar to the PEP-dependent chloroplast gene RUBISCO LARGE SUBUNIT (OsRbcL), and OsPAP3 protein localized to chloroplast nucleoids where PEP functions. Analysis of loss-of-function and gain-of-function mutants showed that the expression of OsPAP3 is tightly linked to chloroplast gene expression and chloroplast biogenesis in rice. Homozygous knockout mutants of OsPAP3 had fewer chloroplasts than wild type, whereas plants overexpressing OsPAP3 had more chloroplasts. Also, OsPAP3 knockout suppressed the PEP-dependent expression of chloroplast genes, but OsPAP3 overexpression increased their expression. These findings indicate that OsPAP3 regulates chloroplast biogenesis in rice by controlling the PEP-dependent expression of chloroplast genes. More importantly, data from 3 seasons of field cultivation revealed that the overexpression of OsPAP3 improves rice grain yield by approximately 25%, largely due to increased tiller formation. Collectively, these observations suggest that OsPAP3 regulates rice growth and productivity by promoting chloroplast development.

Characterization of the complete plastome sequence of Korean endemic, Cardamine glechomifolia H.Lév. (Brassicaceae, Brassicales)

In this study, we report the complete plastome sequence of Cardamine glechomifolia H.Lév. 1913 (NCBI acc. no. OP894664). This plastome shows typical quadripartite structure. The plastome size is 154,307 bp, which consists of 84,015 bp large single-copy (LSC), 17,690 bp small single-copy (SSC), and 26,301 bp inverted repeat (IR) regions. The plastome contains 112 genes, including 78 protein-coding, 30 tRNA, and four rRNA genes. The infA gene is pseudogenized. Sixteen genes contain one intron and two genes (clpP and ycf3) have two introns. The phylogenomic analysis conducted in our study reveals that the genus Cardamine, which encompasses C. glechomifolia, exhibits three distinct clades. In order to elucidate the interrelationship among the three clades, it is imperative to conduct additional investigations by augmenting the number of Cardamine samples.

Hierarchy of codon usage frequencies from codon-anticodon interaction in the crystal basis model

Analyzing the codon usage frequencies of a specimen of 20 plants, for which the codon-anticodon pattern is known, we have remarked that the hierarchy of the usage frequencies present an almost "universal" behavior. Searching to explain this behavior, we assume that the codon usage probability results from the sum of two contributions: the first dominant term is an almost "universal" one and it depends on the codon-anticodon interaction; the second term is a local one, i.e. depends on the biological species. The codon-anticodon interaction is written as a spin-spin plus a z-spin term in the formalism of the crystal basis model. From general considerations, in particular from the choice of the signs and some constraints on the parameters defining the interaction, we are able to explain most of the observed data.

Interaction of PALE CRESS with PAP2/pTAC2 and PAP3/pTAC10 affects the accumulation of plastid-encoded RNA polymerase complexes in Arabidopsis

The transcription of photosynthesis genes in chloroplasts is largely mediated by the plastid-encoded RNA polymerase (PEP), which resembles prokaryotic-type RNA polymerases, but with plant-specific accessory subunits known as plastid transcriptionally active chromosome proteins (pTACs) or PEP-associated proteins (PAPs). However, whether additional factors are involved in the biogenesis of PEP complexes remains unknown. Here, we investigated the function of an essential gene, PALE CRESS (PAC), in the accumulation of PEP complexes in chloroplasts. We established that an Arabidopsis leaf variegation mutant, variegated 6-1 (var6-1), is a hypomorphic allele of PAC. Unexpectedly, we revealed that a fraction of VAR6/PAC is associated with thylakoid membranes, where it interacts with PEP complexes. The accumulation of PEP complexes is defective in both var6-1 and the null allele var6-2. Further protein interaction assays confirmed that VAR6/PAC interacts directly with the PAP2/pTAC2 and PAP3/pTAC10 subunits of PEP complexes. Moreover, we generated viable hypomorphic alleles of the essential gene PAP2/pTAC2, and revealed a genetic interaction between PAC and PAP2/pTAC2 in photosynthesis gene expression and PEP complex accumulation. Our findings establish that VAR6/PAC affects PEP complex accumulation through interactions with PAP2/pTAC2 and PAP3/pTAC10, and provide new insights into the accumulation of PEP and chloroplast development.

Complete chloroplast genomes of Cerastium alpinum, C. arcticum and C. nigrescens: genome structures, comparative and phylogenetic analysis

The genus Cerastium includes about 200 species that are mostly found in the temperate climates of the Northern Hemisphere. Here we report the complete chloroplast genomes of Cerastium alpinum, C. arcticum and C. nigrescens. The length of cp genomes ranged from 147,940 to 148,722 bp. Their quadripartite circular structure had the same gene organization and content, containing 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. Repeat sequences varied from 16 to 23 per species, with palindromic repeats being the most frequent. The number of identified SSRs ranged from 20 to 23 per species and they were mainly composed of mononucleotide repeats containing A/T units. Based on Ka/Ks ratio values, most genes were subjected to purifying selection. The newly sequenced chloroplast genomes were characterized by a high frequency of RNA editing, including both C to U and U to C conversion. The phylogenetic relationships within the genus Cerastium and family Caryophyllaceae were reconstructed based on the sequences of 71 protein-coding genes. The topology of the phylogenetic tree was consistent with the systematic position of the studied species. All representatives of the genus Cerastium were gathered in a single clade with C. glomeratum sharing the least similarity with the others.

The tomato chloroplast stromal proteome compendium elucidated by leveraging a plastid protein-localization prediction Atlas

Tomato (Solanum lycopersicum) is a model species for studying fruit development, wounding, herbivory, and pathogen attack. Despite tomato's world-wide economic importance and the role of chloroplasts as metabolic hubs and integrators of environmental cues, little is known about the stromal proteome of tomato. Using a high-yielding protocol for chloroplast and stromal protein isolation, MudPIT nano-LC-MS/MS analyses, a robust in-house protein database (the Atlas) for predicting the plastid localization of tomato proteins, and rigorous selection criteria for inclusion/exclusion in the stromal proteome, we identified 1,278 proteins of the tomato stromal proteome. We provide one of the most robust stromal proteomes available to date with empirical evidence for 545 and 92 proteins not previously described for tomato plastids and the Arabidopsis stroma, respectively. The relative abundance of tomato stromal proteins was determined using the exponentially modified protein abundance index (emPAI). Comparison of the abundance of tomato and Arabidopsis stromal proteomes provided evidence for the species-specific nature of stromal protein homeostasis. The manual curation of the tomato stromal proteome classified proteins into ten functional categories resulting in an accessible compendium of tomato chloroplast proteins. After curation, only 91 proteins remained as unknown, uncharacterized or as enzymes with unknown functions. The curation of the tomato stromal proteins also indicated that tomato has a number of paralogous proteins, not present in Arabidopsis, which accumulated to different levels in chloroplasts. As some of these proteins function in key metabolic pathways or in perceiving or transmitting signals critical for plant adaptation to biotic and abiotic stress, these data suggest that tomato may modulate the bidirectional communication between chloroplasts and nuclei in a novel manner. The stromal proteome provides a fertile ground for future mechanistic studies in the field of tomato chloroplast-nuclear signaling and are foundational for our goal of elucidating the dynamics of the stromal proteome controlled by the solanaceous-specific, stromal, and wound-inducible leucine aminopeptidase A of tomato.

The complete chloroplast genome of Salix lindleyana (salicaceae), a plateau plant species

Salix lindleyana Wallich ex Andersson 1851 is a species of genus Salix which mainly grows on mountains above 3000 m at sea level in Qinghai-Tibetan Plateau (including the Himalayas and Hengduan Mountains). To determine its phylogenetic position within Salix, we reconstructed S. lindleyana complete chloroplast (cp) genome sequence by de novo assembly using whole-genome sequencing data. The completed chloroplast genome was 155,304 bp, with a total GC content of 36.7%. It had a very typical tetrad structure, including a large single-copy (LSC) region of 84,539 bp, a small single-copy (SSC) region of 16,161 bp, and two inverted repeats (IR) regions of 27,302 bp. A total of 132 functional genes were distributed in the chloroplast genome, including 87 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. Phylogenetic analysis showed that S. lindleyana was clustered with Salix dasyclados Wimmer 1849 and Salix variegata Franchet 1887. The complete chloroplast genome of S. lindleyana provides potential genetic resources for further phylogenetic studies.

Phylogenomics and plastome evolution of Indigofera (Fabaceae)

Introduction

Indigofera L. is the third largest genus in Fabaceae and includes economically important species that are used for indigo dye-producing, medicinal, ornamental, and soil and water conservation. The genus is taxonomically difficult due to the high level of overlap in morphological characters of interspecies, fewer reliability states for classification, and extensive adaptive evolution. Previous characteristic-based taxonomy and nuclear ITS-based phylogenies have contributed to our understanding of Indigofera taxonomy and evolution. However, the lack of chloroplast genomic resources limits our comprehensive understanding of the phylogenetic relationships and evolutionary processes of Indigofera.

Methods

Here, we newly assembled 18 chloroplast genomes of Indigofera. We performed a series of analyses of genome structure, nucleotide diversity, phylogenetic analysis, species pairwise Ka/Ks ratios, and positive selection analysis by combining with allied species in Papilionoideae.

Results and discussion

The chloroplast genomes of Indigofera exhibited highly conserved structures and ranged in size from 157,918 to 160,040 bp, containing 83 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Thirteen highly variable regions were identified, of which trnK-rbcL, ndhF-trnL, and ycf1 were considered as candidate DNA barcodes for species identification of Indigofera. Phylogenetic analysis using maximum likelihood (ML) and Bayesian inference (BI) methods based on complete chloroplast genome and protein-coding genes (PCGs) generated a well-resolved phylogeny of Indigofera and allied species. Indigofera monophyly was strongly supported, and four monophyletic lineages (i.e., the Pantropical, East Asian, Tethyan, and Palaeotropical clades) were resolved within the genus. The species pairwise Ka/Ks ratios showed values lower than 1, and 13 genes with significant posterior probabilities for codon sites were identified in the positive selection analysis using the branch-site model, eight of which were associated with photosynthesis. Positive selection of accD suggested that Indigofera species have experienced adaptive evolution to selection pressures imposed by their herbivores and pathogens. Our study provided insight into the structural variation of chloroplast genomes, phylogenetic relationships, and adaptive evolution in Indigofera. These results will facilitate future studies on species identification, interspecific and intraspecific delimitation, adaptive evolution, and the phylogenetic relationships of the genus Indigofera.

Peptide chain release factor DIG8 regulates plant growth by affecting ROS-mediated sugar transportation in Arabidopsis thaliana

Chloroplasts have important roles in photosynthesis, stress sensing and retrograde signaling. However, the relationship between chloroplast peptide chain release factor and ROS-mediated plant growth is still unclear. In the present study, we obtained a loss-of-function mutant dig8 by EMS mutation. The dig8 mutant has few lateral roots and a pale green leaf phenotype. By map-based cloning, the DIG8 gene was located on AT3G62910, with a point mutation leading to amino acid substitution in functional release factor domain. Using yeast-two-hybrid and BiFC, we confirmed DIG8 protein was characterized locating in chloroplast by co-localization with plastid marker and interacting with ribosome-related proteins. Through observing by transmission electron microscopy, quantifying ROS content and measuring the transport efficiency of plasmodesmata in dig8 mutant, we found that abnormal thylakoid stack formation and chloroplast dysfunction in the dig8 mutant caused increased ROS activity leading to callose deposition and lower PD permeability. A local sugar supplement partially alleviated the growth retardation phenotype of the mutant. These findings shed light on chloroplast peptide chain release factor-affected plant growth by ROS stress.

A Single Nucleotide Variation of CRS2 Affected the Establishment of Photosynthetic System in Rice

Chloroplasts are essential sites for plant photosynthesis, and the biogenesis of the photosynthetic complexes involves the interaction of nuclear genes and chloroplast genes. In this study, we identified a rice pale green leaf mutant, crs2. The crs2 mutant showed different degrees of low chlorophyll phenotypes at different growth stages, especially at the seedling stage. Fine mapping and DNA sequencing of crs2 revealed a single nucleotide substitution (G4120A) in the eighth exons of CRS2, causing a G-to-R mutation of the 229th amino acid of CRS2 (G229R). The results of complementation experiments confirmed that this single-base mutation in crs2 is responsible for the phenotype of the crs2 mutant. CRS2 encodes a chloroplast RNA splicing 2 protein localized in the chloroplast. Western blot results revealed an abnormality in the abundance of the photosynthesis-related protein in crs2. However, the mutation of CRS2 leads to the enhancement of antioxidant enzyme activity, which could reduce ROS levels. Meanwhile, with the release of Rubisco activity, the photosynthetic performance of crs2 was improved. In summary, the G229R mutation in CRS2 causes chloroplast protein abnormalities and affects photosystem performance in rice; the above findings facilitate the elucidation of the physiological mechanism of chloroplast proteins affecting photosynthesis.

Convergent evolution and horizontal gene transfer in Arctic Ocean microalgae

Microbial communities in the world ocean are affected strongly by oceanic circulation, creating characteristic marine biomes. The high connectivity of most of the ocean makes it difficult to disentangle selective retention of colonizing genotypes (with traits suited to biome specific conditions) from evolutionary selection, which would act on founder genotypes over time. The Arctic Ocean is exceptional with limited exchange with other oceans and ice covered since the last ice age. To test whether Arctic microalgal lineages evolved apart from algae in the global ocean, we sequenced four lineages of microalgae isolated from Arctic waters and sea ice. Here we show convergent evolution and highlight geographically limited HGT as an ecological adaptive force in the form of PFAM complements and horizontal acquisition of key adaptive genes. Notably, ice-binding proteins were acquired and horizontally transferred among Arctic strains. A comparison with Tara Oceans metagenomes and metatranscriptomes confirmed mostly Arctic distributions of these IBPs. The phylogeny of Arctic-specific genes indicated that these events were independent of bacterial-sourced HGTs in Antarctic Southern Ocean microalgae.

A comparative study across the parasitic plants of Cuscuta subgenus Grammica (Convolvulaceae) reveals a possible loss of the plastid genome in its section Subulatae

Most species in Cuscuta subgenus Grammica retain many photosynthesis-related plastid genes, generally under purifying selection. A group of holoparasitic species in section Subulatae may have lost their plastid genomes entirely. The c. 153 species of plants belonging to Cuscuta subgenus Grammica are all obligate stem parasites. However, some have completely lost the ability to conduct photosynthesis while others retain photosynthetic machinery and genes. The plastid genome that primarily encodes key photosynthesis genes functions as a bellwether for how reliant plants are on primary production. This research assembles and analyses 17 plastomes across Cuscuta subgenus Grammica with the aim of characterizing the state of the plastome in each of its sections. By comparing the structure and content of plastid genomes across the subgenus, as well as by quantifying the selection acting upon each gene, we reconstructed the patterns of plastome change within the phylogenetic context for this group. We found that species in 13 of the 15 sections that comprise Grammica retain the bulk of plastid photosynthesis genes and are thus hemiparasitic. The complete loss of photosynthesis can be traced to two clades: the entire section Subulatae and a complex of three species within section Ceratophorae. We were unable to recover any significant plastome sequences from section Subulatae, suggesting that plastomes in these species are either drastically reduced or lost entirely.

Complete chloroplast genome molecular structure, comparative and phylogenetic analyses of Sphaeropteris lepifera of Cyatheaceae family: a tree fern from China

Sphaeropteris lepifera is a tree fern in the Cyatheaceae, a family that has played an important role in the evolution of plant systems. This study aimed to analyze the complete chloroplast genome of S. lepifera and compared it with previously published chloroplast genomes Cyatheaceae family. The chloroplast genome of S. lepifera comprised 162,114 bp, consisting of a large single copy (LSC) region of 86,327 bp, a small single copy (SSC) region of 27,731 bp and a pair of inverted repeats (IRa and IRb) of 24,028 bp each. The chloroplast genome encoded 129 genes, comprising 32 transfer RNAs, 8 ribosomal RNAs, and 89 protein-coding genes. Comparison of the genomes of 7 Cyatheaceae plants showed that the chloroplast genome of S. lepifera was missing the gene trnV-UAC. Expansion of the SSC region led to the difference in the chloroplast genome size of S. lepifera. Eight genes, atpI, ccsA, petA, psaB, rpl16, rpoA, rpoC1, and ycf2 have high nucleic acid diversity and can be regarded as potential molecular markers. The genes trnG-trnR and atpB were suitable for DNA barcodes between different communities of S. lepifera. The S. lepifera groups in Zhejiang Province probably diffused from Pingtan and Ningde, Fujian. The results will provide a basis for species identification, biological studies, and endangerment mechanism of S. lepifera.

A comparison of 25 complete chloroplast genomes between sister mangrove species Kandelia obovata and Kandelia candel geographically separated by the South China Sea

In 2003, Kandelia obovata was identified as a new mangrove species differentiated from Kandelia candel. However, little is known about their chloroplast (cp) genome differences and their possible ecological significance. In this study, 25 whole cp genomes, with seven samples of K. candel from Malaysia, Thailand, and Bangladesh and 18 samples of K. obovata from China, were sequenced for comparison. The cp genomes of both species encoded 128 genes, namely 83 protein-coding genes, 37 tRNA genes, and eight rRNA genes, but the cp genome size of K. obovata was ~2 kb larger than that of K. candle due to the presence of more and longer repeat sequences. Of these, tandem repeats and simple sequence repeats exhibited great differences. Principal component analysis based on indels, and phylogenetic tree analyses constructed with homologous protein genes from the single-copy genes, as well as 38 homologous pair genes among 13 mangrove species, gave strong support to the separation of the two species within the Kandelia genus. Homologous genes ndhD and atpA showed intraspecific consistency and interspecific differences. Molecular dynamics simulations of their corresponding proteins, NAD(P)H dehydrogenase chain 4 (NDH-D) and ATP synthase subunit alpha (ATP-A), predicted them to be significantly different in the functions of photosynthetic electron transport and ATP generation in the two species. These results suggest that the energy requirement was a pivotal factor in their adaptation to differential environments geographically separated by the South China Sea. Our results also provide clues for future research on their physiological and molecular adaptation mechanisms to light and temperature.

Anterograde signaling controls plastid transcription via sigma factors separately from nuclear photosynthesis genes

Light initiates chloroplast biogenesis in Arabidopsis by eliminating PHYTOCHROME-INTERACTING transcription FACTORs (PIFs), which in turn de-represses nuclear photosynthesis genes, and synchronously, generates a nucleus-to-plastid (anterograde) signal that activates the plastid-encoded bacterial-type RNA polymerase (PEP) to transcribe plastid photosynthesis genes. However, the identity of the anterograde signal remains frustratingly elusive. The main challenge has been the difficulty to distinguish regulators from the plethora of necessary components for plastid transcription and other essential chloroplast functions, such as photosynthesis. Here, we show that the genome-wide induction of nuclear photosynthesis genes is insufficient to activate the PEP. PEP inhibition is imposed redundantly by multiple PIFs and requires PIF3's activator activity. Among the nuclear-encoded components of the PEP holoenzyme, we identify four light-inducible, PIF-repressed sigma factors as anterograde signals. Together, our results elucidate that light-dependent inhibition of PIFs activates plastid photosynthesis genes via sigma factors as anterograde signals in parallel with the induction of nuclear photosynthesis genes.

PlantRNA 2.0: an updated database dedicated to tRNAs of photosynthetic eukaryotes

PlantRNA (http://plantrna.ibmp.cnrs.fr/) is a comprehensive database of transfer RNA (tRNA) gene sequences retrieved from fully annotated nuclear, plastidial and mitochondrial genomes of photosynthetic organisms. In the first release (PlantRNA 1.0), tRNA genes from 11 organisms were annotated. In this second version, the annotation was implemented to 51 photosynthetic species covering the whole phylogenetic tree of photosynthetic organisms, from the most basal group of Archeplastida, the glaucophyte Cyanophora paradoxa, to various land plants. tRNA genes from lower photosynthetic organisms such as streptophyte algae or lycophytes as well as extremophile photosynthetic species such as Eutrema parvulum were incorporated in the database. As a whole, about 37 000 tRNA genes were accurately annotated. In the frame of the tRNA genes annotation from the genome of the Rhodophyte Chondrus crispus, non-canonical splicing sites in the D- or T-regions of tRNA molecules were identified and experimentally validated. As for PlantRNA 1.0, comprehensive biological information including 5'- and 3'-flanking sequences, A and B box sequences, region of transcription initiation and poly(T) transcription termination stretches, tRNA intron sequences and tRNA mitochondrial import are included.

Ascorbate peroxidase postcold regulation of chloroplast NADPH dehydrogenase activity controls cold memory

Exposure of Arabidopsis (Arabidopsis thaliana) to 4°C imprints a cold memory that modulates gene expression in response to a second (triggering) stress stimulus applied several days later. Comparison of plastid transcriptomes of cold-primed and control plants directly before they were exposed to the triggering stimulus showed downregulation of several subunits of chloroplast NADPH dehydrogenase (NDH) and regulatory subunits of ATP synthase. NDH is, like proton gradient 5 (PGR5)-PGR5-like1 (PGRL1), a thylakoid-embedded, ferredoxin-dependent plastoquinone reductase that protects photosystem I and stabilizes ATP synthesis by cyclic electron transport (CET). Like PGRL1A and PGRL1B transcript levels, ndhA and ndhD transcript levels decreased during the 24-h long priming cold treatment. PGRL1 transcript levels were quickly reset in the postcold phase, but expression of ndhA remained low. The transcript abundances of other ndh genes decreased within the next days. Comparison of thylakoid-bound ascorbate peroxidase (tAPX)-free and transiently tAPX-overexpressing or tAPX-downregulating Arabidopsis lines demonstrated that ndh expression is suppressed by postcold induction of tAPX. Four days after cold priming, when tAPX protein accumulation was maximal, NDH activity was almost fully lost. Lack of the NdhH-folding chaperonin Crr27 (Cpn60β4), but not lack of the NDH activity modulating subunits NdhM, NdhO, or photosynthetic NDH subcomplex B2 (PnsB2), strengthened priming regulation of zinc finger of A. thaliana 10, which is a nuclear-localized target gene of the tAPX-dependent cold-priming pathway. We conclude that cold-priming modifies chloroplast-to-nucleus stress signaling by tAPX-mediated suppression of NDH-dependent CET and that plastid-encoded NdhH, which controls subcomplex A assembly, is of special importance for memory stabilization.

Sorting of mitochondrial and plastid heteroplasmy in Arabidopsis is extremely rapid and depends on MSH1 activity

The fate of new mitochondrial and plastid mutations depends on their ability to persist and spread among the numerous organellar genome copies within a cell (heteroplasmy). The extent to which heteroplasmies are transmitted across generations or eliminated through genetic bottlenecks is not well understood in plants, in part because their low mutation rates make these variants so infrequent. Disruption of MutS Homolog 1 (MSH1), a gene involved in plant organellar DNA repair, results in numerous de novo point mutations, which we used to quantitatively track the inheritance of single nucleotide variants in mitochondrial and plastid genomes in Arabidopsis. We found that heteroplasmic sorting (the fixation or loss of a variant) was rapid for both organelles, greatly exceeding rates observed in animals. In msh1 mutants, plastid variants sorted faster than those in mitochondria and were typically fixed or lost within a single generation. Effective transmission bottleneck sizes (N) for plastids and mitochondria were N ∼ 1 and 4, respectively. Restoring MSH1 function further increased the rate of heteroplasmic sorting in mitochondria (N ∼ 1.3), potentially because of its hypothesized role in promoting gene conversion as a mechanism of DNA repair, which is expected to homogenize genome copies within a cell. Heteroplasmic sorting also favored GC base pairs. Therefore, recombinational repair and gene conversion in plant organellar genomes can potentially accelerate the elimination of heteroplasmies and bias the outcome of this sorting process.

Characterization of Codonopsis pilosula subsp. tangshen plastome and comparative analysis of Codonopsis species

Codonopsis pilosula subsp. tangshen is one of the most important medicinal herbs used in traditional Chinese medicine. Correct identification of materials from C. pilosula subsp. tangshen is critical to ensure the efficacy and safety of the associated medicines. Traditional DNA molecular markers could distinguish Codonopsis species well, so we need to develop super or specific molecular markers. In this study, we reported the plastome of Codonopsis pilosula subsp. tangshen (Oliv.) D.Y. Hong conducted phylogenomic and comparative analyses in the Codonopsis genus for the first time. The entire length of the Codonopsis pilosula subsp. tangshen plastome was 170,672 bp. There were 108 genes in the plastome, including 76 protein-coding genes, 28 transfer RNA (tRNA), and four ribosomal RNA (rRNA) genes. Comparative analysis indicated that Codonopsis pilosula subsp. tangshen had an unusual large inversion in the large single-copy (LSC) region compared with the other three Codonopsis species. And there were two dispersed repeat sequences at both ends of the inverted regions, which might mediate the generation of this inversion. We found five hypervariable regions among the four Codonopsis species. PCR amplification and Sanger sequencing experiments demonstrated that two hypervariable regions could distinguish three medicinal Codonopsis species. Results obtained from this study will support taxonomic classification, discrimination, and molecular evolutionary studies of Codonopsis species.

High-resolution map of plastid-encoded RNA polymerase binding patterns demonstrates a major role of transcription in chloroplast gene expression

Plastids contain their own genomes, which are transcribed by two types of RNA polymerases. One of those enzymes is a bacterial-type, multi-subunit polymerase encoded by the plastid genome. The plastid-encoded RNA polymerase (PEP) is required for efficient expression of genes encoding proteins involved in photosynthesis. Despite the importance of PEP, its DNA binding locations have not been studied on the genome-wide scale at high resolution. We established a highly specific approach to detect the genome-wide pattern of PEP binding to chloroplast DNA using plastid chromatin immunoprecipitation-sequencing (ptChIP-seq). We found that in mature Arabidopsis thaliana chloroplasts, PEP has a complex DNA binding pattern with preferential association at genes encoding rRNA, tRNA, and a subset of photosynthetic proteins. Sigma factors SIG2 and SIG6 strongly impact PEP binding to a subset of tRNA genes and have more moderate effects on PEP binding throughout the rest of the genome. PEP binding is commonly enriched on gene promoters, around transcription start sites. Finally, the levels of PEP binding to DNA are correlated with levels of RNA accumulation, which demonstrates the impact of PEP on chloroplast gene expression. Presented data are available through a publicly available Plastid Genome Visualization Tool (Plavisto) at https://plavisto.mcdb.lsa.umich.edu/.

AtNusG, a chloroplast nucleoid protein of bacterial origin linking chloroplast transcriptional and translational machineries, is required for proper chloroplast gene expression in Arabidopsis thaliana

In Escherichia coli, transcription-translation coupling is mediated by NusG. Although chloroplasts are descendants of endosymbiotic prokaryotes, the mechanism underlying this coupling in chloroplasts remains unclear. Here, we report transcription-translation coupling through AtNusG in chloroplasts. AtNusG is localized in chloroplast nucleoids and is closely associated with the chloroplast PEP complex by interacting with its essential component PAP9. It also comigrates with chloroplast ribosomes and interacts with their two components PRPS5 (uS5c) and PRPS10 (uS10c). These data suggest that the transcription and translation machineries are coupled in chloroplasts. In the atnusg mutant, the accumulation of chloroplast-encoded photosynthetic gene transcripts, such as psbA, psbB, psbC and psbD, was not obviously changed, but that of their proteins was clearly decreased. Chloroplast polysomic analysis indicated that the decrease in these proteins was due to the reduced efficiency of their translation in this mutant, leading to reduced photosynthetic efficiency and enhanced sensitivity to cold stress. These data indicate that AtNusG-mediated coupling between transcription and translation in chloroplasts ensures the rapid establishment of photosynthetic capacity for plant growth and the response to environmental changes. Therefore, our study reveals a conserved mechanism of transcription-translation coupling between chloroplasts and E. coli, which perhaps represents a regulatory mechanism of chloroplast gene expression. This study provides insights into the underlying mechanisms of chloroplast gene expression in higher plants.

PAP8/pTAC6 Is Part of a Nuclear Protein Complex and Displays RNA Recognition Motifs of Viral Origin

Chloroplast biogenesis depends on a complex transcriptional program involving coordinated expression of plastid and nuclear genes. In particular, photosynthesis-associated plastid genes are expressed by the plastid-encoded polymerase (PEP) that undergoes a structural rearrangement during chloroplast formation. The prokaryotic-type core enzyme is rebuilt into a larger complex by the addition of nuclear-encoded PEP-associated proteins (PAP1 to PAP12). Among the PAPs, some have been detected in the nucleus (PAP5 and PAP8), where they could serve a nuclear function required for efficient chloroplast biogenesis. Here, we detected PAP8 in a large nuclear subcomplex that may include other subunits of the plastid-encoded RNA polymerase. We have made use of PAP8 recombinant proteins in Arabidopsis thaliana to decouple its nucleus- and chloroplast-associated functions and found hypomorphic mutants pointing at essential amino acids. While the origin of the PAP8 gene remained elusive, we have found in its sequence a micro-homologous domain located within a large structural homology with a rhinoviral RNA-dependent RNA polymerase, highlighting potential RNA recognition motifs in PAP8. PAP8 in vitro RNA binding activity suggests that this domain is functional. Hence, we propose that the acquisition of PAPs may have occurred during evolution by different routes, including lateral gene transfer.

The complete chloroplast genome sequencing analysis revealed an unusual IRs reduction in three species of subfamily Zygophylloideae

Tetraena mongolica, Zygophyllum xanthoxylon, and Z. fabago are three typical dryland plants with important ecological values in subfamily Zygophylloideae of Zygophyllaceae. Studies on the chloroplast genomes of them are favorable for understanding the diversity and phylogeny of Zygophyllaceae. Here, we sequenced and assembled the whole chloroplast genomes of T. mongolica, Z. xanthoxylon, and Z. fabago, and performed comparative genomic and phylogenetic analysis. The total size, structure, gene content and orders of these three chloroplast genomes were similar, and the three chloroplast genomes exhibited a typical quadripartite structure with a large single-copy region (LSC; 79,696–80,291 bp), a small single-copy region (SSC; 16,462–17,162 bp), and two inverted repeats (IRs; 4,288–4,413 bp). A total of 107 unique genes were identified from the three chloroplast genomes, including 70 protein-coding genes, 33 tRNAs, and 4 rRNAs. Compared with other angiosperms, the three chloroplast genomes were significantly reduced in overall length due to an unusual 16–24 kb shrinkage of IR regions and loss of the 11 genes which encoded subunits of NADH dehydrogenase. Genome-wide comparisons revealed similarities and variations between the three species and others. Phylogenetic analysis based on the three chloroplast genomes supported the opinion that Zygophyllaceae belonged to Zygophyllales in Fabids, and Z. xanthoxylon and Z. fabago belonged to Zygophyllum. The genome-wide comparisons revealed the similarity and variations between the chloroplast genomes of the three Zygophylloideae species and other plant species. This study provides a valuable molecular biology evidence for further studies of phylogenetic status of Zygophyllaceae.

Organellar Evolution: A Path from Benefit to Dependence

Eukaryotic organelles supposedly evolved from their bacterial ancestors because of their benefits to host cells. However, organelles are quite often retained, even when the beneficial metabolic pathway is lost, due to something other than the original beneficial function. The organellar function essential for cell survival is, in the end, the result of organellar evolution, particularly losses of redundant metabolic pathways present in both the host and endosymbiont, followed by a gradual distribution of metabolic functions between the organelle and host. Such biological division of metabolic labor leads to mutual dependence of the endosymbiont and host. Changing environmental conditions, such as the gradual shift of an organism from aerobic to anaerobic conditions or light to dark, can make the original benefit useless. Therefore, it can be challenging to deduce the original beneficial function, if there is any, underlying organellar acquisition. However, it is also possible that the organelle is retained because it simply resists being eliminated or digested untill it becomes indispensable.

Chloroplast pH Homeostasis for the Regulation of Photosynthesis

The pH of various chloroplast compartments, such as the thylakoid lumen and stroma, is light-dependent. Light illumination induces electron transfer in the photosynthetic apparatus, coupled with proton translocation across the thylakoid membranes, resulting in acidification and alkalization of the thylakoid lumen and stroma, respectively. Luminal acidification is crucial for inducing regulatory mechanisms that protect photosystems against photodamage caused by the overproduction of reactive oxygen species (ROS). Stromal alkalization activates enzymes involved in the Calvin-Benson-Bassham (CBB) cycle. Moreover, proton translocation across the thylakoid membranes generates a proton gradient (ΔpH) and an electric potential (ΔΨ), both of which comprise the proton motive force (pmf) that drives ATP synthase. Then, the synthesized ATP is consumed in the CBB cycle and other chloroplast metabolic pathways. In the dark, the pH of both the chloroplast stroma and thylakoid lumen becomes neutral. Despite extensive studies of the above-mentioned processes, the molecular mechanisms of how chloroplast pH can be maintained at proper levels during the light phase for efficient activation of photosynthesis and other metabolic pathways and return to neutral levels during the dark phase remain largely unclear, especially in terms of the precise control of stromal pH. The transient increase and decrease in chloroplast pH upon dark-to-light and light-to-dark transitions have been considered as signals for controlling other biological processes in plant cells. Forward and reverse genetic screening approaches recently identified new plastid proteins involved in controlling ΔpH and ΔΨ across the thylakoid membranes and chloroplast proton/ion homeostasis. These proteins have been conserved during the evolution of oxygenic phototrophs and include putative photosynthetic protein complexes, proton transporters, and/or their regulators. Herein, we summarize the recently identified protein players that control chloroplast pH and influence photosynthetic efficiency in plants.

Complete Chloroplast Genome Analysis of Two Important Medicinal Alpinia Species: Alpinia galanga and Alpinia kwangsiensis

Most Alpinia species are valued as foods, ornamental plants, or plants with medicinal properties. However, morphological characteristics and commonly used DNA barcode fragments are not sufficient for accurately identifying Alpinia species. Difficulties in species identification have led to confusion in the sale and use of Alpinia for medicinal use. To mine resources and improve the molecular methods for distinguishing among Alpinia species, we report the complete chloroplast (CP) genomes of Alpinia galanga and Alpinia kwangsiensis species, obtained via high-throughput Illumina sequencing. The CP genomes of A. galanga and A. kwangsiensis exhibited a typical circular tetramerous structure, including a large single-copy region (87,565 and 87,732 bp, respectively), a small single-copy region (17,909 and 15,181 bp, respectively), and a pair of inverted repeats (27,313 and 29,705 bp, respectively). The guanine-cytosine content of the CP genomes is 36.26 and 36.15%, respectively. Furthermore, each CP genome contained 133 genes, including 87 protein-coding genes, 38 distinct tRNA genes, and 8 distinct rRNA genes. We identified 110 and 125 simple sequence repeats in the CP genomes of A. galanga and A. kwangsiensis, respectively. We then combined these data with publicly available CP genome data from four other Alpinia species (A. hainanensis, A. oxyphylla, A. pumila, and A. zerumbet) and analyzed their sequence characteristics. Nucleotide diversity was analyzed based on the alignment of the complete CP genome sequences, and five candidate highly variable site markers (trnS-trnG, trnC-petN, rpl32-trnL, psaC-ndhE, and ndhC-trnV) were found. Twenty-eight complete CP genome sequences belonging to Alpinieae species were used to construct phylogenetic trees. The results fully demonstrated the phylogenetic relationship among the genera of the Alpinieae, and further proved that Alpinia is a non-monophyletic group. The complete CP genomes of the two medicinal Alpinia species provides lays the foundation for the use of CP genomes in species identification and phylogenetic analyses of Alpinia species.

A Novel Chloroplast Protein RNA Processing 8 Is Required for the Expression of Chloroplast Genes and Chloroplast Development in Arabidopsis thaliana

Chloroplast development involves the coordinated expression of both plastids- and nuclear-encoded genes in higher plants. However, the underlying mechanism still remains largely unknown. In this study, we isolated and characterized an Arabidopsis mutant with an albino lethality phenotype named RNA processing 8 (rp8). Genetic complementation analysis demonstrated that the gene AT4G37920 (RP8) was responsible for the mutated phenotype. The RP8 gene was strongly expressed in photosynthetic tissues at both transcription and translation protein levels. The RP8 protein is localized in the chloroplast and associated with the thylakoid. Disruption of the RP8 gene led to a defect in the accumulation of the rpoA mature transcript, which reduced the level of the RpoA protein, and affected the transcription of PEP-dependent genes. The abundance of the chloroplast rRNA, including 23S, 16S, 4.5S, and 5S rRNA, were reduced in the rp8 mutant, respectively, and the amounts of chloroplast ribosome proteins, such as, PRPS1(uS1c), PRPS5(uS5c), PRPL2 (uL2c), and PRPL4 (uL4c), were substantially decreased in the rp8 mutant, which indicated that knockout of RP8 seriously affected chloroplast translational machinery. Accordingly, the accumulation of photosynthetic proteins was seriously reduced. Taken together, these results indicate that the RP8 protein plays an important regulatory role in the rpoA transcript processing, which is required for the expression of chloroplast genes and chloroplast development in Arabidopsis.

Comparative Analysis of Complete Chloroplast Genomes of 13 Species in Epilobium, Circaea, and Chamaenerion and Insights Into Phylogenetic Relationships of Onagraceae

The evening primrose family, Onagraceae, is a well defined family of the order Myrtales, comprising 22 genera widely distributed from boreal to tropical areas. In this study, we report and characterize the complete chloroplast genome sequences of 13 species in Circaea, Chamaenerion, and Epilobium using a next-generation sequencing method. We also retrieved chloroplast sequences from two other Onagraceae genera to characterize the chloroplast genome of the family. The complete chloroplast genomes of Onagraceae encoded an identical set of 112 genes (with exclusion of duplication), including 78 protein-coding genes, 30 transfer RNAs, and four ribosomal RNAs. The chloroplast genomes are basically conserved in gene arrangement across the family. However, a large segment of inversion was detected in the large single copy region of all the samples of Oenothera subsect. Oenothera. Two kinds of inverted repeat (IR) region expansion were found in Oenothera, Chamaenerion, and Epilobium samples. We also compared chloroplast genomes across the Onagraceae samples in some features, including nucleotide content, codon usage, RNA editing sites, and simple sequence repeats (SSRs). Phylogeny was inferred by the chloroplast genome data using maximum-likelihood (ML) and Bayesian inference methods. The generic relationship of Onagraceae was well resolved by the complete chloroplast genome sequences, showing potential value in inferring phylogeny within the family. Phylogenetic relationship in Oenothera was better resolved than other densely sampled genera, such as Circaea and Epilobium. Chloroplast genomes of Oenothera subsect. Oenothera, which are biparental inheritated, share a syndrome of characteristics that deviate from primitive pattern of the family, including slightly expanded inverted repeat region, intron loss in clpP, and presence of the inversion.

A study of RNA-editing in Populus trichocarpa nuclei revealed acquisition of RNA-editing on the endosymbiont-derived genes, and a preference for intracellular remodeling genes in adaptation to endosymbiosis

RNA-editing is a post-transcriptional modification that can diversify genome-encoded information by modifying individual RNA bases. In contrast to the well-studied RNA-editing in organelles, little is known about nuclear RNA-editing in higher plants. We performed a genome-wide study of RNA-editing in Populus trichocarpa nuclei using the RNA-seq data generated from the sequenced poplar genotype, 'Nisqually-1'. A total of 24,653 nuclear RNA-editing sites present in 8,603 transcripts were identified. Notably, RNA-editing in P. trichocarpa nuclei tended to occur on endosymbiont-derived genes. We then scrutinized RNA-editing in a cyanobacterial strain closely related to chloroplast. No RNA-editing sites were identified therein, implying that RNA-editing of these endosymbiont-derived genes was acquired after endosymbiosis. Gene ontology enrichment analysis of all the edited genes in P. trichocarpa nuclei demonstrated that nuclear RNA-editing was primarily focused on genes involved in intracellular remodeling processes, which suggests that RNA-editing plays contributing roles in organellar establishment during endosymbiosis. We built a coexpression network using all C-to-U edited genes and then decomposed it to obtain 18 clusters, six of which contained a conserved core motif, A/G-C-A/G. Such a short core motif not only attracted the RNA-editing machinery but also enabled large numbers of sites to be targeted though further study is necessary to verify this finding.

The genetic and epigenetic landscape of the Arabidopsis centromeres

Centromeres attach chromosomes to spindle microtubules during cell division and, despite this conserved role, show paradoxically rapid evolution and are typified by complex repeats. We used long-read sequencing to generate the Col-CEN Arabidopsis thaliana genome assembly that resolves all five centromeres. The centromeres consist of megabase-scale tandemly repeated satellite arrays, which support CENTROMERE SPECIFIC HISTONE H3 (CENH3) occupancy and are densely DNA methylated, with satellite variants private to each chromosome. CENH3 preferentially occupies satellites that show the least amount of divergence and occur in higher-order repeats. The centromeres are invaded by ATHILA retrotransposons, which disrupt genetic and epigenetic organization. Centromeric crossover recombination is suppressed, yet low levels of meiotic DNA double-strand breaks occur that are regulated by DNA methylation. We propose that Arabidopsis centromeres are evolving through cycles of satellite homogenization and retrotransposon-driven diversification.

The Arabidopsis PeptideAtlas: Harnessing worldwide proteomics data to create a comprehensive community proteomics resource

We developed a resource, the Arabidopsis PeptideAtlas (www.peptideatlas.org/builds/arabidopsis/), to solve central questions about the Arabidopsis thaliana proteome, such as the significance of protein splice forms and post-translational modifications (PTMs), or simply to obtain reliable information about specific proteins. PeptideAtlas is based on published mass spectrometry (MS) data collected through ProteomeXchange and reanalyzed through a uniform processing and metadata annotation pipeline. All matched MS-derived peptide data are linked to spectral, technical, and biological metadata. Nearly 40 million out of ∼143 million MS/MS (tandem MS) spectra were matched to the reference genome Araport11, identifying ∼0.5 million unique peptides and 17,858 uniquely identified proteins (only isoform per gene) at the highest confidence level (false discovery rate 0.0004; 2 non-nested peptides ≥9 amino acid each), assigned canonical proteins, and 3,543 lower-confidence proteins. Physicochemical protein properties were evaluated for targeted identification of unobserved proteins. Additional proteins and isoforms currently not in Araport11 were identified that were generated from pseudogenes, alternative start, stops, and/or splice variants, and small Open Reading Frames; these features should be considered when updating the Arabidopsis genome. Phosphorylation can be inspected through a sophisticated PTM viewer. PeptideAtlas is integrated with community resources including TAIR, tracks in JBrowse, PPDB, and UniProtKB. Subsequent PeptideAtlas builds will incorporate millions more MS/MS data.

The pentatricopeptide repeat protein EMB1270 interacts with CFM2 to splice specific group II introns in Arabidopsis chloroplasts

Chloroplast biogenesis requires the coordinated expression of chloroplast and nuclear genes. Here, we show that EMB1270, a plastid-localized pentatricopeptide repeat (PPR) protein, is required for chloroplast biogenesis in Arabidopsis thaliana. Knockout of EMB1270 led to embryo arrest, whereas a mild knockdown mutant of EMB1270 displayed a virescent phenotype. Almost no photosynthetic proteins accumulated in the albino emb1270 knockout mutant. By contrast, in the emb1270 knockdown mutant, the levels of ClpP1 and photosystem I (PSI) subunits were significantly reduced, whereas the levels of photosystem II (PSII) subunits were normal. Furthermore, the splicing efficiencies of the clpP1.2, ycf3.1, ndhA, and ndhB plastid introns were dramatically reduced in both emb1270 mutants. RNA immunoprecipitation revealed that EMB1270 associated with these introns in vivo. In an RNA electrophoretic mobility shift assay (REMSA), a truncated EMB1270 protein containing the 11 N-terminal PPR motifs bound to the predicted sequences of the clpP1.2, ycf3.1, and ndhA introns. In addition, EMB1270 specifically interacted with CRM Family Member 2 (CFM2). Given that CFM2 is known to be required for splicing the same plastid RNAs, our results suggest that EMB1270 associates with CFM2 to facilitate the splicing of specific group II introns in Arabidopsis.