RNA-Seq with a novel glabrous-ZM24fl reveals some key lncRNAs and the associated targets in fiber initiation of cotton

BACKGROUND

Cotton fiber is an important natural resource for textile industry and an excellent model for cell biology study. Application of glabrous mutant cotton and high-throughput sequencing facilitates the identification of key genes and pathways for fiber development and cell differentiation and elongation. LncRNA is a type of ncRNA with more than 200 nt in length and functions in the ways of chromatin modification, transcriptional and post-transcriptional modification, and so on. However, the detailed lncRNA and associated mechanisms for fiber initiation are still unclear in cotton.

RESULTS

In this study, we used a novel glabrous mutant ZM24fl, which is endowed with higher somatic embryogenesis, and functions as an ideal receptor for cotton genetic transformation. Combined with the high-throughput sequencing, fatty acid pathway and some transcription factors such as MYB, ERF and bHLH families were identified the important roles in fiber initiation; furthermore, 3,288 lncRNAs were identified, and some differentially expressed lncRNAs were also analyzed. From the comparisons of ZM24_0 DPA vs ZM24_-2 DPA and fl_0 DPA vs ZM24_0 DPA, one common lncRNA MSTRG 2723.1 was found that function upstream of fatty acid metabolism, MBY25-mediating pathway, and pectin metabolism to regulate fiber initiation. In addition, other lncRNAs MSTRG 3390.1, MSTRG 48719.1, and MSTRG 31176.1 were also showed potential important roles in fiber development; and the co-expression analysis between lncRNAs and targets showed the distinct models of different lncRNAs and complicated interaction between lncRNAs in fiber development of cotton.

CONCLUSIONS

From the above results, a key lncRNA MSTRG 2723.1 was identified that might mediate some key genes transcription of fatty acid metabolism, MYB25-mediating pathway, and pectin metabolism to regulate fiber initiation of ZM24 cultivar. Co-expression analysis implied that some other important lncRNAs (e.g., MSTRG 3390.1, MSTRG 48719.1, and MSTRG 31176.1) were also showed the different regulatory model and interaction between them, which proposes some valuable clues for the lncRNAs associated mechanisms in fiber development.

Co-expression network and comparative transcriptome analysis for fiber initiation and elongation reveal genetic differences in two lines from upland cotton CCRI70 RIL population

Upland cotton is the most widely planted for natural fiber around the world, and either lint percentage (LP) or fiber length (FL) is the crucial component tremendously affecting cotton yield and fiber quality, respectively. In this study, two lines MBZ70-053 and MBZ70-236 derived from G. hirsutum CCRI70 recombinant inbred line (RIL) population presenting different phenotypes in LP and FL traits were chosen to conduct RNA sequencing on ovule and fiber samples, aiming at exploring the differences of molecular and genetic mechanisms during cotton fiber initiation and elongation stages. As a result, 249/128, 369/206, 4296/1198 and 3547/2129 up-/down- regulated differentially expressed genes (DGEs) in L2 were obtained at -3, 0, 5 and 10 days post-anthesis (DPA), respectively. Seven gene expression profiles were discriminated using Short Time-series Expression Miner (STEM) analysis; seven modules and hub genes were identified using weighted gene co-expression network analysis. The DEGs were mainly enriched into energetic metabolism and accumulating as well as auxin signaling pathway in initiation and elongation stages, respectively. Meanwhile, 29 hub genes were identified as 14-3-3ω , TBL35, GhACS, PME3, GAMMA-TIP, PUM-7, etc., where the DEGs and hub genes revealed the genetic and molecular mechanisms and differences during cotton fiber development.

Characterization of cotton ARF factors and the role of GhARF2b in fiber development

BACKGROUND

Cotton fiber is a model system for studying plant cell development. At present, the functions of many transcription factors in cotton fiber development have been elucidated, however, the roles of auxin response factor (ARF) genes in cotton fiber development need be further explored.

RESULTS

Here, we identify auxin response factor (ARF) genes in three cotton species: the tetraploid upland cotton G. hirsutum, which has 73 ARF genes, and its putative extent parental diploids G. arboreum and G. raimondii, which have 36 and 35 ARFs, respectively. Ka and Ks analyses revealed that in G. hirsutum ARF genes have undergone asymmetric evolution in the two subgenomes. The cotton ARFs can be classified into four phylogenetic clades and are actively expressed in young tissues. We demonstrate that GhARF2b, a homolog of the Arabidopsis AtARF2, was preferentially expressed in developing ovules and fibers. Overexpression of GhARF2b by a fiber specific promoter inhibited fiber cell elongation but promoted initiation and, conversely, its downregulation by RNAi resulted in fewer but longer fiber. We show that GhARF2b directly interacts with GhHOX3 and represses the transcriptional activity of GhHOX3 on target genes.

CONCLUSION

Our results uncover an important role of the ARF factor in modulating cotton fiber development at the early stage.

GaHD1, a candidate gene for the Gossypium arboreum SMA-4 mutant, promotes trichome and fiber initiation by cellular H2O2 and Ca2+ signals

Cotton fibers are initiated from the epidermal cells of the ovule before or on the day of anthesis. Gossypium arboreum SMA-4 mutant contains recessive mutation (sma-4(ha)) and has the phenotypes of fibreless seeds and glabrous stems. In this study, fine mapping and alternative splicing analysis indicated a nucleotide substitution (AG → AC) at splicing site in a homeodomain-leucine zipper IV family gene (GaHD1) might cause gene A3S (Alternative 3' splicing) mistake, suggested that GaHD1 was the candidate gene of sma-4(ha). Many genes related to the fiber initiation are identified to be differentially expressed in the mutant which could result in the blocked fiber initiation signals such as H₂O_2, or Ca in the mutant. Further comparative physiological analysis of H₂O₂ production and Ca²⁺ flux in the SMA-4 and wide type cotton confirmed that H₂O₂ and Ca were important fiber initiation signals and regulated by GaHD1. The in vitro ovule culture of the mutant with hormones recovered the fibered phenotype coupled with the restoration of these signals. Overexpressing of GaHD1 in Arabidopsis increased trichome densities on the sepal, leaf, and stem tissues while transient silencing of the GaHD1 gene in G. arboreum reduced the trichome densities. These phenotypes indicated that GaHD1 is the candidate gene of SMA-4 with a crucial role in acting upstream molecular switch of signal transductions for cotton trichome and fiber initiations.

Transcriptome profiling of Gossypium arboreum during fiber initiation and the genome-wide identification of trihelix transcription factors

Cotton fiber initiation is the first step in fiber development, and it determines the yield. Here, genome-wide transcriptome profiling of Gossypium arboreum was performed to determine the molecular basis of cotton fiber initiation. A comparison of the transcriptomes of fiber-bearing ovules at -0.5, 0, 0.5, 1, 1.5, 2, 2.5 and 3 d post-anthesis detected 12,049 differentially expressed genes that mainly participated in ribosome, carbon metabolism and amino acid biosynthesis pathways. Genes encoding alcohol dehydrogenase 1 and hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferase, involving in fatty acid degradation and flavonoid biosynthesis, were enriched. Furthermore, 1049 differentially expressed transcription factors were identified. Among these, 17 were trihelix family transcription factors, which play important roles in plant development and responses to biotic and abiotic stresses. In total, 52 full-length trihelix genes, named as GaGTs, were identified in G. arboreum and located in 12 of the 13 cotton chromosomes. Transcriptomic data and a quantitative real-time PCR analysis indicated that several GaGTs were significantly induced during fiber initiation in G. arboreum. Thus, the genome-wide comprehensive analysis of gene expression in G. arboreum fiber initiation will serve as a useful resource for unraveling the functions of specific genes. The phylogenetic relationships and expression analyses of the G. arboreum trihelix genes established a solid foundation for future comprehensive functional analyses of the GaGTs.

Identification and profiling of upland cotton microRNAs at fiber initiation stage under exogenous IAA application

Background

Cotton is the most essential textile crop worldwide, and phytohormones are critical for cotton fiber development. One example is the role of auxin in fiber initiation, but we know little molecular basis. MicroRNAs (miRNAs) have a significant function in cotton development; nevertheless their role in fiber initiation remains unclear. Here, exogenous IAA was applied to cotton plant before anthesis. Utilizing small RNA sequencing, the mechanism underlying miRNA-mediated regulation of fiber initiation under exogenous IAA treatment was investigated.

Results

With exogenous IAA application, the endogenous IAA and GA contents of IAA treated (IT) ovules were higher than control (CK) ovules at the fiber initiation stage, while endogenous ABA content was lower in IT than CK. Using scanning electron microscopy, we found the fiber number and size were significantly promoted in IT at 0 DPA. Fiber quality analysis showed that fiber length, uniformity, strength, elongation, and micronaire of IT were higher than CK, though not statistically significant, while lint percent was significantly higher in IT. We generated six small RNA libraries using − 3, 0, and 3 DPA ovules of IT and CK, and identified 58 known miRNAs and 83 novel miRNAs together with the target genes. The differential expressed miRNAs number between IT and CK at − 3, 0, 3 DPA was 34, 16 and 24, respectively. Gene ontology and KEGG pathway enrichment analyses for the target genes of the miRNAs expressed in a differential manner showed that they were significantly enriched in 30 terms and 8 pathways. QRT-PCR for those identified miRNAs and the target genes related to phytohormones and fiber development was performed, and results suggested a potential role of these miRNAs in fiber initiation.

Conclusions

The exogenous IAA application affected the relative phytohormone contents in ovule and promoted fiber initiation in cotton. Identification and profiling of miRNAs and their targets at the fiber initiation stage provided insights for miRNAs’ regulation function of fiber initiation. These findings not only shed light on the regulatory network of fiber growth but also offer clues for cotton fiber amelioration strategies in cotton.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5760-8) contains supplementary material, which is available to authorized users.

In vivo imaging of Ca2+ accumulation during cotton fiber initiation using fluorescent indicator YC3.60

Non-tip-focused Ca ²⁺ gradient indicated by genetically expressing a FRET-based calcium sensor YC3.60 was established in spherical expanding cotton fibers, which is vital for cotton fiber initiation. Cotton fiber is a single cell elongated from ovule epidermis. It is not only the most important natural fiber used in the textile industry but also an ideal model for studying cell differentiation and elongation. Before linear cell growth, cotton fibers undergo spherical expansion at the beginning of initiation. Ca²⁺, as an important secondary messenger, plays a central role in polarized cell growth including cotton fiber elongation. However, the role of Ca²⁺ in fiber initiation is far from well understood. In this paper, through ovule culture we demonstrate that Ca²⁺ is crucial for fiber initiation. Using transgenic cotton expressing the fluorescent Ca²⁺ indicator YC3.60, we show cellular and intracellular distribution of Ca²⁺ in cotton ovule epidermis and fiber cells. In the initiating fiber cell, Ca²⁺ accumulated mainly at the base of the cell, while in the fast elongating cell, the Ca²⁺ was enriched in the tip region. This cellular distribution of Ca²⁺ reported by YC3.60 was confirmed by the staining with a Ca²⁺-sensitive dye fluo-3/AM. Compared to the fluorescent dye staining, the YC3.60 system can reveal more detailed information on the intracellular distribution without photobleaching. Taken together, our data suggest that Ca²⁺ plays an important role in spherical expansion of cotton fiber initials.