The plastid clpP1 protease gene is essential for plant development

Maize requires Embryo defective27 for embryogenesis and seedling development

The essential role of plastid translation in embryogenesis has been established in many plants, but a retrograde signal triggered by defective plastid translation machinery that may leads to embryogenesis arrest remains unknown. In this study, we characterized an embryo defective27 (emb27) mutant in maize (Zea mays), and cloning indicates that Emb27 encodes the plastid ribosomal protein S13. The null mutant emb27-1 conditions an emb phenotype with arrested embryogenesis; however, the leaky mutant emb27-2 exhibits normal embryogenesis but an albino seedling-lethal phenotype. The emb27-1/emb27-2 trans-heterozygotes display varying phenotypes from emb to normal seeds but albino seedlings. Analysis of the Emb27 transcription levels in these mutants revealed that the Emb27 expression level in the embryo corresponds with the phenotypic expression of the emb27 mutants. In the W22 genetic background, an Emb27 transcription level higher than 6% of the wild-type level renders normal embryogenesis, whereas lower than that arrests embryogenesis. Mutation of Emb27 reduces the level of plastid 16S rRNA and the accumulation of the plastid-encoded proteins. As a secondary effect, splicing of several plastid introns was impaired in emb27-1 and 2 other plastid translation-defective mutants, emb15 and emb16, suggesting that plastome-encoded factors are required for the splicing of these introns, such as Maturase K (MatK). Our results indicate that EMB27 is essential for plastid protein translation, embryogenesis, and seedling development in maize and reveal an expression threshold of Emb27 for maize embryogenesis.

Silencing an ATP-Dependent Caseinolytic Protease Proteolytic Subunit Gene Enhances the Resistance of Rice to Nilaparvata lugens

The ATP-dependent caseinolytic protease (Clp) system has been reported to play an important role in plant growth, development, and defense against pathogens. However, whether the Clp system is involved in plant defense against herbivores remains largely unclear. We explore the role of the Clp system in rice defenses against brown planthopper (BPH) Nilaparvata lugens by combining chemical analysis, transcriptome, and molecular analyses, as well as insect bioassays. We found the expression of a rice Clp proteolytic subunit gene, OsClpP6, was suppressed by infestation of BPH gravid females and mechanical wounding. Silencing OsClpP6 enhanced the level of BPH-induced jasmonic acid (JA), JA-isoleucine (JA-Ile), and ABA, which in turn promoted the production of BPH-elicited rice volatiles and increased the resistance of rice to BPH. Field trials showed that silencing OsClpP6 decreased the population densities of BPH and WBPH. We also observed that silencing OsClpP6 decreased chlorophyll content in rice leaves at early developmental stages and impaired rice root growth and seed setting rate. These findings demonstrate that an OsClpP6-mediated Clp system in rice was involved in plant growth-defense trade-offs by affecting the biosynthesis of defense-related signaling molecules in chloroplasts. Moreover, rice plants, after recognizing BPH infestation, can enhance rice resistance to BPH by decreasing the Clp system activity. The work might provide a new way to breed rice varieties that are resistant to herbivores.

Characterization of Firmiana danxiaensis plastomes and comparative analysis of Firmiana: insight into its phylogeny and evolution

BACKGROUND

Firmiana danxiaensis is a critically endangered and ecologically important tree currently only found in four locations in Danxia or Karst habitats in northern Guangdong Province, China. The specialized habitat preference makes it an ideal model species for study of adaptive evolution. Meanwhile, the phylogenetic relationships of F. danxiaensis in four locations under two landforms are unclear. Therefore, we sequenced its complete chloroplast (cp.) genomes and conducted comprehensive interspecific and intrageneric plastome studies.

RESULTS

The F. danxiaensis plastomes in four locations showed a typical quadripartite and circular structure that ranged from 160,832 to 161,206 bp in size, with 112 unique genes encoded. Comparative genomics showed that the plastomes of F. danxiaensis were relatively conserved with high similarity of genome organization, gene number, GC content and SSRs. While the genomes revealed higher biased codon preferences in Karst habitat than those in Danxia habitats. Eighteen and 11 divergent hotpots were identified at interspecific and intrageneric levels for species identification and further phylogenetic studies. Seven genes (clpP, accD, ccsA, ndhH, rpl20, rpoC2, and rps4) were under positive selection and may be related to adaptation. Phylogenetic analysis revealed that F. danxiaensis is sister to F. major and F. simplex. However, the interspecific relationships are not consistent with the habitat types.

CONCLUSIONS

The characteristics and interspecific relationship of F. danxiaensis plastomes provide new insights into further integration of geographical factors, environmental factors, and genetic variations on the genomic study of F. danxiaensis. Together, our study will contribute to the study of species identification, population genetics, and conservation biology of F. danxiaensis.

Genome-wide identification, characterization, and expression analysis of CCT transcription factors in poplar

The CCT [CONSTANS (CO), CO-like, and TIMING OF CAB EXPRESSION1 (TOC1)] gene family is involved in photoperiodic flowering and adaptation to different environments. In this study, 39 CCT family genes from the poplar genome were identified and characterized, including 18 COL, 7 PRR, and 14 CMF TFs. Phylogenetics analysis showed that the PtrCCT gene family could be classified into five classes (Classes I-V) that have close relationships with Arabidopsis thaliana. Eight pairs of PtrCCTs had collinear relationships through interchromosomal synteny analysis in poplar, suggesting segmental duplication played a vital role in the expansion of the poplar CCT gene family. Besides, synteny analyses of the CCT members among poplar and different species provided more clues for PtrCCT gene family evolution. Cis-acting elements in the promoters of PtrCCTs predicted their involvement in light responses, hormone responses, biotic/abiotic stress responses, and plant growth and development. Eight members of the PpnCCT gene family were differentially expressed in the apical buds and leaves of triploid poplar compared to diploids. We then focused on PpnCCT39 upregulated in triploid poplars and showed that PpnCCT39 was localized in the nucleus, chloroplast, and cytoplasm and could interact with CLPP1 in the chloroplast. Overexpression of PpnCCT39 in poplar increased chlorophyll contents and enhanced photosynthetic rate. This study provided comprehensive information for the CCT gene family and set up a basis for its function identification in poplar.

Comparative chloroplast genome analysis of seven extant Citrullus species insight into genetic variation, phylogenetic relationships, and selective pressure

Citrullus ecirrhosus, Citrullus rehmii, and Citrullus naudinianus are three important related wild species of watermelon in the genus Citrullus, and their morphological differences are clear, however, their chloroplast genome differences remain unknown. This study is the first to assemble, analyze, and publish the complete chloroplast genomes of C. ecirrhosus, C. rehmii, and C. naudinianus. A comparative analysis was then conducted among the complete chloroplast genomes of seven extant Citrullus species, and the results demonstrated that the average genome sizes of Citrullus is 157,005 bp, a total of 130-133 annotated genes were identified, including 8 rRNA, 37 tRNA and 85-88 protein-encoding genes. Their gene content, order, and genome structure were similar. However, noncoding regions were more divergent than coding regions, and rps16-trnQ was a hypervariable fragment. Thirty-four polymorphic SSRs, 1,271 SNPs and 234 INDELs were identified. Phylogenetic trees revealed a clear phylogenetic relationship of Citrullus species, and the developed molecular markers (SNPs and rps16-trnQ) could be used for taxonomy in Citrullus. Three genes (atpB, clpP1, and rpoC2) were identified to undergo selection and would promote the environmental adaptation of Citrullus.

Chloroplast proteostasis: A story of birth, life, and death

Protein homeostasis (proteostasis) is a dynamic balance of protein synthesis and degradation. Because of the endosymbiotic origin of chloroplasts and the massive transfer of their genetic information to the nucleus of the host cell, many protein complexes in the chloroplasts are constituted from subunits encoded by both genomes. Hence, the proper function of chloroplasts relies on the coordinated expression of chloroplast- and nucleus-encoded genes. The biogenesis and maintenance of chloroplast proteostasis are dependent on synthesis of chloroplast-encoded proteins, import of nucleus-encoded chloroplast proteins from the cytosol, and clearance of damaged or otherwise undesired "old" proteins. This review focuses on the regulation of chloroplast proteostasis, its interaction with proteostasis of the cytosol, and its retrograde control over nuclear gene expression. We also discuss significant issues and perspectives for future studies and potential applications for improving the photosynthetic performance and stress tolerance of crops.

Revealing the Complete Chloroplast Genome of an Andean Horticultural Crop, Sweet Cucumber (Solanum muricatum), and Its Comparison with Other Solanaceae Species

Sweet cucumber (Solanum muricatum) sect. Basarthrum is a neglected horticultural crop native to the Andean region. It is naturally distributed very close to other two Solanum crops of high importance, potatoes, and tomatoes. To date, molecular tools for this crop remain undetermined. In this study, the complete sweet cucumber chloroplast (cp) genome was obtained and compared with seven Solanaceae species. The cp genome of S. muricatum was 155,681 bp in length and included a large single copy (LSC) region of 86,182 bp and a small single-copy (SSC) region of 18,360 bp, separated by a pair of inverted repeats (IR) regions of 25,568 bp. The cp genome possessed 87 protein-coding genes (CDS), 37 transfer RNA (tRNA) genes, eight ribosomal RNA (rRNA) genes, and one pseudogene. Furthermore, 48 perfect microsatellites were identified. These repeats were mainly located in the noncoding regions. Whole cp genome comparative analysis revealed that the SSC and LSC regions showed more divergence than IR regions. Similar to previous studies, our phylogenetic analysis showed that S. muricatum is a sister species to members of sections Petota + Lycopersicum + Etuberosum. We expect that this first sweet cucumber chloroplast genome will provide potential molecular markers and genomic resources to shed light on the genetic diversity and population studies of S. muricatum, which will allow us to identify varieties and ecotypes. Finally, the features and the structural differentiation will provide us with information about the genes of interest, generating tools for the most precise selection of the best individuals of sweet cucumber, in less time and with fewer resources.

The plastid-encoded protein Orf2971 is required for protein translocation and chloroplast quality control

Photosynthesis and the biosynthesis of many important metabolites occur in chloroplasts. In these semi-autonomous organelles, the chloroplast genome encodes approximately 100 proteins. The remaining chloroplast proteins, close to 3,000, are encoded by nuclear genes whose products are translated in the cytosol and imported into chloroplasts. However, there is still no consensus on the composition of the protein import machinery including its motor proteins and on how newly imported chloroplast proteins are refolded. In this study, we have examined the function of orf2971, the largest chloroplast gene of Chlamydomonas reinhardtii. The depletion of Orf2971 causes the accumulation of protein precursors, partial proteolysis and aggregation of proteins, increased expression of chaperones and proteases, and autophagy. Orf2971 interacts with the TIC (translocon at the inner chloroplast envelope) complex, catalyzes ATP (adenosine triphosphate) hydrolysis, and associates with chaperones and chaperonins. We propose that Orf2971 is intimately connected to the protein import machinery and plays an important role in chloroplast protein quality control.

Complete chloroplast genome of Lilium ledebourii (Baker) Boiss and its comparative analysis: lights into selective pressure and adaptive evolution

Lilium ledebourii (Baker) Boiss is a rare species, which exhibits valuable traits. However, before its genetic diversity and evolutionary were uncovered, its wild resources were jeopardized. Moreover, some ambiguities in phylogenetic relationships of this genus remain unresolved. Therefore, obtaining the whole chloroplast sequences of L. ledebourii and its comparative analysis along with other Lilium species is crucial and pivotal to understanding the evolution of this genus as well as the genetic populations. A multi-scale genome-level analysis, especially selection pressure, was conducted. Detailed third‑generation sequencing and analysis revealed a whole chloroplast genome of 151,884 bp, with an ordinary quadripartite and protected structure comprising 37.0% GC. Overall, 113 different genes were recognized in the chloroplast genome, consisting of 30 distinct tRNA genes, four distinct ribosomal RNAs genes, and 79 unique protein-encoding genes. Here, 3234 SSRs and 2053 complex repeats were identified, and a comprehensive analysis was performed for IR expansion and contraction, and codon usage bias. Moreover, genome-wide sliding window analysis revealed the variability of rpl32-trnL-ccsA, petD-rpoA, ycf1, psbI-trnS-trnG, rps15-ycf1, trnR, trnT-trnL, and trnP-psaJ-rpl33 were higher among the 48 Lilium cp genomes, displaying higher variability of nucleotide in SC regions. Following 1128 pairwise comparisons, ndhB, psbJ, psbZ, and ycf2 exhibit zero synonymous substitution, revealing divergence or genetic restriction. Furthermore, out of 78 protein-coding genes, we found that accD and rpl36 under positive selection: however, at the entire-chloroplast protein scale, the Lilium species have gone through a purifying selection. Also, a new phylogenetic tree for Lilium was rebuilt, and we believe that the Lilium classification is clearer than before. The genetic resources provided here will aid future studies in species identification, population genetics, and Lilium conservation.

Comparative Plastome Analysis of Three Amaryllidaceae Subfamilies: Insights into Variation of Genome Characteristics, Phylogeny, and Adaptive Evolution

In the latest APG IV classification system, Amaryllidaceae is placed under the order of Asparagus and includes three subfamilies: Agapanthoideae, Allioideae, and Amaryllidoideae, which include many economically important crops. With the development of molecular phylogeny, research on the phylogenetic relationship of Amaryllidaceae has become more convenient. However, the current comparative analysis of Amaryllidaceae at the whole chloroplast genome level is still lacking. In this study, we sequenced 18 Allioideae plastomes and combined them with publicly available data (a total of 41 plastomes), including 21 Allioideae species, 1 Agapanthoideae species, 14 Amaryllidoideae species, and 5 Asparagaceae species. Comparative analyses were performed including basic characteristics of genome structure, codon usage, repeat elements, IR boundary, and genome divergence. Phylogenetic relationships were detected using single-copy genes (SCGs) and ribosomal internal transcribed spacer sequences (ITS), and the branch-site model was also employed to conduct the positive selection analysis. The results indicated that all Amaryllidaceae species showed a highly conserved typical tetrad structure. The GC content and five codon usage indexes in Allioideae species were lower than those in the other two subfamilies. Comparison analysis of Bayesian and ML phylogeny based on SCGs strongly supports the monophyly of three subfamilies and the sisterhood among them. Besides, positively selected genes (PSGs) were detected in each of the three subfamilies. Almost all genes with significant posterior probabilities for codon sites were associated with self-replication and photosynthesis. Our study investigated the three subfamilies of Amaryllidaceae at the whole chloroplast genome level and suggested the key role of selective pressure in the adaptation and evolution of Amaryllidaceae.

Favorable QTLs from Oryza longistaminata improve rice drought resistance

BACKGROUND

Drought is the major abiotic stress to rice grain production under unpredictable changing climatic environments. Wild rice of O. longistaminata show diverse responses and strong tolerance to stress environments. In order to identify whether the O. longistaminata can improve the rice drought resistance or not, a BIL population of 143 BC₂F₂₀ lines derived from the cross between the cultivar rice 9311 and O. longistaminata were assessed under stress of 20% PEG6000.

RESULTS

In total, 28 QTLs related to drought resistance based on eight agronomic traits of seedlings were identified. Of which, thirteen QTLs including two QTLs for leaf drying, one QTL for leaf rolling, one QTL for leaf number, five QTLs for dry weight of root, two QTLs for dry weight of shoot, one QTL for maximum root length and two QTLs for maximum shoot length were derived from O. longistaminata. What's more, qDWR8.1 for dry weight of root was repeatedly detected and fine-mapped to an interval about 36.2 Kb. The unique allele of MH08g0242800 annotated as ATP-dependent Clp protease proteolytic subunit from O. longistaminata was suggested as the candidate gene for drought resistance. Further, six representative BIL lines were stably characterized showing significantly stronger drought resistance than 9311 based on principle component analysis, they each contained 2 ~ 5 QTLs including qDWR8.1 from O. longistaminata.

CONCLUSIONS

Together, our results indicate that the QTLs from O. longistaminata can effectively enhance the drought tolerance of rice, showing great potential value in breeding of elite rice varieties, which will lay a novel insight into the genetic network for drought tolerance of rice.

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

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

Rapid sequence evolution is associated with genetic incompatibilities in the plastid Clp complex

KEY MESSAGE

Replacing the native clpP1 gene in the Nicotiana plastid genome with homologs from different donor species showed that the extent of genetic incompatibilities depended on the rate of sequence evolution. The plastid caseinolytic protease (Clp) complex plays essential roles in maintaining protein homeostasis and comprises both plastid-encoded and nuclear-encoded subunits. Despite the Clp complex being retained across green plants with highly conserved protein sequences in most species, examples of extremely accelerated amino acid substitution rates have been identified in numerous angiosperms. The causes of these accelerations have been the subject of extensive speculation but still remain unclear. To distinguish among prevailing hypotheses and begin to understand the functional consequences of rapid sequence divergence in Clp subunits, we used plastome transformation to replace the native clpP1 gene in tobacco (Nicotiana tabacum) with counterparts from another angiosperm genus (Silene) that exhibits a wide range in rates of Clp protein sequence evolution. We found that antibiotic-mediated selection could drive a transgenic clpP1 replacement from a slowly evolving donor species (S. latifolia) to homoplasmy but that clpP1 copies from Silene species with accelerated evolutionary rates remained heteroplasmic, meaning that they could not functionally replace the essential tobacco clpP1 gene. These results suggest that observed cases of rapid Clp sequence evolution are a source of epistatic incompatibilities that must be ameliorated by coevolutionary responses between plastid and nuclear subunits.

The Intron Retention Variant CsClpP3m Is Involved in Leaf Chlorosis in Some Tea Cultivars

Tea products made from chlorotic or albino leaves are very popular for their unique flavor. Probing into the molecular mechanisms underlying the chlorotic leaf phenotype is required to better understand the formation of these tea cultivars and aid in future practical breeding. In this study, transcriptional alterations of multiple subunit genes of the caseinolytic protease complex (Clp) in the chlorotic tea cultivar 'Yu-Jin-Xiang' (YJX) were found. Cultivar YJX possessed the intron retention variant of ClpP3, named as CsClpP3m, in addition to the non-mutated ClpP3. The mutated variant results in a truncated protein containing only 166 amino acid residues and lacks the catalytic triad S182-H206-D255. Quantitative analysis of two CsClpP3 variants in different leaves with varying degrees of chlorosis in YJX and analyses of different chlorotic tea cultivars revealed that the transcript ratios of CsClpP3m over CsClpP3 were negatively correlated with leaf chlorophyll contents. The chlorotic young leaf phenotype was also generated in the transgenic tobacco by suppressing ClpP3 using the RNAi method; complementation with non-mutated CsClpP3 rescued the wild-type phenotype, whereas CsClpP3m failed to complement. Taken together, CsClpP3m is involved in leaf chlorosis in YJX and some other tea cultivars in a dose-dependent manner, likely resulting from the failure of Clp complex assembly due to the truncated sequence of CsClpP3m. Our data shed light on the mechanisms controlling leaf chlorosis in tea plants.

Unlocking the Complete Chloroplast Genome of a Native Tree Species from the Amazon Basin, Capirona (Calycophyllum Spruceanum, Rubiaceae), and Its Comparative Analysis with Other Ixoroideae Species

Capirona (Calycophyllum spruceanum Benth.) belongs to subfamily Ixoroideae, one of the major lineages in the Rubiaceae family, and is an important timber tree. It originated in the Amazon Basin and has widespread distribution in Bolivia, Peru, Colombia, and Brazil. In this study, we obtained the first complete chloroplast (cp) genome of capirona from the department of Madre de Dios located in the Peruvian Amazon. High-quality genomic DNA was used to construct libraries. Pair-end clean reads were obtained by PE 150 library and the Illumina HiSeq 2500 platform. The complete cp genome of C. spruceanum has a 154,480 bp in length with typical quadripartite structure, containing a large single copy (LSC) region (84,813 bp) and a small single-copy (SSC) region (18,101 bp), separated by two inverted repeat (IR) regions (25,783 bp). The annotation of C. spruceanum cp genome predicted 87 protein-coding genes (CDS), 8 ribosomal RNA (rRNA) genes, 37 transfer RNA (tRNA) genes, and one pseudogene. A total of 41 simple sequence repeats (SSR) of this cp genome were divided into mononucleotides (29), dinucleotides (5), trinucleotides (3), and tetranucleotides (4). Most of these repeats were distributed in the noncoding regions. Whole chloroplast genome comparison with the other six Ixoroideae species revealed that the small single copy and large single copy regions showed more divergence than inverted regions. Finally, phylogenetic analyses resolved that C. spruceanum is a sister species to Emmenopterys henryi and confirms its position within the subfamily Ixoroideae. This study reports for the first time the genome organization, gene content, and structural features of the chloroplast genome of C. spruceanum, providing valuable information for genetic and evolutionary studies in the genus Calycophyllum and beyond.

Plant protease as regulator and signaling molecule for enhancing environmental stress-tolerance

Proteases are ubiquitous in prokaryotes and eukaryotes. Plant proteases are key regulators of various physiological processes, including protein homeostasis, organelle development, senescence, seed germination, protein processing, environmental stress response, and programmed cell death. Proteases are involved in the breakdown of peptide bonds resulting in irreversible posttranslational modification of the protein. Proteases act as signaling molecules that specifically regulate cellular function by cleaving and triggering receptor molecules. Peptides derived from proteolysis regulate ROS signaling under oxidative stress in the plant. It degrades misfolded and abnormal proteins into amino acids to repair the cell damage and regulates the biological process in response to environmental stress. Proteases modulate the biogenesis of phytohormones which control plant growth, development, and environmental stresses. Protein homeostasis, the overall balance between protein synthesis and proteolysis, is required for plant growth and development. Abiotic and biotic stresses are major factors that negatively impact cellular survivability, biomass production, and reduced crop yield potentials. Therefore, the identification of various stress-responsive proteases and their molecular functions may elucidate valuable information for the development of stress-resilient crops with higher yield potentials. However, the understanding of molecular mechanisms of plant protease remains unexplored. This review provides an overview of proteases related to development, signaling, and growth regulation to acclimatize environmental stress in plants.

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.

The chloroplast genomes of four Bupleurum (Apiaceae) species endemic to Southwestern China, a diversity center of the genus, as well as their evolutionary implications and phylogenetic inferences

Background

As one of the largest genera in Apiaceae, Bupleurum L. is well known for its high medicinal value. The genus has frequently attracted the attention of evolutionary biologist and taxonomist for its distinctive characteristics in the Apiaceae family. Although some chloroplast genomes data have been now available, the changes in the structure of chloroplast genomes and selective pressure in the genus have not been fully understood. In addition, few of the species are endemic to Southwest China, a distribution and diversity center of Chinese Bupleurum. Endemic species are key components of biodiversity and ecosystems, and investigation of the chloroplast genomes features of endemic species in Bupleurum will be helpful to develop a better understanding of evolutionary process and phylogeny of the genus. In this study, we analyzed the sequences of whole chloroplast genomes of 4 Southwest China endemic Bupleurum species in comparison with the published data of 17 Bupleurum species to determine the evolutionary characteristics of the genus and the phylogenetic relationships of Asian Bupleurum.

Results

The complete chloroplast genome sequences of the 4 endemic Bupleurum species are 155,025 bp to 155,323 bp in length including a SSC and a LSC region separated by a pair of IRs. Comparative analysis revealed an identical chloroplast gene content across the 21 Bupleurum species, including a total of 114 unique genes (30 tRNA genes, 4 rRNA genes and 80 protein-coding genes). Chloroplast genomes of the 21 Bupleurum species showed no rearrangements and a high sequence identity (96.4–99.2%). They also shared a similar tendency of SDRs and SSRs, but differed in number (59–83). In spite of their high conservation, they contained some mutational hotspots, which can be potentially exploited as high-resolution DNA barcodes for species discrimination. Selective pressure analysis showed that four genes were under positive selection. Phylogenetic analysis revealed that the 21 Bupleurum formed two major clades, which are likely to correspond to their geographical distribution.

Conclusions

The chloroplast genome data of the four endemic Bupleurum species provide important insights into the characteristics and evolution of chloroplast genomes of this genu, and the phylogeny of Bupleurum.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08008-z.

Comparative Chloroplast Genomes of Zosteraceae Species Provide Adaptive Evolution Insights Into Seagrass

Seagrasses are marine flowering plants found in tropical and sub-tropical areas that live in coastal regions between the sea and land. All seagrass species evolved from terrestrial monocotyledons, providing the opportunity to study plant adaptation to sea environments. Here, we sequenced the chloroplast genomes (cpGenomes) of three Zostera species, then analyzed and compared their cpGenome structures and sequence variations. We also performed a phylogenetic analysis using published seagrass chloroplasts and calculated the selection pressure of 17 species within seagrasses and nine terrestrial monocotyledons, as well as estimated the number of shared genes of eight seagrasses. The cpGenomes of Zosteraceae species ranged in size from 143,877 bp (Zostera marina) to 152,726 bp (Phyllospadix iwatensis), which were conserved and displayed similar structures and gene orders. Additionally, we found 17 variable hotspot regions as candidate DNA barcodes for Zosteraceae species, which will be helpful for studying the phylogenetic relationships and interspecies differences between seagrass species. Interestingly, nine genes had positive selection sites, including two ATP subunit genes (atpA and atpF), two ribosome subunit genes (rps4 and rpl20), two DNA-dependent RNA polymerase genes (rpoC1 and rpoC2), as well as accD, clpP, and ycf2. These gene regions may have played key roles in the seagrass adaptation to diverse environments. The Branch model analysis showed that seagrasses had a higher rate of evolution than terrestrial monocotyledons, suggesting that seagrasses experienced greater environmental pressure. Moreover, a branch-site model identified positively selected sites (PSSs) in ccsA, suggesting their involvement in the adaptation to sea environments. These findings are valuable for further investigations on Zosteraceae cpGenomes and will serve as an excellent resource for future studies on seagrass adaptation to sea environments.

Marker-Free Transplastomic Plants by Excision of Plastid Marker Genes Using Directly Repeated DNA Sequences

Excision of marker genes using DNA direct repeats makes use of the efficient native homologous recombination pathway present in the plastids of algae and plants. The method is simple, efficient, and widely applicable to plants and green algae. Marker excision frequency is dependent on the length and number of directly repeated sequences. When two repeats are used a repeat size of greater than 600 bp promotes efficient excision of the marker gene. A wide variety of sequences can be used to make the direct repeats. Only a single round of transformation is required and there is no requirement to introduce site-specific recombinases by retransformation or sexual crosses. Selection is used to maintain the marker and ensure homoplasmy of transgenic plastid genomes (plastomes). Release of selection allows the accumulation of marker-free plastomes generated by marker excision, which is a spontaneous and unidirectional process. Cytoplasmic sorting allows the segregation of cells with marker-free transgenic plastids. The marker-free shoots resulting from direct repeat mediated excision of marker genes have been isolated by vegetative propagation of shoots in the T₀ generation. Alternatively, accumulation of marker-free plastomes during growth, development and flowering of T₀ plants allows for the collection of seeds that give rise to a high proportion of marker-free T₁ seedlings. The procedure enables precise plastome engineering involving insertion of transgenes, point mutations and deletion of genes without the inclusion of any extraneous DNA. The simplicity and convenience of direct repeat excision facilitates its widespread use to isolate marker-free crops.

OsNBL1, a Multi-Organelle Localized Protein, Plays Essential Roles in Rice Senescence, Disease Resistance, and Salt Tolerance

BACKGROUND

Plant senescence is a complicated process involving multiple regulations, such as temperature, light, reactive oxygen species (ROS), endogenous hormone levels, and diseases. Although many such genes have been characterized to understand the process of leaf senescence, there still remain many unknowns, and many more genes need to be characterized.

RESULTS

We identified a rice mutant nbl1 with a premature leaf senescence phenotype. The causative gene, OsNBL1, encodes a small protein with 94 amino acids, which is conserved in monocot, as well as dicot plants. Disruption of OsNBL1 resulted in accelerated dark-induced leaf senescence, accompanied by a reduction in chlorophyll content and up-regulation of several senescence-associated genes. Notably, the nbl1 mutant was more susceptible to rice blast and bacterial blight but more tolerant to sodium chloride. Several salt-induced genes, including HAK1, HAK5, and three SNAC genes, were also up-regulated in the nbl1 mutant. Additionally, the nbl1 mutant was more sensitive to salicylic acid. Plants overexpressing OsNBL1 showed delayed dark-induced senescence, consistent with a higher chlorophyll content compared to wild-type plants. However, the overexpression plants were indistinguishable from the wild-types for resistance to the rice blast disease. OsNBL1 is a multi-organelle localized protein and interacts with OsClpP6, which is associated with senescence.

CONCLUSIONS

We described a novel leaf senescence mutant nbl1 in rice. It is showed that OsNBL1, a multi-organelle localized protein which interacts with a plastidic caseinolytic protease OsClpP6, is essential for controlling leaf senescence, disease resistance, and salt tolerance.

Plastid Genomes of Flowering Plants: Essential Principles

The plastid genome (plastome ) has proved a valuable source of data for evaluating evolutionary relationships among angiosperms. Through basic and applied approaches, plastid transformation technology offers the potential to understand and improve plant productivity, providing food, fiber, energy, and medicines to meet the needs of a burgeoning global population. The growing genomic resources available to both phylogenetic and biotechnological investigations is allowing novel insights and expanding the scope of plastome research to encompass new species. In this chapter, we present an overview of some of the seminal and contemporary research that has contributed to our current understanding of plastome evolution and attempt to highlight the relationship between evolutionary mechanisms and the tools of plastid genetic engineering.

A new method for functional analysis of plastid EMBRYO-DEFECTIVE PPR genes by efficiently constructing cosuppression lines in Arabidopsis

BACKGROUND

Pentatricopeptide-repeat proteins (PPRs) characterized by tandem arrays of a degenerate 35-amino-acid repeat (PPR motif) can bind a single strand RNA and regulate organelle gene expression at the post-transcriptional level, including RNA cleavage, splicing, editing and stability etc. PPRs are conserved in all eukaryotes and extremely expanded in higher plants. Many knockout mutants of PPR genes are embryonically lethal. These genes are named EMB PPRs and functional analysis of them is hindered by the difficulty in obtaining their knockout mutants.

RESULTS

Here, we report a new method for functional analysis of plastid EMB PPRs by efficiently constructing their cosuppression lines in Arabidopsis. When we overexpressed a mutated full length or truncated coding sequence (CDS) of EMB PPRs, such as EMB2279, EMB2654 and EMB976 (all belong to the P family PPRs) in the wild-type (WT) background, a large portion of T₁ plants displayed chlorosis phenotypes, which are similar to those of the weak allele mutants, knockdown lines or partially complementary lines. RT-PCR analysis showed that overexpression of the truncated EMB PPRs led to significant and specific downregulation of their corresponding endogenous mRNAs. However, when these EMB PPRs were overexpressed in the Post transcriptional Gene Silencing (PTGS) deficient mutant, RNA-dependent RNA polymerase 6 (rdr6), none of the T₁ plants displayed chlorosis phenotypes. These results indicate that the chlorosis phenotype results from post transcriptional silencing of the corresponding endogenous gene (also known as sense cosuppression).

CONCLUSIONS

Overexpression of an appropriately truncated EMB PPR CDS in WT leads to gene silencing in a RDR6-dependent manner, and this method can be employed to study the unknown function of EMB PPR genes. By this method, we showed that EMB976 is required for splicing of chloroplast clpP1 intron 2 and ycf3 intron 1.

Protective Roles of Cytosolic and Plastidal Proteasomes on Abiotic Stress and Pathogen Invasion

Protein malfunction is typically caused by abiotic stressors. To ensure cell survival during conditions of stress, it is important for plant cells to maintain proteins in their respective functional conformation. Self-compartmentalizing proteases, such as ATP-dependent Clp proteases and proteasomes are designed to act in the crowded cellular environment, and they are responsible for degradation of misfolded or damaged proteins within the cell. During different types of stress conditions, the levels of misfolded or orphaned proteins that are degraded by the 26S proteasome in the cytosol and nucleus and by the Clp proteases in the mitochondria and chloroplasts increase. This allows cells to uphold feedback regulations to cellular-level signals and adjust to altered environmental conditions. In this review, we summarize recent findings on plant proteolytic complexes with respect to their protective functions against abiotic and biotic stressors.

Insights Into Chloroplast Genome Evolution Across Opuntioideae (Cactaceae) Reveals Robust Yet Sometimes Conflicting Phylogenetic Topologies

Chloroplast genomes (plastomes) are frequently treated as highly conserved among land plants. However, many lineages of vascular plants have experienced extensive structural rearrangements, including inversions and modifications to the size and content of genes. Cacti are one of these lineages, containing the smallest plastome known for an obligately photosynthetic angiosperm, including the loss of one copy of the inverted repeat (∼25 kb) and the ndh gene suite, but only a few cacti from the subfamily Cactoideae have been sufficiently characterized. Here, we investigated the variation of plastome sequences across the second-major lineage of the Cactaceae, the subfamily Opuntioideae, to address (1) how variable is the content and arrangement of chloroplast genome sequences across the subfamily, and (2) how phylogenetically informative are the plastome sequences for resolving major relationships among the clades of Opuntioideae. Our de novo assembly of the Opuntia quimilo plastome recovered an organelle of 150,347 bp in length with both copies of the inverted repeat and the presence of all the ndh gene suite. An expansion of the large single copy unit and a reduction of the small single copy unit was observed, including translocations and inversion of genes, as well as the putative pseudogenization of some loci. Comparative analyses among all clades within Opuntioideae suggested that plastome structure and content vary across taxa of this subfamily, with putative independent losses of the ndh gene suite and pseudogenization of genes across disparate lineages, further demonstrating the dynamic nature of plastomes in Cactaceae. Our plastome dataset was robust in resolving three tribes with high support within Opuntioideae: Cylindropuntieae, Tephrocacteae and Opuntieae. However, conflicting topologies were recovered among major clades when exploring different assemblies of markers. A plastome-wide survey for highly informative phylogenetic markers revealed previously unused regions for future use in Sanger-based studies, presenting a valuable dataset with primers designed for continued evolutionary studies across Cactaceae. These results bring new insights into the evolution of plastomes in cacti, suggesting that further analyses should be carried out to address how ecological drivers, physiological constraints and morphological traits of cacti may be related with the common rearrangements in plastomes that have been reported across the family.

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.

The Mars1 kinase confers photoprotection through signaling in the chloroplast unfolded protein response

In response to proteotoxic stress, chloroplasts communicate with the nuclear gene expression system through a chloroplast unfolded protein response (cpUPR). We isolated Chlamydomonas reinhardtii mutants that disrupt cpUPR signaling and identified a gene encoding a previously uncharacterized cytoplasmic protein kinase, termed Mars1-for mutant affected in chloroplast-to-nucleus retrograde signaling-as the first known component in cpUPR signal transmission. Lack of cpUPR induction in MARS1 mutant cells impaired their ability to cope with chloroplast stress, including exposure to excessive light. Conversely, transgenic activation of cpUPR signaling conferred an advantage to cells undergoing photooxidative stress. Our results indicate that the cpUPR mitigates chloroplast photodamage and that manipulation of this pathway is a potential avenue for engineering photosynthetic organisms with increased tolerance to chloroplast stress.

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

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

Absence of photosynthetic state transitions in alien chloroplasts

MAIN CONCLUSION

The absence of state transitions in a Nt(Hn) cybrid is due to a cleavage of the threonine residue from the misprocessed N-terminus of the LHCII polypeptides. The cooperation between the nucleus and chloroplast genomes is essential for plant photosynthetic fitness. The rapid and specific interactions between nucleus-encoded and chloroplast-encoded proteins are under intense investigation with potential for applications in agriculture and renewable energy technology. Here, we present a novel model for photosynthesis research in which alien henbane (Hyoscyamus niger) chloroplasts function on the nuclear background of a tobacco (Nicotiana tabacum). The result of this coupling is a cytoplasmic hybrid (cybrid) with inhibited state transitions-a mechanism responsible for balancing energy absorption between photosystems. Protein analysis showed differences in the LHCII composition of the cybrid plants. SDS-PAGE analysis revealed a novel banding pattern in the cybrids with at least one additional 'LHCII' band compared to the wild-type parental species. Proteomic work suggested that the N-terminus of at least some of the cybrid Lhcb proteins was missing. These findings provide a mechanistic explanation for the lack of state transitions-the N-terminal truncation of the Lhcb proteins in the cybrid included the threonine residue that is phosphorylated/dephosphorylated in order to trigger state transitions and therefore crucial energy balancing mechanism in plants.

Extreme variation in rates of evolution in the plastid Clp protease complex

Eukaryotic cells represent an intricate collaboration between multiple genomes, even down to the level of multi-subunit complexes in mitochondria and plastids. One such complex in plants is the caseinolytic protease (Clp), which plays an essential role in plastid protein turnover. The proteolytic core of Clp comprises subunits from one plastid-encoded gene (clpP1) and multiple nuclear genes. TheclpP1 gene is highly conserved across most green plants, but it is by far the fastest evolving plastid-encoded gene in some angiosperms. To better understand these extreme and mysterious patterns of divergence, we investigated the history ofclpP1 molecular evolution across green plants by extracting sequences from 988 published plastid genomes. We find thatclpP1 has undergone remarkably frequent bouts of accelerated sequence evolution and architectural changes (e.g. a loss of introns andRNA-editing sites) within seed plants. AlthoughclpP1 is often assumed to be a pseudogene in such cases, multiple lines of evidence suggest that this is rarely true. We applied comparative native gel electrophoresis of chloroplast protein complexes followed by protein mass spectrometry in two species within the angiosperm genusSilene, which has highly elevated and heterogeneous rates ofclpP1 evolution. We confirmed thatclpP1 is expressed as a stable protein and forms oligomeric complexes with the nuclear-encoded Clp subunits, even in one of the most divergentSilene species. Additionally, there is a tight correlation between amino acid substitution rates inclpP1 and the nuclear-encoded Clp subunits across a broad sampling of angiosperms, suggesting continuing selection on interactions within this complex.

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 Analyses of Chloroplast Genomes of Cucurbitaceae Species: Lights into Selective Pressures and Phylogenetic Relationships

Cucurbitaceae is the fourth most important economic plant family with creeping herbaceous species mainly distributed in tropical and subtropical regions. Here, we described and compared the complete chloroplast genome sequences of ten representative species from Cucurbitaceae. The lengths of the ten complete chloroplast genomes ranged from 155,293 bp (C. sativus) to 158,844 bp (M. charantia), and they shared the most common genomic features. 618 repeats of three categories and 813 microsatellites were found. Sequence divergence analysis showed that the coding and IR regions were highly conserved. Three protein-coding genes (accD, clpP, and matK) were under selection and their coding proteins often have functions in chloroplast protein synthesis, gene transcription, energy transformation, and plant development. An unconventional translation initiation codon of psbL gene was found and provided evidence for RNA editing. Applying BI and ML methods, phylogenetic analysis strongly supported the position of Gomphogyne, Hemsleya, and Gynostemma as the relatively original lineage in Cucurbitaceae. This study suggested that the complete chloroplast genome sequences were useful for phylogenetic studies. It would also determine potential molecular markers and candidate DNA barcodes for coming studies and enrich the valuable complete chloroplast genome resources of Cucurbitaceae.

Positive Selection Driving Cytoplasmic Genome Evolution of the Medicinally Important Ginseng Plant Genus Panax

Panax L. (the ginseng genus) is a shade-demanding group within the family Araliaceae and all of its species are of crucial significance in traditional Chinese medicine. Phylogenetic and biogeographic analyses demonstrated that two rounds of whole genome duplications accompanying with geographic and ecological isolations promoted the diversification of Panax species. However, contributions of the cytoplasmic genomes to the adaptive evolution of Panax species remained largely uninvestigated. In this study, we sequenced the chloroplast and mitochondrial genomes of 11 accessions belonging to seven Panax species. Our results show that heterogeneity in nucleotide substitution rate is abundant in both of the two cytoplasmic genomes, with the mitochondrial genome possessing more variants at the total level but the chloroplast showing higher sequence polymorphisms at the genic regions. Genome-wide scanning of positive selection identified five and 12 genes from the chloroplast and mitochondrial genomes, respectively. Functional analyses further revealed that these selected genes play important roles in plant development, cellular metabolism and adaptation. We therefore conclude that positive selection might be one of the potential evolutionary forces that shaped nucleotide variation pattern of these Panax species. In particular, the mitochondrial genes evolved under stronger selective pressure compared to the chloroplast genes.

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

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

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.

Recombination-dependent replication and gene conversion homogenize repeat sequences and diversify plastid genome structure

PREMISE OF THE STUDY

There is a misinterpretation in the literature regarding the variable orientation of the small single copy region of plastid genomes (plastomes). The common phenomenon of small and large single copy inversion, hypothesized to occur through intramolecular recombination between inverted repeats (IR) in a circular, single unit-genome, in fact, more likely occurs through recombination-dependent replication (RDR) of linear plastome templates. If RDR can be primed through both intra- and intermolecular recombination, then this mechanism could not only create inversion isomers of so-called single copy regions, but also an array of alternative sequence arrangements.

METHODS

We used Illumina paired-end and PacBio single-molecule real-time (SMRT) sequences to characterize repeat structure in the plastome of Monsonia emarginata (Geraniaceae). We used OrgConv and inspected nucleotide alignments to infer ancestral nucleotides and identify gene conversion among repeats and mapped long (>1 kb) SMRT reads against the unit-genome assembly to identify alternative sequence arrangements.

RESULTS

Although M. emarginata lacks the canonical IR, we found that large repeats (>1 kilobase; kb) represent ∼22% of the plastome nucleotide content. Among the largest repeats (>2 kb), we identified GC-biased gene conversion and mapping filtered, long SMRT reads to the M. emarginata unit-genome assembly revealed alternative, substoichiometric sequence arrangements.

CONCLUSION

We offer a model based on RDR and gene conversion between long repeated sequences in the M. emarginata plastome and provide support that both intra-and intermolecular recombination between large repeats, particularly in repeat-rich plastomes, varies unit-genome structure while homogenizing the nucleotide sequence of repeats.

Functional diversity of Arabidopsis organelle-localized RNA-recognition motif-containing proteins

RNA-Binding Proteins (RBPs) play key roles in plant gene expression and regulation. RBPs contain a variety of RNA-binding motifs, the most abundant and most widespread one in eukaryotes is the RNA recognition motif (RRM). Many nucleus-encoded RRM-containing proteins are transported into chloroplasts and/or mitochondria, and participate in various RNA-related processes in plant organelles. Loss of these proteins can have a detrimental effect on some critical processes such as photosynthesis and respiration, sometimes leading to lethality. Progress has been made in the last few years in understanding the function of particular organelle-localized RRM-containing proteins. Members of the Organelle RRM protein (ORRM, some also characterized as Glycine-Rich RNA-Binding Proteins) family and the Chloroplast RiboNucleoProtein (cpRNP) family, are involved in various types of RNA metabolism, including RNA editing, RNA stability and RNA processing. Organelle-localized RRM proteins also function in plant development and stress responses, in some conditions acting as protein or RNA chaperones. There has been recent progress in characterizing the function of organelle-localized RRM proteins in RNA-related processes and how RRM proteins contribute to the normal growth and development of plants. WIREs RNA 2017, 8:e1420. doi: 10.1002/wrna.1420 For further resources related to this article, please visit the WIREs website.

The RNA recognition motif protein CP33A is a global ligand of chloroplast mRNAs and is essential for plastid biogenesis and plant development

Chloroplast RNA metabolism depends on a multitude of nuclear-encoded RNA-binding proteins (RBPs). Most known chloroplast RBPs address specific RNA targets and RNA-processing functions. However, members of the small chloroplast ribonucleoprotein family (cpRNPs) play a global role in processing and stabilizing chloroplast RNAs. Here, we show that the cpRNP CP33A localizes to a distinct sub-chloroplastic domain and is essential for chloroplast development. The loss of CP33A yields albino seedlings that exhibit aberrant leaf development and can only survive in the presence of an external carbon source. Genome-wide RNA association studies demonstrate that CP33A associates with all chloroplast mRNAs. For a given transcript, quantification of CP33A-bound versus free RNAs demonstrates that CP33A associates with the majority of most mRNAs analyzed. Our results further show that CP33A is required for the accumulation of a number of tested mRNAs, and is particularly relevant for unspliced and unprocessed precursor mRNAs. Finally, CP33A fails to associate with polysomes or to strongly co-precipitate with ribosomal RNA, suggesting that it defines a ribodomain that is separate from the chloroplast translation machinery. Collectively, these findings suggest that CP33A contributes to globally essential RNA processes in the chloroplasts of higher plants.

Analysis of Arabidopsis Accessions Hypersensitive to a Loss of Chloroplast Translation

Natural accessions of Arabidopsis (Arabidopsis thaliana) differ in their ability to tolerate a loss of chloroplast translation. These differences can be attributed in part to variation in a duplicated nuclear gene (ACC2) that targets homomeric acetyl-coenzyme A carboxylase (ACCase) to plastids. This functional redundancy allows limited fatty acid biosynthesis to occur in the absence of heteromeric ACCase, which is encoded in part by the plastid genome. In the presence of functional ACC2, tolerant alleles of several nuclear genes, not yet identified, enhance the growth of seedlings and embryos disrupted in chloroplast translation. ACC2 knockout mutants, by contrast, are hypersensitive. Here we describe an expanded search for hypersensitive accessions of Arabidopsis, evaluate whether all of these accessions are defective in ACC2, and characterize genotype-to-phenotype relationships for homomeric ACCase variants identified among 855 accessions with sequenced genomes. Null alleles with ACC2 nonsense mutations, frameshift mutations, small deletions, genomic rearrangements, and defects in RNA splicing are included among the most sensitive accessions examined. By contrast, most missense mutations affecting highly conserved residues failed to eliminate ACC2 function. Several accessions were identified where sensitivity could not be attributed to a defect in either ACC2 or Tic20-IV, the chloroplast membrane channel required for ACC2 uptake. Overall, these results underscore the central role of ACC2 in mediating Arabidopsis response to a loss of chloroplast translation, highlight future applications of this system to analyzing chloroplast protein import, and provide valuable insights into the mutational landscape of an important metabolic enzyme that is highly conserved throughout eukaryotes.

The E3 ligase AtCHIP positively regulates Clp proteolytic subunit homeostasis

The caseinolytic peptidase (Clp) core proteins are essential for plant growth and development, especially for chloroplast function. Antisense or overexpression of ClpP4, which is one of the Clp core subunits, causes chlorotic phenotypes in Arabidopsis. An E3 ligase gene, AtCHIP, has previously been found to ubiquitylate ClpP4 in vitro. ClpP4 antisense and overexpressing plants that also overexpressed AtCHIP were constructed to explore the effect of AtCHIP on ClpP4. Overexpression of AtCHIP was found to rescue the chlorotic phenotypes of both ClpP4 antisense and overexpressing plants. The unbalanced levels of Clp core proteins in ClpP4 antisense and overexpressing plants with overexpression of AtCHIP were similar to wild-type levels, suggesting that AtCHIP regulates Clp core proteins. The results also show that AtCHIP can interact with ClpP3 and ClpP5 in yeast and ubiquitylate ClpP3 and ClpP5 in vitro. This suggests that AtCHIP is directly related to ClpP3 and ClpP5. Given these results, the inference is that through selective degradation of Clp subunits, AtCHIP could positively regulate homeostasis of Clp proteolytic subunits and maximize the production of functional chloroplasts. Similar results were obtained from transgenic tobacco plants, suggesting that regulation of the Clp protease by AtCHIP is conserved.

The Complete Sequence of the Acacia ligulata Chloroplast Genome Reveals a Highly Divergent clpP1 Gene

Legumes are a highly diverse angiosperm family that include many agriculturally important species. To date, 21 complete chloroplast genomes have been sequenced from legume crops confined to the Papilionoideae subfamily. Here we report the first chloroplast genome from the Mimosoideae, Acacia ligulata, and compare it to the previously sequenced legume genomes. The A. ligulata chloroplast genome is 158,724 bp in size, comprising inverted repeats of 25,925 bp and single-copy regions of 88,576 bp and 18,298 bp. Acacia ligulata lacks the inversion present in many of the Papilionoideae, but is not otherwise significantly different in terms of gene and repeat content. The key feature is its highly divergent clpP1 gene, normally considered essential in chloroplast genomes. In A. ligulata, although transcribed and spliced, it probably encodes a catalytically inactive protein. This study provides a significant resource for further genetic research into Acacia and the Mimosoideae. The divergent clpP1 gene suggests that Acacia will provide an interesting source of information on the evolution and functional diversity of the chloroplast Clp protease complex.

Complete chloroplast genome of the multifunctional crop globe artichoke and comparison with other Asteraceae

With over 20,000 species, Asteraceae is the second largest plant family. High-throughput sequencing of nuclear and chloroplast genomes has allowed for a better understanding of the evolutionary relationships within large plant families. Here, the globe artichoke chloroplast (cp) genome was obtained by a combination of whole-genome and BAC clone high-throughput sequencing. The artichoke cp genome is 152,529 bp in length, consisting of two single-copy regions separated by a pair of inverted repeats (IRs) of 25,155 bp, representing the longest IRs found in the Asteraceae family so far. The large (LSC) and the small (SSC) single-copy regions span 83,578 bp and 18,641 bp, respectively. The artichoke cp sequence was compared to the other eight Asteraceae complete cp genomes available, revealing an IR expansion at the SSC/IR boundary. This expansion consists of 17 bp of the ndhF gene generating an overlap between the ndhF and ycf1 genes. A total of 127 cp simple sequence repeats (cpSSRs) were identified in the artichoke cp genome, potentially suitable for future population studies in the Cynara genus. Parsimony-informative regions were evaluated and allowed to place a Cynara species within the Asteraceae family tree. The eight most informative coding regions were also considered and tested for “specific barcode” purpose in the Asteraceae family. Our results highlight the usefulness of cp genome sequencing in exploring plant genome diversity and retrieving reliable molecular resources for phylogenetic and evolutionary studies, as well as for specific barcodes in plants.

Accumulation of high contents of free amino acids in the leaves of Nicotiana benthamiana by the co-suppression of NbClpC1 and NbClpC2 genes

KEY MESSAGE

We report the significant increase of the content of free amino acids in Nicotiana benthamiana by the co-suppression of the ClpC1 and ClpC2 genes, which are translated to be the chaperonic part in the Clp protease at plastids. Clp protease with ClpC1 and ClpC2 proteins as the chaperonic part degrades denatured or improperly folded protein in plastids. Nicotiana benthamiana ClpC1 and ClpC2 genes (NbClpC1 and NbClpC2: NbClpC1/C2) share 93% similarities; therefore, co-suppression of the NbClpC1/C2 was possible using a single virus-induced silencing vector. Co-suppression of NbClpC1/C2 resulted in a pleiotropic phenotype including disappearance of apical dominance and formation of chlorotic leaves. NbClpC1/C2 co-suppressed leaves accumulated 11.9-fold more free amino acids than the GFP-silenced leaves. The co-suppression of NbClpC1/C2 did not change the expression levels of some selected genes in the biosynthetic pathways for the free amino acids, but reduced the total protein amounts to 32.5%, indicating that co-suppression affected the incorporation of free amino acids in proteins during translation. The loosely packed mesophyll cells and abnormal vascular bundles in the leaves suggested structural problems associated with translocation of free amino acids to sink tissues. NbClpC1/C2 co-suppression can offer a novel strategy for accumulation of free amino acids though it results in stunted growth.

Exploring the limits for reduction of plastid genomes: a case study of the mycoheterotrophic orchids Epipogium aphyllum and Epipogium roseum

The question on the patterns and limits of reduction of plastid genomes in nonphotosynthetic plants and the reasons of their conservation is one of the intriguing topics in plant genome evolution. Here, we report sequencing and analysis of plastid genome in nonphotosynthetic orchids Epipogium aphyllum and Epipogium roseum, which, with sizes of 31 and 19 kbp, respectively, represent the smallest plastid genomes characterized by now. Besides drastic reduction, which is expected, we found several unusual features of these “minimal” plastomes: Multiple rearrangements, highly biased nucleotide composition, and unprecedentedly high substitution rate. Only 27 and 29 genes remained intact in the plastomes of E. aphyllum and E. roseum—those encoding ribosomal components, transfer RNAs, and three additional housekeeping genes (infA, clpP, and accD). We found no signs of relaxed selection acting on these genes. We hypothesize that the main reason for retention of plastid genomes in Epipogium is the necessity to translate messenger RNAs (mRNAs) of accD and/or clpP proteins which are essential for cell metabolism. However, these genes are absent in plastomes of several plant species; their absence is compensated by the presence of a functional copy arisen by gene transfer from plastid to the nuclear genome. This suggests that there is no single set of plastid-encoded essential genes, but rather different sets for different species and that the retention of a gene in the plastome depends on the interaction between the nucleus and plastids.

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

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

Organization, function and substrates of the essential Clp protease system in plastids

Intra-plastid proteolysis is essential in plastid biogenesis, differentiation and plastid protein homeostasis (proteostasis). We provide a comprehensive review of the Clp protease system present in all plastid types and we draw lessons from structural and functional information of bacterial Clp systems. The Clp system plays a central role in plastid development and function, through selective removal of miss-folded, aggregated, or otherwise unwanted proteins. The Clp system consists of a tetradecameric proteolytic core with catalytically active ClpP and inactive ClpR subunits, hexameric ATP-dependent chaperones (ClpC,D) and adaptor protein(s) (ClpS1) enhancing delivery of subsets of substrates. Many structural and functional features of the plastid Clp system are now understood though extensive reverse genetics analysis combined with biochemical analysis, as well as large scale quantitative proteomics for loss-of-function mutants of Clp core, chaperone and ClpS1 subunits. Evolutionary diversification of Clp system across non-photosynthetic and photosynthetic prokaryotes and organelles is illustrated. Multiple substrates have been suggested based on their direct interaction with the ClpS1 adaptor or screening of different loss-of-function protease mutants. The main challenge is now to determine degradation signals (degrons) in Clp substrates and substrate delivery mechanisms, as well as functional interactions of Clp with other plastid proteases. This article is part of a Special Issue entitled: Chloroplast Biogenesis.

Conditional repression of essential chloroplast genes: Evidence for new plastid signaling pathways

The development of a repressible chloroplast gene expression system in Chlamydomonas reinhardtii has opened the door for studying the role of essential chloroplast genes. This approach has been used to analyze three chloroplast genes of this sort coding for the α subunit of RNA polymerase (rpoA), a ribosomal protein (rps12) and the catalytic subunit of the ATP-dependent ClpP protease (clpP1). Depletion of the three corresponding proteins leads to growth arrest and cell death. Shutdown of chloroplast transcription and translation increases the abundance of a set of plastid transcripts that includes mainly those involved in transcription, translation and proteolysis and reveals multiple regulatory feedback loops in the chloroplast gene circuitry. Depletion of ClpP profoundly affects plastid protein homeostasis and elicits an autophagy-like response with extensive cytoplasmic vacuolization of cells. It also triggers changes in chloroplast and nuclear gene expression resulting in increased abundance of chaperones, proteases, ubiquitin-related proteins and proteins involved in lipid trafficking and thylakoid biogenesis. These features are hallmarks of an unfolded protein response in the chloroplast and raise new questions on plastid protein homeostasis and plastid signaling. This article is part of a Special Issue entitled: Chloroplast Biogenesis.

Natural variation in sensitivity to a loss of chloroplast translation in Arabidopsis

Mutations that eliminate chloroplast translation in Arabidopsis (Arabidopsis thaliana) result in embryo lethality. The stage of embryo arrest, however, can be influenced by genetic background. To identify genes responsible for improved growth in the absence of chloroplast translation, we examined seedling responses of different Arabidopsis accessions on spectinomycin, an inhibitor of chloroplast translation, and crossed the most tolerant accessions with embryo-defective mutants disrupted in chloroplast ribosomal proteins generated in a sensitive background. The results indicate that tolerance is mediated by ACC2, a duplicated nuclear gene that targets homomeric acetyl-coenzyme A carboxylase to plastids, where the multidomain protein can participate in fatty acid biosynthesis. In the presence of functional ACC2, tolerance is enhanced by a second locus that maps to chromosome 5 and heightened by additional genetic modifiers present in the most tolerant accessions. Notably, some of the most sensitive accessions contain nonsense mutations in ACC2, including the "Nossen" line used to generate several of the mutants studied here. Functional ACC2 protein is therefore not required for survival in natural environments, where heteromeric acetyl-coenzyme A carboxylase encoded in part by the chloroplast genome can function instead. This work highlights an interesting example of a tandem gene duplication in Arabidopsis, helps to explain the range of embryo phenotypes found in Arabidopsis mutants disrupted in essential chloroplast functions, addresses the nature of essential proteins encoded by the chloroplast genome, and underscores the value of using natural variation to study the relationship between chloroplast translation, plant metabolism, protein import, and plant development.

A pair of homoeolog ClpP5 genes underlies a virescent yellow-like mutant and its modifier in maize

Gene-background interaction is a commonly observed phenomenon in many species, but the molecular mechanisms of such an interaction is less well understood. Here we report the cloning of a maize mutant gene and its modifier. A recessive mutant with a virescent yellow-like (vyl) phenotype was identified in an ethyl methanesulfonate-mutagenized population derived from the maize inbred line B73. Homozygous mutant maize plants exhibited a yellow leaf phenotype after emergence but gradually recovered and became indistinguishable from wild-type plants after approximately 2 weeks. Taking the positional cloning approach, the Chr.9_ClpP5 gene, one of the proteolytic subunits of the chloroplast Clp protease complex, was identified and validated as the candidate gene for vyl. When introgressed by backcross into the maize inbred line PH09B, the mutant phenotype of vyl lasted much longer in the greenhouse and was lethal in the field, implying the presence of a modifier(s) for vyl. A major modifier locus was identified on chromosome 1, and a paralogous ClpP5 gene was isolated and confirmed as the candidate for the vyl-modifier. Expression of Chr.1_ClpP5 is induced significantly in B73 by the vyl mutation, while the expression of Chr.1_ClpP5 in PH09B is not responsive to the vyl mutation. Moreover, expression and sequence analysis suggests that the PH09B Chr.1_ClpP5 allele is functionally weaker than the B73 allele. We propose that functional redundancy between duplicated paralogous genes is the molecular mechanism for the interaction between vyl and its modifier.