Organelle Genes Research Articles

Mitochondrial Splicing Efficiency Is Lower in Holoparasites Than in Free-Living Plants.

Mitochondria play a crucial role in eukaryotic organisms, housing their own genome with genes vital for oxidative phosphorylation. Coordination between nuclear and mitochondrial genomes is pivotal for organelle gene expression. Splicing, editing and processing of mitochondrial transcripts are regulated by nuclear-encoded factors. Splicing efficiency (SEf) of the many group II introns present in plant mitochondrial genes is critical for mitochondrial function since a splicing defect or splicing deficiency can severely impact plant growth and development. This study investigates SEf in free-living and holoparasitic plants, focusing on 25 group II introns from 15 angiosperm species. Our comparative analyses reveal distinctive splicing patterns with holoparasites exhibiting significantly lower SEf, potentially linked to their unique evolutionary trajectory. Given the preponderance of horizontal gene transfer (HGT) in parasitic plants, we investigated the effect of HGT on SEf, such as the presence of foreign introns or foreign nuclear-encoded splicing factors. Contrary to expectations, the SEf reductions do not correlate with HGT events, suggesting that other factors are at play, such as the loss of photosynthesis or the transition to a holoparasitic lifestyle. The findings of this study broaden our understanding of the molecular evolution in parasitic plants and shed light on the multifaceted factors influencing organelle gene expression.

The Organelle in the Ointment: improved detection of cryptic mitochondrial reads resolves many unknown sequences in cross-species microbiome analyses

Abstract The genomes of mitochondria and chloroplasts contain ribosomal RNA (rRNA) genes, reflecting their ancestry as free-living bacteria. These organellar rRNAs are often amplified in microbiome studies of animals and plants. If identified, they can be discarded, merely reducing sequencing depth. However, we identify certain high-abundance organeller RNAs not identified by common pipelines, which may compromise statistical analysis of microbiome structure and diversity. We quantified this by reanalyzing 7459 samples from seven 16S rRNA studies, including microbiomes from 927 unique animal genera. We find that under-annotation of cryptic mitochondrial and chloroplast reads affects multiple of these large-scale cross-species microbiome comparisons, and varies between host species, biasing comparisons. We offer a straightforward solution: supplementing existing taxonomies with diverse organelle rRNA sequences. This resolves up to 97% of unique unclassified sequences in some entire studies as mitochondrial (14% averaged across all studies), without increasing false positive annotations in mitochondria-free mock communities. Improved annotation decreases the proportion of unknown sequences by ≥10-fold in 2262 of 7459 samples (30%), spanning 5 of 7 major studies examined. We recommend leveraging organelle sequence diversity to better identify organelle gene sequences in microbiome studies, and provide code, data resources and tutorials that implement this approach.

PPR21 is involved in the splicing of nad2 introns via interacting with PPR–small MutS-related 1 and small PPR protein 2 and is essential to maize seed development

Pentatricopeptide repeat (PPR) proteins are a large group of eukaryote-specific RNA-binding proteins that play pivotal roles in plant organelle gene expression. Here, we report the function of PPR21 in mitochondrial intron splicing and its role in maize kernel development. PPR21 is a typical P-type PPR protein targeted to mitochondria. The ppr21 mutants are arrested in embryogenesis and endosperm development, leading to embryo lethality. Null mutations of PPR21 reduce the splicing efficiency of nad2 intron 1, 2, and 4 and impair the assembly and activity of mitochondrial complex I. Previous studies show that the P-type PPR protein EMP12 is required for the splicing of identical introns. However, our protein interaction analyses reveal that PPR21 does not interact with EMP12. Instead, both PPR21 and EMP12 interact with the small MutS-related (SMR) domain–containing PPR protein 1 (PPR-SMR1) and the short P-type PPR protein 2 (SPR2). PPR-SMR1 interacts with SPR2, and both proteins are required for the splicing of many introns in mitochondria, including nad2 intron 1, 2, and 4. These results suggest that a PPR21–(PPR-SMR1/SPR2)–EMP12 complex is involved in the splicing of nad2 introns in maize mitochondria.

Transcriptomic Analysis Reveals Adaptive Evolution and Conservation Implications for the Endangered Magnolia lotungensis.

Magnolia lotungensis is an extremely endangered endemic tree in China. To elucidate the genetic basis of M. lotungensis, we performed a comprehensive transcriptome analysis using a sample integrating the plant's bark, leaves, and flowers. De novo transcriptome assembly yielded 177,046 transcripts and 42,518 coding sequences. Notably, we identified 796 species-specific genes enriched in organelle gene regulation and defense responses. A codon usage bias analysis revealed that mutation bias appears to be the primary driver of selection in shaping the species' genetic architecture. An evolutionary analysis based on dN/dS values of paralogous and orthologous gene pairs indicated a predominance of purifying selection, suggesting strong evolutionary constraints on most genes. A comparative transcriptomic analysis with Magnolia sinica identified approximately 1000 ultra-conserved genes, enriched in essential cellular processes such as transcriptional regulation, protein synthesis, and genome stability. Interestingly, only a limited number of 511 rapidly evolving genes under positive selection were detected compared to M. sinica and Magnolia kuangsiensis. These genes were enriched in metabolic processes associated with adaptation to specific environments, potentially limiting the species' ability to expand its range. Our findings contribute to understanding the genetic architecture of M. lotungensis and suggest that an insufficient number of adaptive genes contribute to its endangered status.

Malaria parasites require a divergent heme oxygenase for apicoplast gene expression and biogenesis.

Malaria parasites have evolved unusual metabolic adaptations that specialize them for growth within heme-rich human erythrocytes. During blood-stage infection, Plasmodium falciparum parasites internalize and digest abundant host hemoglobin within the digestive vacuole. This massive catabolic process generates copious free heme, most of which is biomineralized into inert hemozoin. Parasites also express a divergent heme oxygenase (HO)-like protein (PfHO) that lacks key active-site residues and has lost canonical HO activity. The cellular role of this unusual protein that underpins its retention by parasites has been unknown. To unravel PfHO function, we first determined a 2.8 Å-resolution X-ray structure that revealed a highly α-helical fold indicative of distant HO homology. Localization studies unveiled PfHO targeting to the apicoplast organelle, where it is imported and undergoes N-terminal processing but retains most of the electropositive transit peptide. We observed that conditional knockdown of PfHO was lethal to parasites, which died from defective apicoplast biogenesis and impaired isoprenoid-precursor synthesis. Complementation and molecular-interaction studies revealed an essential role for the electropositive N-terminus of PfHO, which selectively associates with the apicoplast genome and enzymes involved in nucleic acid metabolism and gene expression. PfHO knockdown resulted in a specific deficiency in levels of apicoplast-encoded RNA but not DNA. These studies reveal an essential function for PfHO in apicoplast maintenance and suggest that Plasmodium repurposed the conserved HO scaffold from its canonical heme-degrading function in the ancestral chloroplast to fulfill a critical adaptive role in organelle gene expression.

Ecological Predictors of Organelle Genome Evolution: Phylogenetic Correlations with Taxonomically Broad, Sparse, Unsystematized Data.

Comparative analysis of variables across phylogenetically linked observations can reveal mechanisms and insights in evolutionary biology. As the taxonomic breadth of the sample of interest increases, challenges of data sparsity, poor phylogenetic resolution, and complicated evolutionary dynamics emerge. Here, we investigate a cross-eukaryotic question where all these problems exist: which organismal ecology features are correlated with gene retention in mitochondrial and chloroplast DNA (organelle DNA or oDNA). Through a wide palette of synthetic control studies, we first characterize the specificity and sensitivity of a collection of parametric and non-parametric phylogenetic comparative approaches to identify relationships in the face of such sparse and awkward datasets. This analysis is not directly focused on oDNA, and so provides generalizable insights into comparative approaches with challenging data. We then combine and curate ecological data coupled to oDNA genome information across eukaryotes, including a new semi-automated approach for gathering data on organismal traits from less systematized open-access resources including encyclopedia articles on species and taxa. The curation process also involved resolving several issues with existing datasets, including enforcing the clade-specificity of several ecological features and fixing incorrect annotations. Combining this unique dataset with our benchmarked comparative approaches, we confirm support for several known links between organismal ecology and organelle gene retention, identify several previously unidentified relationships constituting possible ecological contributors to oDNA genome evolution, and provide support for a recently hypothesized link between environmental demand and oDNA retention. We, with caution, discuss the implications of these findings for organelle evolution and of this pipeline for broad comparative analyses in other fields.

Comparative analysis of the organelle genomes of Aconitum carmichaelii revealed structural and sequence differences and phylogenetic relationships

In this study, we conducted an assembly and analysis of the organelle genomes of Aconitum carmichaelii. Our investigation encompassed the examination of organelle genome structures, gene transfer events, and the environmental selection pressures affecting A. carmichaelii. The results revealed distinct evolutionary patterns in the organelle genomes of A. carmichaelii. Especially, the plastome exhibited a more conserved structure but a higher nucleotide substitution rate (NSR), while the mitogenome displayed a more complex structure with a slower NSR. Through homology analysis, we identified several instances of unidirectional protein-coding genes (PCGs) transferring from the plastome to the mitogenome. However, we did not observe any events which genes moved from the mitogenome to the plastome. Additionally, we observed multiple transposable element (TE) fragments in the organelle genomes, with both organelles showing different preferences for the type of nuclear TE insertion. Divergence time estimation suggested that rapid differentiation occurred in Aconitum species approximately 7.96 million years ago (Mya). This divergence might be associated with the reduction in CO2 levels and the significant uplift of the Qinghai-Tibet Plateau (QTP) during the late Miocene. Selection pressure analysis indicated that the dN/dS values of both organelles were less than 1, suggested that organelle PCGs were subject to purification selection. However, we did not detect any positively selected genes (PSGs) in Subg. Aconitum and Subg. Lycoctonum. This observation further supports the idea that stronger negative selection pressure on organelle genes in Aconitum results in a more conserved amino acid sequence. In conclusion, this study contributes to a deeper understanding of organelle evolution in Aconitum species and provides a foundation for future research on the genetic mechanisms underlying the structure and function of the Aconitum plastome and mitogenome.

Comparative analysis of chloroplast and mitochondrial genomes of sweet potato provides evidence of gene transfer

The increasing number of plant mitochondrial DNA genomes (mtDNA) sequenced reveals the extent of transfer from both chloroplast DNA genomes (cpDNA) and nuclear DNA genomes (nDNA). This study created a library and assembled the chloroplast and mitochondrial genomes of the leafy sweet potato better to understand the extent of mitochondrial and chloroplast gene transfer. The full-length chloroplast genome of the leafy sweet potato (OM808940) is 161,387 bp, with 132 genes annotated, including 87 protein-coding genes, 8 rRNA genes, and 37 tRNA genes. The mitochondrial genome (OM808941) was 269,578 bp in length and contained 69 functional genes, including 39 protein-coding genes, 6 rRNA genes, and 24 tRNA genes. 68 SSR loci were found in the leafy sweet potato organelle genome, including 54 in the chloroplast genome and 14 in the mitochondria genome. In the sweet potato mitochondrial genome, most genes have RNA editing sites, and the conversion ratio from hydrophilic amino acids to hydrophobic amino acids is the highest, reaching 47.12%. Horizontal transfer occurs in the sweet potato organelle genome and nuclear genome. 40 mitochondrial genome segments share high homology with 14 chloroplast genome segments, 33 of which may be derived from chloroplast genome horizontal transfer. 171 mitochondrial genome sequences come from the horizontal transfer of nuclear genome. The phylogenetic analysis of organelle genes revealed that the leafy sweet potato was closely related to the tetraploid wild species Ipomoea tabascana and the wild diploid species Ipomoea trifida.

Cytonuclear interplay in auto- and allopolyploids: amultifaceted perspective from the Festuca-Lolium complex.

Restoring cytonuclear stoichiometry is necessary after whole-genome duplication (WGD) and interspecific/intergeneric hybridization in plants. We investigated this phenomenon in auto- and allopolyploids of the Festuca-Lolium complex providing insights into the mechanisms governing cytonuclear interactions in early polyploid and hybrid generations. Our study examined the main processes potentially involved in restoring the cytonuclear balance after WGD comparing diploids and new and well-established autopolyploids. We uncovered that both the number of chloroplasts and the number of chloroplast genome copies were significantly higher in the newly established autopolyploids and grew further in more established autopolyploids. The increase in the copy number of the chloroplast genome exceeded the rise in the number of chloroplasts and fully compensated for the doubling of the nuclear genome. In addition, changes in nuclear and organelle gene expression were insignificant. Allopolyploid Festuca × Lolium hybrids displayed potential structural conflicts in parental protein variants within the cytonuclear complexes. While biased maternal allele expression has been observed in numerous hybrids, our results suggest that its role in cytonuclear stabilization in the Festuca × Lolium hybrids is limited. This study provides insights into the restoration of the cytonuclear stoichiometry, yet it emphasizes the need for future research to explore post-transcriptional regulation and its impact on cytonuclear gene expression stoichiometry. Our findings may enhance the understanding of polyploid plant evolution, with broader implications for the study of cytonuclear interactions in diverse biological contexts.

Alteration of Mitochondrial Transcript Expression in Arabidopsis thaliana Using a Custom-Made Library of Pentatricopeptide Repeat Proteins.

Pentatricopeptide repeat (PPR) proteins are considered a potential tool for manipulating organelle gene expression in plants because they can recognise a wide range of different RNA sequences, and the molecular basis for this sequence recognition is partially known and understood. A library of redesigned PPR proteins related to restorer-of-fertility proteins was created and transformed into plants in order to target mitochondrial transcripts. Ninety different variants tested in vivo showed a wide range of phenotypes. One of these lines, which displayed slow growth and downward curled leaves, showed a clear reduction in complex V. The phenotype was due to a specific cleavage of atp1 transcripts induced by a modified PPR protein from the library, validating the use of this library as a source of mitochondrial 'mutants'. This study is a step towards developing specific RNA targeting tools using PPR proteins that can be aimed at desired targets.

Transcriptome-Based Study on the Phylogeny and Hybridization of Marattialean Ferns (Marattiaceae).

Marattiaceae is a phylogenetically isolated family of tropical eusporangiate ferns including six genera with more than one-hundred species. In Marattiaceae, monophyly of genera has been well-supported phylogenetically. However, the phylogenetic relationships among them were elusive and controversial. Here, a dataset of 26 transcriptomes (including 11 newly generated) were used to assess single-copy nuclear genes and to obtain the organelle gene sequences. Through phylotranscriptomic analysis, the phylogeny and hybridization events of Marattiaceae were explored and a robust phylogenomic framework for the evolution of Marattiaceae was provided. Using both concatenation- and coalescent-based phylogenies, the gene-tree discordance, incomplete lineage sorting (ILS) simulations, and network inference were examined. Except the low support with mitochondrial genes of Marattiaceae, nuclear genes and chloroplast genes strongly supported a sister relationship between Marattiaceae and leptosporangiate ferns. At the genus level, all phylogenetic analysis based on nuclear genes datasets recovered five genera in Marattiaceae as monophyletic with strong support. Danaea and Ptisana were the first two diverged clades in turn. Christensenia was a sister clade to the clade Marattia + Angiopteris s.l. In Angiopteris s.l., three clades (Angiopteris s.s., the Archangiopteris group, and An. sparsisora) were well identified with maximum support. The Archangiopteris group was derived from Angiopteris s.s. at ca. 18 Ma. The putative hybrid species An. sparsisora between Angiopteris s.s. and the Archangiopteris group was verified by the species network analyses and the maternal plastid genes. This study will improve our understanding for using the phylotranscriptomic method to explore phylogeny and investigate hybridization events for difficult taxa in ferns.

Transcriptional activity of mitochondrial genes in intraspecific and interspecific sunflower hybrids

Relevance. The genetic structure of plant cells implies the coordinated work of three genomes: nuclear, plastid, and mitochondrial. Hybridization between genetically heterogeneous parents can lead to changes in the established nuclear-cytoplasmic balance, which in turn can affect the level and consistency of their gene expression. Changes in the transcriptional activity of organelle genes (in particular, mitochondria) during distant (interspecific) hybridization remain poorly understood. Results. The present study employed the qPCR technique to evaluate the transcriptional activity level of the mitochondrial genes atp1, atp4, atp6, atp9, nad3, nad6, cox1, and cox3 in intra- and interspecific sunflower hybrids and their parental forms from the VIR collection. According to the analyzed transcriptional activity of mitochondrial genes, they can be divided into three groups: genes with a relatively high level of expression – atp1, atp6, and nad6, those with a medium level of expression – atp4, cox1, cox3, and genes with a low level of expression – atp9 and nad3. Comparative analysis showed no significant difference (P<0.05) between maternal lines and hybrids. However, the expression of the nad6 gene in the case of Helianthus argophyllus (Torr. & A. Gray) was 2.6 times higher than in the cultivated sunflower lines. Conclusion. The absence of substantial changes in the expression of mitochondrial genes both in intra- and interspecific hybrids indicates the lack of significant changes in the regulation of nuclear-cytoplasmic interactions in these hybrids.

Cellular and environmental dynamics influence species-specific extents of organelle gene retention.

Mitochondria and plastids rely on many nuclear-encoded genes, but retain small subsets of the genes they need to function in their own organelle DNA (oDNA). Different species retain different numbers of oDNA genes, and the reasons for these differences are not completely understood. Here, we use a mathematical model to explore the hypothesis that the energetic demands imposed by an organism's changing environment influence how many oDNA genes it retains. The model couples the physical biology of cell processes of gene expression and transport to a supply-and-demand model for the environmental dynamics to which an organism is exposed. The trade-off between fulfilling metabolic and bioenergetic environmental demands, and retaining genetic integrity, is quantified for a generic gene encoded either in oDNA or in nuclear DNA. Species in environments with high-amplitude, intermediate-frequency oscillations are predicted to retain the most organelle genes, whereas those in less dynamic or noisy environments the fewest. We discuss support for, and insight from, these predictions with oDNA data across eukaryotic taxa, including high oDNA gene counts in sessile organisms exposed to day-night and intertidal oscillations (including plants and algae) and low counts in parasites and fungi.

Combination treatment of hepatitis C virus-associated hepatocellular carcinoma by simultaneously blocking genes in multiple organelles via functionally engineered graphene oxide

Chemical reagents are widely used to block disease development, including tumor progression and increase in abnormal cell populations; however, their consistent usage is limited owing to drug-related side effects and resistance. Although strategies to enhance the intracellular delivery of anticancer drugs, even at low concentrations, and decrease their undesired side effects have been developed, practical approaches for liver cancer are still limited. Combination therapy has become a solution for improving chemotherapeutic efficacy. Here, we designed a dual-therapeutic platform by engineering graphene oxide (GO) as a multi-functional co-delivery carrier that targets multiple organelles in cells. It delivers two therapeutic reagents as a single vehicle while simultaneously blocking genes in the nucleus and cytoplasm. We adopted a xenograft mouse model of hepatitis C virus infection to verify the therapeutic efficacy of this combinatorial approach. Our strategy successfully suppressed tumor growth by enhancing the intracellular accumulation of anticancer drugs even at one-tenth of the conventional dosage of drugs. Simultaneously, the dual-treat method showed a twofold therapeutic efficiency compared with single-drug treatment. By facilitating a combination of individual pharmacological effects, the present functional GO-based delivery system could be a promising multimodal therapeutic platform owing to its high stability, biocompatibility, and multifunctionality.

Chloroplast Genome Annotation Tools: Prolegomena to the Identification of Inverted Repeats.

The development of next-generation sequencing technology and the increasing amount of sequencing data have brought the bioinformatic tools used in genome assembly into focus. The final step of the process is genome annotation, which works on assembled genome sequences to identify the location of genome features. In the case of organelle genomes, specialized annotation tools are used to identify organelle genes and structural features. Numerous annotation tools target chloroplast sequences. Most chloroplast DNA genomes have a quadripartite structure caused by two copies of a large inverted repeat. We investigated the strategies of six annotation tools (Chloë, Chloroplot, GeSeq, ORG.Annotate, PGA, Plann) for identifying inverted repeats and analyzed their success using publicly available complete chloroplast sequences of taxa belonging to the asterid and rosid clades. The annotation tools use two different approaches to identify inverted repeats, using existing general search tools or implementing stand-alone solutions. The chloroplast sequences studied show that there are different types of imperfections in the assembled data and that each tool performs better on some sequences than the others.

MISF2 Encodes an Essential Mitochondrial Splicing Cofactor Required for nad2 mRNA Processing and Embryo Development in Arabidopsis thaliana.

Mitochondria play key roles in cellular energy metabolism in eukaryotes. Mitochondria of most organisms contain their own genome and specific transcription and translation machineries. The expression of angiosperm mtDNA involves extensive RNA-processing steps, such as RNA trimming, editing, and the splicing of numerous group II-type introns. Pentatricopeptide repeat (PPR) proteins are key players in plant organelle gene expression and RNA metabolism. In the present analysis, we reveal the function of the MITOCHONDRIAL SPLICING FACTOR 2 gene (MISF2, AT3G22670) and show that it encodes a mitochondria-localized PPR protein that is crucial for early embryo development in Arabidopsis. Molecular characterization of embryo-rescued misf2 plantlets indicates that the splicing of nad2 intron 1, and thus respiratory complex I biogenesis, are strongly compromised. Moreover, the molecular function seems conserved between MISF2 protein in Arabidopsis and its orthologous gene (EMP10) in maize, suggesting that the ancestor of MISF2/EMP10 was recruited to function in nad2 processing before the monocot–dicot divergence ~200 million years ago. These data provide new insights into the function of nuclear-encoded factors in mitochondrial gene expression and respiratory chain biogenesis during plant embryo development.

The Melatonin Receptor CAND2/PMTR1 Is Involved in the Regulation of Mitochondrial Gene Expression under Photooxidative Stress.

Melatonin is a signaling molecule that mediates multiple stress-dependent reactions. Under photooxidative stress conditions generating intensive ROS production, exogenous melatonin (50 µM) contributed to maintaining the expression of mitochondrial encoded genes and up-regulation of RNA-polymerase genes RPOTm and RPOTmp, operating through the CAND2 receptor and α-subunit of the heterotrimeric G protein GPA1 coupled with CAND2. Unlike wild-type plants, mutants with defective CAND2 and GPA1 genes exhibited no decrease in the alternative pathway of leaf respiration, as well as the activity of an alternative oxidase, and the expression of the AOX1a gene. At the same time, the protective effect of exogenous melatonin on some physiological indicators did not depend on the receptor and was associated with the direct antioxidant function of the regulator. Thus, melatonin under photooxidative stress conditions can act as an antioxidant and as a hormone capable of regulating the expression of nuclear and organelle genes through the components of melatonin signal perception.

Organelle genome assembly uncovers the dynamic genome reorganization and cytoplasmic male sterility associated genes in tomato

To identify cytoplasmic male sterility (CMS)-associated genes in tomato, we determined the genome sequences of mitochondria and chloroplasts in three CMS tomato lines derived from independent asymmetric cell fusions, their nuclear and cytoplasmic donors, and male fertile weedy cultivated tomato and wild relatives. The structures of the CMS mitochondrial genomes were highly divergent from those of the nuclear and cytoplasmic donors, and genes of the donors were mixed up in these genomes. On the other hand, the structures of CMS chloroplast genomes were moderately conserved across the donors, but CMS chloroplast genes were unexpectedly likely derived from the nuclear donors. Comparative analysis of the structures and contents of organelle genes and transcriptome analysis identified three genes that were uniquely present in the CMS lines, but not in the donor or fertile lines. RNA-sequencing analysis indicated that these three genes transcriptionally expressed in anther, and identified different RNA editing levels in one gene, orf265, that was partially similar to ATP synthase subunit 8, between fertile and sterile lines. The orf265 was a highly potential candidate for CMS-associated gene. This study suggests that organelle reorganization mechanisms after cell fusion events differ between mitochondria and chloroplasts, and provides insight into the development of new F1 hybrid breeding programs employing the CMS system in tomato.

The role of PPR proteins in posttranscriptional regulation of organelle components in plants.

Pentatricopeptide repeat (PPR) proteins constitute one of the largest protein families in land plants. They are sequence-specific RNA-binding proteins and play key roles in posttranscriptional processes within organelles. Their combined actions have profound effects on chloroplast photosynthetic electron transport chain and mitochondrial respiratory chain, affecting photosynthesis and respiration respectively, and ultimately on yield, fertility, and grain quality. Over the past decade, much has been learned about the molecular functions of these proteins on plant growth and development. However, due to the large size of this protein family, the functions of most members remain largely unknown. Here, we summarize the molecular mechanisms of PPR proteins functions on organelle genes, and effects on development of organelles and plants. Problems that need to be resolved are also identified. This article will provide a theoretical basis for understanding the functions of PPR protein family and genetic improvements of grain yield and quality.

PPR proteins in the tea plant (Camellia sinensis) and their potential roles in the leaf color changes

Pentatricopeptide repeat (PPR) family protein is one of the largest protein families in plants, which is involved in many biological processes, including transcription, RNA stabilization, RNA cleavage, RNA splicing, RNA editing, and translation of organelle genes, thereby affecting the expression of the organelle genes. However, the corresponding PPR protein in tea plants has not been investigated. In this study, 578 CsPPR genes were identified in tea plants (Camellia sinensis); the PPR motifs, types of PPR proteins, subcellular locations, gene structure, conserved motifs, phylogenetic relationships, and GO annotation were analyzed. Based on motif structure and evolutionary relationship, these genes were divided into the P subfamily and PLS subfamily. GO annotation indicated that most PPR proteins had binding activity, including RNA and DNA binding. Transcriptome and expression profile analyses suggested that 22 CsPPR genes might be related to the etiolation and albino process of tea plants. SNP analysis found that 28 SNP loci were related to the albino of HB1, and 16 SNP loci were related to the etiolation of HJYa. 12 CsPPR proteins that may be combined with chloroplast genes were found to use the PPR codes. According to the expression profile and SNP analysis, CsPPR178, CsPPR332, and CsPPR368 might have a close relationship with the albino and etiolation of tea leaves. This study provides comprehensive insights into the CsPPR proteins, which may contribute to further study of the role of CsPPR proteins on the albino and etiolation of tea leaves.