MicroRNAs in Control of Plant Development

Comparative genome-wide analysis of circular RNAs in Brassica napus L.: target-site versus non-target-site resistance to herbicide stress

Circular RNAs (circRNAs), a class of non-coding RNA molecules, are recognized for their unique functions; however, their responses to herbicide stress in Brassica napus remain unclear. In this study, the role of circRNAs in response to herbicide treatment was investigated in two rapeseed cultivars: MH33, which confers non-target-site resistance (NTSR), and EM28, which exhibits target-site resistance (TSR). The genome-wide circRNA profiles of herbicide-stressed and non-stressed seedlings were analyzed. The findings indicate that NTSR seedlings exhibited a greater abundance of circRNAs, shorter lengths of circRNAs and their parent genes, and more diverse functions of parent genes compared with TSR seedlings. Compared to normal-growth plants, the herbicide-stressed group exhibited similar trends in the number of circRNAs, functions of parent genes, and differentially expressed circRNAs as observed in NTSR seedlings. In addition, a greater number of circRNAs that function as competing microRNA (miRNA) sponges were identified in the herbicide stress and NTSR groups compared to the normal-growth and TSR groups, respectively. The differentially expressed circRNAs were validated by qPCR. The differntially expressed circRNA-miRNA networks were predicted, and the mRNAs targeted by these miRNAs were annotated. Our results suggest that circRNAs play a crucial role in responding to herbicide stress, exhibiting distinct responses between NTSR and TSR in rapeseed. These findings offer valuable insights into the mechanisms underlying herbicide resistance in rapeseed.

Comparative analysis of miRNA expression profiles in flowering and non-flowering tissue of Crocus sativus L

Crocus sativus is a valuable plant due to the presence of apocarotenoids in its stigma. Considerable work has been done in the past to understand the apocarotenoid biosynthetic pathway in saffron. However, the reports on understanding the regulation of flowering at the post-transcriptional level are meagre. The study aimed to discover the candidate miRNAs, target genes, transcription factors (TFs), and apocarotenoid biosynthetic pathway genes associated with the regulation and transition of flowering in C. sativus. In the present investigation, miRNA profiling was performed in flowering and non-flowering corms of saffron, along with expression analysis of apocarotenoid genes and transcription factors involved in the synthesis of secondary metabolites. Significant modulation in the expression of miR156, miR159, miR166, miR172, miR395, miR396, miR399, and miR408 gene families was observed. We obtained 36 known miRNAs (26 in flowering and 10 in non-flowering) and 64 novel miRNAs (40 in flowering and 24 in non-flowering) unique to specific tissues in our analysis. TFs, including CsMADS and CsMYb, showed significant modulation in expression in flowering tissue, followed by CsHB. Additionally, the miRNAs were predicted to be involved in carbohydrate metabolism, phytohormone signalling, regulation of flower development, and response to stress, cold, and defence. The comprehensive study has enhanced our understanding of the regulatory machinery comprising factors like phytohormones, abiotic stress, apocarotenoid genes, transcription factors, and miRNAs responsible for the synthesis of apocarotenoids and developmental processes during and after flowering.

MicroRNA1432 regulates rice drought stress tolerance by targeting the CALMODULIN-LIKE2 gene

Due to climate change, drought has become a major threat to rice (Oryza sativa L.) growth and yield worldwide. Understanding the genetic basis of drought tolerance in rice is therefore of great importance. Here, we identified a microRNA, miR1432, which regulates rice drought tolerance by targeting the CALMODULIN-LIKE2 (OsCaML2) gene. Mutation of MIR1432 or suppression of miR1432 expression significantly impaired seed germination and seedling growth under drought-stress conditions. Molecular analysis demonstrated that miR1432 affected rice drought tolerance by directly targeting OsCaML2, which encodes an EF-hand chiral calcium-binding protein. Overexpression of a miR1432-resistant form of OsCaML2 (OEmCaML2) phenocopied the mir1432 mutant and miR1432 suppression plants. Furthermore, the suppression of miR1432 severely affected the expression of genes involved in responses to stimulation, metabolism and signal transduction, especially the mitogen-activated protein kinase (MAPK) pathway and hormone transduction pathway in rice under drought stress. Thus, our findings show that the miR1432-OsCaML2 module plays an important role in the regulation of rice drought tolerance, suggesting its potential utilization in developing molecular breeding strategies that improve crop drought tolerance.

MicroRNA Profiling Revealed the Mechanism of Enhanced Cold Resistance by Grafting in Melon (Cucumis melo L.)

Grafting is widely used to improve the resistance to abiotic stresses in cucurbit plants, but the effect and molecular mechanism of grafting on cold stress are still unknown in melon. In this study, phenotypic characteristics, physiological indexes, small-RNA sequencing and expression analyses were performed on grafted plants with pumpkin rootstock (PG) and self-grafted plants (SG) to explore the mechanism of changed cold tolerance by grafting in melon. Compared with SG plants, the cold tolerance was obviously enhanced, the malondialdehyde (MDA) content was significantly decreased and the activities of antioxidant enzymes (superoxide dismutase, SOD; catalase, CAT; peroxidase, POD) were significantly increased in PG plants. Depend on differentially expressed miRNA (DEM) identification and expression pattern analyses, cme-miR156b, cme-miR156f and chr07_30026 were thought to play a key role in enhancing low-temperature resistance resulting from grafting. Subsequently, 24, 37 and 17 target genes of cme-miR156b, cme-miR156f and chr07_30026 were respectively predicted, and 21 target genes were co-regulated by cme-miR156b and cme-miR156f. Among these 57 unique target genes, the putative promoter of 13 target genes contained the low-temperature responsive (LTR) cis-acting element. The results of qRT-PCR indicated that six target genes (MELO3C002370, MELO3C009217, MELO3C018972, MELO3C016713, MELO3C012858 and MELO3C000732) displayed the opposite expression pattern to their corresponding miRNAs. Furthermore, MELO3C002370, MELO3C016713 and MELO3C012858 were significantly downregulated in cold-resistant cultivars and upregulated in cold-sensitive varieties after cold stimulus, and they acted as the key negative regulators of low-temperature response in melon. This study revealed three key miRNAs and three putative target genes involved in the cold tolerance of melon and provided a molecular basis underlying how grafting improved the low-temperature resistance of melon plants.

MicroRNA397 promotes rice flowering by regulating the photorespiration pathway

The precise timing of flowering plays a pivotal role in ensuring successful plant reproduction and seed production. This process is intricately governed by complex genetic networks that integrate internal and external signals. This study delved into the regulatory function of microRNA397 (miR397) and its target gene LACCASE-15 (OsLAC15) in modulating flowering traits in rice (Oryza sativa). Overexpression of miR397 led to earlier heading dates, decreased number of leaves on the main stem, and accelerated differentiation of the spikelet meristem. Conversely, overexpression of OsLAC15 resulted in delayed flowering and prolonged vegetative growth. Through biochemical and physiological assays, we uncovered that miR397-OsLAC15 had a profound impact on carbohydrate accumulation and photosynthetic assimilation, consequently enhancing the photosynthetic intensity in miR397-overexpressing rice plants. Notably, we identified that OsLAC15 is at least partially localized within the peroxisome organelle, where it regulates the photorespiration pathway. Moreover, we observed that a high CO2 concentration could rescue the late flowering phenotype in OsLAC15-overexpressing plants. These findings shed valuable insights into the regulatory mechanisms of miR397-OsLAC15 in rice flowering and provided potential strategies for developing crop varieties with early flowering and high-yield traits through genetic breeding.

FvemiR160-FveARF18A-FveAP1/FveFUL module regulates flowering time in woodland strawberry

Flowering is an indicator of plant transformation from vegetative to reproductive growth. miR160 has been shown to have a significant effect on the growth and development of fruits, leaves, and roots of plants or their stress response to environment, but the participation of miR160 in regulating flowering time in plants is unclear. In this study, we found that two FvemiR160s (FvemiR160a/FvemiR160b) mature sequences in strawberry (Fragaria vesca) were consistent. It was displayed that the miR160 mature sequence is highly conserved in various species, and the miR160 mature sequence formed by the 5' arm of the MIR160 precursor was more conserved. Three FveARFs in woodland strawberry were negatively regulated by FvemiR160a, among which FveARF18A was the most significant. Phylogenetic analysis indicated that FvemiR160 is closely related to apple (Malus domestica), grape (Vitis vinifera), and Arabidopsis thaliana, while FveARF18A is closely related to RcARF18. Subsequently, we demonstrated that FvemiR160a can target cutting FveARF18A to negatively regulate its expression by RLM-5' RACE, cleavage site mutation, and GFP fluorescence assay. Moreover, we observed that FveMIR160a overexpressed plants have advanced flowering, while mFveARF18A overexpressed plants have delayed flowering. We also verified that FveARF18A negatively regulates the expression of FveAP1 and FveFUL by binding their promoters by yeast one-hybrid, LUC, and GUS assay, and FveAP1 and FveFUL transgenic Arabidopsis showed early flowering phenotype. In addition, the expression level of FvemiR160a was decreased obviously while that of FveARF18A was increased obviously by MeJA, GA and IAA. In conclusion, our study reveals the important role of the FvemiR160-FveARF18A-FveAP1/FveFUL module in the flowering process of woodland strawberry and provides a new pathway for studying flowering.

General Analysis of Heat Shock Factors in the Cymbidium ensifolium Genome Provided Insights into Their Evolution and Special Roles with Response to Temperature

Heat shock factors (HSFs) are the key regulators of heat stress responses and play pivotal roles in tissue development and the temperature-induced regulation of secondary metabolites. In order to elucidate the roles of HSFs in Cymbidium ensifolium, we conducted a genome-wide identification of CeHSF genes and predicted their functions based on their structural features and splicing patterns. Our results revealed 22 HSF family members, with each gene containing more than one intron. According to phylogenetic analysis, 59.1% of HSFs were grouped into the A subfamily, while subfamily HSFC contained only two HSFs. And the HSF gene families were differentiated evolutionarily between plant species. Two tandem repeats were found on Chr02, and two segmental duplication pairs were observed on Chr12, Chr17, and Chr19; this provided evidence for whole-genome duplication (WGD) events in C. ensifolium. The core region of the promoter in most CeHSF genes contained cis-acting elements such as AP2/ERF and bHLH, which were associated with plant growth, development, and stress responses. Except for CeHSF11, 14, and 19, each of the remaining CeHSFs contained at least one miRNA binding site. This included binding sites for miR156, miR393, and miR319, which were responsive to temperature and other stresses. The HSF gene family exhibited significant tissue specificity in both vegetative and floral organs of C. ensifolium. CeHSF13 and CeHSF15 showed relatively significant expression in flowers compared to other genes. During flower development, CeHSF15 exhibited markedly elevated expression in the early stages of flower opening, implicating critical regulatory functions in organ development and floral scent-related regulations. During the poikilothermic treatment, CeHSF14 was upregulated over 200-fold after 6 h of heat treatment. CeHSF13 and CeHSF14 showed the highest expression at 6 h of low temperature, while the expression of CeHSF15 and CeHSF21 continuously decreased at a low temperature. The expression patterns of CeHSFs further confirmed their role in responding to temperature stress. Our study may help reveal the important roles of HSFs in plant development and metabolic regulation and show insight for the further molecular design breeding of C. ensifolium.

The update and transport of aluminum nanoparticles in plants and their biochemical and molecular phototoxicity on plant growth and development: A systematic review

As aluminum nanoparticles (Al-NPs) are widely used in our daily life and various industries, Al-NPs has been becoming an emerging pollution in the environment. The impact of this NP has been attracting more and more attention from the scientific communities. In this review, we systematically summarized the interactions, uptake, and transport of Al-NPs in the plant system. Al-NPs can enter plants through different pathways and accumulate in various tissues, leading to alter plant growth and development. Al-NPs also affected root, shoot, and leaf characteristics as well as changing nutrient uptake and distribution and inducing oxidative stress via excess reactive radical generation, thereby impairing plant defense systems. Additionally, Al-NPs altered gene expression, which involved in various signaling pathways and metabolic processes in plants, that further altered plants susceptible or tolerant to stressors. The review also emphasized the effects of Al-NP size, surface charge, concentration, and exposure duration on plant growth and development. In the future, more research should be focused on mechanisms underlying Al-NPs phytotoxicity and potential risk to humans and off-target species.

High CO2 adaptation mechanisms revealed in the miR156-regulated flowering time pathway

Elevated CO2 concentrations have been observed to accelerate flowering time in Arabidopsis through the action of a highly conserved regulatory network controlled by miR156 and miR172. However, the network's robustness to the impact of increasing CO2 concentrations on flowering time remains poorly understood. In this study, we investigate this question by conducting a comprehensive analysis of the global landscape of network dynamics, including quantifying the probabilities associated with juvenile and flowering states and assessing the speed of the transition between them. Our findings reveal that a CO2 concentration range of 400-800ppm only mildly advances flowering time, contrasting with the dramatic changes from 200 to 300ppm. Notably, the feedback regulation of miR156 by squamosal promoter binding protein-like proteins (SPLs) plays a substantial role in mitigating the effects of increasing CO2 on flowering time. Intriguingly, we consistently observe a correlation between delayed flowering time and increased variance in flowering time, and vice versa, suggesting that this might be an intrinsic adaptation mechanism embedded within the network. To gain a deeper understanding of this network's dynamics, we identified the sensitive features within the feedback loops of miR156 SPLs and miR172-APETALA2 family proteins (AP2s), with the latter proving to be the most sensitive. Strikingly, our study underscores the indispensability of all feedback regulations in maintaining both juvenile and adult states as well as the transition time between them. Together, our research provides the first physical basis in plant species, aiding in the elucidation of novel regulatory mechanisms and the robustness of the miRNAs-regulated network in response to increasing CO2, therefore influencing the control of flowering time. Moreover, this study provides a promising strategy for engineering plant flowering time to enhance their adaptation and resilience.

Heat-responsive microRNAs participate in regulating the pollen fertility stability of CMS-D2 restorer line under high-temperature stress

Anther development and pollen fertility of cytoplasmic male sterility (CMS) conditioned by Gossypium harknessii cytoplasm (CMS-D2) restorer lines are susceptible to continuous high-temperature (HT) stress in summer, which seriously hinders the large-scale application of "three-line" hybrids in production. Here, integrated small RNA, transcriptome, degradome, and hormone profiling was performed to explore the roles of microRNAs (miRNAs) in regulating fertility stability in mature pollens of isonuclear alloplasmic near-isogenic restorer lines NH and SH under HT stress at two environments. A total of 211 known and 248 novel miRNAs were identified, of which 159 were differentially expressed miRNAs (DEMs). Additionally, 45 DEMs in 39 miRNA clusters (PmCs) were also identified, and most highly expressed miRNAs were significantly induced in SH under extreme HT, especially four MIR482 and six MIR6300 family miRNAs. PmC28 was located in the fine-mapped interval of the Rf1 gene and contained two DEMs, gra-miR482_L-2R + 2 and gma-miR2118a-3p_R + 1_1ss18TG. Transcriptome sequencing identified 6281 differentially expressed genes, of which heat shock protein (HSP)-related genes, such as HSP70, HSP22, HSP18.5-C, HSP18.2 and HSP17.3-B, presented significantly reduced expression levels in SH under HT stress. Through integrating multi-omics data, we constructed a comprehensive molecular network of miRNA-mRNA-gene-KEGG containing 35 pairs of miRNA/target genes involved in regulating the pollen development in response to HT, among which the mtr-miR167a_R + 1, tcc-miR167c and ghr-miR390a, tcc-miR396c_L-1 and ghr-MIR169b-p3_1ss6AG regulated the pollen fertility by influencing ARF8 responsible for the auxin signal transduction, ascorbate and aldarate metabolism, and the sugar and lipid metabolism and transport pathways, respectively. Further combination with hormone analysis revealed that HT-induced jasmonic acid signaling could activate the expression of downstream auxin synthesis-related genes and cause excessive auxin accumulation, followed by a cascade of auxin signal transduction, ultimately resulting in pollen abortion. The results provide a new understanding of how heat-responsive miRNAs regulate the stability of fertility restoration for CMS-D2 cotton under heat stress.

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.

Acclimation response and management strategies to combat heat stress in wheat for sustainable agriculture: A state-of-the-art review

Unpredicted variability in climate change on the planet is associated with frequent extreme high-temperature events impacting crop yield globally. Wheat is an economically and nutritionally important crop that fulfils global food requirements and each degree rise in temperature results in ∼6% of its yield reduction. Thus, understanding the impact of climate change, especially the terminal heat stress on global wheat production, becomes critically important for policymakers, crop breeders, researchers and scientists to ensure global food security. This review describes how wheat perceives heat stress and induces stress adaptation events by its morpho-physiological, phenological, molecular, and biochemical makeup. Temperature above a threshold level in crop vicinity leads to irreversible injuries, viz. destruction of cellular membranes and enzymes, generation of active oxygen species, redox imbalance, etc. To cope with these changes, wheat activates its heat tolerance mechanisms characterized by hoarding up soluble carbohydrates, signalling molecules, and heat tolerance gene expressions. Being vulnerable to heat stress, increasing wheat production without delay seeks strategies to mitigate the detrimental effects and provoke the methods for its sustainable development. Thus, to ensure the crop's resilience to stress and increasing food demand, this article circumscribes the integrated management approaches to enhance wheat's performance and adaptive capacity besides its alleviating risks of increasing temperature anticipated with climate change. Implementing these integrated strategies in the face of risks from rising temperatures will assist us in producing sustainable wheat with improved yield.

Profiling of conserved orthologs and miRNAs for understanding their role in salt tolerance mechanism of Sesuvium portulacastrum L

BACKGROUND

Sesuvium portulacastrum is a facultative halophyte capable of thriving in a saline environment. Despite molecular studies conducted to unravel its salt adaptation mechanism, there is a paucity of information on the role of salt-responsive orthologs and microRNAs (miRNAs) in this halophyte. Here, we searched the orthology to identify salt-responsive orthologs and miRNA targets of Sesuvium using the Arabidopsis genome.

METHODS

The relative fold change of orthologs, conserved miRNAs, and miRNA targets of Sesuvium was analyzed under 100 mM (LS) and 250 mM NaCl (HS) treatment at 24 h using qRT-PCR. The comparison between the expression of Sesuvium orthologs and Arabidopsis orthologs (Arabidopsis eFP browser database) was used to identify differentially expressed genes.

RESULTS

Upon salt treatment, we found that SpCIPK3 (1.95-fold in LS and 2.90-fold in HS) in Sesuvium roots, and SpNHX7 (1.61-fold in LS and 6.39-fold in HS) and, SpSTPK2 (2.54-fold in LS and 7.65-fold in HS) in Sesuvium leaves were upregulated in a salt concentration-specific manner. In Arabidopsis, these genes were either downregulated or did not show significant variation, implicating its significance in the halophytic nature of Sesuvium. Furthermore, miRNAs like miR394a, miR396a, and miR397a exhibited a negative correlation with their targets-Frigida interacting protein 1, Cysteine proteinases superfamily protein, and Putative laccase, respectively under different salt treatments.

CONCLUSION

The study revealed that the high salt tolerance in Sesuvium is associated with distinct transcriptional reprogramming, hence, to gain holistic mechanistic insights, global-scale profiling is required.

MicroRNA-1 targets ribosomal protein genes to regulate the growth, development and reproduction of Schistosoma japonicum

Eggs laid by mature female schistosomes are primarily responsible for the pathogenesis of schistosomiasis and critical for transmission. Consequently, elucidating the mechanism of sexual maturation as well as egg production may lead to new strategies for the control of schistosomiasis. MicroRNAs (miRNAs) are involved in multiple biological processes including reproduction in many organisms, yet their roles have not been well characterized in schistosomes. Here, we investigated microRNA-1 (miR-1), which was downregulated gradually in both male and female Schistosoma japonicum after they reached sexually maturity. The expression of miR-1, as shown with quantitative reverse transcription PCR (qRT-PCR), was lower in the reproductive organs of adult females compared with the somatic tissues. Overexpression of miR-1 in adult worms destroyed the morphological architecture of reproductive organs and reduced the subsequent oviposition, which may be due to the activation of apoptosis pathways. Through in silico analysis, 34 potential target genes of miR-1 were identified, including five ribosomal protein genes, called rp-s13, rp-l7ae, rp-l14, rp-l11 and rp-s24e. In vitro dual-luciferase reporter gene assays and miRNA overexpression experiments further validated that these ribosomal protein genes were directly regulated by miR-1. In contrast to the gene expression of miR-1, qRT-PCR and in situ hybridization experiments demonstrated these ribosomal protein genes were enriched in the sexual organs of adult females. Using RNA interference to silence the ribosomal protein genes in different developmental stages in a mouse model system, we demonstrated that these miR-1 target genes not only participated in the reproductive development of S. japonicum, but also were required for the growth and survival of the parasite in the early developmental stages. Taken together, our data suggested that miR-1 may affect the growth, reproduction and oviposition of S. japonicum by targeting the ribosomal protein genes, which provides insights for exploration of new anti-schistosome strategies.

Genome-wide identification of the longan R2R3-MYB gene family and its role in primary and lateral root

R2R3-MYB is an important transcription factor family that regulates plant growth and development. Root development directly affects the absorption of water and nutrients by plants. Therefore, to understand the regulatory role of R2R3-MYB transcription factor family in root development of longan, this study identified the R2R3-MYB gene family members at the genome-wide level, and analyzed their phylogenetic characteristics, physical and chemical properties, gene structure, chromosome location and tissue expression. The analysis identified 124 R2R3-MYB family members in the longan genome. Phylogenetic analysis divided these members into 22 subfamilies, and the members of the unified subfamily had similar motifs and gene structures. The result of qRT-PCR showed that expression levels of DlMYB33, DlMYB34, DlMYB59, and DlMYB77 were significantly higher in main roots than in lateral as opposed to those of DlMYB35, DlMYB69, DlMYB70, and DlMYB83, which were significantly lower. SapBase database prediction and miRNAs sequencing results showed that 34 longan miRNAs could cleave R2R3-MYB, including 17 novel miRNAs unique to longan. The qRT-PCR and subcellular localization experiments of DlMYB92 and DlMYB98 showed that DlMYB92 is a key factor that regulates transcription in the nucleus and participates in the regulation of longan lateral root development. Longan also has a conserved miRNA-MYB-lateral root development regulation mechanism. This study provides a reference for further research on the transcriptional regulation of the miRNA-R2R3-MYB module in the root development of longan.

Cold-responsive transcription factors in Arabidopsis and rice: A regulatory network analysis using array data and gene co-expression network

Plant growth and development can be influenced by cold stress. Responses of plants to cold are regulated in part by transcription factors (TFs) and microRNAs, which their determination would be necessary in comprehension of the corresponding molecular cues. Here, transcriptomes of Arabidopsis and rice were analyzed to computationally determine TFs and microRNAs that are differentially responsive to cold treatment, and their co-expression networks were established. Among 181 Arabidopsis and 168 rice differentially expressed TF genes, 37 (26 novel) were up- and 16 (8 novel) were downregulated. Common TF encoding genes were from ERF, MYB, bHLH, NFY, bZIP, GATA, HSF and WRKY families. NFY A4/C2/A10 were the significant hub TFs in both plants. Phytohormone responsive cis-elements such as ABRE, TGA, TCA and LTR were the common cis-elements in TF promoters. Arabidopsis had more responsive TFs compared to rice possibly due to its greater adaptation to ranges geographical latitudes. Rice had more relevant miRNAs probably because of its bigger genome size. The interacting partners and co-expressed genes were different for the common TFs so that of the downstream regulatory networks and the corresponding metabolic pathways. Identified cold-responsive TFs in (A + R) seemed to be more engaged in energy metabolism esp. photosynthesis, and signal transduction, respectively. At post-transcriptional level, miR5075 showed to target many identified TFs in rice. In comparison, the predictions showed that identified TFs are being targeted by diverse groups of miRNAs in Arabidopsis. Novel TFs, miRNAs and co-expressed genes were introduced as cold-responsive markers that can be harnessed in future studies and development of crop tolerant varieties.

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.

Regulation of the regulators: Transcription factors controlling biosynthesis of plant secondary metabolites during biotic stresses and their regulation by miRNAs

Biotic stresses threaten to destabilize global food security and cause major losses to crop yield worldwide. In response to pest and pathogen attacks, plants trigger many adaptive cellular, morphological, physiological, and metabolic changes. One of the crucial stress-induced adaptive responses is the synthesis and accumulation of plant secondary metabolites (PSMs). PSMs mitigate the adverse effects of stress by maintaining the normal physiological and metabolic functioning of the plants, thereby providing stress tolerance. This differential production of PSMs is tightly orchestrated by master regulatory elements, Transcription factors (TFs) express differentially or undergo transcriptional and translational modifications during stress conditions and influence the production of PSMs. Amongst others, microRNAs, a class of small, non-coding RNA molecules that regulate gene expression post-transcriptionally, also play a vital role in controlling the expression of many such TFs. The present review summarizes the role of stress-inducible TFs in synthesizing and accumulating secondary metabolites and also highlights how miRNAs fine-tune the differential expression of various stress-responsive transcription factors during biotic stress.

MicroRNA Expression Profiles in Response to Phytophthora infestans and Oidium neolycopersici and Functional Identification of sly-miR397 in Tomato

Late blight and powdery mildew are two widespread tomato diseases caused by Phytophthora infestans and Oidium neolycopersici, respectively, which reduce the quantity and quality of tomato. MicroRNAs (miRNAs) play critical roles in tomato resistance to various pathogens. Investigating the function of miRNAs is of great significance in controlling tomato diseases. To identify potential miRNAs involved in the interaction of tomato with P. infestans or O. neolycopersici, we analyzed the expression profiles of small RNAs in tomato leaves infected with these two pathogens using RNA-seq technology. A total of 330 and 288 miRNAs exhibited differences in expression levels after exposure to P. infestans and O. neolycopersici, respectively. One hundred and forty-six commonly differentially expressed (DE) miRNAs responsive to P. infestans and O. neolycopersici infestation were detected, including 10 commonly known conserved DE miRNAs and 136 novel miRNAs. Among these known DE miRNAs, sly-miR397 was strongly downregulated in response to P. infestans or O. neolycopersici infection. Silencing of sly-miR397 resulted in enhanced tolerance to the pathogens, whereas overexpression of sly-miR397 showed increased susceptibility. Furthermore, changes in sly-miR397 expression could also affect expression levels of pathogenesis-related genes and reactive oxygen species-scavenging genes, leading to altered necrotic cells and H₂O₂ levels. In addition, the number of lateral branches significantly changed in transgenic plants. Taken together, our results provide potential miRNA resources for further research of miRNA-disease associations and indicates that sly-miR397 acts as a negative regulator of disease resistance and influences lateral branch development in tomato.

Mdm-miR858 targets MdMYB9 and MdMYBPA1 to participate anthocyanin biosynthesis in red-fleshed apple

Anthocyanins are important secondary metabolites in plants. They are important for human health because of their antioxidant activities and because their dietary intake reduces the incidence of cardiovascular and cerebrovascular diseases and tumors. The biosynthesis of anthocyanins and its regulation in fruits and vegetables is a global research hotspot. Compared with cultivated apples, the red-fleshed apple is a relatively new and popular commodity in the market. Previous studies on red-fleshed apples have focused on the basis for the high anthocyanin content and the transcriptional regulation of anthocyanin synthesis. In the present study, we focused on the mechanism of microRNA-mediated post-transcriptional regulation of anthocyanin synthesis in red-fleshed apples. We identified a microRNA (miRNA), designated mdm-miR858, that is specifically expressed in the flesh of apple fruit. The expression level of miR858 was significantly lower in red-fleshed apples than in white-fleshed apples. The overexpression of mdm-miR858 significantly inhibited anthocyanin accumulation, whereas the silencing of mdm-miR858 promoted anthocyanin synthesis in STTM858 transgenic apple calli. Further analyses showed that mdm-miR858 targets the transcription factor genes MdMYB9 and MdMYBPA1 to participate anthocyanin accumulation in apple. Our results also show that MdHY5, a transcription factor in the light signaling pathway, can bind to the promoter of mdm-miR858 to inhibit its transcription, thereby regulating anthocyanin synthesis. Based on our results, we describe a novel HY5-miR858-MYB loop involved in the modulation of anthocyanin biosynthesis. These findings provide new information about how plant miRNAs regulate anthocyanin anabolism and provide a basis for breeding new anthocyanin-rich, red-fleshed apple varieties.

Tweaking the Small Non-Coding RNAs to Improve Desirable Traits in Plant

Plant transcriptome contains an enormous amount of non-coding RNAs (ncRNAs) that do not code for proteins but take part in regulating gene expression. Since their discovery in the early 1990s, much research has been conducted to elucidate their function in the gene regulatory network and their involvement in plants' response to biotic/abiotic stresses. Typically, 20-30 nucleotide-long small ncRNAs are a potential target for plant molecular breeders because of their agricultural importance. This review summarizes the current understanding of three major classes of small ncRNAs: short-interfering RNAs (siRNAs), microRNA (miRNA), and transacting siRNAs (tasiRNAs). Furthermore, their biogenesis, mode of action, and how they have been utilized to improve crop productivity and disease resistance are discussed here.

Involvement of miRNAs in Metabolic Herbicide Resistance to Bispyribac-Sodium in Echinochloa crus-galli (L.) P. Beauv

Several mechanisms involved in weed herbicide resistance are unknown, particularly those acting at the epigenetic level, such as the capacity of small-non-coding RNAs (sncRNAs) to target messenger RNAs of genes involved in herbicide detoxification. The transcription of these sncRNAs is stimulated by epigenetic factors, thereby affecting gene expression. This study was carried out in order to evaluate, for the first time in Echinochloa crus-galli (L.) P. Beauv. (barnyardgrass), the capacity of miRNAs to regulate the expression of genes associated with bispyribac-sodium detoxification. The expression profiles of eight miRNAs with a high degree of complementarity (≥80%) with mRNAs of genes involved in herbicide detoxification (CYP450, GST and eIF4B) were determined by qRT-PCR before and after herbicide spraying. Five of the miRNAs studied (gra-miR7487c, gma-miR396f, gra-miR8759, osa-miR395f, ath-miR847) showed an increased expression after herbicide application in both susceptible and resistant biotypes. All the miRNAs, except gra-miR8759, were more highly expressed in the herbicide-resistant biotypes. In specimens with increased expression of miRNAs, we observed reduced expression of the target genes. The remaining three miRNAs (ata-miR166c-5p, ath-miR396b-5p and osa-miR5538) showed no over-expression after herbicide treatment, and no difference in expression was recorded between susceptible and resistant biotypes. Our results represent a first overview of the capacity of miRNAs to regulate the expression of genes involved in bispyribac-sodium detoxification in the genus Echinochloa. Further research is required to identify novel miRNAs and target genes to develop more focused and sustainable strategies of weed control.

MdZAT5 regulates drought tolerance via mediating accumulation of drought-responsive miRNAs and mRNAs in apple

Drought limits apple yield and fruit quality. However, the molecular mechanism of apple in response to drought is not well known. Here, we report a Cys2/His2 (C2H2)-type zinc-finger protein, MdZAT5, that positively regulates apple drought tolerance by regulating drought-responsive RNAs and microRNAs (miRNAs). DNA affinity purification and sequencing and yeast-one hybrid analysis identified the binding motifs of MdZAT5, T/ACACT/AC/A/G. Chromatin immunoprecipitation quantitative polymerase chain reaction (ChIP-qPCR) and electrophoretic mobility shift assays (EMSAs) showed that MdZAT5 directly binds to the promoters of the drought-responsive genes including MdRHA2a, MdLEA14, MdTPX1, and MdCAT3, and activates their expression under drought stress. MdZAT5 interacts with and directly targets HYPONASTIC LEAVES1 (MdHYL1). MdZAT5 may facilitate the interaction of MdHYL1 with pri-miRNAs or MdDCL1 by activating MdHYL1 expression, thereby regulating the biogenesis of drought-responsive miRNAs. Genetic dissection showed that MdHYL1 is essential for MdZAT5-mediated drought tolerance and miRNA biogenesis. In addition, ChIP-qPCR and EMSA revealed that MdZAT5 binds directly to the promoters of some MIR genes including Mdm-miR171i and Mdm-miR172c, and modulates their transcription. Taken together, our findings improve our understanding of the molecular mechanisms of drought response in apple and provide a candidate gene for the breeding of drought-tolerant cultivars.

Temperature-sensitive male sterility in rice determined by the roles of AGO1d in reproductive phasiRNA biogenesis and function

Phased secondary siRNAs (phasiRNAs) are broadly present in the reproductive tissues of flowering plants, with spatial-temporal specificity. However, the ARGONAUTE (AGO) proteins associated with phasiRNAs and their miRNA triggers remain elusive. Here, through histological and high-throughput sequencing analyses, we show that rice AGO1d, which is specifically expressed in anther wall cells before and during meiosis, associates with both miR2118 and miR2275 to mediate phasiRNA biogenesis. AGO1d preferentially binds to miR2118-triggered 21-nucleotide (nt) phasiRNAs with a 5'-terminal uridine, suggesting a dual role in phasiRNA biogenesis and function. Depletion of AGO1d causes a reduction of 21- and 24-nt phasiRNAs and temperature-sensitive male sterility. At lower temperatures, anthers of the ago1d mutant predominantly show excessive tapetal cells with little starch accumulation during pollen formation, possibly caused by the dysregulation of cell metabolism. These results uncover an essential role of AGO1d in rice anther development at lower temperatures and demonstrate coordinative roles of AGO proteins during reproductive phasiRNA biogenesis and function.

Roles of non-coding RNAs in the hormonal and nutritional regulation in nodulation and nitrogen fixation

Symbiotic nitrogen fixation is an important component in the nitrogen cycle and is a potential solution for sustainable agriculture. It is the result of the interactions between the plant host, mostly restricted to legume species, and the rhizobial symbiont. From the first encounter between the host and the symbiont to eventual successful nitrogen fixation, there are delicate processes involved, such as nodule organogenesis, rhizobial infection thread progression, differentiation of the bacteroid, deregulation of the host defense systems, and reallocation of resources. All these processes are tightly regulated at different levels. Recent evidence revealed that non-coding RNAs (ncRNAs), including microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs), participate in these processes by controlling the transcription and translation of effector genes. In general, ncRNAs are functional transcripts without translation potential and are important gene regulators. MiRNAs, negative gene regulators, bind to the target mRNAs and repress protein production by causing the cleavage of mRNA and translational silencing. LncRNAs affect the formation of chromosomal loops, DNA methylation, histone modification, and alternative splicing to modulate gene expression. Both lncRNAs and circRNAs could serve as target mimics of miRNA to inhibit miRNA functions. In this review, we summarized and discussed the current understanding of the roles of ncRNAs in legume nodulation and nitrogen fixation in the root nodule, mainly focusing on their regulation of hormone signal transduction, the autoregulation of nodulation (AON) pathway and nutrient homeostasis in nodules. Unraveling the mediation of legume nodulation by ncRNAs will give us new insights into designing higher-performance leguminous crops for sustainable agriculture.

Copper microRNAs modulate the formation of giant feeding cells induced by the root knot nematode Meloidogyne incognita in Arabidopsis thaliana

Root-knot nematodes (RKNs) are root endoparasites that induce the dedifferentiation of a few root cells and the reprogramming of their gene expression to generate giant hypermetabolic feeding cells. We identified two microRNA families, miR408 and miR398, as upregulated in Arabidopsis thaliana and Solanum lycopersicum roots infected by RKNs. In plants, the expression of these two conserved microRNA families is known to be activated by the SPL7 transcription factor in response to copper starvation. By combining functional approaches, we deciphered the network involving these microRNAs, their regulator and their targets. MIR408 expression was located within nematode-induced feeding cells like its regulator SPL7 and was regulated by copper. Moreover, infection assays with mir408 and spl7 knockout mutants or lines expressing targets rendered resistant to cleavage by miR398 demonstrated the essential role of the SPL7/MIR408/MIR398 module in the formation of giant feeding cells. Our findings reveal how perturbation of plant copper homeostasis, via the SPL7/MIR408/MIR398 module, modulates the development of nematode-induced feeding cells.

Parasitic plant small RNA analyses unveil parasite-specific signatures of microRNA retention, loss, and gain

Parasitism is a successful life strategy that has evolved independently in several families of vascular plants. The genera Cuscuta and Orobanche represent examples of the two profoundly different groups of parasites: one parasitizing host shoots and the other infecting host roots. In this study, we sequenced and described the overall repertoire of small RNAs from Cuscuta campestris and Orobanche aegyptiaca. We showed that C. campestris contains a number of novel microRNAs (miRNAs) in addition to a conspicuous retention of miRNAs that are typically lacking in other Solanales, while several typically conserved miRNAs seem to have become obsolete in the parasite. One new miRNA appears to be derived from a horizontal gene transfer event. The exploratory analysis of the miRNA population (exploratory due to the absence of a full genomic sequence for reference) from the root parasitic O. aegyptiaca also revealed a loss of a number of miRNAs compared to photosynthetic species from the same order. In summary, our study shows partly similar evolutionary signatures in the RNA silencing machinery in both parasites. Our data bear proof for the dynamism of this regulatory mechanism in parasitic plants.

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.

4D genetic networks reveal the genetic basis of metabolites and seed oil-related traits in 398 soybean RILs

Background

The yield and quality of soybean oil are determined by seed oil-related traits, and metabolites/lipids act as bridges between genes and traits. Although there are many studies on the mode of inheritance of metabolites or traits, studies on multi-dimensional genetic network (MDGN) are limited.

Results

In this study, six seed oil-related traits, 59 metabolites, and 107 lipids in 398 recombinant inbred lines, along with their candidate genes and miRNAs, were used to construct an MDGN in soybean. Around 175 quantitative trait loci (QTLs), 36 QTL-by-environment interactions, and 302 metabolic QTL clusters, 70 and 181 candidate genes, including 46 and 70 known homologs, were previously reported to be associated with the traits and metabolites, respectively. Gene regulatory networks were constructed using co-expression, protein–protein interaction, and transcription factor binding site and miRNA target predictions between candidate genes and 26 key miRNAs. Using modern statistical methods, 463 metabolite–lipid, 62 trait–metabolite, and 89 trait–lipid associations were found to be significant. Integrating these associations into the above networks, an MDGN was constructed, and 128 sub-networks were extracted. Among these sub-networks, the gene–trait or gene–metabolite relationships in 38 sub-networks were in agreement with previous studies, e.g., oleic acid (trait)–GmSEI–GmDGAT1a–triacylglycerol (16:0/18:2/18:3), gene and metabolite in each of 64 sub-networks were predicted to be in the same pathway, e.g., oleic acid (trait)–GmPHS–d-glucose, and others were new, e.g., triacylglycerol (16:0/18:1/18:2)–GmbZIP123–GmHD-ZIPIII-10–miR166s–oil content.

Conclusions

This study showed the advantages of MGDN in dissecting the genetic relationships between complex traits and metabolites. Using sub-networks in MGDN, 3D genetic sub-networks including pyruvate/threonine/citric acid revealed genetic relationships between carbohydrates, oil, and protein content, and 4D genetic sub-networks including PLDs revealed the relationships between oil-related traits and phospholipid metabolism likely influenced by the environment. This study will be helpful in soybean quality improvement and molecular biological research.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13068-022-02191-1.

Expression Profiling and MicroRNA Regulatory Networks of Homeobox Family Genes in Sugarcane Saccharum spontaneum L

Homeobox (HB) genes play important roles in plant growth and development processes, particularly in the formation of lateral organs. Thus, they could influence leaf morphogenesis and biomass formation in plants. However, little is known about HBs in sugarcane, a crucial sugar crop, due to its complex genetic background. Here, 302 allelic sequences for 104 HBs were identified and divided into 13 subfamilies in sugarcane Saccharum spontaneum. Comparative genomics revealed that whole-genome duplication (WGD)/segmental duplication significantly promoted the expansion of the HB family in S. spontaneum, with SsHB26, SsHB63, SsHB64, SsHB65, SsHB67, SsHB95, and SsHB96 being retained from the evolutionary event before the divergence of dicots and monocots. Based on the analysis of transcriptome and degradome data, we speculated that SsHB15 and SsHB97 might play important roles in regulating sugarcane leaf morphogenesis, with miR166 and SsAGO10 being involved in the regulation of SsHB15 expression. Moreover, subcellular localization and transcriptional activity detection assays demonstrated that these two genes, SsHB15 and SsHB97, were functional transcription factors. This study demonstrated the evolutionary relationship and potential functions of SsHB genes and will enable the further investigation of the functional characterization and the regulatory mechanisms of SsHBs.

A comparative study of microRNAs in different stages of Eimeria tenella

Apicomplexan parasites have divergent biogenesis machinery for small RNA generation. Analysis has shown that parasites in Plasmodium and Cryptosporidium as well as many species in Leishmania or Trypanosoma do not have a complete machinery in small RNA biogenesis. Recently, the miRNA-generating system of Toxoplasma has been identified as plant/fungal-like and its miRNAome has been elucidated. However, the microRNA (miRNA) expression profiles and their potential regulatory functions in different stages of Eimeria tenella remain largely unknown. In this study, we characterized the RNA silencing machinery of E. tenella and investigated the miRNA population distribution at different life stages by high-throughput sequencing. We characterized the expression of miRNAs in the unsporulated oocyst, sporulated oocyst and schizogony stages, obtaining a total of 392 miRNAs. We identified 58 differentially expressed miRNAs between USO (unsporulated oocysts) and SO (sporulated oocysts) that were significantly enriched for their potential target genes in the regulation of gene expression and chromatin binding, suggesting an epigenetic modulation of sporulating by these miRNAs. In comparing miRNA expression at endogenous and exogenous developmental stages, twenty-four miRNAs were identified differently expressed. Those were mainly associated with the regulation of genes with protein kinase activity, suggesting control of protein phosphorylation. This is the first study about the evolution of miRNA biogenesis system and miRNA control of gene expression in Eimeria species. Our data may lead to functional insights into of the regulation of gene expression during parasite life cycle in apicomplexan parasites.

Biogenesis, conservation, and function of miRNA in liverworts

MicroRNAs (miRNAs) are small non-coding endogenous RNA molecules, 18-24 nucleotides long, that control multiple gene regulatory pathways via post-transcriptional gene silencing in eukaryotes. To develop a comprehensive picture of the evolutionary history of miRNA biogenesis and action in land plants, studies on bryophyte representatives are needed. Here, we review current understanding of liverwort MIR gene structure, miRNA biogenesis, and function, focusing on the simple thalloid Pellia endiviifolia and the complex thalloid Marchantia polymorpha. We review what is known about conserved and non-conserved miRNAs, their targets, and the functional implications of miRNA action in M. polymorpha and P. endiviifolia. We note that most M. polymorpha miRNAs are encoded within protein-coding genes and provide data for 23 MIR gene structures recognized as independent transcriptional units. We identify M. polymorpha genes involved in miRNA biogenesis that are homologous to those identified in higher plants, including those encoding core microprocessor components and other auxiliary and regulatory proteins that influence the stability, folding, and processing of pri-miRNAs. We analyzed miRNA biogenesis proteins and found similar domain architecture in most cases. Our data support the hypothesis that almost all miRNA biogenesis factors in higher plants are also present in liverworts, suggesting that they emerged early during land plant evolution.

Moderate Soil Drying-Induced Alternative Splicing Provides a Potential Novel Approach for the Regulation of Grain Filling in Rice Inferior Spikelets

Poor grain filling of inferior spikelets, especially in some large-panicle rice varieties, is becoming a major limitation in breaking the ceiling of rice production. In our previous studies, we proved that post-anthesis moderate soil drying (MD) was an effective way to promote starch synthesis and inferior grain filling. As one of the most important regulatory processes in response to environmental cues and at different developmental stages, the function of alternative splicing (AS) has not yet been revealed in regulating grain filling under MD conditions. In this study, AS events at the most active grain-filling stage were identified in inferior spikelets under well-watered control (CK) and MD treatments. Of 16,089 AS events, 1840 AS events involving 1392 genes occurred differentially between the CK and MD treatments, many of which function on spliceosome, ncRNA metabolic process, starch, and sucrose metabolism, and other functions. Some of the splicing factors and starch synthesis-related genes, such as SR protein, hnRNP protein, OsAGPL2, OsAPS2, OsSSIVa, OsSSIVb, OsGBSSII, and OsISA1 showed differential AS changes under MD treatment. The expression of miR439f and miR444b was reduced due to an AS event which occurred in the intron where miRNAs were located in the MD-treated inferior spikelets. On the contrary, OsAGPL2, an AGPase encoding gene, was alternatively spliced, resulting in different transcripts with or without the miR393b binding site, suggesting a potential mechanism for miRNA-mediated gene regulation on grain filling of inferior spikelets in response to MD treatment. This study provides some new insights into the function of AS on the MD-promoted grain filling of inferior spikelets, and potential application in agriculture to increase rice yields by genetic approaches.

Differential Responses of Wheat (Triticum aestivum L.) and Cotton (Gossypium hirsutum L.) to Nitrogen Deficiency in the Root Morpho-Physiological Characteristics and Potential MicroRNA-Mediated Mechanisms

Understanding the mechanism of crop response to nitrogen (N) deficiency is very important for developing sustainable agriculture. In addition, it is unclear if the microRNA-mediated mechanism related to root growth complies with a common mechanism in monocots and dicots under N deficiency. Therefore, the root morpho-physiological characteristics and microRNA-mediated mechanisms were studied under N deficiency in wheat (Triticum aestivum L.) and cotton (Gossypium hirsutum L.). For both crops, shoot dry weight, plant dry weight and total leaf area as well as some physiological traits, i.e., the oxygen consuming rate in leaf and root, the performance index based on light energy absorption were significantly decreased after 8 days of N deficiency. Although N deficiency did not significantly impact the root biomass, an obvious change on the root morphological traits was observed in both wheat and cotton. After 8 days of treatment with N deficiency, the total root length, root surface area, root volume of both crops showed an opposite trend with significantly decreasing in wheat but significantly increasing in cotton, while the lateral root density was significantly increased in wheat but significantly decreased in cotton. At the same time, the seminal root length in wheat and the primary root length in cotton were increased after 8 days of N deficiency treatment. Additionally, the two crops had different root regulatory mechanisms of microRNAs (miRNAs) to N deficiency. In wheat, the expressions of miR167, miR319, miR390, miR827, miR847, and miR165/166 were induced by N treatment; these miRNAs inhibited the total root growth but promoted the seminal roots growth and lateral root formation to tolerate N deficiency. In cotton, the expressions of miR156, miR167, miR171, miR172, miR390, miR396 were induced and the expressions of miR162 and miR393 were inhibited; which contributed to increasing in the total root length and primary root growth and to decreasing in the lateral root formation to adapt the N deficiency. In conclusion, N deficiency significantly affected the morpho-physiological characteristics of roots that were regulated by miRNAs, but the miRNA-mediated mechanisms were different in wheat and cotton.

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.

Boundary domain genes were recruited to suppress bract growth and promote branching in maize

Grass inflorescence development is diverse and complex and involves sophisticated but poorly understood interactions of genes regulating branch determinacy and leaf growth. Here, we use a combination of transcript profiling and genetic and phylogenetic analyses to investigate tasselsheath1 (tsh1) and tsh4, two maize genes that simultaneously suppress inflorescence leaf growth and promote branching. We identify a regulatory network of inflorescence leaf suppression that involves the phase change gene tsh4 upstream of tsh1 and the ligule identity gene liguleless2 (lg2). We also find that a series of duplications in the tsh1 gene lineage facilitated its shift from boundary domain in nongrasses to suppressed inflorescence leaves of grasses. Collectively, these results suggest that the boundary domain genes tsh1 and lg2 were recruited to inflorescence leaves where they suppress growth and regulate a nonautonomous signaling center that promotes inflorescence branching, an important component of yield in cereal grasses.

Regulation of Phenolic Compound Production by Light Varying in Spectral Quality and Total Irradiance

Photosynthetically active radiation (PAR) is an important environmental cue inducing the production of many secondary metabolites involved in plant oxidative stress avoidance and tolerance. To examine the complex role of PAR irradiance and specific spectral components on the accumulation of phenolic compounds (PheCs), we acclimated spring barley (Hordeum vulgare) to different spectral qualities (white, blue, green, red) at three irradiances (100, 200, 400 µmol m⁻² s⁻¹). We confirmed that blue light irradiance is essential for the accumulation of PheCs in secondary barley leaves (in UV-lacking conditions), which underpins the importance of photoreceptor signals (especially cryptochrome). Increasing blue light irradiance most effectively induced the accumulation of B-dihydroxylated flavonoids, probably due to the significantly enhanced expression of the F3′H gene. These changes in PheC metabolism led to a steeper increase in antioxidant activity than epidermal UV-A shielding in leaf extracts containing PheCs. In addition, we examined the possible role of miRNAs in the complex regulation of gene expression related to PheC biosynthesis.

Regulatory role of microRNAs (miRNAs) in the recent development of abiotic stress tolerance of plants

MicroRNAs (miRNAs) are a distinct groups of single-stranded non-coding, tiny regulatory RNAs approximately 20-24 nucleotides in length. miRNAs negatively influence gene expression at the post-transcriptional level and have evolved considerably in the development of abiotic stress tolerance in a number of model plants and economically important crop species. The present review aims to deliver the information on miRNA-mediated regulation of the expression of major genes or Transcription Factors (TFs), as well as genetic and regulatory pathways. Also, the information on adaptive mechanisms involved in plant abiotic stress responses, prediction, and validation of targets, computational tools, and databases available for plant miRNAs, specifically focus on their exploration for engineering abiotic stress tolerance in plants. The regulatory function of miRNAs in plant growth, development, and abiotic stresses consider in this review, which uses high-throughput sequencing (HTS) technologies to generate large-scale libraries of small RNAs (sRNAs) for conventional screening of known and novel abiotic stress-responsive miRNAs adds complexity to regulatory networks in plants. The discoveries of miRNA-mediated tolerance to multiple abiotic stresses, including salinity, drought, cold, heat stress, nutritional deficiency, UV-radiation, oxidative stress, hypoxia, and heavy metal toxicity, are highlighted and discussed in this review.

Global DNA Methylation and mRNA-miRNA Variations Activated by Heat Shock Boost Early Microspore Embryogenesis in Cabbage (Brassica oleracea)

Microspore culture, a type of haploid breeding, is extensively used in the cultivation of cruciferous crops such as cabbage. Heat shock (HS) treatment is essential to improve the embryo rate during the culture process; however, its molecular role in boosting early microspore embryogenesis (ME) remains unknown. Here we combined DNA methylation levels, miRNAs, and transcriptome profiles in isolated microspores of cabbage ‘01-88’ under HS (32 °C for 24 h) and normal temperature (25 °C for 24 h) to investigate the regulatory roles of DNA methylation and miRNA in early ME. Global methylation levels were significantly different in the two pre-treatments, and 508 differentially methylated regions (DMRs) were identified; 59.92% of DMRs were correlated with transcripts, and 39.43% of miRNA locus were associated with methylation levels. Significantly, the association analysis revealed that 31 differentially expressed genes (DEGs) were targeted by methylation and miRNA and were mainly involved in the reactive oxygen species (ROS) response and abscisic acid (ABA) signaling, indicating that HS induced DNA methylation, and miRNA might affect ME by influencing ROS and ABA. This study revealed that DNA methylation and miRNA interfered with ME by modulating key genes and pathways, which could broaden our understanding of the molecular regulation of ME induced by HS pre-treatment.

Cotton miR393-TIR1 Module Regulates Plant Defense Against Verticillium dahliae via Auxin Perception and Signaling

Plant auxin is essential in plant growth and development. However, the molecular mechanisms of auxin involvement in plant immunity are unclear. Here, we addressed the function of the cotton (Gossypium hirsutum) miR393-TIR1 module in plant defense against Verticillium dahliae infection via auxin perception and signaling. GhTIR1 was directedly cleaved by ghr-miR393 according to mRNA degradome data, 5'-RACE analysis, and a GUS reporter assay. Ghr-miR393 knockdown significantly increased plant susceptibility to V. dahliae compared to the control, while ghr-miR393 overexpression and GhTIR1 knockdown significantly increased plant resistance. External indole-3-acetic acid (IAA) application significantly enhanced susceptibility to V. dahliae in ghr-miR393 knockdown and control plants compared to mock treatment, and only slightly increased susceptibility in overexpressing ghr-miR393 and GhTIR1-silenced plants. Application of external PEO-IAA (an auxin antagonist) had a contrary trend with IAA application. Based on yeast two-hybrid and bimolecular fluorescence complementation assays, GhTIR1 interacted with GhIAA14 in the nucleus, and GhIAA14 knockdown reduced plant resistance to V. dahliae infection. The results suggested that the ghr-miR393-GhTIR1 module regulates plant defense via auxin perception and signaling. Additionally, simultaneous knockdown of GhTIR1 and GhICS1 significantly increased plant susceptibility to V. dahliae compared to the control, indicating that salicylic acid (SA) accumulation is vital for the ghr-miR393-GhTIR1 module to regulates plant resistance. Transcriptome data also demonstrated that GhTIR1 knockdown significantly downregulated expression of auxin-related genes and upregulated expression of SA-related genes. Overall, the ghr-miR393-GhTIR1 module participates in plant response to V. dahliae infection via IAA perception and signaling partially depending on the SA defense pathway.

The Characters of Non-Coding RNAs and Their Biological Roles in Plant Development and Abiotic Stress Response

Plant growth and development are greatly affected by the environment. Many genes have been identified to be involved in regulating plant development and adaption of abiotic stress. Apart from protein-coding genes, more and more evidence indicates that non-coding RNAs (ncRNAs), including small RNAs and long ncRNAs (lncRNAs), can target plant developmental and stress-responsive mRNAs, regulatory genes, DNA regulatory regions, and proteins to regulate the transcription of various genes at the transcriptional, posttranscriptional, and epigenetic level. Currently, the molecular regulatory mechanisms of sRNAs and lncRNAs controlling plant development and abiotic response are being deeply explored. In this review, we summarize the recent research progress of small RNAs and lncRNAs in plants, focusing on the signal factors, expression characters, targets functions, and interplay network of ncRNAs and their targets in plant development and abiotic stress responses. The complex molecular regulatory pathways among small RNAs, lncRNAs, and targets in plants are also discussed. Understanding molecular mechanisms and functional implications of ncRNAs in various abiotic stress responses and development will benefit us in regard to the use of ncRNAs as potential character-determining factors in molecular plant breeding.

Time-sequential detection of quantitative trait loci and candidate genes underlying the dynamic growth of Salix suchowensis

Elucidating the genetic factors underlying long-term biological processes remains challenging since the relevant genes and their effects may vary across different developmental stages. In this study, we carried out a large-scale field trial of the progeny of an F1 full-sib pedigree of Salix suchowensis and measured plant height and ground diameter periodically over a time course of 240 days. With the obtained data, we characterized plant growth rhythms and performed time-sequential analyses of quantitative trait loci underlying the dynamic growth of the plants. The dynamic mapping of quantitative trait loci revealed that stem height and ground diameter were under the control of four quantitative trait loci, and the effects of these quantitative trait loci varied greatly throughout the growth process, in which two quantitative trait loci were found to exert a pleiotropic effect determining the correlation between stem height and ground diameter. The analysis of candidate genes in the target genetic intervals showed that the pleiotropic effect of the two quantitative trait loci arises from the colocalization of genes with independent effects on stem height and ground diameter. Further examination of the expression patterns of the candidate genes indicated that height and circumference growth involve different activities of leaf and cambium tissues. This study provides unprecedented information to help us understand the dynamic growth of plants and presents an applicable strategy for elucidating the genetic mechanism underlying a long-term biological process by using plant growth as an example.

An Integration of MicroRNA and Transcriptome Sequencing Analysis Reveal Regulatory Roles of miRNAs in Response to Chilling Stress in Wild Rice

A chromosome single segment substitution line (CSSL) DC90, which was generated by introgressing CTS-12, a locus derived from common wild rice (Oryza rufipogon Griff.), into the 9311 (Oryza sativa L. ssp. indica) background, exhibits a chilling tolerance phenotype under chilling stress. Here, an integration of microRNA (miRNA) deep sequencing and transcriptomic sequencing analysis was performed to explore the expression profiles of miRNAs and their target genes mediated by CTS-12 under chilling stress, and to reveal the possible regulatory mechanisms of miRNAs that are involved in chilling tolerance. Integration analysis revealed that a number of differentially expressed miRNAs (DEMs) and putative target genes with different expression patterns and levels were identified in 9311 and DC90 under chilling stress. KEGG enrichment analysis revealed that the target genes that are regulated by chilling-induced miRNAs are involved in the regulation of various biological processes/pathways, including protein biosynthesis, redox process, photosynthetic process, and chloroplast development in two genotypes. CRISPR/Cas9 editing of the target genes of the key DEMs in a chilling tolerant rice variety Zhonghua 11 (ZH11) found that LOC_Os11g48020 (OsGL1-11), one of the putative target genes of osa-miR1846a/b-5p and encoding a wax synthesis protein, is correlated with a chilling stress tolerance phenotype, implying osa-miR1846a/b-5p/OsGL1-11 plays an important role in CTS-12-mediated chilling stress tolerance regulatory pathway(s). Therefore, we speculate that the CTS-12 may regulate the key miRNA target genes in response to chilling stress by differential regulation of miRNAs in wild rice, thereby resulting in the variation of chilling tolerance phenotype between 9311 and DC90.

Nanomaterials coupled with microRNAs for alleviating plant stress: a new opening towards sustainable agriculture

Plant growth and development is influenced by their continuous interaction with the environment. Their cellular machinery is geared to make rapid changes for adjusting the morphology and physiology to withstand the stressful changes in their surroundings. The present scenario of climate change has however intensified the occurrence and duration of stress and this is getting reflected in terms of yield loss. A number of breeding and molecular strategies are being adopted to enhance the performance of plants under abiotic stress conditions. In this context, the use of nanomaterials is gaining momentum. Nanotechnology is a versatile field and its application has been demonstrated in almost all the existing fields of science. In the agriculture sector, the use of nanoparticles is still limited, even though it has been found to increase germination and growth, enhance physiological and biochemical activities and impact gene expression. In this review, we have summarized the use and role of nanomaterial and small non-coding RNAs in crop improvement while highlighting the potential of nanomaterial assisted eco-friendly delivery of small non-coding RNAs as an innovative strategy for mitigating the effect of abiotic stress.

Genome-Wide Identification of Cotton (Gossypium spp.) Glycerol-3-Phosphate Dehydrogenase (GPDH) Family Members and the Role of GhGPDH5 in Response to Drought Stress

Glycerol-3-phosphate dehydrogenase (GPDH) is a key enzyme in plant glycerol synthesis and metabolism, and plays an important role in plant resistance to abiotic stress. Here, we identified 6, 7, 14 and 14 GPDH genes derived from Gossypium arboreum, Gossypium raimondii, Gossypium barbadense and Gossypium hirsutum, respectively. Phylogenetic analysis assigned these genes into three classes, and most of the genes within the family were expanded by whole-genome duplication (WGD) and segmental duplications. Moreover, determination of the nonsynonymous substitution rate/synonymous substitution rate (Ka/Ks) ratio showed that the GPDH had an evolutionary preference for purifying selection. Transcriptome data revealed that GPDH genes were more active in the early stages of fiber development. Additionally, numerous stress-related cis-elements were identified in the potential promoter region. Then, a protein-protein-interaction (PPI) network of GPDH5 in G. hirsutum was constructed. In addition, we predicted 30 underlying miRNAs in G. hirsutum. Functional validation results indicated that silencing GhGPDH5 diminished drought tolerance in the upland cotton TM-1 line. In summary, this study provides a fundamental understanding of the GPDH gene family in cotton, GhGPDH5 exerts a positive effect during drought stress and is potentially involved in stomatal closure movements.

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.

Genome-wide investigation of microRNAs and expression profiles during rhizome development in ginger (Zingiber officinale Roscoe)

BACKGROUND

MicroRNAs (miRNAs) are endogenous, non-coding small functional RNAs that govern the post-transcriptional regulatory system of gene expression and control the growth and development of plants. Ginger is an herb that is well-known for its flavor and medicinal properties. The genes involved in ginger rhizome development and secondary metabolism have been discovered, but the genome-wide identification of miRNAs and their overall expression profiles and targets during ginger rhizome development are largely unknown. In this study, we used BGISEQ-500 technology to perform genome-wide identification of miRNAs from the leaf, stem, root, flower, and rhizome of ginger during three development stages.

RESULTS

In total, 104 novel miRNAs and 160 conserved miRNAs in 28 miRNA families were identified. A total of 181 putative target genes for novel miRNAs and 2772 putative target genes for conserved miRNAs were predicted. Transcriptional factors were the most abundant target genes of miRNAs, and 17, 9, 8, 4, 13, 8, 3 conserved miRNAs and 5, 7, 4, 5, 5, 15, 9 novel miRNAs showed significant tissue-specific expression patterns in leaf, stem, root, flower, and rhizome. Additionally, 53 miRNAs were regarded as rhizome development-associated miRNAs, which mostly participate in metabolism, signal transduction, transport, and catabolism, suggesting that these miRNAs and their target genes play important roles in the rhizome development of ginger. Twelve candidate miRNA target genes were selected, and then, their credibility was confirmed using qRT-PCR. As the result of qRT-PCR analysis, the expression of 12 candidate target genes showed an opposite pattern after comparison with their miRNAs. The rhizome development system of ginger was observed to be governed by miR156, miR319, miR171a_2, miR164, and miR529, which modulated the expression of the SPL, MYB, GRF, SCL, and NAC genes, respectively.

CONCLUSION

This is a deep genome-wide investigation of miRNA and identification of miRNAs involved in rhizome development in ginger. We identified 52 rhizome-related miRNAs and 392 target genes, and this provides an important basis for understanding the molecular mechanisms of the miRNA target genes that mediate rhizome development in ginger.