Identification and characterization of salt responsive miRNA-SSR markers in rice (Oryza sativa)

Development of miRNA-SSR and target-SSR markers from yield-associate genes and their applicability in the assessment of genetic diversity and association mapping in rice (Oryza sativa L.)

The gene-derived functional markers are considered effective to use in marker-assisted breeding and genetic diversity analysis. As of now, no functional markers have been identified from miRNAs regulating yield traits. The miRNAs play a key role as regulators in controlling the candidate genes involved in grain yield improvement in rice. In this study, 13 miRNA-SSR and their target gene SSR markers were mined from 29 yield-responsive miRNA along with their 29 target genes in rice. The validation of these markers showed that four miRNA-SSRs and one target gene SSR markers had shown polymorphism among 120 diverse rice genotypes. The PIC values ranged from 0.25 (OsARF18-SSR) to 0.72 (miR408-SSR, miR172b-SSR, and miR396f-SSR) with an average value of 0.57. These polymorphic markers grouped 120 rice genotypes into 3 main clusters based on the levels of high genetic diversity. These markers also showed significant association with key yield traits. Among all, miR172b-SSR showed a strong association with plant height in two seasons. This investigation suggests that this new class of molecular markers has great potential in the characterization of rice germplasm by genetic diversity and population structure and in marker-assisted breeding for the development of high-yielding varieties.

Supplementary information

The online version contains supplementary material available at 10.1007/s11032-024-01462-z.

Identification of self- and pathogen-targeted miRNAs from resistant and susceptible Theobroma cacao variety to black pod disease

Cacao (Theobroma cacao) is a highly valuable crop with growing demand in the global market. However, cacao farmers often face challenges posed by black pod disease caused by Phytophthora spp., with P. palmivora being the most dominant. Regulations of various gene expressions influence plant resistance to pathogens. One mechanism involves targeting the mRNA of virulence genes in the invading pathogens, suppressing their infection. However, resistance also could be suppressed by plant-derived miRNAs that target their own defence genes. The objective of this study is to identify differentially expressed miRNAs in black pod-resistant and susceptible cacao varieties and to predict their targets in T. cacao and P. palmivora transcripts. Extracted miRNA from resistant and susceptible varieties of T. Cacao was sequenced, identified, and matched to host and pathogen mRNA. In total, 54 known miRNAs from 40 miRNA families and 67 novel miRNAs were identified. Seventeen miRNAs were differentially expressed in susceptible variety compared to resistant one, with 9 miRNAs upregulated and 8 miRNAs downregulated. In T. cacao transcripts, the upregulated miRNAs were predicted to target several genes, including defence genes. The suppression of these defense genes can lead to a reduction in plant resistance against pathogen infection. In P. palmivora transcripts, the upregulated miRNAs were predicted to target several genes, including P. palmivora effector genes. In the future, limiting expression of miRNAs that target T. cacao's defence genes and applying miRNAs that target P. palmivora effector genes hold promise for enhancing cacao plant resistance against P. palmivora infection.

Physiological and Molecular Responses of Pyrus pyraster Seedlings to Salt Treatment Analyzed by miRNA and Cytochrome P450 Gene-Based Markers

Physiological and molecular marker-based changes were studied in the tissues of two-year-old Pyrus pyraster (L.) Burgsd. seedlings under salt treatment. For 60 days, 5 mL of 100 mM NaCl solution was applied to each plant per day to a cumulative volume of 300 mL in the substrate. In response to osmotic stress, the seedlings increased their water use efficiency (WUE) on day 20 of regular NaCl application and maintained a stable net photosynthetic rate (An) per unit area. Under conditions of increasing salinity, the young plants maintained a balanced water regime of the leaf tissues (Ψwl). The seedlings invested mass to their root growth (R/S), retained a substantial portion (72%) of Na+ ions in the roots, and protected their leaves against intoxication and damage. A significant decrease in the leaf gas exchange parameters (gs, E, An) was manifested on day 60 of the experiment when the cumulative NaCl intake was 300 mL per plant. The variability in the reactions of the seedlings to salinity is related to the use of open-pollinated progeny (54 genotypes) in the experiment. Lus-miR168 showed tissue- and genotype-specific genome responses to the applied stress. Polymorphic miRNA-based loci were mostly detected in the root samples on the 20th and 35th days of the experiment. The cumulative effect of the salt treatment was reflected in the predominance of polymorphic loci in the leaves. We can confirm that miRNA-based markers represent a sensitive detection tool for plant stress response on an individual level. The screening and selection of the optimal type of miRNA for this type of research is crucial. The cytochrome P450-Based Analog (PBA) techniques were unable to detect polymorphism among the control and treated seedlings, except for the primer pair CYP2BF+R, where, in the roots of the stressed plant, insertions in the amplicons were obtained. The expression ratios of cytochrome P450 in the salt-stressed plants were higher in the roots in the case of 20/100 mL and in the leaves with higher doses. The observed physiological and molecular responses to salinity reflect the potential of P. pyraster seedlings in adaptation to osmotic and ionic stress.

High Daytime Temperature Responsive MicroRNA Profiles in Developing Grains of Rice Varieties with Contrasting Chalkiness

High temperature impairs starch biosynthesis in developing rice grains and thereby increases chalkiness, affecting the grain quality. Genome encoded microRNAs (miRNAs) fine-tune target transcript abundances in a spatio-temporal specific manner, and this mode of gene regulation is critical for a myriad of developmental processes as well as stress responses. However, the role of miRNAs in maintaining rice grain quality/chalkiness during high daytime temperature (HDT) stress is relatively unknown. To uncover the role of miRNAs in this process, we used five contrasting rice genotypes (low chalky lines Cyp, Ben, and KB and high chalky lines LaGrue and NB) and compared the miRNA profiles in the R6 stage caryopsis samples from plants subjected to prolonged HDT (from the onset of fertilization through R6 stage of caryopsis development). Our small RNA analysis has identified approximately 744 miRNAs that can be grouped into 291 families. Of these, 186 miRNAs belonging to 103 families are differentially regulated under HDT. Only two miRNAs, Osa-miR444f and Osa-miR1866-5p, were upregulated in all genotypes, implying that the regulations greatly varied between the genotypes. Furthermore, not even a single miRNA was commonly up/down regulated specifically in the three tolerant genotypes. However, three miRNAs (Osa-miR1866-3p, Osa-miR5150-3p and canH-miR9774a,b-3p) were commonly upregulated and onemiRNA (Osa-miR393b-5p) was commonly downregulated specifically in the sensitive genotypes (LaGrue and NB). These observations suggest that few similarities exist within the low chalky or high chalky genotypes, possibly due to high genetic variation. Among the five genotypes used, Cypress and LaGrue are genetically closely related, but exhibit contrasting chalkiness under HDT, and thus, a comparison between them is most relevant. This comparison revealed a general tendency for Cypress to display miRNA regulations that could decrease chalkiness under HDT compared with LaGrue. This study suggests that miRNAs could play an important role in maintaining grain quality in HDT-stressed rice.

Role of non-coding RNAs against salinity stress in Oryza species: Strategies and challenges in analyzing miRNAs, tRFs and circRNAs

Salinity is an imbalanced concentration of mineral salts in the soil or water that causes yield loss in salt-sensitive crops. Rice plant is vulnerable to soil salinity stress at seedling and reproductive stages. Different non-coding RNAs (ncRNAs) post-transcriptionally regulate different sets of genes during different developmental stages under varying salinity tolerance levels. While microRNAs (miRNAs) are well known small endogenous ncRNAs, tRNA-derived RNA fragments (tRFs) are an emerging class of small ncRNAs derived from tRNA genes with a demonstrated regulatory role, like miRNAs, in humans but unexplored in plants. Circular RNA (circRNA), another ncRNA produced by back-splicing events, acts as target mimics by preventing miRNAs from binding with their target mRNAs, thereby reducing the miRNA's action upon its target. Same may hold true between circRNAs and tRFs. Hence, the work done on these ncRNAs was reviewed and no reports were found for circRNAs and tRFs under salinity stress in rice, either at seedling or reproductive stages. Even the reports on miRNAs are restricted to seedling stage only, in spite of severe effects on rice crop production due to salt stress during reproductive stage. Moreover, this review sheds light on strategies to predict and analyze these ncRNAs in an effective manner.

Suppression of microRNA168 enhances salt tolerance in rice (Oryza sativa L.)

BACKGROUND

Rice is a salt-sensitive crop. Complex gene regulatory cascades are likely involved in salinity stress in rice roots. microRNA168 (miR168) is a conserved miRNA among different plant species. It in-directly regulates the expression of all miRNAs by targeting gene ARGONAUTE1(AGO1). Short Tandem Target Mimic (STTM) technology is an ideal approach to study miRNA functions by in-activating mature miRNA in plants.

RESULTS

In this study, rice miR168 was inactivated by STTM. The T3 generation seedlings of STTM168 exhibited significantly enhanced salt resistance. Direct target genes of rice miR168 were obtained by in silico prediction and further confirmed by degradome-sequencing. PINHEAD (OsAGO1), which was previously suggested to be a plant abiotic stress response regulator. RNA-Seq was performed in root samples of 150mM salt-treated STTM168 and control seedlings. Among these screened 481 differentially expressed genes within STTM168 and the control, 44 abiotic stress response related genes showed significant difference, including four known salt-responsive genes.

CONCLUSION

Based on sequencing and qRT-PCR, a "miR168-AGO1-downstream" gene regulation model was proposed to be responsible for rice salt stress response. The present study proved miR168-AGO1 cascade to play important role in rice salinity stress responding, as well as to be applied in agronomic improvement in further.

In silico analysis of key regulatory networks related to microfibril angle in Populus trichocarpa Hook

Dissection of regulatory network that control wood structure is highly challenging in functional genomics. Nevertheless, due to the availability of genomic, transcriptomic and proteomic sequences, a large amount of information is available for use in achieving this goal. MicroRNAs, which compose a class of small non-coding RNA molecules that inhibit protein translation by targeting mRNA cleavage sites and thus regulate a wide variety of developmental and physiological processes in plants, are important parts of this regulatory network. These findings and the availability of sequence information have made it possible to carry out an in silico analysis to predict and annotate miRNAs and their target genes associated with an important factor affecting wood rigidity, microfibril angle (MFA), throughout the Populus trichocarpa Hook. genome. Our computational approach revealed miRNAs and their targets via ESTs, sequences putatively associated with microfibril angle. In total, 250 miRNAs were identified as RNA molecules with roles in the silencing and post-transcriptional regulation of the expression of nine genes. We found SHY2, IAA4 (ATAUX2–11), BZIP60, AP2, MYB15, ABI3, MYB17, LAF1 and MYB28 as important nodes in a network with possible role in MFA determination. Other co-expressed genes putatively involved in this regulatory system were also identified by construction of a co-expression network. The candidate genes from this study may help unravel the regulatory networks putatively linked to microfibril angle.

Understanding natural genetic variation for grain phytic acid content and functional marker development for phytic acid-related genes in rice

BACKGROUND

The nutritional value of rice can be improved by developing varieties with optimum levels of grain phytic acid (PA). Artificial low-PA mutants with impaired PA biosynthesis have been developed in rice through induced mutagenesis. However, low-PA mutant stocks with drastically reduced grain PA content have poor breeding potential, and their use in rice breeding is restricted due to their detrimental pleiotropic effects, which include decreased seed viability, low grain weight, and low seed yield. Therefore, it is necessary to take advantage of the natural variation in grain PA content in order to reduce the PA content to an ideal level without compromising the crop's agronomic performance. Natural genetic diversity in grain PA content has not been thoroughly examined among elite genetic stocks. Additionally, given grain PA content as a quantitative trait driven by polygenes, DNA marker-assisted selection may be required for manipulation of such a trait; however, informative DNA markers for PA content have not yet been identified in rice. Here we investigated and dissected natural genetic variation and genetic variability components for grain PA content in rice varieties cultivated in Eastern and North-Eastern India during the last 50 years. We developed novel gene-based markers for the low-PA-related candidate genes in rice germplasm, and their allelic diversity and association with natural variation in grain PA content were studied.

RESULTS

A wide (0.3-2.8%), significant variation for grain PA content, with decade-wise and ecology-wise differences, was observed among rice varieties. Significant genotype x environment interaction suggested polygenic inheritance. The novel candidate gene-based markers detected 43 alleles in the rice varieties. The new markers were found highly informative as indicated by PIC values (0.11-0.65; average: 0.34) and coverage of total diversity. Marker alleles developed from two putative transporter genes viz., SPDT and OsPT8 were significantly associated with grain PA variation assayed on the panel. A 201 bp allele at the 3' UTR of SPDT gene was negatively associated with grain PA content and explained 7.84% of the phenotypic variation. A rare allele in the coding sequence of OsPT8 gene was positively associated with grain PA content which explained phenotypic variation of 18.49%.

CONCLUSION

Natural variation in grain PA content is substantial and is mostly controlled by genetic factors. The unique DNA markers linked with PA content have significant potential as genomic resources for the development of low-PA rice varieties through genomics-assisted breeding procedures.

Development and verification of SSR markers from drought stress-responsive miRNAs in Dongxiang wild rice (Oryza rufipogon Griff.)

Rice production worldwide has continued to decline due to various environmental stresses, with drought stress being a prominent factor, as rice is a semi-aquatic plant. Thus, development of drought stress-resistant rice varieties is of great importance for rice production. In our previous study, we found that microRNAs (miRNAs) play a crucial role in the response to drought stress in Dongxiang wild rice (DXWR) (Oryza rufipogon Griff.). Developing drought stress-responsive miRNA-based single sequence repeat (SSR) markers for DXWR will be of great value for the efficient identification and utilization of miRNA genes to breed drought stress-resistant rice varieties. In this study, ninety-nine novel SSR markers were developed based on the drought stress-responsive miRNAs of DXWR. These markers were distributed in all 12 rice chromosomes, and most were in chromosomes 2 and 6, with di- and tri-nucleotides being the most abundant repeat motifs. Twelve out of fourteen synthesized markers displayed high levels of genetic diversity in the genomes of three populations of DXWR and 40 modern rice varieties worldwide. The number of alleles per locus ranged from 2 to 7, with an average of 4.67; the genetic diversity index ranged from 0.21 to 0.76, with an average of 0.58; and the polymorphism information content value ranged from 0.18 to 0.72, with an average of 0.53. These novel molecular markers developed from the drought stress-responsive miRNAs of DXWR could be additional tools for mapping elite miRNA genes and breeding drought stress-resistant rice varieties.

Development and characterization of nitrogen and phosphorus use efficiency responsive genic and miRNA derived SSR markers in wheat

Among all the nutrients, nitrogen (N) and phosphorous (P) are the most limiting factors reducing wheat production and productivity world-wide. These macronutrients are directly applied to soil in the form of fertilizers. However, only 30-40% of these applied fertilizers are utilized by crop plants, while the rest is lost through volatilization, leaching, and surface run off. Therefore, to overcome the deficiency of N and P, it becomes necessary to improve their use efficiency. Marker-assisted selection (MAS) combined with traditional plant breeding approaches is considered best to improve the N and P use efficiency (N/PUE) of wheat varieties. In this study, we developed and evaluated a total of 98 simple sequence repeat (SSR) markers including 66 microRNAs and 32 gene-specific SSRs on a panel of 10 (N and P efficient/deficient) wheat genotypes. Out of these, 35 SSRs were found polymorphic and have been used for the study of genetic diversity and population differentiation. A set of two SSRs, namely miR171a and miR167a were found candidate markers able to discriminate contrasting genotypes for N/PUE, respectively. Therefore, these two markers could be used as functional markers for characterization of wheat germplasm for N and P use efficiency. Target genes of these miRNAs were found to be highly associated with biological processes (24 GO terms) as compared to molecular function and cellular component and shows differential expression under various P starving conditions and abiotic stresses.

Deciphering the role of miRNA in reprogramming plant responses to drought stress

Drought is the most prevalent environmental stress that affects plants' growth, development, and crop productivity. However, plants have evolved adaptive mechanisms to respond to the harmful effects of drought. They reprogram their: transcriptome, proteome, and metabolome that alter their cellular and physiological processes and establish cellular homeostasis. One of the crucial regulatory processes that govern this reprogramming is post-transcriptional regulation by microRNAs (miRNAs). miRNAs are small non-coding RNAs, involved in the downregulation of the target mRNA via translation inhibition/mRNA degradation/miRNA-mediated mRNA decay/ribosome drop off/DNA methylation. Many drought-inducible miRNAs have been identified and characterized in plants. Their main targets are regulatory genes that influence growth, development, osmotic stress tolerance, antioxidant defense, phytohormone-mediated signaling, and delayed senescence during drought stress. Overexpression of drought-responsive miRNAs (Osa-miR535, miR160, miR408, Osa-miR393, Osa-miR319, and Gma-miR394) in certain plants has led to tolerance against drought stress indicating their vital role in stress mitigation. Similarly, knock down (miR166/miR398c) or deletion (miR169 and miR827) of miRNAs has also resulted in tolerance to drought stress. Likewise, engineered Arabidopsis plants with miR165, miR166 using short tandem target mimic strategy, exhibited drought tolerance. Since miRNAs regulate the expression of an array of drought-responsive genes, they can act as prospective targets for genetic manipulations to enhance drought tolerance in crops and achieve sustainable agriculture. Further investigations toward functional characterization of diverse miRNAs, and understanding stress-responses regulated by these miRNAs and their utilization in biotechnological applications is highly recommended.

Novel miRNA-SSRs for Improving Seed Hardness Trait of Pomegranate (Punica granatum L.)

Present research discovered novel miRNA-SSRs for seed type trait from 761 potential precursor miRNA sequences of pomegranate. SSR mining and BLASTx of the unique sequences identified 69 non-coding pre-miRNA sequences, which were then searched for BLASTn homology against Dabenzi genome. Sixty three true pri-miRNA contigs encoding 213 pre-miRNAs were predicted. Analysis of the resulting sequences enabled discovery of SSRs within pri-miRNA (227) and pre-miRNA sequences (79). A total of 132 miRNA-SSRs were developed for seed type trait from 63 true pri-miRNAs, of which 46 were specific to pre-miRNAs. Through ePCR, 123 primers were validated and mapped on eight Tunisia chromosomes. Further, 80 SSRs producing specific amplicons were ePCR-confirmed on multiple genomes i.e. Dabenzi, Taishanhong, AG2017 and Tunisia, yielding a set of 63 polymorphic SSRs (polymorphism information content ≥0.5). Of these, 32 miRNA-SSRs revealed higher polymorphism level (89.29%) when assayed on six pomegranate genotypes. Furthermore, target prediction and network analysis suggested a possible association of miRNA-SSRs i.e. miRNA_SH_SSR69, miRNA_SH_SSR36, miRNA_SH_SSR103, miRNA_SH_SSR35 and miRNA_SH_SSR53 with seed type trait. These miRNA-SSRs would serve as important genomic resource for rapid and targeted improvement of seed type trait of pomegranate.

Genome wide study of cysteine rich receptor like proteins in Gossypium sp

Cysteine-rich receptor-like-kinases (CRKs), a transmembrane subfamily of receptor-like kinase, play crucial roles in plant adaptation. As such cotton is the major source of fiber for the textile industry, but environmental stresses are limiting its growth and production. Here, we have performed a deep computational analysis of CRKs in five Gossypium species, including G. arboreum (60 genes), G. raimondii (74 genes), G. herbaceum (65 genes), G. hirsutum (118 genes), and G. barbadense (120 genes). All identified CRKs were classified into 11 major classes and 43 subclasses with the finding of several novel CRK-associated domains including ALMT, FUSC_2, Cript, FYVE, and Pkinase. Of these, DUF26_DUF26_Pkinase_Tyr was common and had elevated expression under different biotic and abiotic stresses. Moreover, the 35 land plants comparison identified several new CRKs domain-architectures. Likewise, several SNPs and InDels were observed in CLCuD resistant G. hirsutum. The miRNA target side prediction and their expression profiling in different tissues predicted miR172 as a major CRK regulating miR. The expression profiling of CRKs identified multiple clusters with co-expression under certain stress conditions. The expression analysis under CLCuD highlighted the role of GhCRK057, GhCRK059, GhCRK058, and GhCRK081 in resistant accession. Overall, these results provided primary data for future potential functional analysis as well as a reference study for other agronomically important crops.

NCodR: A multi-class support vector machine classification to distinguish non-coding RNAs in Viridiplantae

Non-coding RNAs (ncRNAs) are major players in the regulation of gene expression. This study analyses seven classes of ncRNAs in plants using sequence and secondary structure-based RNA folding measures. We observe distinct regions in the distribution of AU content along with overlapping regions for different ncRNA classes. Additionally, we find similar averages for minimum folding energy index across various ncRNAs classes except for pre-miRNAs and lncRNAs. Various RNA folding measures show similar trends among the different ncRNA classes except for pre-miRNAs and lncRNAs. We observe different k-mer repeat signatures of length three among various ncRNA classes. However, in pre-miRs and lncRNAs, a diffuse pattern of k-mers is observed. Using these attributes, we train eight different classifiers to discriminate various ncRNA classes in plants. Support vector machines employing radial basis function show the highest accuracy (average F1 of ~96%) in discriminating ncRNAs, and the classifier is implemented as a web server, NCodR.

Identification of novel drought-responsive miRNA regulatory network of drought stress response in common vetch (Vicia sativa)

Drought is among the most important natural disasters with severe effects on animals and plants. MicroRNAs are a class of noncoding RNAs that play a crucial role in plant growth, development, and response to stress factors, including drought. However, the microRNAs in drought responses in common vetch (Vicia sativa), an annual herbaceous leguminous plant commonly used for forage by including it in mixed seeding during winter and spring, have not been characterized. To explore the microRNAs’ response to drought in common vetch, we sequenced 10 small RNA (sRNA) libraries by the next-generation sequencing technology. We obtained 379 known miRNAs belonging to 38 families and 47 novel miRNAs. The two groups had varying numbers of differentially expressed miRNAs: 85 in the comparison group D5 vs C5 and 38 in the comparison group D3 vs C3. Combined analysis of mRNA and miRNA in the same samples under drought treatment identified 318 different target genes of 123 miRNAs. Functional annotation of the target genes revealed that the miRNAs regulate drought-responsive genes, such as leucine-rich repeat receptor-like kinase-encoding genes (LRR-RLKs), ABC transporter G family member 1 (ABCG1), and MAG2-interacting protein 2 (MIP2). The genes were involved in various pathways, including cell wall biosynthesis, reactive oxygen removal, and protein transport. The findings in this study provide new insights into the miRNA-mediated regulatory networks of drought stress response in common vetch.

Targeted designing functional markers revealed the role of retrotransposon derived miRNAs as mobile epigenetic regulators in adaptation responses of pistachio

We developed novel miRNA-based markers based on salt responsive miRNA sequences to detect polymorphisms in miRNA sequences and locations. The validation of 76 combined miRNA + miRNA and miRNA + ISSR markers in the three extreme pistachio populations led to the identification of three selected markers that could link salt tolerance phenotype to genotype and divided pistachio genotypes and Pistacia species into three clusters. This novel functional marker system, in addition to more efficient performance, has higher polymorphisms than previous miRNA-based marker systems. The functional importance of the target gene of five miRNAs in the structure of the three selected markers in regulation of different genes such as ECA2, ALA10, PFK, PHT1;4, PTR3, KUP2, GRAS, TCP, bHLH, PHD finger, PLATZ and genes involved in developmental, signaling and biosynthetic processes shows that the polymorphism associated with these selected miRNAs can make a significant phenotypic difference between salt sensitive and tolerant pistachio genotypes. The sequencing results of selected bands showed the presence of conserved miRNAs in the structure of the mitochondrial genome. Further notable findings of this study are that the sequences of PCR products of two selected markers were annotated as Gypsy and Copia retrotransposable elements. The transposition of retrotransposons with related miRNAs by increasing the number of miRNA copies and changing their location between nuclear and organellar genomes can affect the regulatory activity of these molecules. These findings show the crucial role of retrotransposon-derived miRNAs as mobile epigenetic regulators between intracellular genomes in regulating salt stress responses as well as creating new and tolerant phenotypes for adaptation to environmental conditions.

Current status of microRNA-mediated regulation of drought stress responses in cereals

Drought is one of the most important abiotic stress factors impeding crop productivity. With the uncovering of their role as potential regulators of gene expression, microRNAs (miRNAs) have been recognized as new targets for developing stress resistance. MicroRNAs are small noncoding RNAs whose abundance is significantly altered under stress conditions. Interestingly, plant miRNAs predominantly targets transcription factors (TFs), and some of which are also the most critical drought-responsive genes that in turn could regulate the expression of numerous loci with drought-adaptive potential. The phytohormone ABA plays important roles in regulating stomatal conductance and in initiating an adaptive response to drought stress. miRNAs are implicated in regulating ABA-(abscisic acid) and non-ABA-mediated drought resistance pathways. For instance, miR159-MYB module and miR169-NFYA module participates in an ABA-dependent pathway, whereas several other ABA-independent miRNA-target modules (miR156-SPL; miR393-TIR1; miR160-ARF10, ARF16, ARF17; miR167-ARF6 and ARF8; miR390/TAS3siRNA-ARF2, ARF3, ARF4) collectively regulate drought responses in plants. Overall, miRNA-mediated drought response manifests diverse molecular, biochemical and physiological processes. Because of their immense role in controlling gene expression, miRNA manipulation has significant potential to augment plant tolerance to drought stress. This review compiles the current understanding of drought-responsive miRNAs in major cereals. Also, potential miRNA manipulation strategies currently in use along with the challenges and future perspectives are discussed.

Improving rice salt tolerance by precision breeding in a new era

Rice is a premier staple food that constitutes the bulk of the daily diet of the majority of people in Asia. Agricultural productivity must be boosted to support this huge demand for rice. However, production is jeopardized by soil salinity. Advances in whole-genome sequencing, marker-assisted breeding strategies, and targeted mutagenesis have substantially improved the toolbox of today's breeders. Given that salinity has a major influence on rice at both the seedling and reproductive stages, understanding and manipulating this trait will have an enormous impact on sustainable production. This article summarizes recent developments in the understanding of the mechanisms of salt tolerance and how state-of-the-art tools such as RNA guided CRISPR endonuclease technology including targeted mutagenesis or base and prime editing can help in gene discovery and functional analysis as well as in transferring favorable alleles into elite breeding material to speed the breeding of salt-tolerant rice cultivars.

Development and use of miRNA-derived SSR markers for the study of genetic diversity, population structure, and characterization of genotypes for breeding heat tolerant wheat varieties

Heat stress is an important abiotic factor that limits wheat production globally, including south-east Asia. The importance of micro (mi) RNAs in gene expression under various biotic and abiotic stresses is well documented. Molecular markers, specifically simple sequence repeats (SSRs), play an important role in the wheat improvement breeding programs. Given the role of miRNAs in heat stress-induced transcriptional regulation and acclimatization, the development of miRNA-derived SSRs would prove useful in studying the allelic diversity at the heat-responsive miRNA-genes in wheat. In the present study, efforts have been made to identify SSRs from 96 wheat heat-responsive miRNA-genes and their characterization using a panel of wheat genotypes with contrasting reactions (tolerance/susceptible) to heat stress. A set of 13 miRNA-derived SSR markers were successfully developed as an outcome. These miRNA-SSRs are located on 11 different common wheat chromosomes (2A, 3A, 3B, 3D, 4D, 5A, 5B, 5D, 6A, 6D, and 7A). Among 13 miRNA-SSRs, seven were polymorphic on a set of 37 selected wheat genotypes. Within these polymorphic SSRs, three makers, namely HT-169j, HT-160a, and HT-160b, were found promising as they could discriminate heat-tolerant and heat-susceptible genotypes. This is the first report of miRNA-SSR development in wheat and their deployment in genetic diversity and population structure studies and characterization of trait-specific germplasm. The study suggests that this new class of molecular makers has great potential in the marker-assisted breeding (MAB) programs targeted at improving heat tolerance and other adaptability or developmental traits in wheat and other crops.

Identification and development of novel salt-responsive candidate gene based SSRs (cg-SSRs) and MIR gene based SSRs (mir-SSRs) in bread wheat (Triticum aestivum)

Salt stress adversely affects the global wheat production and productivity. To improve salinity tolerance of crops, identification of robust molecular markers is highly imperative for development of salt-tolerant cultivars to mimic yield losses under saline conditions. In this study, we mined 171 salt-responsive genes (including 10 miRNAs) from bread wheat genome using the sequence information of functionally validated salt-responsive rice genes. Salt-stress, tissue and developmental stage-specific expression analysis of RNA-seq datasets revealed the constitutive as well as the inductive response of salt-responsive genes in different tissues of wheat. Fifty-four genotypes were phenotyped for salt stress tolerance. The stress tolerance index of the genotypes ranged from 0.30 to 3.18. In order to understand the genetic diversity, candidate gene based SSRs (cg-SSRs) and MIR gene based SSRs (miR-SSRs) were mined from 171 members of salt-responsive genes of wheat and validated among the contrasting panels of 54 tolerant as well as susceptible wheat genotypes. Among 53 SSR markers screened, 10 cg-SSRs and 8 miR-SSRs were found to be polymorphic. Polymorphic information content between the wheat genotypes ranged from 0.07 to 0.67, indicating the extant of wide genetic variation among the salt tolerant and susceptible genotypes at the DNA level. The genetic diversity analysis based on the allelic data grouped the wheat genotypes into three separate clusters of which single group encompassing most of the salt susceptible genotypes and two of them containing salt tolerance and moderately salt tolerance wheat genotypes were in congruence with penotypic data. Our study showed that both salt-responsive genes and miRNAs based SSRs were more diverse and can be effectively used for diversity analysis. This study reports the first extensive survey on genome-wide analysis, identification, development and validation of salt-responsive cg-SSRs and miR-SSRs in wheat. The information generated in the present study on genetic divergence among genotypes having a differential response to salt will help in the selection of suitable lines as parents for developing salt tolerant cultivars in wheat.

Development and validation of heat-responsive candidate gene and miRNA gene based SSR markers to analysis genetic diversity in wheat for heat tolerance breeding

Being a major staple food crop of the world, wheat provides nutritional food security to the global populations. Heat stress is a major abiotic stress that adversely affects wheat production throughout the world including Indo-Gangatic Plains (IGP) where four wheat growing countries viz., India, Bangladesh, Nepal and Pakistan produce 42% of the total wheat production. Therefore, identification of heat stress responsive molecular markers is imperative to marker assisted breeding programs. Information about trait specific gene based SSRs is available but there is lack of information on SSRs from non-coding regions. In the present study, we developed 177 heat-responsive gene-based SSRs (cg-SSR) and MIR gene-based SSR (miRNA-SSR) markers from wheat genome for assessing genetic diversity analysis of thirty- six contrasting wheat genotypes for heat tolerance. Of the 177 SSR loci, 144 yielded unambiguous and repeatable amplicons, however, thirty-seven were found polymorphic among the 36 wheat genotypes. The polymorphism information content (PIC) of primers used in this study ranged from 0.03-0.73, with a mean of 0.35. Number of alleles produced per primer varied from 2 to 6, with a mean of 2.58. The UPGMA dendrogram analysis grouped all wheat genotypes into four clusters. The markers developed in this study has potential application in the MAS based breeding programs for developing heat tolerant wheat cultivars and genetic diversity analysis of wheat germplasm. Identification of noncoding region based SSRs will be fruitful for identification of trait specific wheat germplasm.

Identification of lncRNAs in response to infection by Plasmodiophora brassicae in Brassica napus and development of lncRNA-based SSR markers

Clubroot resistance in spring canola has been introgressed from different Brassica sources; however, molecular mechanism underlying this resistance, especially the involvement of long non-coding RNAs (lncRNAs), is yet to be understood. We identified 464 differentially expressed (DE) lncRNAs from the roots of clubroot-resistant canola, carrying resistance on chromosome BnaA03, and susceptible canola lines challenged with Plasmodiophora brassicae pathotype 3. Pathway enrichment analysis showed that most of the target genes regulated by these DE lncRNAs belonged to plant-pathogen interaction and hormone signaling, as well as primary and secondary metabolic pathways. Comparative analysis of these lncRNAs with 530 previously reported DE lncRNAs, identified using resistance located on BnaA08, detected 12 lncRNAs that showed a similar trend of upregulation in both types of resistant lines; these lncRNAs probably play a fundamental role in clubroot resistance. We identified SSR markers within 196 DE lncRNAs. Genotyping of two DH populations carrying resistance on BnaA03 identified a marker capable of detecting the resistance in 98% of the DH lines. To our knowledge, this is the first report of the identification of SSRs within lncRNAs responsive to P. brassicae infection, demonstrating the potential use of lncRNAs in the breeding of Brassica crops.

Genetic variability analysis of partially salt tolerant local and inbred rice (Oryza sativa L.) through molecular markers

Random Amplified Polymorphic DNA (RAPD) analysis was performed to assess the genetic variability in sixteen selected germplasms of rice, Oryza sativa L. using eight decamer RAPD primers. The data obtained from this investigation reveals a high level of polymorphism between cultivars. The primers produced a total of 255 bands of which all 255 bands were polymorphic indicating 100% polymorphism. The size of the amplified bands ranged from 220 bp to 2290 bp. The number of polymorphic fragments ranged from 24 to 49 with an average of 32 polymorphic fragments for each primer. The primer OPX04 produced the maximum number (49) of polymorphic bands while the OPB04 and OPB17 produced the minimum number (24) of polymorphic bands. The polymorphic information content (PIC) values ranged from 0.6616 to 0.8845 with an average of 0.832. The highest PIC value (0.8845) was obtained for primer OPL03. The RAPD data was analyzed to determine the pair-wise genetic similarity coefficients which ranged from 0.00 to 0.83. The BRRIdhan 23 and the BRRIdhan 41 varieties were the closest genotypes with the highest similarity index of 83%. This was followed by 77% similarity between a pair of cultivars Kalamona and Horkuch. On the other hand, 100% dissimilarity was seen between BRRIdhan 53, BRRIdhan 50, BRRIdhan 10, BRRIdhan 70, BRRIdhan 54, BRRIdhan 40, BRRIdhan 23, BRRIdhan 47, BRRIdhan 41 and Dadsail respectively and between BRRIdhan 53 and Horkuch; indicating a high level of variability between paired genotypes. Cluster analysis was performed using Unweighted Paired Group of Arithmetic Means (UPGMA). The UPGMA dendrogram resolved the selected rice cultivars into four clusters.

Properties of Ginkgo biloba L.: Antioxidant Characterization, Antimicrobial Activities, and Genomic MicroRNA Based Marker Fingerprints

The aim of this study was to characterize extracts from the leaves of Ginkgo biloba L. from selected Slovakian localities in terms of the content of bioactive constituents, antioxidants and their antimicrobial properties. The results indicated that the content of antioxidants was sample-specific, and this specificity was statistically significant. Ginkgo biloba L. from the locality of Košice had the best activity determined by the free radical scavenging activity (DPPH) (1.545 mg Trolox equivalent antioxidant capacity (TEAC)/g fresh matter (FM)) as well as the molybdenum-reducing antioxidant power (35.485 mg TEAC/g FM) methods. The highest content of total polyphenols (2.803 mg gallic acid equivalent (GAE)/g FM) and flavonoids (4.649 μg quercetin equivalent (QE)/g FM) was also detected in this sample. All samples of G. biloba leaf extracts showed significant antimicrobial activity against one or more of the examined bacterial species, and Staphylococcus aureus subsp. aureus CCM 2461 was found to be the most susceptible (minimal inhibition concentration MIC50 and MIC90 values of 64.2 and 72.2 µg/mL, respectively). Based on the results it was concluded that Ginkgo biloba L. extracts can be used as antimicrobial and antioxidant additives. Selected miRNA-based molecular markers were used to examine the environmental adaptability of Ginkgo biloba L. An almost-complete genotype clustering pattern based on locality was determined in the analysis that involved a species-specific gb-miR5261 marker. Morphologically specific exemplar, cv. Ohatsuki, was excluded.

Genome wide identification, characterization and validation of novel miRNA-based SSR markers in pomegranate (Punica granatum L.)

A total of 17,439 mature miRNAs (~ 21 nt) earlier generated through RNA seq in the pomegranate were used for in silico analysis. After complexity reduction, a total of 1922 representative mature miRNAs were selected and used as query sequences against pomegranate genome to retrieve 2540 homologous contigs with flanking regions (~ 800). By using pre-miRNA prediction web server, a total of 1028 true contigs harbouring pri-miRNAs encoding 1162 pre-miRNAs were identified. Survey of these sequences for SSRs yielded a total of 1358 and 238 SSRs specific to pri-miRNA and pre-miRNAs, respectively. Of these, primer pairs were designed for 897 pri-miRNA and 168 pre-miRNA SSRs. In pri-miRNA sequences, hexa-nucleotides repeats were found to be most abundant (44.18%) followed by mono- (18.41%) and di-nucleotide (17.01%), which is also observed in pre-miRNA sequences. Further, a set of 51 randomly selected pre-miRNA-SSRs was examined for marker polymorphism. The experimental validation of these markers on eight pomegranate genotypes demonstrated 92.15% polymorphism. Utility of these functional markers was confirmed via examination of genetic diversity of 18 pomegranate genotypes using 15 miRNA-SSRs. Further, potential application of miRNA-SSRs for discovery of trait specific candidate genes was showed by validating 51 mature miRNA against publically available 2047 EST sequences of pomegranate by target and network analysis. In summary, the current study offers novel functional molecular markers for pomegranate genetic improvement.

Plant microRNAs: biogenesis, gene silencing, web-based analysis tools and their use as molecular markers

MicroRNAs (miRNAs) are tiny (20-24 nt bp) regulatory non-protein-coding RNA molecules that have been extensively characterized and found important for many physiological and developmental processes. The miss-expression of miRNAs leads to various defects in plants. MicroRNAs repress gene expression by directing mRNA degradation or translational arrest. Several proteins such as PP43A, HYL1, DCL, HST are indispensable role players in promoting miRNA biogenesis in plants. During miRNA biogenesis, lariat RNAs are produced as by-products of pre-mRNA splicing which have a negative role in regulation of miRNA homeostasis. By acting as a decoy and by sequestering to the dicing complex, lariat RNA can prevent the processing of miRNAs. A number of bioinformatic tools with different methodologies are available to identify and validate miRNAs and their targets. Many miRNAs have been reported in different crops for different traits; however, no reports are available on their use in plant breeding. Recently, researchers have developed trait specific miRNA-based molecular markers (miRNA-SSRs/SNP) for many quantitative traits in different plant species. In the future, these molecular markers can be used for plant breeding programs. In this review, a comprehensive up-to-date information is provided on the bioinformatic tools used for analysis of plant miRNAs and their targets, the number of miRNAs, their biogenesis, gene silencing mechanism and miRNA-based molecular markers.

Analysis of genetic diversity and population structure using SSR markers and validation of a Cleavage Amplified Polymorphic Sequences (CAPS) marker involving the sodium transporter OsHKT1;5 in saline tolerant rice (Oryza sativa L.) landraces

Naturally evolved saline tolerant rice landraces found along the coastline of India are a valuable genomic resource to explore the complex, polygenic nature of salinity tolerance. In the present study, a set of 28 genome wide SSR markers, 11 salt responsive genic SSR markers and 8 Saltol QTL linked SSR markers were used to estimate genetic relatedness and population structure within a collection of 47 rice landraces (including a tolerant and 2 sensitive checks) originating from geographically divergent coastal regions of India. All three marker types identified substantial genetic variation among the landraces, as evident from their higher PIC values (0.53 for genomic SSRs, 0.43 for Genic SSRs and 0.59 for Saltol SSRs). The markers RM431, RM484 (Genomic SSRs), OsCAX (D), OsCAX (T) (Genic SSRs) and RM562 (Saltol SSR) were identified as good candidates to be used in breeding programs for improving salinity tolerance in rice. STRUCTURE analysis divided the landraces into five distinct populations, with classification correlating with their geographical locations. Principal coordinate and hierarchical cluster analyses (UPGMA and neighbor joining) are in close agreement with STRUCTURE results. AMOVA analysis indicated a higher magnitude of genetic differentiation within individuals of groups (58%), than among groups (42%). We also report the development and validation of a new Cleavage Amplified Polymorphic Sequence (CAPS) marker (OsHKT1;5V395) that targets a codon in the sodium transporter gene OsHKT1;5 (Saltol/SKC1 locus) that is associated with sodium transport rates in the above rice landraces. The CAPS marker was found to be present in all landraces except in IR29, Kamini, Gheus, Matla 1 and Matla 2. Significant molecular genetic diversity established among the analyzed salt tolerant rice landraces will aid in future association mapping; the CAPS marker, OsHKT1;5V395 can be used to map rice landraces for the presence of the SNP (Single Nucleotide Polymorphism) associated with increased sodium transport rates and concomitant salinity tolerance in rice.

MicroRNA-guided regulation of heat stress response in wheat

Background

With rising global temperature, understanding plants’ adaptation to heat stress has implications in plant breeding. MicroRNAs (miRNAs) are small, non-coding, regulatory RNAs guiding gene expression at the post-transcriptional level. In this study, small RNAs and the degradome (parallel analysis of RNA ends) of leaf tissues collected from control and heat-stressed wheat plants immediately at the end of the stress period and 1 and 4 days later were analysed.

Results

Sequencing of 24 small RNA libraries produced 55.2 M reads while 404 M reads were obtained from the corresponding 24 PARE libraries. From these, 202 miRNAs were ascertained, of which mature miRNA evidence was obtained for 104 and 36 were found to be differentially expressed after heat stress. The PARE analysis identified 589 transcripts targeted by 84 of the ascertained miRNAs. PARE sequencing validated the targets of the conserved members of miRNA156, miR166 and miR393 families as squamosa promoter-binding-like, homeobox leucine-zipper and transport inhibitor responsive proteins, respectively. Heat stress responsive miRNA targeted superoxide dismutases and an array of homeobox leucine-zipper proteins, F-box proteins and protein kinases. Query of miRNA targets to interactome databases revealed a predominant association of stress responses such as signalling, antioxidant activity and ubiquitination to superoxide dismutases, F-box proteins, pentatricopeptide repeat-containing proteins and mitochondrial transcription termination factor-like proteins.

Conclusion

The interlaced data set generated in this study identified and validated heat stress regulated miRNAs and their target genes associated with thermotolerance. Such accurate identification and validation of miRNAs and their target genes are essential to develop novel regulatory gene-based breeding strategies.

Electronic supplementary material

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

Advances in understanding salt tolerance in rice

This review presents a comprehensive overview of the recent research on rice salt tolerance in the areas of genomics, proteomics, metabolomics and chemical genomics. Salinity is one of the major constraints in rice cultivation globally. Traditionally, rice is a glycophyte except for a few genotypes that have been widely used in salinity tolerance breeding of rice. Both seedling and reproductive stages of rice are considered to be the salt-susceptible stages; however, research efforts have been biased towards improving the understanding of seedling-stage salt tolerance. An extensive literature survey indicated that there have been very few attempts to develop reproductive stage-specific salt tolerance in rice probably due to the lack of salt-tolerant phenotypes at the reproductive stage. Recently, the role of DNA methylation, genome duplication and codon usage bias in salinity tolerance of rice have been studied. Furthermore, the study of exogenous salt stress alleviants in rice has opened up another potential avenue for understanding and improving its salt tolerance. There is a need to not only generate additional genomic resources in the form of salt-responsive QTLs and molecular markers and to characterize the genes and their upstream regulatory regions, but also to use them to gain deep insights into the mechanisms useful for developing tolerant varieties. We analysed the genomic locations of diverse salt-responsive genomic resources and found that rice chromosomes 1-6 possess the majority of these salinity-responsive genomic resources. The review presents a comprehensive overview of the recent research on rice salt tolerance in the areas of genomics, proteomics, metabolomics and chemical genomics, which should help in understanding the molecular basis of salinity tolerance and its more effective improvement in rice.

Alfalfa response to heat stress is modulated by microRNA156

Heat stress and extreme temperatures negatively affect plant development by disrupting regular cellular and biochemical functions, ultimately leading to reduced crop production. Alfalfa (Medicago sativa) is an important forage crop grown worldwide as forage for livestock feed. Limiting the effects of abiotic stress by developing alfalfa cultivars that are stress tolerant would help mitigate losses to crop production. Members of the microRNA156 (miR156) family regulate the Squamosa Promoter-Binding Protein-Like (SPL) genes that in turn impact plant growth and development by regulating downstream genes in response to various abiotic stresses. In this study, alfalfa with miR156 overexpression and SPL13 RNAi knockdown show increased tolerance to heat stress (40°C). Transgenic plants show high water potential and increased non-enzymatic antioxidant content under heat stress. Moreover, anthocyanin content and chlorophyll abundance were increased under stress. Expression of some important transcription factors and downstream genes involved in abiotic stress response were altered in miR156-overexpressing genotypes under heat. Taken together, our results demonstrate that the miR156/SPL13 network contributes to improving heat stress tolerance in alfalfa.

miRNomes involved in imparting thermotolerance to crop plants

Thermal stress is one of the challenges to crop plants that negatively impacts crop yield. To overcome this ever-growing problem, utilization of regulatory mechanisms, especially microRNAs (miRNAs), that provide efficient and precise regulation in a targeted manner have been found to play determining roles. Besides their roles in plant growth and development, many recent studies have shown differential regulation of several miRNAs during abiotic stresses including heat stress (HS). Thus, understanding the underlying mechanism of miRNA-mediated gene expression during HS will enable researchers to exploit this regulatory mechanism to address HS responses. This review focuses on the miRNAs and regulatory networks that were involved in physiological, metabolic and morphological adaptations during HS in plant, specifically in crops. Illustrated examples including, the miR156-SPL, miR169-NF-YA5, miR395-APS/AST, miR396-WRKY, etc., have been discussed in specific relation to the crop plants. Further, we have also discussed the available plant miRNA databases and bioinformatics tools useful for miRNA identification and study of their regulatory role in response to HS. Finally, we have briefly discussed the future prospects about the miRNA-related mechanisms of HS for improving thermotolerance in crop plants.

Regulation mechanism of microRNA in plant response to abiotic stress and breeding

microRNAs (miRNAs) in plants are a class of small RNAs consisting of approximately 21-24 nucleotides. The mature miRNA binds to the target mRNA through the formation of a miRNA-induced silencing complex (MIRISC), and cleaves or inhibits translation, thereby achieving negative regulation of the target gene. Based on miRNA plays an important role in regulating plant gene expression, studies on the prediction, identification, function and evolution of plant miRNAs have been carried out. In addition, many researches prove that miRNAs are also involved in many kinds of abiotic and biotic stress, under abiotic stress, plants can express some miRNA, and act on stress-related target genes, which can make plants adapt to stress in physiological response. In this review, the synthetic pathway and mechanism of plant miRNA are briefly described, and we discuss the biological functions and regulatory mechanisms of miRNAs responding to abiotic stresses including low temperature, salt, drought stress and breeding to lay the foundation for further exploring the mechanism of action of miRNAs in stress resistance of plant. And analyze its utilization prospects in plant stress resistance research.

Genetic Variation and Association Analysis of the SSR Markers Linked to the Major Drought-Yield QTLs of Rice

Drought is one of the major abiotic stresses, which hampers the production of rice worldwide. Informative molecular markers are valuable tools for improving the drought tolerance in various varieties of rice. The present study was conducted to evaluate the informative simple sequence repeat (SSR) markers in a diverse set of rice genotypes. The genetic diversity analyses of the 83 studied rice genotypes were performed using 34 SSR markers closely linked to the major quantitative trait loci (QTLs) of grain yield under drought stress (qDTYs). In general, our results indicated high levels of polymorphism. In addition, we screened these rice genotypes at the reproductive stage under both drought stress and nonstressful conditions. The results of the regression analysis demonstrated a significant relationship between 11 SSR marker alleles and the plant paddy weight under stressful conditions. Under the nonstressful conditions, 16 SSR marker alleles showed a significant correlation with the plant paddy weight. Finally, four markers (RM279, RM231, RM166, and RM231) demonstrated a significant association with the plant paddy weight under both stressful and nonstressful conditions. These informative-associated alleles may be useful for improving the crop yield under both drought stress and nonstressful conditions in breeding programs.

Discovery of microRNA-target modules of African rice (Oryza glaberrima) under salinity stress

Oryza glaberrima is the second edible rice in the genus Oryza. It is grown in the African countries. miRNAs are regulatory molecules that are involved in every domains of gene expression including salinity stress response. Although several miRNAs have been reported from various species of Oryza, yet none of them are from this species. Salt treated (200 mM NaCl for 48 h) and control smallRNA libraries of RAM-100, a salt tolerant genotype, each with 2 replications generated 150 conserve and 348 novel miRNAs. We also used smallRNAseq data of NCBI of O. glaberrima to discover additional 246 known miRNAs. Totally, 29 known and 32 novel miRNAs were differentially regulated under salinity stress. Gene ontology and KEGG analysis indicated several targets were involved in vital biological pathways of salinity stress tolerance. Expression of selected miRNAs as indicated by Illumina data were found to be coherent with real time-PCR analysis. However, target gene expression was inversely correlated with their corresponding miRNAs. Finally based upon present results as well as existing knowledge of literature, we proposed the miRNA-target modules that were induced by salinity stress. Therefore, the present findings provide valuable information about miRNA-target networks in salinity adaption of O. glaberrima.

Regulation of low temperature stress in plants by microRNAs

Low temperature is one of the most common environmental stresses that seriously affect the growth and development of plants. However, plants have the plasticity in their defence mechanisms enabling them to tolerate and, sometimes, even survive adverse environmental conditions. MicroRNAs (miRNAs) are small non-coding RNAs, approximately 18-24 nucleotides in length, and are being increasingly recognized as regulators of gene expression at the post-transcriptional level and have the ability to influence a broad range of biological processes. There is growing evidence in the literature that reprogramming of gene expression mediated through miRNAs is a major defence mechanism in plants enabling them to respond to stresses. To date, numerous studies have established the importance of miRNA-based regulation of gene expression under low temperature stress. Individual miRNAs can modulate the expression of multiple mRNA targets, and, therefore, the manipulation of a single miRNA has the potential to affect multiple biological processes. Numerous functional studies have attempted to identify the miRNA-target interactions and have elaborated the role of several miRNAs in cold-stress regulation. This review summarizes the current understanding of miRNA-mediated modulation of the expression of key genes as well as genetic and regulatory pathways, involved in low temperature stress responses in plants.

Genome-Wide Development of MicroRNA-Based SSR Markers in Medicago truncatula with Their Transferability Analysis and Utilization in Related Legume Species

Microsatellite (simple sequence repeats, SSRs) marker is one of the most widely used markers in marker-assisted breeding. As one type of functional markers, MicroRNA-based SSR (miRNA-SSR) markers have been exploited mainly in animals, but the development and characterization of miRNA-SSR markers in plants are still limited. In the present study, miRNA-SSR markers for Medicago truncatula (M. truncatula) were developed and their cross-species transferability in six leguminous species was evaluated. A total of 169 primer pairs were successfully designed from 130 M. truncatula miRNA genes, the majority of which were mononucleotide repeats (70.41%), followed by dinucleotide repeats (14.20%), compound repeats (11.24%) and trinucleotide repeats (4.14%). Functional classification of SSR-containing miRNA genes showed that all targets could be grouped into three Gene Ontology (GO) categories: 17 in biological process, 11 in molecular function, and 14 in cellular component. The miRNA-SSR markers showed high transferability in other six leguminous species, ranged from 74.56% to 90.53%. Furthermore, 25 Mt-miRNA-SSR markers were used to evaluate polymorphisms in 20 alfalfa accessions, and the polymorphism information content (PIC) values ranged from 0.39 to 0.89 with an average of 0.71, the allele number per marker varied from 3 to 18 with an average of 7.88, indicating a high level of informativeness. The present study is the first time developed and characterized of M. truncatula miRNA-SSRs and demonstrated their utility in transferability, these novel markers will be valuable for genetic diversity analysis, marker-assisted selection and genotyping in leguminous species.

Comprehensive survey and evolutionary analysis of genome-wide miRNA genes from ten diploid Oryza species

BACKGROUND

MicroRNAs (miRNAs) are non-coding RNAs that play versatile roles in post-transcriptional gene regulation. Although much is known about their biogenesis, and gene regulation very little is known about their evolutionary relation among the closely related species.

RESULT

All the orthologous miRNA genes of Oryza sativa (japonica) from 10 different Oryza species were identified, and the evolutionary changes among these genes were analysed. Significant differences in the expansion of miRNA gene families were observed across the Oryza species. Analysis of the nucleotide substitution rates indicated that the mature sequences show the least substitution rates among the different regions of miRNA genes, and also show a very much less substitution rates as compared to that of all protein-coding genes across the Oryza species. Evolution of miRNA genes was also found to be contributed by transposons. A non-neutral selection was observed at 80 different miRNA loci across Oryza species which were estimated to have lost ~87% of the sequence diversity during the domestication. The phylogenetic analysis revealed that O. longistaminata diverged first among the AA-genomes, whereas O. brachyantha and O. punctata appeared as the eminent out-groups. The miR1861 family organised into nine distinct compact clusters in the studied Oryza species except O. brachyantha. Further, the expression analysis showed that 11 salt-responsive miRNAs were differentially regulated between O. coarctata and O. glaberrima.

CONCLUSION

Our study provides the evolutionary dynamics in the miRNA genes of 10 different Oryza species which will support more investigations about the structural and functional organization of miRNA genes of Oryza species.

Genome-wide analysis of DUF221 domain-containing gene family in Oryza species and identification of its salinity stress-responsive members in rice

DUF221 domain-containing genes (DDP genes) play important roles in developmental biology, hormone signalling transduction, and responses to abiotic stress. Therefore to understand their structural and evolutionary relationship, we did a genome-wide analysis of this important gene family in rice. Further, through comparative genomics, DDP genes from Oryza sativa subsp. (indica), nine different wild species of rice and Arabidopsis were also identified. We also found an expansion of the DDP gene families in rice and Arabidopsis which is due to the segmental duplication events in some of the gene family members. In general, a highly purifying selection was found acting on all the deduced paralogous and orthologous DDP gene pairs. The data from microarray and subsequent qRT-PCR analysis revealed that although several OsDDPs were differentially regulated under salinity stress, yet OsDDP6 was upregulated at all the developmental stages in salt tolerant rice genotype, FL478. Interestingly, OsDDP6 was found to be involved in proline metabolism pathway as indicated by protein network analysis. The diverse gene structures, varied transmembrane topologies and the differential expression patterns implied the functional diversity in DDP genes. Therefore, the comprehensive evolutionary analysis of DDP genes from different Oryza species and Arabidopsis performed in this study will provide the basis for further functional validation studies vis-à-vis DDP genes of rice and other plant species.

Extrapolative microRNA precursor based SSR mining from tea EST database in respect to agronomic traits

Tea (Camellia sinensis, (L.) Kuntze) is considered as most popular drink across the world and it is widely consumed beverage for its several health-benefit characteristics. These positive traits primarily rely on its regulatory networks of different metabolic pathways. Development of microsatellite markers from the conserved genomic regions are being worthwhile for reviewing the genetic diversity of closely related species or self-pollinated species. Although several SSR markers have been reported, in tea, the trait-specific Simple Sequence Repeat (SSR) markers, leading to be useful in marker assisted breeding technique, are yet to be identified. Micro RNAs are short, non-coding RNA molecules, involved in post transcriptional mode of gene regulation and thus effects on related phenotype. Present study deals with identification of the microsatellite motifs within the reported and predicted miRNA precursors that are effectively followed by designing of primers from SSR flanking regions in order to PCR validation. In addition to the earlier reports, two new miRNAs are predicting here from tea expressed tag sequence database. Furthermore, 18 SSR motifs are found to be in 13 of all 33 predicted miRNAs. Trinucleotide motifs are most abundant among all followed by dinucleotides. Since, miRNA based SSR markers are evidenced to have significant role on genetic fingerprinting study, these outcomes would pave the way in developing novel markers for tagging tea specific agronomic traits as well as substantiating non-conventional breeding program.

Identification of microRNAs in Response to Drought in Common Wild Rice (Oryza rufipogon Griff.) Shoots and Roots

Background

Drought is the most important factor that limits rice production in drought-prone environments. Plant microRNAs (miRNAs) are involved in biotic and abiotic stress responses. Common wild rice (Oryza rufipogon Griff.) contains abundant drought-resistant genes, which provide an opportunity to explore these excellent resources as contributors to improve rice resistance, productivity, and quality.

Results

In this study, we constructed four small RNA libraries, called CL and CR from PEG6000-free samples and DL and DR from PEG6000-treated samples, where ‘R’ indicates the root tissue and ‘L’ indicates the shoot tissue. A total of 200 miRNAs were identified to be differentially expressed under the drought-treated conditions (16% PEG6000 for 24 h), and the changes in the miRNA expression profile of the shoot were distinct from those of the root. At the miRNA level, 77 known miRNAs, which belong to 23 families, including 40 up-regulated and 37 down-regulated in the shoot, and 85 known miRNAs in 46 families, including 65 up-regulated and 20 down-regulated in the root, were identified as differentially expressed. In addition, we predicted 26 new miRNA candidates from the shoot and 43 from the root that were differentially expressed during the drought stress. The quantitative real-time PCR analysis results were consistent with high-throughput sequencing data. Moreover, 88 miRNAs that were differentially-expressed were predicted to match with 197 targets for drought-stress.

Conclusion

Our results suggest that the miRNAs of O. rufipogon are responsive to drought stress. The differentially expressed miRNAs that are tissue-specific under drought conditions could play different roles in the regulation of the auxin pathway, the flowering pathway, the drought pathway, and lateral root formation. Thus, the present study provides an account of tissue-specific miRNAs that are involved in the drought adaption of O. rufipogon.

Rice Improvement Through Genome-Based Functional Analysis and Molecular Breeding in India

Rice is one of the main pillars of food security in India. Its improvement for higher yield in sustainable agriculture system is also vital to provide energy and nutritional needs of growing world population, expected to reach more than 9 billion by 2050. The high quality genome sequence of rice has provided a rich resource to mine information about diversity of genes and alleles which can contribute to improvement of useful agronomic traits. Defining the function of each gene and regulatory element of rice remains a challenge for the rice community in the coming years. Subsequent to participation in IRGSP, India has continued to contribute in the areas of diversity analysis, transcriptomics, functional genomics, marker development, QTL mapping and molecular breeding, through national and multi-national research programs. These efforts have helped generate resources for rice improvement, some of which have already been deployed to mitigate loss due to environmental stress and pathogens. With renewed efforts, Indian researchers are making new strides, along with the international scientific community, in both basic research and realization of its translational impact.

A Genome-Wide Survey of the Microsatellite Content of the Globe Artichoke Genome and the Development of a Web-Based Database

The recently acquired genome sequence of globe artichoke (Cynara cardunculus var. scolymus) has been used to catalog the genome’s content of simple sequence repeat (SSR) markers. More than 177,000 perfect SSRs were revealed, equivalent to an overall density across the genome of 244.5 SSRs/Mbp, but some 224,000 imperfect SSRs were also identified. About 21% of these SSRs were complex (two stretches of repeats separated by <100 nt). Some 73% of the SSRs were composed of dinucleotide motifs. The SSRs were categorized for the numbers of repeats present, their overall length and were allocated to their linkage group. A total of 4,761 perfect and 6,583 imperfect SSRs were present in 3,781 genes (14.11% of the total), corresponding to an overall density across the gene space of 32,5 and 44,9 SSRs/Mbp for perfect and imperfect motifs, respectively. A putative function has been assigned, using the gene ontology approach, to the set of genes harboring at least one SSR. The same search parameters were applied to reveal the SSR content of 14 other plant species for which genome sequence is available. Certain species-specific SSR motifs were identified, along with a hexa-nucleotide motif shared only with the other two Compositae species (sunflower (Helianthus annuus) and horseweed (Conyza canadensis)) included in the study. Finally, a database, called “Cynara cardunculus MicroSatellite DataBase” (CyMSatDB) was developed to provide a searchable interface to the SSR data. CyMSatDB facilitates the retrieval of SSR markers, as well as suggested forward and reverse primers, on the basis of genomic location, genomic vs genic context, perfect vs imperfect repeat, motif type, motif sequence and repeat number. The SSR markers were validated via an in silico based PCR analysis adopting two available assembled transcriptomes, derived from contrasting globe artichoke accessions, as templates.

The genome sequence of the outbreeding globe artichoke constructed de novo incorporating a phase-aware low-pass sequencing strategy of F1 progeny

Globe artichoke (Cynara cardunculus var. scolymus) is an out-crossing, perennial, multi-use crop species that is grown worldwide and belongs to the Compositae, one of the most successful Angiosperm families. We describe the first genome sequence of globe artichoke. The assembly, comprising of 13,588 scaffolds covering 725 of the 1,084 Mb genome, was generated using ~133-fold Illumina sequencing data and encodes 26,889 predicted genes. Re-sequencing (30×) of globe artichoke and cultivated cardoon (C. cardunculus var. altilis) parental genotypes and low-coverage (0.5 to 1×) genotyping-by-sequencing of 163 F₁ individuals resulted in 73% of the assembled genome being anchored in 2,178 genetic bins ordered along 17 chromosomal pseudomolecules. This was achieved using a novel pipeline, SOILoCo (Scaffold Ordering by Imputation with Low Coverage), to detect heterozygous regions and assign parental haplotypes with low sequencing read depth and of unknown phase. SOILoCo provides a powerful tool for de novo genome analysis of outcrossing species. Our data will enable genome-scale analyses of evolutionary processes among crops, weeds, and wild species within and beyond the Compositae, and will facilitate the identification of economically important genes from related species.