Identification of cotton microRNAs and their targets

Identification and Characterization of Evolutionary Conserved Muskmelon Non-coding miRNAs and Their Target Proteins

MicroRNAs (miRNAs) are short, endogenously encoded small RNAs, 18-26 nucleotides in length, which can posttranscriptionally regulate gene expression through translation inhibition or endonucleolytic cleavage. The muskmelon is one of the most widely cultivated fruits in the Cucurbitaceae family. Despite its significance, only 120 miRNAs from different families have been reported in muskmelon. In this study, we aimed to expand this knowledge base by predicting 40 new miRNAs in muskmelon using a spectrum of genomic-based tools. Precursor and mature sequences were obtained from microRNA registry database as reference and analyzed via the basic local alignment search tool (Blastn) for ESTs identification. After removing the non-coding sequences, the remaining candidate sequences were analyzed using MFOLD to generate secondary structures for the newly predicted miRNAs. Additionally, the predicted muskmelon miRNAs were validated using a set of five randomly chosen primers and RT-PCR. Through gene ontology (GO) analysis, we identified 89 targets associated with newly predicted muskmelon miRNAs. Transcription factor-coding genes play a crucial role in plant growth and development. Additionally, the miR4249 has been found to have the same targets in muskmelon that have been linked to cell signaling and transcription factors. The identified targets are integral for diverse biological processes including plant growth, development, metabolism, aging, disease resistance, and resistance to environmental stresses, such as salt, cold, and oxidative stress. As a result, the outcomes of this study demonstrate that this mechanism not only contributes to the production of a higher quality crop but also enhances overall production.

Genome-wide association study for boll weight in Gossypium hirsutum races

High yield has always been an essential target in almost all of the cotton breeding programs. Boll weight (BW) is a key component of cotton yield. Numerous linkage mapping and genome-wide association studies (GWAS) have been performed to understand the genetic mechanism of BW, but information on the markers/genes controlling BW remains limited. In this study, we conducted a GWAS for BW using 51,268 high-quality single-nucleotide polymorphisms (SNPs) and 189 Gossypium hirsutum accessions across five different environments. A total of 55 SNPs significantly associated with BW were detected, of which 29 and 26 were distributed in the A and D subgenomes, respectively. Five SNPs were simultaneously detected in two environments. For TM5655, TM8662, TM36371, and TM50258, the BW grouped by alleles of each SNP was significantly different. The ± 550 kb regions around these four key SNPs contained 262 genes. Of them, Gh_A02G1473, Gh_A10G1765, and Gh_A02G1442 were expressed highly at 0 to 1 days post-anthesis (dpa), - 3 to 0 dpa, and - 3 to 0 dpa in ovule of TM-1, respectively. They were presumed as the candidate genes for fiber cell differentiation, initiation, or elongation based on gene annotation of their homologs. Overall, these results supplemented valuable information for dissecting the genetic architecture of BW and might help to improve cotton yield through molecular marker-assisted selection breeding and molecular design breeding.

Advances in understanding epigenetic regulation of plant trichome development: a comprehensive review

Plant growth and development are controlled by a complex gene regulatory network, which is currently a focal point of research. It has been established that epigenetic factors play a crucial role in plant growth. Trichomes, specialized appendages that arise from epidermal cells, are of great significance in plant growth and development. As a model system for studying plant development, trichomes possess both commercial and research value. Epigenetic regulation has only recently been implicated in the development of trichomes in a limited number of studies, and microRNA-mediated post-transcriptional regulation appears to dominate in this context. In light of this, we have conducted a review that explores the interplay between epigenetic regulations and the formation of plant trichomes, building upon existing knowledge of hormones and transcription factors in trichome development. Through this review, we aim to deepen our understanding of the regulatory mechanisms underlying trichome formation and shed light on future avenues of research in the field of epigenetics as it pertains to epidermal hair growth.

Integrated miRNA, target mRNA, and metabolome profiling of Tinospora cordifolia with reference to berberine biosynthesis

MicroRNAs play a central role in gene regulation and emerge as novel targets for secondary metabolites improvement in plants. The crops thus can be improved through knowledge obtained by the study of miRNAs because of their conserved nature in gene regulation. The present study has been carried out on Tinospora cordifolia (T. cordifolia), because of its illimitable application for the treatment of various diseases. This plant has tremendous medicinal properties, yet unexplored at the molecular level, and has not received much recognition in the scientific field. Thus, here computational analysis was performed to identify T. cordifolia miRNAs using EST database. Using these miRNAs, we predicted their targets which were found to be associated with the regulation of diverse gene networks including 433 berberine biosynthesis genes in T. cordifolia. Further, selected miRNAs were validated and their expression was detected in different T. cordifolia tissues followed by expression analysis of their target mRNAs. These data were then compared with the metabolic profile of T. cordifolia with an emphasis on therapeutically important compound berberine. In this study, we did simultaneous miRNA/target gene expression and metabolome analysis which opens a new way for initiating new proposition and prioritization of miRNAs/genes/metabolites for targeted follow‑up metabolic engineering experimentations.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03342-9.

An overview of salinity stress, mechanism of salinity tolerance and strategies for its management in cotton

Salinity stress is one of the primary threats to agricultural crops resulting in impaired crop growth and development. Although cotton is considered as reasonably salt tolerant, it is sensitive to salt stress at some critical stages like germination, flowering, boll formation, resulting in reduced biomass and fiber production. The mechanism of partial ion exclusion (exclusion of Na+ and/or Cl-) in cotton appears to be responsible for the pattern of uptake and accumulation of harmful ions (Na+ and Cl) in tissues of plants exposed to saline conditions. Maintaining high tissue K+/Na+ and Ca2+/Na+ ratios has been proposed as a key selection factor for salt tolerance in cotton. The key adaptation mechanism in cotton under salt stress is excessive sodium exclusion or compartmentation. Among the cultivated species of cotton, Egyptian cotton (Gossypium barbadense L.) exhibit better salt tolerance with good fiber quality traits as compared to most cultivated cotton and it can be used to improve five quality traits and transfer salt tolerance into Upland or American cotton (Gossypium hirsutum L.) by interspecific introgression. Cotton genetic studies on salt tolerance revealed that the majority of growth, yield, and fiber traits are genetically determined, and controlled by quantitative trait loci (QTLs). Molecular markers linked to genes or QTLs affecting key traits have been identified, and they could be utilized as an indirect selection criterion to enhance breeding efficiency through marker-assisted selection (MAS). Transfer of genes for compatible solute, which are an important aspect of ion compartmentation, into salt-sensitive species is, theoretically, a simple strategy to improve tolerance. The expression of particular stress-related genes is involved in plant adaptation to environmental stressors. As a result, enhancing tolerance to salt stress can be achieved by marker assisted selection added with modern gene editing tools can boost the breeding strategies that defend and uphold the structure and function of cellular components. The intent of this review was to recapitulate the advancements in salt screening methods, tolerant germplasm sources and their inheritance, biochemical, morpho-physiological, and molecular characteristics, transgenic approaches, and QTLs for salt tolerance in cotton.

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.

Analysis of miRNAs responsive to long-term calcium deficiency in tef (Eragrostis tef (Zucc.) Trotter)

MicroRNAs (miRNAs) play an important role in growth, development, stress resilience, and epigenetic modifications of plants. However, the effect of calcium (Ca²⁺) deficiency on miRNA expression in the orphan crop tef (Eragrostis tef) remains unknown. In this study, we analyzed expression of miRNAs in roots and shoots of tef in response to Ca²⁺ treatment. miRNA-seq followed by bioinformatic analysis allowed us to identify a large number of small RNAs (sRNAs) ranging from 17 to 35 nt in length. A total of 1380 miRNAs were identified in tef experiencing long-term Ca²⁺ deficiency while 1495 miRNAs were detected in control plants. Among the miRNAs identified in this study, 161 miRNAs were similar with those previously characterized in other plant species and 348 miRNAs were novel, while the remaining miRNAs were uncharacterized. Putative target genes and their functions were predicted for all the known and novel miRNAs that we identified. Based on gene ontology (GO) analysis, the predicted target genes are known to have various biological and molecular functions including calcium uptake and transport. Pairwise comparison of differentially expressed miRNAs revealed that some miRNAs were specifically enriched in roots or shoots of low Ca²⁺-treated plants. Further characterization of the miRNAs and their targets identified in this study may help in understanding Ca²⁺ deficiency responses in tef and related orphan crops.

MiRNA fine tuning for crop improvement: using advance computational models and biotechnological tools

MiRNAs modulate target genes expression at post-transcriptional levels, by reducing spatial abundance of mRNAs. MiRNAs regulats plant metabolism, and emerged as regulators of plant stress responses. Which make miRNAs promising candidates for fine tuning to affectively alter crop stress tolerance and other important traits. With recent advancements in the computational biology and biotechnology miRNAs structure and target prediction is possible resulting in pin point editing; miRNA modulation can be done by up or down regulating miRNAs using recently available biotechnological tools (CRISPR Cas9, TALENS and RNAi). In this review we have focused on miRNA biogenesis, miRNA roles in plant development, plant stress responses and roles in signaling pathways. Additionally we have discussed latest computational prediction models for miRNA to target gene interaction and biotechnological systems used recently for miRNA modulation. We have also highlighted setbacks and limitations in the way of miRNA modulation; providing entirely a new direction for improvement in plant genomics primarily focusing miRNAs.

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.

Genome-Wide Analysis and Expression Profile of Superoxide Dismutase (SOD) Gene Family in Rapeseed (Brassica napus L.) under Different Hormones and Abiotic Stress Conditions

Superoxide dismutase (SOD) is an important enzyme that acts as the first line of protection in the plant antioxidant defense system, involved in eliminating reactive oxygen species (ROS) under harsh environmental conditions. Nevertheless, the SOD gene family was yet to be reported in rapeseed (Brassica napus L.). Thus, a genome-wide investigation was carried out to identify the rapeseed SOD genes. The present study recognized 31 BnSOD genes in the rapeseed genome, including 14 BnCSDs, 11 BnFSDs, and six BnMSDs. Phylogenetic analysis revealed that SOD genes from rapeseed and other closely related plant species were clustered into three groups based on the binding domain with high bootstrap values. The systemic analysis exposed that BnSODs experienced segmental duplications. Gene structure and motif analysis specified that most of the BnSOD genes displayed a relatively well-maintained exon-intron and motif configuration within the same group. Moreover, we identified five hormones and four stress- and several light-responsive cis-elements in the promoters of BnSODs. Thirty putative bna-miRNAs from seven families were also predicted, targeting 13 BnSODs. Gene ontology annotation outcomes confirm the BnSODs role under different stress stimuli, cellular oxidant detoxification processes, metal ion binding activities, SOD activity, and different cellular components. Twelve BnSOD genes exhibited higher expression profiles in numerous developmental tissues, i.e., root, leaf, stem, and silique. The qRT-PCR based expression profiling showed that eight genes (BnCSD1, BnCSD3, BnCSD14, BnFSD4, BnFSD5, BnFSD6, BnMSD2, and BnMSD10) were significantly up-regulated under different hormones (ABA, GA, IAA, and KT) and abiotic stress (salinity, cold, waterlogging, and drought) treatments. The predicted 3D structures discovered comparable conserved BnSOD protein structures. In short, our findings deliver a foundation for additional functional investigations on the BnSOD genes in rapeseed breeding programs.

In Silico Identification and Functional Characterization of Conserved miRNAs in the Genome of Cryptosporidium parvum

Cryptosporidium parvum, a predominant causal agent of a fatal zoonotic protozoan diarrhoeal disease called cryptosporidiosis, bears a worldwide public health concern for childhood mortality and poses a key threat to the dairy and water industries. MicroRNAs (miRNAs), small but powerful posttranscriptional gene silencing RNA molecules, regulate a variety of molecular, biological, and cellular processes in animals and plants. As to the present date, there is a paucity of information regarding miRNAs of C. parvum; hence, this study was used to identify miRNAs in the organism using a comprehensible expressed sequence tag-based homology search approach consisting of a series of computational screening process from the identification of putative miRNA candidates to the functional annotation of the important gene targets in C. parvum. The results revealed a conserved miRNA that targeted 487 genes in the model organism (Drosophila melanogaster) and 85 genes in C. parvum, of which 11 genes had direct involvements in several crucial virulence factors such as environmental oocyst protection, excystation, locomotion, adhesion, invasion, stress protection, intracellular growth, and survival. Besides, 20 genes showed their association with various major pathways dedicated for the ribosomal biosynthesis, DNA repair, transportation, protein production, gene expression, cell cycle, cell proliferation, development, immune response, differentiation, and nutrient metabolism of the organism in the host. Thus, this study provides a strong evidence of great impact of identified miRNA on the biology, virulence, and pathogenesis of C. parvum. Furthermore, the study suggests that the detected miRNA could be a potential epigenomic tool for controlling the protozoon through silencing those virulent and pathway-related target genes.

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.

Utilization of Transcriptome, Small RNA, and Degradome Sequencing to Provide Insights Into Drought Stress and Rewatering Treatment in Medicago ruthenica

Drought is a major limiting factor in foraging grass yield and quality. Medicago ruthenica (M. ruthenica) is a high-quality forage legume with drought resistance, cold tolerance, and strong adaptability. In this study, we integrated transcriptome, small RNA, and degradome sequencing in identifying drought response genes, microRNAs (miRNAs), and key miRNA-target pairs in M. ruthenica under drought and rewatering treatment conditions. A total of 3,905 genes and 50 miRNAs (45 conserved and 5 novel miRNAs) were significantly differentially expressed in three test conditions (CK: control, DS: plants under drought stress, and RW: plants rewatering after drought stress). The degradome sequencing (AllenScore < 4) analysis revealed that 104 miRNAs (11 novel and 93 conserved miRNAs) were identified with 263 target transcripts, forming 296 miRNA-target pairs in three libraries. There were 38 differentially expressed targets from 16 miRNAs in DS vs. CK, 31 from 11 miRNAs in DS vs. RW, and 6 from 3 miRNAs in RW vs. CK; 21, 18, and 3 miRNA-target gene pairs showed reverse expression patterns in DS vs. CK, DS vs. RW, and RW vs. CK comparison groups, respectively. These findings provide valuable information for further functional characterization of genes and miRNAs in response to abiotic stress, in general, and drought stress in M. ruthenica, and potentially contribute to drought resistance breeding of forage in the future.

Non-coding RNAs having strong positive interaction with mRNAs reveal their regulatory nature during flowering in a wild relative of pigeonpea (Cajanus scarabaeoides)

In higher plants, flower development is a result of crosstalk between many factors like photoperiod, vernalization, hormone concentration, epigenetic modification etc. and is also regulated by non-coding RNAs (ncRNAs). In the present study, we are reporting the involvement of long non-coding RNAs (lncRNAs) and miRNAs during the process of flower development in Cajanus scarabaeoides, an important wild relative of pigeonpea. The transcriptome of floral and leaf tissues revealed a total of 1672 lncRNAs and 57 miRNAs being expressed during flower development. Prediction analysis of identified lncRNAs showed that 1593 lncRNAs were targeting 3420 mRNAs and among these, 98 were transcription factors (TFs) belonging to 48 groups. All the identified 57 miRNAs were novel, suggesting their genera specificity. Prediction of the secondary structure of lncRNAs and miRNAs followed by interaction analysis revealed that 199 lncRNAs could interact with 47 miRNAs where miRNAs were acting in the root of interaction. Gene Ontology of the ncRNAs and their targets showed the potential role of lncRNAs and miRNAs in the flower development of C. scarabaeoides. Among the identified interactions, 17 lncRNAs were endogenous target mimics (eTMs) for miRNAs that target flowering-related transcription factors. Expression analysis of identified transcripts revealed that higher expression of Csa-lncRNA_1231 in the bud sequesters Csa-miRNA-156b by indirectly mimicking the miRNA and leading to increased expression of flower-specific SQUAMOSA promoter-binding protein-like (SPL-12) TF indicating their potential role in flower development. The present study will help in understanding the molecular regulatory mechanism governing the induction of flowering in C. scarabaeoides.

Machine learning improves our knowledge about miRNA functions towards plant abiotic stresses

During the last two decades, human has increased his knowledge about the role of miRNAs and their target genes in plant stress response. Biotic and abiotic stresses result in simultaneous tissue-specific up/down-regulation of several miRNAs. In this study, for the first time, feature selection algorithms have been used to investigate the contribution of individual plant miRNAs in Arabidopsis thaliana response towards different levels of several abiotic stresses including drought, salinity, cold, and heat. Results of information theory-based feature selection revealed that miRNA-169, miRNA-159, miRNA-396, and miRNA-393 had the highest contributions to plant response towards drought, salinity, cold, and heat, respectively. Furthermore, regression models, i.e., decision tree (DT), support vector machines (SVMs), and Naïve Bayes (NB) were used to predict the plant stress by having the plant miRNAs' concentration. SVM with Gaussian kernel was capable of predicting plant stress (R2 = 0.96) considering miRNA concentrations as input features. Findings of this study prove the performance of machine learning as a promising tool to investigate some aspects of miRNAs' contribution to plant stress responses that have been undiscovered until today.

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.

A Comparative Whole Genome Sequence Analysis Leads to Identification of Repeat-Associated Evolutionarily Conserved miRNAs in Bombyx mori (Lepidoptera: Bombycidae)

MicroRNAs (miRNAs) are evolutionarily conserved small noncoding RNAs, which play important regulatory roles in various biological processes. In this study, we have developed a computational approach for detecting conserved miRNAs based on comparison of whole genome sequences of closely related species by considering various key features of experimentally validated miRNAs. By applying this approach, we have identified 34 new miRNAs from Bombyx mori (L.), which are also conserved in Drosophila melanogaster (Diptera: Drosophilidae) and Anopheles gambiae (Diptera: Culicidae). Most of these miRNAs were associated with repeat region of the genome. We did an expression analysis of the 34 newly predicted miRNAs and found that 30 of these miRNAs are expressing in different tissues of B. mori. Besides, we have also predicted the putative targets of these miRNAs in B. mori based on several known characteristic features of miRNA::mRNA duplexes and found that these targets include diverse range of functions, suggesting multiple layers of gene regulation of various important biological processes.

Vascular Calcification and not Arrhythmia in Idiopathic Atrial Fibrillation Associates with Sex Differences in Diabetic Microvascular Injury miRNA Profiles

BACKGROUND

Atrial Fibrillation (AF) in patients without concomitant cardiovascular pathophysiological disease, is called idiopathic Atrial Fibrillation (iAF). Nonetheless, iAF patients have often times subclinical coronary (micro) vascular dysfunction and, particularly in women, a higher prevalence of subsequent cardiovascular comorbidities. Previously, we identified a plasma miRNA association with diabetes and microvascular injury in Diabetic Nephropathy (DN) patients. Therefore, in this study we assessed whether plasma levels of these diabetic, microvascular injury associated miRNAs reflect microvascular integrity in iAF patients, associated with the presence of paroxysmal arrhythmia or instead are determined by concealed coronary artery disease.

METHODS

Circulating levels of a pre-selected set of diabetic, (micro) vascular injury associated miRNAs, were measured in 59 iAF patients compared to 176 Sinus Rhythm (SR) controls. Furthermore, the presence of coronary artery and aortic calcification in each patient was assessed using Cardiac Computed Tomography Angiography (CCTA).

RESULTS

Paroxysmal arrhythmia in iAF patients did not result in significant miRNA expression profile differences in iAF patients compared to SR controls. Nonetheless, coronary artery calcification (CAC) was associated with higher levels of miRNAs-103, -125a-5p, -221 and -223 in men. In women, CAC was associated with higher plasma levels of miRNA-27a and miRNA-126 and correlated with Agatston scores. Within the total population, ascending Aortic Calcification (AsAC) patients displayed increased plasma levels of miRNA-221, while women, in particular, demonstrated a Descending Aorta Calcification (DAC) associated increase in miRNA-212 levels.

CONCLUSIONS

Diabetic microvascular injury associated miRNAs in iAF are associated with subclinical coronary artery disease in a sex-specific way and confirm the notion that biological sex identifies iAF subgroups that may require dedicated clinical care.

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.

Identification and expression analysis of a microRNA cluster derived from pre-ribosomal RNA in Papaver somniferum L. and Papaver bracteatum L

Opium poppy (Papaver somniferum L.) is one of the ancient medical crops, which produces several important alkaloids such as morphine, noscapine, sanguinarine and codeine. MicroRNAs are endogenous non-coding RNAs that play important regulatory roles in plant diverse biological processes. Many plant miRNAs are encoded as single transcriptional units, in contrast to animal miRNAs, which are often clustered. Herein, using computational approaches, a total of 22 miRNA precursors were identified, which five of them were located as a clustered in pre-ribosomal RNA. Afterward, the transcript level of the precursor and the mature of clustered miRNAs in two species of the Papaveraceae family, i.e. P. somniferum L. and P. bracteatum L, were quantified by RT-PCR. With respect to obtained results, these clustered miRNAs were expressed differentially in different tissues of these species. Moreover, using target prediction and Gene Ontology (GO)-based on functional classification indicated that these miRNAs might play crucial roles in various biological processes as well as metabolic pathways. In this study, we discovered the clustered miRNA derived from pre-rRNA, which may shed some light on the importance of miRNAs in the plant kingdom.

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

Cytoplasmic male sterility (CMS) lines provide crucial material to harness heterosis for crop plants, which serves as an important strategy for hybrid seed production. However, the molecular mechanism remains obscure. Although microRNAs (miRNAs) play important roles in vegetative growth and reproductive growth, there are few reports on miRNAs regulating the development of male sterility in Upland cotton. In present study, 12 small RNA libraries were constructed and sequenced for two development stages of flower buds from a CMS line and its maintainer line. Based on the results, 256 novel miRNAs were allocated to 141 new miRNA families, and 77 known miRNAs belonging to 54 conserved miRNA families were identified as well. Comparative analysis revealed that 61 novel and 10 conserved miRNAs were differentially expressed. Further transcriptome analysis identified 232 target genes for these miRNAs, which participated in cellular developmental process, cell death, pollen germination, and sexual reproduction. In addition, expression patterns of typical miRNA and the negatively regulated target genes, such as PPR, ARF, AP2, and AFB, were verified by qRT-PCR in cotton flower buds. These targets were previously reported to be related to reproduction development and male sterility, suggesting that miRNAs might act as regulators of CMS occurrence. Some miRNAs displayed specific expression profiles in special developmental stages of CMS line and its fertile hybrid (F₁). Present study offers new information on miRNAs and their related target genes in exploiting CMS mechanism, and revealing the miRNA regulatory networks in Upland cotton.

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.

Biocomputational identification and validation of novel microRNAs predicted from bubaline whole genome shotgun sequences

MicroRNAs (miRNAs) are small (19-25 base long), non-coding RNAs that regulate post-transcriptional gene expression by cleaving targeted mRNAs in several eukaryotes. The miRNAs play vital roles in multiple biological and metabolic processes, including developmental timing, signal transduction, cell maintenance and differentiation, diseases and cancers. Experimental identification of microRNAs is expensive and lab-intensive. Alternatively, computational approaches for predicting putative miRNAs from genomic or exomic sequences rely on features of miRNAs viz. secondary structures, sequence conservation, minimum free energy index (MFEI) etc. To date, not a single miRNA has been identified in bubaline (Bubalus bubalis), which is an economically important livestock. The present study aims at predicting the putative miRNAs of buffalo using comparative computational approach from buffalo whole genome shotgun sequencing data (INSDC: AWWX00000000.1). The sequences were blasted against the known mammalian miRNA. The obtained miRNAs were then passed through a series of filtration criteria to obtain the set of predicted (putative and novel) bubaline miRNA. Eight miRNAs were selected based on lowest E-value and validated by real time PCR (SYBR green chemistry) using RNU6 as endogenous control. The results from different trails of real time PCR shows that out of selected 8 miRNAs, only 2 (hsa-miR-1277-5p; bta-miR-2285b) are not expressed in bubaline PBMCs. The potential target genes based on their sequence complementarities were then predicted using miRanda. This work is the first report on prediction of bubaline miRNA from whole genome sequencing data followed by experimental validation. The finding could pave the way to future studies in economically important traits in buffalo.

In-silico based identification and functional analyses of miRNAs and their targets in Cowpea (Vigna unguiculata L.)

Cowpea (Vigna unguiculata L.) is an important leguminous plant and a good diet due to presence of carbohydrate and high protein contents. Currently, only few cowpea microRNAs (miRNAs) are reported. This study is intended to identify and functionally analyze new miRNAs and their targets in cowpea. An in-silico based homology search approach was applied and a total of 46 new miRNAs belonging to 45 families were identified and functionally annotated from the cowpea expressed sequence tags (ESTs). All these potential miRNAs are reported here for the first time in cowpea. The 46 new miRNAs were also observed with stable hairpin structures with minimum free energy, ranging from -10 to -132 kcal mol^-1 with an average of -40 kcal mol^-1. The length of new cowpea miRNAs are ranged from 18 to 26 nt with an average of 21 nt. The cowpea miRNA-vun-mir4414, is found as pre-miRNA cluster for the first time in cowpea. Furthermore, a set of 138 protein targets were also identified for these newly identified 46 cowpea miRNAs. These targets have significant role in various biological processes, like metabolism, transcription regulation as transcription factor, cell transport, signal transduction, growth & development and structural proteins. These findings are the significant basis to utilize and manage this important leguminous plant-cowpea for better nutritional properties and tolerance for biotic and abiotic stresses.

Genome-wide characterization and expression analyses of superoxide dismutase (SOD) genes in Gossypium hirsutum

BACKGROUND

Superoxide dismutases (SODs) are a key antioxidant enzyme family, which have been implicated in protecting plants against the toxic effects of reactive oxygen species. Despite current studies have shown that the gene family are involved in plant growth and developmental processes and biotic and abiotic stress responses, little is known about its functional role in upland cotton.

RESULTS

In the present study, we comprehensively analyzed the characteristics of the SOD gene family in upland cotton (Gossypium hirsutum). Based on their conserved motifs, 18 GhSOD genes were identified and phylogenetically classified into five subgroups which corroborated their classifications based on gene-structure patterns and subcellular localizations. The GhSOD sequences were distributed at different densities across 12 of the 26 chromosomes. The conserved domains, gene family evolution cis-acting elements of promoter regions and miRNA-mediated posttranscriptional regulation were predicted and analyzed. In addition, the expression pattern of 18 GhSOD genes were tested in different tissues/organs and developmental stages, and different abiotic stresses and abscisic acid, which indicated that the SOD gene family possessed temporal and spatial specificity expression specificity and may play important roles in reactive oxygen species scavenging caused by various stresses in upland cotton.

CONCLUSIONS

This study describes the first genome-wide analysis of the upland cotton SOD gene family, and the results will help establish a foundation for the further cloning and functional verification of the GhSOD gene family during stress responses, leading to crop improvement.

microRNAs as reference genes for quantitative PCR in cotton

Cotton (Gossypium hirsutum) is the most important non-food plant in the world. Studies concerning the fiber quality and plant fitness of cotton at molecular level depend on high sensitive and reproducible gene-expression assays. However, only a few reports have described genes suitable for normalizing gene expression data. In this study, we report for the first time that microRNAs (miRNAs) are reliable reference genes (RGs) for cotton gene expression normalization in quantitative real-time reverse transcription (RT)-PCR. The stability of cotton miRNAs was assayed in root, stem, leaf and flower samples from three different cultivars [FiberMax (FM966), Delta Opal (DO) and Cedro] and under conditions of biotic stress caused by infection with Cotton leafroll dwarf virus (CLRDV). The stability of mRNAs already described as reference genes in cotton was also assessed. The geNorm, NormFinder, BestKeeper and ΔCt algorithms were used to select the best reference genes. In 8 of the 12 sets tested, miRNAs (miR172, 168 and 390) were found to be the best RGs. To validate the best selected RGs, miR159, miR164, miR2118, miR2910, miR3476, GhDCL2 and GhDCL4 expression levels were evaluated under biotic stress conditions, and miR164 and a putative myo-inositol oxigenase gene (GhMIOX) were assessed in leaves and flowers. The RGs selected in this work proved to be excellent reference genes in the two cases studied. Our results support the use of miRNAs as reference genes for miRNA and protein-coding genes.

Functional Roles of microRNAs in Agronomically Important Plants-Potential as Targets for Crop Improvement and Protection

MicroRNAs (miRNAs) are a class of small non-coding RNAs that have recently emerged as important regulators of gene expression, mainly through cleavage and/or translation inhibition of the target mRNAs during or after transcription. miRNAs play important roles by regulating a multitude of biological processes in plants which include maintenance of genome integrity, development, metabolism, and adaptive responses toward environmental stresses. The increasing population of the world and their food demands requires focused efforts for the improvement of crop plants to ensure sustainable food production. Manipulation of mRNA transcript abundance via miRNA control provides a unique strategy for modulating differential plant gene expression and miRNAs are thus emerging as the next generation targets for genetic engineering for improvement of the agronomic properties of crops. However, a deeper understanding of its potential and the mechanisms involved will facilitate the design of suitable strategies to obtain the desirable traits with minimum trade-offs in the modified crops. In this regard, this review highlights the diverse roles of conserved and newly identified miRNAs in various food and industrial crops and recent advances made in the uses of miRNAs to improve plants of agronomically importance so as to significantly enhance crop yields and increase tolerance to various environmental stress agents of biotic-or abiotic origin.

Mining, identification and function analysis of microRNAs and target genes in peanut (Arachis hypogaea L.)

In the present investigation, a total of 60 conserved peanut (Arachis hypogaea L.) microRNA (miRNA) sequences, belonging to 16 families, were identified using bioinformatics methods. There were 392 target gene sequences, identified from 58 miRNAs with Target-align software and BLASTx analyses. Gene Ontology (GO) functional analysis suggested that these target genes were involved in mediating peanut growth and development, signal transduction and stress resistance. There were 55 miRNA sequences, verified employing a poly (A) tailing test, with a success rate of up to 91.67%. Twenty peanut target gene sequences were randomly selected, and the 5' rapid amplification of the cDNA ends (5'-RACE) method were used to validate the cleavage sites of these target genes. Of these, 14 (70%) peanut miRNA targets were verified by means of gel electrophoresis, cloning and sequencing. Furthermore, functional analysis and homologous sequence retrieval were conducted for target gene sequences, and 26 target genes were chosen as the objects for stress resistance experimental study. Real-time fluorescence quantitative PCR (qRT-PCR) technology was applied to measure the expression level of resistance-associated miRNAs and their target genes in peanut exposed to Aspergillus flavus (A. flavus) infection and drought stress, respectively. In consequence, 5 groups of miRNAs & targets were found accorded with the mode of miRNA negatively controlling the expression of target genes. This study, preliminarily determined the biological functions of some resistance-associated miRNAs and their target genes in peanut.

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.

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.

Sequence data mining in search of hookworm (Necator americanus) microRNAs

The new world hookworm, Necator americanus is a soil-transmitted nematode responsible for Necatoriasis (a type of helminthiasis) in hosts such as humans, dogs, and cats. N. americanus genome and transcriptome has been sequenced and a draft assembly analysis has been published highlighting protein coding genes and possible drug target proteins. Hookworm microRNA identification, annotations and their public release is yet to be attempted. The same is evident from lack of hookworm miRNA information in related popular public nucleotide sequence repositories such as miRBase, GenBank, WormBase etc. Therefore, in the present study we addressed these issues using EST and assembled transcript sequence information of hookworm. Using computational approaches, we identified three miRNAs precursor sequences and their mature forms. We also identified their potential targets from hookworm ESTs and transcripts, and from human transcriptome. Overall, the results indicate presence of nematode specific miRNA homologs in N. americanus and shades light on their putative targets in worm itself and the human host.

Identification and Characterization of miRNAs in Chondrus crispus by High-Throughput Sequencing and Bioinformatics Analysis

Chondrus crispus, an economically and medicinally important red alga, is a medicinally active substance and important for anti-tumor research. In this study, 117 C. crispus miRNAs (108 conserved and 9 novel) were identified from 2,416,181 small-RNA reads using high-throughput sequencing and bioinformatics methods. According to the BLAST search against the miRBase database, these miRNAs belonged to 110 miRNA families. Sequence alignment combined with homology searching revealed both the conservation and diversity of predicted potential miRNA families in different plant species. Four and 19 randomly selected miRNAs were validated by northern blotting and stem-loop quantitative real-time reverse transcription polymerase chain reaction detection, respectively. The validation rates (75% and 94.7%) demonstrated that most of the identified miRNAs could be credible. A total of 160 potential target genes were predicted and functionally annotated by Gene Ontology analysis and Kyoto Encyclopedia of Genes and Genomes analysis. We also analyzed the interrelationship of miRNAs, miRNA-target genes and target genes in C. crispus by constructing a Cytoscape network. The 117 miRNAs identified in our study should supply large quantities of information that will be important for red algae small RNA research.

Identification and Profiling of microRNAs Expressed in Elongating Cotton Fibers Using Small RNA Deep Sequencing

Plant microRNAs (miRNAs) have been shown to play essential roles in the regulation of gene expression. In this study, small RNA deep sequencing was applied to explore novel miRNAs expressed in elongating cotton fibers. A total of 46 novel and 96 known miRNAs, primarily derived from the corresponding specific loci in genome of Gossypium arboreum, were identified. 64 miRNAs were shown to be differentially expressed during the fiber elongation process; 16 were predicted to be novel miRNAs while the remaining 48 belong to known miRNA families. Furthermore, RLM-5' RACE (RNA ligase-mediated rapid amplification of 5'-cDNA ends) experiments identified the targets of eight important miRNAs, and the expression levels of these target genes were confirmed to be negatively correlated with the expression patterns of their corresponding miRNAs. We propose a potential functional network mediated through these eight miRNAs to illustrate their important functions in fiber elongation. Our study provides novel insights into the dynamic profiles of these miRNAs and a basis for investigating the regulatory mechanisms involved in the elongation of cotton fibers.

Transcriptome-wide identification of Rauvolfia serpentina microRNAs and prediction of their potential targets

MicroRNAs (miRNAs) are small non-coding RNAs of ∼ 19-24 nucleotides (nt) in length and considered as potent regulators of gene expression at transcriptional and post-transcriptional levels. Here we report the identification and characterization of 15 conserved miRNAs belonging to 13 families from Rauvolfia serpentina through in silico analysis of available nucleotide dataset. The identified mature R. serpentina miRNAs (rse-miRNAs) ranged between 20 and 22nt in length, and the average minimal folding free energy index (MFEI) value of rse-miRNA precursor sequences was found to be -0.815 kcal/mol. Using the identified rse-miRNAs as query, their potential targets were predicted in R. serpentina and other plant species. Gene Ontology (GO) annotation showed that predicted targets of rse-miRNAs include transcription factors as well as genes involved in diverse biological processes such as primary and secondary metabolism, stress response, disease resistance, growth, and development. Few rse-miRNAs were predicted to target genes of pharmaceutically important secondary metabolic pathways such as alkaloids and anthocyanin biosynthesis. Phylogenetic analysis showed the evolutionary relationship of rse-miRNAs and their precursor sequences to homologous pre-miRNA sequences from other plant species. The findings under present study besides giving first hand information about R. serpentina miRNAs and their targets, also contributes towards the better understanding of miRNA-mediated gene regulatory processes in plants.

MicroRNA and mRNA expression profiling analysis revealed the regulation of plant height in Gossypium hirsutum

BACKGROUND

Dwarf cottons are more resistant to damage from wind and rain and associated with stable, increased yields, and also desirable source for breeding the machine harvest varieties. In an effort to uncover the transcripts and miRNA networks involved in plant height, the transcriptome and small RNA sequencing were performed based on dwarf mutant Ari1327 (A1), tall-culm mutant Ari3697 (A3) and wild type Ari971 (A9) in Gossypium hirsutum.

METHODS

The stem apexes of wild-type upland cotton (Ari971) and its dwarf mutant (Ari1327) and tall-culm mutant (Ari3697) at the fifth true leaf stage were extracted for RNA, respectively. Transcriptome and small RNA libraries were constructed and subjected to next generation sequencing.

RESULTS

The transcriptome sequencing analysis showed that the enriched pathways of top 3 differentially expressed genes (DEGs) were categorized as carotenoid biosynthesis, plant-pathogen interaction and plant hormone signal transduction in both A1-A9 and A3-A9. The ABA and IAA related factors were differentially expressed in the mutants. Importantly, we found the lower expressed SAUR and elevated expressed GH3, and ABA related genes such as NCED and PP2C maybe relate to reduced growth of the plant height in Ari1327 which was consistent with the higher auxin and ABA content in this mutant. Furthermore, miRNA160 targeted to the auxin response factor (ARF) and miRNA166 (gma-miR166u and gma-miR166h-3p) targeted to ABA responsive element binding factor were related to the mutation in cotton. We have noticed that the cell growth related factors (smg7 targeted by gra-miR482 and 6 novel miRNAs and pectate-lyases targeted by osa-miR159f), the redox reactions related factors (Cytochrome P450 targeted by miR172) and MYB genes targeted by miR828, miR858 and miR159 were also involved in plant height of the cotton mutants. A total of 226 conserved miRNAs representing 32 known miRNA families were obtained, and 38 novel miRNAs corresponding to 23 unique RNA sequences were identified. Total 531 targets for 211 conserved miRNAs were obtained. Using PAREsnip, 27 and 29 miRNA/target conserved interactions were validated in A1-A9 and A3-A9, respectively. Furthermore, miRNA160, miRNA858 and miRNA172 were validated to be up-regulated in A1-A9 but down-regulated in A3-A9, whereas miRNA159 showed the opposite regulation.

CONCLUSIONS

This comprehensive interaction of the transcriptome and miRNA at tall-culm and dwarf mutant led to the discovery of regulatory mechanisms in plant height. It also provides the basis for in depth analyses of dwarf mutant genes for further breeding of dwarf cotton.

Aquaporins in Plants

Aquaporins are membrane channels that facilitate the transport of water and small neutral molecules across biological membranes of most living organisms. In plants, aquaporins occur as multiple isoforms reflecting a high diversity of cellular localizations, transport selectivity, and regulation properties. Plant aquaporins are localized in the plasma membrane, endoplasmic reticulum, vacuoles, plastids and, in some species, in membrane compartments interacting with symbiotic organisms. Plant aquaporins can transport various physiological substrates in addition to water. Of particular relevance for plants is the transport of dissolved gases such as carbon dioxide and ammonia or metalloids such as boron and silicon. Structure-function studies are developed to address the molecular and cellular mechanisms of plant aquaporin gating and subcellular trafficking. Phosphorylation plays a central role in these two processes. These mechanisms allow aquaporin regulation in response to signaling intermediates such as cytosolic pH and calcium, and reactive oxygen species. Combined genetic and physiological approaches are now integrating this knowledge, showing that aquaporins play key roles in hydraulic regulation in roots and leaves, during drought but also in response to stimuli as diverse as flooding, nutrient availability, temperature, or light. A general hydraulic control of plant tissue expansion by aquaporins is emerging, and their role in key developmental processes (seed germination, emergence of lateral roots) has been established. Plants with genetically altered aquaporin functions are now tested for their ability to improve plant tolerance to stresses. In conclusion, research on aquaporins delineates ever expanding fields in plant integrative biology thereby establishing their crucial role in plants.

Role of aquaporins in determining transpiration and photosynthesis in water-stressed plants: crop water-use efficiency, growth and yield

The global shortage of fresh water is one of our most severe agricultural problems, leading to dry and saline lands that reduce plant growth and crop yield. Here we review recent work highlighting the molecular mechanisms allowing some plant species and genotypes to maintain productivity under water stress conditions, and suggest molecular modifications to equip plants for greater production in water-limited environments. Aquaporins (AQPs) are thought to be the main transporters of water, small and uncharged solutes, and CO2 through plant cell membranes, thus linking leaf CO2 uptake from the intercellular airspaces to the chloroplast with water loss pathways. AQPs appear to play a role in regulating dynamic changes of root, stem and leaf hydraulic conductivity, especially in response to environmental changes, opening the door to using AQP expression to regulate plant water-use efficiency. We highlight the role of vascular AQPs in regulating leaf hydraulic conductivity and raise questions regarding their role (as well as tonoplast AQPs) in determining the plant isohydric threshold, growth rate, fruit yield production and harvest index. The tissue- or cell-specific expression of AQPs is discussed as a tool to increase yield relative to control plants under both normal and water-stressed conditions.

New Insights on Coffea miRNAs: Features and Evolutionary Conservation

Small RNAs influence the gene expression at the post-transcriptional level by guiding messenger RNA (mRNA) cleavage, translational repression, and chromatin modifications. In addition to model plants, the microRNAs (miRNAs) have been identified in different crop species. In this work, we developed a specific pipeline to search for coffee miRNA homologs on expressed sequence tags (ESTs) and genome survey sequences (GSS) databases. As a result, 36 microRNAs were identified and a total of 616 and 362 potential targets for Coffea arabica and Coffea canephora, respectively. The evolutionary analyses of these molecules were performed by comparing the primary and secondary structures of precursors and mature miRNAs with their orthologs. Moreover, using a stem-loop RT-PCR assay, we evaluated the accumulation of mature miRNAs in genomes with different ploidy levels, detecting an increase in the miRNAs accumulation according to the ploidy raising. Finally, a 5' RACE (Rapid Amplification of cDNA Ends) assay was performed to verify the regulation of auxin responsive factor 8 (ARF8) by MIR167 in coffee plants. The great variety of target genes indicates the functional plasticity of these molecules and reinforces the importance of understanding the RNAi-dependent regulatory mechanisms. Our results expand the study of miRNAs and their target genes in this crop, providing new challenges to understand the biology of these species.

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.