Long non-coding RNAs and their potential functions in Ligon-lintless-1 mutant cotton during fiber development

Exploring the regulatory role of non-coding RNAs in fiber development and direct regulation of GhKCR2 in the fatty acid metabolic pathway in upland cotton

Cotton fiber holds immense importance as the primary raw material for the textile industry. Consequently, comprehending the regulatory mechanisms governing fiber development is pivotal for enhancing fiber quality. Our study aimed to construct a regulatory network of competing endogenous RNAs (ceRNAs) and assess the impact of non-coding RNAs on gene expression throughout fiber development. Through whole transcriptome data analysis, we identified differentially expressed genes (DEGs) regulated by non-coding RNA (ncRNA) that were predominantly enriched in phenylpropanoid biosynthesis and the fatty acid elongation pathway. This analysis involved two contrasting phenotypic materials (J02-508 and ZRI015) at five stages of fiber development. Additionally, we conducted a detailed analysis of genes involved in fatty acid elongation, including KCS, KCR, HACD, ECR, and ACOT, to unveil the factors contributing to the variation in fatty acid elongation between J02-508 and ZRI015. Through the integration of histochemical GUS staining, dual luciferase assay experiments, and correlation analysis of expression levels during fiber development stages for lncRNA MSTRG.44818.23 (MST23) and GhKCR2, we elucidated that MST23 positively regulates GhKCR2 expression in the fatty acid elongation pathway. This identification provides valuable insights into the molecular mechanisms underlying fiber development, emphasizing the intricate interplay between non-coding RNAs and protein-coding genes.

An analysis of lncRNAs related to fiber quality and the discovery of their target genes in a Gossypium hirsutum line with Gossypium mustelinum introgression

KEY MESSAGE

Analysis of fiber quality lncRNAs and their target genes from a pair of Gossypium mustelinum near-isogenic lines provide new prospects for improving the fiber quality of Upland cotton. Long noncoding RNAs (lncRNAs) are an important part of genome transcription and play roles in a wide range of biological processes in plants. In this research, a pair of near-isogenic cotton lines, namely, a Gossypium mustelinum introgression line (IL9) with outstanding fiber quality and its recurrent Upland cotton parent (PD94042), were used as the experimental materials. Cotton fibers were selected for lncRNA sequencing at 17 and 21 days post-anthesis. A total of 2693 differentially expressed genes were identified. In total, 5841 lncRNAs were ultimately screened, from which 163 differentially expressed lncRNAs were identified. Target genes of the lncRNAs were predicted by two different methods: cis and trans. Some of the target genes were related to cell components, membrane components, plant hormone signal transduction and catalytic metabolism, and the results indicated that there might also be important effects on the development of fiber. Four differentially expressed target genes related to fiber quality (Gomus.D05G015100, Gomus.A05G281300, Gomus.A12G023400 and Gomus.A10G226800) were screened through gene function annotation, and the functions of these four genes were verified through virus-induced gene silencing (VIGS). Compared to the negative controls, plants in which any of these four genes were silenced showed significant reductions in fiber strength. In addition, the plants in which the Gomus.A12G023400 gene was silenced showed a significant reduction in fiber uniformity, whereas the plants in which Gomus.A05G281300 was silenced showed a significant increase in fiber fineness as measured via micronaire. Our results showed that these genes play different roles during fiber development, impacting fiber quality.

Long non-coding RNAs and their potential function in response to postharvest senescence of Sparassis latifolia during cold storage

Long non-coding RNAs (lncRNAs) have been shown to play crucial roles in response to aging processes. However, how lncRNAs regulate postharvest senescence of Sparassis latifolia (S. latifolia) with oriented polypropylene (OPP) film packing during cold storage remains unclear. In this study, we performed RNA-seq using the fruiting bodies of S. latifolia stored at 4 ℃ for 0, 8, 16 and 24 days after harvest, and profiled the lncRNA and mRNA transcriptome, respectively. In total, 1003 putative lncRNAs were identified, and there were 495, 483 and 162 differentially expressed (DE) lncRNAs, and 3680, 3941 and 1870 differentially expressed mRNAs after 8, 16 and 24 days of storage, respectively, compared to 0 day of storage. Target genes of differentially expressed lncRNAs were found to significantly associate with carbon and energy metabolism, response to abiotic stimulus, amino acid biosynthesis and metabolism, and protein synthesis and transcription. In addition, DE-lncRNA-mRNA co-expression networks in response to aging stress were also constructed. Taken together, these results confirm the regulatory role of lncRNAs in postharvest senescence of S. latifolia and will facilitate for improving preservation method.

VirhostlncR: A comprehensive database to explore lncRNAs and their targets in viral infections

Long non-coding-RNAs (lncRNAs) are an expanding set of cis-/trans-regulatory RNA genes that outnumber the protein-coding genes. Although being increasingly discovered, the functional role of the majority of lncRNAs in diverse biological conditions is undefined. Increasing evidence supports the critical role of lncRNAs in the emergence, regulation, and progression of various viral infections including influenza, hepatitis, coronavirus, and human immunodeficiency virus. Hence, the identification of signature lncRNAs would facilitate focused analysis of their functional roles accounting for their targets and regulatory mechanisms associated with infections. Towards this, we compiled 2803 lncRNAs identified to be modulated by 33 viral strains in various mammalian cell types and are provided through the resource named VirhostlncR (http://ciods.in/VirhostlncR/). The information on each of the viral strains, their multiplicity of infection, duration of infection, host cell name and cell types, fold change of lncRNA expression, and their specific identification methods are integrated into VirhostlncR. Based on the current datasets, we report 150 lncRNAs including differentiation antagonizing non-protein coding RNA (DANCR), metastasis-associated lung adenocarcinoma transcript 1 (MALAT1), maternally expressed gene 3 (MEG3), nuclear paraspeckle assembly transcript 1 (NEAT1), and plasmacytoma variant translocation 1 (PVT1) to be perturbed by two or more viruses. Analysis of viral protein interactions with human transcription factors (TFs) or TF-containing protein complexes identified that distinct viruses can transcriptionally regulate many of these lncRNAs through multiple protein complexes. Together, we believe that the current dataset will enable priority selection of lncRNAs for identification of their targets and serve as an effective platform for the analysis of noncoding RNA-mediated regulations in viral infections.

Long non-coding RNA and microRNA landscape of two major domesticated cotton species

Allotetraploid cotton plants Gossypium hirsutum and Gossypium barbadense have been widely cultivated for their natural, renewable textile fibres. Even though ncRNAs in domesticated cotton species have been extensively studied, systematic identification and annotation of lncRNAs and miRNAs expressed in various tissues and developmental stages under various biological contexts are limited. This influences the comprehension of their functions and future research on these cotton species. Here, we report high confidence lncRNAs and miRNA collection from G. hirsutum accession and G. barbadense accession using large-scale RNA-seq and small RNA-seq datasets incorporated into a user-friendly database, CoNCRAtlas. This database provides a wide range and depth of lncRNA and miRNA annotation based on the systematic integration of extensive annotations such as expression patterns derived from transcriptome data analysis in thousands of samples, as well as multi-omics annotations. We assume this comprehensive resource will accelerate evolutionary and functional studies in ncRNAs and inform future breeding programs for cotton improvement. CoNCRAtlas is accessible at http://www.nipgr.ac.in/CoNCRAtlas/.

Transcriptome, Ectopic Expression and Genetic Population Analysis Identify Candidate Genes for Fiber Quality Improvement in Cotton

Comparative transcriptome analysis of fiber tissues between Gossypium barbadense and Gossypium hirsutum could reveal the molecular mechanisms underlying high-quality fiber formation and identify candidate genes for fiber quality improvement. In this study, 759 genes were found to be strongly upregulated at the elongation stage in G. barbadense, which showed four distinct expression patterns (I-IV). Among them, the 346 genes of group IV stood out in terms of the potential to promote fiber elongation, in which we finally identified 42 elongation-related candidate genes by comparative transcriptome analysis between G. barbadense and G. hirsutum. Subsequently, we overexpressed GbAAR3 and GbTWS1, two of the 42 candidate genes, in Arabidopsis plants and validated their roles in promoting cell elongation. At the secondary cell wall (SCW) biosynthesis stage, 2275 genes were upregulated and exhibited five different expression profiles (I-V) in G. barbadense. We highlighted the critical roles of the 647 genes of group IV in SCW biosynthesis and further picked out 48 SCW biosynthesis-related candidate genes by comparative transcriptome analysis. SNP molecular markers were then successfully developed to distinguish the SCW biosynthesis-related candidate genes from their G. hirsutum orthologs, and the genotyping and phenotyping of a BC₃F₅ population proved their potential in improving fiber strength and micronaire. Our results contribute to the better understanding of the fiber quality differences between G. barbadense and G. hirsutum and provide novel alternative genes for fiber quality improvement.

Investigating nicotine pathway-related long non-coding RNAs in tobacco

Long non-coding RNAs (lncRNAs) are transcripts longer than 200 bp with low or no protein-coding ability, which play essential roles in various biological processes in plants. Tobacco is an ideal model plant for studying nicotine biosynthesis and metabolism, and there is little research on lncRNAs in this field. Therefore, how to take advantage of the mature tobacco system to profoundly investigate the lncRNAs involved in the nicotine pathway is intriguing. By exploiting 549 public RNA-Seq datasets of tobacco, 30,212 lncRNA candidates were identified, including 24,084 large intervening non-coding RNAs (lincRNAs), 5,778 natural antisense transcripts (NATs) and 350 intronic non-coding RNAs (incRNAs). Compared with protein-coding genes, lncRNAs have distinct properties in terms of exon number, sequence length, A/U content, and tissue-specific expression pattern. lincRNAs showed an asymmetric evolutionary pattern, with a higher proportion (68.71%) expressed from the Nicotiana sylvestris (S) subgenome. We predicted the potential cis/trans-regulatory effects on protein-coding genes. One hundred four lncRNAs were detected as precursors of 30 known microRNA (miRNA) family members, and 110 lncRNAs were expected to be the potential endogenous target mimics for 39 miRNAs. By combining the results of weighted gene co-expression network analysis with the differentially expressed gene analysis of topping RNA-seq data, we constructed a sub-network containing eight lncRNAs and 25 nicotine-related coding genes. We confirmed that the expression of seven lncRNAs could be affected by MeJA treatment and may be controlled by the transcription factor NtMYC2 using a quantitative PCR assay and gene editing. The results suggested that lncRNAs are involved in the nicotine pathway. Our findings further deepened the understanding of the features and functions of lncRNAs and provided new candidates for regulating nicotine biosynthesis in tobacco.

The Conservation of Long Intergenic Non-Coding RNAs and Their Response to Verticillium dahliae Infection in Cotton

Long intergenic non-coding RNAs (lincRNAs) have been demonstrated to be vital regulators of diverse biological processes in both animals and plants. While many lincRNAs have been identified in cotton, we still know little about the repositories and conservativeness of lincRNAs in different cotton species or about their role in responding to biotic stresses. Here, by using publicly available RNA-seq datasets from diverse sources, including experiments of Verticillium dahliae (Vd) infection, we identified 24,425 and 17,713 lincRNAs, respectively, in Gossypium hirsutum (Ghr) and G. barbadense (Gba), the two cultivated allotetraploid cotton species, and 6933 and 5911 lincRNAs, respectively, in G. arboreum (Gar) and G. raimondii (Gra), the two extant diploid progenitors of the allotetraploid cotton. While closely related subgenomes, such as Ghr_At and Gba_At, tend to have more conserved lincRNAs, most lincRNAs are species-specific. The majority of the synthetic and transcribed lincRNAs (78.2%) have a one-to-one orthologous relationship between different (sub)genomes, although a few of them (0.7%) are retained in all (sub)genomes of the four species. The Vd responsiveness of lincRNAs seems to be positively associated with their conservation level. The major functionalities of the Vd-responsive lincRNAs seem to be largely conserved amongst Gra, Ghr, and Gba. Many Vd-responsive Ghr-lincRNAs overlap with Vd-responsive QTL, and several lincRNAs were predicted to be endogenous target mimicries of miR482/2118, with a pair being highly conserved between Ghr and Gba. On top of the confirmation of the feature characteristics of the lincRNAs previously reported in cotton and other species, our study provided new insights into the conservativeness and divergence of lincRNAs during cotton evolution and into the relationship between the conservativeness and Vd responsiveness of lincRNAs. The study also identified candidate lincRNAs with a potential role in disease response for functional characterization.

Long intergenic non-protein coding RNA 847 promotes laryngeal squamous cell carcinoma progression through the microRNA-181a-5p/zinc finger E-box binding homeobox 2 axis

The present study is targeted at investigating the effects of long intergenic non-protein coding RNA 847 (LINC00847) on the malignant biological behaviors of laryngeal squamous cell carcinoma (LSCC) cells, and the mechanisms. Quantitative real-time PCR and Western blotting were conducted for detecting the expressions of LINC00847, microRNA-181a-5p (miR-181a-5p) and zinc finger E-box binding homeobox 2 (ZEB2) in LSCC cell lines and tissue samples. BrdU, cell counting kit-8, scratch wound healing, Transwell and flow cytometry assays were utilized for detecting cell proliferation, migration, invasion, and cell cycle progression. Dual-luciferase reporter gene, RNA binding protein immunoprecipitation (RIP), and RNA pull-down assays were utilized to investigate the interaction among LINC00847, miR-181a-5p, and ZEB2. The subcellular location of LINC00847 was determined by RNA fluorescence in situ hybridization (RNA-FISH) assay. Tumor growth was evaluated using a xenograft model of nude mice. It was revealed that LINC00847 expression was increased in LSCC tissues, and its high expression was associated with lymph node metastasis and poor differentiation. LINC00847 was mainly located in the cytoplasm of LSCC cells, and LINC00847 overexpression promoted LSCC cell proliferation, migration, invasion, and accelerated the cell cycle progression while knocking down LINC00847 had the opposite effects in vitro and inhibited the tumor growth in vivo. LINC00847 directly targeted miR-181a-5p and negatively modulated miR-181a-5p expression. ZEB2 was a target gene of miR-181a-5p, and was positively and indirectly modulated by LINC00847. Our data suggest that LINC00847 promotes LSCC progression by regulating the miR-181a-5p/ZEB2 axis.

Identification of long non-coding RNAs in Verticillium dahliae following inoculation of cotton

Long non-coding RNAs (lncRNAs) play important roles in diverse biological processes. However, these functions have not been assessed in Verticillium dahliae, a soil-borne fungal pathogen that causes devastating wilt diseases in many crops. The discovery and identity of novel lncRNAs and their association with virulence may contribute to an increased understanding of the regulation of virulence in V. dahliae. Here, we identified a total of 352 lncRNAs in V. dahliae. The lncRNAs were transcribed from all V. dahliae chromosomes, typically with shorter open reading frames, lower GC content, and fewer exons than protein-coding genes. In addition, 308 protein-coding genes located within 10 kb upstream and 10 kb downstream of lncRNAs were identified as neighboring genes, and which were considered as potential targets of lncRNA. These neighboring genes encode products involved in development, stress responses, and pathogenicity of V. dahliae, such as transcription factors (TF), kinase, and members of the secretome. Furthermore, 47 lncRNAs were significantly differentially expressed in V. dahliae following inoculation of susceptible cotton (Gossyoiumhisutum) cultivar Junmian No.1, suggesting that lncRNAs may be involved in the regulation of virulence in V. dahliae. Moreover, correlations in expression patterns between lncRNA and their neighboring genes were detected. Expression of lncRNA012077 and its neighboring gene was up-regulated 6 h following inoculation of cotton, while the expression of lncRNA007722 was down-regulated at 6 h but up-regulated at 24 h, in a pattern opposite to that of its neighboring gene. Overexpression of lncRNA012077 in wild-type strain (Vd991) enhanced its virulence on cotton while overexpression of lncRNA009491 reduced virulence. Identification of novel lncRNAs and their association with virulence may provide new targets for disease control.

Functional examination of lncRNAs in allotetraploid Gossypium hirsutum

BACKGROUND

An evolutionary model using diploid and allotetraploid cotton species identified 80 % of non-coding transcripts in allotetraploid cotton as being uniquely activated in comparison with its diploid ancestors. The function of the lncRNAs activated in allotetraploid cotton remain largely unknown.

RESULTS

We employed transcriptome analysis to examine the relationship between the lncRNAs and mRNAs of protein coding genes (PCGs) in cotton leaf tissue under abiotic stresses. LncRNA expression was preferentially associated with that of the flanking PCGs. Selected highly-expressed lncRNA candidates (n = 111) were subjected to a functional screening pilot test in which virus-induced gene silencing was integrated with abiotic stress treatment. From this low-throughput screen, we obtained candidate lncRNAs relating to plant height and tolerance to drought and other abiotic stresses.

CONCLUSIONS

Low-throughput screen is an effective method to find functional lncRNA for further study. LncRNAs were more active in abiotic stresses than PCG expression, especially temperature stress. LncRNA XLOC107738 may take a cis-regulatory role in response to environmental stimuli. The degree to which lncRNAs are constitutively expressed may impact expression patterns and functions on the individual gene level rather than in genome-wide aggregate.

Noncoding RNAs in Medicinal Plants and their Regulatory Roles in Bioactive Compound Production

BACKGROUND

Noncoding RNAs (ncRNAs), such as microRNAs (miRNAs), small interfering RNAs (siRNAs) and long noncoding RNAs (lncRNAs), play significant regulatory roles in plant development and secondary metabolism and are involved in plant response to biotic and abiotic stresses. They have been intensively studied in model systems and crops for approximately two decades and massive amount of information have been obtained. However, for medicinal plants, ncRNAs, particularly their regulatory roles in bioactive compound biosynthesis, are just emerging as a hot research field.

OBJECTIVE

This review aims to summarize current knowledge on herbal ncRNAs and their regulatory roles in bioactive compound production.

RESULTS

So far, scientists have identified thousands of miRNA candidates from over 50 medicinal plant species and 11794 lncRNAs from Salvia miltiorrhiza, Panax ginseng, and Digitalis purpurea. Among them, more than 30 miRNAs and five lncRNAs have been predicted to regulate bioactive compound production.

CONCLUSION

The regulation may achieve through various regulatory modules and pathways, such as the miR397-LAC module, the miR12112-PPO module, the miR156-SPL module, the miR828-MYB module, the miR858-MYB module, and other siRNA and lncRNA regulatory pathways. Further functional analysis of herbal ncRNAs will provide useful information for quality and quantity improvement of medicinal plants.

Understanding the role of phytohormones in cotton fiber development through omic approaches; recent advances and future directions

Cotton is among the most important fiber crops for the textile-based industry, thanks to its cellulose-rich mature fibers. The fiber initiation and elongation are one of the best models for deciphering mechanisms of single-cell differentiation and growth, that also target of fiber development programs. During the last couple of decades, high yielding omics approaches (genomics, transcriptomics, and proteomics), have helped in the identification of several genes and gene products involved in fiber development along with functional relationship to phytohormones. For example, MYB transcription factor family and Sus gene family have been evidenced by controlling cotton fiber initiation. Most importantly, the biosynthesis, responses, and transporting of phytohormones is documented to participate in the initiation of cotton fibers. Herein, in this review, the reliable genetic evidence by manipulating the above genes in cotton have been summarized to describe the relationships among key phytohormones, transcription factors, proteins, and downstream fiber growth-related genes such as Sus. The effect of other important factors such as ROS, fatty acid metabolism, and actin (globular multi-functional proteins) over fiber development has also been discussed. The challenges and deficiencies in the research of cotton fiber development have been mentioned along with a future perspective to discover new crucial genes using multiple omics analysis.

The Ligon lintless -2 Short Fiber Mutation Is Located within a Terminal Deletion of Chromosome 18 in Cotton

Extreme elongation distinguishes about one-fourth of cotton (Gossypium sp.) seed epidermal cells as "lint" fibers, useful for the textile industry, from "fuzz" fibers (<5 mm). Ligon lintless-2 (Li ₂ ), a dominant mutation that results in no lint fiber but normal fuzz fiber, offers insight into pathways and mechanisms that differentiate spinnable cotton from its progenitors. A genetic map developed using 1,545 F2 plants showed that marker CISP15 was 0.4 cM from Li ₂ , and "dominant" simple sequence repeat (SSR) markers (i.e. with null alleles in the Li ₂ genotype) SSR7 and SSR18 showed complete linkage with Li ₂ Nonrandom distribution of markers with null alleles suggests that the Li ₂ phenotype results from a 176- to 221-kb deletion of the terminal region of chromosome 18 that may have been masked in prior pooled-sample mapping strategies. The deletion includes 10 genes with putative roles in fiber development. Two Glycosyltransferase Family 1 genes showed striking expression differences during elongation of wild-type versus Li ₂ fiber, and virus-induced silencing of these genes in the wild type induced Li ₂ -like phenotypes. Further, at least 7 of the 10 putative fiber development genes in the deletion region showed higher expression in the wild type than in Li ₂ mutants during fiber development stages, suggesting coordinated regulation of processes in cell wall development and cell elongation, consistent with the hypothesis that some fiber-related quantitative trait loci comprise closely spaced groups of functionally diverse but coordinately regulated genes.

Investigation of the EIL/EIN3 Transcription Factor Gene Family Members and Their Expression Levels in the Early Stage of Cotton Fiber Development

The ethylene-insensitive3-like/ethylene-insensitive3 (EIL/EIN3) protein family can serve as a crucial factor for plant growth and development under diverse environmental conditions. EIL/EIN3 protein is a form of a localized nuclear protein with DNA-binding activity that potentially contributes to the intricate network of primary and secondary metabolic pathways of plants. In light of recent research advances, next-generation sequencing (NGS) and novel bioinformatics tools have provided significant breakthroughs in the study of the EIL/EIN3 protein family in cotton. In turn, this paved the way to identifying and characterizing the EIL/EIN3 protein family. Hence, the high-throughput, rapid, and cost-effective meta sequence analyses have led to a remarkable understanding of protein families in addition to the discovery of novel genes, enzymes, metabolites, and other biomolecules of the higher plants. Therefore, this work highlights the recent advance in the genomic-sequencing analysis of higher plants, which has provided a plethora of function profiles of the EIL/EIN3 protein family. The regulatory role and crosstalk of different metabolic pathways, which are apparently affected by these transcription factor proteins in one way or another, are also discussed. The ethylene hormone plays an important role in the regulation of reactive oxygen species in plants under various environmental stress circumstances. EIL/EIN3 proteins are the key ethylene-signaling regulators and play important roles in promoting cotton fiber developmental stages. However, the function of EIL/EIN3 during initiation and early elongation stages of cotton fiber development has not yet been fully understood. The results provided valuable information on cotton EIL/EIN3 proteins, as well as a new vision into the evolutionary relationships of this gene family in cotton species.