Study on the causes of changes in colour during Hibiscus syriacus flowering based on transcriptome and metabolome analyses.

BACKGROUND: The flower colour of H. syriacus 'Qiansiban' transitions from fuchsia to pink-purple and finally to pale purple, thereby enhancing the ornamental value of the cultivars. However, the molecular mechanism underlying this change in flower colour in H. syriacus has not been elucidated. In this study, the transcriptomic data of H. syriacus 'Qiansiban' at five developmental stages were analysed to investigate the impact of flavonoid components on flower colour variation. Additionally, five cDNA libraries were constructed from H. syriacus 'Qiansiban' during critical blooming stages, and the transcriptomes were sequenced to investigate the molecular mechanisms underlying changes in flower colouration. RESULTS: High-performance liquid chromatography‒mass spectrometry detected five anthocyanins in H. syriacus 'Qiansiban', with malvaccin-3-O-glucoside being the predominant compound in the flowers of H. syriacus at different stages, followed by petunigenin-3-O-glucoside. The levels of these five anthocyanins exhibited gradual declines throughout the flowering process. In terms of the composition and profile of flavonoids and flavonols, a total of seven flavonoids were identified: quercetin-3-glucoside, luteolin-7-O-glucoside, Santianol-7-O-glucoside, kaempferol-O-hexosyl-C-hexarbonoside, apigenin-C-diglucoside, luteolin-3,7-diglucoside, and apigenin-7-O-rutinoside. A total of 2,702 DEGs were identified based on the selected reference genome. Based on the enrichment analysis of differentially expressed genes, we identified 9 structural genes (PAL, CHS, FLS, DRF, ANS, CHI, F3H, F3'5'H, and UFGT) and 7 transcription factors (3 MYB, 4 bHLH) associated with flavonoid biosynthesis. The qRT‒PCR results were in good agreement with the high-throughput sequencing data. CONCLUSION: This study will establish a fundamental basis for elucidating the mechanisms underlying alterations in the flower pigmentation of H. syriacus.

A Breakthrough in Kitavirids: Genetic Variability, Reverse Genetics, Koch's Postulates, and Transmission of Hibiscus Green Spot Virus 2.

Hibiscus green spot virus 2 (HGSV-2), a member of the genus Higrevirus (family Kitaviridae), is a positive-stranded RNA virus associated with leprosis-like symptoms in citrus and green spots on leaves in hibiscus. HGSV-2 has only been reported in Hawaii, and while it is speculated that mites in the genus Brevipalpus might be responsible for its transmission, proper transmission assays have yet to be conducted. This study characterizes additional citrus and hibiscus isolates of HGSV-2 collected from two Hawaiian Islands. We constructed an infectious cDNA clone from a hibiscus isolate of HGSV-2 collected on Oahu and demonstrated its ability to infect several experimental hosts, including Phaseolus vulgaris, Nicotiana tabacum, and N. benthamiana, as well as natural hosts, Citrus reticulata and Hibiscus arnottianus. Bacilliform virions with varied sizes of 33 to 120 nm (length) and 14 to 70 nm (diameter) were observed in partially purified preparations obtained from agroinoculated leaves. Virus progeny from the infectious cDNA clone was found to be infectious after mechanical transmission to N. benthamiana and to cause local lesions. Finally, an isoline colony of the mite Brevipalpus azores had vector competence to transmit a citrus isolate of HGSV-2 collected from Maui to citrus and hibiscus plants, demonstrating the mite-borne nature of HGSV-2. The infectious cDNA clone developed in this study is the first reverse-genetics system for a kitavirid and will be fundamental to better characterize basic biology of HGSV-2 and its interactions with host plants and mite vectors.

Physiological and transcriptional mechanisms associated with cadmium stress tolerance in Hibiscus syriacus L.

BACKGROUND: Cadmium (Cd) pollution of soils is a global concern because its accumulation in plants generates severe growth retardation and health problems. Hibiscus syriacus is an ornamental plant that can tolerate various abiotic stresses, including Cd stress. Therefore, it is proposed as a plant material in Cd-polluted areas. However, the molecular mechanisms of H. syriacus tolerance to Cd are not yet understood. RESULTS: This study investigated the physiological and transcriptional response of "Hongxing", a Cd2+-tolerant H. syriacus variety, grown on a substrate containing higher concentration of Cd (400 mg/kg). The Cd treatment induced only 28% of plant mortality, but a significant decrease in the chlorophyll content was observed. Malondialdehyde content and activity of the antioxidant enzymes catalase, peroxidase, and superoxide dismutase were significantly increased under Cd stress. Transcriptome analysis identified 29,921 differentially expressed genes (DEGs), including 16,729 down-regulated and 13,192 up-regulated genes, under Cd stress. Functional enrichment analyses assigned the DEGs mainly to plant hormone signal transduction, transport, nucleosome and DNA processes, mitogen-activated protein kinase signaling pathway, antioxidant process, fatty acid metabolism, and biosynthesis of secondary metabolites. Many MYB, EP2/ERF, NAC, WRKY family genes, and genes containing metal binding domains were up-regulated, implying that they are essential for the Cd-stress response in H. syriacus. The most induced genes were filtered out, providing valuable resources for future studies. CONCLUSIONS: Our findings provide insights into the molecular responses to Cd stress in H. syriacus. Moreover, this study offers comprehensive and important resources for future studies toward improving the plant Cd tolerance and its valorization in phytoremediation.

Spatiotemporal process of long-distance seed dispersal in a pantropically distributed sea hibiscus group.

Long-distance dispersal (LDD) of seeds plays an essential role in the migration of plants to a new habitat and maintaining gene flow among geographically isolated populations. Pantropical plants with sea-drifted seeds, which have one of the largest distributions in all flowering plants, have achieved their global distribution by LDD. However, the spatiotemporal processes to achieve the wide distribution and the role of LDD in it have not yet been elucidated. In this study, we conducted phylogenomic analyses on the plastome, genome-wide nuclear SNP, and low-copy gene data of Hibiscus tiliaceus and its relatives. The dated phylogeny suggested that global expansion started approximately 4 million years ago (Ma), and species diversification occurred 1 Ma. Plastome phylogeny confirmed the nonmonophyly of the haplotypes in the two widely distributed coastal species, H. tiliaceus and H. pernambucensis. In contrast, genome-wide nuclear SNP phylogenies demonstrated clear genetic segregation among species and/or geographical regions. Ancestral polymorphisms in chloroplast genomes shared among widely distributed species have remained below the range of rapid expansion and speciation of marginal populations. This study demonstrated that the LDD of sea-drifted seeds contributed to the rapid expansion and pantropical distribution of sea hibiscus in the last few million years, and adaptation to local environment or isolation by regional effect after LDD promoted speciation, suppressing gene flow.

Genetic assessment improves conservation efforts for the critically endangered oceanic island endemic Hibiscus liliiflorus.

Hibiscus liliiflorus, endemic to the Indian Ocean island of Rodrigues, is one of the rarest plant species in the world; only 2 wild individuals remain. Previously, when 4 wild individuals remained, the Mauritian Wildlife Foundation (MWF) in Rodrigues propagated cuttings of them in their nursery, then planted seedlings produced in the nursery into 3 outplanted populations on the island. Our goals were to: 1) assess whether all 4 original wild genotypes are represented in the MWF nursery; 2) determine whether ex situ living collections at international botanical gardens maintain unique genotypes of H. liliiflorus; 3) assess whether nursery individuals have crossed or self-fertilized to produce seed and quantify their relative contributions to outplanted populations; and 4) provide recommendations for future conservation actions. We used a 2b-RADseq approach to produce 2,711 genome-wide single nucleotide polymorphisms (SNPs) from 98 samples. Genotype identity analysis, principal component analysis, and model-based clustering in STRUCTURE found 4 genotypes extant in Rodrigues but no unique genotypes in ex situ botanic garden collections. Only 3 genotypes are represented in the MWF nursery; the one remaining genotype is represented by an extant wild individual. Parentage analysis showed that seeds produced in the MWF nursery resulted from both self-fertilization and crossing between genotypes, a result supported by internal relatedness and hybrid index calculations. Each outplanted population is dominated by a subset of parental genotypes, and we propose actions to balance the parental contributions to outplanted populations. Our study highlights how genetic assessments of ex situ conservation projects help conserve critically endangered species.

Genome-Wide Analysis of the MADS-Box Gene Family in Hibiscus syriacus and Their Role in Floral Organ Development.

Hibiscus syriacus belongs to the Malvaceae family, and is a plant with medicinal, edible, and greening values. MADS-box transcription factor is a large family of regulatory factors involved in a variety of biological processes in plants. Here, we performed a genome-wide characterization of MADS-box proteins in H. syriacus and investigated gene structure, phylogenetics, cis-acting elements, three-dimensional structure, gene expression, and protein interaction to identify candidate MADS-box genes that mediate petal developmental regulation in H. syriacus. A total of 163 candidate MADS-box genes were found and classified into type I (Mα, Mß, and Mγ) and type II (MIKC and Mδ). Analysis of cis-acting elements in the promoter region showed that most elements were correlated to plant hormones. The analysis of nine HsMADS expressions of two different H. syriacus cultivars showed that they were differentially expressed between two type flowers. The analysis of protein interaction networks also indicated that MADS proteins played a crucial role in floral organ identification, inflorescence and fruit development, and flowering time. This research is the first to analyze the MADS-box family of H. syriacus and provides an important reference for further study of the biological functions of the MADS-box, especially in flower organ development.

Integrated Methylome and Transcriptome Analysis Provides Insights into the DNA Methylation Underlying the Mechanism of Cytoplasmic Male Sterility in Kenaf (Hibiscus cannabinus L.).

Cytoplasmic male sterility (CMS) is widely exploited in hybrid seed production. Kenaf is an important fiber crop with high heterosis. The molecular mechanism of kenaf CMS remains unclear, particularly in terms of DNA methylation. Here, using the anthers of a kenaf CMS line (P3A) and its maintainer line (P3B), comparative physiological, DNA methylation, and transcriptome analyses were performed. The results showed that P3A had considerably lower levels of IAA, ABA, photosynthetic products and ATP contents than P3B. DNA methylome analysis revealed 650 differentially methylated genes (DMGs) with 313 up- and 337 down methylated, and transcriptome analysis revealed 1788 differentially expressed genes (DEGs) with 558 up- and 1230 downregulated genes in P3A compared with P3B. Moreover, 45 genes were characterized as both DEGs and DMGs, including AUX,CYP, BGL3B, SUS6, AGL30 and MYB21. Many DEGs may be regulated by related DMGs based on methylome and transcriptome studies. These DEGs were involved in carbon metabolism, plant hormone signal transduction, the TCA cycle and the MAPK signaling pathway and were shown to be important for CMS in kenaf. These results provide new insights into the epigenetic mechanism of CMS in kenaf and other crops.

Genome-Wide Identification and Characterization of NAC Family in Hibiscus hamabo Sieb. et Zucc. under Various Abiotic Stresses.

NAC transcription factor is one of the largest plant gene families, participating in the regulation of plant biological and abiotic stresses. In this study, 182 NAC proteins (HhNACs) were identified based on genomic datasets of Hibiscus hamabo Sieb. et Zucc (H. hamabo). These proteins were divided into 19 subfamilies based on their phylogenetic relationship, motif pattern, and gene structure analysis. Expression analysis with RNA-seq revealed that most HhNACs were expressed in response to drought and salt stress. Research of quantitative real-time PCR analysis of nine selected HhNACs supported the transcriptome data's dependability and suggested that HhNAC54 was significantly upregulated under multiple abiotic stresses. Overexpression of HhNAC54 in Arabidopsis thaliana (A. thaliana) significantly increased its tolerance to salt. This study provides a basis for a comprehensive analysis of NAC transcription factor and insight into the abiotic stress response mechanism in H. hamabo.

Nucleo-cytoplasmic interactions affect the 5' terminal transcription of mitochondrial genes between the isonuclear CMS line UG93A and its maintainer line UG93B of kenaf (Hibiscus cannabinus).

Gene expression and translation in plant mitochondria remain poorly understood due to the complicated transcription of its mRNA. In this study, we report the 5' and 3' RNA extremities and promoters of five mitochondrial genes, atp1, atp4, atp6, atp9, and cox3. The results reveal that four genes (atp1, atp4, atp6, and cox3) are transcribed from multiple initiation sites but with a uniform transcript at the 3' end, indicating that heterogeneity of the 5' end is a common feature in the transcription of kenaf mitochondrial genes. Furthermore, we found that the transcription initiation sites of these four genes are significantly different in UG93A, UG93B, and the F1 hybrid. These data indicate that nuclear loci and unknown transcription factors within the mitochondria of different cytoplasmic types may be involved in mitochondrial transcription. Promoter architecture analysis showed that the promoter core sequences are conserved in the kenaf mitochondrial genome but are highly divergent, suggesting that these elements are essential for the promoter activity of mitochondrial genes in kenaf. Our results reveal that the heterogeneity of the 5' end and uniformity at the 3' end are common transcriptional features of mitochondrial genes. These data provide essential information for understanding the transcription of mitochondrial genes in kenaf and can be used as a reference for other plants.

Distribution of FISH oligo-5S rDNA and oligo-(AGGGTTT)3 in Hibiscus mutabilis L.

Hibiscus exhibits high variation in chromosome number both within and among species. The Hibiscus mutabilis L. karyotype was analyzed in detail using fluorescence in situ hybridization (FISH) with oligonucleotide probes for (AG3T3)3 and 5S rDNA, which were tested here for the first time. In total, 90 chromosomes were counted in prometaphase and metaphase, and all exhibited similarly intense (AG3T3)3 signals at both ends. (AG3T3)3 showed little variation and thus did not allow discrimination among H. mutabilis chromosomes, but its location at both ends confirmed the integrity of each chromosome, thus contributing to accurate counting of the numerous, small chromosomes. Oligo-5S rDNA marked the proximal/distal regions of six chromosomes: weak signals on chromosomes 7 and 8, slightly stronger signals on chromosomes 15 and 16, and very strong signals on chromosomes 17 and 18. Therefore, 5S rDNA could assist in chromosome identification in H. mutabilis. Metaphase chromosome lengths ranged from 3.00 to 1.18 µm, indicating small chromosomes. The ratios of longest to shortest chromosome length in prometaphase and metaphase were 2.58 and 2.54, respectively, indicating karyotype asymmetry in H. mutabilis. These results provide an exact chromosome number and a physical map, which will be useful for genome assembly and contribute to molecular cytogenetics in the genus Hibiscus.

A comprehensive integrated transcriptome and metabolome analyses to reveal key genes and essential metabolic pathways involved in CMS in kenaf.

KEY MESSAGE: Numbers of critical genes and pathways were found from the levels of transcriptome and metabolome, which were useful information for understanding of kenaf CMS mechanism. Cytoplasmic male sterility (CMS) is a maternally inherited trait in higher plants that leads to the inability to produce or release functional pollen. However, there is lack of comprehensive studies to reveal the molecular basis of CMS occurrence in kenaf. Herein, we performed transcriptome and UPLC-MS-based metabolome analyses in the anthers of a CMS (UG93A) and its maintainer (UG93B) to sort out essential genes and metabolites responding to CMS in kenaf. Transcriptome characterized 7769 differentially expressed genes (DEGs) between these two materials, and pathway enrichment analysis indicated that these DEGs were involved mainly in pentose and glucuronate interconversions, starch and sucrose metabolism, taurine and hypotaurine metabolism. In the metabolome assay, a total of 116 significantly different metabolites (SDMs) were identified between the CMS and its maintainer line, and these SDMs were involved in eight KEGG pathways, including flavone and flavonol biosynthesis, glycerophospholipid metabolism, flavonoid biosynthesis, glycosylphosphatidylinositol-anchor biosynthesi. Integrated analyses of transcriptome and metabolome showed that 50 genes had strong correlation coefficient values (R2 > 0.9) with ten metabolites enriched in six pathways; notably, most genes and metabolites of flavonoid biosynthesis pathways and flavone and flavonol biosynthesis pathways involved in flavonoids biosynthetic pathways were downregulated in CMS compared to those in maintainer. Taken together, the decreased accumulation of flavonoids resulted from the compromised biosynthesis pathways coupled with energy deficiency in the anthers may contribute largely to CMS in UG93A of kenaf.

Chloroplast genome of Hibiscus rosa-sinensis (Malvaceae): Comparative analyses and identification of mutational hotspots.

Previous studies to resolve phylogenetic and taxonomic discrepancies of Hibiscus remained inconclusive. Here, we report chloroplast genome sequence of Hibiscus rosa-sinensis. Hibiscus rosa-sinensis chloroplast genome was 160,951 bp, comprising of large single copy (89,509 bp) and small single copy (20,246 bp) regions, separated by IRa and IRb (25,598 bp each). The genome contained 130 genes including 85 protein-coding genes, 37 transfer RNAs and 8 ribosomal RNAs. Comparative analyses of chloroplast genomes revealed similar structure among 12 species within family Malvaceae. Evolutionary rates of 77 protein-coding genes showed 95% similarities. Analyses of codon usage, amino acid frequency, putative RNA editing sites, and repeats showed a great extent of similarities between Hibiscus rosa-sinensis and Hibiscus syriacus. We identified 30 mutational hotpots including psbZ-trnG, trnK-rps16, trnD-trnY, trnW-trnP, rpl33-rps18, petG-trnW, trnS-trnG, trnH-psbA, atpB-rbcL, and rpl32-trnL that might be used as polymorphic and robust markers to resolve phylogenetic discrepancies in genus Hibiscus.

Transcriptome Analysis of Salt Stress in Hibiscus hamabo Sieb. et Zucc Based on Pacbio Full-Length Transcriptome Sequencing.

Hibiscus hamabo Sieb. et Zucc is an important semi-mangrove plant with great morphological features and strong salt resistance. In this study, by combining single molecule real time and next-generation sequencing technologies, we explored the transcriptomic changes in the roots of salt stressed H. hamabo. A total of 94,562 unigenes were obtained by clustering the same isoforms using the PacBio RSII platform, and 2269 differentially expressed genes were obtained under salt stress using the Illumina platform. There were 519 differentially expressed genes co-expressed at each treatment time point under salt stress, and these genes were found to be enriched in ion signal transduction and plant hormone signal transduction. We used Arabidopsis thaliana (L.) Heynh. transformation to confirm the function of the HhWRKY79 gene and discovered that overexpression enhanced salt tolerance. The full-length transcripts generated in this study provide a full characterization of the transcriptome of H. hamabo and may be useful in mining new salt stress-related genes specific to this species, while facilitating the understanding of the salt tolerance mechanisms.

Genome-wide Analysis of Basic Helix-Loop-Helix Family Genes and Expression Analysis in Response to Drought and Salt Stresses in Hibiscus hamabo Sieb. et Zucc.

The basic helix-loop-helix (bHLH) family of transcription factors is one of the most significant and biggest in plants. It is involved in the regulation of both growth and development, as well as stress response. Numerous members of the bHLH family have been found and characterized in woody plants in recent years. However, no systematic study of the bHLH gene family has been published for Hibiscus hamabo Sieb. et Zucc. In this research, we identified 162 bHLH proteins (HhbHLHs) from the genomic and transcriptomic datasets of H. hamabo, which were phylogenetically divided into 19 subfamilies. According to a gene structural study, the number of exon-introns in HhbHLHs varied between zero and seventeen. MEME research revealed that the majority of HhbHLH proteins contained three conserved motifs, 1, 4, and 5. The examination of promoter cis-elements revealed that the majority of HhbHLH genes had several cis-elements involved in plant growth and development and abiotic stress responses. In addition, the overexpression of HhbHLH2 increased salt and drought stress tolerance in Arabidopsis.

Reactive Oxygen Species Accumulation Strongly Allied with Genetic Male Sterility Convertible to Cytoplasmic Male Sterility in Kenaf.

Male sterility (MS) plays a key role in the hybrid breed production of plants. Researchers have focused on the association between genetic male sterility (GMS) and cytoplasmic male sterility (CMS) in kenaf. In this study, P9BS (a natural GMS mutant of the kenaf line P9B) and male plants of P9B were used as parents in multiple backcross generations to produce P9SA, a CMS line with stable sterility, to explore the molecular mechanisms of the association between GMS and CMS. The anthers of the maintainer (P9B), GMS (P9BS), and CMS (P9SA) lines were compared through phenotypic, cell morphological, physiological, biochemical observations, and transcriptome analysis. Premature degradation of the tapetum was observed at the mononuclear stage in P9BS and P9SA, which also had lower activity of reactive oxygen species (ROS) scavenging enzymes compared with P9B. Many coexpressed differentially expressed genes were related to ROS balance, including ATP synthase, electron chain transfer, and ROS scavenging processes were upregulated in P9B. CMS plants had a higher ROS accumulation than GMS plants. The MDA content in P9SA was 3.2 times that of P9BS, and therefore, a higher degree of abortion occurred in P9SA, which may indicate that the conversion between CMS and GMS is related to intracellular ROS accumulation. Our study adds new insights into the natural transformation of GMS and CMS in plants in general and kenaf in particular.

The genome of kenaf (Hibiscus cannabinus L.) provides insights into bast fibre and leaf shape biogenesis.

Kenaf is an annual crop that is widely cultivated as a source of bast (phloem) fibres, the phytoremediation of heavy metal-contaminated farmlands and textile-relevant compounds. Leaf shape played a unique role in kenaf improvement, due to the inheritance as a single locus and the association with fibre development in typical lobed-leaf varieties. Here we report a high-quality genome assembly and annotation for var. 'Fuhong 952' with 1078 Mbp genome and 66 004 protein-coding genes integrating single-molecule real-time sequencing, a high-density genetic map and high-throughput chromosome conformation capture techniques. Gene mapping assists the identification of a homeobox transcription factor LATE MERISTEM IDENTITY 1 (HcLMI1) gene controlling lobed-leaf. Virus-induced gene silencing (VIGS) of HcLMI1 in a lobed-leaf variety was critical to induce round (entire)-like leaf formation. Candidate genes involved in cell wall formation were found in quantitative trait loci (QTL) for fibre yield and quality-related traits. Comparative genomic and transcriptome analyses revealed key genes involved in bast fibre formation, among which there are twice as many cellulose synthase A (CesA) genes due to a recent whole-genome duplication after divergence from Gossypium. Population genomic analysis showed two recent population bottlenecks in kenaf, suggesting domestication and improvement process have led to an increase in fibre biogenesis and yield. This chromosome-scale genome provides an important framework and toolkit for sequence-directed genetic improvement of fibre crops.

Comparative profile analysis reveals differentially expressed microRNAs regulate anther and pollen development in kenaf cytoplasmic male sterility line.

Cytoplasmic male sterility (CMS) is advantageous in extensive crop breeding and represents a perfect model for understanding anther and pollen development research. MicroRNAs (miRNAs) play key roles in regulating various biological processes. However, the miRNA-mediated regulatory network in kenaf CMS occurrence remains largely unknown. In the present study, a comparative deep sequencing approach was used to investigate the miRNAs and their roles in regulating anther and pollen development during CMS occurrence. We identified 283 known and 46 new candidate miRNAs in kenaf anther. A total of 67 differentially expressed miRNAs (DEMs) were discovered between CMS and its maintainer line. Among them, 40 and 27 miRNAs were up- and downregulated, respectively. These 67 DEMs were predicted to target 189 genes. Validation of DEMs and putative target genes were confirmed by using real-time quantitative PCR. In addition, a potential miRNA-mediated regulatory network, which mainly involves the auxin signaling pathway, signal transduction, glycolysis and energy metabolism, gene expression, transmembrane transport, protein modification and metabolism, and floral development, that mediates anther development during CMS occurrence was proposed. Taken together, our findings provide a better understanding of the molecular mechanism of miRNA regulation in pollen development and CMS occurrence in kenaf.

Comparative acetylomic analysis reveals differentially acetylated proteins regulating anther and pollen development in kenaf cytoplasmic male sterility line.

Cytoplasmic male sterility (CMS) is widely used in plant breeding and represents a perfect model to understand cyto-nuclear interactions and pollen development research. Lysine acetylation in proteins is a dynamic and reversible posttranslational modification (PTM) that plays an important roles in diverse cell processes and signaling. However, studies addressing acetylation PTM regarding to anther and pollen development in CMS background are largely lacking. To reveal the possible mechanism of kenaf (Hibiscus cannabinus L.) CMS and pollen development, we performed a label-free-based comparative acetylome analysis in kenaf anther of a CMS line and wild-type (Wt). Using whole transcriptome unigenes of kenaf as the reference genome, we identified a total of 1204 Kac (lysin acetylation) sites on 1110 peptides corresponding to 672 unique proteins. Futher analysis showed 56 out of 672 proteins were differentially acetylated between CMS and Wt line, with 13 and 43 of those characterized up- and downregulated, respectively. Thirty-eight and 82 proteins were detected distinctively acetylated in CMS and Wt lines, respectively. And evaluation of the acetylomic and proteomic results indicated that the most significantly acetylated proteins were not associated with abundant changes at the protein level. Bioinformatics analysis demonstrated that many of these proteins were involved in various biological processes which may play key roles in pollen development, inculding tricarboxylic acid (TCA) cycle and energy metabolism, protein folding, protein metabolism, cell signaling, gene expression regulation. Taken together, our results provide insight into the CMS molecular mechanism and pollen development in kenaf from a protein acetylation perspective.

Comparative Transcriptome Analysis Reveals the Cause for Accumulation of Reactive Oxygen Species During Pollen Abortion in Cytoplasmic Male-Sterile Kenaf Line 722HA.

Cytoplasmic male sterility (CMS) is a maternally inherited trait used for hybrid production in plants, a novel kenaf CMS line 722HA was derived from the thermo-sensitive male-sterile mutant 'HMS' by recurrent backcrossing with 722HB. The line 722HA has great potential for hybrid breeding in kenaf. However, the underlying molecular mechanism that controls pollen abortion in 722HA remains unclear, thus limiting the full utilization of this line. To understand the possible mechanism governing pollen abortion in 722HA, cytological, transcriptomic, and biochemical analyses were carried out to compare the CMS line 722HA and its maintainer line 722HB. Cytological observations of the microspore development revealed premature degradation of the tapetum at the mononuclear stage, which resulted in pollen dysfunction. The k-means clustering analysis of differentially expressed genes (DEGs) revealed that these genes are related to processes associated with the accumulation of reactive oxygen species (ROS), including electron transport chain, F1F0-ATPase proton transport, positive regulation of superoxide dismutase (SOD), hydrogen peroxide catabolic, and oxidation-reduction. Biochemical analysis indicated that ROS-scavenging capability was lower in 722HA than in 722HB, resulting in an accumulation of excess ROS, which is consistent with the transcriptome results. Taken together, these results demonstrate that excessive ROS accumulation may affect the normal development of microspores. Our study provides new insight into the molecular mechanism of pollen abortion in 722HA and will promote further studies of kenaf hybrids.

Comparative phosphoproteomic analysis reveals differentially phosphorylated proteins regulate anther and pollen development in kenaf cytoplasmic male sterility line.

Cytoplasmic male sterility (CMS) is widely used in plant breeding and represents a perfect model to understand cyto-nuclear interactions and pollen development research. Protein phosphorylation is ubiquitous and is involved in the regulation of diverse cellular processes. To reveal the possible mechanism of CMS and pollen development in kenaf, we performed an iTRAQ-based comparative phosphoproteome analysis in the anthers of a CMS line and wild-type plant (Wt). Whole transcriptome unigenes of kenaf as the reference genome, we identified a total of 3045 phosphorylated sites on 1640 peptides corresponding to 974 unique proteins. 292 of the peptides which corresponding to 247 unique proteins were differentially phosphorylated (fold change ≥ 1.20 with P value< 0.05) between these two materials. 113 and 134 proteins were characterized as up-regulated or down-regulated phosphorylated, respectively. An evaluation of the phosphoproteome and proteomic results indicated that the most significantly phosphorylated proteins were not associated with abundant changes at the protein level. Bioinformatics analysis demonstrated that many of these proteins were involved in various biological processes which may play key roles in pollen development, including carbohydrate metabolism, energy metabolism, transport, gene expression regulation, signal transduction, and cell cycle control. Our results provide insight into the CMS mechanism and pollen development in kenaf from a protein phosphorylation perspective.