Exploration of miRNAs and target genes of cytoplasmic male sterility line in cotton during flower bud development

The roles of non-coding RNAs in male reproductive development and abiotic stress responses during this unique process in flowering plants

Successful male reproductive development is the guarantee for sexual reproduction of flowering plants. Male reproductive development is a complicated and multi-stage process that integrates physiological processes and adaptation and tolerance to a myriad of environmental stresses. This well-coordinated process is governed by genetic and epigenetic machineries. Non-coding RNAs (ncRNAs) play pleiotropic roles in the plant growth and development. The identification, characterization and functional analysis of ncRNAs and their target genes have opened a new avenue for comprehensively revealing the regulatory network of male reproductive development and its response to environmental stresses in plants. This review briefly addresses the types, origin, biogenesis and mechanisms of ncRNAs in plants, highlights important updates on the roles of ncRNAs in regulating male reproductive development and emphasizes the contribution of ncRNAs, especially miRNAs and lncRNAs, in responses to abiotic stresses during this unique process in flowering plants.

The applications of CRISPR/Cas-mediated microRNA and lncRNA editing in plant biology: shaping the future of plant non-coding RNA research

MAIN CONCLUSION

CRISPR/Cas technology has greatly facilitated plant non-coding RNA (ncRNA) biology research, establishing itself as a promising tool for ncRNA functional characterization and ncRNA-mediated plant improvement. Throughout the last decade, the promising genome editing tool clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated proteins (Cas; CRISPR/Cas) has allowed unprecedented advances in the field of plant functional genomics and crop improvement. Even though CRISPR/Cas-mediated genome editing system has been widely used to elucidate the biological significance of a number of plant protein-coding genes, this technology has been barely applied in the functional analysis of those non-coding RNAs (ncRNAs) that modulate gene expression, such as microRNAs (miRNAs) and long non-coding RNAs (lncRNAs). Nevertheless, compelling findings indicate that CRISPR/Cas-based ncRNA editing has remarkable potential for deciphering the biological roles of ncRNAs in plants, as well as for plant breeding. For instance, it has been demonstrated that CRISPR/Cas tool could overcome the challenges associated with other approaches employed in functional genomic studies (e.g., incomplete knockdown and off-target activity). Thus, in this review article, we discuss the current status and progress of CRISPR/Cas-mediated ncRNA editing in plant science in order to provide novel prospects for further assessment and validation of the biological activities of plant ncRNAs and to enhance the development of ncRNA-centered protocols for crop improvement.

Insights to Gossypium defense response against Verticillium dahliae: the Cotton Cancer

The soil-borne pathogen Verticillium dahliae, also referred as "The Cotton Cancer," is responsible for causing Verticillium wilt in cotton crops, a destructive disease with a global impact. To infect cotton plants, the pathogen employs multiple virulence mechanisms such as releasing enzymes that degrade cell walls, activating genes that contribute to virulence, and using protein effectors. Conversely, cotton plants have developed numerous defense mechanisms to combat the impact of V. dahliae. These include strengthening the cell wall by producing lignin and depositing callose, discharging reactive oxygen species, and amassing hormones related to defense. Despite the efforts to develop resistant cultivars, there is still no permanent solution to Verticillium wilt due to a limited understanding of the underlying molecular mechanisms that drive both resistance and pathogenesis is currently prevalent. To address this challenge, cutting-edge technologies such as clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9), host-induced gene silencing (HIGS), and gene delivery via nano-carriers could be employed as effective alternatives to control the disease. This article intends to present an overview of V. dahliae virulence mechanisms and discuss the different cotton defense mechanisms against Verticillium wilt, including morphophysiological and biochemical responses and signaling pathways including jasmonic acid (JA), salicylic acid (SA), ethylene (ET), and strigolactones (SLs). Additionally, the article highlights the significance of microRNAs (miRNAs), circular RNAs (circRNAs), and long non-coding RNAs (lncRNAs) in gene expression regulation, as well as the different methods employed to identify and functionally validate genes to achieve resistance against this disease. Gaining a more profound understanding of these mechanisms could potentially result in the creation of more efficient strategies for combating Verticillium wilt in cotton crops.

Identification of fertility restoration candidate genes from a restorer line R186 for Gossypium harknessii cytoplasmic male sterile cotton

BACKGROUND

The utilization of heterosis based on three-line system is an effective strategy in crop breeding. However, cloning and mechanism elucidation of restorer genes for cytoplasmic male sterility (CMS) in upland cotton have yet been realized.

RESULTS

This research is based on CMS line 2074A with the cytoplasm from Gossypium harknessii (D_2-2) and restorer line R186. The offspring of 2074A × R186 were used to conduct genetic analysis. The fertility mechanism of 2074A can be speculated to be governed by multiple genes, since neither the single gene model nor the double genes model could be used. The bulked segregant analysis (BSA) for (2074A × R186) F₂ determined the genetic interval of restorer genes on a region of 4.30 Mb on chromosome D05 that contains 77 annotated genes. Four genes were identified as candidates for fertility restoration using the RNA-seq data of 2074A, 2074B, and R186. There are a number of large effect variants in the four genes between 2074A and R186 that could cause amino acid changes. Evolutionary analysis and identity analysis revealed that GH_D05G3183, GH_D05G3384, and GH_D05G3490 have high identity with their homologs in D_2-2, respectively. Tissue differential expression analysis revealed that the genes GH_D05G3183, GH_D05G3384, and GH_D05G3490 were highly expressed in the buds of the line R186. The predicted results demonstrated that GH_D05G3183, GH_D05G3384 and GH_D05G3490 might interact with GH_A02G1295 to regulate orf610a in mitochondria.

CONCLUSION

Our study uncovered candidate genes for fertility restoration in the restorer line R186 and predicted the possible mechanism for restoring the male fertility in 2074A. This research provided valuable insight into the nucleoplasmic interactions.

Whole-transcriptome analyses identify key differentially expressed mRNAs, lncRNAs, and miRNAs associated with male sterility in watermelon

Male sterility is a valuable trait for watermelon breeding, as watermelon hybrids exhibit obvious heterosis. However, the underlying regulatory mechanism is still largely unknown, especially regarding the related non-coding genes. In the present study, approximately 1035 differentially expressed genes (DEGs), as well as 80 DE-lncRNAs and 10 DE-miRNAs, were identified, with the overwhelming majority down-regulated in male-sterile floral buds. Enrichment analyses revealed that the general phenylpropanoid pathway as well as its related metabolisms was predicted to be altered in a mutant compared to its fertile progenitor. Meanwhile, the conserved genetic pathway DYT1-TDF1-AMS-MS188-MS1, as well as the causal gene ClAMT1 for the male-sterile mutant Se18, was substantially disrupted during male reproductive development. In addition, some targets of the key regulators AMS and MS188 in tapetum development were also down-regulated at a transcriptional level, such as ABCG26 (Cla004479), ACOS5 (Cla022956), CYP703A2 (Cla021151), PKSA (Cla021099), and TKPR1 (Cla002563). Considering lncRNAs may act as functional endogenous target mimics of miRNAs, competitive endogenous RNA networks were subsequently constructed, with the most complex one containing three DE-miRNAs, two DE-lncRNAs, and 21 DEGs. Collectively, these findings not only contribute to a better understanding of genetic regulatory networks underlying male sterility in watermelon, but also provide valuable candidates for future research.

Plant non-coding RNAs function in pollen development and male sterility

Male sterility is classified as either cytoplasmic male sterility (CMS) or genic male sterility (GMS). Generally, CMS involves mitochondrial genomes interacting with the nuclear genome, while GMS is caused by nuclear genes alone. Male sterility is regulated by multilevel mechanisms in which non-coding RNAs (ncRNAs), including microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and phased small interfering RNAs (phasiRNAs), which have been proven to be critical elements. The development of high-throughput sequencing technology offers new opportunities to evaluate the genetic mechanism of ncRNAs in plant male sterility. In this review, we summarize the critical ncRNAs that regulate gene expression in ways dependent on or independent of hormones, which involve the differentiation of the stamen primordia, degradation of the tapetum, formation of microspores, and the release of pollen. In addition, the key mechanisms of the miRNA-lncRNA-mRNA interaction networks mediating male sterility in plants are elaborated. We present a different perspective on exploring the ncRNA-mediated regulatory pathways that control CMS in plants and create male-sterile lines through hormones or genome editing. A refined understanding of the ncRNA regulatory mechanisms in plant male sterility for the development of new sterile lines would be conducive to improve hybridization breeding.

Homology-based identification of candidate genes for male sterility editing in upland cotton (Gossypium hirsutum L.)

Upland cotton (Gossypium hirsutum L.) accounts for more than 90% of the world's cotton production, providing natural material for the textile and oilseed industries worldwide. One strategy for improving upland cotton yields is through increased adoption of hybrids; however, emasculation of cotton flowers is incredibly time-consuming and genetic sources of cotton male sterility are limited. Here we review the known biochemical modes of plant nuclear male sterility (NMS), often known as plant genetic male sterility (GMS), and characterized them into four groups: transcriptional regulation, splicing, fatty acid transport and processing, and sugar transport and processing. We have explored protein sequence homology from 30 GMS genes of three monocots (maize, rice, and wheat) and three dicots (Arabidopsis, soybean, and tomato). We have analyzed evolutionary relationships between monocot and dicot GMS genes to describe the relative similarity and relatedness of these genes identified. Five were lowly conserved to their source species, four unique to monocots, five unique to dicots, 14 highly conserved among all species, and two in the other category. Using this source, we have identified 23 potential candidate genes within the upland cotton genome for the development of new male sterile germplasm to be used in hybrid cotton breeding. Combining homology-based studies with genome editing may allow for the discovery and validation of GMS genes that previously had no diversity observed in cotton and may allow for development of a desirable male sterile mutant to be used in hybrid cotton production.

Integrated Analysis of Small RNA, Transcriptome, and Degradome Sequencing Reveals the MiR156, MiR5488 and MiR399 are Involved in the Regulation of Male Sterility in PTGMS Rice

A photoperiod- and thermo-sensitive genic male sterile (PTGMS) line is the basic material for two-hybrid rice and is an important genetic breeding resource. Peiai64S (PA64S) is an important germplasm resource of PTGMS rice, and it has been applied to two-line hybrid rice systems in China. Pollen fertility in PA64S is regulated by the temperature and photoperiod, but the mechanism of the fertility transition is unclear. In this study, we obtained the male fertile plant PA64S(F) and the male sterile plant PA64S(S) by controlling different temperatures under long light conditions and used the male fertile and sterile plants to investigate the role of microRNAs (miRNAs) in regulating male fertility in rice. We performed the small RNA library sequencing of anthers from PA64S(S) and PA64S(F). A total of 196 miRNAs were identified-166 known miRNAs among 27 miRNA families and 30 novel miRNAs. In the transcriptome analysis, the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of differentially expressed genes revealed significant enrichment in the synthesis and metabolism of fatty acids and some secondary metabolism pathways such as fatty acid metabolism and phenylalanine metabolism. With a comprehensive analysis of miRNA, transcriptome, and degradome sequencing, we identified that 13 pairs of miRNA/target genes regulated male fertility in rice by responding to temperature change, among which the miR156, miR5488, and miR399 affect the male fertility of PA64S by influencing SPLs, the lignin synthesis of anther walls, and the flavonoid metabolism pathway. The results provide a new understanding of PTGMS rice, which will help us better understand the potential regulatory mechanisms of male sterility in the future.

Molecular characterization of miRNA genes and their expression in Dimocarpus longan Lour

A genome-wide analysis of longan miRNA genes was conducted, and full-length pri-miRNA transcripts were cloned. Bioinformatics and expression analyses contributed to the functional characterization of longan miRNA genes. MicroRNAs are important for the post-transcriptional regulation of target genes. However, little is known about the transcription and regulation of miRNA genes in longan (Dimocarpus longan Lour.). In this study, 80 miRNA precursors (pre-miRNA) were predicted, and their secondary structure, size, conservation, and diversity were analyzed. Furthermore, the full-length cDNA sequences of 13 longan primary miRNAs (pri-miRNAs) were amplified by RLM-RACE and SMART-RACE and analyzed, which revealed that longan pri-miRNA transcripts have multiple transcription start sites (TSSs) and the downstream pre-miRNAs are polymorphic. Accordingly, the longan pri-miRNAs and protein-encoding genes may have similar transcriptional specificities. An analysis of the longan miRNA gene promoter elements indicated that the three most abundant cis-acting elements were light-responsive, stress-responsive, and hormone-responsive elements. A quantitative real-time PCR assay elucidated the potential spatial and temporal expression patterns of longan pre-miRNAs during the early stages of somatic embryogenesis (SE) and in different longan organs/tissues. This is the first report regarding the molecular characterization of miRNA genes and their expression profiles in longan. The generated data may serve as a foundation for future research aimed at clarifying the longan miRNA gene functions.

Embryogenic Calli Induction and Salt Stress Response Revealed by RNA-Seq in Diploid Wild Species Gossypium sturtianum and Gossypium raimondii

Wild cotton species can contribute to a valuable gene pool for genetic improvement, such as genes related to salt tolerance. However, reproductive isolation of different species poses an obstacle to produce hybrids through conventional breeding. Protoplast fusion technology for somatic cell hybridization provides an opportunity for genetic manipulation and targeting of agronomic traits. Transcriptome sequencing analysis of callus under salt stress is conducive to study salt tolerance genes. In this study, calli were induced to provide materials for extracting protoplasts and also for screening salt tolerance genes. Calli were successfully induced from leaves of Gossypium sturtianum (C₁ genome) and hypocotyls of G. raimondii (D₅ genome), and embryogenic calli of G. sturtianum and G. raimondii were induced on a differentiation medium with different concentrations of 2, 4-D, KT, and IBA, respectively. In addition, embryogenic calli were also induced successfully from G. raimondii through suspension cultivation. Transcriptome sequencing analysis was performed on the calli of G. raimondii and G. sturtianum, which were treated with 200 mM NaCl at 0, 6, 12, 24, and 48 h, and a total of 12,524 genes were detected with different expression patterns under salt stress. Functional analysis showed that 3,482 genes, which were differentially expressed in calli of G. raimondii and G. sturtianum, were associated with biological processes of nucleic acid binding, plant hormone (such as ABA) biosynthesis, and signal transduction. We demonstrated that DEGs or TFs which related to ABA metabolism were involved in the response to salt stress, including xanthoxin dehydrogenase genes (ABA2), sucrose non-fermenting 1-related protein kinases (SnRK2), NAM, ATAT1/2, and CUC2 transcription factors (NAC), and WRKY class of zinc-finger proteins (WRKY). This research has successfully induced calli from two diploid cotton species and revealed new genes responding to salt stress in callus tissue, which will lay the foundation for protoplast fusion for further understanding of salt stress responses in cotton.

Profiling of MicroRNAs Involved in Mepiquat Chloride-Mediated Inhibition of Internode Elongation in Cotton (Gossypium hirsutum L.) Seedlings

Mepiquat chloride (MC) is the most important plant growth retardant that is widely used in cotton (Gossypium hirsutum L.) production to suppress excessive vegetative growth and improve plant architecture. MicroRNAs (miRNAs) are important gene expression regulators that control plant growth and development. However, miRNA-mediated post-transcriptional regulation in MC-induced growth inhibition remains unclear. In this study, the dynamic expression profiles of miRNAs responsive to MC in cotton internodes were investigated. A total of 508 known miRNAs belonging to 197 families and five novel miRNAs were identified. Among them, 104 miRNAs were differentially expressed at 48, 72, or 96 h post MC treatment compared with the control (0 h); majority of them were highly conserved miRNAs. The number of differentially expressed miRNAs increased with time after treatment. The expression of 14 known miRNAs was continuously suppressed, whereas 12 known miRNAs and one novel miRNA were continuously induced by MC. The expression patterns of the nine differentially expressed miRNAs were verified using qRT-PCR. The targets of the known and novel miRNAs were predicted. Four conserved and six novel targets were validated using the RLM-5' RACE assay. This study revealed that miRNAs play crucial regulatory roles in the MC-induced inhibition of internode elongation. It can improve our understanding of post-transcriptional gene regulation in MC-mediated growth inhibition and could potentially facilitate the breeding of dwarf cotton.

Identification and Investigation of miRNAs From Gastrodia elata Blume and Their Potential Function

Gastrodia elata Blume (G. elata) is a valuable traditional Chinese medicine with neuroprotection, anti-inflammatory, and immune regulatory functions. MicroRNAs (miRNA) is a kind of endogenous noncoding small RNAs that plays distinctly important roles for gene regulation of organisms. So far, the research on G. elata is mainly focused on the pharmacological functions of the natural chemical ingredients, and the function of G. elata miRNA remains unknown. In this study, 5,718 known miRNAs and 38 novel miRNAs were identified by high-throughput sequencing from G. elata. Based on GO and KEGG analysis, we found that the human genes possibly regulated by G. elata miRNAs were related to the cell cycle, immune regulation, intercellular communication, etc. Furthermore, two novel miRNAs as Gas-miR01 and Gas-miR02 have stable and high expression in the medicinal tissues of G. elata. Further bioinformatics prediction showed that both Gas-miR01 and Gas-miR02 could target Homo sapiens A20 gene, furthermore, the dual-luciferase reporter gene assay and Western Blotting verified the interaction of Gas-miR01 or Gas-miR02 with A20. These evidences suggested that G. elata-unique miRNAs might be involved in certain physiological processes. The animal experiment showed that Gas-miR01 and Gas-miR02 could be detected in some tissues of mice by intragastric administration; meanwhile, the A20 expression in some tissues of mice was downregulated. These results supported for the functional study of G. elata miRNAs.

Mitochondrial proteomic analysis reveals that proteins relate to oxidoreductase activity play a central role in pollen fertility in cotton

Cotton (Gossypium hirsutum L.) is an important economic crop. Cytoplasm male sterility (CMS) has been used to develop hybrid system and to produce hybrid seeds in cotton, but the molecular mechanism of CMS remains unclear. Mitochondria are semi-autonomous organelles, which play an important role in the reproduction of flowering plants. Male sterility has been proved associated with mitochondrial dysfunction in plants. In present study, a new strategy of proteomic sequencing data-independent acquisition (DIA) was used to analysis protein abundance across CMS lines 2074A (cytoplasm of Gossypium harknessii, D_2-2) and 2074S (cytoplasm of G. hirsutum, AD₁), and their maintainer 2074B. Comparing with transcriptome results showed that there is little consistence between proteome and transcriptome. A total of 2095 protein species were identified in three materials, and 186 and 161 differentially proteins were detected in the comparisons of 2074A vs 2074B, and 2074S vs 2074B, respectively. Among them, 49 and 50 proteins were specific existed in anther, and mainly participated in oxidoreductase activity, carbohydrate metabolism, fatty acid metabolism, cell aging, wax or cutin deposition and signal transduction. Gh_A07G0770 and Gh_D05G1908 were specific up-regulated in sterility lines, and the other genes Gh_D08G1196, Gh_D12G1971, Gh_A11G1250, Gh_D08G0388 were down-regulated, which presented similar expression tendency verified by qRT-PCR, transcriptome and proteome results. These six genes related to lipid synthesis, response to oxidative stress and cell aging, suggested them being involved in CMS occurrence. Using virus-induced gene silencing (VIGS) system, sterility obtained demonstrated the silencing Gh_A11G1250 in maintainer 2074B led to partial anthers abortion. Gh_A11G1250 encoded a mitochondrial localization of peroxisomal-like protein, participated in response to reactive oxygen species (ROS). Twenty-two proteins interacting with Gh_A11G1250 mainly related to chlorophyll biosynthetic process, photoperiodism and flowering, which showed different expression pattern between the male sterile line 2074A and maintainer 2074B. This novel research based on mitochondrial proteomics comparison confirmed that DAPs related to oxidative stress are critical to pollen abortion. BIOLOGICAL SIGNIFICANCE: Cytoplasm male sterility (CMS) system is utilized widely for hybrid production in cotton. However, the genetic and molecular mechanisms of CMS still need to be further elucidated. Up till now, fewer comprehensive comparisons of the mitochondrial proteomes from cotton CMS line and maintainer line have been reported. In this study, we performed a novel comparison of mitochondrial protein profiles in two CMS lines and their common maintainer line. Based on our results, we found a potential protein related to oxidative stress led to the anthers abortion. These results accumulate data to interpret the molecular mechanisms of CMS in cotton.

Sex Determination Using RNA-Sequencing Analyses in Early Prenatal Pig Development

Sexual dimorphism is a relevant factor in animal science, since it can affect the gene expression of economically important traits. Eventually, the interest in the prenatal phase in a transcriptome study may not comprise the period of development in which male and female conceptuses are phenotypically divergent. Therefore, it would be interesting if sex differentiation could be performed using transcriptome data, with no need for extra techniques. In this study, the sex of pig conceptuses (embryos at 25 days-old and fetuses at 35 days-old) was determined by reads counts per million (CPM) of Y chromosome-linked genes that were discrepant among samples. Thus, ten genes were used: DDX3Y, KDM5D, ZFY, EIF2S3Y, EIF1AY, LOC110255320, LOC110257894, LOC396706, LOC100625207, and LOC110255257. Conceptuses that presented reads CPM sum for these genes (ΣCPM_chrY) greater than 400 were classified as males and those with ΣCPM_chrY below 2 were classified as females. It was demonstrated that the sex identification can be performed at early stages of pig development from RNA-sequencing analysis of genes mapped on Y chromosome. Additionally, these results reinforce that sex determination is a mechanism conserved across mammals, highlighting the importance of using pigs as an animal model to study sex determination during human prenatal development.

Identification of Differentially Expressed miRNAs between a Wheat K-type Cytoplasmic Male Sterility Line and Its Near-Isogenic Restorer Line

K-type cytoplasmic male sterility (KCMS) lines were ideal material for three-line hybrid wheat system due to the major role in hybrid wheat production. In this study, the morphology of developing microspore and mature pollen was compared between a KCMS line and its near-isogenic restorer line (KCMS-NIL). The most striking difference is that the microspore was unable to develop into tricellular pollen in the KCMS line. MicroRNA plays vital roles in flowering and gametophyte development. Small RNA sequencing identified a total of 274 known and 401 novel miRNAs differentially expressed between two lines or two developmental stages. Most of miRNAs with high abundance were differentially expressed at the uninucleate stage, and their expression level recovered or remained at the binucleate stage. Further degradome sequencing identified target genes which were mainly enriched in transcription regulation, phytohormone signaling and RNA degradation pathways. Combining with the transcriptome data, a correlation was found between the abnormal anther development, such as postmeiotic mitosis cessation, deformative pollen wall and the chromosome condensation of the vegetative cell, and the alterations in the related miRNA and their targets expression profiles. According to the correlation and pathway analysis, we propose a hypothetic miRNA-mediated network for the control of KCMS restoration.

Transcriptome analysis identified aberrant gene expression in pollen developmental pathways leading to CGMS in cotton (Gossypium hirsutum L.)

Male sterility (induced or natural) is a potential tool for commercial hybrid seed production in different crops. Despite numerous endeavors to understand the physiological, hereditary, and molecular cascade of events governing CMS in cotton, the exact biological process controlling sterility and fertility reconstruction remains obscure. During current study, RNA-Seq using Ion Torrent S5 platform is carried out to identify 'molecular portraits' in floral buds among the Cytoplasmic Genic Male Sterility (CGMS) line, its near-isogenic maintainer, and restorer lines. A total of 300, 438 and 455 genes were differentially expressed in CGMS, Maintainer, and Restorer lines respectively. The functional analysis using AgriGo revealed suppression in the pathways involved in biogenesis and metabolism of secondary metabolites which play an important role in pollen and anther maturation. Enrichment analysis showed dearth related to pollen and anther's development in sterile line, including anomalous expression of genes and transcription factors that have a role in the development of the reproductive organ, abnormal cytoskeleton formation, defects in cell wall formation. The current study found aberrant expression of DYT1, AMS and cytochrome P450 genes involved in tapetum formation, pollen development, pollen exine and anther cuticle formation associated to male sterility as well as fertility restoration of CGMS. In the current study, more numbers of DEGs were found on Chromosome D05 and A05 as compared to other chromosomes. Expression pattern analysis of fourteen randomly selected genes using qRT-PCR showed high concurrence with gene expression profile of RNA-Seq analysis accompanied by a strong correlation of 0.82. The present study provides an important support for future studies in identifying interaction between cyto-nuclear molecular portraits, to accelerate functional genomics and molecular breeding related to cytoplasmic male sterility studies in cotton.