Functioning of PPR Proteins in Organelle RNA Metabolism and Chloroplast Biogenesis

Genetic tapestry of Capsicum fruit colors: a comparative analysis of four cultivated species

KEY MESSAGE

Genome-wide association study of color spaces across the four cultivated Capsicum spp. revealed a shared set of genes influencing fruit color, suggesting mechanisms and pathways across Capsicum species are conserved during the speciation. Notably, Cytochrome P450 of the carotenoid pathway, MYB transcription factor, and pentatricopeptide repeat-containing protein are the major genes responsible for fruit color variation across the Capsicum species. Peppers (Capsicum spp.) rank among the most widely consumed spices globally. Fruit color, serving as a determinant for use in food colorants and cosmeceuticals and an indicator of nutritional contents, significantly influences market quality and price. Cultivated Capsicum species display extensive phenotypic diversity, especially in fruit coloration. Our study leveraged the genetic variance within four Capsicum species (Capsicum baccatum, Capsicum chinense, Capsicum frutescens, and Capsicum annuum) to elucidate the genetic mechanisms driving color variation in peppers and related Solanaceae species. We analyzed color metrics and chromatic attributes (Red, Green, Blue, L*, a*, b*, Luminosity, Hue, and Chroma) on samples cultivated over six years (2015-2021). We resolved genomic regions associated with fruit color diversity through the sets of SNPs obtained from Genotyping by Sequencing (GBS) and genome-wide association study (GWAS) with a Multi-Locus Mixed Linear Model (MLMM). Significant SNPs with FDR correction were identified, within the Cytochrome P450, MYB-related genes, Pentatricopeptide repeat proteins, and ABC transporter family were the most common among the four species, indicating comparative evolution of fruit colors. We further validated the role of a pentatricopeptide repeat-containing protein (Chr01:31,205,460) and a cytochrome P450 enzyme (Chr08:45,351,919) via competitive allele-specific PCR (KASP) genotyping. Our findings advance the understanding of the genetic underpinnings of Capsicum fruit coloration, with developed KASP assays holding potential for applications in crop breeding and aligning with consumer preferences. This study provides a cornerstone for future research into exploiting Capsicum's diverse fruit color variation.

Identification of Chromosomal Regions and Candidate Genes for Round leaf Locus in Cucumis melo L

Leaf morphology plays a crucial role in plant classification and provides a significant model for studying plant diversity while directly impacting photosynthetic efficiency. In the case of melons, leaf shape not only influences production and classification but also represents a key genetic trait that requires further exploration. In this study, we utilized forward genetics to pinpoint a recessive locus, dubbed Cmrl (Round leaf), which is responsible for regulating melon leaf shape. Through bulked segregant analysis sequencing and extensive evaluation of a two-year F2 population, we successfully mapped the Cmrl locus to a 537.07 kb region on chromosome 8 of the melon genome. Subsequent genetic fine-mapping efforts, leveraging a larger F2 population encompassing 1322 plants and incorporating F2:3 phenotypic data, further refined the locus to an 80.27 kb interval housing five candidate genes. Promoter analysis and coding sequence cloning confirmed that one of these candidates, MELO3C019152.2 (Cmppr encoding a pentatricopeptide repeat-containing family protein, Cmppr), stands out as a strong candidate gene for the Cmrl locus. Notably, comparisons of Cmrl expressions across various stages of leaf development and different leaf regions suggest a pivotal role of Cmrl in the morphogenesis of melon leaves.

OsPGL3A encodes a DYW-type pentatricopeptide repeat protein involved in chloroplast RNA processing and regulated chloroplast development

The chloroplast serves as the primary site of photosynthesis, and its development plays a crucial role in regulating plant growth and morphogenesis. The Pentatricopeptide Repeat Sequence (PPR) proteins constitute a vast protein family that function in the post-transcriptional modification of RNA within plant organelles. In this study, we characterized mutant of rice with pale green leaves (pgl3a). The chlorophyll content of pgl3a at the seedling stage was significantly reduced compared to the wild type (WT). Transmission electron microscopy (TEM) and quantitative PCR analysis revealed that pgl3a exhibited aberrant chloroplast development compared to the wild type (WT), accompanied by significant alterations in gene expression levels associated with chloroplast development and photosynthesis. The Mutmap analysis revealed that a single base deletionin the coding region of Os03g0136700 in pgl3a. By employing CRISPR/Cas9 mediated gene editing, two homozygous cr-pgl3a mutants were generated and exhibited a similar phenotype to pgl3a, thereby confirming that Os03g0136700 was responsible for pgl3a. Consequently, it was designated as OsPGL3A. OsPGL3A belongs to the DYW-type PPR protein family and is localized in chloroplasts. Furthermore, we demonstrated that the RNA editing efficiency of rps8-182 and rpoC2-4106, and the splicing efficiency of ycf3-1 were significantly decreased in pgl3a mutants compared to WT. Collectively, these results indicate that OsPGL3A plays a crucial role in chloroplast development by regulating the editing and splicing of chloroplast genes in rice.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01468-7.

The dicot homolog of maize PPR103 carries a C-terminal DYW domain and may have a role in C-to-U editing of some chloroplast RNA transcripts

In plants, cytidine-to-uridine (C-to-U) editing is a crucial step in processing mitochondria- and chloroplast-encoded transcripts. This editing requires nuclear-encoded proteins including members of the pentatricopeptide (PPR) family, especially PLS-type proteins carrying the DYW domain. IPI1/emb175/PPR103 is a nuclear gene encoding a PLS-type PPR protein essential for survival in Arabidopsis thaliana and maize. Arabidopsis IPI1 was identified as likely interacting with ISE2, a chloroplast-localized RNA helicase associated with C-to-U RNA editing in Arabidopsis and maize. Notably, while the Arabidopsis and Nicotiana IPI1 orthologs possess complete DYW motifs at their C-termini, the maize homolog, ZmPPR103, lacks this triplet of residues which are essential for editing. In this study we examined the function of IPI1 in chloroplast RNA processing in N. benthamiana to gain insight into the importance of the DYW domain to the function of the EMB175/PPR103/ IPI1 proteins. Structural predictions suggest that evolutionary loss of residues identified as critical for catalyzing C-to-U editing in other members of this class of proteins, were likely to lead to reduced or absent editing activity in the Nicotiana and Arabidopsis IPI1 orthologs. Virus-induced gene silencing of NbIPI1 led to defects in chloroplast ribosomal RNA processing and changes to stability of rpl16 transcripts, revealing conserved function with its maize ortholog. NbIPI1-silenced plants also had defective C-to-U RNA editing in several chloroplast transcripts, a contrast from the finding that maize PPR103 had no role in editing. The results indicate that in addition to its role in transcript stability, NbIPI1 may contribute to C-to-U editing in N. benthamiana chloroplasts.

A mitochondrial pentatricopeptide repeat protein enhances cold tolerance by modulating mitochondrial superoxide in rice

Cold stress affects rice growth and productivity. Defects in the plastid-localized pseudouridine synthase OsPUS1 affect chloroplast ribosome biogenesis, leading to low-temperature albino seedlings and accumulation of reactive oxygen species (ROS). Here, we report an ospus1-1 suppressor, sop10. SOP10 encodes a mitochondria-localized pentatricopeptide repeat protein. Mutations in SOP10 impair intron splicing of the nad4 and nad5 transcripts and decrease RNA editing efficiency of the nad2, nad6, and rps4 transcripts, resulting in deficiencies in mitochondrial complex I, thus decrease ROS generation and rescuing the albino phenotype. Overexpression of different compartment-localized superoxide dismutases (SOD) genes in ospus1-1 reverses the ROS over-accumulation and albino phenotypes to various degrees, with Mn-SOD reversing the best. Mutation of SOP10 in indica rice varieties enhances cold tolerance with lower ROS levels. We find that the mitochondrial superoxide plays a key role in rice cold responses, and identify a mitochondrial superoxide modulating factor, informing efforts to improve rice cold tolerance.

Comparative Root Transcriptome Profiling and Gene Regulatory Network Analysis between Eastern and Western Carrot (Daucus carota L.) Cultivars Reveals Candidate Genes for Vascular Tissue Patterning

Carrot (Daucus carota L.) is a highly consumed vegetable rich in carotenoids, known for their potent antioxidant, anti-inflammatory, and immune-protecting properties. While genetic and molecular studies have largely focused on wild and Western carrot cultivars (cvs), little is known about the evolutionary interactions between closely related Eastern and Western cvs. In this study, we conducted comparative transcriptome profiling of root tissues from Eastern (UHSBC-23-1) and Western (UHSBC-100) carrot cv. to better understand differentially expressed genes (DEGs) associated with storage root development and vascular cambium (VC) tissue patterning. Through reference-guided TopHat mapping, we achieved an average mapping rate of 73.87% and identified a total of 3544 DEGs (p < 0.05). Functional annotation and gene ontology classification revealed 97 functional categories, including 33 biological processes, 19 cellular components, 45 metabolic processes, and 26 KEGG pathways. Notably, Eastern cv. exhibited enrichment in cell wall, plant-pathogen interaction, and signal transduction terms, while Western cv. showed dominance in photosynthesis, metabolic process, and carbon metabolism terms. Moreover, constructed gene regulatory network (GRN) for both cvs. obtained orthologs with 1222 VC-responsive genes of Arabidopsis thaliana. In Western cv, GRN revealed VC-responsive gene clusters primarily associated with photosynthetic processes and carbon metabolism. In contrast, Eastern cv. exhibited a higher number of stress-responsive genes, and transcription factors (e.g., MYB15, WRKY46, AP2/ERF TF connected via signaling pathways with NAC036) were identified as master regulators of xylem vessel differentiation and secondary cell wall thickening. By elucidating the comparative transcriptome profiles of Eastern and Western cvs. for the first time, our study provides valuable insights into the differentially expressed genes involved in root development and VC tissue patterning. The identification of key regulatory genes and their roles in these processes represents a significant advancement in our understanding of the evolutionary relations and molecular mechanisms underlying secondary growth of carrot and regulation by vascular cambium.

Unravelling the genetic framework associated with grain quality and yield-related traits in maize (Zea mays L.)

Maize serves as a crucial nutrient reservoir for a significant portion of the global population. However, to effectively address the growing world population's hidden hunger, it is essential to focus on two key aspects: biofortification of maize and improving its yield potential through advanced breeding techniques. Moreover, the coordination of multiple targets within a single breeding program poses a complex challenge. This study compiled mapping studies conducted over the past decade, identifying quantitative trait loci associated with grain quality and yield related traits in maize. Meta-QTL analysis of 2,974 QTLs for 169 component traits (associated with quality and yield related traits) revealed 68 MQTLs across different genetic backgrounds and environments. Most of these MQTLs were further validated using the data from genome-wide association studies (GWAS). Further, ten MQTLs, referred to as breeding-friendly MQTLs (BF-MQTLs), with a significant phenotypic variation explained over 10% and confidence interval less than 2 Mb, were shortlisted. BF-MQTLs were further used to identify potential candidate genes, including 59 genes encoding important proteins/products involved in essential metabolic pathways. Five BF-MQTLs associated with both quality and yield traits were also recommended to be utilized in future breeding programs. Synteny analysis with wheat and rice genomes revealed conserved regions across the genomes, indicating these hotspot regions as validated targets for developing biofortified, high-yielding maize varieties in future breeding programs. After validation, the identified candidate genes can also be utilized to effectively model the plant architecture and enhance desirable quality traits through various approaches such as marker-assisted breeding, genetic engineering, and genome editing.

Protein-Protein Interactions and Quantitative Phosphoproteomic Analysis Reveal Potential Mitochondrial Substrates of Protein Phosphatase 2A-B'ζ Holoenzyme

Protein phosphatase 2A (PP2A) is a heterotrimeric conserved serine/threonine phosphatase complex that includes catalytic, scaffolding, and regulatory subunits. The 3 A subunits, 17 B subunits, and 5 C subunits that are encoded by the Arabidopsis genome allow 255 possible PP2A holoenzyme combinations. The regulatory subunits are crucial for substrate specificity and PP2A complex localization and are classified into the B, B', and B" non-related families in land plants. In Arabidopsis, the close homologs B'η, B'θ, B'γ, and B'ζ are further classified into a subfamily of B' called B'η. Previous studies have suggested that mitochondrial targeted PP2A subunits (B'ζ) play a role in energy metabolism and plant innate immunity. Potentially, the PP2A-B'ζ holoenzyme is involved in the regulation of the mitochondrial succinate/fumarate translocator, and it may affect the enzymes involved in energy metabolism. To investigate this hypothesis, the interactions between PP2A-B'ζ and the enzymes involved in the mitochondrial energy flow were investigated using bimolecular fluorescence complementation in tobacco and onion cells. Interactions were confirmed between the B'ζ subunit and the Krebs cycle proteins succinate/fumarate translocator (mSFC1), malate dehydrogenase (mMDH2), and aconitase (ACO3). Additional putative interacting candidates were deduced by comparing the enriched phosphoproteomes of wild type and B'ζ mutants: the mitochondrial regulator Arabidopsis pentatricopeptide repeat 6 (PPR6) and the two metabolic enzymes phosphoenolpyruvate carboxylase (PPC3) and phosphoenolpyruvate carboxykinase (PCK1). Overall, this study identifies potential PP2A substrates and highlights the role of PP2A in regulating energy metabolism in mitochondria.

Characterization of the Liriodendron chinense Pentatricopeptide Repeat (PPR) Gene Family and Its Role in Osmotic Stress Response

The Pentatricopeptide repeat (PPR) superfamily is a large gene family in plants that regulates organelle RNA metabolism, which is important for plant growth and development. However, a genome-wide analysis of the PPR gene family and its response to abiotic stress has not been reported for the relict woody plant Liriodendron chinense. In this paper, we identified 650 PPR genes from the L. chinense genome. A phylogenetic analysis showed that the LcPPR genes could roughly be divided into the P and PLS subfamilies. We found that 598 LcPPR genes were widely distributed across 19 chromosomes. An intraspecies synteny analysis indicated that duplicated genes from segmental duplication contributed to the expansion of the LcPPR gene family in the L. chinense genome. In addition, we verified the relative expression of Lchi03277, Lchi06624, Lchi18566, and Lchi23489 in the roots, stems, and leaves and found that all four genes had the highest expression in the leaves. By simulating a drought treatment and quantitative reverse transcription PCR (qRT-PCR) analysis, we confirmed the drought-responsive transcriptional changes in four LcPPR genes, two of which responded to drought stress independent of endogenous ABA biosynthesis. Thus, our study provides a comprehensive analysis of the L. chinense PPR gene family. It contributes to research into their roles in this valuable tree species' growth, development, and stress resistance.

A mitochondrion-associated PPR protein, WBG1, regulates grain chalkiness in rice

Rice kernel quality has vital commercial value. Grain chalkiness deteriorates rice's appearance and palatability. However, the molecular mechanisms that govern grain chalkiness remain unclear and may be regulated by many factors. In this study, we identified a stable hereditary mutant, white belly grain 1 (wbg1), which has a white belly in its mature grains. The grain filling rate of wbg1 was lower than that of the wild type across the whole filling period, and the starch granules in the chalky part were oval or round and loosely arranged. Map-based cloning showed that wbg1 was an allelic mutant of FLO10, which encodes a mitochondrion-targeted P-type pentatricopeptide repeat protein. Amino acid sequence analysis found that two PPR motifs present in the C-terminal of WBG1 were lost in wbg1. This deletion reduced the splicing efficiency of nad1 intron 1 to approximately 50% in wbg1, thereby partially reducing the activity of complex I and affecting ATP production in wbg1 grains. Furthermore, haplotype analysis showed that WBG1 was associated with grain width between indica and japonica rice varieties. These results suggested that WBG1 influences rice grain chalkiness and grain width by regulating the splicing efficiency of nad1 intron 1. This deepens understanding of the molecular mechanisms governing rice grain quality and provides theoretical support for molecular breeding to improve rice quality.

The dicot homolog of maize PPR103 carries a C-terminal DYW domain and is required for C-to-U editing of chloroplast RNA transcripts

In plants, cytidine-to-uridine (C-to-U) editing is a crucial step in processing mitochondria and chloroplast-encoded transcripts. This editing requires nuclear-encoded proteins including members of the pentatricopeptide (PPR) family, especially PLS-type proteins carrying the DYW domain. IPI1/emb175/PPR103 is a nuclear gene encoding a PLS-type PPR protein essential for survival in Arabidopsis thaliana and maize. Arabidopsis IPI1 was identified as likely interacting with ISE2, a chloroplast-localized RNA helicase associated with C-to-U RNA editing in Arabidopsis and maize. Notably, while the Arabidopsis and Nicotiana IPI1 homologs possess complete DYW motifs at their C-termini, the maize homolog, ZmPPR103, lacks this triplet of residues which are essential for editing. We examined the function of ISE2 and IPI1 in chloroplast RNA processing in N. benthamiana. A combination of deep sequencing and Sanger sequencing revealed C-to-U editing at 41 sites in 18 transcripts, with 34 sites conserved in the closely related N. tabacum. Virus induced gene silencing of NbISE2 or NbIPI1 led to defective C-to-U revealed that they have overlapping roles at editing a site in the rpoB transcript but have distinct roles in editing other transcripts. This finding contrasts with maize ppr103 mutants that showed no defects in editing. The results indicate that NbISE2 and NbIPI1 are important for C-to-U editing in N. benthamiana chloroplasts, and they may function in a complex to edit specific sites while having antagonistic effects on editing others. That NbIPI1, carrying a DYW domain, is involved in organelle C-to-U RNA editing supports previous work showing that this domain catalyzes RNA editing.

KRGG1 function in RNA editing in Trypanosoma brucei

Mitochondrial gene expression in trypanosomes requires numerous multiprotein complexes that are unique to kinetoplastids. Among these, the most well characterized are RNA editing catalytic complexes (RECCs) that catalyze the guide RNA (gRNA)-specified insertion and deletion of uridines during mitochondrial mRNA maturation. This post-transcriptional resequencing of mitochondrial mRNAs can be extensive, involving dozens of different gRNAs and hundreds of editing sites with most of the mature mRNA sequences resulting from the editing process. Proper coordination of the editing with the cognate gRNAs is attributed to RNA editing substrate-binding complexes (RESCs), which are also required for RNA editing. Although the precise mechanism of RESC function is less well understood, their affinity for binding both editing substrates and products suggests that these complexes may provide a scaffold for RECC catalytic processing. KRGG1 has been shown to bind RNAs, and although affinity purification co-isolates RESC complexes, its role in RNA editing remains uncertain. We show here that KRGG1 is essential in BF parasites and required for normal editing. KRGG1 repression results in reduced amounts of edited A6 mRNA and increased amounts of edited ND8 mRNA. Sequence and structure analysis of KRGG1 identified a region of homology with RESC6, and both proteins have predicted tandem helical repeats that resemble ARM/HEAT motifs. The ARM/HEAT-like region is critical for function as exclusive expression of mutated KRGG1 results in growth inhibition and disruption of KRGG1 association with RESCs. These results indicate that KRGG1 is critical for RNA editing and its specific function is associated with RESC activity.

OsPPR11 encoding P-type PPR protein that affects group II intron splicing and chloroplast development

OsPPR11 belongs to the P-type PPR protein family and can interact with OsCAF2 to regulate Group II intron splicing and affect chloroplast development in rice. Pentatricopeptide repeat (PPR) proteins participate in chloroplasts or mitochondria group II introns splicing in plants. The PPR protein family contains 491 members in rice, but most of their functions are unknown. In this study, we identified a nuclear gene encoding the P-type PPR protein OsPPR11 in chloroplasts. The qRT-PCR analysis demonstrated that OsPPR11 was expressed in all plant tissues, but leaves had the highest expression. The osppr11 mutants had yellowing leaves and a lethal phenotype that inhibited chloroplast development and photosynthesis-related gene expression and reduced photosynthesis-related protein accumulation in seedlings. Moreover, photosynthetic complex accumulation decreased significantly in osppr11 mutants. The OsPPR11 is required for ndhA, and ycf3-1 introns splicing and interact with CRM family protein OsCAF2, suggesting that these two proteins may form splicing complexes to regulate group II introns splicing. Further analysis revealed that OsCAF2 interacts with OsPPR11 through the N-terminus. These results indicate that OsPPR11 is essential for chloroplast development and function by affecting group II intron splicing in rice.

OsPPR11 encoding P-type PPR protein that affects group II intron splicing and chloroplast development

OsPPR11 belongs to the P-type PPR protein family and can interact with OsCAF2 to regulate Group II intron splicing and affect chloroplast development in rice. Pentatricopeptide repeat (PPR) proteins participate in chloroplasts or mitochondria group II introns splicing in plants. The PPR protein family contains 491 members in rice, but most of their functions are unknown. In this study, we identified a nuclear gene encoding the P-type PPR protein OsPPR11 in chloroplasts. The qRT-PCR analysis demonstrated that OsPPR11 was expressed in all plant tissues, but leaves had the highest expression. The osppr11 mutants had yellowing leaves and a lethal phenotype that inhibited chloroplast development and photosynthesis-related gene expression and reduced photosynthesis-related protein accumulation in seedlings. Moreover, photosynthetic complex accumulation decreased significantly in osppr11 mutants. The OsPPR11 is required for ndhA, and ycf3-1 introns splicing and interact with CRM family protein OsCAF2, suggesting that these two proteins may form splicing complexes to regulate group II introns splicing. Further analysis revealed that OsCAF2 interacts with OsPPR11 through the N-terminus. These results indicate that OsPPR11 is essential for chloroplast development and function by affecting group II intron splicing in rice.

Pentatricopeptide Protein PTCD2 Regulates COIII Translation in Mitochondria of the HeLa Cell Line

Protein biosynthesis in mitochondria is tightly coupled with assembly of inner membrane complexes and therefore must be coordinated with cytosolic translation of the mRNAs corresponding to the subunits which are encoded in the nucleus. Molecular mechanisms underlying the regulation of mitochondrial translation remain unclear despite recent advances in structural biology. Until now, only one translational regulator of protein biosynthesis in mammalian mitochondria is known-protein TACO1, which regulates translation of COI mRNA. Here we describe the function of pentatricopeptide-containing protein PTCD2 as a translational regulator of another mitochondrially encoded subunit of cytochrome c oxidase-COIII in the HeLa cell line. Deletion of the PTCD2 gene leads to significant decrease in COIII translation efficiency and impairment in CIV activity. Additionally, we show that PTCD2 protein is partially co-sedimentates with associated mitochondrial ribosome and associates with mitochondrial ribosome proteins in pull-down assays. These data allow concluding that PTCD2 is a specific translational regulator of COIII which attracts the mRNA to the mitochondrial ribosome.

Bulk segregant analysis-sequencing and RNA-Seq analyses reveal candidate genes associated with albino phenotype in Brassica napus

Inheritable albino mutants are excellent models for exploring the mechanism of chloroplast biogenesis and development. However, only a few non-lethal albino mutations have been reported to date in Brassica species. Here, we describe a resynthesized Brassica napus mutant, whose leaf, stem, and silique tissues showed an inheritable albino phenotype under field conditions after the bud stage but green phenotype in the greenhouse during the whole growing season, indicating that the albino phenotype depends on environmental conditions. Compared with the green leaves of the field-grown wild-type (GL) and greenhouse-grown mutant (WGL) plants, white leaves of the field-grown mutant (WL) showed significantly lower chlorophyll contents and structural defects in chloroplasts. Genetic analysis revealed that the albino phenotype of WL is recessive and is controlled by multiple genes. Bulk segregant analysis-sequencing (BSA-Seq) indicated that the candidate regions responsible for the albino phenotype spanned a total physical distance of approximately 49.68 Mb on chromosomes A03, A07, A08, C03, C04, C06, and C07. To gain insights into the molecular mechanisms that control chloroplast development in B. napus, we performed transcriptome (RNA-Seq) analysis of GL, WGL, and WL samples. GO and KEGG enrichment analyses suggested that differentially expressed genes (DEGs) associated with leaf color were significantly enriched in photosynthesis, ribosome biogenesis and chlorophyll metabolism. Further analysis indicated that DEGs involved in chloroplast development and chlorophyll metabolism were likely the main factors responsible for the albino phenotype in B. napus. A total of 59 DEGs were screened in the candidate regions, and four DEGs (BnaC03G0522600NO, BnaC07G0481600NO, BnaC07G0497800NO, and BnaA08G0016300NO) were identified as the most likely candidates responsible for the albino phenotype. Altogether, this study provides clues for elucidating the molecular mechanisms underlying chloroplast development in B. napus.

Evolution of mitochondrial RNA editing in extant gymnosperms

To unveil the evolution of mitochondrial RNA editing in gymnosperms, we characterized mitochondrial genomes (mitogenomes), plastid genomes, RNA editing sites, and pentatricopeptide repeat (PPR) proteins from 10 key taxa representing four of the five extant gymnosperm clades. The assembled mitogenomes vary in gene content due to massive gene losses in Gnetum and Conifer II clades. Mitochondrial gene expression levels also vary according to protein function, with the most highly expressed genes involved in the respiratory complex. We identified 9132 mitochondrial C-to-U editing sites, as well as 2846 P-class and 8530 PLS-class PPR proteins. Regains of editing sites were demonstrated in Conifer II rps3 transcripts whose corresponding mitogenomic sequences lack introns due to retroprocessing. Our analyses reveal that non-synonymous editing is efficient and results in more codons encoding hydrophobic amino acids. In contrast, synonymous editing, although performed with variable efficiency, can increase the number of U-ending codons that are preferentially utilized in gymnosperm mitochondria. The inferred loss-to-gain ratio of mitochondrial editing sites in gymnosperms is 2.1:1, of which losses of non-synonymous editing are mainly due to genomic C-to-T substitutions. However, such substitutions only explain a small fraction of synonymous editing site losses, indicating distinct evolutionary mechanisms. We show that gymnosperms have experienced multiple lineage-specific duplications in PLS-class PPR proteins. These duplications likely contribute to accumulated RNA editing sites, as a mechanistic correlation between RNA editing and PLS-class PPR proteins is statistically supported.

Screening and Verification of Photosynthesis and Chloroplast-Related Genes in Mulberry by Comparative RNA-Seq and Virus-Induced Gene Silencing

Photosynthesis is one of the most important factors in mulberry growth and production. To study the photosynthetic regulatory network of mulberry we sequenced the transcriptomes of two high-yielding (E1 and E2) and one low-yielding (H32) mulberry genotypes at two-time points (10:00 and 12:00). Re-annotation of the mulberry genome based on the transcriptome sequencing data identified 22,664 high-quality protein-coding genes with a BUSCO-assessed completeness of 93.4%. A total of 6587 differentially expressed genes (DEGs) were obtained in the transcriptome analysis. Functional annotation and enrichment revealed 142 out of 6587 genes involved in the photosynthetic pathway and chloroplast development. Moreover, 3 out of 142 genes were further examined using the VIGS technique; the leaves of MaCLA1- and MaTHIC-silenced plants were markedly yellowed or even white, and the leaves of MaPKP2-silenced plants showed a wrinkled appearance. The expression levels of the ensiled plants were reduced, and the levels of chlorophyll b and total chlorophyll were lower than those of the control plants. Co-expression analysis showed that MaCLA1 was co-expressed with CHUP1 and YSL3; MaTHIC was co-expressed with MaHSP70, MaFLN1, and MaEMB2794; MaPKP2 was mainly co-expressed with GH9B7, GH3.1, and EDA9. Protein interaction network prediction revealed that MaCLA1 was associated with RPE, TRA2, GPS1, and DXR proteins; MaTHIC was associated with TH1, PUR5, BIO2, and THI1; MaPKP2 was associated with ENOC, LOS2, and PGI1. This study offers a useful resource for further investigation of the molecular mechanisms involved in mulberry photosynthesis and preliminary insight into the regulatory network of photosynthesis.

Chloroplastic pentatricopeptide repeat proteins (PPR) in albino plantlets of Agave angustifolia Haw. reveal unexpected behavior

BACKGROUND

Pentatricopeptide repeat (PPR) proteins play an essential role in the post-transcriptional regulation of genes in plastid genomes. Although important advances have been made in understanding the functions of these genes, there is little information available on chloroplastic PPR genes in non-model plants and less in plants without chloroplasts. In the present study, a comprehensive and multifactorial bioinformatic strategy was applied to search for putative PPR genes in the foliar and meristematic tissues of green and albino plantlets of the non-model plant Agave angustifolia Haw.

RESULTS

A total of 1581 PPR transcripts were identified, of which 282 were chloroplastic. Leaf tissue in the albino plantlets showed the highest levels of expression of chloroplastic PPRs. The search for hypothetical targets of 12 PPR sequences in the chloroplast genes of A. angustifolia revealed their action on transcripts related to ribosomes and translation, photosystems, ATP synthase, plastid-encoded RNA polymerase and RuBisCO.

CONCLUSIONS

Our results suggest that the expression of PPR genes depends on the state of cell differentiation and plastid development. In the case of the albino leaf tissue, which lacks functional chloroplasts, it is possible that anterograde and retrograde signaling networks are severely compromised, leading to a compensatory anterograde response characterized by an increase in the expression of PPR genes.

A cytosolic pentatricopeptide repeat protein is essential for tapetal plastid development by regulating OsGLK1 transcript levels in rice

Most plant pentatricopeptide repeat (PPR) proteins localize to and function inside plastids and mitochondria. However, the function of PPRs that only localize to the cytoplasm remains unknown. Here, we demonstrated that the rice (Oryza sativa) PPR protein CYTOPLASM-LOCALIZED PPR1 (OsCPPR1) contributes to pollen development and localizes to the cytoplasm. Knocking down OsCPPR1 led to abnormal plastid development in tapetal cells, prolonged tapetal programmed cell death (PCD) and tapetum degradation, and significantly reduced pollen fertility. Transcriptome analysis revealed that the transcript level of OsGOLDEN-LIKE1 (OsGLK1), which encodes a transcription factor that regulates plastid development and maintenance, was significantly higher in the OsCPPR1 knockdown plants compared to wild-type plants. We further determined that OsCPPR1 downregulates OsGLK1 transcription by directly binding to the single-stranded regions of OsGLK1 mRNAs. Overexpression of OsGLK1 resulted in abnormal tapetum and plastid development, similar to that seen in OsCPPR1 knockdown plants, and suppression of OsGLK1 partially restored pollen fertility in the OsCPPR1 knockdown plants. We therefore conclude that OsCPPR1 suppresses OsGLK1 in the regulation of plastid development and PCD in the tapetum. Our work revealed novel functions for a cytosolic PPR, demonstrating the diverse roles of PPRs in plants and identifying a new regulatory mechanism for regulating pollen development in rice.

An Insight Into Pentatricopeptide-Mediated Chloroplast Necrosis via microRNA395a During Rhizoctonia solani Infection

Sheath blight (ShB) disease, caused by Rhizoctonia solani, is one of the major biotic stress-oriented diseases that adversely affect the rice productivity worldwide. However, the regulatory mechanisms are not understood yet comprehensively. In the current study, we had investigated the potential roles of miRNAs in economically important indica rice variety Pusa Basmati-1 upon R. solani infection by carrying out in-depth, high-throughput small RNA sequencing with a total data size of 435 million paired-end raw reads from rice leaf RNA samples collected at different time points. Detailed data analysis revealed a total of 468 known mature miRNAs and 747 putative novel miRNAs across all the libraries. Target prediction and Gene Ontology functional analysis of these miRNAs were found to be unraveling various cellular, molecular, and biological functions by targeting various plant defense-related genes. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) was performed to validate the miRNAs and their putative target genes. Out of the selected miRNA-specific putative target genes, miR395a binding and its cleavage site on pentatricopeptide were determined by 5’ RACE-PCR. It might be possible that R. solani instigated chloroplast degradation by modulating the pentatricopeptide which led to increased susceptibility to fungal infection.

Genome-Wide Identification and Characterization of RNA/DNA Differences Associated with Drought Response in Wheat

RNA/DNA difference (RDD) is a post-transcriptional RNA modification to enrich genetic information, widely involved in regulating diverse biological processes in eukaryotes. RDDs in the wheat nuclear genome, especially those associated with drought response or tolerance, were not well studied up to now. In this study, we investigated the RDDs related to drought response based on the RNA-seq data of drought-stressed and control samples in wheat. In total, 21,782 unique RDDs were identified, of which 265 were found to be drought-induced, representing the first drought-responsive RDD landscape in the wheat nuclear genome. The drought-responsive RDDs were located in 69 genes, of which 35 were differentially expressed under drought stress. Furthermore, the effects of RNA/DNA differences were investigated, showing that they could result in changes of RNA secondary structure, miRNA-target binding as well as protein conserved domains in the RDD-containing genes. In particular, the A to C mutation in TraesCS2A02G053100 (orthology to OsRLCK) led to the loss of tae-miR9657b-5p targeting, indicating that RNA/DNA difference might mediate miRNA to regulate the drought-response process. This study reported the first drought-responsive RDDs in the wheat nuclear genome. It sheds light on the roles of RDD in drought tolerance, and may also contribute to wheat genetic improvement based on epi-transcriptome methods.

Exploiting sterility and fertility variation in cytoplasmic male sterile vegetable crops

Cytoplasmic male sterility (CMS) has long been used to economically produce hybrids that harness growth vigor through heterosis. Yet, how CMS systems operate within commercially viable seed production strategies in various economically important vegetable crops, and their underlying molecular mechanisms, are often overlooked details that could expand the utility of CMS as a cost-effective and stable system. We provide here an update on the nature of cytoplasmic-nuclear interplay for pollen sterility and fertility transitions in vegetable crops, based on the discovery of components of nuclear fertility restoration and reversion determinants. Within plant CMS systems, pollen fertility can be rescued by the introduction of nuclear fertility restorer genes (Rfs), which operate by varied mechanisms to countermand the sterility phenotype. By understanding these systems, it is now becoming feasible to achieve fertility restoration with Rfs designed for programmable CMS-associated open reading frames (ORFs). Likewise, new opportunities exist for targeted disruption of CMS-associated ORFs by mito-TALENs in crops where natural Rfs have not been readily identified, providing an alternative approach to recovering fertility of cytoplasmic male sterile lines in crops. Recent findings show that facultative gynodioecy, as a reproductive strategy, can coordinate the sterility and fertility transition in response to environmental cues and/or metabolic signals that reflect ecological conditions of reproductive isolation. This information is important to devising future systems that are more inherently stable.

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

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

Plastid Deficient 1 Is Essential for the Accumulation of Plastid-Encoded RNA Polymerase Core Subunit β and Chloroplast Development in Arabidopsis

Plastid-encoded RNA polymerase (PEP)-dependent transcription is an essential process for chloroplast development and plant growth. It is a complex event that is regulated by numerous nuclear-encoded proteins. In order to elucidate the complex regulation mechanism of PEP activity, identification and characterization of PEP activity regulation factors are needed. Here, we characterize Plastid Deficient 1 (PD1) as a novel regulator for PEP-dependent gene expression and chloroplast development in Arabidopsis. The PD1 gene encodes a protein that is conserved in photoautotrophic organisms. The Arabidopsis pd1 mutant showed albino and seedling-lethal phenotypes. The plastid development in the pd1 mutant was arrested. The PD1 protein localized in the chloroplasts, and it colocalized with nucleoid protein TRXz. RT-quantitative real-time PCR, northern blot, and run-on analyses indicated that the PEP-dependent transcription in the pd1 mutant was dramatically impaired, whereas the nuclear-encoded RNA polymerase-dependent transcription was up-regulated. The yeast two-hybrid assays and coimmunoprecipitation experiments showed that the PD1 protein interacts with PEP core subunit β (PEP-β), which has been verified to be essential for chloroplast development. The immunoblot analysis indicated that the accumulation of PEP-β was barely detected in the pd1 mutant, whereas the accumulation of the other essential components of the PEP complex, such as core subunits α and β′, were not affected in the pd1 mutant. These observations suggested that the PD1 protein is essential for the accumulation of PEP-β and chloroplast development in Arabidopsis, potentially by direct interaction with PEP-β.

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

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

Molecular Analysis Uncovers the Mechanism of Fertility Restoration in Temperature-Sensitive Polima Cytoplasmic Male-Sterile Brassica napus

Temperature-sensitive male sterility is a heritable agronomic trait affected by genotype-environment interactions. In rapeseed (Brassica napus), Polima (pol) temperature-sensitive cytoplasmic male sterility (TCMS) is commonly used for two-line breeding, as the fertility of pol TCMS lines can be partially restored at certain temperatures. However, little is known about the underlying molecular mechanism that controls fertility restoration. Therefore, we aimed to investigate the fertility conversion mechanism of the pol TCMS line at two different ambient temperatures (16 °C and 25 °C). Our results showed that the anthers developed and produced vigorous pollen at 16 °C but not at 25 °C. In addition, we identified a novel co-transcript of orf224-atp6 in the mitochondria that might lead to fertility conversion of the pol TCMS line. RNA-seq analysis showed that 1637 genes were significantly differentially expressed in the fertile flowers of 596-L when compared to the sterile flower of 1318 and 596-H. Detailed analysis revealed that differentially expressed genes were involved in temperature response, ROS accumulation, anther development, and mitochondrial function. Single-molecule long-read isoform sequencing combined with RNA sequencing revealed numerous genes produce alternative splicing transcripts at high temperatures. Here, we also found that alternative oxidase, type II NAD(P)H dehydrogenases, and transcription factor Hsfs might play a crucial role in male fertility under the low-temperature condition. RNA sequencing and bulked segregant analysis coupled with whole-genome sequencing identified the candidate genes involved in the post-transcriptional modification of orf224. Overall, our study described a putative mechanism of fertility restoration in a pol TCMS line controlled by ambient temperature that might help utilise TCMS in the two-line breeding of Brassica crops.

Organellar Introns in Fungi, Algae, and Plants

Introns are ubiquitous in eukaryotic genomes and have long been considered as ‘junk RNA’ but the huge energy expenditure in their transcription, removal, and degradation indicate that they may have functional significance and can offer evolutionary advantages. In fungi, plants and algae introns make a significant contribution to the size of the organellar genomes. Organellar introns are classified as catalytic self-splicing introns that can be categorized as either Group I or Group II introns. There are some biases, with Group I introns being more frequently encountered in fungal mitochondrial genomes, whereas among plants Group II introns dominate within the mitochondrial and chloroplast genomes. Organellar introns can encode a variety of proteins, such as maturases, homing endonucleases, reverse transcriptases, and, in some cases, ribosomal proteins, along with other novel open reading frames. Although organellar introns are viewed to be ribozymes, they do interact with various intron- or nuclear genome-encoded protein factors that assist in the intron RNA to fold into competent splicing structures, or facilitate the turn-over of intron RNAs to prevent reverse splicing. Organellar introns are also known to be involved in non-canonical splicing, such as backsplicing and trans-splicing which can result in novel splicing products or, in some instances, compensate for the fragmentation of genes by recombination events. In organellar genomes, Group I and II introns may exist in nested intronic arrangements, such as introns within introns, referred to as twintrons, where splicing of the external intron may be dependent on splicing of the internal intron. These nested or complex introns, with two or three-component intron modules, are being explored as platforms for alternative splicing and their possible function as molecular switches for modulating gene expression which could be potentially applied towards heterologous gene expression. This review explores recent findings on organellar Group I and II introns, focusing on splicing and mobility mechanisms aided by associated intron/nuclear encoded proteins and their potential roles in organellar gene expression and cross talk between nuclear and organellar genomes. Potential application for these types of elements in biotechnology are also discussed.

The Chloroplast Epitranscriptome: Factors, Sites, Regulation, and Detection Methods

Modifications in nucleic acids are present in all three domains of life. More than 170 distinct chemical modifications have been reported in cellular RNAs to date. Collectively termed as epitranscriptome, these RNA modifications are often dynamic and involve distinct regulatory proteins that install, remove, and interpret these marks in a site-specific manner. Covalent nucleotide modifications-such as methylations at diverse positions in the bases, polyuridylation, and pseudouridylation and many others impact various events in the lifecycle of an RNA such as folding, localization, processing, stability, ribosome assembly, and translational processes and are thus crucial regulators of the RNA metabolism. In plants, the nuclear/cytoplasmic epitranscriptome plays important roles in a wide range of biological processes, such as organ development, viral infection, and physiological means. Notably, recent transcriptome-wide analyses have also revealed novel dynamic modifications not only in plant nuclear/cytoplasmic RNAs related to photosynthesis but especially in chloroplast mRNAs, suggesting important and hitherto undefined regulatory steps in plastid functions and gene expression. Here we report on the latest findings of known plastid RNA modifications and highlight their relevance for the post-transcriptional regulation of chloroplast gene expression and their role in controlling plant development, stress reactions, and acclimation processes.

Plant RNA Binding Proteins as Critical Modulators in Drought, High Salinity, Heat, and Cold Stress Responses: An Updated Overview

Plant abiotic stress responses are tightly regulated by different players at multiple levels. At transcriptional or post-transcriptional levels, several RNA binding proteins (RBPs) regulate stress response genes through RNA metabolism. They are increasingly recognized as critical modulators of a myriad of biological processes, including stress responses. Plant RBPs are heterogeneous with one or more conservative RNA motifs that constitute canonical/novel RNA binding domains (RBDs), which can bind to target RNAs to determine their regulation as per the plant requirements at given environmental conditions. Given its biological significance and possible consideration as a potential tool in genetic manipulation programs to improve key agronomic traits amidst frequent episodes of climate anomalies, studies concerning the identification and functional characterization of RBP candidate genes are steadily mounting. This paper presents a comprehensive overview of canonical and novel RBPs and their functions in major abiotic stresses including drought, heat, salt, and cold stress conditions. To some extent, we also briefly describe the basic motif structure of RBPs that would be useful in forthcoming studies. Additionally, we also collected RBP genes that were modulated by stress, but that lacked functional characterization, providing an impetus to conduct further research.

An Analytical Review of the Structural Features of Pentatricopeptide Repeats: Strategic Amino Acids, Repeat Arrangements and Superhelical Architecture

Tricopeptide repeats are common in natural proteins, and are exemplified by 34- and 35-residue repeats, known respectively as tetratricopeptide repeats (TPRs) and pentatricopeptide repeats (PPRs). In both classes, each repeat unit forms an antiparallel bihelical structure, so that multiple such units in a polypeptide are arranged in a parallel fashion. The primary structures of the motifs are nonidentical, but amino acids of similar properties occur in strategic positions. The focus of the present work was on PPR, but TPR, its better-studied cousin, is often included for comparison. The analyses revealed that critical amino acids, namely Gly, Pro, Ala and Trp, were placed at distinct locations in the higher order structure of PPR domains. While most TPRs occur in repeats of three, the PPRs exhibited a much greater diversity in repeat numbers, from 1 to 30 or more, separated by spacers of various sequences and lengths. Studies of PPR strings in proteins showed that the majority of PPR units are single, and that the longer tandems (i.e., without space in between) occurred in decreasing order. The multi-PPR domains also formed superhelical vortices, likely governed by interhelical angles rather than the spacers. These findings should be useful in designing and understanding the PPR domains.