Enrichment of a set of microRNAs during the cotton fiber development

Computational identification and experimental validation of novel Saccharum officinarum microRNAs along with their targets through RT-PCR approach

Various metabolic and cell signaling processes impact the functions of sugarcane plant cells. MicroRNAs (miRNAs) play critical regulatory roles in enhancing yield and providing protection against various stressors. This study seeks to identify and partially characterize several novel miRNAs in sugarcane using in silico tools, while also offering a preliminary assessment of their functions. This was accomplished by predicting novel conserved miRNAs in sugarcane plants using a variety of genomics-based techniques like BLASTn, MFOLD, psRNA Target, sequence logo, Weblogo, primer-3, etc. and annotated using miRBase and NCBI. For validation, RT-PCR method was used along with agarose gel. After the preparation of fourteen randomly chosen primers, they were validated by RT-PCR. Accordingly, they contain fifty specific targeted proteins with substantial targets in the structural, transcriptional protein, etc. Furthermore, the sof-miR5025a directs the heat repeat protein while the voltage-dependent anion is governed by sof-miR8005a. Similarly, the sof-miR7768b and sof-miR6249b monitor the pathogenesis-related protein and zinc finger, C2H2 type protein, which assist as transcription factors. Thus, the novel sugarcane miRNAs target a wide range of important genes help regulate the environment for sugarcane to generate a higher-quality crop.

Differential miRNA expression and regulatory mechanisms in pigmentation and fiber development of white and brown cotton (Gossypium hirsutum)

Cotton (Gossypium hirsutum) is a major global natural fiber crop used in the textile industry. Although white colored cotton remains the most popular form of cultivated cotton, colored varieties could replace chemically dyed fibers and provide more environmental friendly alternatives. In order to investigate the role of miRNAs in fiber color, we selected white and brown cotton varieties for comparative investigations. Through small RNA sequencing, a number of known miRNA families were discovered (74 in white cotton and 61 in brown cotton, with 44 shared) in which 11 miRNA families were significantly elevated in brown cotton variety. Functional enrichment and network analysis of target genes of these miRNAs revealed their regulatory role in secondary metabolite biosynthesis pathway, particularly the flavonoids pathway, which are known to be associated with fiber coloration. Pigmentation and developmental-related miRNA members such as miR396e-5p, miR167l, and miR1446 were also significantly enriched. Real-time PCR results suggest the regulatory role of miRNAs in these two cotton varieties. Furthermore, 30 and 25 novel miRNAs were also identified in white and brown cotton, respectively. Our findings also show miRNAs associated with fiber coloration and development through the intricate networks of miRNA and targets. Understanding these systems may provide novel insights on improving the fiber color and quality.

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.

Co-localization and analysis of miR477b with fiber length quantitative trait loci in cotton

Cotton is an important cash crop for the textile industry. However, the understanding of natural genetic variation of fiber elongation in relation to miRNA is lacking. A miRNA gene (miR477b) was found to co-localize with a previously mapped fiber length (FL) quantitative trait locus (QTL). The miR477b was differentially expressed during fiber elongation between two backcross inbred lines (BILs) differing in FL and its precursor sequences. Bioinformatics and qRT-PCR analysis were further used to analyse the miRNA genes, which could produce mature miR477b. Cotton plants with virus-induced gene silencing (VIGS) constructs to over-express the allele of miR477b from the BIL with longer fibers had significantly longer fibers as compared with negative control plants, while the VIGS plants with suppressed miRNA expression had significantly shorter fibers. The expression level of the target gene (DELLA) and related genes (RDL1 and EXPA1 for DELLA through HOX3 protein) in the two BILs and/or the VIGS plants were generally congruent, as expected. This report represents one of the first comprehensive studies to integrate QTL linkage mapping and physical mapping of small RNAs with both small and mRNA transcriptome analysis, followed by VIGS, to identify candidate small RNA genes affecting the natural variation of fiber elongation in cotton.

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 of conserved miRNAs and their targets in Jatropha curcas: an in silico approach

BACKGROUND

MicroRNAs (miRNAs) are small endogenous RNAs with an approximate length of 18-22 nucleotides and involved in the regulation of gene expression in transcriptional or post-transcriptional levels. They were found to be associated with leaf morphogenesis, flowering time, vegetative phase change, and response to environmental cues in plants, where they act as a critical regulatory factor. The nature of high conservancy of plant miRNAs within the plant species made it possible to detect the conserved miRNAs by computational approaches. Expressed Sequence Tags (EST) based comparative genomic approaches provide advantages over wet lab approaches as it is convenient, easy to carry out and less time consuming. EST-based in silico approach can unravel new conserved miRNAs in plants, even when the complete genome sequence is not available.

RESULTS

To identify the novel miRNAs, a total of 46,865 ESTs from Jatropha curcas were searched for homology to all available 6746 mature miRNAs of plant eudicotyledons. Finally, we ended up with 12 novel miRNAs in Jatropha that range from 18 to 19 nucleotides where their respective precursor miRNAs had 54.11-71.76% (A + U) content. The putative miRNAs belong to 12 individual miRNA family and most of them have higher (A + U) content ranging from 47.36 to 77.77% than their respective miRNA homologs. Many of the target genes by the newly identified miRNAs were associated with plant growth and development, stress response, defense and hormone signaling, and oil synthesis pathways.

CONCLUSION

These findings have the potential to speed up miRNA identification and expand our understanding of miRNA functions in J. curcas.

Analysis of the MIR396 gene family and the role of MIR396b in regulating fiber length in cotton

Cotton fiber is one of the most important natural raw materials in the world textile industry. Improving fiber yield and quality has always been the main goal. MicroRNAs, as typical small noncoding RNAs, could affect fiber length during different stages of fiber development. Based on differentially expressed microRNA in the two interspecific backcross inbred lines (BILs) with a significant difference in fiber length, we identified the miR396 gene family in the two tetraploid cotton genomes and found MIR396b_D13 as the functional precursor to produce mature miR396 during the fiber elongation stage. Among 46 target genes regulated by miR396b, the GROWTH-REGULATING FACTOR 5 gene (GRF5, Gh_A10G0492) had a differential expression level in the two BILs during fiber elongation stage. The expression patterns indicated that the miR396b-GRF5 regulatory module has a critical role in fiber development. Furthermore, virus-induced gene silencing (VIGS) of miR396b significantly produced longer fiber than the wild type, and the expression level of GRF5 showed the reverse trends of the miR396b expression level. The analysis of co-expression network for the GRF5 gene suggested that a cytochrome P450 gene functions as an allene oxide synthase (Gh_D06G0089, AOS), which plays a critical role in jasmonate biosynthetic pathway. In conclusion, our results revealed that the miR396b-GRF5 module has a critical role in fiber development. These findings provide a molecular foundation for fiber quality improvement in the future.

Genome-wide analysis of IQD proteins and ectopic expression of watermelon ClIQD24 in tomato suggests its important role in regulating fruit shape

The plant-specific IQ67 domain (IQD) is the largest class of calmodulin targets found in plants, and plays an important role in many biological processes, especially fruit development processes. However, the functional role of IQD proteins in the development of watermelon (Citrullus lanatus) shape remains unknown, as the IQD protein family in watermelon has not been systematically characterized. Herein, we elucidated the gene structures, chromosomal locations, evolutionary divergence, and functions of 35 IQD genes in the watermelon genome. The transcript profiles and quantitative real-time PCR analysis at different stages of fruit development showed that the ClIQD24 gene was highly expressed on 0 days after pollination. Furthermore, we found that the ectopic overexpression of ClIQD24 promoted tomato fruit elongation, thereby revealing the significance of ClIQD24 in the progression of watermelon shape. Our study will serve as a reference for further investigations on the molecular mechanisms underlying watermelon fruit shape formation.

Genome-Wide Identification and Characterization of CPR5 Genes in Gossypium Reveals Their Potential Role in Trichome Development

Trichomes protect plants against insects, microbes, herbivores, and abiotic damages and assist seed dispersal. The function of CPR5 genes have been found to be involved in the trichome development but the research on the underlying genetic and molecular mechanisms are extremely limited. Herein, genome wide identification and characterization of CPR5 genes was performed. In total, 26 CPR5 family members were identified in Gossypium species. Phylogenetic analysis, structural characteristics, and synteny analysis of CPR5s showed the conserved evolution relationships of CPR5. The promoter analysis of CPR5 genes revealed hormone, stress, and development-related cis-elements. Gene ontology (GO) enrichment analysis showed that the CPR5 genes were largely related to biological regulation, developmental process, multicellular organismal process. Protein-protein interaction analysis predicted several trichome development related proteins (SIM, LGO, and GRL) directly interacting with CPR5 genes. Further, nine putative Gossypium-miRNAs were also identified, targeting Gossypium CPR5 genes. RNA-Seq data of G. arboreum (with trichomes) and G. herbaceum (with no trichomes) was used to perform the co-expression network analysis. GheCPR5.1 was identified as a hub gene in a co-expression network analysis. RT-qPCR of GheCPR5.1 gene in different tissues suggests that this gene has higher expressions in the petiole and might be a key candidate involved in the trichome development. Virus induced gene silencing of GheCPR5.1 (Ghe02G17590) confirms its role in trichome development and elongation. Current results provide proofs of the possible role of CPR5 genes and provide preliminary information for further studies of GheCPR5.1 functions in trichome development.

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.

In silico identification of conserved miRNAs in the genome of fibre biogenesis crop Corchorus capsularis

Corchorus capsularis, commonly known as jute occupies the leading position in the production of natural fibre alongside lower environmental threat. Small noncoding ~21 to 24 nucleotides long microRNAs play significant roles in regulating the gene expression as well as different functions in cellular growth and development. Here, the study adopted a comprehensive in silico approach to identify and characterize the conserved miRNAs in the genome of C. capsularis including functional annotation of specific gene targets. Expressed Sequence Tags (ESTs) based homology search of 3350 known miRNAs of dicotyledons were allowed against 763 non-redundant ESTs of jute genome, resulted in the prediction of 5 potential miRNA candidates belonging five different miRNA families (miR1536, miR9567-3p, miR4391, miR11300, and miR8689). The putative miRNAs were composed of 18 nucleotides having a range of -0.49 to -1.56 MFEI values and 55%–61% of (A + U) content in their pre-miRNAs. A total of 1052 gene targets of putative miRNAs were identified and their functions were extensively analyzed. Most of the gene targets were involved in plant growth, cell cycle regulation, organelle synthesis, developmental process and environmental responses. Five gene targets, namely, NAC Domain Containing Protein, WRKY DNA binding protein, 3-dehydroquinate synthase, S-adenosyl-L-Met–dependent methyl transferase and Vascular-related NAC-Domain were found to be involved in the lignin biosynthesis, phenylpropanoid pathways and secondary wall formation. The present study might accelerate the more miRNA discovery, strengthening the complete understanding of miRNAs association in the cellular basis of lignin biosynthesis towards the production of high standard jute products.

A Malvaceae-specific miRNA targeting the newly duplicated GaZIP1L to regulate Zn2+ ion transporter capacity in cotton ovules

MicroRNAs (miRNAs) play critical roles in regulating gene expression in plants, yet their functions underlying cultivated diploid Gossypium arboreum cotton ovule development are largely unknown. Here, we acquired small RNA profiles from G. arboreum ovules and fibers collected at different growth stages, and identified 46 novel miRNAs that accounted for 23.7% of all miRNAs in G. arboreum reported in the latest plant sRNA database. Through analysis of 84 (including 38 conserved) differentially expressed G. arboreum miRNAs, we detected 215 putative protein-coding genes in 26 biological processes as their potential targets. A Malvaceae-specific novel miRNA named gar-miRN44 was found to likely regulate cotton ovule growth by targeting to a newly duplicated Zn²⁺ ion transporter gene GaZIP1L. During cotton ovule development, gar-miRN44 transcript level decreased sharply after 10 to 15 days post-anthesis (DPA), while that of the GaZIP1L increased significantly, with a concomitant increase of Zn²⁺ ion concentration in late ovule developmental stages. Molecular dynamics simulation and ion absorption analysis showed that GaZIP1L has stronger Zn²⁺ ion binding ability than the original GaZIP1, indicating that the newly evolved GaZIP1L may be more suitable for maintaining high Zn²⁺ ion transport capacity that is likely required for cotton ovule growth via enhanced cellulose synthase activities. Our systematic miRNA profiling in G. arboreum and characterization of gar-miRN44 not only contribute to the understanding of miRNA function in cotton, but also provide potential targets for plant breeding.

Identification, characterization and expression profiling of circular RNAs in the early cotton fiber developmental stages

Circular RNA is one of the endogenous non-coding RNAs with a covalently closed loop structure and largely involved in regulating gene expression. However, the abundance of circular RNAs and their regulatory functions during the early stages of fiber development are still not known. In this work, we conducted high-throughput sequencing of the Ligonlintless-1 and its wild-type at 0 DPA, 8 DPA and stem. A total of 2811 circular RNAs were identified and unevenly distributed across cotton chromosomes. We found 34, 142 and 27 circular RNAs were differentially expressed between Ligonlintless-1 and wild-type at 0 DPA, 8 DPA and stem, respectively. Both circular RNA-microRNA-mRNA network and MeJA treatment results in Ligonlintless-1 and wild-type might provide a strong indication of four circular RNAs and ghr_miR169b being important biological molecular associating with fiber development. The results provide new insight into the putative molecular function of circular RNAs in the regulation of fiber development.

Genome-Wide Characterization and Analysis of CIPK Gene Family in Two Cultivated Allopolyploid Cotton Species: Sequence Variation, Association with Seed Oil Content, and the Role of GhCIPK6

Calcineurin B-like protein-interacting protein kinases (CIPKs), as key regulators, play an important role in plant growth and development and the response to various stresses. In the present study, we identified 80 and 78 CIPK genes in the Gossypium hirsutum and G. barbadense, respectively. The phylogenetic and gene structure analysis divided the cotton CIPK genes into five groups which were classified into an exon-rich clade and an exon-poor clade. A synteny analysis showed that segmental duplication contributed to the expansion of Gossypium CIPK gene family, and purifying selection played a major role in the evolution of the gene family in cotton. Analyses of expression profiles showed that GhCIPK genes had temporal and spatial specificity and could be induced by various abiotic stresses. Fourteen GhCIPK genes were found to contain 17 non-synonymous single nucleotide polymorphisms (SNPs) and co-localized with oil or protein content quantitative trait loci (QTLs). Additionally, five SNPs from four GhCIPKs were found to be significantly associated with oil content in one of the three field tests. Although most GhCIPK genes were not associated with natural variations in cotton oil content, the overexpression of the GhCIPK6 gene reduced the oil content and increased C18:1 and C18:1+C18:1d6 in transgenic cotton as compared to wild-type plants. In addition, we predicted the potential molecular regulatory mechanisms of the GhCIPK genes. In brief, these results enhance our understanding of the roles of CIPK genes in oil synthesis and stress responses.

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.

Analysis of the MIR160 gene family and the role of MIR160a_A05 in regulating fiber length in cotton

The MIR160 family in Gossypium hirsutum and G. barbadense was characterized, and miR160a_A05 was found to increase cotton-fiber length by downregulating its target gene (ARF17) and several GH3 genes. Cotton fiber is the most important raw material for the textile industry. MicroRNAs are involved in regulating cotton-fiber development, but a role in fiber elongation has not been demonstrated. In this study, miR160a was found to be differentially expressed in elongating fibers between two interspecific (between Gossypium hirsutum and G. barbadense) backcross inbred lines (BILs) with different fiber lengths. The gene MIR160 colocalized with a previously mapped fiber-length quantitative trait locus. Its target gene ARF17 was differentially expressed between the two BILs during fiber elongation, but in the inverse fashion. Bioinformatics was used to analyze the MIR160 family in both G. hirsutum and G. barbadense. Moreover, qRT-PCR analysis identified MIR160a as the functional MIR160 gene encoding the miR160a precursor during fiber elongation. Using virus-induced gene silencing and overexpression, overexpressed MIR160a_A05 resulted in significantly longer fibers compared with wild type, whereas suppression of miR160 resulted in significantly shorter fibers. Expression levels of the target gene auxin-response factor 17 (ARF17) and related genes GH3 in the two BILs and/or the virus-infected plants demonstrated similar changes in response to modulation of miR160a level. Finally, overexpression or suppression of miR160 increased or decreased, respectively, the cellular level of indole-3-acetic acid, which is involved in fiber elongation. These results describe a specific regulatory mechanism for fiber elongation in cotton that can be utilized for future crop improvement.

Comparative Proteomic Analysis of Molecular Differences between Leaves of Wild-Type Upland Cotton and Its Fuzzless-Lintless Mutant

Fuzzless-lintless mutant (fl) ovules of upland cotton have been used to investigate cotton fiber development for decades. However, the molecular differences of green tissues between fl and wild-type (WT) cotton were barely reported. Here, we found that gossypol content, the most important secondary metabolite of cotton leaves, was higher in Gossypium hirsutum L. cv Xuzhou-142 (Xu142) WT than in fl. Then, we performed comparative proteomic analysis of the leaves from Xu142 WT and its fl. A total of 4506 proteins were identified, of which 103 and 164 appeared to be WT- and fl-specific, respectively. In the 4239 common-expressed proteins, 80 and 74 were preferentially accumulated in WT and fl, respectively. Pathway enrichment analysis and protein–protein interaction network analysis of both variety-specific and differential abundant proteins showed that secondary metabolism and chloroplast-related pathways were significantly enriched. Quantitative real-time PCR confirmed that the expression levels of 12 out of 16 selected genes from representative pathways were consistent with their protein accumulation patterns. Further analyses showed that the content of chlorophyll a in WT, but not chlorophyll b, was significantly increased compared to fl. This work provides the leaf proteome profiles of Xu142 and its fl mutant, indicating the necessity of further investigation of molecular differences between WT and fl leaves.

Role of MicroRNAs and small RNAs in regulation of developmental processes and agronomic traits in Gossypium species

Small RNAs (sRNAs) are short, non-coding, 17-24 nucleotides long RNA molecules that play vital roles in regulating gene expression in every known organism investigated to date including cotton (Gossypium ssp.). These tiny RNA molecules target diverse categories of genes from different bioliogical and metabolic processes and have been reported in the three domains of life. Small RNAs, including miRNAs, are involved in ovule and fiber development, biotic and abiotic stresses, fertility, and other biochemical processes in cotton species. Also, sRNAs are the critical components in RNA interference pathway. In this article, we have reviewed the research efforts related to the isolation and characterization of miRNAs using molecular and genomic approaches. The progress made in understanding the functional roles of miRNAs in regulation, alteration, and inactivation of fundamental plant processes and traits of importance in cotton are presented here.

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

BACKGROUND

Long non-coding RNAs (LncRNAs) are part of genes, which are not translated into proteins and play a vital role in plant growth and development. Nevertheless, the presence of LncRNAs and how they functions in Ligon-lintless-1 mutant during the early cessation of cotton fiber development are still not well understood. In order to investigate the function of LncRNAs in cotton fiber development, it is necessary and important to identify LncRNAs and their potential roles in fiber cell development.

RESULTS

In this work, we identified 18,333 LncRNAs, with the proportion of long intergenic noncoding RNAs (LincRNAs) (91.5%) and anti-sense LncRNAs (8.5%), all transcribed from Ligon-lintless-1 (Li1) and wild-type (WT). Expression differences were detected between Ligon-lintless-1 and wild-type at 0 and 8 DPA (day post anthesis). Pathway analysis and Gene Ontology based on differentially expressed LncRNAs on target genes, indicated fatty acid biosynthesis and fatty acid elongation being integral to lack of fiber in mutant cotton. The result of RNA-seq and RT-qPCR clearly singles out two potential LncRNAs, LNC_001237 and LNC_017085, to be highly down-regulated in the mutant cotton. The two LncRNAs were found to be destabilized or repressed by ghr-miR2950. Both RNA-seq analysis and RT-qPCR results in Ligon-lintless-1 mutant and wild-type may provide strong evidence of LNC_001237, LNC_017085 and ghr-miR2950 being integral molecular elements participating in various pathways of cotton fiber development.

CONCLUSION

The results of this study provide fundamental evidence for the better understanding of LncRNAs regulatory role in the molecular pathways governing cotton fiber development. Further research on designing and transforming LncRNAs will help not only in the understanding of their functions but will also in the improvement of fiber quality.

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.

Overexpression of miRNA in Cotton via Agrobacterium-Mediated Transformation

microRNAs (miRNAs) are an extensive class of newly identified endogenous small regulatory molecules. Many studies show that miRNAs play a critical role in almost all biological and metabolic progress through targeting protein-coding genes for mRNA cleavage or translation inhibition. Many miRNAs are also identified from cotton using computational and/or experimental approaches, including the next-generation deep sequencing technology. However, their functions are unclear. In this chapter, we describe a detailed method for overexpressing miRNA, miR156 as an example, in cotton using Agrobacterium-mediated genetic transformation. This provides an approach to investigate the function and regulatory mechanism of miRNAs in cotton.

The miRNAome of ramie (Boehmeria nivea L.): identification, expression, and potential roles of novel microRNAs in regulation of cadmium stress response

BACKGROUND

MicroRNAs (miRNAs) regulate numerous crucial abiotic stress processes in plants. However, information is limited on their involvement in cadmium (Cd) stress response and tolerance mechanisms in plants, including ramie (Boehmeria nivea L.) that produces a number of economic valuable as an important natural fibre crop and an ideal crop for Cd pollution remediation.

RESULTS

Four small RNA libraries of Cd-stressed and non-stressed leaves and roots of ramie were constructed. Using small RNA-sequencing, 73 novel miRNAs were identified. Genome-wide expression analysis revealed that a set of miRNAs was differentially regulated in response to Cd stress. In silico target prediction identified 426 potential miRNA targets that include several uptake or transport factors for heavy metal ions. The reliability of small RNA sequencing and the relationship between the expression levels of miRNAs and their target genes were confirmed by quantitative PCR (q-PCR). We showed that the expression patterns of miRNAs obtained by q-PCR were consistent with those obtained from small RNA sequencing. Moreover, we demonstrated that the expression of six randomly selected target genes was inversely related to that of their corresponding miRNAs, indicating that the miRNAs regulate Cd stress response in ramie.

CONCLUSIONS

This study enriches the number of Cd-responsive miRNAs and lays a foundation for the elucidation of the miRNA-mediated regulatory mechanism in ramie during Cd stress.

In Silico Identification of Conserved MiRNAs from Physcomitrella patens ESTs and their Target Characterization

Background:

MicroRNAs (miRNAs) are groups of small non-protein-coding endogenous single stranded RNAs with approximately 18-24 nucleotides in length. High evolutionary sequence conservation of miRNAs among plant species and availability of powerful computational tools allow identification of new orthologs and paralogs.

Methods:

New conserved miRNAs in P. patens were found by EST-based homology search approaches. All candidates were screened according to a series of miRNA filtering criteria. Unigene, DFCI Gene Index (PpspGI) databases and psRNATarget algorithm were applied to identify target transcripts using P. patens putative conserved miRNA sequences.

Results:

Nineteen conserved P. patens miRNAs were identified. The sequences were homologous to known reference plant mature miRNA from 10 miRNA families. They could be folded into the typical miRNA secondary structures. RepeatMasker algorithm demonstrated that ppt-miR2919e and pptmiR1533 had simple sequence repeats in their sequences. Target sites (49 genes) were identified for 7 out of 19 miRNAs. GO and KEGG analysis of targets indicated the involvement of some in important multiple biological and metabolic processes.

Conclusion:

The majority of the registered miRNAs in databases were predicted by computational approaches while many more have remained unknown. Due to the conserved nature of miRNAs in plant species from closely to distantly related, homology search-based approaches between plants species could lead to the identification of novel miRNAs in other plant species providing baseline information for further search about the biological functions and evolution of miRNAs.

Identification of Expressed miRNAs in Human Rheumatoid Arthritis Using Computational Approach - Discovery of a New miR-7167 from Human

BACKGROUND

Rheumatoid Arthritis (RA) is a chronic inflammatory and autoimmune disease leading to bones and joints destruction. It is one of the major causes of lifetime disability and mortality among humans in the developing and developed countries. It was evident that epigenetic dysregulation is related to the pathogenesis of RA. MicroRNAs (miRNAs) are small non-coding RNAs that are epigenetic regulators for diverse biological processes and also provided novel molecular insights in the formation of arthritis.

OBJECTIVE

The influences of miRNAs in the alteration of gene regulation during the pathogenesis of arthritis were exposed in recent years.

METHOD

The computational approach to identify miRNA through EST-based homology is more powerful, economical and time-efficient. In this study, we applied EST-based homology search to identify miRNAs responsible for the development of arthritis in human beings.

RESULTS

Our study on 36519 ESTs in human RA condition revealed the expression of four miRNAs, HSA-miR-198, HSA-miR-4647, has-miR-7167-5p and has-miR-7167-3p. The present study is the first report about has-miR-7167 that was homologous to Macaca mulatta.

CONCLUSION

The predicted targets of these identified miRNAs revealed many biological functions in the pathogenesis of RA. Further elaborated studies on these miRNAs will help to understand their function in the development of RA and the use of miRNAs as therapeutic targets in the future.

New insights into tomato microRNAs

Cultivated tomato, Solanum lycopersicum, is one of the most common fruits in the global food industry. Together with the wild tomato Solanum pennellii, it is widely used for developing better cultivars. MicroRNAs affect mRNA regulation, inhibiting its translation and/or promoting its degradation. Important proteins involved in these processes are ARGONAUTE and DICER. This study aimed to identify and characterize the genes involved in the miRNA processing pathway, miRNA molecules and target genes in both species. We validated the presence of pathway genes and miRNA in different NGS libraries and 6 miRNA families using quantitative RT-PCR. We identified 71 putative proteins in S. lycopersicum and 108 in S. pennellii likely involved in small RNAs processing. Of these, 29 and 32 participate in miRNA processing pathways, respectively. We identified 343 mature miRNAs, 226 pre-miRNAs in 87 families, including 192 miRNAs, which were not previously identified, belonging to 38 new families in S. lycopersicum. In S. pennellii, we found 388 mature miRNAs and 234 pre-miRNAs contained in 85 families. All miRNAs found in S. pennellii were unpublished, being identified for the first time in our study. Furthermore, we identified 2471 and 3462 different miRNA target in S. lycopersicum and S. pennellii, respectively.

Genome-wide characterization, identification, and expression analysis of the WD40 protein family in cotton

WD40 repeat proteins are largely distributed across the plant kingdom and play an important role in diverse biological activities. In this work, we performed genome-wide identification, characterization, and expression level analysis of WD40 genes in cotton. A total of 579, 318, and 313 WD40 genes were found in Gossypium hirsutum, G. arboreum, and G. raimondii, respectively. Based on phylogenetic tree analyses, WD40 genes were divided into 11 groups with high similarities in exon/intron features and protein domains within the group. Expression analysis of WD40 genes showed differential expression at different stages of cotton fiber development (0 and 8 DPA) and cotton stem. A number of miRNAs were identified to target WD40 genes that are significantly involved in cotton fiber development during the initiation and elongation stages. These include miR156, miR160, miR162, miR164, miR166, miR167, miR169, miR171, miR172, miR393, miR396, miR398, miR2950, and miR7505. The findings provide a stronger indication of WD40 gene function and their involvement in the regulation of cotton fiber development during the initiation and elongation stages.

Identification of microRNAs involved in drought stress responses in early-maturing cotton by high-throughput sequencing

Drought stress is one of the most important abiotic stresses. Cotton is classified as drought tolerant crop but the regulatory mechanism is unknown. MicroRNAs (miRNAs) have been implicated important roles in stress responses in many plants. However, the study of miRNAs in cotton responsive to drought stress is limited, especially in early-maturing cotton. In this study, we performed deep sequencing of small RNAs to identify known and novel miRNAs involved in the regulation of drought stress and understand the expression profile of miRNAs in early-maturing cotton. Three cotton small RNA libraries: non-stressed Shizao1 (early-maturing cotton variety) library (NSS), drought-stressed Shizao1 library (DSS) and non-stressed Jimian958 (medium-maturing cotton variety) library (NSJ) were constructed for deep sequencing. As a result, we identified a total of 64 known and 67 novel miRNAs in the 3 libraries and 88 of them were dramatically differentially expressed (greater than twofold) during drought stress. In addition, we found the expression of 41 miRNAs increased or reduced more than twofold in early-maturing cotton variety compared with that in medium-maturing cotton variety. Our results significantly increased the number of miRNAs in cotton and revealed for the first time the expression profile of miRNAs for early-maturing cotton.

Computational Identification of MicroRNAs and Their Transcript Target(s) in Field Mustard (Brassica rapa L.)

BACKGROUND

Micro RNAs (miRNAs) are a pivotal part of non-protein-coding endogenous small RNA molecules that regulate the genes involved in plant growth and development, and respond to biotic and abiotic environmental stresses posttranscriptionally.

OBJECTIVE

In the present study, we report the results of a systemic search for identification of new miRNAs in B. rapa using homology-based ESTs (Expressed Sequence Tags) analysis and considering a series of fi ltration criteria.

MATERIALS AND METHODS

Plant mature miRNA sequences were searched in non-protein coding ESTs registered in NCBI EST database. Zuker RNA folding algorithm was used to generate the secondary structures of the ESTs. Potential sequences were candidate as miRNA genes and characterized evolutionarily only and if only they fi t some described criteria. Also, the web tool psRNATarget was applied to predict candidate B. rapa miRNA targets.

RESULTS

In this study, 10 novel miRNAs from B. rapa belonging to 6 miRNA families were identified using EST-based homology analysis by considering a series of fi ltration criteria. All potent miRNAs appropriate fold back structure. Several potential targets with known/unknown functions for these novel miRNAs were identified. The target genes mainly encode transcription factors, enzymes, DNA binding proteins, disease resistance proteins, carrier proteins and other biological processes.

CONCLUSIONS

MicroRNA having diverse functions in plant species growth, development and evolution by posttranscriptionally regulating the levels of specific transcriptome so by effecting on their translation products. Research in miRNA led to the identification of many miRNAs and their regulating genes from diverse plant species.

Genome-wide identification and characterization of microRNAs differentially expressed in fibers in a cotton phytochrome A1 RNAi line

Cotton fiber is an important commodity throughout the world. Fiber property determines fiber quality and commercial values. Previous studies showed that silencing phytochrome A1 gene (PHYA1) by RNA interference in Upland cotton (Gossypium hirsutum L. cv. Coker 312) had generated PHYA1 RNAi lines with simultaneous improvements in fiber quality (longer, stronger and finer fiber) and other key agronomic traits. Characterization of the altered molecular processes in these RNAi genotypes and its wild-type controls is a great interest to better understand the PHYA1 RNAi phenotypes. In this study, a total of 77 conserved miRNAs belonging to 61 families were examined in a PHYA1 RNAi line and its parental Coker 312 genotype by using multiplex sequencing. Of these miRNAs, seven (miR7503, miR7514, miR399c, miR399d, miR160, miR169b, and miR2950) were found to be differentially expressed in PHYA1 RNAi cotton. The target genes of these differentially expressed miRNAs were involved in the metabolism and signaling pathways of phytohormones, which included Gibberellin, Auxin and Abscisic Acid. The expression of several MYB transcription factors was also affected by miRNAs in RNAi cotton. In addition, 35 novel miRNAs (novel miR1-novel miR35) were identified in fibers for the first time in this study. Target genes of vast majority of these novel miRNAs were also predicted. Of these, nine novel miRNAs (novel-miR1, 2, 16, 19, 26, 27, 28, 31 and 32) were targeted to cytochrome P450-like TATA box binding protein (TBP). The qRT-PCR confirmed expression levels of several differentially regulated miRNAs. Expression patterns of four miRNAs-targets pairs were also examined via RNA deep sequencing. Together, the results imply that the regulation of miRNA expression might confer to the phenotype of the PHYA1 RNAi line(s) with improved fiber quality.

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.

Identification and profiling of conserved and novel microRNAs involved in oil and oleic acid production during embryogenesis in Carya cathayensis Sarg

MicroRNAs (miRNAs) are important regulators of plant development and fruit formation. Mature embryos of hickory (Carya cathayensis Sarg.) nuts contain more than 70% oil (comprising 90% unsaturated fatty acids), along with a substantial amount of oleic acid. To understand the roles of miRNAs involved in oil and oleic acid production during hickory embryogenesis, three small RNA libraries from different stages of embryogenesis were constructed. Deep sequencing of these three libraries identified 95 conserved miRNAs with 19 miRNA*s, 7 novel miRNAs (as well as their corresponding miRNA*s), and 26 potentially novel miRNAs. The analysis identified 15 miRNAs involved in oil and oleic acid production that are differentially expressed during embryogenesis in hickory. Among them, nine miRNA sequences, including eight conserved and one novel, were confirmed by qRT-PCR. In addition, 145 target genes of the novel miRNAs were predicted using a bioinformatic approach. Our results provide a framework for better understanding the roles of miRNAs during embryogenesis in hickory.

Identification and expression analysis of cellular and viral microRNAs in CyHV3-infected KCF-1 cells

MicroRNAs (miRNAs) are small non-coding RNAs with approximately 22 nucleotides (nt) that are encoded by a diverse range of metazoan eukaryotes, plants and viruses. CyHV-3 (cyprinid herpesvirus-3) is a member of the Alloherpesviridae virus family and has caused severe economic losses for the common carp and koi carp fishery industries. In this study, a total of 15,987,652 clean reads were generated from a cDNA library of CyHV-3-infected KCF-1 (koi caudal fin) cells using high-throughput sequencing technology. Following annotation and secondary structure prediction, 28 miRNAs were identified as novel candidate miRNAs encoded by common carp (Cyprinus carpio), and seven miRNAs were shown to be encoded by CyHV-3. Next, 19 host miRNAs and seven viral miRNAs were validated by stem-loop real-time PCR. Northern blot analysis confirmed the presence of 14 host miRNAs and five CyHV-3-encoded novel miRNAs. The results of this study expand the knowledge of common carp and CyHV-3 microRNAs and provide a useful theoretical foundation for further study of CyHV-3.

Transcriptome and Degradome of microRNAs and Their Targets in Response to Drought Stress in the Plants of a Diploid and Its Autotetraploid Paulownia australis

MicroRNAs (miRNAs) are small, non-coding RNAs that play vital roles in plant growth, development, and stress response. Increasing numbers of studies aimed at discovering miRNAs and analyzing their functions in plants are being reported. In this study, we investigated the effect of drought stress on the expression of miRNAs and their targets in plants of a diploid and derived autotetraploid Paulownia australis. Four small RNA (sRNA) libraries and four degradome libraries were constructed from diploid and autotetraploid P. australis plants treated with either 75% or 25% relative soil water content. A total of 33 conserved and 104 novel miRNAs (processing precision value > 0.1) were identified, and 125 target genes were identified for 36 of the miRNAs by using the degradome sequencing. Among the identified miRNAs, 54 and 68 were differentially expressed in diploid and autotetraploid plants under drought stress (25% relative soil water content), respectively. The expressions of miRNAs and target genes were also validated by quantitative real-time PCR. The results showed that the relative expression trends of the randomly selected miRNAs were similar to the trends predicted by Illumina sequencing. And the correlations between miRNAs and their target genes were also analyzed. Furthermore, the functional analysis showed that most of these miRNAs and target genes were associated with plant development and environmental stress response. This study provided molecular evidence for the possible involvement of certain miRNAs in the drought response and/or tolerance in P. australis, and certain level of differential expression between diploid and autotetraploid plants.

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.

Integrative transcriptome, proteome, phosphoproteome and genetic mapping reveals new aspects in a fiberless mutant of cotton

To investigate the molecular mechanisms of fiber initiation in cotton (Gossypium spp.), an integrated approach combining transcriptome, iTRAQ-based proteome and genetic mapping was taken to compare the ovules of the Xuzhou 142 wild type (WT) with its fuzzless-lintless (fl) mutant at −3 and 0 day post-anthesis. A total of 1,953 mRNAs, 187 proteins, and 131 phosphoproteins were differentially expressed (DE) between WT and fl, and the levels of transcripts and their encoded proteins and phosphoproteins were highly congruent. A functional analysis suggested that the abundance of proteins were mainly involved in amino sugar, nucleotide sugar and fatty acid metabolism, one carbon pool for folate metabolism and flavonoid biosynthesis. qRT-PCR, Western blotting, and enzymatic assays were performed to confirm the regulation of these transcripts and proteins. A molecular mapping located the lintless gene li3 in the fl mutant on chromosome 26 for the first time. A further in-silico physical mapping of DE genes with sequence variations between fl and WT identified one and four candidate genes in the li3 and n2 regions, respectively. Taken together, the transcript abundance, phosphorylation status of proteins at the fiber initiation stage and candidate genes have provided insights into regulatory processes underlying cotton fiber initiation.

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.

A cotton miRNA is involved in regulation of plant response to salt stress

The present study functionally identified a new microRNA (microRNA ovual line 5, miRNVL5) with its target gene GhCHR from cotton (Gossypium hirsutum). The sequence of miRNVL5 precursor is 104 nt long, with a well developed secondary structure. GhCHR contains two DC1 and three PHD Cys/His-rich domains, suggesting that GhCHR encodes a zinc-finger domain-containing transcription factor. miRNVL5 and GhCHR express at various developmental stages of cotton. Under salt stress (50-400 mM NaCl), miRNVL5 expression was repressed, with concomitant high expression of GhCHR in cotton seedlings. Ectopic expression of GhCHR in Arabidopsis conferred salt stress tolerance by reducing Na(+) accumulation in plants and improving primary root growth and biomass. Interestingly, Arabidopsis constitutively expressing miRNVL5 showed hypersensitivity to salt stress. A GhCHR orthorlous gene At2g44380 from Arabidopsis that can be cleaved by miRNVL5 was identified by degradome sequencing, but no confidential miRNVL5 homologs in Arabidopsis have been identified. Microarray analysis of miRNVL5 transgenic Arabidopsis showed six downstream genes (CBF1, CBF2, CBF3, ERF4, AT3G22920, and AT3G49200), which were induced by salt stress in wild-type but repressed in miRNVL5-expressing Arabidopsis. These results indicate that miRNVL5 is involved in regulation of plant response to salt stress.

Unveiling the Micronome of Cassava (Manihot esculenta Crantz)

MicroRNAs (miRNAs) are an important class of endogenous non-coding single-stranded small RNAs (21-24 nt in length), which serve as post-transcriptional negative regulators of gene expression in plants. Despite the economic importance of Manihot esculenta Crantz (cassava) only 153 putative cassava miRNAs (from multiple germplasm) are available to date in miRBase (Version 21), and identification of a number of miRNAs from the cassava EST database have been limited to comparisons with Arabidopsis. In this study, mature sequences of all known plant miRNAs were used as a query for homologous searches against cassava EST and GSS databases, and additional identification of novel and conserved miRNAs were gleaned from next generation sequencing (NGS) of two cassava landraces (T200 from southern Africa and TME3 from West Africa) at three different stages post explant transplantation and acclimatization. EST and GSS derived data revealed 259 and 32 miRNAs in cassava, and one of the miRNA families (miR2118) from previous studies has not been reported in cassava. NGS data collectively displayed expression of 289 conserved miRNAs in leaf tissue, of which 230 had not been reported previously. Of the 289 conserved miRNAs identified in T200 and TME3, 208 were isomiRs. Thirty-nine novel cassava-specific miRNAs of low abundance, belonging to 29 families, were identified. Thirty-eight (98.6%) of the putative new miRNAs identified by NGS have not been previously reported in cassava. Several miRNA targets were identified in T200 and TME3, highlighting differential temporal miRNA expression between the two cassava landraces. This study contributes to the expanding knowledge base of the micronome of this important crop.

MicroRNA Profiles in Spontaneous Decidualized Menstrual Endometrium and Early Pregnancy Decidua with Successfully Implanted Embryos

To comparatively analyze the human microRNA (miRNA) profiles between spontaneous decidualized menstrual endometrium and early pregnancy decidua by an in-depth sequencing of miRNAs. The specific miRNAs expressed at conception might be involved in pregnancy establishment and expression of let-7f-5p and let-7g-5p was experimentally up-regulated or inhibited to assess the effect on the expression of IGF2BP-1 and IGF2R in vitro, respectively. Samples of endometria and deciduas were obtained from 25 women who suffered from tubal or male factor subfertility and from 35 early pregnant women who underwent pregnancy termination at 6–8 weeks gestation were irrespectively collected and comparatively analyzed by miRNA sequencing and differential expression of known and novel miRNAs was analyzed using bioinformatics. The 2042 miRNA expression was analyzed in the study and the differential expression of six miRNAs was validated by qRT-PCR. The expression of four miRNAs in decidua samples was down-regulated (miR-34c, miR-92a, miR-181a-5p, and miR-191), whereas the expression of miR-10a-5p and let-7f-5p was significantly up-regulated. The expression of IGF2BP-1 and IGF2R declined and increased with overexpression and inhibition of let-7f-5p and let-7g-5p, respectively. Changes in the expression of particular miRNAs might play a role in the physiology of decidualization following successful embryo implantation, ultimately resulting in continuous decidualization.

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.

Identification of microRNAs differentially expressed involved in male flower development

Hickory (Carya cathayensis Sarg.) is one of the most economically important woody trees in eastern China, but its long flowering phase delays yield. Our understanding of the regulatory roles of microRNAs (miRNAs) in male flower development in hickory remains poor. Using high-throughput sequencing technology, we have pyrosequenced two small RNA libraries from two male flower differentiation stages in hickory. Analysis of the sequencing data identified 114 conserved miRNAs that belonged to 23 miRNA families, five novel miRNAs including their corresponding miRNA*s, and 22 plausible miRNA candidates. Differential expression analysis revealed 12 miRNA sequences that were upregulated in the later (reproductive) stage of male flower development. Quantitative real-time PCR showed similar expression trends as that of the deep sequencing. Novel miRNAs and plausible miRNA candidates were predicted using bioinformatic analysis methods. The miRNAs newly identified in this study have increased the number of known miRNAs in hickory, and the identification of differentially expressed miRNAs will provide new avenues for studies into miRNAs involved in the process of male flower development in hickory and other related trees.