Uncovering the protein translocon at the chloroplast inner envelope membrane

Removal of the large inverted repeat from the plastid genome reveals gene dosage effects and leads to increased genome copy number

The chloroplast genomes of most plants and algae contain a large inverted repeat (IR) region that separates two single-copy regions and harbours the ribosomal RNA operon. We have addressed the functional importance of the IR region by removing an entire copy of the 25.3-kb IR from the tobacco plastid genome. Using plastid transformation and subsequent selectable marker gene elimination, we precisely excised the IR, thus generating plants with a substantially reduced plastid genome size. We show that the lack of the IR results in a mildly reduced plastid ribosome number, suggesting a gene dosage benefit from the duplicated presence of the ribosomal RNA operon. Moreover, the IR deletion plants contain an increased number of plastid genomes, suggesting that genome copy number is regulated by measuring total plastid DNA content rather than by counting genomes. Together, our findings (1) demonstrate that the IR can enhance the translation capacity of the plastid, (2) reveal the relationship between genome size and genome copy number, and (3) provide a simplified plastid genome structure that will facilitate future synthetic biology applications.

Polyploid plants take cytonuclear perturbations in stride

Hybridization in plants is often accompanied by nuclear genome doubling (allopolyploidy), which has been hypothesized to perturb interactions between nuclear and organellar (mitochondrial and plastid) genomes by creating imbalances in the relative copy number of these genomes and producing genetic incompatibilities between maternally derived organellar genomes and the half of the allopolyploid nuclear genome from the paternal progenitor. Several evolutionary responses have been predicted to ameliorate these effects, including selection for changes in protein sequences that restore cytonuclear interactions; biased gene retention/expression/conversion favoring maternal nuclear gene copies; and fine-tuning of relative cytonuclear genome copy numbers and expression levels. Numerous recent studies, however, have found that evolutionary responses are inconsistent and rarely scale to genome-wide generalities. The apparent robustness of plant cytonuclear interactions to allopolyploidy may reflect features that are general to allopolyploids such as the lack of F2 hybrid breakdown under disomic inheritance, and others that are more plant-specific, including slow sequence divergence in organellar genomes and preexisting regulatory responses to changes in cell size and endopolyploidy during development. Thus, cytonuclear interactions may only rarely act as the main barrier to establishment of allopolyploid lineages, perhaps helping to explain why allopolyploidy is so pervasive in plant evolution.

A comparative analysis of plastome evolution in autotrophic Piperales

PREMISE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

The translocon protein FtsHi1 is an ATP-dependent DNA/RNA helicase that prevents R-loop accumulation in chloroplasts

R-loops, three-stranded nucleic acid structures consisting of a DNA: RNA hybrid and displaced single-stranded DNA, play critical roles in gene expression and genome stability. How R-loop homeostasis is integrated into chloroplast gene expression remains largely unknown. We found an unexpected function of FtsHi1, an inner envelope membrane-bound AAA-ATPase in chloroplast R-loop homeostasis of Arabidopsis thaliana. Previously, this protein was shown to function as a component of the import motor complex for nuclear-encoded chloroplast proteins. However, this study provides evidence that FtsHi1 is an ATP-dependent helicase that efficiently unwinds both DNA-DNA and DNA-RNA duplexes, thereby preventing R-loop accumulation. Over-accumulation of R-loops could impair chloroplast transcription but not necessarily genome integrity. The dual function of FtsHi1 in both protein import and chloroplast gene expression may be important to coordinate the biogenesis of nuclear- and chloroplast-encoded subunits of multi-protein photosynthetic complexes. This study suggests a mechanical link between protein import and R-loop homeostasis in chloroplasts of higher plants.

Plastid genome of Passiflora tripartita var. mollissima (poro-poro) from Huánuco, Peru

Passiflora tripartita var. mollissima, known locally as poro-poro, is an important native fruit used in traditional Peruvian medicine with relevant agro-industrial and pharmaceutical potential for its antioxidant capacity for human health. However, to date, only a few genetic data are available, which limits exploring its genetic diversity and developing new genetic studies for its improvement. We report the poro-poro plastid genome to expand the knowledge of its molecular markers, evolutionary studies, molecular pathways, and conservation genetics. The complete chloroplast (cp) genome is 163,451 bp in length with a typical quadripartite structure, containing a large single-copy region of 85,525 bp and a small single-copy region of 13,518 bp, separated by a pair of inverted repeat regions (IR) of 32,204 bp, and the overall GC content was 36.87%. This cp genome contains 128 genes (110 genes were unique and 18 genes were found duplicated in each IR region), including 84 protein-coding genes, 36 transfer RNA-coding genes, eight ribosomal RNA-coding genes, and 13 genes with introns (11 genes with one intron and two genes with two introns). The inverted repeat region boundaries among species were similar in organization, gene order, and content, with a few revisions. The phylogenetic tree reconstructed based on single-copy orthologous genes and maximum likelihood analysis demonstrates poro-poro is most closely related to Passiflora menispermifolia and Passiflora oerstedii. In summary, our study constitutes a valuable resource for studying molecular evolution, phylogenetics, and domestication. It also provides a powerful foundation for conservation genetics research and plant breeding programs. To our knowledge, this is the first report on the plastid genome of Passiflora tripartita var. mollissima from Peru.

Light controls mesophyll-specific post-transcriptional splicing of photoregulatory genes by AtPRMT5

Light plays a central role in plant growth and development, providing an energy source and governing various aspects of plant morphology. Previous study showed that many polyadenylated full-length RNA molecules within the nucleus contain unspliced introns (post-transcriptionally spliced introns, PTS introns), which may play a role in rapidly responding to changes in environmental signals. However, the mechanism underlying post-transcriptional regulation during initial light exposure of young, etiolated seedlings remains elusive. In this study, we used FLEP-seq2, a Nanopore-based sequencing technique, to analyze nuclear RNAs in Arabidopsis (Arabidopsis thaliana) seedlings under different light conditions and found numerous light-responsive PTS introns. We also used single-nucleus RNA sequencing (snRNA-seq) to profile transcripts in single nucleus and investigate the distribution of light-responsive PTS introns across distinct cell types. We established that light-induced PTS introns are predominant in mesophyll cells during seedling de-etiolation following exposure of etiolated seedlings to light. We further demonstrated the involvement of the splicing-related factor A. thaliana PROTEIN ARGININE METHYLTRANSFERASE 5 (AtPRMT5), working in concert with the E3 ubiquitin ligase CONSTITUTIVE PHOTOMORPHOGENIC 1 (COP1), a critical repressor of light signaling pathways. We showed that these two proteins orchestrate light-induced PTS events in mesophyll cells and facilitate chloroplast development, photosynthesis, and morphogenesis in response to ever-changing light conditions. These findings provide crucial insights into the intricate mechanisms underlying plant acclimation to light at the cell-type level.

Plastome phylogenomics and morphological traits analyses provide new insights into the phylogenetic position, species delimitation and speciation of Triplostegia (Caprifoliaceae)

BACKGROUND

The genus Triplostegia contains two recognized species, T. glandulifera and T. grandiflora, but its phylogenetic position and species delimitation remain controversial. In this study, we assembled plastid genomes and nuclear ribosomal DNA (nrDNA) cistrons sampled from 22 wild Triplostegia individuals, each from a separate population, and examined these with 11 recently published Triplostegia plastomes. Morphological traits were measured from herbarium specimens and wild material, and ecological niche models were constructed.

RESULTS

Triplostegia is a monophyletic genus within the subfamily Dipsacoideae comprising three monophyletic species, T. glandulifera, T. grandiflora, and an unrecognized species Triplostegia sp. A, which occupies much higher altitude than the other two. The new species had previously been misidentified as T. glandulifera, but differs in taproot, leaf, and other characters. Triplotegia is an old genus, with stem age 39.96 Ma, and within it T. glandulifera diverged 7.94 Ma. Triplostegia grandiflora and sp. A diverged 1.05 Ma, perhaps in response to Quaternary climate fluctuations. Niche overlap between Triplostegia species was positively correlated with their phylogenetic relatedness.

CONCLUSIONS

Our results provide new insights into the species delimitation of Triplostegia, and indicate that a taxonomic revision of Triplostegia is needed. We also identified that either rpoB-trnC or ycf1 could serve as a DNA barcode for Triplostegia.

Herbgenomics meets Papaveraceae: a promising -omics perspective on medicinal plant research

Herbal medicines were widely used in ancient and modern societies as remedies for human ailments. Notably, the Papaveraceae family includes well-known species, such as Papaver somniferum and Chelidonium majus, which possess medicinal properties due to their latex content. Latex-bearing plants are a rich source of diverse bioactive compounds, with applications ranging from narcotics to analgesics and relaxants. With the advent of high-throughput technologies and advancements in sequencing tools, an opportunity exists to bridge the knowledge gap between the genetic information of herbs and the regulatory networks underlying their medicinal activities. This emerging discipline, known as herbgenomics, combines genomic information with other -omics studies to unravel the genetic foundations, including essential gene functions and secondary metabolite biosynthesis pathways. Furthermore, exploring the genomes of various medicinal plants enables the utilization of modern genetic manipulation techniques, such as Clustered Regularly-Interspaced Short Palindromic Repeats (CRISPR/Cas9) or RNA interference. This technological revolution has facilitated systematic studies of model herbs, targeted breeding of medicinal plants, the establishment of gene banks and the adoption of synthetic biology approaches. In this article, we provide a comprehensive overview of the recent advances in genomic, transcriptomic, proteomic and metabolomic research on species within the Papaveraceae family. Additionally, it briefly explores the potential applications and key opportunities offered by the -omics perspective in the pharmaceutical industry and the agrobiotechnology field.

Chloroplast protein translocation pathways and ubiquitin-dependent regulation at a glance

Chloroplasts conduct photosynthesis and numerous metabolic and signalling processes that enable plant growth and development. Most of the ∼3000 proteins in chloroplasts are nucleus encoded and must be imported from the cytosol. Thus, the protein import machinery of the organelle (the TOC-TIC apparatus) is of fundamental importance for chloroplast biogenesis and operation. Cytosolic factors target chloroplast precursor proteins to the TOC-TIC apparatus, which drives protein import across the envelope membranes into the organelle, before various internal systems mediate downstream routing to different suborganellar compartments. The protein import system is proteolytically regulated by the ubiquitin-proteasome system (UPS), enabling centralized control over the organellar proteome. In addition, the UPS targets a range of chloroplast proteins directly. In this Cell Science at a Glance article and the accompanying poster, we present mechanistic details of these different chloroplast protein targeting and translocation events, and of the UPS systems that regulate chloroplast proteins.

Genomic Analysis of Plastid-Nuclear Interactions and Differential Evolution Rates in Coevolved Genes across Juglandaceae Species

The interaction between the nuclear and chloroplast genomes in plants is crucial for preserving essential cellular functions in the face of varying rates of mutation, levels of selection, and modes of transmission. Despite this, identifying nuclear genes that coevolve with chloroplast genomes at a genome-wide level has remained a challenge. In this study, we conducted an evolutionary rate covariation analysis to identify candidate nuclear genes coevolving with chloroplast genomes in Juglandaceae. Our analysis was based on 4,894 orthologous nuclear genes and 76 genes across seven chloroplast partitions in nine Juglandaceae species. Our results indicated that 1,369 (27.97%) of the nuclear genes demonstrated signatures of coevolution, with the Ycf1/2 partition yielding the largest number of hits (765) and the ClpP1 partition yielding the fewest (13). These hits were found to be significantly enriched in biological processes related to leaf development, photoperiodism, and response to abiotic stress. Among the seven partitions, AccD, ClpP1, MatK, and RNA polymerase partitions and their respective hits exhibited a narrow range, characterized by dN/dS values below 1. In contrast, the Ribosomal, Photosynthesis, Ycf1/2 partitions and their corresponding hits, displayed a broader range of dN/dS values, with certain values exceeding 1. Our findings highlight the differences in the number of candidate nuclear genes coevolving with the seven chloroplast partitions in Juglandaceae species and the correlation between the evolution rates of these genes and their corresponding chloroplast partitions.

Characterization of the complete chloroplast genome of Gypsophila huashanensis Y. W. Tsui & D. Q. Lu, an endemic herb species in China

Gypsophila huashanensis Y. W. Tsui & D. Q. Lu (Caryophyllaceae) is an endemic herb species to the Qinling Mountains in China. In this study, we characterized its whole plastid genome using the Illumina sequencing platform. The complete plastid genome of G. huashanensis is 152,457 bp in length, including a large single-copy DNA region of 83,476 bp, a small single-copy DNA region of 17,345 bp, and a pair of inverted repeat DNA sequences of 25,818 bp. The genome contains 130 genes comprising 85 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Evolutionary analysis showed that the non-coding regions of Caryophyllaceae exhibit a higher level of divergence than the exon regions. Gene site selection analysis suggested that 11 coding protein genes (accD, atpF, ndhA, ndhB, petB, petD, rpoCl, rpoC2, rps16, ycfl, and ycf2) have some sites under protein sequence evolution. Phylogenetic analysis showed that G. huashanensis is most closely related to the congeneric species G. oldhamiana. These results are very useful for studying phylogenetic evolution and species divergence in the family Caryophyllaceae.

The journey of preproteins across the chloroplast membrane systems

The photosynthetic capacity of chloroplasts is vital for autotrophic growth in algae and plants. The origin of the chloroplast has been explained by the endosymbiotic theory that proposes the engulfment of a cyanobacterium by an ancestral eukaryotic cell followed by the transfer of many cyanobacterial genes to the host nucleus. As a result of the gene transfer, the now nuclear-encoded proteins acquired chloroplast targeting peptides (known as transit peptides; transit peptide) and are translated as preproteins in the cytosol. Transit peptides contain specific motifs and domains initially recognized by cytosolic factors followed by the chloroplast import components at the outer and inner envelope of the chloroplast membrane. Once the preprotein emerges on the stromal side of the chloroplast protein import machinery, the transit peptide is cleaved by stromal processing peptidase. In the case of thylakoid-localized proteins, cleavage of the transit peptides may expose a second targeting signal guiding the protein to the thylakoid lumen or allow insertion into the thylakoid membrane by internal sequence information. This review summarizes the common features of targeting sequences and describes their role in routing preproteins to and across the chloroplast envelope as well as the thylakoid membrane and lumen.

The First Complete Chloroplast Genome of Cordia monoica: Structure and Comparative Analysis

Cordia monoica is a member of the Boraginaceae family. This plant is widely distributed in tropical regions and has a great deal of medical value as well as economic importance. In the current study, the complete chloroplast (cp) genome of C. monoica was sequenced, assembled, annotated, and reported. This circular chloroplast genome had a size of 148,711 bp, with a quadripartite structure alternating between a pair of repeated inverted regions (26,897-26,901 bp) and a single copy region (77,893 bp). Among the 134 genes encoded by the cp genome, there were 89 protein-coding genes, 37 transfer RNA (tRNA) genes, and 8 ribosomal RNA (rRNA) genes. A total of 1387 tandem repeats were detected, with the hexanucleotides class making up 28 percent of the repeats. Cordia monoica has 26,303 codons in its protein-coding regions, and leucine amino acid was the most frequently encoded amino acid in contrast to cysteine. In addition, 12 of the 89 protein-coding genes were found to be under positive selection. The phyloplastomic taxonomical clustering of the Boraginaceae species provides further evidence that chloroplast genome data are reliable not only at family level but also in deciphering the phylogeny at genus level (e.g., Cordia).

Plastids: diving into their diversity, their functions, and their role in plant development

Plastids are a group of essential, heterogenous semi-autonomous organelles characteristic of plants that perform photosynthesis and a diversity of metabolic pathways that impact growth and development. Plastids are remarkably dynamic and can interconvert in response to specific developmental and environmental cues, functioning as a central metabolic hub in plant cells. By far the best studied plastid is the chloroplast, but in recent years the combination of modern techniques and genetic analyses has expanded our current understanding of plastid morphological and functional diversity in both model and non-model plants. These studies have provided evidence of an unexpected diversity of plastid subtypes with specific characteristics. In this review, we describe recent findings that provide insights into the characteristics of these specialized plastids and their functions. We concentrate on the emerging evidence that supports the model that signals derived from particular plastid types play pivotal roles in plant development, environmental, and defense responses. Furthermore, we provide examples of how new technologies are illuminating the functions of these specialized plastids and the overall complexity of their differentiation processes. Finally, we discuss future research directions such as the use of ectopic plastid differentiation as a valuable tool to characterize factors involved in plastid differentiation. Collectively, we highlight important advances in the field that can also impact future agricultural and biotechnological improvement in plants.

Understanding protein import in diverse non-green plastids

The spectacular diversity of plastids in non-green organs such as flowers, fruits, roots, tubers, and senescing leaves represents a Universe of metabolic processes in higher plants that remain to be completely characterized. The endosymbiosis of the plastid and the subsequent export of the ancestral cyanobacterial genome to the nuclear genome, and adaptation of the plants to all types of environments has resulted in the emergence of diverse and a highly orchestrated metabolism across the plant kingdom that is entirely reliant on a complex protein import and translocation system. The TOC and TIC translocons, critical for importing nuclear-encoded proteins into the plastid stroma, remain poorly resolved, especially in the case of TIC. From the stroma, three core pathways (cpTat, cpSec, and cpSRP) may localize imported proteins to the thylakoid. Non-canonical routes only utilizing TOC also exist for the insertion of many inner and outer membrane proteins, or in the case of some modified proteins, a vesicular import route. Understanding this complex protein import system is further compounded by the highly heterogeneous nature of transit peptides, and the varying transit peptide specificity of plastids depending on species and the developmental and trophic stage of the plant organs. Computational tools provide an increasingly sophisticated means of predicting protein import into highly diverse non-green plastids across higher plants, which need to be validated using proteomics and metabolic approaches. The myriad plastid functions enable higher plants to interact and respond to all kinds of environments. Unraveling the diversity of non-green plastid functions across the higher plants has the potential to provide knowledge that will help in developing climate resilient crops.

Multiple domestication events explain the origin of Gossypium hirsutum landraces in Mexico

Several Mesoamerican crops constitute wild-to-domesticated complexes generated by multiple initial domestication events, and continuous gene flow among crop populations and between these populations and their wild relatives. It has been suggested that the domestication of cotton (Gossypium hirsutum) started in the northwest of the Yucatán Peninsula, from where it spread to other regions inside and outside of Mexico. We tested this hypothesis by assembling chloroplast genomes of 23 wild, landraces, and breeding lines (transgene-introgressed and conventional). The phylogenetic analysis showed that the evolutionary history of cotton in Mexico involves multiple events of introgression and genetic divergence. From this, we conclude that Mexican landraces arose from multiple wild populations. Our results also revealed that their structural and functional chloroplast organizations had been preserved. However, genetic diversity decreases as a consequence of domestication, mainly in transgene-introgressed (TI) individuals (π = 0.00020, 0.00001, 0.00016, 0, and 0, of wild, TI-wild, landraces, TI-landraces, and breeding lines, respectively). We identified homologous regions that differentiate wild from domesticated plants and indicate a relationship among the samples. A decrease in genetic diversity associated with transgene introgression in cotton was identified for the first time, and our outcomes are therefore relevant to both biosecurity and agrobiodiversity conservation.

Architecture of chloroplast TOC-TIC translocon supercomplex

Chloroplasts rely on the translocon complexes in the outer and inner envelope membranes (the TOC and TIC complexes, respectively) to import thousands of different nuclear-encoded proteins from the cytosol^1-4. Although previous studies indicated that the TOC and TIC complexes may assemble into larger supercomplexes^5-7, the overall architectures of the TOC-TIC supercomplexes and the mechanism of preprotein translocation are unclear. Here we report the cryo-electron microscopy structure of the TOC-TIC supercomplex from Chlamydomonas reinhardtii. The major subunits of the TOC complex (Toc75, Toc90 and Toc34) and TIC complex (Tic214, Tic20, Tic100 and Tic56), three chloroplast translocon-associated proteins (Ctap3, Ctap4 and Ctap5) and three newly identified small inner-membrane proteins (Simp1-3) have been located in the supercomplex. As the largest protein, Tic214 traverses the inner membrane, the intermembrane space and the outer membrane, connecting the TOC complex with the TIC proteins. An inositol hexaphosphate molecule is located at the Tic214-Toc90 interface and stabilizes their assembly. Four lipid molecules are located within or above an inner-membrane funnel formed by Tic214, Tic20, Simp1 and Ctap5. Multiple potential pathways found in the TOC-TIC supercomplex may support translocation of different substrate preproteins into chloroplasts.

Regulation of alternative splicing by retrograde and light signals converges to control chloroplast proteins

Retrograde signals sent by chloroplasts control transcription in the nucleus. These signals antagonistically converge with light signals to coordinate the expression of genes involved in chloroplast functioning and seedling development. Although significant advances have been made in understanding the molecular interplay between light and retrograde signals at the transcriptional level, little is known about their interconnection at the post-transcriptional level. By using different publicly available datasets, this study addresses the influence of retrograde signaling on alternative splicing and defines the molecular and biological functions of this regulation. These analyses revealed that alternative splicing mimics transcriptional responses triggered by retrograde signals at different levels. First, both molecular processes similarly depend on the chloroplast-localized pentatricopeptide-repeat protein GUN1 to modulate the nuclear transcriptome. Secondly, as described for transcriptional regulation, alternative splicing coupled with the nonsense-mediated decay pathway effectively downregulates expression of chloroplast proteins in response to retrograde signals. Finally, light signals were found to antagonistically control retrograde signaling-regulated splicing isoforms, which consequently generates opposite splicing outcomes that likely contribute to the opposite roles these signals play in controlling chloroplast functioning and seedling development.

Understanding protein translocation across chloroplast membranes: Translocons and motor proteins

Subcellular organelles in eukaryotes are surrounded by lipid membranes. In an endomembrane system, vesicle trafficking is the primary mechanism for the delivery of organellar proteins to specific organelles. However, organellar proteins for chloroplasts, mitochondria, the nucleus, and peroxisomes that are translated in the cytosol are directly imported into their target organelles. Chloroplasts are a plant-specific organelle with outer and inner envelope membranes, a dual-membrane structure that is similar to mitochondria. Interior chloroplast proteins translated by cytosolic ribosomes are thus translocated through TOC and TIC complexes (translocons in the outer and inner envelope of chloroplasts, respectively), with stromal ATPase motor proteins playing a critical role in pulling pre-proteins through these import channels. Over the last three decades, the identity and function of TOC/TIC components and stromal motor proteins have been actively investigated, which has shed light on the action mechanisms at a molecular level. However, there remains some disagreement over the exact composition of TIC complexes and genuine stromal motor proteins. In this review, we discuss recent findings on the mechanisms by which proteins are translocated through TOC/TIC complexes and discuss future prospects for this field of research.

Genome-wide association analysis for drought tolerance and associated traits in faba bean (Vicia faba L.)

Faba bean (Vicia faba L.) is an important high protein legume adapted to diverse climatic conditions with multiple benefits for the overall sustainability of the cropping systems. Plant-based protein demand is being expanded and faba bean is a good candidate to cover this need. However, the crop is very sensitive to abiotic stresses, especially drought, which severely affects faba bean yield and development worldwide. Therefore, identifying genes associated with drought stress tolerance is a major challenge in faba bean breeding. Although the faba bean response to drought stress has been widely studied, the molecular approaches to improve drought tolerance in this crop are still limited. Here we built on recent genomic advances such as the development of the first high-density SNP genotyping array, to conduct a genome-wide association study (GWAS) using thousands of genetic polymorphisms throughout the entire faba bean genome. A worldwide collection of 100 faba bean accessions was grown under control and drought conditions and 10 morphological, phenological and physiological traits were evaluated to identify single nucleotide polymorphism (SNP) markers associated with drought tolerance. We identified 29 SNP markers significantly correlated with these traits under drought stress conditions. The flanking sequences were blasted to the Medicago truncatula reference genomes in order to annotate potential candidate genes underlying the causal variants. Three of the SNPs for chlorophyll content after the stress, correspond to uncharacterized proteins indicating the presence of novel genes associated with drought tolerance in faba bean. The significance of stress-inducible signal transducers provides valuable information on the possible mechanisms underlying the faba bean response to drought stress, thus providing a foundation for future marker-assisted breeding in the crop.

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.

Structure of a TOC-TIC supercomplex spanning two chloroplast envelope membranes

The TOC and TIC complexes are essential translocons that facilitate the import of the nuclear genome-encoded preproteins across the two envelope membranes of chloroplast, but their exact molecular identities and assembly remain unclear. Here, we report a cryoelectron microscopy structure of TOC-TIC supercomplex from Chlamydomonas, containing a total of 14 identified components. The preprotein-conducting pore of TOC is a hybrid β-barrel co-assembled by Toc120 and Toc75, while the potential translocation path of TIC is formed by transmembrane helices from Tic20 and YlmG, rather than a classic model of Tic110. A rigid intermembrane space (IMS) scaffold bridges two chloroplast membranes, and a large hydrophilic cleft on the IMS scaffold connects TOC and TIC, forming a pathway for preprotein translocation. Our study provides structural insights into the TOC-TIC supercomplex composition, assembly, and preprotein translocation mechanism, and lays a foundation to interpret the evolutionary conservation and diversity of this fundamental translocon machinery.

Ten Plastomes of Crassula (Crassulaceae) and Phylogenetic Implications

The genus Crassula is the second-largest genus in the family Crassulaceae, with about 200 species. As an acknowledged super-barcode, plastomes have been extensively utilized for plant evolutionary studies. Here, we first report 10 new plastomes of Crassula. We further focused on the structural characterizations, codon usage, aversion patterns, and evolutionary rates of plastomes. The IR junction patterns-IRb had 110 bp expansion to rps19-were conservative among Crassula species. Interestingly, we found the codon usage patterns of matK gene in Crassula species are unique among Crassulaceae species with elevated ENC values. Furthermore, subgenus Crassula species have specific GC-biases in the matK gene. In addition, the codon aversion motifs from matK, pafI, and rpl22 contained phylogenetic implications within Crassula. The evolutionary rates analyses indicated all plastid genes of Crassulaceae were under the purifying selection. Among plastid genes, ycf1 and ycf2 were the most rapidly evolving genes, whereas psaC was the most conserved gene. Additionally, our phylogenetic analyses strongly supported that Crassula is sister to all other Crassulaceae species. Our findings will be useful for further evolutionary studies within the Crassula and Crassulaceae.

Presequence translocase-associated motor subunits of the mitochondrial protein import apparatus are dual-targeted to mitochondria and plastids

The import and assembly of most of the mitochondrial proteome is regulated by protein translocases located within the mitochondrial membranes. The Presequence Translocase-Associated Motor (PAM) complex powers the translocation of proteins across the inner membrane and consists of Hsp70, the J-domain containing co-chaperones, Pam16 and Pam18, and their associated proteins Tim15 and Mge1. In Arabidopsis, multiple orthologues of Pam16, Pam18, Tim15 and Mge1 have been identified and a mitochondrial localization has been confirmed for most. As the localization of Pam18-1 has yet to be determined and a plastid localization has been observed for homologues of Tim15 and Mge1, we carried out a comprehensive targeting analysis of all PAM complex orthologues using multiple in vitro and in vivo methods. We found that, Pam16 was exclusively targeted to the mitochondria, but Pam18 orthologues could be targeted to both the mitochondria and plastids, as observed for the PAM complex interacting partner proteins Tim15 and Mge1.

Chloroplast protein import machinery and quality control

Most chloroplast proteins are nucleus-encoded, translated on cytoplasmic ribosomes as precursor proteins, and imported into chloroplasts through TOC and TIC, the translocons of the outer and inner chloroplast envelope membranes. While the composition of the TOC complex is well established, there is still some controversy about the importance of a recently identified TIC complex consisting of Tic20, Tic214, Tic100, and Tic56. TOC and TIC form a supercomplex with a protein channel at the junction of the outer and inner envelope membranes through which preproteins are pulled into the stroma by the ATP-powered Ycf2 complex consisting of several FtsH-like ATPases and/or by chloroplast Hsp proteins. Several components of the TOC/TIC system are moonlighting proteins with additional roles in chloroplast gene expression and metabolism. Chaperones and co-chaperones, associated with TOC and TIC on the cytoplasmic and stromal side of the chloroplast envelope, participate in the unfolding and folding of the precursor proteins and act together with the ubiquitin-proteasome system in protein quality control. Chloroplast protein import is also intimately linked with retrograde signaling, revealing altogether an unsuspected complexity in the regulation of this process.

Tic12, a 12-kDa essential component of the translocon at the inner envelope membrane of chloroplasts in Arabidopsis

The complexes translocon at the outer envelope membrane of chloroplasts and translocon at the inner envelope membrane of chloroplasts (TIC) mediate preprotein translocation across the chloroplast outer and inner envelope membranes, respectively. Tic20, Tic56, Tic100, and Tic214 form a stable one-megadalton TIC whose function is essential for Arabidopsis thaliana and Chlamydomonas reinhardtii. Tic20 plays a central role in preprotein translocation by forming a protein-conducting channel. Tic56, Tic100, and Tic214 are also indispensable for TIC function, but whether other components are required for this process remain unclear. Here, we demonstrate that a 12-kDa protein named Tic12 is part of the TIC in A. thaliana and participates in preprotein translocation across the inner envelope membrane. Tic12 was tightly associated with the TIC but disassociated under high-salt conditions in combination with Triton X-100. Site-specific UV crosslinking experiments revealed that Tic12 and Tic20 directly interact with the transit peptide of a translocating preprotein. The tic12 null mutants are albino and seedling lethal, similar to the other tic null mutants. Tic12 and Tic20 were also involved in preprotein translocation in (Pisum sativum) pea chloroplasts. Thus, Tic12 is an essential constituent that forms the functional core together with Tic20 in the one-megadalton TIC.

Complete chloroplast genomes of two medicinal Swertia species: the comparative evolutionary analysis of Swertia genus in the Gentianaceae family

The complete chloroplast genome of Swertia kouitchensis has been sequenced and assembled, compared with that of S. bimaculata to determine the evolutionary relationships among species of the Swertia in the Gentianaceae family. Swertia kouitchensis and S. bimaculata are from the Gentianaceae family. The complete chloroplast genome of S. kouitchensis was newly assembled, annotated, and analyzed by Illumina Hiseq 2500 platform. The chloroplast genomes of the two species encoded a total of 133, 134 genes, which included 88-89 protein-coding genes, 37 transfer RNA (tRNA) genes, and 8 ribosomal RNA genes. One intron was contained in each of the eight protein-coding genes and eight tRNA-coding genes, whereas two introns were found in two genes (ycf3 and clpP). The most abundant codon of the two species was for isoleucine, and the least abundant codon was for cysteine. The number of microsatellite repeat sequences was twenty-eight and thirty-two identified in the chloroplast genomes of S. kouitchensis and S. bimaculata, respectively. A total of 1127 repeat sequences were identified in all the 23 Swertia chloroplast genomes, and they fell into four categories. Furthermore, five divergence hotspot regions can be applied to discriminate these 23 Swertia species through genomes comparison. One pair of genus-specific DNA barcodes primer has been accurately identified. Therefore, the diverse regions cloned by a specific primer may become an effective and powerful molecular marker for the identification of Swertia genus. Moreover, four genes (ccsA, ndhK, rpoC1, and rps12) were positive selective pressure. The phylogenetic tree showed that the 23 Swertia species were clustered into a large clade including four evident subbranches, whereas the two species of S. kouitchensis and S. bimaculata were separately clustered into the diverse but correlated species group.

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.

A distinct class of GTP-binding proteins mediates chloroplast protein import in Rhodophyta

Chloroplast protein import is mediated by translocons named TOC and TIC on the outer and inner envelope membranes, respectively. Translocon constituents are conserved among green lineages, including plants and green algae. However, it remains unclear whether Rhodophyta (red algae) share common chloroplast protein import mechanisms with the green lineages. We show that in the rhodophyte Cyanidioschyzon merolae, plastome-encoded Tic20_pt localized to the chloroplast envelope and was transiently associated with preproteins during import, suggesting its conserved function as a TIC constituent. Besides plastome-encoded FtsH_pt and several chaperones, a class of GTP (guanosine 5′-triphosphate)-binding proteins distinct from the Toc34/159 GTPase family associated transiently with preproteins. This class of proteins resides mainly in the cytosol and shows sequence similarities with Sey1/RHD3, required for endoplasmic reticulum membrane fusion, and with the periplastid-localized import factor PPP1, previously identified in the Apicomplexa and diatoms. These GTP-binding proteins, named plastid targeting factor for protein import 1 (PTF1) to PTF3, may act as plastid targeting factors in Rhodophyta.

Mutations in the chloroplast inner envelope protein TIC100 impair and repair chloroplast protein import and impact retrograde signaling

Chloroplast biogenesis requires synthesis of proteins in the nucleocytoplasm and the chloroplast itself. Nucleus-encoded chloroplast proteins are imported via multiprotein translocons in the organelle's envelope membranes. Controversy exists around whether a 1-MDa complex comprising TIC20, TIC100, and other proteins constitutes the inner membrane TIC translocon. The Arabidopsis thaliana cue8 virescent mutant is broadly defective in plastid development. We identify CUE8 as TIC100. The tic100cue8 mutant accumulates reduced levels of 1-MDa complex components and exhibits reduced import of two nucleus-encoded chloroplast proteins of different import profiles. A search for suppressors of tic100cue8 identified a second mutation within the same gene, tic100soh1, which rescues the visible, 1 MDa complex-subunit abundance, and chloroplast protein import phenotypes. tic100soh1 retains but rapidly exits virescence and rescues the synthetic lethality of tic100cue8 when retrograde signaling is impaired by a mutation in the GENOMES UNCOUPLED 1 gene. Alongside the strong virescence, changes in RNA editing and the presence of unimported precursor proteins show that a strong signaling response is triggered when TIC100 function is altered. Our results are consistent with a role for TIC100, and by extension the 1-MDa complex, in the chloroplast import of photosynthetic and nonphotosynthetic proteins, a process which initiates retrograde signaling.

Organellar transcripts dominate the cellular mRNA pool across plants of varying ploidy levels

Mitochondrial and plastid functions depend on coordinated expression of proteins encoded by genomic compartments that have radical differences in copy number of organellar and nuclear genomes. In polyploids, doubling of the nuclear genome may add challenges to maintaining balanced expression of proteins involved in cytonuclear interactions. Here, we use ribo-depleted RNA sequencing (RNA-seq) to analyze transcript abundance for nuclear and organellar genomes in leaf tissue from four different polyploid angiosperms and their close diploid relatives. We find that even though plastid genomes contain <1% of the number of genes in the nuclear genome, they generate the majority (69.9 to 82.3%) of messenger RNA (mRNA) transcripts in the cell. Mitochondrial genes are responsible for a much smaller percentage (1.3 to 3.7%) of the leaf mRNA pool but still produce much higher transcript abundances per gene compared to nuclear genome. Nuclear genes encoding proteins that functionally interact with mitochondrial or plastid gene products exhibit mRNA expression levels that are consistently more than 10-fold lower than their organellar counterparts, indicating an extreme cytonuclear imbalance at the RNA level despite the predominance of equimolar interactions at the protein level. Nevertheless, interacting nuclear and organellar genes show strongly correlated transcript abundances across functional categories, suggesting that the observed mRNA stoichiometric imbalance does not preclude coordination of cytonuclear expression. Finally, we show that nuclear genome doubling does not alter the cytonuclear expression ratios observed in diploid relatives in consistent or systematic ways, indicating that successful polyploid plants are able to compensate for cytonuclear perturbations associated with nuclear genome doubling.

Comparative Analysis of the Complete Chloroplast Genomes in Allium Section Bromatorrhiza Species (Amaryllidaceae): Phylogenetic Relationship and Adaptive Evolution

With the development of molecular sequencing approaches, many taxonomic and phylogenetic problems of the genus Allium L. have been solved; however, the phylogenetic relationships of some subgenera or sections, such as section Bromatorrhiza, remain unresolved, which has greatly impeded our full understanding of the species relationships among the major clades of Allium. In this study, the complete chloroplast (cp) genomes of nine species in the Allium sect. Bromatorrhiza were determined using the Illumina paired-end sequencing, the NOVOPlasty de novo assembly strategy, and the PGA annotation method. The results showed that the cp genome exhibited high conservation and revealed a typical circular tetrad structure. Among the sect. Bromatorrhiza species, the gene content, SSRs, codon usage, and RNA editing site were similar. The genome structure and IR regions’ fluctuation were investigated while genes, CDSs, and non-coding regions were extracted for phylogeny reconstruction. Evolutionary rates (Ka/Ks values) were calculated, and positive selection analysis was further performed using the branch-site model. Five hypervariable regions were identified as candidate molecular markers for species authentication. A clear relationship among the sect. Bromatorrhiza species were detected based on concatenated genes and CDSs, respectively, which suggested that sect. Bromatorrhiza is monophyly. In addition, there were three genes with higher Ka/Ks values (rps2, ycf1, and ycf2), and four genes (rpoC2, atpF, atpI, and rpl14) were further revealed to own positive selected sites. These results provide new insights into the plastome component, phylogeny, and evolution of Allium species.

Characterization of the complete chloroplast genome sequences of six Dalbergia species and its comparative analysis in the subfamily of Papilionoideae (Fabaceae)

Dalbergia spp. are numerous and widely distributed in pantropical areas in Asia, Africa and America, and most of the species have important economic and ecological value as precious timber. In this study, we determined and characterized six complete chloroplast genomes of Dalbergia species (Dalbergia obtusifolia, D. hupeana, D. mimosoides, D. sissoo, D. hancei, D. balansae), which displayed the typical quadripartite structure of angiosperms. The sizes of the genomes ranged from 155,698 bp (D. hancei) to 156,419 bp (D. obtusifolia). The complete chloroplast genomes of Dalbergia include 37 tRNA genes, eight rRNA genes and 84 protein-coding genes. We analysed the sequence diversity of Dalberigia chloroplast genomes coupled with previous reports. The results showed 12 noncoding regions (rps16-accD, trnR-UCU-trnG-UCC, ndhE-ndhG, trnG-UCC-psbZ, rps8-rpl14, trnP-UGG-psaJ, ndhH-rps15, trnQ-UUG-rps16, trnS-GCU-psbI, rps12-clpP, psbA-trnK-UUU, trnK-UUU-intron), and four coding regions (rps16, ycf1, rps15 and ndhF) showed many nucleotide variations that could be used as potential molecular markers. Based on a site-specific model, we analysed the selective pressure of chloroplast genes in Dalbergia species. Twenty-two genes with positively selected sites were detected, involving the photosynthetic system (ndhC, adhD, ndhF, petB, psaA, psaB, psbB, psbC, psbK and rbcL), self-replication category of genes (rpoA, rpoC2, rps3, rps12 and rps18) and others (accD, ccsA, cemA, clpP, matK, ycf1 and ycf2). Additionally, we identified potential RNA editing sites that were relatively conserved in the genus Dalbergia. Furthermore, the comparative analysis of cp genomes of Dalbergieae species indicated that the boundary of IRs/SSC was highly variable, which resulted in the size variation of cp genomes. Finally, phylogenetic analysis showed an inferred phylogenetic tree of Papilionoideae species with high bootstrap support and suggested that Amorpheae was the sister of the clade Dalbergieae. Moreover, three genera of the Pterocarpus clade showed a nested evolutionary relationship. These complete cp genomes provided valuable information for understanding the genetic variation and phylogenetic relationship of Dalbergia species with their relatives.

Comparative proteomic analysis on chloroplast proteins provides new insights into the effects of low temperature in sugar beet

BACKGROUND

Low temperature, which is one of the main environmental factors that limits geographical distribution and sucrose yield, is a common abiotic stress during the growth and development of sugar beet. As a regulatory hub of plant response to abiotic stress, activity in the chloroplasts is related to many molecular and physiological processes, particularly in response to low temperature stress.

RESULTS

The contents of chlorophyll (Chl) and malondialdehyde (MDA), relative electrical conductivity (REL), and superoxide dismutase (SOD) activity were measured. The results showed that sugar beet could manage low temperature stress by regulating the levels of Chl, REL and MDA, and the activity of SOD. The physiological responses indicated that sugar beets respond positively to low temperature treatments and are not significantly damaged. Moreover, to determine the precise time to response low temperature in sugar beet, well-known abiotic stresses-responsive transcript factor family, namely DEHYDRATION RESPONSIVE ELEMENT BINDING PROTEIN (DREB), was selected as the marker gene. The results of phylogenetic analyses showed that BvDREBA1 and BvDREBA4 were in the same branch as the cold- and drought-responsive AtDREB gene. In addition, the expression of BvDREBs reached its maximum level at 24 h after low temperature by RNA-Seq and qRT-PCR analysis. Furthermore, the changes in chloroplast proteome after low temperature at 24 h were detected using a label-free technique. A total of 416 differentially expressed proteins were identified. GO enrichment analysis showed that 16 GO terms were significantly enriched, particularly chloroplast stroma, chloroplast envelope, and chloroplast thylakoid membrane. It is notable that the transport of photosynthetic proteins (BvLTD and BvTOC100), the formation of starch granules (BvPU1, BvISA3, and BvGWD3) and the scavenging of reactive oxygen species (BvCu/Zn-SOD, BvCAT, BvPrx, and BvTrx) were the pathways used by sugar beets to respond to low temperatures at an early stage.

CONCLUSIONS

These results provide a preliminarily analysis of how chloroplasts of sugar beet respond to low temperature stress at the translational level and provide a theoretical basis for breeding low temperature resistant varieties of sugar beet.

Genetic Mapping of the Gmpgl3 Mutant Reveals the Function of GmTic110a in Soybean Chloroplast Development

The generation of oxygen and organic matter in plants mainly depends on photosynthesis, which directly affects plant growth and development. The chloroplast is the main organelle in which photosynthesis occurs. In this study, a Glycine max pale green leaf 3-1 (Gmpgl3-1) mutant was isolated from the soybean mutagenized population. The Gmpgl3-1 mutant presented with decreased chlorophyll contents, reduced chloroplast stroma thylakoids, reduced yields, and decreased numbers of pods per plant. Bulked segregant analysis (BSA) together with map-based cloning revealed a single-nucleotide non-synonymous mutation at the 341st nucleotide of the first exon of the chloroplast development-related GmTic110a gene. The phenotype of the knockout plants was the same as that of the mutant. The GmTic110a gene was highly expressed in the leaves at various developmental stages, and its protein was localized to the inner chloroplast membrane. Split luciferase complementation assays and coimmunoprecipitation (co-IP) experiments revealed that GmTic110a interacted with GmTic20, GmTic40a, and GmTic40b in tobacco leaves. These results indicated that the GmTic110a gene plays an important role in chloroplast development.

Genetic analysis of chlorophyll synthesis and degradation regulated by BALANCE of CHLOROPHYLL METABOLISM

Chlorophyll (Chl) serves a number of essential functions, capturing and converting light energy as a component of photosystem supercomplexes. Chl degradation during leaf senescence is also required for adequate degeneration of chloroplasts and salvaging of nutrients from senescent leaves. In this study, we performed genetic analysis to determine the functions of BALANCE of CHLOROPHYLL METABOLISM1 (BCM1) and BCM2, which control Chl levels by regulating synthesis and degradation, and STAY-GREEN (SGR)1 (also known as NON-YELLOWING1 [NYE1]) and SGR2, which encode Mg-dechelatase and catalyze Chl a degradation in Arabidopsis (Arabidopsis thaliana). Analysis of bcm1 bcm2 revealed that both BCM1 and BCM2 are involved in the regulation of Chl levels in presenescent leaves and Chl degradation in senescing leaves. Analysis of bcm1 bcm2 nye1 nye2 suggested that BCMs repress Chl-degrading activity in both presenescent and senescing leaves by regulating SGR activity. Furthermore, transactivation analysis and chromatin immunoprecipitation (ChIP) assay revealed that GOLDEN2-LIKE1 (GLK1), a central transcription factor regulating the expression of genes encoding photosystem-related proteins, such as light-harvesting Chl a/b-binding proteins (LHCPs), directly regulates the transcription of BCM1. LHCPs are stabilized by Chl binding, suggesting that GLKs control the amount of LHCP through transcriptional and post-translational regulation via BCM-mediated Chl-level regulation. Meanwhile, we generated a mutant of the BCM ortholog in lettuce (Lactuca sativa) by genome editing and found that it showed an early yellowing phenotype, but only a slight reduction in Chl in presenescent leaves. Thus, this study revealed a conserved but slightly diversified regulation of Chl and LHCP levels via the GLK-BCM pathway in eudicots.

Molecular evolution and phylogenetic relationships of Ligusticum (Apiaceae) inferred from the whole plastome sequences

Background

The genus Ligusticum belongs to Apiaceae, and its taxonomy has long been a major difficulty. A robust phylogenetic tree is the basis of accurate taxonomic classification of Ligusticum. We herein used 26 (including 14 newly sequenced) plastome-scale data to generate reliable phylogenetic trees to explore the phylogenetic relationships of Chinese Ligusticum.

Results

We found that these plastid genomes exhibited diverse plastome characteristics across all four currently identified clades in China, while the plastid protein-coding genes were conserved. The phylogenetic analyses by the concatenation and coalescent methods obtained a more robust molecular phylogeny than prior studies and showed the non-monophyly of Chinese Ligusticum. In the concatenation-based phylogeny analyses, the two datasets yielded slightly different topologies that may be primarily due to the discrepancy in the number of variable sites.

Conclusions

Our plastid phylogenomics analyses emphasized that the current circumscription of the Chinese Ligusticum should be reduced, and the taxonomy of Ligusticum urgently needs revision. Wider taxon sampling including the related species of Ligusticum will be necessary to explore the phylogenetic relationships of this genus. Overall, our study provided new insights into the taxonomic classification of Ligusticum and would serve as a framework for future studies on taxonomy and delimitation of Ligusticum from the perspective of the plastid genome.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02010-z.

Complete Chloroplast Genome of Cnidium monnieri (Apiaceae) and Comparisons with Other Tribe Selineae Species

Cnidium monnieri is an economically important traditional Chinese medicinal plant. In this study, the complete chloroplast (cp) genome of C. monnieri was determined using the Illumina paired-end sequencing, the GetOrganelle de novo assembly strategy, as well as the GeSeq annotation method. Our results showed that the cp genome was 147,371 bp in length with 37.4% GC content and included a large single-copy region (94,361 bp) and a small single-copy region (17,552 bp) separated by a pair of inverted repeat regions (17,729 bp). A total of 129 genes were contained in the cp genome, including 85 protein-coding genes, 36 tRNA genes, and eight rRNA genes. We also investigated codon usage, RNA editing, repeat sequences, simple sequence repeats (SSRs), IR boundaries, and pairwise Ka/Ks ratios. Four hypervariable regions (trnD-trnY-trnE-trnT, ycf2, ndhF-rpl32-trnL, and ycf1) were identified as candidate molecular markers for species authentication. The phylogenetic analyses supported non-monophyly of Cnidium and C. monnieri located in tribe Selineae based on the cp genome sequences and internal transcribed spacer (ITS) sequences. The incongruence of the phylogenetic position of C. monnieri between ITS and cpDNA phylogenies suggested that C. monnieri might have experienced complex evolutions with hybrid and incomplete lineage sorting. All in all, the results presented herein will provide plentiful chloroplast genomic resources for studies of the taxonomy, phylogeny, and species authentication of C. monnieri. Our study is also conducive to elucidating the phylogenetic relationships and taxonomic position of Cnidium.

Comparative plastome analysis of Musaceae and new insights into phylogenetic relationships

Background

Musaceae is an economically important family consisting of 70-80 species. Elucidation of the interspecific relationships of this family is essential for a more efficient conservation and utilization of genetic resources for banana improvement. However, the scarcity of herbarium specimens and quality molecular markers have limited our understanding of the phylogenetic relationships in wild species of Musaceae. Aiming at improving the phylogenetic resolution of Musaceae, we analyzed a comprehensive set of 49 plastomes for 48 species/subspecies representing all three genera of this family.

Results

Musaceae plastomes have a relatively well-conserved genomic size and gene content, with a full length ranging from 166,782 bp to 172,514 bp. Variations in the IR borders were found to show phylogenetic signals to a certain extent in Musa. Codon usage bias analysis showed different preferences for the same codon between species and three genera and a common preference for A/T-ending codons. Among the two genes detected under positive selection (dN/dS > 1), ycf2 was indicated under an intensive positive selection. The divergent hotspot analysis allowed the identification of four regions (ndhF-trnL, ndhF, matK-rps16, and accD) as specific DNA barcodes for Musaceae species.

Bayesian and maximum likelihood phylogenetic analyses using full plastome resulted in nearly identical tree topologies with highly supported relationships between species. The monospecies genus Musella is sister to Ensete, and the genus Musa was divided into two large clades, which corresponded well to the basic number of n = x = 11 and n = x =10/9/7, respectively. Four subclades were divided within the genus Musa. A dating analysis covering the whole Zingiberales indicated that the divergence of Musaceae family originated in the Palaeocene (59.19 Ma), and the genus Musa diverged into two clades in the Eocene (50.70 Ma) and then started to diversify from the late Oligocene (29.92 Ma) to the late Miocene. Two lineages (Rhodochlamys and Australimusa) radiated recently in the Pliocene /Pleistocene periods.

Conclusions

The plastome sequences performed well in resolving the phylogenetic relationships of Musaceae and generated new insights into its evolution. Plastome sequences provided valuable resources for population genetics and phylogenetics at lower taxon.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08454-3.

Chloroplast Protein Tic55 Involved in Dark-Induced Senescence through AtbHLH/AtWRKY-ANAC003 Controlling Pathway of Arabidopsis thaliana

The chloroplast comprises the outer and inner membranes that are composed of the translocon protein complexes Toc and Tic (translocon at the outer/inner envelope membrane of chloroplasts), respectively. Tic55, a chloroplast Tic protein member, was shown to be not vital for functional protein import in Arabidopsis from previous studies. Instead, Tic55 was revealed to be a dark-induced senescence-related protein in our earlier study. To explore whether Tic55 elicits other biological functions, a tic55-II knockout mutant (SALK_086048) was characterized under different stress treatments. Abiotic stress conditions, such as cold, heat, and high osmotic pressure, did not cause visible effects on tic55-II mutant plant, when compared to the wild type (WT). In contrast, senescence was induced in the individually darkened leaves (IDLs), resulting in the differential expression of the senescence-related genes PEROXISOME DEFECTIVE 1 (PED1), BLUE COPPER-BINDING PROTEIN (BCB), SENESCENCE 1 (SEN1), and RUBISCO SMALL SUBUNIT GENE 2B (RBCS2B). The absence of Tic55 in tic55-II knockout mutant inhibited expression of the senescence-related genes PED1, BCB, and SEN1 at different stages of dark adaptation, while causing stimulation of RBCS2B gene expression at an early stage of dark response. Finally, yeast one-hybrid assays located the ANAC003 promoter region with cis-acting elements are responsible for binding to the different AtbHLH proteins, thereby causing the transactivation of an HIS3 reporter gene. ANAC003 was shown previously as a senescence-related protein and its activation would lead to expression of senescence-associated genes (SAGs), resulting in plant senescence. Thus, we propose a hypothetical model in which three signaling pathways may be involved in controlling the expression of ANAC003, followed by expression of SAGs that in turn leads to leaf senescence in Arabidopsis by this study and previous data.

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.

Genome structure and diversity among Cynanchum wilfordii accessions

BACKGROUND

Cynanchum wilfordii (Cw) and Cynanchum auriculatum (Ca) have long been used in traditional medicine and as functional food in Korea and China, respectively. They have diverse medicinal functions, and many studies have been conducted, including pharmaceutical efficiency and metabolites. Especially, Cw is regarded as the most famous medicinal herb in Korea due to its menopausal symptoms relieving effect. Despite the high demand for Cw in the market, both species are cultivated using wild resources with rare genomic information.

RESULTS

We collected 160 Cw germplasm from local areas of Korea and analyzed their morphological diversity. Five Cw and one Ca of them, which were morphologically diverse, were sequenced, and nuclear ribosomal DNA (nrDNA) and complete plastid genome (plastome) sequences were assembled and annotated. We investigated the genomic characteristics of Cw as well as the genetic diversity of plastomes and nrDNA of Cw and Ca. The Cw haploid nuclear genome was approximately 178 Mbp. Karyotyping revealed the juxtaposition of 45S and 5S nrDNA on one of 11 chromosomes. Plastome sequences revealed 1226 interspecies polymorphisms and 11 Cw intraspecies polymorphisms. The 160 Cw accessions were grouped into 21 haplotypes based on seven plastome markers and into 108 haplotypes based on seven nuclear markers. Nuclear genotypes did not coincide with plastome haplotypes that reflect the frequent natural outcrossing events.

CONCLUSIONS

Cw germplasm had a huge morphological diversity, and their wide range of genetic diversity was revealed through the investigation with 14 molecular markers. The morphological and genomic diversity, chromosome structure, and genome size provide fundamental genomic information for breeding of undomesticated Cw plants.

Comparative chloroplast genome analyses of Paraboea (Gesneriaceae): Insights into adaptive evolution and phylogenetic analysis

Paraboea (Gesneriaceae) distributed in the karst areas of South and Southwest China and Southeast Asia, is an ideal genus to study the phylogeny and adaptive evolution of karst plants. In this study, the complete chloroplast genomes of twelve Paraboea species were sequenced and analyzed. Twelve chloroplast genomes ranged in size from 153166 to 154245 bp. Each chloroplast genome had a typical quartile structure, and relatively conserved type and number of gene components, including 131 genes which are composed of 87 protein coding genes, 36 transfer RNAs and 8 ribosomal RNAs. A total of 600 simple sequence repeats and 389 non-overlapped sequence repeats were obtained from the twelve Paraboea chloroplast genomes. We found ten divergent regions (trnH-GUG-psbA, trnM-CAU, trnC-GCA, atpF-atpH, ycf1, trnK-UUU-rps16, rps15, petL, trnS-GCU-trnR-UCU and psaJ-rpl33) among the 12 Paraboea species to be potential molecular markers. In the phylogenetic tree of 31 Gesneriaceae plants including twelve Paraboea species, all Paraboea species clustered in a clade and confirmed the monophyly of Paraboea. Nine genes with positive selection sites were detected, most of which were related to photosynthesis and protein synthesis, and might played crucial roles in the adaptability of Paraboea to diverse karst environments. These findings are valuable for further study of the phylogeny and karst adaptability of Gesneriaceae plants.

Different rates of pollen and seed gene flow cause branch-length and geographic cytonuclear discordance within Asian butternuts

Topological cytonuclear discordance is commonly observed in plant phylogenetic and phylogeographic studies, yet few studies have attempted to detect two other forms of cytonuclear discordance (branch length and geographical) and to uncover the causes of the discordance. We used the whole nuclear and chloroplast genome data from 80 individual Asian butternuts to reveal the pattern and processes of cytonuclear discordance. Our findings indicate that the chloroplast genome had substantially deeper divergence (branch-length discordance) and a steeper cline in the contact zone (geographic discordance) compared with the nuclear genome. After various hypothesis have been tested, the results suggest that incomplete lineage sorting, positive selection and cytonuclear incompatibility are probably insufficient to explain this pattern. However, isolation-by-distance analysis and gene flow estimation point to a much higher level of gene flow by pollen compared with by seeds, which may have slowed down lineage divergence and mediated wider contact for nuclear genome compared with the chloroplast genome. Altogether, this study highlights a critical role of sex-biased dispersal in causing discordance between the nuclear and plastid genome of Asian butternuts. Given its ubiquity among plants, asymmetric gene flow should be given a high priority in future studies of cytonuclear discordance.

Extensive genomic rearrangements mediated by repetitive sequences in plastomes of Medicago and its relatives

BACKGROUND

Although plastomes are highly conserved with respect to gene content and order in most photosynthetic angiosperms, extensive genomic rearrangements have been reported in Fabaceae, particularly within the inverted repeat lacking clade (IRLC) of Papilionoideae. Two hypotheses, i.e., the absence of the IR and the increased repeat content, have been proposed to affect the stability of plastomes. However, this is still unclear for the IRLC species. Here, we aimed to investigate the relationships between repeat content and the degree of genomic rearrangements in plastomes of Medicago and its relatives Trigonella and Melilotus, which are nested firmly within the IRLC.

RESULTS

We detected abundant repetitive elements and extensive genomic rearrangements in the 75 newly assembled plastomes of 20 species, including gene loss, intron loss and gain, pseudogenization, tRNA duplication, inversion, and a second independent IR gain (IR ~ 15 kb in Melilotus dentata) in addition to the previous first reported cases in Medicago minima. We also conducted comparative genomic analysis to evaluate plastome evolution. Our results indicated that the overall repeat content is positively correlated with the degree of genomic rearrangements. Some of the genomic rearrangements were found to be directly linked with repetitive sequences. Tandem repeated sequences have been detected in the three genes with accelerated substitution rates (i.e., accD, clpP, and ycf1) and their length variation could be explained by the insertions of tandem repeats. The repeat contents of the three localized hypermutation regions around these three genes with accelerated substitution rates are also significantly higher than that of the remaining plastome sequences.

CONCLUSIONS

Our results suggest that IR reemergence in the IRLC species does not ensure their plastome stability. Instead, repeat-mediated illegitimate recombination is the major mechanism leading to genome instability, a pattern in agreement with recent findings in other angiosperm lineages. The plastome data generated herein provide valuable genomic resources for further investigating the plastome evolution in legumes.

Highly Diverse Shrub Willows (Salix L.) Share Highly Similar Plastomes

Plastome phylogenomics is used in a broad range of studies where single markers do not bear enough information. Phylogenetic reconstruction in the genus Salix is difficult due to the lack of informative characters and reticulate evolution. Here, we use a genome skimming approach to reconstruct 41 complete plastomes of 32 Eurasian and North American Salix species representing different lineages, different ploidy levels, and separate geographic regions. We combined our plastomes with published data from Genbank to build a comprehensive phylogeny of 61 samples (50 species) using RAxML (Randomized Axelerated Maximum Likelihood). Additionally, haplotype networks for two observed subclades were calculated, and 72 genes were tested to be under selection. The results revealed a highly conserved structure of the observed plastomes. Within the genus, we observed a variation of 1.68%, most of which separated subg. Salix from the subgeneric Chamaetia/Vetrix clade. Our data generally confirm previous plastid phylogenies, however, within Chamaetia/Vetrix phylogenetic results represented neither taxonomical classifications nor geographical regions. Non-coding DNA regions were responsible for most of the observed variation within subclades and 5.6% of the analyzed genes showed signals of diversifying selection. A comparison of nuclear restriction site associated DNA (RAD) sequencing and plastome data on a subset of 10 species showed discrepancies in topology and resolution. We assume that a combination of (i) a very low mutation rate due to efficient mechanisms preventing mutagenesis, (ii) reticulate evolution, including ancient and ongoing hybridization, and (iii) homoplasy has shaped plastome evolution in willows.

Rampant Nuclear Transfer and Substitutions of Plastid Genes in Passiflora

Gene losses in plastid genomes (plastomes) are often accompanied by functional transfer to the nucleus or substitution of an alternative nuclear-encoded gene. Despite the highly conserved gene content in plastomes of photosynthetic land plants, recent gene loss events have been documented in several disparate angiosperm clades. Among these lineages, Passiflora lacks several essential ribosomal genes, rps7, rps16, rpl20, rpl22, and rpl32, the two largest plastid genes, ycf1 and ycf2, and has a highly divergent rpoA. Comparative transcriptome analyses were performed to determine the fate of the missing genes in Passiflora. Putative functional transfers of rps7, rpl22, and rpl32 to nucleus were detected, with the nuclear transfer of rps7, representing a novel event in angiosperms. Plastid-encoded rps7 was transferred into the intron of a nuclear-encoded plastid-targeted thioredoxin m-type gene, acquiring its plastid transit peptide (TP). Plastid rpl20 likely experienced a novel substitution by a duplicated, nuclear-encoded mitochondrial-targeted rpl20 that has a similar gene structure. Additionally, among rosids, evidence for a third independent transfer of rpl22 in Passiflora was detected that gained a TP from a nuclear gene containing an organelle RNA recognition motif. Nuclear transcripts representing rpoA, ycf1, and ycf2 were not detected. Further analyses suggest that the divergent rpoA remains functional and that the gene is under positive or purifying selection in different clades. Comparative analyses indicate that alternative translocon and motor protein complexes may have substituted for the loss of ycf1 and ycf2 in Passiflora.

The FtsHi Enzymes of Arabidopsis thaliana: Pseudo-Proteases with an Important Function

FtsH metalloproteases found in eubacteria, animals, and plants are well-known for their vital role in the maintenance and proteolysis of membrane proteins. Their location is restricted to organelles of endosymbiotic origin, the chloroplasts, and mitochondria. In the model organism Arabidopsis thaliana, there are 17 membrane-bound FtsH proteases containing an AAA+ (ATPase associated with various cellular activities) and a Zn2+ metalloprotease domain. However, in five of those, the zinc-binding motif HEXXH is either mutated (FtsHi1, 2, 4, 5) or completely missing (FtsHi3), rendering these enzymes presumably inactive in proteolysis. Still, homozygous null mutants of the pseudo-proteases FtsHi1, 2, 4, 5 are embryo-lethal. Homozygous ftshi3 or a weak point mutant in FTSHi1 are affected in overall plant growth and development. This review will focus on the findings concerning the FtsHi pseudo-proteases and their involvement in protein import, leading to consequences in embryogenesis, seed growth, chloroplast, and leaf development and oxidative stress management.