Small RNA sequencing identifies miRNA roles in ovule and fibre development

Integrated microRNA and whole-transcriptome sequencing reveals the involvement of small and long non-coding RNAs in the fiber growth of ramie plant

BACKGROUND

MicroRNAs (miRNAs) and long non-coding RNAs (lncRNAs) are the two main types of non-coding RNAs that play crucial roles in plant growth and development. However, their specific roles in the fiber growth of ramie plant (Boehmeria nivea L. Gaud) remain largely unknown.

METHODS

In this study, we performed miRNA and whole-transcriptome sequencing of two stem bark sections exhibiting different fiber growth stages to determine the expression profiles of miRNAs, lncRNAs, and protein-encoding genes.

RESULTS

Among the identified 378 miRNAs and 6,839 lncRNAs, 88 miRNAs and 1,288 lncRNAs exhibited differential expression. Bioinformatics analysis revealed that 29 and 228 differentially expressed protein-encoding genes were targeted by differentially expressed miRNAs and lncRNAs, respectively, constituting eight putative competing endogenous RNA networks. lncR00022274 exhibited downregulated expression in barks with growing fibers. It also had an antisense overlap with the MYB gene, BntWG10016451, whose overexpression drastically increased the xylem fiber number and secondary wall thickness of fibers in the stems of transgenic Arabidopsis, suggesting the potential association of lncR00022274-BntWG10016451 expression with fiber growth.

CONCLUSIONS

These findings provide insights into the roles of ncRNAs in the regulation of fiber growth in ramie, which can be used for the biotechnological improvement of its fiber yield and quality in the future.

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.

The cotton miR171a-SCL6 module mediates plant resistance through regulating GhPR1 expression

Plants have developed intricate defense mechanisms in response to fluctuating environmental cues, including the use of microRNA (miRNA) as post-transcriptional regulators. However, the specific mechanisms through which miRNA contributes to disease resistance remain largely elusive. While the miR171-SCLs have been investigated in an eclectic array of plants, there has been a notable scarcity of research specifically focused on cotton (Gossypium hirsutum). In our previous miRNA-sequencing analysis, we found that ghr-miR171a displayed a differential response to infections by Verticillium dahliae. In this study, we further investigated the function of the miR171a-SCL6 module in cotton during V. dahliae infection. The ghr-miR171a was confirmed to direct the cleavage of GhSCL6 mRNA in the post-transcriptional process, as evidenced by 5' RLM-RACE, β-glucuronidase (GUS) histochemical staining and enzyme activity assay. Interestingly, we found that overexpressing ghr-miR171a reduced cotton plants' resistance to V. dahliae, while suppressing ghr-miR171a increased the plants' defense capacity. The GhSCL6 protein, when fused with green fluorescent protein (GFP), localizes in the cell nucleus, indicating its potential role in gene regulation. This was further corroborated by yeast two-hybrid assays, which verified GhSCL6's transcriptional activation ability. Through quantitative reverse transcriptase PCR (qRT-PCR), luciferase (LUC) fluorescence, and yeast one-hybrid assays, we found that GhSCL6 binds to the GT-box element of the GhPR1 promoter, activating its expression and thereby enhancing plant disease resistance. Taken together, our findings demonstrate that the cotton miR171a-SCL6 module regulates Verticillium wilt resistance in plants through the post-transcriptional process. This insight may offer new perspectives for disease resistance strategies in cotton.

Insights to Gossypium defense response against Verticillium dahliae: the Cotton Cancer

The soil-borne pathogen Verticillium dahliae, also referred as "The Cotton Cancer," is responsible for causing Verticillium wilt in cotton crops, a destructive disease with a global impact. To infect cotton plants, the pathogen employs multiple virulence mechanisms such as releasing enzymes that degrade cell walls, activating genes that contribute to virulence, and using protein effectors. Conversely, cotton plants have developed numerous defense mechanisms to combat the impact of V. dahliae. These include strengthening the cell wall by producing lignin and depositing callose, discharging reactive oxygen species, and amassing hormones related to defense. Despite the efforts to develop resistant cultivars, there is still no permanent solution to Verticillium wilt due to a limited understanding of the underlying molecular mechanisms that drive both resistance and pathogenesis is currently prevalent. To address this challenge, cutting-edge technologies such as clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9), host-induced gene silencing (HIGS), and gene delivery via nano-carriers could be employed as effective alternatives to control the disease. This article intends to present an overview of V. dahliae virulence mechanisms and discuss the different cotton defense mechanisms against Verticillium wilt, including morphophysiological and biochemical responses and signaling pathways including jasmonic acid (JA), salicylic acid (SA), ethylene (ET), and strigolactones (SLs). Additionally, the article highlights the significance of microRNAs (miRNAs), circular RNAs (circRNAs), and long non-coding RNAs (lncRNAs) in gene expression regulation, as well as the different methods employed to identify and functionally validate genes to achieve resistance against this disease. Gaining a more profound understanding of these mechanisms could potentially result in the creation of more efficient strategies for combating Verticillium wilt in cotton crops.

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.

Whole-transcriptome analyses of Sorghum leaves identify key mRNAs and ncRNAs associated with GA3-mediated alleviation of salt stress

Sorghum has recently attracted much attention for its tolerance in high salt environment. However, the effect and regulatory mechanism of the gibberellic acid (GA₃)-mediated alleviation of salt stress in sorghum remains unclear. Herein, we reported that a GA₃ concentration of 50 mg/L is optimal for sorghum ("Jitian 3") development under salt stress. We conducted a whole-transcriptome analysis between GA₃-treated and control sorghum leaves under salt stress, and we identified 1002 differentially expressed (DE)-messenger RNAs (mRNAs), 81 DE-long non-coding RNAs (lncRNAs), 7 DE-circular RNAs (circRNAs), and 26 DE-microRNA (miRNAs) in sorghum following GA₃ treatment. We also identified a majority of DE-mRNAs and non-coding RNAs (ncRNAs) targets that serve essential roles in phenylpropanoid biosynthesis and plant hormone networks. In addition, we generated a competitive endogenous RNA (ceRNA)-miRNA-target gene network, and 3 circRNAs (circRNA_2746, circRNA_6515, circRNA_5622), 4 lncRNAs (XR_002450182.1, XR_002452422.1, XR_002448510.1, XR_002448296.1) and 4 genes (LOC8056546, LOC8062245, LOC8061469, LOC8071960) probably act as valuable candidates for the regulation of the GA₃-mediated alleviation of salt stress in sorghum. Our findings uncovered potential mRNA and non-coding RNAs that contribute to GA₃ regulation, thus offering a basis for the future investigation of underlying mechanisms of salt stress in sorghum.

RAD gene family analysis in cotton provides some key genes for flowering and stress tolerance in upland cotton G. hirsutum

Background

RADIALIS (RAD), belongs to the MYB gene family and regulates a variety of functions including floral dorsoventral asymmetry in Antirrhinum majus and development of fruit proteins in Solanum lycopersicum. RAD genes contain an SNF2_N superfamily domain. Here, we comprehensively identified 68 RAD genes from six different species including Arabidopsis and five species of cotton.

Results

Phylogenetic analysis classified RAD genes into five groups. Gene structure, protein motifs and conserved amino acid residues indicated that GhRAD genes were highly conserved during the evolutionary process. Chromosomal location information showed that GhRAD genes were distributed unevenly on different chromosomes. Collinearity and selection pressure analysis indicated RAD gene family expansion in G. hirsutum and G. barbadense with purifying selection pressure. Further, various growth and stress related promotor cis-acting elements were observed. Tissue specific expression level indicated that most GhRAD genes were highly expressed in roots and flowers (GhRAD2, GhRAD3, GhRAD4 and GhRAD11). Next, GhRAD genes were regulated by phytohormonal stresses (JA, BL and IAA). Moreover, Ghi-miRN1496, Ghi-miR1440, Ghi-miR2111b, Ghi-miR2950a, Ghi-miR390a, Ghi-miR390b and Ghi-miR7495 were the miRNAs targeting most of GhRAD genes.

Conclusions

Our study revealed that RAD genes are evolutionary conserved and might be involved in different developmental processes and hormonal stress response. Data presented in our study could be used as the basis for future studies of RAD genes in cotton.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08248-z.

Genome-wide identification and analysis of Oleosin gene family in four cotton species and its involvement in oil accumulation and germination

BACKGROUND

Cotton is not only a major textile fiber crop but also a vital oilseed, industrial, and forage crop. Oleosins are the structural proteins of oil bodies, influencing their size and the oil content in seeds. In addition, the degradation of oleosins is involved in the mobilization of lipid and oil bodies during seed germination. However, comprehensive identification and the systematic analysis of the Oleosin gene (OLEOs) family have not been conducted in cotton.

RESULTS

An in-depth analysis has enabled us to identify 25 and 24 OLEOs in tetraploid cotton species G. hirsutum and G. barbadense, respectively, while 12 and 13 OLEOs were identified in diploid species G. arboreum and G. raimondii, respectively. The 74 OLEOs were further clustered into three lineages according to the phylogenetic tree. Synteny analysis revealed that most of the OLEOs were conserved and that WGD or segmental duplications might drive their expansion. The transmembrane helices in GhOLEO proteins were predicted, and three transmembrane models were summarized, in which two were newly proposed. A total of 24 candidate miRNAs targeting GhOLEOs were predicted. Three highly expressed oil-related OLEOs, GH_A07G0501 (SL), GH_D10G0941 (SH), and GH_D01G1686 (U), were cloned, and their subcellular localization and function were analyzed. Their overexpression in Arabidopsis increased seed oil content and decreased seed germination rates.

CONCLUSION

We identified OLEO gene family in four cotton species and performed comparative analyses of their relationships, conserved structure, synteny, and gene duplication. The subcellular localization and function of three highly expressed oil-related OLEOs were detected. These results lay the foundation for further functional characterization of OLEOs and improving seed oil content.

Exploring diverse roles of micro RNAs in banana: Current status and future prospective

Micro RNAs (miRNAs) are 20-24 nucleotides long non-coding RNA sequences identified and characterized in multiple plant and animal systems. miRNAs play multifarious roles ranging from plant development to stress tolerance by synchronizing physiological processes at the level of transcription and translation. Banana is a major horticultural crop with colossal production worldwide. Despite the recent encouraging developments, the information on functions of miRNAs in banana plants is still in its infancy. The available literature pertaining to miRNAs in banana plants hints towards their contribution as master regulators in crucial physiological processes for instance abiotic stress responses, pathogenic defence response, fruit ripening and so on. This review is focused on biogenesis of miRNAs, their identification and deciphering their respective roles in banana plants with special emphasis on abiotic stress responses, plant immune responses, fruit ripening and storage. Based on the prior reports, we identified a few miRNAs with prospective roles in stress tolerance and illustrated the potential applications of miRNAs in banana crop improvement utilizing recent biotechnological tools such as CRISPR cas9, RNAi and the nano particle based foliar spray of miRNAs. The review briefly explained the future directions in banana research with a special emphasis on miRNA regulatory networks and agronomic traits improvement. Finally, future domains in miRNA research in plants and their possible applications towards crop improvement in agriculture are described briefly.

The emerging role of small RNAs in ovule development, a kind of magic

In plants, small RNAs have been recognized as key genetic and epigenetic regulators of development. Small RNAs are usually 20 to 30 nucleotides in length and they control, in a sequence specific manner, the transcriptional or post-transcriptional expression of genes. In this review, we present a comprehensive overview of the most recent findings about the function of small RNAs in ovule development, including megasporogenesis and megagametogenesis, both in sexual and apomictic plants. We discuss recent studies on the role of miRNAs, siRNAs and trans-acting RNAs (ta-siRNAs) in early female germline differentiation. The mechanistic complexity and unique regulatory features are reviewed, and possible directions for future research are provided.

Influence of virus-host interactions on plant response to abiotic stress

Environmental factors play a significant role in controlling growth, development and defense responses of plants. Changes in the abiotic environment not only significantly alter the physiological and molecular pathways in plants, but also result in attracting the insect pests that carry a payload of viruses. Invasion of plants by viruses triggers the RNA silencing based defense mechanism in plants. In counter defense the viruses have gained the ability to suppress the host RNA silencing activities. A new paradigm has emerged, with the recognition that plant viruses also have the intrinsic capacity to modulate host plant response to environmental cues, in an attempt to favour their own survival. Thus, plant-virus interactions provide an excellent system to understand the signals in crosstalk between biotic (virus) and abiotic stresses. In this review, we have summarized the basal plant defense responses to pathogen invasion while emphasizing on the role of RNA silencing as a front line of defense response to virus infection. The emerging knowledge indicates overlap between RNA silencing with the innate immune responses during antiviral defense. The suppressors of RNA silencing serve as Avr proteins, which can be recognized by the host R proteins. The defense signals also function in concert with the phytohormones to influence plant responses to abiotic stresses. The current evidence on the role of virus induced host tolerance to abiotic stresses is also discussed.

MicroRNAs play important roles in regulating the rapid growth of the Phyllostachys edulis culm internode

Moso bamboo (Phyllostachys edulis) is a fast-growing species with uneven growth and lignification from lower to upper segments within one internode. MicroRNAs (miRNAs) play a vital role in post-transcriptional regulation in plants. However, how miRNAs regulate fast growth in bamboo internodes is poorly understood. In this study, one moso bamboo internode was divided during early rapid growth into four segments called F4 (bottom) to F1 (upper) and these were then analysed for transcriptomes, miRNAs and degradomes. The F4 segment had a higher number of actively dividing cells as well as a higher content of auxin (IAA), cytokinin (CK) and gibberellin (GA) compared with the F1 segment. RNA-seq analysis showed DNA replication and cell division-associated genes highly expressed in F4 rather than in F1. In total, 63 miRNAs (DEMs) were identified as differentially expressed between F4 and F1. The degradome and the transcriptome indicated that many downstream transcription factors and hormonal responses genes were modulated by DEMs. Several miR-target interactions were further validated by tobacco co-infiltration. Our findings give new insights into miRNA-mediated regulatory pathways in bamboo, and will contribute to a comprehensive understanding of the molecular mechanisms governing rapid growth.

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 identification, gene cloning, subcellular location and expression analysis of SPL gene family in P. granatum L

BACKGROUNDS

Pomegranate is an excellent tree species with nutritional, medicinal, ornamental and ecological values. Studies have confirmed that SPL factors play an important role in floral transition and flower development.

RESULTS

Used bioinformatics methods, 15 SPL (SQUAMOSA promoter-binding protein-like) genes were identified and analyzed from the 'Taishanhong' pomegranate (P. granatum L.) genome. Phylogenetic analysis showed that PgSPLs were divided into six subfamilies (G1 ~ G6). PgSPL promoter sequences contained multiple cis-acting elements associated with abiotic stress or hormonal response. Based on the transcriptome data, expression profiles of different tissues and different developmental stages showed that PgSPL genes had distinct temporal and spatial expression characteristics. The expression analysis of miR156 in small RNA sequencing results showed that miR156 negatively regulated the expression of target genes. qRT-PCR analysis showed that the expression levels of PgSPL2, PgSPL3, PgSPL6, PgSPL11 and PgSPL14 in leaves were significantly higher than those in buds and stems (p < 0.05). The expression levels of PgSPL5, PgSPL12 and PgSPL13 in flower buds were significantly higher than that in leaves and stems (p < 0.05). The full-length of coding sequence of PgSPL5 and PgSPL13 were obtained by homologous cloning technology. The full length of PgSPL5 is 1020 bp, and PgSPL13 is 489 bp, which encodes 339 and 162 amino acids, respectively. Further investigation revealed that PgSPL5 and PgSPL13 proteins were located in the nucleus. Exogenous plant growth regulator induction experiments showed that PgSPL5 was up-regulated in leaves and stems. PgSPL13 was up-regulated in leaves and down-regulated in stems. When sprayed with 6-BA, IBA and PP333 respectively, PgSPL5 and PgSPL13 were up-regulated most significantly at P2 (bud vertical diameter was 5.1 ~ 12.0 mm) stage of bisexual and functional male flowers.

CONCLUSIONS

Our findings suggested that PgSPL2, PgSPL3, PgSPL6, PgSPL11 and PgSPL14 played roles in leaves development of pomegranate. PgSPL5, PgSPL12 and PgSPL13 played roles in pomegranate flower development. PgSPL5 and PgSPL13 were involved in the response process of different plant hormone signal transduction in pomegranate development. This study provided a robust basis for further functional analyses of SPL genes in pomegranate.

Non-coding RNAs in the interaction between rice and Meloidogyne graminicola

Background

Root knot nematodes (RKN) are plant parasitic nematodes causing major yield losses of widely consumed food crops such as rice (Oryza sativa). Because non-coding RNAs, including small interfering RNAs (siRNA), microRNAs (miRNAs) and long non-coding RNAs (lncRNAs), are key regulators of various plant processes, elucidating their regulation during this interaction may lead to new strategies to improve crop protection. In this study, we aimed to identify and characterize rice siRNAs, miRNAs and lncRNAs responsive to early infection with RKN Meloidogyne graminicola (Mg), based on sequencing of small RNA, degradome and total RNA libraries from rice gall tissues compared with uninfected root tissues.

Results

We found 425 lncRNAs, 3739 siRNAs and 16 miRNAs to be differentially expressed between both tissues, of which a subset was independently validated with RT-qPCR. Functional prediction of the lncRNAs indicates that a large part of their potential target genes code for serine/threonine protein kinases and transcription factors. Differentially expressed siRNAs have a predominant size of 24 nts, suggesting a role in DNA methylation. Differentially expressed miRNAs are generally downregulated and target transcription factors, which show reduced degradation according to the degradome data.

Conclusions

To our knowledge, this work is the first to focus on small and long non-coding RNAs in the interaction between rice and Mg, and provides an overview of rice non-coding RNAs with the potential to be used as a resource for the development of new crop protection strategies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07735-7.

Comparative analysis of sRNAs, degradome and transcriptomics in sweet sorghum reveals the regulatory roles of miRNAs in Cd accumulation and tolerance

Key miRNAs including sbi-miR169p/q, sbi-miR171g/j, sbi-miR172a/c/d, sbi-miR172e, sbi-miR319a/b, sbi-miR396a/b, miR408, sbi-miR5384, sbi-miR5565e and nov_23 were identified to function in the regulation of Cd accumulation and tolerance. As an energy plant, sweet sorghum shows great potential in the phytoremediation of Cd-contaminated soils. However, few studies have focused on the regulatory roles of miRNAs and their targets under Cd stress. In this study, comparative analysis of sRNAs, degradome and transcriptomics was conducted in high-Cd accumulation (H18) and low-Cd accumulation (L69) genotypes of sweet sorghum. A total of 38 conserved and 23 novel miRNAs with differential expressions were identified under Cd stress or between H18 and L69, and 114 target genes of 41 miRNAs were validated. Furthermore, 25 miRNA-mRNA pairs exhibited negatively correlated expression profiles and sbi-miR172e together with its target might participate in the distinct Cd tolerance between H18 and L69 as well as sbi-miR172a/c/d. Additionally, two groups of them: miR169p/q-nov_23 and miR408 were focused through the co-expression analysis, which might be involved in Cd uptake and tolerance by regulating their targets associated with transmembrane transportation, cytoskeleton activity, cell wall construction and ROS (reactive oxygen species) homeostasis. Further experiments exhibited that cell wall components of H18 and L69 were different when exposed to cadmium, which might be regulated by miR169p/q, miR171g/j, miR319a/b, miR396a/b, miR5384 and miR5565e through their targets. Through this study, we aim to reveal the potential miRNAs involved in sweet sorghum in response to Cd stress and provide references for developing high-Cd accumulation or high Cd-resistant germplasm of sweet sorghum that can be used in phytoremediation.

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.

MicroRNAs: Potential Targets for Developing Stress-Tolerant Crops

Crop yield is challenged every year worldwide by changing climatic conditions. The forecasted climatic scenario urgently demands stress-tolerant crop varieties to feed the ever-increasing global population. Molecular breeding and genetic engineering approaches have been frequently exploited for developing crops with desired agronomic traits. Recently, microRNAs (miRNAs) have emerged as powerful molecules, which potentially serve as expression markers during stress conditions. The miRNAs are small non-coding endogenous RNAs, usually 20-24 nucleotides long, which mediate post-transcriptional gene silencing and fine-tune the regulation of many abiotic- and biotic-stress responsive genes in plants. The miRNAs usually function by specifically pairing with the target mRNAs, inducing their cleavage or repressing their translation. This review focuses on the exploration of the functional role of miRNAs in regulating plant responses to abiotic and biotic stresses. Moreover, a methodology is also discussed to mine stress-responsive miRNAs from the enormous amount of transcriptome data available in the public domain generated using next-generation sequencing (NGS). Considering the functional role of miRNAs in mediating stress responses, these molecules may be explored as novel targets for engineering stress-tolerant crop varieties.

Functional Exploration of Chaperonin (HSP60/10) Family Genes and their Abiotic Stress-induced Expression Patterns in Sorghum bicolor

Background

Sorghum, the C4 dry-land cereal, important for food, fodder, feed and fuel, is a model crop for abiotic stress tolerance with smaller genome size, genetic diversity, and bio-energy traits. The heat shock proteins/chaperonin 60s (HSP60/Cpn60s) assist the plastid proteins, and participate in the folding and aggregation of proteins. However, the functions of HSP60s in abiotic stress tolerance in Sorghum remain unclear.

Methods

Genome-wide screening and in silico characterization of SbHSP60s were carried out along with tissue and stress-specific expression analysis.

Results

A total of 36 HSP60 genes were identified in Sorghum bicolor. They were subdivided into 2 groups, the HSP60 and HSP10 co-chaperonins encoded by 30 and 6 genes, respectively. The genes are distributed on all the chromosomes, chromosome 1 being the hot spot with 9 genes. All the HSP60s were found hydrophilic and highly unstable. The HSP60 genes showed a large number of introns, the majority of them with more than 10. Among the 12 paralogs, only 1 was tandem and the remaining 11 segmental, indicating their role in the expansion of SbHSP60s. Majority of the SbHSP60 genes expressed uniformly in leaf while a moderate expression was observed in the root tissues, with the highest expression displayed by SbHSP60-1. From expression analysis, SbHSP60-3 for drought, SbHSP60-9 for salt, SbHSP60-9 and 24 for heat and SbHSP60-3, 9 and SbHSP10-2 have been found implicated for cold stress tolerance and appeared as the key regulatory genes.

Conclusion

This work paves the way for the utilization of chaperonin family genes for achieving abiotic stress tolerance in plants.

Small RNA and degradome deep sequencing reveals important roles of microRNAs in cotton (Gossypium hirsutum L.) response to root-knot nematode Meloidogyne incognita infection

Investigation of cotton response to nematode infection will allow us to better understand the cotton immune defense mechanism and design a better biotechnological approach for efficiently managing pest nematodes in cotton. In this study, we firstly treated cotton by root knot nematode (RKN, Meloidogyne incognita) infections, then we employed the high throughput deep sequencing technology to sequence and genome-widely identify all miRNAs in cotton; finally, we analyzed the functions of these miRNAs in cotton response to RKN infections. A total of 266 miRNAs, including 193 known and 73 novel miRNAs, were identified by deep sequencing technology, which belong to 67 conserved and 66 novel miRNA families, respectively. A majority of identified miRNA families only contain one miRNA; however, miR482 family contains 14 members and some others contain 2-13 members. Certain miRNAs were specifically expressed in RKN-infected cotton roots and others were completely inhibited by RKN infection. A total of 50 miRNAs were differentially expressed after RKN infection, in which 28 miRNAs were up-regulated and 22 were inhibited by RKN treatment. Based on degradome sequencing, 87 gene targets were identified to be targeted by 57 miRNAs. These miRNA-targeted genes are involved in the interaction of cotton plants and nematode infection. Based on GO (gene ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis, 466 genes from all 636 miRNA targets were mapped to 6340 GO terms, 181 genes from 228 targets of differentially expressed miRNAs were mapped to 1588 GO terms. The GO terms were then categorized into the three main GO classes: biological processes, cellular components, and molecular functions. The targets of differentially expressed miRNAs were enriched in 43 GO terms, including 22 biological processes, 10 cellular components, and 11 molecular functions (p < 0.05). Many identified processes were associated with organism responses to the environmental stresses, including regulation of nematode larval development, response to nematode, and response to flooding. Our results will enhance the study and application of developing new cotton cultivars for nematode resistance.

Profiling of MicroRNAs Involved in Mepiquat Chloride-Mediated Inhibition of Internode Elongation in Cotton (Gossypium hirsutum L.) Seedlings

Mepiquat chloride (MC) is the most important plant growth retardant that is widely used in cotton (Gossypium hirsutum L.) production to suppress excessive vegetative growth and improve plant architecture. MicroRNAs (miRNAs) are important gene expression regulators that control plant growth and development. However, miRNA-mediated post-transcriptional regulation in MC-induced growth inhibition remains unclear. In this study, the dynamic expression profiles of miRNAs responsive to MC in cotton internodes were investigated. A total of 508 known miRNAs belonging to 197 families and five novel miRNAs were identified. Among them, 104 miRNAs were differentially expressed at 48, 72, or 96 h post MC treatment compared with the control (0 h); majority of them were highly conserved miRNAs. The number of differentially expressed miRNAs increased with time after treatment. The expression of 14 known miRNAs was continuously suppressed, whereas 12 known miRNAs and one novel miRNA were continuously induced by MC. The expression patterns of the nine differentially expressed miRNAs were verified using qRT-PCR. The targets of the known and novel miRNAs were predicted. Four conserved and six novel targets were validated using the RLM-5' RACE assay. This study revealed that miRNAs play crucial regulatory roles in the MC-induced inhibition of internode elongation. It can improve our understanding of post-transcriptional gene regulation in MC-mediated growth inhibition and could potentially facilitate the breeding of dwarf cotton.

Role of phasiRNAs from two distinct phasing frames of GhMYB2 loci in cis- gene regulation in the cotton genome

BACKGROUND

Phased small interfering RNA (phasiRNA) is primarily derived from the 22-nt miRNA targeting loci. GhMYB2, a gene with potential roles in cotton fiber cell fate determination, is a target gene of miR828 and miR858 in the generation of phasiRNAs.

RESULTS

In the presented work, through the evaluation of phasing scores and phasiRNA distribution pattern, we found that phasiRNAs from GhMYB2 were derived from the 3' cleavage fragments of 22-nt miR828 and 21-nt miR858 respectively. These two miRNA targeting sites initiated two phasing frames on transcripts of one locus. By means of RNA ligase-mediated rapid amplification of cDNA ends (RLM-RACE), we further demonstrated that phasiRNAs derived from the two phasing frames played a role in cis-regulation of GhMYB2. The phasiRNAs derived from GhMYB2 were expressed in the somatic tissues, especially in anther and hypocotyl. We further employed our previous small RNA sequencing data as well as the degradome data of cotton fiber bearing ovules, anthers, hypocotyls and embryogenic calli tissues published in public databases, to validate the expression, phasing pattern and functions of phasiRNAs.

CONCLUSIONS

The presenting research provide insights of the molecular mechanism of phasiRNAs in regulation of GhMYB2 loci.

Engineering Multiple Abiotic Stress Tolerance in Canola, Brassica napus

Impacts of climate change like global warming, drought, flooding, and other extreme events are posing severe challenges to global crop production. Contribution of Brassica napus towards the oilseed industry makes it an essential component of international trade and agroeconomics. Consequences from increasing occurrences of multiple abiotic stresses on this crop are leading to agroeconomic losses making it vital to endow B. napus crop with an ability to survive and maintain yield when faced with simultaneous exposure to multiple abiotic stresses. For an improved understanding of the stress sensing machinery, there is a need for analyzing regulatory pathways of multiple stress-responsive genes and other regulatory elements such as non-coding RNAs. However, our understanding of these pathways and their interactions in B. napus is far from complete. This review outlines the current knowledge of stress-responsive genes and their role in imparting multiple stress tolerance in B. napus. Analysis of network cross-talk through omics data mining is now making it possible to unravel the underlying complexity required for stress sensing and signaling in plants. Novel biotechnological approaches such as transgene-free genome editing and utilization of nanoparticles as gene delivery tools are also discussed. These can contribute to providing solutions for developing climate change resilient B. napus varieties with reduced regulatory limitations. The potential ability of synthetic biology to engineer and modify networks through fine-tuning of stress regulatory elements for plant responses to stress adaption is also highlighted.

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.

Small RNA and mRNA Sequencing Reveal the Roles of microRNAs Involved in Pomegranate Female Sterility

Female sterility is a key factor restricting plant reproduction. Our previous studies have revealed that pomegranate female sterility mainly arose from the abnormality of ovule development. MicroRNAs (miRNAs) play important roles in ovule development. However, little is known about the roles of miRNAs in female sterility. In this study, a combined high-throughput sequencing approach was used to investigate the miRNAs and their targeted transcripts involved in female development. A total of 103 conserved and 58 novel miRNAs were identified. Comparative profiling indicated that the expression of 43 known miRNAs and 14 novel miRNAs were differentially expressed between functional male flowers (FMFs) and bisexual flowers (BFs), 30 known miRNAs and nine novel miRNAs showed significant differences among different stages of BFs, and 20 known miRNAs and 18 novel miRNAs exhibited remarkable expression differences among different stages of FMFs. Gene ontology (GO) analyses of 144 predicted targets of differentially expressed miRNAs indicated that the “reproduction process” and “floral whorl development” processes were significantly enriched. The miRNA–mRNA interaction analyses revealed six pairs of candidate miRNAs and their targets associated with female sterility. Interestingly, pg-miR166a-3p was accumulated, whereas its predicted targets (Gglean012177.1 and Gglean013966.1) were repressed in functional male flowers (FMFs), and the interaction between pg-miR166a-3p and its targets (Gglean012177.1 and Gglean013966.1) were confirmed by transient assay. A. thaliana transformed with 35S-pre-pg-miR166a-3p verified the role of pg-miR166a-3p in ovule development, which indicated pg-miR166a-3p’s potential role in pomegranate female sterility. The results provide new insights into molecular mechanisms underlying the female sterility at the miRNA level.

Small RNA sequencing identifies cucumber miRNA roles in waterlogging-triggered adventitious root primordia formation

The formation of adventitious roots (ARs) is a key morphological adaptation of cucumber (Cucumis sativus L.) to waterlogging stress. MicroRNAs (miRNAs) constitute a group of non-coding small RNAs (sRNA) that play crucial roles in regulating diverse biological processes, including waterlogging acclimation. However, which specific miRNAs and how they are involved in waterlogging-triggered de novo AR primordia formation are not fully known. Here, Illumina sRNA sequencing was applied to sequence six sRNA libraries generated from the waterlogging-tolerant cucumber Zaoer-N after 48 h of waterlogging and the control. A total of 358 cucumber miRNAs, 312 known and 46 novel, were obtained. Among them, 23 were differentially expressed, with 10 and 13 being up- and downregulated, respectively. A qPCR expression study confirmed that the identified differentially expressed miRNAs were credible. A total of 657 putative miRNA target genes were predicted for the 23 miRNAs using an in silico approach. A gene ontology enrichment analysis revealed that target genes functioning in cell redox homeostasis, cytoskeleton, photosynthesis and cell growth were over-represented. In total, 58 of the 657 target genes showed inverse expression patterns compared with their respective miRNAs through a combined analysis of sRNA- and RNA-sequencing-based transcriptome datasets using the same experimental design. The target gene annotation included a peroxidase, a GDSL esterases/lipase and two heavy metal-associated isoprenylated plant proteins. Our results provide an important framework for understanding the unique miRNA patterns seen in responses to waterlogging and the miRNA-mediated formation of de novo AR primordia in cucumber.

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.

Electronic supplementary material

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

The effect of DNA methylation on the miRNA expression pattern in lipopolysaccharide-induced inflammatory responses in human dental pulp cells

Endodontic infection is a widespread oral problem. DNA methylation is a key epigenetic modification that plays important roles in various inflammatory responses, but its role in dental pulp inflammation is poorly understood. In this study, we assessed the expression of DNA methyltransferases (DNMTs) in human dental pulp cells (hDPCs) during lipopolysaccharide (LPS)-induced inflammation and found that DNMT3B mRNA expression was reduced and DNMT1 mRNA and protein levels decreased significantly. Pretreatment with the DNMT inhibitor 5-Aza-2'-deoxycytidine (5-Aza-CdR) significantly enhanced the expression of the inflammatory cytokines IL-6 and IL-8 in LPS-stimulated hDPCs, indicating that DNA methylation may play a role in hDPC inflammation. Studies have reported that some microRNAs (miRNAs) are involved in dental pulp infection. DNA methylation can modulate the inflammatory response by regulating miRNA expression, but this phenomenon has not yet been reported in pulp inflammation. The present study used next-generation sequencing to examine the effect of 5-Aza-CdR on the miRNA expression profile of LPS-treated hDPCs, and the results showed that 5-Aza-CdR pretreatment changed the miRNA expression pattern in hDPCs during inflammation. Among the changed miRNAs, miR-146a-5p, which is a pulp inflammation-related miRNA, demonstrated the most noticeably altered expression. miR-146a-5p could be induced by LPS in hDPCs, and 5-Aza-CdR preincubation or DNMT1 knockdown markedly increased its expression level. However, no significant difference was found in the methylation pattern of the MIR146A promoter with 5-Aza-CdR pretreatment or DNMT1 knockdown in LPS-stimulated hDPCs. These results indicate that DNA methylation may regulate the LPS-induced inflammatory response by changing the miRNA expression in hDPCs.

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.

Whirly (Why) transcription factors in tomato (Solanum lycopersicum L.): genome-wide identification and transcriptional profiling under drought and salt stresses

Whirly (Why) transcription factor (TFs) constitute one of the important TF families which plays essential roles in plant metabolism to cope with environmental stresses. In the present study, Why genes were identified at genome-wide scale in tomato (Solanum lycopersicum), and bioinformatics analyses were implemented. Validation of Why genes expressions under drought and salt stresses were also performed using RT-qPCR. The analyses revealed the presence of two Why genes in tomato genome, SlWhy1 (Solyc05g007100.2.1) and SlWhy2 (Solyc11g044750.1.1). Both genes contained Whirly transcription factor domain structure (PF08536), and Why proteins were in basic character (pI ≥ 7). While the lengths of the proteins ranged from 268 to 236 amino acid residues for SlWhy1 and SlWhy2 respectively, exon numbers identified in both genes were seven. According to the digital expression data, SlWhy genes are expressed at medium level in different anatomical parts and developmental stages. In the promotor sequence analysis, 13 types of putative TF binding sites were identified, and the highest motif number was 46, found for GATA TF. Gene co-expression analyses revealed that complex networks for SlWhy genes, which are connected with various metabolic pathways. Based on the RT-qPCR data, both SlWhy1 and SlWhy2 genes were up-regulated under salt and drought stresses. 3D structure analyses revealed that SlWhy1 protein had a more diverged structure than SlWhy2 protein, based on their comparisons in Arabidopsis and potato. The results obtained in the present study could be a useful scientific basis for understanding Why genes in tomato and their functions under abiotic stress conditions.

Identification of lncRNAs involved in rice ovule development and female gametophyte abortion by genome-wide screening and functional analysis

Background

As important female reproductive tissues, the rice (Oryza sativa L.) ovule and female gametophyte is significant in terms of their fertility. Long noncoding RNAs (lncRNAs) play important and wide-ranging roles in the growth and development of plants and have become a major research focus in recent years. Therefore, we explored the characterization and expression change of lncRNAs during ovule development and female gametophytic abortion.

Results

In our study, whole-transcriptome strand-specific RNA sequencing (ssRNA-seq) was performed in the ovules of a high-frequency female-sterile rice line (fsv1) and a wild-type rice line (Gui99) at the megaspore mother cell meiosis stage (stage 1), functional megaspore mitosis stage (stage 2) and female gametophyte mature stage (stage 3). By comparing two rice lines, we identified 152, 233, and 197 differentially expressed lncRNAs at the three ovule developmental stages. Functional analysis of the coherent target genes of these differentially expressed lncRNAs indicated that many lncRNAs participate in multiple pathways such as hormone and cellular metabolism and signal transduction. Moreover, there were many differentially expressed lncRNAs acting as the precursors of some miRNAs that are involved in the development of ovules and female gametophytes. In addition, we have found that lncRNAs can act as decoys, competing with mRNAs for binding to miRNAs to maintain the normal expression of genes related to ovule and female gametophyte development.

Conclusion

These results provide important clues for elucidating the female gametophyte abortion mechanism in rice. This study also expands our understanding about the biological functions of lncRNAs and the annotation of the rice genome.

Electronic supplementary material

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

miRLocator: A Python Implementation and Web Server for Predicting miRNAs from Pre-miRNA Sequences

microRNAs (miRNAs) are short, noncoding regulatory RNAs derived from hairpin precursors (pre-miRNAs). In synergy with experimental approaches, computational approaches have become an invaluable tool for identifying miRNAs at the genome scale. We have recently reported a method called miRLocator, which applies machine learning algorithms to accurately predict the localization of most likely miRNAs within their pre-miRNAs. One major strength of miRLocator is the fact that the machine learning-based miRNA prediction model can be automatically trained using a set of miRNAs of particular interest, with informative features extracted from miRNA-miRNA* duplexes and the optimized ratio between positive and negative samples. Here, we present a detailed protocol for miRLocator that performs the training and prediction processes using a python implementation and web interface. The source codes, web interface, and manual documents are freely available to academic users at https://github.com/cma2015/miRLocator .

Molecular Control by Non-coding RNAs During Fruit Development: From Gynoecium Patterning to Fruit Ripening

Fruits are originated from the transition of a quiescent ovary to a fast-growing young fruit. The evolution of reproductive structures such as ovary and fruit has made seed dispersal easier, which is a key process for reproductive success in flowering plants. The complete fruit development and ripening are characterized by a remarkable phenotypic plasticity which is orchestrated by a myriad of genetic factors. In this context, transcriptional regulation by non-coding small (i.e., microRNAs) and long (lncRNAs) RNAs underlies important mechanisms controlling reproductive organ development. These mechanisms may act together and interact with other pathways (i.e., phytohormones) to regulate cell fate and coordinate reproductive organ development. Functional genomics has shown that non-coding RNAs regulate a diversity of developmental reproductive stages, from carpel formation and ovary development to the softening of the ripe/ripened fruit. This layer of transcriptional control has been associated with ovule, seed, and fruit development as well as fruit ripening, which are crucial developmental processes in breeding programs because of their relevance for crop production. The final ripe fruit is the result of a process under multiple levels of regulation, including mechanisms orchestrated by microRNAs and lncRNAs. Most of the studies we discuss involve work on tomato and Arabidopsis. In this review, we summarize non-coding RNA-controlled mechanisms described in the current literature that act coordinating the main steps of gynoecium development/patterning and fruit ripening.

The role of miRNA in somatic embryogenesis

Somatic embryogenesis (SEG) is one of the best techniques for mass production of economically important plants. It is also used for the study of morphology, anatomy, physiology, genetics and molecular mechanism of embryo development. Somatic Embryos (SE) are bipolar structures that develop from a cell other than a gamete or zygote. SEG reflects the unique developmental potential of plant somatic cells, resulting in the transition of the differentiated somatic cells to embryogenic cells to follow the zygotic embryo stages. There are several biochemical and physiological processes that transformed a single somatic cell to a whole plant. SE studies provide insight into cell mechanisms governing the totipotency process in plants. Previously, in vitro studies have suggested the role of various regulatory genes in embryogenic transition that are triggered by plant hormones in response to stress. The omic studies identify the specific genes, transcripts, and proteins required for somatic embryogenesis development. MicroRNAs (miRNAs) are small, 19-24 nucleotides (nt), non-coding small RNA regulatory molecules controlling a large number of biological processes. In addition to their role in SEG, miRNAs play vital role in plant development, secondary metabolite synthesis and metabolism of macromolecules, hormone signal transduction, and tolerance of plants to biotic and abiotic stresses. During last decade several types of miRNAs involved in SEG have been reported. Among these miRNAs, miR156, miR162, miR166a, miR167, miR168, miR171a/b, miR171c, miR393, miR397 and miR398 played very active role during various stages of SEG. In this review, we highlighted the role of these as well as other miRNAs in some economically important plants.

A critical and speculative review on microRNA technology in crop improvement: Current challenges and future directions

MicroRNAs (miRNAs) lie at the center of gene regulation and, as such, have become novel targets for crop improvement including the enhancement of crop quality and yields as well as responses to environmental stresses. There are several major issues related to miRNA technology including the functional analysis of miRNAs and their nomenclature. In this critical and speculative review, we recommend several directions for future plant miRNA research and perspectives. Research on miRNA needs to be extended from merely descriptive studies to functional studies. More genetic tools, such as genome editing, should be developed for miRNA functional study. Obtaining transgenic plants is a bottleneck for plant miRNA functional studies and, hence, more reliable transformation methods need to be developed. We also propose a new terminology approach for miRNA nomenclature. The current miRNA nomenclature is confusing and has mislead much research. Here we suggest to name a miRNA as miR#-5p or -3p, and to name their opposite strand as miR#*-3p or -5p. The advantages of the new nomenclature is that it covers information on the history, relationship, family, and location of an individual miRNA. It recognizes both traditional and new discovery.

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.

MicroRNA-target gene responses to root knot nematode (Meloidogyne incognita) infection in cotton (Gossypium hirsutum L.)

MicroRNAs (miRNAs) are a large class of small regulatory RNA molecules, however no study has been performed to elucidate the role of miRNAs in cotton (Gossypium hirsutum) response to the root knot nematode (RKN, Meloidogyne incognita) infection. We selected 28 miRNAs and 8 miRNA target genes to investigate the miRNA-target gene response to M. incognita infection. Our results show that RKN infection significantly affected the expression of several miRNAs and their targeted genes. After 10 days of RKN infection, expression fold changes on miRNA expressions ranged from down-regulated by 33% to upregulated by 406%; meanwhile the expression levels of miRNA target genes were 45.8% to 231%. Three miRNA-target pairs, miR159-MYB, miR319-TCP4 and miR167-ARF8, showed inverse expression patterns between gene targets and their corresponded miRNAs, suggesting miRNA-mediated gene regulation in cotton roots in response to RKN infection.

Characterization of the late embryogenesis abundant (LEA) proteins family and their role in drought stress tolerance in upland cotton

Background

Late embryogenesis abundant (LEA) proteins are large groups of hydrophilic proteins with major role in drought and other abiotic stresses tolerance in plants. In-depth study and characterization of LEA protein families have been carried out in other plants, but not in upland cotton. The main aim of this research work was to characterize the late embryogenesis abundant (LEA) protein families and to carry out gene expression analysis to determine their potential role in drought stress tolerance in upland cotton. Increased cotton production in the face of declining precipitation and availability of fresh water for agriculture use is the focus for breeders, cotton being the backbone of textile industries and a cash crop for many countries globally.

Results

In this work, a total of 242, 136 and 142 LEA genes were identified in G. hirsutum, G. arboreum and G. raimondii respectively. The identified genes were classified into eight groups based on their conserved domain and phylogenetic tree analysis. LEA 2 were the most abundant, this could be attributed to their hydrophobic character. Upland cotton LEA genes have fewer introns and are distributed in all chromosomes. Majority of the duplicated LEA genes were segmental. Syntenic analysis showed that greater percentages of LEA genes are conserved. Segmental gene duplication played a key role in the expansion of LEA genes. Sixty three miRNAs were found to target 89 genes, such as miR164, ghr-miR394 among others. Gene ontology analysis revealed that LEA genes are involved in desiccation and defense responses. Almost all the LEA genes in their promoters contained ABRE, MBS, W-Box and TAC-elements, functionally known to be involved in drought stress and other stress responses. Majority of the LEA genes were involved in secretory pathways. Expression profile analysis indicated that most of the LEA genes were highly expressed in drought tolerant cultivars Gossypium tomentosum as opposed to drought susceptible, G. hirsutum. The tolerant genotypes have a greater ability to modulate genes under drought stress than the more susceptible upland cotton cultivars.

Conclusion

The finding provides comprehensive information on LEA genes in upland cotton, G. hirsutum and possible function in plants under drought stress.

Electronic supplementary material

The online version of this article (10.1186/s12863-017-0596-1) contains supplementary material, which is available to authorized users.

Characterization of the late embryogenesis abundant (LEA) proteins family and their role in drought stress tolerance in upland cotton

BACKGROUND

Late embryogenesis abundant (LEA) proteins are large groups of hydrophilic proteins with major role in drought and other abiotic stresses tolerance in plants. In-depth study and characterization of LEA protein families have been carried out in other plants, but not in upland cotton. The main aim of this research work was to characterize the late embryogenesis abundant (LEA) protein families and to carry out gene expression analysis to determine their potential role in drought stress tolerance in upland cotton. Increased cotton production in the face of declining precipitation and availability of fresh water for agriculture use is the focus for breeders, cotton being the backbone of textile industries and a cash crop for many countries globally.

RESULTS

In this work, a total of 242, 136 and 142 LEA genes were identified in G. hirsutum, G. arboreum and G. raimondii respectively. The identified genes were classified into eight groups based on their conserved domain and phylogenetic tree analysis. LEA 2 were the most abundant, this could be attributed to their hydrophobic character. Upland cotton LEA genes have fewer introns and are distributed in all chromosomes. Majority of the duplicated LEA genes were segmental. Syntenic analysis showed that greater percentages of LEA genes are conserved. Segmental gene duplication played a key role in the expansion of LEA genes. Sixty three miRNAs were found to target 89 genes, such as miR164, ghr-miR394 among others. Gene ontology analysis revealed that LEA genes are involved in desiccation and defense responses. Almost all the LEA genes in their promoters contained ABRE, MBS, W-Box and TAC-elements, functionally known to be involved in drought stress and other stress responses. Majority of the LEA genes were involved in secretory pathways. Expression profile analysis indicated that most of the LEA genes were highly expressed in drought tolerant cultivars Gossypium tomentosum as opposed to drought susceptible, G. hirsutum. The tolerant genotypes have a greater ability to modulate genes under drought stress than the more susceptible upland cotton cultivars.

CONCLUSION

The finding provides comprehensive information on LEA genes in upland cotton, G. hirsutum and possible function in plants under drought stress.

High throughput deep sequencing reveals the important roles of microRNAs during sweetpotato storage at chilling temperature

Sweetpotato (Impomoea batatas L.) is a globally important economic food crop with a potential of becoming a bioenergy and pharmaceutical crop. Thus, studying the molecular mechanism of tuberous root development and storage is very important. However, not too much progress has been made in this field. In this study, we employed the next generation high-throughput deep sequencing technology to sequence all small RNAs and degradome of sweetpotato for systematically investigating sweetpotato response to chilling stress during storage. A total of 190 known microRNAs (miRNAs) and 191 novel miRNAs were identified, and 428 transcripts were targeted by 184 identified miRNAs. More importantly, we identified 26 miRNAs differentially expressed between chilling stress and control conditions. The expression of these miRNAs and their targets was also confirmed by qRT-PCR. Integrated analysis of small RNAs and degradome sequencing reveals that miRNA-mediated SA signaling, ABA-dependent, and ROS response pathways are involved in sweetpotato root response to chilling stress during storage.

Genome-Wide Analysis of DNA Methylation During Ovule Development of Female-Sterile Rice fsv1

The regulation of female fertility is an important field of rice sexual reproduction research. DNA methylation is an essential epigenetic modification that dynamically regulates gene expression during development processes. However, few reports have described the methylation profiles of female-sterile rice during ovule development. In this study, ovules were continuously acquired from the beginning of megaspore mother cell meiosis until the mature female gametophyte formation period, and global DNA methylation patterns were compared in the ovules of a high-frequency female-sterile line (fsv1) and a wild-type rice line (Gui99) using whole-genome bisulfite sequencing (WGBS). Profiling of the global DNA methylation revealed hypo-methylation, and 3471 significantly differentially methylated regions (DMRs) were observed in fsv1 ovules compared with Gui99. Based on functional annotation and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis of differentially methylated genes (DMGs), we observed more DMGs enriched in cellular component, reproduction regulation, metabolic pathway, and other pathways. In particular, many ovule development genes and plant hormone-related genes showed significantly different methylation patterns in the two rice lines, and these differences may provide important clues for revealing the mechanism of female gametophyte abortion.

Identification and expression analysis of microRNAs during ovule development in rice (Oryza sativa) by deep sequencing

MicroRNA (miRNA) expression profiles during rice ovule development revealed the possible miRNA-mediated regulation between ovule sporophytic tissue and female gametophyte and the involvement of miRNAs in programmed cell death. MiRNAs are 20-24-nucleotide small RNAs that play key roles in the regulation of many growth and developmental processes in plants. Rice ovule development comprises a series of biological events, which are regulated by complex molecular mechanisms. To gain insight into miRNA-mediated regulation of rice ovule development, Illumina sequencing was used to examine the expression of miRNAs from the megaspore mother cell meiosis stage to the fertilized ovule stage. Based on the sequencing data, 486 known and 204 novel miRNAs were identified during rice ovule development. Moreover, 56, 65 and 11 differentially expressed miRNAs between adjacent developmental stages were identified. By analyzing transcriptome and degradome data, we identified 41, 65 and 12 coherent target genes for the differentially expressed miRNAs in ovule development. We found that changes in the expression of plant hormone-related miRNAs may play important roles in embryo sac development, providing evidence for cross-talk communication between sporophytic tissue and the female gametophyte. Additionally, we revealed that miRNAs may be involved in programmed cell death after fertilization. Finally, we constructed miRNA-mediated regulatory networks that are active during rice ovule development.

Propofol regulates the expression of TLR4 through miR‑21 in human umbilical vein endothelial cells

Propofol (2,6-diisopropylphenol) is one of the most commonly used intravenous anesthetics. Anesthetics can regulate the inflammatory process; however, the mechanism remains to be fully elucidated. The present study aimed to investigate whether and how propofol affects the inflammatory reaction in human umbilical vein endothelial cells (HUVECs). The expression levels of Toll‑like receptor 4 (TLR4) and cluster of differentiation 14 (CD14) were determined in HUVECs treated with propofol and lipopolysaccharide (LPS) using western blot and reverse transcription‑quantitative polymerase chain reaction analyses. In addition, whether propofol regulated the expression of TLR4 though microRNA (miR)‑21 was examined. The results showed that LPS promoted the expression levels of TLR4, CD14 and tumor necrosis factor α (TNFα), and suppressed the expression of miR‑21 in HUVECs. Propofol suppressed the expression levels of TLR4, CD14 and TNFα, and upregulated the expression of miR‑21 in a concentration‑dependent manner. miR‑21 downregulated the expression of TLR4 at the mRNA and protein levels, whereas the miR‑21 mimic reversed the effect of LPS on the expression of TLR4. In addition, the miR‑21 inhibitor inhibited the downregulatory effect of propofol on the expression of TLR4. TargetScan analysis showed that TLR4 was included in the list of targets of miR‑21. Fluorescent reporter assays showed that the miR‑21 mimic and propofol treatment reduced the fluorescence intensity in cells transfected with a reporter vector containing the wild‑type TLR4 3'‑untranslated region. Taken together, the results of the present study demonstrated that propofol regulated the expression of TLR4 in HUVECs through miR‑21.

Genome wide identification of cotton (Gossypium hirsutum)-encoded microRNA targets against Cotton leaf curl Burewala virus

Cotton leaf curl Burewala virus (CLCuBV, genus Begomovirus) causes devastating cotton leaf curl disease. Among various known virus controlling strategies, RNAi-mediated one has shown potential to protect host crop plants. Micro(mi) RNAs, are the endogenous small RNAs and play a key role in plant development and stress resistance. In the present study we have identified cotton (Gossypium hirsutum)-encoded miRNAs targeting the CLCuBV. Based on threshold free energy and maximum complementarity scores of host miRNA-viral mRNA target pairs, a number of potential miRNAs were annotated. Among them, ghr-miR168 was selected as the most potent candidate, capable of targeting several vital genes namely C1, C3, C4, V1 and V2 of CLCuBV genome. In addition, ghr-miR395a and ghr-miR395d were observed to target the overlapping transcripts of C1 and C4 genes. We have verified the efficacy of these miRNA targets against CLCuBV following suppression of RNAi-mediated virus control through translational inhibition or cleavage of viral mRNA.

Genetic Subtraction Profiling Identifies Candidate miRNAs Involved in Rice Female Gametophyte Abortion

The female gametophyte is an important participant in the sexual reproduction of plants. The molecular mechanism of its development has received much attention in recent years. As important regulators of gene expression, miRNAs have been certified to play a significant role in many biological processes of plants, including sexual reproduction. In this study, to investigate the potential regulatory effects of miRNAs on rice female gametophyte abortion, we used the high-throughput sequencing method to compare the miRNA transcriptome in ovules of a high frequency female-sterile line (fsv1) and a rice wild-type line (Gui 99) during ovule development. As a result, 522 known miRNAs and 295 novel miRNAs were expressed in the developing ovule of rice, while 100 known miRNAs were significantly differentially expressed between these two rice lines during ovule development. Combining with gene expression information, a total of 627 coherent target genes of these differential expressed known miRNAs between fsv1 and Gui 99 were identified. The functional analyses of these coherent target genes revealed that the coherent target genes of differential expressed known miRNAs between the two rice lines are involved in many biological pathways, such as protein degradation, auxin signal transduction, and transcription factor regulation. These results provide us with important clues to investigate the regulatory roles of miRNAs in rice female gametophyte abortion.