Identification of an miRNA candidate reflects the possible significance of transcribed microsatellites in the hairpin precursors of black pepper

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.

Small RNA sequencing reveals various microRNAs involved in piperine biosynthesis in black pepper (Piper nigrum L.)

Background

Black pepper (Piper nigrum L.), an important and long-cultivated spice crop, is native to South India and grown in the tropics. Piperine is the main pungent and bioactive alkaloid in the berries of black pepper, but the molecular mechanism for piperine biosynthesis has not been determined. MicroRNAs (miRNAs), which are classical endogenous noncoding small RNAs, play important roles in regulating secondary metabolism in many species, but less is known regarding black pepper or piperine biosynthesis.

Results

To dissect the functions of miRNAs in secondary metabolism especially in piperine biosynthesis, 110 known miRNAs, 18 novel miRNAs and 1007 individual targets were identified from different tissues of black pepper by small RNA sequencing. qRT-PCR and 5′-RLM-RACE experiments were conducted to validate the reliability of the sequencing data and predicted targets. We found 3 miRNAs along with their targets including miR166-4CL, miR396-PER and miR397-CCR modules that are involved in piperine biosynthesis.

Conclusion

MiRNA regulation of secondary metabolism is a common phenomenon in plants. Our study revealed new miRNAs that regulate piperine biosynthesis, which are special alkaloids in the piper genus, and they might be useful for future piperine genetic improvement of black pepper.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08154-4.

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.

Elucidating the Functional Role of Predicted miRNAs in Post- Transcriptional Gene Regulation Along with Symbiosis in Medicago truncatula

Background:

microRNAs are small non-coding RNAs which inhibit translational and post-transcriptional processes whereas long non-coding RNAs are found to regulate both transcriptional and post-transcriptional gene expression. Medicago truncatula is a well-known model plant for studying legume biology and is also used as a forage crop. In spite of its importance in nitrogen fixation and soil fertility improvement, little information is available about Medicago non-coding RNAs that play important role in symbiosis.

Objective:

In this study we have tried to understand the role of Medicago ncRNAs in symbiosis and regulation of transcription factors.

Methods:

We have identified novel miRNAs by computational methods considering various parameters like length, MFEI, AU content, SSR signatures and tried to establish an interaction model with their targets obtained through psRNATarget server.

Results:

149 novel miRNAs are predicted along with their 770 target proteins. We have also shown that 51 of these novel miRNAs are targeting 282 lncRNAs.

Conclusion:

In this study role of Medicago miRNAs in the regulation of various transcription factors are elucidated. Knowledge gained from this study will have a positive impact on the nitrogen fixing ability of this important model plant, which in turn will improve the soil fertility.

Network analysis of noncoding RNAs in pepper provides insights into fruit ripening control

Pepper is an important vegetable worldwide and is a model plant for nonclimacteric fleshy fruit ripening. Drastic visual changes and internal biochemical alterations are involved in fruit coloration, flavor, texture, aroma, and palatability to animals during the pepper fruit ripening process. To explore the regulation of bell pepper fruit ripening by noncoding RNAs (ncRNAs), we examined their expression profiles; 43 microRNAs (miRNAs), 125 circular RNAs (circRNAs), 366 long noncoding RNAs (lncRNAs), and 3266 messenger RNAs (mRNAs) were differentially expressed (DE) in mature green and red ripe fruit. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed that the targets of the DE ncRNAs and DE mRNAs included several kinds of transcription factors (TFs) (ERF, bHLH, WRKY, MYB, NAC, bZIP, and ARF), enzymes involved in cell wall metabolism (beta-galactosidase, beta-glucosidase, beta-amylase, chitinase, pectate lyase (PL), pectinesterase (PE) and polygalacturonase (PG)), enzymes involved in fruit color accumulation (bifunctional 15-cis-phytoene synthase, 9-cis-epoxycarotenoid dioxygenase, beta-carotene hydroxylase and carotene epsilon-monooxygenase), enzymes associated with fruit flavor and aroma (glutamate-1-semialdehyde 2,1-aminomutase, anthocyanin 5-aromatic acyltransferase, and eugenol synthase 1) and enzymes involved in the production of ethylene (ET) (ACO1/ACO4) as well as other plant hormones such as abscisic acid (ABA), auxin (IAA), and gibberellic acid (GA). Based on accumulation profiles, a network of ncRNAs and mRNAs associated with bell pepper fruit ripening was developed that provides a foundation for further developing a more refined understanding of the molecular biology of fruit ripening.

Analysis of the Coding and Non-Coding RNA Transcriptomes in Response to Bell Pepper Chilling

Increasing evidence suggests that long non-coding RNAs (lncRNAs), circular RNAs (circRNAs), and microRNAs (miRNAs) have roles during biotic and abiotic stress, though their exact contributions remain unclear. To explore their biological functions in response to chilling in bell pepper, we examined their accumulation profiles by deep sequencing and identified 380 lncRNAs, 36 circRNAs, 18 miRNAs, and 4128 differentially expressed mRNAs in the chilled versus the non-chilled fruit. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed differentially expressed genes and putative ncRNA targets, including transcription factors of multiple classes, such as myeloblastosis (MYB), basic helix-loop-helix (bHLH), and ethylene response factor (ERF) transcription factors (TFs), enzymes involved in bio-oxidation and oxidative phosphorylation (serine/threonine-protein kinase, polyphenol oxidase, catalase, peroxidase, lipoxygenase, and ATPase), and cell wall metabolism-related enzymes (beta-galactosidase, pectate lyase, pectinesterase, and polygalacturonase). On the basis of the accumulation profiles, a network of putatively interacting RNAs associated with bell pepper chilling was developed, which pointed to ncRNAs that could provide the foundation for further developing a more refined understanding of the molecular response to chilling injury.

Abiotic Stress Responsive miRNA-Target Network and Related Markers (SNP, SSR) in Brassica juncea

Abiotic stress is one of the major factors responsible for huge yield loss in crop plants. MicroRNAs play a key role in adaptive responses of plants under abiotic stress conditions through post-transcriptional gene regulations. In present study, 95 potential miRNAs were predicted in Brassica juncea using comparative genomics approach. It was noted that these miRNAs, target several transcription factors (TFs), transporter family proteins, signaling related genes, and protease encoding genes. Nineteen distinct miRNA-target regulatory networks were observed with significant involvement in regulation of transcription, response to stimulus, hormone and auxin mediated signaling pathway related gene ontology (GO) term. The sucrose-starch metabolism and pentose-gluconate interconversion pathways were found significantly enriched for these target genes. Molecular markers such as Simple Sequence Repeats (SSR) and Single Nucleotide Polymorphism (SNPs) were identified on miRNAs (miR-SSRs and miR-SNPs) and their target genes in B. juncea. Notably, one of the miR-SNP (C/T) was found at the 5^th position on mature region of miR2926. This C/T transition led to the distorted and unstable hairpin structure of miR2926, consequently complete loss of target function. Hence, findings from this study will lay a foundation for marker assisted breeding for abiotic stress tolerant varieties of B. juncea.

Transfer RNA Derived Small RNAs Targeting Defense Responsive Genes Are Induced during Phytophthora capsici Infection in Black Pepper (Piper nigrum L.)

Small RNAs derived from transfer RNAs were recently assigned as potential gene regulatory candidates for various stress responses in eukaryotes. In this study, we report on the cloning and identification of tRNA derived small RNAs from black pepper plants in response to the infection of the quick wilt pathogen, Phytophthora capsici. 5'tRFs cloned from black pepper were validated as highly expressed during P. capsici infection. A high-throughput systematic analysis of the small RNAome (sRNAome) revealed the predominance of 5'tRFs in the infected leaf and root. The abundance of 5'tRFs in the sRNAome and the defense responsive genes as their potential targets indicated their regulatory role during stress response in black pepper. The 5'Ala(CGC) tRF mediated cleavage was experimentally mapped at the tRF binding sites on the mRNA targets of Non-expresser of pathogenesis related protein (NPR1), which was down-regulated during pathogen infection. Comparative sRNAome further demonstrated sequence conservation of 5'Ala tRFs across the angiosperm plant groups, and many important genes in the defense response were identified in silico as their potential targets. Our findings uncovered the diversity, differential expression and stress responsive functional role of tRNA-derived small RNAs during Phytophthora infection in black pepper.

Unravelling the complexity of microRNA-mediated gene regulation in black pepper (Piper nigrum L.) using high-throughput small RNA profiling

Analysis of high-throughput small RNA deep sequencing data, in combination with black pepper transcriptome sequences revealed microRNA-mediated gene regulation in black pepper ( Piper nigrum L.). Black pepper is an important spice crop and its berries are used worldwide as a natural food additive that contributes unique flavour to foods. In the present study to characterize microRNAs from black pepper, we generated a small RNA library from black pepper leaf and sequenced it by Illumina high-throughput sequencing technology. MicroRNAs belonging to a total of 303 conserved miRNA families were identified from the sRNAome data. Subsequent analysis from recently sequenced black pepper transcriptome confirmed precursor sequences of 50 conserved miRNAs and four potential novel miRNA candidates. Stem-loop qRT-PCR experiments demonstrated differential expression of eight conserved miRNAs in black pepper. Computational analysis of targets of the miRNAs showed 223 potential black pepper unigene targets that encode diverse transcription factors and enzymes involved in plant development, disease resistance, metabolic and signalling pathways. RLM-RACE experiments further mapped miRNA-mediated cleavage at five of the mRNA targets. In addition, miRNA isoforms corresponding to 18 miRNA families were also identified from black pepper. This study presents the first large-scale identification of microRNAs from black pepper and provides the foundation for the future studies of miRNA-mediated gene regulation of stress responses and diverse metabolic processes in black pepper.

Computational prediction of miRNAs and their targets in Phaseolus vulgaris using simple sequence repeat signatures

BACKGROUND

MicroRNAs (miRNAs) are endogenous, noncoding, short RNAs directly involved in regulating gene expression at the post-transcriptional level. In spite of immense importance, limited information of P. vulgaris miRNAs and their expression patterns prompted us to identify new miRNAs in P. vulgaris by computational methods. Besides conventional approaches, we have used the simple sequence repeat (SSR) signatures as one of the prediction parameter. Moreover, for all other parameters including normalized Shannon entropy, normalized base pairing index and normalized base-pair distance, instead of taking a fixed cut-off value, we have used 99% probability range derived from the available data.

RESULTS

We have identified 208 mature miRNAs in P. vulgaris belonging to 118 families, of which 201 are novel. 97 of the predicted miRNAs in P. vulgaris were validated with the sequencing data obtained from the small RNA sequencing of P. vulgaris. Randomly selected predicted miRNAs were also validated using qRT-PCR. A total of 1305 target sequences were identified for 130 predicted miRNAs. Using 80% sequence identity cut-off, proteins coded by 563 targets were identified. The computational method developed in this study was also validated by predicting 229 miRNAs of A. thaliana and 462 miRNAs of G. max, of which 213 for A. thaliana and 397 for G. max are existing in miRBase 20.

CONCLUSIONS

There is no universal SSR that is conserved among all precursors of Viridiplantae, but conserved SSR exists within a miRNA family and is used as a signature in our prediction method. Prediction of known miRNAs of A. thaliana and G. max validates the accuracy of our method. Our findings will contribute to the present knowledge of miRNAs and their targets in P. vulgaris. This computational method can be applied to any species of Viridiplantae for the successful prediction of miRNAs and their targets.

Microsatellite resources of Eucalyptus: current status and future perspectives

Eucalyptus is the premier paper pulp, short rotation plantation species grown all over the world. Genetic improvement programs integrating molecular marker tools are in progress in many parts of the globe to increase the productivity. Whole genome sequence and expressed sequence tags (ESTs) of the eucalypts paved way for introduction of molecular genetics and breeding in this genus. Different molecular characterization approaches have been used simultaneously in eucalypts, however, microsatellites or simple sequence repeats (SSRs) with their prolific characteristics could occupy a special niche in Eucalyptus genetic improvement. Further, highly informative SSRs were used for the clonal identity, genetic fidelity and in certification of breeder's rights. Eucalyptus genetic linkage maps generated with microsatellite loci were used successfully to identify quantitative trait loci (QTLs) for various economically important traits. Progressively more numbers of microsatellites are being linked to genes associated with adaptive and functional variations, therefore making their utility broader in genetic applications. Availability of common SSR markers across the species provides an opportunity to validate the expression of QTLs across variable genetic backgrounds and accurately compare the position of QTLs in other species. Recent evidences suggest that the presence of SSRs in micro RNAs of plant species play a role in the quantitative trait expression. Similar studies in eucalypts may provide new insights into the genetic architecture of transcript-level variations and post transcriptional gene regulation. This review on eucalypts microsatellites, highlights the availability and characteristics of genomic and eSSRs and their potential in genetic analysis of natural and breeding populations and also discusses the future prospects in population genetics and marker assisted selection.

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

Salinity is an important abiotic stress that affects agricultural production and productivity. It is a complex trait that is regulated by different molecular mechanisms. miRNAs are non-coding RNAs which are highly conserved and regulate gene expression. Simple sequence repeats (SSRs) are robust molecular markers for studying genetic diversity. Although several SSR markers are available now, challenge remains to identify the trait-specific SSRs which can be used for marker assisted breeding. In order to understand the genetic diversity of salt responsive-miRNA genes in rice, SSR markers were mined from 130 members of salt-responsive miRNA genes of rice and validated among the contrasting panels of tolerant as well as susceptible rice genotypes, each with 12 genotypes. Although 12 miR-SSRs were found to be polymorphic, only miR172b-SSR was able to differentiate the tolerant and susceptible genotypes in 2 different groups. It had also been found that miRNA genes were more diverse in susceptible genotypes than the tolerant one (as indicated by polymorphic index content) which might interfere to form the stem-loop structure of premature miRNA and their subsequent synthesis in susceptible genotypes. Thus, we concluded that length variations of the repeats in salt responsive miRNA genes may be responsible for a possible sensitivity to salinity adaptation. This is the first report of characterization of trait specific miRNA derived SSRs in plants.

Revisiting an important component of plant genomes: microsatellites

Microsatellites are some of the most highly variable repetitive DNA tracts in genomes. Few studies focus on whether the characteristic instability of microsatellites is linked to phenotypic effects in plants. We summarise recent data to investigate how microsatellite variations affect gene expression and hence phenotype. We discuss how the basic characteristics of microsatellites may contribute to phenotypic effects. In summary, microsatellites in plants are universal and highly mutable, they coexist and coevolve with transposable elements, and are under selective pressure. The number of motif nucleotides, the type of motif and transposon activity all contribute to the nonrandom generation and decay of microsatellites, and to conservation and distribution biases. Although microsatellites are generated by accident, they mature through responses to environmental change before final decay. This process is mediated by organism adjustment mechanisms, which maintain a balance between birth versus death and growth versus decay in microsatellites. Close relationships also exist between the physical structure, variation and functionality of microsatellites: in most plant species, sequences containing microsatellites are associated with catalytic activity and binding functions, are expressed in the membrane and organelles, and participate in the developmental and metabolic processes. Microsatellites contribute to genome structure and functional plasticity, and may be considered to promote species evolution in plants in response to environmental changes. In conclusion, the generation, loss, functionality and evolution of microsatellites can be related to plant gene expression and functional alterations. The effect of microsatellites on phenotypic variation may be as significant in plants as it is in animals.

De novo transcriptome sequencing reveals a considerable bias in the incidence of simple sequence repeats towards the downstream of 'Pre-miRNAs' of black pepper

Next generation sequencing has an advantageon transformational development of species with limited available sequence data as it helps to decode the genome and transcriptome. We carried out the de novo sequencing using illuminaHiSeq™ 2000 to generate the first leaf transcriptome of black pepper (Piper nigrum L.), an important spice variety native to South India and also grown in other tropical regions. Despite the economic and biochemical importance of pepper, a scientifically rigorous study at the molecular level is far from complete due to lack of sufficient sequence information and cytological complexity of its genome. The 55 million raw reads obtained, when assembled using Trinity program generated 2,23,386 contigs and 1,28,157 unigenes. Reports suggest that the repeat-rich genomic regions give rise to small non-coding functional RNAs. MicroRNAs (miRNAs) are the most abundant type of non-coding regulatory RNAs. In spite of the widespread research on miRNAs, little is known about the hair-pin precursors of miRNAs bearing Simple Sequence Repeats (SSRs). We used the array of transcripts generated, for the in silico prediction and detection of ‘43 pre-miRNA candidates bearing different types of SSR motifs’. The analysis identified 3913 different types of SSR motifs with an average of one SSR per 3.04 MB of thetranscriptome. About 0.033% of the transcriptome constituted ‘pre-miRNA candidates bearing SSRs’. The abundance, type and distribution of SSR motifs studied across the hair-pin miRNA precursors, showed a significant bias in the position of SSRs towards the downstream of predicted ‘pre-miRNA candidates’. The catalogue of transcripts identified, together with the demonstration of reliable existence of SSRs in the miRNA precursors, permits future opportunities for understanding the genetic mechanism of black pepper and likely functions of ‘tandem repeats’ in miRNAs.