Cotton Fiber Development Requires the Pentatricopeptide Repeat Protein GhIm for Splicing of Mitochondrial nad7 mRNA

Genome-Wide Identification of Pentatricopeptide Repeat (PPR) Gene Family and Multi-Omics Analysis Provide New Insights Into the Albinism Mechanism of Kandelia obovata Propagule Leaves

Pentatricopeptide repeat (PPR) gene family constitutes one of the largest gene families in plants, which mainly participate in RNA editing and RNA splicing of organellar RNAs, thereby affecting the organellar development. Recently, some evidence elucidated the important roles of PPR proteins in the albino process of plant leaves. However, the functions of PPR genes in the woody mangrove species have not been investigated. In this study, using a typical true mangrove Kandelia obovata, we systematically identified 298 PPR genes and characterized their general features and physicochemical properties, including evolutionary relationships, the subcellular localization, PPR motif type, the number of introns and PPR motifs, and isoelectric point, and so forth. Furthermore, we combined genome-wide association studies (GWAS) and transcriptome analysis to identify the genetic architecture and potential PPR genes associated with propagule leaves colour variations of K. obovata. As a result, we prioritized 16 PPR genes related to the albino phenotype using different strategies, including differentially expressed genes analysis and genetic diversity analysis. Further analysis discovered two genes of interest, namely Maker00002998 (PLS-type) and Maker00003187 (P-type), which were differentially expressed genes and causal genes detected by GWAS analysis. Moreover, we successfully predicted downstream target chloroplast genes (rps14, rpoC1 and rpoC2) bound by Maker00002998 PPR proteins. The experimental verification of RNA editing sites of rps14, rpoC1, and rpoC2 in our previous study and the verification of interaction between Maker00002998 and rps14 transcript using in vitro RNA pull-down assays revealed that Maker00002998 PPR protein might be involved in the post-transcriptional process of chloroplast genes. Our result provides new insights into the roles of PPR genes in the albinism mechanism of K. obovata propagule leaves.

Comparative Transcriptomic Analysis of Gossypium hirsutum Fiber Development in Mutant Materials (xin w 139) Provides New Insights into Cotton Fiber Development

Cotton is the most widely planted fiber crop in the world, and improving cotton fiber quality has long been a research hotspot. The development of cotton fibers is a complex process that includes four consecutive and overlapping stages, and although many studies on cotton fiber development have been reported, most of the studies have been based on cultivars that are promoted in production or based on lines that are used in breeding. Here, we report a phenotypic evaluation of Gossypium hirsutum based on immature fiber mutant (xin w 139) and wild-type (Xin W 139) lines and a comparative transcriptomic study at seven time points during fiber development. The results of the two-year study showed that the fiber length, fiber strength, single-boll weight and lint percentage of xin w 139 were significantly lower than those of Xin W 139, and there were no significant differences in the other traits. Principal component analysis (PCA) and cluster analysis of the RNA-sequencing (RNA-seq) data revealed that these seven time points could be clearly divided into three different groups corresponding to the initiation, elongation and secondary cell wall (SCW) synthesis stages of fiber development, and the differences in fiber development between the two lines were mainly due to developmental differences after twenty days post anthesis (DPA). Differential expression analysis revealed a total of 5131 unique differentially expressed genes (DEGs), including 290 transcription factors (TFs), between the 2 lines. These DEGs were divided into five clusters. Each cluster functional category was annotated based on the KEGG database, and different clusters could describe different stages of fiber development. In addition, we constructed a gene regulatory network by weighted correlation network analysis (WGCNA) and identified 15 key genes that determined the differences in fiber development between the 2 lines. We also screened seven candidate genes related to cotton fiber development through comparative sequence analysis and qRT-PCR; these genes included three TFs (GH_A08G1821 (bHLH), GH_D05G3074 (Dof), and GH_D13G0161 (C3H)). These results provide a theoretical basis for obtaining an in-depth understanding of the molecular mechanism of cotton fiber development and provide new genetic resources for cotton fiber research.

Attenuated Total Reflection Fourier Transform Infrared Spectroscopy with Soft Independent Modeling of Class Analogy-Principal Component Analysis for Classifying Cotton Fiber Maturity Phenotypes of Cotton Population Composed of Various Genotypes

Maturity is a major fiber trait that affects the processing and performance of cotton fiber. Rapid and accurate identification of fiber maturity phenotypes and genotypes is of importance to breeders. Previous studies showed that either conventional fiber measurements or attenuated total reflection Fourier transform infrared spectroscopy (ATR FT-IR) analysis discriminated the immature fiber (im) phenotype from the wild type (WT) mature fiber phenotype in a segregating F2 population from a cross between two upland cotton lines differing in fiber maturity. However, both conventional fiber property measurement methods and FT-IR analyses with current algorithms could not detect the subtle differences among the WT fibers composed of two different genotypes, WT homozygosity (WT-homo) and WT heterozygosity (WT-hetero). This research explored the FT-IR method, in combination with soft independent modeling of class analogy of principal component analysis (SIMCA-PCA), for the discrimination of WT fiber phenotypes consisting of two different genotypes (WT-homo and WT-hetero). The new approach enabled the detection of IR spectral intensity differences between WT-homo and WT-hetero fibers. Successful classification originated from a distinctive spectral difference in the low-wavenumber region (<700 cm-1) between WT-hetero fibers and WT-homo fibers. This observation emphasized that ATR FT-IR with a SIMCA-PCA approach would be a sensitive tool for classifying the WT fibers demonstrating minor phenotypic differences. The improved sensitivity of the infrared method may provide a way of dissecting genotype-phenotype interactions of cotton fibers rapidly and efficiently.

Genetic Analysis of Mutagenesis That Induces the Photoperiod Insensitivity of Wild Cotton Gossypium hirsutum Subsp. purpurascens

Cotton genus Gossypium L., especially its wild species, is rich in genetic diversity. However, this valuable genetic resource is barely used in cotton breeding programs. In part, due to photoperiod sensitivities, the genetic diversity of Gossypium remains largely untapped. Herein, we present a genetic analysis of morphological, cytological, and genomic changes from radiation-mediated mutagenesis that induced plant photoperiod insensitivity in the wild cotton of Gossypium hirsutum. Several morphological and agronomical traits were found to be highly inheritable using the progeny between the wild-type G. hirsutum subsp. purpurascens (El-Salvador) and its mutant line (Kupaysin). An analysis of pollen mother cells (PMCs) revealed quadrivalents that had an open ring shape and an adjoining type of divergence of chromosomes from translocation complexes. Using 336 SSR markers and 157 F₂ progenies that were grown with parental genotypes and F₁ hybrids in long day and short night conditions, five quantitative trait loci (QTLs) associated with cotton flowering were located on chromosomes At-05, At-11, and Dt-07. Nineteen candidate genes related to the flowering traits were suggested through molecular and in silico analysis. The DNA markers associated with the candidate genes, upon future functional analysis, would provide useful tools in marker-assisted selection (MAS) in cotton breeding programs for early flowering and maturity.

Comparison of Mitochondrial Genomes between a Cytoplasmic Male-Sterile Line and Its Restorer Line for Identifying Candidate CMS Genes in Gossypium hirsutum

As the core of heterosis utilization, cytoplasmic male sterility (CMS) has been widely used in hybrid seed production. Previous studies have shown that CMS is always closely related to the altered programming of mitochondrial genes. To explore candidate CMS genes in cotton (Gossypium hirsutum), sequencing and de novo assembly were performed on the mitochondrial genome of the G. hirsutum CMS line SI3A, with G. harknessii CMS-D2 cytoplasm, and the corresponding G. hirsutum restorer line 0-613-2R. Remarkable variations in genome structure and gene transcripts were detected. The mitochondrial genome of SI3A has three circle molecules, including one main circle and two sub-circles, while 0-613-2R only has one. RNA-seq and RT-qPCR analysis proved that orf606a and orf109a, which have a chimeric structure and transmembrane domain, were highly expressed in abortive anthers of SI3A. In addition, comparative analysis of RNA-seq and full-length transcripts revealed the complex I gene nad4 to be expressed at a lower level in SI3A than in its restorer and that it featured an intron retention splicing pattern. These two novel chimeric ORFs and nad4 are potential candidates that confer CMS character in SI3A. This study provides new insight into the molecular basis of the nuclear–cytoplasmic interaction mechanism, and that putative CMS genes might be important sources for future precise design cross-breeding of cotton.

Evidence for thermosensitivity of the cotton (Gossypium hirsutum L.) immature fiber (im) mutant via hypersensitive stomatal activity

Thickness of cotton fiber, referred to as fiber maturity, is a key determinant of fiber quality, lint yield, and textile performance. The cotton immature fiber (im) mutant has been used to study fiber maturity since its fiber is thinner than the wild type near isogeneic line (NIL), Texas Marker-1 (TM-1). The im phenotype is caused by a single recessive mutation of a pentatricopeptide repeat (PPR) gene that reduces the activity of mitochondrial complex I and up-regulates stress responsive genes. However, the mechanisms altering the stress responses in im mutant are not well understood. Thus, we characterized growth and gas exchange in im and TM-1 under no stress and also investigated their stress responses by comparing gas exchange and transcriptomic profiles under high temperature. Phenotypic differences were detected between the NILs in non-fiber tissues although less pronounced than the variation in fibers. At near optimum temperature (28±3°C), im maintained the same photosynthetic performance as TM-1 by means of greater stomatal conductance. In contrast, under high temperature stress (>34°C), im leaves reduced photosynthesis by decreasing the stomatal conductance disproportionately more than TM-1. Transcriptomic analyses showed that the genes involved in heat stress responses were differentially expressed between the NIL leaves. These results indicate that the im mutant previously reported to have low activity of mitochondrial complex I displays increased thermosensitivity by impacting stomatal conductance. They also support a notion that mitochondrial complex I activity is required for maintenance of optimal photosynthetic performance and acclimation of plants to high temperature stress. These findings may be useful in the future efforts to understand how physiological mechanisms play a role in determining cotton fiber maturity and may influence stress responses in other crops.

Research Progress of PPR Proteins in RNA Editing, Stress Response, Plant Growth and Development

RNA editing is a posttranscriptional phenomenon that includes gene processing and modification at specific nucleotide sites. RNA editing mainly occurs in the genomes of mitochondria and chloroplasts in higher plants. In recent years, pentatricopeptide repeat (PPR) proteins, which may act as trans-acting factors of RNA editing have been identified, and the study of PPR proteins has become a research focus in molecular biology. The molecular functions of these proteins and their physiological roles throughout plant growth and development are widely studied. In this minireview, we summarize the current knowledge of the PPR family, hoping to provide some theoretical reference for future research and applications.