Next Article in Journal
Comparing Hot and Cold Loading in an Integrated Biomass Recovery Operation
Next Article in Special Issue
Transcriptome Sequencing and Differential Expression Analysis Reveal Molecular Mechanisms for Starch Accumulation in Chestnut
Previous Article in Journal
Wood Density Determination by Drilling Chips Extraction in Ten Softwood and Hardwood Species
Previous Article in Special Issue
Update of Genetic Linkage Map and QTL Analysis for Growth Traits in Eucommia ulmoides Oliver
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification and in Silico Characterization of Novel and Conserved MicroRNAs in Methyl Jasmonate-Stimulated Scots Pine (Pinus sylvestris L.) Needles

1
Genetic Resource Centre, Latvian State Forest Research Institute “Silava”, Riga street 111, Salaspils LV-2169, Latvia
2
Norwegian Institute of Bioeconomy Research, Postboks 115, Ås NO-1431, Norway
*
Author to whom correspondence should be addressed.
Forests 2020, 11(4), 384; https://doi.org/10.3390/f11040384
Submission received: 5 March 2020 / Revised: 23 March 2020 / Accepted: 26 March 2020 / Published: 28 March 2020
(This article belongs to the Special Issue Forest Genetics and Tree Improvement)

Abstract

:
MicroRNAs (miRNAs) are non-protein coding RNAs of ~20–24 nucleotides in length that play an important role in many biological and metabolic processes, including the regulation of gene expression, plant growth and developmental processes, as well as responses to stress and pathogens. The aim of this study was to identify and characterize novel and conserved microRNAs expressed in methyl jasmonate-treated Scots pine needles. In addition, potential precursor sequences and target genes of the identified miRNAs were determined by alignment to the Pinus unigene set. Potential precursor sequences were identified using the miRAtool, conserved miRNA precursors were also tested for the ability to form the required stem-loop structure, and the minimal folding free energy indexes were calculated. By comparison with miRBase, 4975 annotated sequences were identified and assigned to 173 miRNA groups, belonging to a total of 60 conserved miRNA families. A total of 1029 potential novel miRNAs, grouped into 34 families were found, and 46 predicted precursor sequences were identified. A total of 136 potential target genes targeted by 28 families were identified. The majority of previously reported highly conserved plant miRNAs were identified in this study, as well as some conserved miRNAs previously reported to be monocot specific. No conserved dicot-specific miRNAs were identified. A number of potential gymnosperm or conifer specific miRNAs were found, shared among a range of conifer species.

1. Introduction

Scots pine is a long-lived organism with a wide distribution in the northern hemisphere. This requires an adaption to a broad range of growing and environmental conditions—likely facilitated by substantial phenotypic plasticity and epigenetic variation [1]. Plants have developed a range of epigenetic mechanisms to deal with biotic and abiotic stress, including DNA methylation, histone modification and expression of various non-coding RNAs (including microRNAs), which influence gene expression and regulation [2]. MicroRNAs (miRNAs) are a class of non-protein coding small RNAs (sRNAs) of ~20–24 nucleotides in length that play an important role in a variety of biological and metabolic processes, primarily through the coordinated action on the post-transcriptional control of gene expression [3]. While miRNAs control essential aspects of plant growth and development, miRNA and sRNA expression are similarly important in the responses to the challenges of stress and pathogens [4,5,6,7]. In plants, miRNA studies have been mainly concentrated in angiosperms and relatively few reports have been made in conifers.
MiRNAs are initially synthesized as primary miRNA transcripts (pri-RNAs), which can be up to several hundred nucleotides in length and contain at least one hairpin stem loop structure. A range of enzymes process these transcripts (DICER-LIKE 1 enzyme, HASTY and others) to generate miRNA precursors (pre-miRNAs) and subsequently produce the mature miRNA molecules [5,8]. In addition, alternative miRNA processing pathways have been described in plant species that involve DICER independent pathways and the processing of other non-coding RNAs [9]. These include miRNAs that are produced from an intron of a protein-coding gene by DICER-independent pre-mRNA splicing machinery (miRtrons), Argonaute RISC catalytic subunit 2 (AGO2) processing, and origins from various other noncoding RNA genes. Plants produce many distinct types of DCL/AGO-associated regulatory sRNAs, from which miRNAs, phased siRNAs and heterochromatic siRNAs are three of the major types of small plant RNAs [10,11].
miRNA precursors can produce both precise excision or distinct miRNA variants known as isoforms (isomiRs). IsomiRs from the same miRNA arm typically differ at their 5′ end, 3′ end, or both, thus abundant isoforms are frequently annotated alongside a canonical miRNA sequence found in databases. IsomiRs in plants can originate from imprecise cleavage by DCL1 (templated isomiRs), which generates variants that show complementarity to their pre-miRNA. Alternatively, isomiRs can be spawned by post-transcriptional modification due to the addition or removal of specific nucleotides to miRNA ends (non-templated isomiRs) [12]. Non-templated modifications most frequently occur at the 3′ end, but less frequently at the 5’ end, and these modifications can influence both the miRNA stability and the efficiency of target repression) [13,14]. One of the mechanisms that can increase the diversity of miRNA action is seed shifting, where the dominant mature miRNA isoform is shifted in complementarity relative to its target by one to several nucleotides in the 5’ or 3’ direction relative to its original position. Therefore, modifications at the 5′ end of the mature product could relocate the seed position and thus change in the seed sequence (called “seed shifting”), thereby altering mRNA target recognition and function [15].
While there are relatively few conifer sequences deposited in specific miRNA databases (e.g., miRBase), there have been a number of reports published about miRNA studies in conifer species. These species include lodgepole pine (Pinus contorta) [16], Eastern white pine (Pinus strobus) [17], larch (Larix leptolepsis) [18,19,20,21,22], Norway spruce (Picea abies) [23,24,25,26], Chinese fir (Cunninghamia lanceolata) [27,28], Chinese yew (Taxus chinensis) [29], sequoia (Sequoia sempervirens) [30] and others. However, there are no publications to date about miRNA studies in Scots pine (Pinus sylvestris L). The majority of miRNA studies in gymnosperms have investigated various developmental stages or plant tissues, e.g., the expression of small RNAs (sRNAs) in Sequoia sempervirens during phase changes, specifically in the juvenile, adult and in vitro propagated plants [30], expression patterns of conserved miRNAs from mature and germinated pollen of Pinus taeda [31], and miRNAs in zygotic embryos and female gametophytes of Pinus taeda [32]. Yakovlev et al. [26] reported on the expression of miRNAs in Norway spruce seedlings derived from plants regenerated after embryogenesis in a cold and warm environment. In addition to developmental processes, miRNAs are involved in regulation of defense responses. Differentially expressed miRNAs target NB-LRR genes in the bark of Norway spruce are produced in response to inoculation with Ceratocystis polonica and wounding [24], while in Pinus taeda miRNAs were also identified in the xylem during the process of fusiform rust gall development [33].
miRNAs are known to modulate the expression of genes during plant defense and the effect of methyl jasmonate on miRNA expression has also been investigated. The action of miRNAs on upon the biosynthesis and perception pathway is of particular importance in trees since jasmonic acid (JA) and its methyl ester, methyl jasmonate (MeJ), are plant signaling molecules that broadly affect gene expression impacting both plant growth and development as well as the response to pathogen attack, wounding and plant responses to abiotic and biotic stresses [34,35]. MeJ, which is synthesized from linolenic acid, and is one of the few plant compounds that are effective at low vapour concentrations. Jasmonates include jasmonic acid, its derivatives and conjugates; the jasmonates and in particular, the active hormone jasmonoyl-isoleucine is known to regulate defenses against chewing herbivores [36]. Remarkably MeJ has been reported to be involved in defense priming, and the induction of anatomical and chemical changes such as the formation of traumatic resin ducts in the xylem, and the synthesis of phenolic and alkaloid compounds in many conifer species [37,38,39]. The interaction between miRNA-controlled gene expression and MeJ biosynthesis and perception is less known in gymnosperm species, but in Taxus chinensis cells, marked changes in miRNA profiles were reported in response to MeJ [29]. Likewise, the mechanisms by which MeJA regulates paclitaxel biosynthesis were also investigated in Taxus × media cells, and this study showed the potential for mRNAs being targeted by miRNAs [40].
The aim of this study was to identify and characterize novel and conserved miRNAs expressed in MeJ-treated Scots pine needles. In addition, potential precursor sequences and target genes of the identified miRNAs were determined, to understand the type of processes regulated by both conserved and novel miRNAs under stress conditions. A combined strategy—high throughput sequencing and computational prediction—was utilized to identify conserved P. sylvestris miRNAs from six sRNA libraries. The obtained mature miRNA sequences were analyzed and filtered based on known characteristics of plant miRNAs, and compared to other plant miRNAs available in databases. Conserved miRNAs are identified and their role is discussed.

2. Materials and Methods

Six one-year-old Scots pine ramets of one clone derived from Latvian Scots pine breeding program (Sm9-III-2) were transferred into growth chambers two weeks prior to the start of the experiment. Three of the ramets were each treated with 5 ml of a 10 mM MeJA/0.1% Tween 80 solution in deionized water, applied with a hand sprayer. Three control ramets were treated with 5 ml of a 0.1% Tween 80 solution in deionized water. The MeJA treated and control ramets were kept in separate growth chambers. After treatment, ramets were covered with plastic bags for 48 hours to allow the volatilization of excess MeJA. Ramets were grown at 17–22 °C under long day conditions (16 h light + 8 h dark). Two weeks after MeJA treatment, needles were collected and immediately frozen in liquid nitrogen. Control samples were harvested from untreated ramets in the separate growth chamber at the same time. Needle samples were stored at −80 °C until RNA extraction. Total RNA was isolated from 100 mg of needles using a phenol: chloroform: isoamyl alcohol extraction protocol [41]. Total RNA and small RNA quality, quantity and integrity number (for total RNA) was verified using the Agilent Technologies 2100 Bioanalyzer with RNA Agilent RNA 6000 Nano Kit and Agilent Small RNA kit. Total RNA preparations were stored at −80 °C. Total RNA samples were enriched for small RNA as outlined in the Ion RNA-Seq Library Preparation guide (Thermofisher Scientific Manual 4476286 revision E) and 6 small RNA non-barcoded libraries were prepared using Ion Total RNA-Seq Kit v2 for Small RNA Libraries (Part No. 4476289, Thermofisher Scientific) according to the manufacturer’s protocol. Each amplified sRNA library was quantified and the quality analyzed using the Agilent Technologies 2100 Bioanalyzer with a High Sensitivity DNA Kit as intended for templating and separate sequencing as individual non-barcoded libraries on separate sequencing chips. Template-positive Ion Sphere™ Particles (ISPs) were prepared with the Ion OneTouch™ 2 Instrument and enriched with the Ion OneTouch™ ES following the manufacturer’s protocol. The ISP enrichment was then assessed using the Qubit® 2.0 Fluorometer and Ion Sphere™ Quality Control Kitl. Enriched ISPs were then loaded onto an Ion 316 chip (Cat. No. 4483188) and sequenced on an Ion Personal Genome Machine® (PGM™) System at the Norwegian Institute of Bioeconomy Research (NIBIO). The sequences were base called on the Ion Torrent Server with version 4.0.2.
The sRNA sequences were analyzed using the CLC Genomics Workbench software version 7.5.1 (QIAGEN). Low quality reads and adapter sequences were removed and sequences were filtered by length for miRNA identification: minimum length 19 nt and maximum length 25 nt [42,43]. As is known, plant miRNAs tend to be 21 or 22 nt in length and, as previously reported by Axtell and Meyers [11], no RNAs < 20 nt or > 24 nt should be annotated as miRNAs, and annotations of 23-or 24-nt miRNAs require extremely strong evidence; then, cut off limits were used from 19–25 nt. Conserved miRNAs were identified by comparison with miRNA sequences from various tree species (Pinus taeda, Pinus densata, Picea abies, Populus trichocarpa, Populus euphratica, Acacia auriculiformis) as well as other plant species (Arabidopsis thaliana, Oryza sativa, Nicotiana tabacum, Vitis vinifera, Zea mays). Mature miRNA and pre-miRNA sequences of these species were obtained from miRBase (v20 and 21) [44,45]. Two mismatches were allowed between Scots pine miRNA sequences and miRNAs obtained from miRBase. Using the CLC genomics Workbench software, sequences were counted and assigned to families by comparison with mature miRNA sequences from miRBase. The parameters for sequence comparison were: additional downstream bases-2, additional upstream bases-2, missing downstream bases-2, missing upstream bases-2. miRNA sequences obtained from published reports, which were not present in miRBase, were compared to sequences in miRBase as well as sequences obtained in this study. The family classification utilised in miRBase was used to categorize miRNAs not published in miRBase. miRNAs and their isoforms, that contained more than 2 mismatches with published miRNAs were considered as potential new miRNAs for further validation.
In an attempt to verify both the conserved and novel mature miRNAs sequences obtained in this study, potential precursor and target gene sequences were identified. To identify potential Scots pine precursor miRNA sequences, both the conserved as well as unannotated miRNA sequences were aligned to the Pinus PGI_v9.0_032811 unigene sequences (available from DFCI Pine Gene Index at ftp://occams.dfci.harvard.edu/pub/bio/tgi/data/Pinus/). The miRA tool [46] was used for identification of pre-miRNAs. Small RNA sequences were mapped to Pinus unigene sequences (allowing two nucleotide mismatches) using the CLC genomics workbench software, and the miRA tool was used for identification of canonical miRNAs, 5’ and 3’ isomiRs and polymorphic isomiRs and precursor sequences. Previously described criteria [46] were used with some modifications: minimum length (in nt) of the double stranded segment within the folded sequence-19; minimum length (inclusive, in nt) of mature/star miRNA-19; maximum length (exclusive, in nt) of mature/star miRNA-25.
In addition, for conserved miRNA sequences, the minimum free-folding energy index (MFEI) [47] was calculated to confirm that the precursor sequences conformed to the requirements for forming the miRNA precursor structures [11]. Sequences with a maximum of two mismatches with the miRNA sequences were identified and taking into account reports that plant pre-miRNAs vary from approximately 80–200 nt in length [47], regions flanking the mapped mature miRNAs (approximately 150 nt downstream and 150 nt upstream) were used to predict folding structures using the Mfold program (http://mfold.rna.albany.edu/?q=mfold) web server [48] and the CLC genomics workbench software. If the length of a sequence was less than 300 nt, the entire available sequence was used as a miRNA precursor sequence. MFE (minimal negative folding free energy, ΔG), AMFE (adjusted MFE), MFEI (minimal folding free energy index), length of sequence, nucleotide percentage (A, U, G, and C), A + U content, G + C content, and number of base pairs were calculated [49]. A sRNA sequence was considered as a potential miRNA candidate only if it met the following criteria [10,11,47,50]: (1) the predicted mature miRNA was allowed to have 0–2 nucleotide mismatches with the best matched known plant mature miRNA and sequence length was between 19 and 25 nucleotides; (2) a RNA sequence could fold into an appropriate stem–loop hairpin secondary structure; (3) the predicted mature miRNA is located on one arm of the hairpin structure; (4) there were less than 6 mismatches in the complementary site (the opposite miRNA* sequence on the other arm); (5) three mismatched positions of nucleotides in asymmetric bulges; (6) the predicted pre-miRNA had a high negative minimal free– folding energy (MFE) from which the negative folding free energies and MFE index (minimal free folding energy index, MFEI) were calculated in order to distinguish potential precursor miRNAs from other small RNAs. The MFEI was calculated using the formula:
[(MFE/length of the RNA sequence) *100]/(G+C) %
Predicted secondary structures of precursor miRNAs have folding free energy indexes (MFEIs) >= 0.85, distinguishing them from other small RNAs such as tRNAs, and rRNAs whose MFEI are between 0.59 and 0.66; 7) 30–70% A + U content [47]. Application of these criteria can significantly reduce false positive identification of potential precursor miRNAs [50].
Identification of potential miRNA target genes was done by searching for complementary regions between the identified miRNAs in this study and by using all the Pinus transcript sequence input using online web server, and the psRNATarget-Plant Small RNA Target Analysis Server as described previously [51]. Potential target genes were annotated according to GO categories using the Blast2GO PRO Plugin and all non-redundant GeneBank CDS translations + PDB + SwissProt sequences as well SwissProt –non-redundant UniProtKB/SwissProt sequences [52].

3. Results

3.1. Sequencing of Scots Pine Small RNA Libraries

Sequencing of the six small RNA libraries yielded approximately 5.8 million reads before trimming. Prior to trimming, the average length of small RNA reads in the control sample libraries was 21.57 nt and 20.76 nt in the MeJA treated sample libraries. After trimming (19–25 nt), 4.5 million reads remained, with an average length 21.50 nt in the control libraries and 21.46 nt MeJA treated libraries. Small RNAs of 21–22 nt length were the most prevalent among the obtained sequences. In total, 1,021,696 unique small RNA sequences were found.

3.2. Identification of Conserved and Novel miRNAs in Scots Pine

To identify conserved miRNAs expressed in Scots pine, all unique small RNA sequences were compared to annotated mature plant miRNAs in miRBase. Sequences from 11 species in miRBase were utilised, of which six are woody species, including three conifer species-Pinus taeda, Pinus densata, Picea abies. Of the 1,021,696 unique small RNA sequences obtained, 4975 potentially conserved miRNA sequences were identified (consisting of 317,195 reads from a total of 4,488,459). Of these, 957 were ambiguously annotated, i.e., a small RNA sequence was similar to the mature regions of two different miRBase sequences (from the same miRNA family). Of the 4975 annotated sequences (Table 1) identified in our data set, 33.7% were identified from Picea abies, 29.3% from Pinus taeda and 11.8% from Pinus densata (Table 1). Only 0.8% were identified from Acacia auriculiformis and 1.4% from Zea mays. Comparing our data with miRBase sequences, we found that 91.7% of annotated Pinus taeda, 75% of annotated Picea abies and 70% of annotated Pinus densata miRNA sequences were also present in the Scots pine sequences. The least represented miRNA database sequences in Scots pine sequences were from Oryza sativa—only 13.9% of all Oryza sativa sequences present in miRBase were found in Scots pine sequences (Table 2).
The 4975 annotated Scots pine miRNAs were assigned to 173 miRNA groups (Supplementary File 1), based on mature miRNA sequences from miRBase, creating consensus sequences for the identified miRNAs. Two mismatches were allowed between identified Scots pine miRNAs and annotated plant miRNAs in miRBase. The additional parameters were: additional upstream or downstream bases-2, missing upstream or downstream bases-2. The conifer miRNA sequences reported in various publications, but which were not in miRBase, were compared to miRNA sequences in miRBase in order to assign them to miRNA families, provided that the sequences were publically available. These 173 consensus sequences belonged to a total of 60 conserved miRNA families (Supplementary File 1) The majority of miRNA families (32 families), contained 1 consensus sequence (1 group), while the miR159 family (including miR319) contained the most groups (18). Most of the miRNA families found only in conifers contained only 1 or 2 groups.
Only four conifer species are represented in miRBase v21-P.taeda (74 miRNAs), P.densata (60 miRNAs), P.abies (81 miRNAs) and C.lanceolata (9 miRNAs). Comparison of miRBase annotated conifer miRNAs and our data identified 34 conifer specific conserved miRNA families (Table 3). Of these, 34 potentially conifer specific miRNA families, 18 were found in Scots pine, 13 families in loblolly pine and 10 families in Sikang pine, while 24 families were found in Norway spruce. No conifer-specific miRNA families were reported in C. lanceolata. Two families—miR950 and miR1311—were identified in four conifer species (P.sylvestris, P.taeda, P.densata and P.abies). Fourteen conserved conifer miRNA families were not identified in pine species, but were found only in Picea abies. miRNA families miR3699, miR3702 and miR3710 were identified only in P.sylvestris and in Picea abies, but were not found in the other conifer species represented in miRBase.
It was found that all novel and conserved miRNA families were found in both sample types- control and with MJA treated, but the different isomiRs were also found in different sample types.

3.3. Identification of Potential miRNA Precursors

The determination of potential novel miRNAs by identification of miRNA precursors, as well as the identification of precursors for conserved miRNAs, was performed using the miRA tool by mapping the mature miRNA sequences to the Pinus PGI_v9.0_032811 unigene sequences. Additional P. sylvestris sequence databases were also analyzed (e.g., expressed P. sylvestris sequences in GenBank and the draft P. sylvestris genome); however, these yielded a smaller number of potential precursor sequences in comparison to the P. sylvestris unigene set. Therefore, only this database was subsequently utilised.
A total of 1029 potential novel miRNAs that had no homology (as described previously) to miRBase v22 annotations were found. They were grouped into 34 families and 46 predicted precursor sequences were identified (Supplementary File 2). The largest family was miR00005 and contains 124 miRNA isoforms; the smallest family was miR00004, with 3 isoforms (Supplementary File 2). Most of novel miRNAs were located on the 3’arm (526 sequences), 216 novel miRNAs were located on the 5’arm, and 287 miRNA sequences were identified as star sequences.
Analyzing the 4975 Scots pine mature miRNA sequences with homology to miRBase miRNAs, 50 potential precursor sequences for 2209 of these mature miRNAs were identified. These 50 sequences were predicted to be precursors of 20 families (Supplementary File 3). The most isomiRs (483), were found for miR950, but the most precursors were found for miR482 and miR950.
Using the miRA tool, only 19 potential precursor sequences for 780 of these mature miRNAs were identified (Supplementary File 4). These 19 sequences were predicted to be precursors of 9 miRNA families-miR396, miR482, miR946, miR949, miR952, miR1312, miR1313, miR1314 and miR3701. All 19 families were also found using manual searching and selection.
A large number of isomiRs that were homologous to each of these families were identified with the miRA tool, ranging from 242 isomiRs for the miR482 family to 23 isomiRs for the miR946 family. Only one potential precursor sequence was found for the miRNA families (miR946, miR952, miR1312, miR1313), while five potential precursor sequences were found for the family miR482. Of the 780 isomiRs, 342 were found on the 3’ arm, 250 on the 5’ arm, but 188 isomiRs were identified as star sequences. Only one precursor identified with the miRA tool (TC159053 for miR952) was not identified by manual searching and selection.
In addition to the miRA tool, potential precursor sequences for conserved miRNAs were analyzed for the ability to form the required stem-loop structure, and the minimal folding free energy indexes were calculated (Supplementary File 3). Predicted pre-miRNA sequences were trimmed in the primary miRNA sequence region until the next bulge or loop after the miRNA* region. Minimum folding free energy indexes ranged from 0.72–1.31, with most being > 0.85 (47 precursor sequences (94%)) (Supplementary File 3), and corroborated with previously reported sequences. Three predicted precursors (for more than 250 isomiRs) had an MFEI less than 0.85. A number of sequences were identified as potential precursors for multiple miRNAs which are in the same family. This is due to the parameters utilized, which allows a maximum of two mismatches between a miRNA and a potential precursor miRNA sequence. Two or more than two potential precursor sequences were identified for miR396, miR482, miR949, miR1314, miR3701. In some cases, the bioinformatic identification of potential precursor sequences was complicated, because the mature miRNA was identical to the sequence in one precursor, but the MFEI was lower than in a different potential precursor with less homology to the mature miRNA sequence. These mismatches between potential precursors and mature miRNAs, as well as between two potential precursors could be due to the sequences being derived from different individuals. Another possibility is that either the true precursor or mature miRNA sequences are not present in the data set or databases.

3.4. miRNA Target Identification

Analysis of conserved miRNA putative target genes identified 119 genes targeted by 58 miRNA families (Supplementary File 3). Seven of these target genes, TC195763, TC162935, DR694512, CF470498, TC170132, TC159334, and TC195763, were targets of two different miRNA families. Ten of the putative gene targets (targeted by eight miRNA families) were of unknown function. Of these eight families, four miRNA families targeted one or more unannotated genes, but four families targeted multiple genes with both unknown and known functions. No putative targets were identified for seven miRNA families (defined as “no result”). However, putative target genes for six of these miRNA families have been described in other publications [53,54,55,56,57,58].
Six of the identified putative target genes (targeted by five miRNA families—miR393, miR950, miR951, miR3623, and miR024) were homologous to TIR/TIR like/NBS-LRR disease-resistance protein genes and an additional three targets (targeted by three miRNA families) were with unknown function, but these miRNAs have been previously described as targeting disease resistance protein genes. Plant NBS-LRR proteins are involved in the detection of pathogen-associated proteins, most often the effector molecules of pathogens which are responsible for virulence.
GO annotation [52] of the identified target genes indicated that the most common functions in the biological process domain (Supplementary File 5, Figure S1) were related to transcription regulation and signal transduction, as well as protein phosphorylation and ubiquitination, and response to stress. The most common GO annotations in the molecular function domain (Supplementary File 5, Figure S2) were binding, including DNA binding, and transcription factor activity.
Analysis of novel miRNA putative target genes identified 136 genes targeted by 28 families (from a total of 34 potential novel miRNA families identified by precursor sequence analysis) (Supplementary File 2) or 217 isoforms (from 1029 miRNA isoforms) with 471 target sites. These results indicated that the same isoforms can target not only different target genes, but also more than one target site within a single target gene, and also that different isoforms from more than one family can target the same target gene (sequences with red color in Supplementary File 3). The largest number of target genes were predicted for family miR00002 with 20 targets. The most target sites (65) were found for family miR00025, where 38 isoforms targeted three target genes and family miR00027, where 15 isoforms targeted three target genes and 58 target sites. Fourteen miRNA families were homologous to resistance-related target genes. Four potential target genes (targeted by five miRNA families or 19 isoforms) were of unknown function. Sixteen target genes were homologous to alcohol dehydrogenases, eight target genes were homologous to phytocyanins and seven were homologous to peroxidases. Using the Swissprot protein database, the most common GO annotations of the identified target genes in the biological process domain (Supplementary File 5, Figures S3 and S4) were response to stimulus (including response to abscisic acid, response to cytokinin, heat acclimation, response to wounding, hypoxia and others), single organisms process and cellular process. The most common GO annotations in the molecular function domain were binding and catalytic activity (Supplementary File 5, Figure S5). Utilizing the larger to non-redundant protein database, the most common GO annotations of the identified target genes in the biological process domain (Supplementary File 5, Figure S6) were metabolic processes, single organism processes, cellular processes and responses to stimuli. However, the most common GO annotations in the molecular function domain were the same-binding and catalytic activity (Supplementary File 5, Figure S7).

4. Discussion

Potential novel miRNAs were determined by identifying precursor sequences and target genes. Most novel isomiRs were located on the 3′ arm of the precursor stem-loop structure, similarly to the conserved miRNA sequences, and corresponding to previous reports [59]. miRNA isoforms can be processed not only by the Dicer enzyme, but also by other Dicer-like enzymes and Dicer independent mechanisms. This suggests that only a subset of small RNAs can be identified as mature miRNAs in concordance with a Dicer mediated miRNA biogensis pathway, but other isoforms might be produced by other Dicer-like enzymes or by Dicer-independent mechanisms, or can be caused by recurrent somatic mutations in Drosha, which can induce changes in miRNA expression, and somatic mutations in Drosha and Dicer1 can impair miRNA biogenesis [60,61].
It has been reported that more than 90% of miRNA precursors had an MFEI greater than 0.85, and no mRNAs, tRNAs, or rRNAs had more MFEI higher than 0.85 [62]. Our data indicate that most minimum folding free energy indexes were > 0.85 and included 94% of all precursors, in corroboration with previously reported data.
Comparing our data, we concluded that only 23 from 58 miRNAs—including isomiRs, from Taxus chinesis [29], which also had treatment with MJA—were also found in our data, from which more than half (15) of the miRNA sequences were found in the same sample (control or MJA-treated) type as in T. chinesis. In some other gymnosperm miRNA studies [25,28,63] more miRNA families were found that have precursors, but in this study, more isomiRs were found for these families, compared to the previous three studies. These differences in identified isomiR number between studies is probably due to the utilized methodology and parameters.
A survey of miRNA and other small RNA sequences across a wide range of species has identified conserved miRNA families [64], and the miRNA families identified from our study is broadly consistent with this report. The five miRNA families classified as ubiquitous (miR156, miR166, miR167, miR168, miR172) were all identified in our data. In addition, the 16 miRNA families classified as present in most taxonomic groups were also present in our data, with the exception of miR4414. None of the miRNA families enriched in dicots were found in our data, and they were also poorly represented in other gymnosperm species (Cycas rumphi, Ginkgo biloba and Picea abies) [64]. All three miRNA families reported to be enriched in gymnosperms (miR536, miR1083, miR1314) were identified in our dataset. Of the 20 miRNA families enriched in monocots, only one (miR1432) was identified in this study. All of the miRNA families conserved in Arabidopsis, Oryza and Populus [65], were also identified within the sequences expressed in Pinus sylvestris. However, the miRNAs that were enriched in monocots or dicots were not identified in P.sylvestris [64]. A number of putative gymnosperm-specific miRNAs were identified; however, in many cases, their sequences were less conserved than the highly conserved plant lineage miRNAs identified. For example, the most highly conserved miRNA sequences were found between P. sylvestris and A. thaliana, while the miRNA families that were common with other gymnosperm species had one or more nucleotide mismatches compared to the P. sylvestris sequences.
Many of the highly conserved miRNAs are found in a range of plant groups (e.g., monocots, dicots and gymnosperms), suggesting a common ancestry and function for these miRNAs. However, some of the P. sylvestris miRNA families were similar to those miRNAs only reported to be present in monocots (i.e., miR1432, miR2275, miR5072, miR5083). In some cases, this could be due to the underrepresentation of a particular miRNA family in miRBase (e.g., miR5072 or miR5083—which are only reported in rice). However, in other cases (miR1432, miR2275) these families are reported to be present in a range of monocot species, including rice, sorghum, maize, Brachypodium and Aegilops). These miRNA families were also reported at least one other gymnosperm species, not only in monocots. Two conserved miRNA families identified in P. sylvestris (miR6024, miR6478) have not been previously reported in any other gymnosperm species.
The miRNA families miR156, miR159/319, miR160, miR166, miR171, and miR408 are reported to be present in all green plants (embryophyta); the miR396 family is present in all vascular plants (tracheophytes); while miR397 and miR398 are present in all seed plants (spermatophytes) [66]. The miR403, miR828 and miR2111 families were reported as eudicot-specific families; however, these miRNA families, except for miR403 (restricted to core eudicot lineages), were also identified in P.sylvestris in this study. Only some families-miR156, miR159/319, miR160, miR166, miR171, miR408, miR390/391, miR395, miR529 and miR536-identified in bryophytes were also found in P.sylvestris, which could be a result of the relatively small number of miRNAs reported in bryophytes.
The higher proportion of perfect matches with angiosperms rather than conifers may be due to the low number of conifer sequences compared to angiosperms and the lack of conserved conifer miRNA sequences in the database, or it may indicate that miRNA sequences are more diverse in conifers than in angiosperms. Conserved miRNA families are relatively easy to identify, not only because they are annotated with a higher confidence from a range of species, but also because they are reported to be more abundant, with the 21 most highly conserved miRNA families accounting for 54–98% of miRNA sequences [64]. Therefore, species or lineage-specific miRNA families may be unrecognized due to inadequate sequence coverage in addition to the difficulties in unambiguously identifying novel miRNAs.

4.1. Comparison of Scots Pine miRNA Families with Published Conifer miRNAs

A literature review identified 22 publications reporting on miRNA studies in conifers, of which seven publications were about studies in five different pine species—P.taeda, P.densata, P.strobus, P.contorta, P.tabuliformis. Of these, miRNA sequences from only four publications are available in miRbase: Lu et al. [33] (Pinus taeda), Wan et al. [67] (Pinus densata), Yakovlev et al. [26] (Picea abies) and Wan et al. [28] (Cunninghamia lanceolata). Comparing our data with previously reported miRNA families, we found that of the 58 conserved miRNA families identified in Scots pine, 56 also were reported in at least one other conifer species (Supplementary File 2). Eighteen miRNA families were found only in conifer species (Pinus taeda, Pinus densata, Picea abies): miR950, miR951, miR952, miR946, miR947, miR949, miR3699, miR3701, miR3702, miR3710, miR3712, miR1309, miR1311, miR1312, miR1313, miR1314, miR1315, and miR1316. The miR6024 and miR6478 families were identified only in Scots pine, but two families, miR5072 and miR5083, were found in Scots pine and in only one additional conifer species. The functional annotation of the target genes of these conserved miRNA families indicates that most are involved in plant growth and development, biotic and abiotic stress response and disease resistance.

4.2. Target Genes

The majority of highly conserved miRNA families (e.g., as reported in Cuperus et al. [66]), which are found in all plants including bryophytes, target transcription factors such as SPL, MYB, HD ZIP III, AP2 and others [68]. In contrast, none of the gymnosperm-specific miRNA target genes were annotated as having DNA-binding functions, suggesting that transcription factors are underrepresented in this group. However, this could also be an artefact due to the small number of gymnosperm-specific miRNAs and target genes. In addition, transcription factors may be underrepresented in the publicly available P. sylvestris gene sequences. In this study, the predicted target genes included stress-related gene families such as Leucine Rich Repeat (LRR) protein genes, protein kinase domain, Toll-Interleukin receptor (TIR) domain, disease resistance protein RPP1-WsB, heat shock proteins, and blue cooper proteins (including uclacyanin, plantacyanin). Transcription factors were also among the identified target genes-CCAAT-binding transcription factors, AP2-related transcription factors, homocysteine S-methyltransferases, MYB-like proteins, regulatory protein GntR, MarR family genes, GAMYB, and zinc finger protein genes. Target genes also included transferases and other enzymes. A similar functional distribution of target genes was also reported in crop species, as almost 50% of the miRNA targets were transcription factors in pathways that are likely important in setting or maintaining the developmental program leading to high quality soybean seeds [69], and 66% of target genes in crops are reported to be transcription factors; however, 11% are lncRNA, 7% are NB-LRR proteins, 2% are pathogen proteins and 24% are other proteins [70].

4.3. Resistance/Stress Genes

TIR-NBS-LRR and other NBS-LRR genes are involved in defense responses to pathogen infection and disease resistance. This is a highly diverse class of genes that can include highly conserved genes as well as lineage- or species-specific genes. They can form gene clusters, and are regulated by miRNAs [71]. The TIR-NBS-LRR gene family has been lost in monocots [72], but is present in gymnosperms. The conserved miR482/2118 family targets the highly-conserved P-loop motif in NBS genes [71], and was also identified in this study. Three main mechanisms of miRNA locus evolution have been proposed: the inverted duplication of a target gene sequence, leading to the formation of the required stem-loop structure [73]; the formation of a stem-loop structure by the self-complementarity of a transcribed sequence [74]; the mutation of existing miRNA sequences [75]. Zhang et al. [71] suggest that the clustering of NBS-LRR genes facilitates the expression of inverted tandem duplications of target genes, thereby facilitating the tandem co-evolution of miRNAs and their target NBS-LRR genes. Of the gymnosperm/conifer specific miRNA families identified in this study, at least two target NBS-LRR genes (miR1311 and miR1312) and the targets of the majority of the others are resistance genes. Conifer genomes are reported to contain a larger proportion of gene families in comparison to angiosperms [76], and contain a large proportion of repeated sequences and retrotransposons [77,78] which can facilitate the formation of inverted genes and other sequences. Further characterization of novel miRNAs in Scot pines will reveal if these unique miRNAs share similar origins.

4.4. Transcription Factors

Many genes are activated in response to stresses at the transcriptional level, and they provide stress tolerance through the production of vital metabolic proteins and also by regulating downstream genes [79]. Transcription factors (TFs) are essential for the regulation of gene expression, and usually belong to members of multigene families [80]. TFs families can evolve by a range of mechanisms such as exon capture, duplication, translocation and mutation [81,82]. The majority of miRNA targets are transcription factors, which regulate plant growth and developments [83,84,85]. The role of miRNAs in the regulation of transcription factors influencing traits such as meristem identity, polarity and flowering has been established; however, even highly conserved miRNA families can be involved in different developmental processes, and their role and mechanisms of action can vary between species, or even between different tissues types [86]. Therefore, the role of these conserved miRNAs must be further investigated in conifers, as these miRNA families may have evolved to play different roles to model angiosperm species.

5. Conclusions

A large number of mature and precursor miRNA sequences have been identified and published in miRBase and other databases, but there are still problems with the correct identification and annotation of miRNAs and the subsequent accuracy of available information [87,88]. The consequences of this are inconsistencies in the naming of miRNA families and isoforms, leading to difficulties in comparisons between studies. It is also difficult to identify true precursor miRNAs based on homology searches with mature miRNA sequences. There are several tools available for the prediction (MiRFinder, Rfam, MIReNA, etc.) of precursor miRNAs, but distinguishing true pre-miRNAs from other hairpin sequences with stem loop structures can be complicated, as these structures are not unique to pre-miRNAs, and many other coding or non-coding RNAs—such as rRNAs, tRNAs and mRNAs—can also form similar hairpin structures.
The majority of highly conserved plant miRNAs were identified in this study, as well as some conserved miRNAs previously reported to be monocot specific. No conserved dicot-specific miRNAs were identified. A number of potential gymnosperm or conifer specific miRNAs were found, shared among a range of conifer species. Potential target genes were identified, of which the targets of highly conserved miRNAs present in most plant families were transcription factors, while the conserved conifer-specific miRNA targets were involved in disease resistance.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4907/11/4/384/s1, Supplementary File 1: Table S1: Conserved miRNAs grouped in groups by consensus sequence, Supplementary File 2: Novel miRNAs, their targets, families and Blast2Go data, Supplementary File 3: Excel file S1: conserved miRNA precursors, targets, Supplementary File 4: Excel file S1: identified conserved miRNA precursors by miRA tool. Supplementary File 5: Figure S1: Conserved miRNA target sequences distribution in category of biological process, Figure S2: Conserved miRNA target sequences distribution in category of molecular function, Figure S3: Novel miRNA target sequences distribution in category of biological process in SwissProt database, Figure S4: Novel miRNA target sequences distribution in category in response to stimulus in SwissProt database, Figure S5: Novel miRNA target sequences distribution in category of molecular function in SwissProt database, Figure S6: Novel miRNA target sequences distribution in category of biological process in non-redundant protein database, Figure S7: Novel miRNA target sequences distribution in category of biological process in non-redundant protein database.

Author Contributions

B.K., D.R., C.G.F. conceived and designed the research, I.Š. cultivated, treated and provided the scots pine ramets, B.K., I.Š. and V.Š. performed the experiments, B.K. and D.R. analyzed and investigated the data and wrote the manuscript (original draft paper), A.V.-S., I.Y. and C.G.F. trained, consulted and assisted with sequencing with the IonTorrent PGM. All authors reviewed and edited the manuscript and approve the published version of the manuscript.

Funding

This project was funded by the Latvian Council of Science grant 284/2012 “Investigation of molecular defense mechanisms in Scots pine (Pinus sylvestris L.).

Conflicts of Interest

The authors declare no completing interests.

Data Availability

Unique small RNA sequences were deposited to the SRA (Short Read Archive, NCBI) with the accession number PRJNA531446 (https://www.ncbi.nlm.nih.gov/sra/?term=PRJNA531446).

References

  1. Bräutigam, K.; Vining, K.J.; Lafon-Placette, C.; Fossdal, C.G.; Mirouze, M.; Marcos, J.G.; Fluch, S.; Fraga, M.F.; Guevara, M.Á.; Abarca, D.; et al. Epigenetic regulation of adaptive responses of forest tree species to the environment. Ecol. Evol. 2013, 3, 399–415. [Google Scholar] [CrossRef] [PubMed]
  2. Rohde, A.; Junttila, O. Remembrances of an embryo: Long-term effects on phenology traits in spruce. New Phytol. 2008, 177, 2–5. [Google Scholar] [CrossRef] [PubMed]
  3. Carrington, J.C.; Ambros, V. Role of microRNAs in plant and animal development. Science 2003, 301, 336–338. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Bartel, D.P. MicroRNAs: Genomics, biogenesis, mechanism, and function. Cell 2004, 116, 281–297. [Google Scholar] [CrossRef] [Green Version]
  5. Borges, F.; Martienssen, R.A. The expanding world of small RNAs in plants. Nat. Rev. Mol. Cell Biol. 2015, 16, 727–741. [Google Scholar] [CrossRef] [Green Version]
  6. Wang, H.; Jiao, X.; Kong, X.; Hamera, S.; Wu, Y.; Chen, X.; Fang, R.; Yan, Y. A signaling cascade from miR444 to RDR1 in rice antiviral RNA silencing pathway. Plant Physiol. 2016, 170, 2365–2377. [Google Scholar] [CrossRef] [Green Version]
  7. Cai, Q.; Qiao, L.; Wang, M.; He, B.; Lin, F.M.; Palmquist, J.; Da Huang, S.; Jin, H. Plants send small RNAs in extracellular vesicles to fungal pathogen to silence virulence genes. Science 2018, 360, 1126–1129. [Google Scholar] [CrossRef] [Green Version]
  8. Neutelings, G.; Fénart, S.; Lucau-Danila, A.; Hawkins, S. Identification and characterization of miRNAs and their potential targets in flax. J. Plant Physiol. 2012, 169, 1754–1766. [Google Scholar] [CrossRef]
  9. Desvignes, T.; Batzel, P.; Berezikov, E.; Eilbeck, K.; Eppig, J.T.; McAndrews, M.S.; Singer, A.; Postlethwait, J.H. MiRNA Nomenclature: A view incorporating genetic origins, biosynthetic pathways, and sequence variants. Trends Genet. 2015, 31, 613–626. [Google Scholar] [CrossRef] [Green Version]
  10. Meyers, B.C.; Axtell, M.J.; Bartel, B.; Bartel, D.P.; Baulcombe, D.; Bowman, J.L.; Cao, X.; Carrington, J.C.; Chen, X.; Green, P.J.; et al. Criteria for annotation of plant microRNAs. Plant Cell 2008, 20, 3186–3190. [Google Scholar] [CrossRef]
  11. Axtell, M.J.; Meyers, B.C. Revisiting criteria for plant microRNA annotation in the Era of big data. Plant Cell 2018, 30, 272–284. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Zhai, J.; Zhao, Y.; Simon, S.A.; Huang, S.; Petsch, K.; Arikit, S.; Pillay, M.; Ji, L.; Xie, M.; Cao, X.; et al. Plant microRNAs display differential 39 truncation and tailing modifications that are ARGONAUTE1 dependent and conserved across species. Plant Cell 2013, 25, 2417–2428. [Google Scholar] [CrossRef] [PubMed]
  13. Newman, M.A.; Mani, V.; Hammond, S.M. Deep sequencing of microRNA precursors reveals extensive 39 end modification. RNA 2011, 17, 1795–1803. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Wyman, S.K.; Knouf, E.C.; Parkin, R.K.; Fritz, B.R.; Lin, D.W.; Dennis, L.M.; Krouse, M.A.; Webster, P.J.; Tewari, M. Post-transcriptional generation of miRNA variants by multiple nucleotidyl transferases contributes to miRNA transcriptome complexity. Genome Res. 2011, 21, 1450–1461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Tan, G.C.; Chan, E.; Molnar, A.; Sarkar, R.; Alexieva, D.; Isa, I.M.; Robinson, S.; Zhang, S.; Ellis, P.; Langford, C.F.; et al. 5′ isomiR variation is of functional and evolutionary importance. Nucleic Acids Res. 2014, 42, 9424–9435. [Google Scholar] [CrossRef]
  16. Morin, R.D.; Aksay, G.; Dolgosheina, E.; Ebhardt, H.A.; Magrini, V.; Mardis, E.R.; Sahinalp, S.C.; Unrau, P.J. Comparative analysis of the small RNA transcriptomes of Pinus contorta and Oryza sativa. Genome Res. 2008, 18, 571–584. [Google Scholar] [CrossRef] [Green Version]
  17. Fei, Y.; Xiao, B.; Yang, M.; Ding, Q.; Tang, W. MicroRNAs, polyamines, and the activities antioxidant enzymes are associated with in vitro rooting in white pine (Pinus strobus L.). Springerplus 2016, 5, 1–11. [Google Scholar] [CrossRef] [Green Version]
  18. Zhang, S.; Zhou, J.; Han, S.; Yang, W.; Li, W.; Wei, H.; Li, X.; Qi, L. Four abiotic stress-induced miRNA families differentially regulated in the embryogenic and non-embryogenic callus tissues of Larix leptolepis. Biochem. Biophys. Res. Commun. 2010, 398, 355–360. [Google Scholar] [CrossRef]
  19. Zhang, J.; Zhang, S.; Han, S.; Wu, T.; Li, X.; Li, W.; Qi, L. Genome-wide identification of microRNAs in larch and stage-specific modulation of 11 conserved microRNAs and their targets during somatic embryogenesis. Planta 2012, 236, 647–657. [Google Scholar] [CrossRef]
  20. Zhang, J.; Zhang, S.; Han, S.; Li, X.; Tong, Z.; Qi, L. Deciphering small noncoding RNAs during the transition from dormant embryo to germinated embryo in larches (Larix leptolepis). PLoS ONE 2013, 8, e81452. [Google Scholar] [CrossRef]
  21. Zhang, J.; Wu, T.; Li, L.; Han, S.; Li, X.; Zhang, S.; Qi, L. Dynamic expression of small RNA populations in larch (Larix leptolepis). Planta 2013, 237, 89–101. [Google Scholar] [CrossRef]
  22. Zhang, J.; Zhang, S.; Li, S.; Han, S.; Wu, T.; Li, X.; Qi, L. A genome-wide survey of microRNA truncation and 3′ nucleotide addition events in larch (Larix leptolepis). Planta 2013, 237, 1047–1056. [Google Scholar] [CrossRef]
  23. Källman, T.; Chen, J.; Gyllenstrand, N.; Lagercrantz, U. A significant fraction of 21-nucleotide small RNA originates from phased degradation of resistance genes in several perennial species 1[C][W][OA]. Plant Physiol. 2013, 162, 741–754. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Fossdal, C.G.; Yaqoob, N.; Krokene, P.; Kvaalen, H.; Solheim, H.; Yakovlev, I.A. Local and systemic changes in expression of resistance genes, NB-LRR genes and their putative microRNAs in norway spruce after wounding and inoculation with the pathogen ceratocystis polonica. BMC Plant Biol. 2012, 12, 1–11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Xia, R.; Xu, J.; Arikit, S.; Meyers, B.C. Extensive families of miRNAs and PHAS loci in Norway spruce demonstrate the origins of complex phasiRNA networks in seed plants. Mol. Biol. Evol. 2015, 32, 2905–2918. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Yakovlev, I.A.; Fossdal, C.G.; Johnsen, Ø. MicroRNAs, the epigenetic memory and climatic adaptation in Norway spruce. New Phytol. 2010, 187, 1154–1169. [Google Scholar] [CrossRef]
  27. Qiu, Z.; Li, X.; Zhao, Y.; Zhang, M.; Wan, Y.; Cao, D.; Lu, S.; Lin, J. Genome-wide analysis reveals dynamic changes in expression of microRNAs during vascular cambium development in Chinese fir, Cunninghamia lanceolata. J. Exp. Bot. 2015, 66, 3041–3054. [Google Scholar] [CrossRef] [Green Version]
  28. Wan, L.C.; Wang, F.; Guo, X.; Lu, S.; Qiu, Z.; Zhao, Y.; Zhang, H.; Lin, J. Identification and characterization of small non-coding RNAs from Chinese fir by high throughput sequencing. BMC Plant Biol. 2012, 12. [Google Scholar] [CrossRef] [Green Version]
  29. Qiu, D.; Pan, X.; Wilson, I.W.; Li, F.; Liu, M.; Teng, W.; Zhang, B. High throughput sequencing technology reveals that the taxoid elicitor methyl jasmonate regulates microRNA expression in Chinese yew (Taxus chinensis). Gene 2009, 436, 37–44. [Google Scholar] [CrossRef]
  30. Chen, Y.T.; Shen, C.H.; Lin, W.D.; Chu, H.A.; Huang, B.L.; Kuo, C.I.; Yeh, K.W.; Huang, L.C.; Chang, I.F. Small RNAs of Sequoia sempervirens during rejuvenation and phase change. Plant Biol. 2013, 15, 27–36. [Google Scholar] [CrossRef]
  31. Quinn, C.R.; Iriyama, R.; Fernando, D.D. Expression patterns of conserved microRNAs in the male gametophyte of loblolly pine (Pinus taeda). Plant Reprod. 2014, 27, 69–78. [Google Scholar] [CrossRef]
  32. Oh, T.J.; Wartell, R.M.; Cairney, J.; Pullman, G.S. Evidence for stage-specific modulation of specific microRNAs (miRNAs) and miRNA processing components in zygotic embryo and female gametophyte of loblolly pine (Pinus taeda). New Phytol. 2008, 179, 67–80. [Google Scholar] [CrossRef] [PubMed]
  33. Lu, S.; Sun, Y.H.; Amerson, H.; Chiang, V.L. MicroRNAs in loblolly pine (Pinus taeda L.) and their association with fusiform rust gall development. Plant J. 2007, 51, 1077–1098. [Google Scholar] [CrossRef]
  34. Staswick, P.E.; Su, W.; Howell, S.H. Methyl jasmonate inhibition of root growth and induction of a leaf protein are decreased in an Arabidopsis thaliana mutant. Proc. Natl. Acad. Sci. USA 1992, 89, 6837–6840. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Wasternack, C. Jasmonates: An update on biosynthesis, signal transduction and action in plant stress response, growth and development. Ann. Bot. 2007, 100, 681–697. [Google Scholar] [CrossRef] [Green Version]
  36. Kang, J.H.; Wang, L.; Giri, A.; Baldwin, I.T. Silencing threonine deaminase and JAR4 in Nicotiana attenuata impairs jasmonic acid-isoleucine-mediated defenses against Manduca sexta. Plant Cell 2006, 18, 3303–3320. [Google Scholar] [CrossRef] [Green Version]
  37. Martin, D.; Tholl, D.; Gershenzon, J.; Bohlmann, J. Methyl jasmonate induces traumatic resin ducts, terpenoid resin biosynthesis, and terpenoid accumulation in developing xylem of Norway spruce stems. Plant Physiol. 2002, 129, 1003–1018. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Hudgins, J.W.; Franceschi, V.R. Methyl jasmonate-induced ethylene production is responsible for conifer phloem defense responses and reprogramming of stem cambial zone for traumatic resin duct formation. Plant Physiol. 2004, 135, 2134–2149. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Zeneli, G.; Krokene, P.; Christiansen, E.; Krekling, T.; Gershenzon, J. Methyl jasmonate treatment of mature Norway spruce (Picea abies) trees increases the accumulation of terpenoid resin components and protects against infection by Ceratocystis polonica, a bark beetle-associated fungus. Tree Physiol. 2006, 26, 977–988. [Google Scholar] [CrossRef] [PubMed]
  40. Sun, G.; Yang, Y.; Xie, F.; Wen, J.-F.; Wu, J.; Wilson, I.W.; Tang, Q.; Liu, H.; Qiu, D. Deep Sequencing Reveals Transcriptome Re-Programming of Taxus × media Cells to the Elicitation with Methyl Jasmonate. PLoS ONE 2013, 8, e62865. [Google Scholar] [CrossRef] [Green Version]
  41. Rubio-Piña, J.A.; Zapata-Pérez, O. Isolation of total RNA from tissues rich in polyphenols and polysaccharides of mangrove plants. Electron. J. Biotechnol. 2011, 14. [Google Scholar] [CrossRef]
  42. Singh, N.; Srivastava, S.; Sharma, A. Identification and analysis of miRNAs and their targets in ginger using bioinformatics approach. Gene 2016, 575, 570–576. [Google Scholar] [CrossRef] [PubMed]
  43. Jike, W.; Sablok, G.; Bertorelle, G.; Li, M.; Varotto, C. In silico identification and characterization of a diverse subset of conserved microRNAs in bioenergy crop Arundo donax L. Sci. Rep. 2018, 8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Kozomara, A.; Griffiths-Jones, S. MiRBase: Integrating microRNA annotation and deep-sequencing data. Nucleic Acids Res. 2011, 39. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Kozomara, A.; Griffiths-Jones, S. MiRBase: Annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 2014, 42. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Evers, M.; Huttner, M.; Dueck, A.; Meister, G.; Engelmann, J.C. miRA: Adaptable novel miRNA identification in plants using small RNA sequencing data. BMC Bioinform. 2015, 16, 1–10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Zhang, B.; Pan, X.; Wang, Q.; Cobb, G.P.; Anderson, T.A. Computational identification of microRNAs and their targets. Comput. Biol. Chem. 2006, 30, 395–407. [Google Scholar] [CrossRef]
  48. Zuker, M. Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31, 3406–3415. [Google Scholar] [CrossRef]
  49. Zhang, B.; Pan, X.; Cobb, G.P.; Anderson, T.A. Plant microRNA: A small regulatory molecule with big impact. Dev. Biol. 2006, 289, 3–16. [Google Scholar] [CrossRef]
  50. Ambros, V.; Bartel, B.; Bartel, D.P.; Burge, C.B.; Carrington, J.C.; Chen, X.; Dreyfuss, G.; Eddy, S.R.; Griffiths-Jones, S.; Marshall, M.; et al. A uniform system for microRNA annotation. RNA 2003, 9, 277–279. [Google Scholar] [CrossRef] [Green Version]
  51. Dai, X.; Zhao, P.X. PsRNATarget: A plant small RNA target analysis server. Nucleic Acids Res. 2011, 39, 155–159. [Google Scholar] [CrossRef] [Green Version]
  52. Conesa, A.; Götz, S. Blast2GO: A comprehensive suite for functional analysis in plant genomics. Int. J. Plant Genom. 2008, 2008. [Google Scholar] [CrossRef]
  53. Liu, X.; Huang, J.; Wang, Y.; Khanna, K.; Xie, Z.; Owen, H.A.; Zhao, D. The role of floral organs in carpels, an Arabidopsis loss-of-function mutation in MicroRNA160a, in organogenesis and the mechanism regulating its expression. Plant J. 2010, 62, 416–428. [Google Scholar] [CrossRef] [PubMed]
  54. Chen, X. microRNA biogenesis and function in plants. FEBS Lett. 2005, 579, 5923–5931. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Aida, M.; Ishida, T.; Fukaki, H.; Fujisawa, H.; Tasaka, M. Genes involved in organ separation in Arabidopsis: An analysis of the cup-shaped cotyledon mutant. Plant Cell 1997, 9, 841–857. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Gray, W.M.; Kepinski, S.; Rouse, D.; Leyser, O.; Estelle, M. Auxin regulates SCFTIR1-dependent degradation of AUX/IAA proteins. Nature 2001, 414, 271–276. [Google Scholar] [CrossRef] [PubMed]
  57. Xie, Z.; Kasschau, K.D.; Carrington, J.C. Negative feedback regulation of Dicer-Like1 in Arabidopsis by microRNA-guided mRNA degradation. Curr. Biol. 2003, 13, 784–789. [Google Scholar] [CrossRef] [Green Version]
  58. Khraiwesh, B.; Zhu, J.K.; Zhu, J. Role of miRNAs and siRNAs in biotic and abiotic stress responses of plants. Biochim. Biophys. Acta Gene Regul. Mech. 2012, 1819, 137–148. [Google Scholar] [CrossRef] [Green Version]
  59. Guo, L.; Chen, F. A challenge for miRNA: Multiple isomiRs in miRNAomics. Gene 2014, 544, 1–7. [Google Scholar] [CrossRef]
  60. Rakheja, D.; Chen, K.S.; Liu, Y.; Shukla, A.A.; Schmid, V.; Chang, T.C.; Khokhar, S.; Wickiser, J.E.; Karandikar, N.J.; Malter, J.S.; et al. Somatic mutations in DROSHA and DICER1 impair microRNA biogenesis through distinct mechanisms in Wilms tumours. Nat. Commun. 2014, 2. [Google Scholar] [CrossRef]
  61. Torrezan, G.T.; Ferreira, E.N.; Nakahata, A.M.; Barros, B.D.F.; Castro, M.T.M.; Correa, B.R.; Krepischi, A.C.V.; Olivieri, E.H.R.; Cunha, I.W.; Tabori, U.; et al. Recurrent somatic mutation in DROSHA induces microRNA profile changes in Wilms tumour. Nat. Commun. 2014, 5. [Google Scholar] [CrossRef] [Green Version]
  62. Zhang, B.H.; Pan, X.P.; Cox, S.B.; Cobb, G.P.; Anderson, T.A. Evidence that miRNAs are different from other RNAs. Cell. Mol. Life Sci. 2006, 63, 246–254. [Google Scholar] [CrossRef] [PubMed]
  63. Yakovlev, I.A.; Fossdal, C.G. In silico analysis of small RNAs suggest roles for novel and conserved miRNAs in the formation of epigenetic memory in somatic embryos of Norway spruce. Front. Physiol. 2017, 8. [Google Scholar] [CrossRef] [PubMed]
  64. Chávez Montes, R.A.; De Fátima Rosas-Cárdenas, F.; De Paoli, E.; Accerbi, M.; Rymarquis, L.A.; Mahalingam, G.; Marsch-Martínez, N.; Meyers, B.C.; Green, P.J.; De Folter, S. Sample sequencing of vascular plants demonstrates widespread conservation and divergence of microRNAs. Nat. Commun. 2014, 5, 3722. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Jones-Rhoades, M.W.; Bartel, D.P.; Bartel, B. MicroRNAs and their regulatory roles in plants. Annu. Rev. Plant Biol. 2006, 57, 19–53. [Google Scholar] [CrossRef]
  66. Cuperus, J.T.; Fahlgren, N.; Carrington, J.C. Evolution and functional diversification of MIRNA genes. Plant Cell 2011, 23, 431–442. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Wan, L.C.; Zhang, H.; Lu, S.; Zhang, L.; Qiu, Z.; Zhao, Y.; Zeng, Q.Y.; Lin, J. Transcriptome-wide identification and characterization of miRNAs from Pinus densata. BMC Genom. 2012, 13. [Google Scholar] [CrossRef] [Green Version]
  68. Samad, A.F.A.; Sajad, M.; Nazaruddin, N.; Fauzi, I.A.; Murad, A.M.A.; Zainal, Z.; Ismail, I. MicroRNA and transcription factor: Key players in plant regulatory network. Front. Plant Sci. 2017, 8, 565. [Google Scholar] [CrossRef] [Green Version]
  69. Shamimuzzaman, M.; Vodkin, L. Identification of soybean seed developmental stage-specific and tissue-specific miRNA targets by degradome sequencing. BMC Genom. 2012, 13, 310. [Google Scholar] [CrossRef] [Green Version]
  70. Tang, J.; Chu, C. MicroRNAs in crop improvement: Fine-tuners for complex traits. Nat. Plants 2017, 3, 17077. [Google Scholar] [CrossRef]
  71. Zhang, Y.; Xia, R.; Kuang, H.; Meyers, B.C. The diversification of plant NBS-LRR defense genes directs the evolution of MicroRNAs that target them. Mol. Biol. Evol. 2016, 33, 2692–2705. [Google Scholar] [CrossRef] [Green Version]
  72. Meyers, B.C.; Morgante, M.; Michelmore, R.W. TIR-X and TIR-NBS proteins: Two new families related to disease resistance TIR-NBS-LRR proteins encoded in Arabidopsis and other plant genomes. Plant J. 2002, 32, 77–92. [Google Scholar] [CrossRef] [PubMed]
  73. Allen, E.; Xie, Z.; Gustafson, A.M.; Sung, G.H.; Spatafora, J.W.; Carrington, J.C. Evolution of microRNA genes by inverted duplication of target gene sequences in Arabidopsis thaliana. Nat. Genet. 2004, 36, 1282–1290. [Google Scholar] [CrossRef] [PubMed]
  74. De Felippes, F.F.; Schneeberger, K.; Dezulian, T.; Huson, D.H.; Weigel, D. Evolution of Arabidopsis thaliana microRNAs from random sequences. RNA 2008, 14, 2455–2459. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Xia, R.; Meyers, B.C.; Liu, Z.; Beers, E.P.; Ye, S.; Liu, Z. MicroRNA superfamilies descended from miR390 and their roles in secondary small interfering RNA biogenesis in eudicots. Plant Cell 2013, 25, 1555–1572. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Kinlaw, C.S.; Neale, D.B. Complex gene families in pine genomes. Trends Plant Sci. 1997, 2, 356–359. [Google Scholar] [CrossRef]
  77. Nystedt, B.; Street, N.R.; Wetterbom, A.; Zuccolo, A.; Lin, Y.C.; Scofield, D.G.; Vezzi, F.; Delhomme, N.; Giacomello, S.; Alexeyenko, A.; et al. The Norway spruce genome sequence and conifer genome evolution. Nature 2013, 497, 579–584. [Google Scholar] [CrossRef] [Green Version]
  78. Neale, D.B.; McGuire, P.E.; Wheeler, N.C.; Stevens, K.A.; Crepeau, M.W.; Cardeno, C.; Zimin, A.V.; Puiu, D.; Pertea, G.M.; Sezen, U.U.; et al. The Douglas-Fir genome sequence reveals specialization of the photosynthetic apparatus in Pinaceae. G3 (Bethesda) 2017, 7, 3157–3167. [Google Scholar] [CrossRef] [Green Version]
  79. Kavar, T.; Maras, M.; Kidrič, M.; Šuštar-Vozlič, J.; Meglič, V. Identification of genes involved in the response of leaves of Phaseolus vulgaris to drought stress. Mol. Breed. 2008, 21, 159–172. [Google Scholar] [CrossRef]
  80. Salih, H.; Gong, W.; He, S.; Sun, G.; Sun, J.; Du, X. Genome-wide characterization and expression analysis of MYB transcription factors in Gossypium hirsutum. BMC Genet. 2016, 17, 129. [Google Scholar] [CrossRef] [Green Version]
  81. Edger, P.P.; Pires, J.C. Gene and genome duplications: The impact of dosage-sensitivity on the fate of nuclear genes. Chromosome Res. 2009, 17, 699–717. [Google Scholar] [CrossRef] [Green Version]
  82. Sharma, N.; Bhalla, P.L.; Singh, M.B. Transcriptome-wide profiling and expression analysis of transcription factor families in a liverwort, Marchantia polymorpha. BMC Genom. 2013, 14, 915. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Li, C.; Zhang, A.B. MicroRNAs in control of plant development; MicroRNAs in control of plant development. J. Cell. Physiol. 2016, 231, 303–313. [Google Scholar] [CrossRef] [PubMed]
  84. Shriram, V.; Kumar, V.; Devarumath, R.M.; Khare, T.S.; Wani, S.H. Micrornas as potential targets for abiotic stress tolerance in plants. Front. Plant Sci. 2016, 7, 817. [Google Scholar] [CrossRef] [PubMed]
  85. Shu, Y.; Liu, Y.; Li, W.; Song, L.; Zhang, J.; Guo, C. Genome-wide investigation of MicroRNAs and their targets in response to freezing stress in Medicago sativa L., based on high-throughput sequencing. G3 (Bethesda) 2016, 6, 755–765. [Google Scholar] [CrossRef] [Green Version]
  86. D’Ario, M.; Griffiths-Jones, S.; Kim, M. Small RNAs: Big impact on plant development. Trends Plant Sci. 2017, 22, 1056–1068. [Google Scholar] [CrossRef] [Green Version]
  87. Meng, Y.; Shao, C.; Wang, H.; Chen, M. Are all the miRBase-registered microRNAs true? A structure- and expression-based re-examination in plants. RNA Biol. 2012, 9, 249–253. [Google Scholar] [CrossRef] [Green Version]
  88. Taylor, R.S.; Tarver, J.E.; Foroozani, A.; Donoghue, P.C.J. MicroRNA annotation of plant genomes—Do it right or not at all. BioEssays 2017, 39, 1–6. [Google Scholar] [CrossRef]
Table 1. Conserved small RNAs identified from miRBase.
Table 1. Conserved small RNAs identified from miRBase.
AnnotationCountPercentage
Annotated49750.5%
Acacia auriculiformis410.8%
Arabidopsis thaliana4579.2%
Oryza sativa3076.2%
Picea abies167633.7%
Pinus taeda145929.3%
Pinus densata58611.8%
Populus euphratica20.0%
Populus trichocarpa1623.3%
Nicotiana tabacum821.6%
Vitis vinifera1322.7%
Zea mays711.4%
Unannotated1,016,72199.5%
Total1,021,696100.0%
Table 2. Comparison of the number of P. sylvestris small RNAs with matching miRNA sequences in miRBase.
Table 2. Comparison of the number of P. sylvestris small RNAs with matching miRNA sequences in miRBase.
SpeciesNo. of Sequences in miRBaseNo. of Matching SequencesPercentage Found
Acacia auriculiformis7457.1%
Arabidopsis thaliana2989030.2%
Oryza sativa5928213.9%
Picea abies403075.0%
Pinus taeda363391.7%
Pinus densata302170.0%
Populus euphratica4125.0%
Populus trichocarpa3527320.7%
Nicotiana tabacum1623018.5%
Vitis vinifera1632716.6%
Zea mays1723721.5%
Table 3. miRNAs found only in conifers based on miRBase and Scots pine data.
Table 3. miRNAs found only in conifers based on miRBase and Scots pine data.
miRNA FamilyP.sylvestrisP.taedaP.densataP.abiesC.lanceolata
miR946+++--
miR948-+---
miR947++-+-
miR949++---
miR950++++-
miR952+++--
miR951++-+-
miR1309++---
miR1311++++-
miR1312+++--
miR1313+++--
miR1314+++--
miR1315++---
miR1316++---
miR3693---+-
miR3694---+-
miR3695---+-
miR3696---+-
miR3697---+-
miR3698---+-
miR3699+--+-
miR3700---+-
miR3701+-++-
miR3702+--+-
miR3703---+-
miR3704--++-
miR3705---+-
miR3706---+-
miR3707---+-
miR3708---+-
miR3709---+-
miR3710+--+-
miR3711---+-
miR3712+-++-

Share and Cite

MDPI and ACS Style

Krivmane, B.; Šņepste, I.; Šķipars, V.; Yakovlev, I.; Fossdal, C.G.; Vivian-Smith, A.; Ruņģis, D. Identification and in Silico Characterization of Novel and Conserved MicroRNAs in Methyl Jasmonate-Stimulated Scots Pine (Pinus sylvestris L.) Needles. Forests 2020, 11, 384. https://doi.org/10.3390/f11040384

AMA Style

Krivmane B, Šņepste I, Šķipars V, Yakovlev I, Fossdal CG, Vivian-Smith A, Ruņģis D. Identification and in Silico Characterization of Novel and Conserved MicroRNAs in Methyl Jasmonate-Stimulated Scots Pine (Pinus sylvestris L.) Needles. Forests. 2020; 11(4):384. https://doi.org/10.3390/f11040384

Chicago/Turabian Style

Krivmane, Baiba, Ilze Šņepste, Vilnis Šķipars, Igor Yakovlev, Carl Gunnar Fossdal, Adam Vivian-Smith, and Dainis Ruņģis. 2020. "Identification and in Silico Characterization of Novel and Conserved MicroRNAs in Methyl Jasmonate-Stimulated Scots Pine (Pinus sylvestris L.) Needles" Forests 11, no. 4: 384. https://doi.org/10.3390/f11040384

APA Style

Krivmane, B., Šņepste, I., Šķipars, V., Yakovlev, I., Fossdal, C. G., Vivian-Smith, A., & Ruņģis, D. (2020). Identification and in Silico Characterization of Novel and Conserved MicroRNAs in Methyl Jasmonate-Stimulated Scots Pine (Pinus sylvestris L.) Needles. Forests, 11(4), 384. https://doi.org/10.3390/f11040384

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop