Genome-wide analysis of small RNA and novel microRNA discovery during fiber and seed initial development in Gossypium hirsutum. L

Coordination of floral and fiber development in cotton (Gossypium) by hormone- and flavonoid-signalling associated regulatory miRNAs

Various plant development activities and stress responses are tightly regulated by various microRNAs (miRNA) and their target genes, or transcription factors in a spatiotemporal manner. Here, to exemplify how flowering-associated regulatory miRNAs synchronize their expression dynamics during floral and fiber development in cotton, constitutive expression diminution transgenic lines of auxin-signaling regulatory Gh-miR167 (35S-MIM167) were developed through target mimicry approach. 'Moderate' (58% to 80%)- and 'high' (> 80%)-Gh-miR167 diminution mimic lines showed dosage-dependent developmental deformities in anther development, pollen maturation, and fruit (= boll) formation. Cross pollination of 'moderate' 35S-MIM167 mimic lines with wild type (WT) plant partially restored boll formation and emergence of fiber initials on the ovule surface. Gh-miR167 diminution favored organ-specific transcription biases in miR159, miR166 as well as miR160, miR164, and miR172 along with their target genes during anther and petal development, respectively. Similarly, accumulative effect of percent Gh-miR167 diminution, cross regulation of its target ARF6/8 genes, and temporal mis-expression of hormone signaling- and flavonoid biosynthesis-associated regulatory miRNAs at early fiber initiation stage caused irregular fiber formation. Spatial and temporal transcription proportions of regulatory miRNAs were also found crucial for the execution of hormone- and flavonoid-dependent progression of floral and fiber development. These observations discover how assorted regulatory genetic circuits get organized in response to Gh-miR167 diminution and converge upon ensuing episodes of floral and fiber development in cotton.

Identification and analysis of oil candidate genes reveals the molecular basis of cottonseed oil accumulation in Gossypium hirsutum L

Based on the integration of QTL-mapping and regulatory network analyses, five high-confidence stable QTL regions, six candidate genes and two microRNAs that potentially affect the cottonseed oil content were discovered. Cottonseed oil is increasingly becoming a promising target for edible oil with its high content of unsaturated fatty acids. In this study, a recombinant inbred line (RIL) cotton population was constructed to detect quantitative trait loci (QTLs) for the cottonseed oil content. A total of 39 QTLs were detected across eight different environments, of which five QTLs were stable. Forty-three candidate genes potentially involved in carbon metabolism, fatty acid synthesis and triacylglycerol biosynthesis processes were further obtained in the stable QTL regions. Transcriptome analysis showed that nineteen of these candidate genes expressed during  the developing cottonseed ovules and may affect the cottonseed oil content. Besides, transcription factor (TF) and microRNA (miRNA) co-regulatory network analyses based on the nineteen candidate genes suggested that six genes, two core miRNAs (ghr-miR2949b and ghr-miR2949c), and one TF GhHSL1 were considered to be closely associated with the cottonseed oil content. Moreover, four vital genes were validated by quantitative real-time PCR (qRT-PCR). These results provide insights into the oil accumulation mechanism in developing cottonseed ovules through the construction of a detailed oil accumulation model.

Global expression dynamics and miRNA evolution profile govern floral/fiber architecture in the modern cotton (Gossypium)

Majority of differentially expressed miRNAs with functional attributes have been recruited independently and parallelly during allopolyploidy followed by the millennia of human selection of both domesticated G. hirsutum and G. barbadense. The genus Gossypium is a marvelous evolutionary model for studying allopolyploidy and morpho-evolution of long-spinnable fibers from the ancestral wild-fuzz. Many genes, transcription factors, and notably, the regulatory miRNAs essentially govern such remarkable modern fiber phenotypes. To comprehend the impact of allopolyploidy on the evolutionary selection of transcriptional dynamicity of key miRNAs, comparative transcriptome profiling of vegetative and fiber tissues of domesticated diploid G. arboreum (A₂) and allopolyploid cotton species G. hirsutum (AD₁), and G. barbadense (AD₂) identified > 300 differentially expressed miRNAs (DEmiRs) within or between corresponding tissues of A₂, AD₁ and AD₂ species. Up to 49% and 32% DEmiRs were up- and down-regulated at fiber initiation stage of AD₁ and AD₂ species, respectively, whereas 50% and 18% DEmiRs were up- and down-regulated at fiber elongation stage of both the allopolyploid species. Interestingly, A-subgenome-specific DEmiRs exhibit expression dominance in the allopolyploid genetic backgrounds. Comparative spatio-temporal expression analyses of AD₁ and AD₂ species discovered that a majority of DEmiRs were recruited independently under millennia of human selection during domestication. Functional annotations of these DEmiRs revealed selection of associated molecular functions such as hormone-signaling, calcium-signaling and reactive oxygen species (ROS) signaling during fiber initiation and elongation. To validate the functional attributes of annotated DEmiRs, we demonstrated for the first time that the target-mimicry-based constitutive diminution of auxin-signaling associated miR167 directly affected the differentiation of floral and fiber tissues of transgenic cotton. These results strongly suggested that the evolutionarily favored DEmiRs including miR167 were involved in the transcriptional regulation of numerous genes during cotton evolution for enhanced fiber-associated agronomic traits.

Identification and characterization of microRNAs in phloem and xylem from ramie (Boehmeria nivea)

Ramie (Boehmeria nivea) is a widely cropped species in southern China due to its high economic value of natural fiber for industry. Development of phloem and xylem is key evidence for generating fiber. However, the MicroRNA (miRNA) profiles of phloem and xylem in ramie have not been reported yet. miRNA belong to a small RNA family which has been recognized as an important regulator for various biological processes. In the present study, we aimed to identify differently expressed miRNAs between phloem and xylem in adult ramie. The results showed that 137 and 122 unique conserved miRNAs were identified from phloem and xylem libraries, respectively. Meanwhile, 4 novel miRNAs were identified from ramie by miRDeep2. Of these miRNAs, 77 conserved miRNAs in ramie were differentially expressed. Among the differentially expressed miRNAs, 44 miRNAs and 33 miRNAs were up-regulated and down-regulated in phloem compared to that in xylem, respectively. The functions of differentially expressed miRNAs were associated with regulating the development and differentiation of phloem and xylem. The present study provides a glance of miRNA profiles for further understanding of miRNA role in ramie development.

Integration of proteomic and transcriptomic profiles reveals multiple levels of genetic regulation of salt tolerance in cotton

BACKGROUND

Salinity is a major abiotic stress that limits upland cotton growth and reduces fibre production worldwide. To reveal genetic regulation via transcript and protein levels after salt stress, we comprehensively analysed the global changes in mRNA, miRNA, and protein profiles in response to salt stress in two contrasting salt-tolerant cotton genotypes.

RESULTS

In the current study, proteomic and mRNA-seq data were combined to reveal that some genes are differentially expressed at both the proteomic and mRNA levels. However, we observed no significant change in mRNA corresponding to most of the strongly differentially abundant proteins. This finding may have resulted from global changes in alternative splicing events and miRNA levels under salt stress conditions. Evidence was provided indicating that several salt stress-responsive proteins can alter miRNAs and modulate alternative splicing events in upland cotton. The results of the stringent screening of the mRNA-seq and proteomic data between the salt-tolerant and salt-sensitive genotypes identified 63 and 85 candidate genes/proteins related to salt tolerance after 4 and 24 h of salt stress, respectively, between the tolerant and sensitive genotype. Finally, we predicted an interaction network comprising 158 genes/proteins and then discovered that two main clusters in the network were composed of ATP synthase (CotAD_74681) and cytochrome oxidase (CotAD_46197) in mitochondria. The results revealed that mitochondria, as important organelles involved in energy metabolism, play an essential role in the synthesis of resistance proteins during the process of salt exposure.

CONCLUSION

We provided a plausible schematic for the systematic salt tolerance model; this schematic reveals multiple levels of gene regulation in response to salt stress in cotton and provides a list of salt tolerance-related genes/proteins. The information here will facilitate candidate gene discovery and molecular marker development for salt tolerance breeding in cotton.

MicroRNA expression profiles during cotton (Gossypium hirsutum L) fiber early development

The role of microRNAs (miRNAs) during cotton fiber development remains unclear. Here, a total of 54 miRNAs belonging to 39 families were selected to characterize miRNA regulatory mechanism in eight different fiber development stages in upland cotton cv BM-1. Among 54 miRNAs, 18 miRNAs were involved in cotton fiber initiation and eight miRNAs were related to fiber elongation and secondary wall biosynthesis. Additionally, 3,576 protein-coding genes were candidate target genes of these miRNAs, which are potentially involved in cotton fiber development. We also investigated the regulatory network of miRNAs and corresponding targets in fiber initiation and elongation, and secondary wall formation. Our Gene Ontology-based term classification and KEGG-based pathway enrichment analyses showed that the miRNA targets covered 220 biological processes, 67 molecular functions, 45 cellular components, and 10 KEGG pathways. Three of ten KEGG pathways were involved in lignan synthesis, cell elongation, and fatty acid biosynthesis, all of which have important roles in fiber development. Overall, our study shows the potential regulatory roles of miRNAs in cotton fiber development and the importance of miRNAs in regulating different cell types. This is helpful to design miRNA-based biotechnology for improving fiber quality and yield.

Differential expression of microRNAs during fiber development between fuzzless-lintless mutant and its wild-type allotetraploid cotton

Cotton is one of the most important textile crops but little is known how microRNAs regulate cotton fiber development. Using a well-studied cotton fiberless mutant Xu-142-fl, we compared 54 miRNAs for their expression between fiberless mutant and its wildtype. In wildtype Xu-142, 26 miRNAs are involved in cotton fiber initiation and 48 miRNAs are related to primary wall synthesis and secondary wall thickening. Thirty three miRNAs showed different expression in fiber initiation between Xu-142 and Xu-142-fl. These miRNAs potentially target 723 protein-coding genes, including transcription factors, such as MYB, ARF, and LRR. ARF18 was newly predicted targets of miR160a, and miR160a was expressed at higher level in −2DPA of Xu-142-fl compared with Xu-142. Furthermore, the result of Gene Ontology-based term classification (GO), EuKaryotic Orthologous Groups (KOG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis shows that miRNA targets were classified to 222 biological processes, 64 cellular component and 42 molecular functions, enriched in 22 KOG groups, and classified into 28 pathways. Together, our study provides evidence for better understanding of miRNA regulatory roles in the process of fiber development, which is helpful to increase fiber yield and improve fiber quality.

Small RNA sequencing and degradome analysis of developing fibers of short fiber mutants Ligon-lintles-1 (Li 1 ) and -2 (Li 2 ) revealed a role for miRNAs and their targets in cotton fiber elongation

BACKGROUND

The length of cotton fiber is an important agronomic trait that directly affects the quality of yarn and fabric. Understanding the molecular basis of fiber elongation would provide a means for improvement of fiber length. Ligon-lintless-1 (Li 1 ) and -2 (Li 2 ) are monogenic and dominant mutations that result in an extreme reduction in the length of lint fiber on mature seeds. In a near-isogenic state with wild type cotton these two short fiber mutants provide an effective model system to study the mechanisms of fiber elongation. Plant miRNAs regulate many aspects of growth and development. However, the mechanism underlying the miRNA-mediated regulation of fiber development is largely unknown.

RESULTS

Small RNA libraries constructed from developing fiber cells of the short fiber mutants Li 1 and Li 2 and their near-isogenic wild type lines were sequenced. We identified 24 conservative and 147 novel miRNA families with targets that were detected through degradome sequencing. The distribution of the target genes into functional categories revealed the largest set of genes were transcription factors. Expression profiles of 20 miRNAs were examined across a fiber developmental time course in wild type and short fiber mutations. We conducted correlation analysis between miRNA transcript abundance and the length of fiber for 11 diverse Upland cotton lines. The expression patterns of 4 miRNAs revealed significant negative correlation with fiber lengths of 11 cotton lines.

CONCLUSIONS

Our results suggested that the mutations have changed the regulation of miRNAs expression during fiber development. Further investigations of differentially expressed miRNAs in the Li 1 and Li 2 mutants will contribute to better understanding of the regulatory mechanisms of cotton fiber development. Four miRNAs negatively correlated with fiber length are good candidates for further investigations of miRNA regulation of important genotype dependent fiber traits. Thus, our results will contribute to further studies on the role of miRNAs in cotton fiber development and will provide a tool for fiber improvement through molecular breeding.

RNA Interference for Functional Genomics and Improvement of Cotton (Gossypium sp.)

RNA interference (RNAi), is a powerful new technology in the discovery of genetic sequence functions, and has become a valuable tool for functional genomics of cotton (Gossypium sp.). The rapid adoption of RNAi has replaced previous antisense technology. RNAi has aided in the discovery of function and biological roles of many key cotton genes involved in fiber development, fertility and somatic embryogenesis, resistance to important biotic and abiotic stresses, and oil and seed quality improvements as well as the key agronomic traits including yield and maturity. Here, we have comparatively reviewed seminal research efforts in previously used antisense approaches and currently applied breakthrough RNAi studies in cotton, analyzing developed RNAi methodologies, achievements, limitations, and future needs in functional characterizations of cotton genes. We also highlighted needed efforts in the development of RNAi-based cotton cultivars, and their safety and risk assessment, small and large-scale field trials, and commercialization.

MicroRNAs in cotton: an open world needs more exploration

This paper reviews the progress and current problems in the field of cotton microRNAs. Cotton is not only one of the most important crops in terms of fiber usage and economic value, but also a model species for investigating cell wall and cellulose biosynthesis as well plant polyploidization. Compared with model plant species, such as Arabidopsis and rice, the research on cotton microRNAs (miRNAs) is lagging, although great progress has been made in the past decade. Since the first reports on identifying miRNAs in cotton in 2007, hundreds of miRNAs have been identified using an in silico comparative genome-based approach and direct cloning. Next-generation deep sequencing has opened the door for cotton miRNA research. In cotton, miRNAs are associated with many biological and metabolic processes, including fiber initiation and development, floral development, embryogenesis, and response to biotic and abiotic stresses. However, the majority of current research is focused on miRNA identification. Although several targets have been predicted using computational approaches and degradome sequencing, more functional studies should be performed in the next couple of years to elucidate the roles of miRNAs in cotton fiber development and response to different environmental stresses using transgenic technology. This paper reviews the history, identification, and function of cotton miRNAs as well as future directions for this research.

Small RNA sequencing identifies miRNA roles in ovule and fibre development

MicroRNAs (miRNAs) have been found to be differentially expressed during cotton fibre development. However, which specific miRNAs and how they are involved in fibre development is unclear. Here, using deep sequencing, 65 conserved miRNA families were identified and 32 families were differentially expressed between leaf and ovule. At least 40 miRNAs were either leaf or ovule specific, whereas 62 miRNAs were shared in both leaf and ovule. qRT-PCR confirmed these miRNAs were differentially expressed during fibre early development. A total of 820 genes were potentially targeted by the identified miRNAs, whose functions are involved in a series of biological processes including fibre development, metabolism and signal transduction. Many predicted miRNA-target pairs were subsequently validated by degradome sequencing analysis. GO and KEGG analyses showed that the identified miRNAs and their targets were classified to 1027 GO terms including 568 biological processes, 324 molecular functions and 135 cellular components and were enriched to 78 KEGG pathways. At least seven unique miRNAs participate in trichome regulatory interaction network. Eleven trans-acting siRNA (tasiRNA) candidate genes were also identified in cotton. One has never been found in other plant species and two of them were derived from MYB and ARF, both of which play important roles in cotton fibre development. Sixteen genes were predicted to be tasiRNA targets, including sucrose synthase and MYB2. Together, this study discovered new miRNAs in cotton and offered evidences that miRNAs play important roles in cotton ovule/fibre development. The identification of tasiRNA genes and their targets broadens our understanding of the complicated regulatory mechanism of miRNAs in cotton.

Bioinformatics analysis of small RNAs in pima (Gossypium barbadense L.)

Small RNAs (sRNAs) are ~20 to 24 nucleotide single-stranded RNAs that play crucial roles in regulation of gene expression. In plants, sRNAs are classified into microRNAs (miRNAs), repeat-associated siRNAs (ra-siRNAs), phased siRNAs (pha-siRNAs), cis and trans natural antisense transcript siRNAs (cis- and trans-nat siRNAs). Pima (Gossypium barbadense L.) is one of the most economically important fiber crops, producing the best and longest spinnable fiber. Although some miRNAs are profiled in Pima, little is known about siRNAs, the largest subclass of plant sRNAs. In order to profile these gene regulators in Pima, a comprehensive analysis of sRNAs was conducted by mining publicly available sRNA data, leading to identification of 678 miRNAs, 3,559,126 ra-siRNAs, 627 pha-siRNAs, 136,600 cis-nat siRNAs and 79,994 trans-nat siRNAs. The 678 miRNAs, belonging to 98 conserved and 402 lineage-specific families, were produced from 2,138 precursors, of which 297 arose from introns, exons, or intron/UTR-exon junctions of protein-coding genes. Ra-siRNAs were produced from various repeat loci, while most (97%) were yielded from retrotransposons, especially LTRs (long terminal repeats). The genes encoding auxin-signaling-related proteins, NBS-LRRs and transcription factors were major sources of pha-siRNAs, while two conserved TAS3 homologs were found as well. Most cis-NATs in Pima overlapped in enclosed and convergent orientations, while a few hybridized in divergent and coincided orientations. Most cis- and trans-nat siRNAs were produced from overlapping regions. Additionally, characteristics of length and the 5’-first nucleotide of each sRNA class were analyzed as well. Results in this study created a valuable molecular resource that would facilitate studies on mechanism of controlling gene expression.

Deep sequencing reveals important roles of microRNAs in response to drought and salinity stress in cotton

Drought and salinity are two major environmental factors adversely affecting plant growth and productivity. However, the regulatory mechanism is unknown. In this study, the potential roles of small regulatory microRNAs (miRNAs) in cotton response to those stresses were investigated. Using next-generation deep sequencing, a total of 337 miRNAs with precursors were identified, comprising 289 known miRNAs and 48 novel miRNAs. Of these miRNAs, 155 miRNAs were expressed differentially. Target prediction, Gene Ontology (GO)-based functional classification, and Kyoto Encyclopedia of Genes and Genomes (KEGG)-based functional enrichment show that these miRNAs might play roles in response to salinity and drought stresses through targeting a series of stress-related genes. Degradome sequencing analysis showed that at least 55 predicted target genes were further validated to be regulated by 60 miRNAs. CitationRank-based literature mining was employed to determinhe the importance of genes related to drought and salinity stress. The NAC, MYB, and MAPK families were ranked top under the context of drought and salinity, indicating their important roles for the plant to combat drought and salinity stress. According to target prediction, a series of cotton miRNAs are associated with these top-ranked genes, including miR164, miR172, miR396, miR1520, miR6158, ghr-n24, ghr-n56, and ghr-n59. Interestingly, 163 cotton miRNAs were also identified to target 210 genes that are important in fibre development. These results will contribute to cotton stress-resistant breeding as well as understanding fibre development.