Next Article in Journal
Molecular Cloning, Tissue Distribution and Antiviral Immune Response of Duck Src
Next Article in Special Issue
Micro RNA-175 Targets Claudin-1 to Inhibit Madin–Darby Canine Kidney Cell Adhesion
Previous Article in Journal
Molecular Manipulation of the miR160/AUXIN RESPONSE FACTOR Expression Module Impacts Root Development in Arabidopsis thaliana
Previous Article in Special Issue
Revealing Differential RNA Editing Specificity of Human ADAR1 and ADAR2 in Schizosaccharomyces pombe
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of the miR482 Gene Family in Plants

College of Agriculture and Forestry Sciences, Linyi University, Linyi 276000, China
*
Author to whom correspondence should be addressed.
Genes 2024, 15(8), 1043; https://doi.org/10.3390/genes15081043
Submission received: 19 May 2024 / Revised: 5 July 2024 / Accepted: 5 August 2024 / Published: 8 August 2024
(This article belongs to the Special Issue RNAs in Biology)

Abstract

:
MicroRNA482 (miR482) is a conserved microRNA family in plants, playing critical regulatory roles in different biological activities. Though the members of the miR482 gene family have been identified in plants, a systematic study has not been reported yet. In the present research, 140 mature sequences generated by 106 precursors were used for molecular characterization, phylogenetic analysis, and target gene prediction, and the competing endogenous RNA (ceRNA) network mediated by miR482 was summarized. The length of mature sequences ranged from 17 nt to 25 nt, with 22 nt being the most abundant, and the start and end of the mature sequences had a preference for uracil (U). By sequence multiplex comparison, it was found that the mature sequences of 5p were clustered into one group, and others were clustered into the other group. Phylogenetic analysis revealed that the 140 mature sequences were categorized into six groups. Meanwhile, all the precursor sequences formed a stable hairpin structure, and the 106 precursors were divided into five groups. However, the expression of miR482 varied significantly between different species and tissues. In total, 149 target genes were predicted and their functions focused on single-organism process, cellular process, and cell and cell part. The ceRNA network of miR482 in tomato, cotton, and peanut was summarized based on related publications. In conclusion, this research will provide a foundation for further understanding of the miR482 gene family.

1. Introduction

MicroRNAs (miRNAs) are endogenous small non-coding RNAs (ncRNA), 17–25 nts in length, which are widely distributed in both plants and animals [1]. They were first reported in the nematode Caenorhabditis elegans as early as 1993 [2]. However, the first miRNAs in plants were reported in 2002, with miRNAs discovered in Arabidopsis thaliana [3]. Since then, with the development and application of high-throughput sequencing, more and more miRNAs have been reported in radish [4], wild soybean [5], tomato [6], and so on. So far, miRNAs from 271 species, including nearly 200 plant genomes, have been studied.
miRNAs are regarded as gene expression regulators, which play key roles at the post-transcriptional level by binding to target mRNAs and combining with argonaute proteins to achieve the blockage of mRNA translation or degradation of target mRNAs [7,8]. Additionally, in some cases, miRNA could induce the generation of secondary phased siRNAs (phasiRNAs) from the target genes after the cleavage [9]. miRNAs are involved in various biological processes, ranging from plant growth and development to biotic and abiotic stresses [10,11]. In particular, their function in the responses to different stress has attracted more attention. For example, apple miR160 can improve drought resistance by promoting the development of roots and rhizomes [12]. The knockdown lines of soybean miR166 displays reduced gibberellin content and plant dwarfing [13]. miR398 regulates the expression of target genes (the Cu/Zn superoxide dismutase gene family) to adapt to a variety of adverse stresses [14]. miR482 inhibits the nucleotide-binding site–leucine-rich repeat (NBS-LRR) genes, which are involved in resistance to pathogens [15].
Long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) are other kinds of non-coding RNAs. By interacting with miRNA response elements (MREs), they can act as competing endogenous RNAs (ceRNAs) to regulate the expression of target genes with miRNAs [16,17,18]. This regulatory model reveals the complex ceRNAs network formed by different kinds of non-coding RNAs mediated by miRNAs. Currently, the regulatory mechanism of ceRNAs is a hotspot for further research on miRNAs.
miR482 is a class of miRNAs distributed in various plants. The first plant that miR482 was identified in was Populus trichocarpa, while it was absent from Arabidopsis, and it was found that miR482 was conserved between rice and Populus by searching the rice genome [19]. It is worth mentioning that miR2118 has highly similar sequences to miR482, and together with miR482, they are now known as the miR482/2118 superfamily. However, there are differences between their beginning sequences; miR482 expands more widely in dicots while miR2118 expands in monocots, and they have different functions: miR482 is mainly responsible for anti-disease functions, and miR2118 is mainly related to growth and development [20]. Afterward, it was successively reported in soybean [21], apple [22], tomato [23], and other plants. The regulation modes of miR482 vary among different plants, making its function interesting. For example, miR482 in black cottonwood targets the gene encoding a resistance protein [19]. It also directly targets the NBS-LRR gene in tomato, engaging in the late blight response process [24]. In cotton, miR482 targets the NLR gene, playing a role in the response to dahlia yellow mosaic disease [25]. Additionally, miR482/2118 in litchi cleaves a long non-coding RNA gene, LcTASL1, and induces the production of phasiRNAs, which is involved in GA signaling [26]. These findings suggest that miR482 plays a positive role in regulating multiple stress responses [24,27,28]. Therefore, research on the molecular features, evolutionary patterns, target gene functions, and ceRNAs network of the miR482 family in plants is very meaningful, and could lay a theoretical foundation for further understanding of the genetic evolutionary structure and functional studies of the plant miR482 family.

2. Materials and Methods

2.1. Data Preparation of miR482 Family Members for Bioinformatic Analysis

For the identification of the miR482 family in plants, the mature and precursor sequences of the miR482 gene family in plants were respectively downloaded from the miRBase database server v. 22.0 (http://www.mirbase.org/, accessed on 10 August 2023) [29]. The sequences were used for the following bioinformatic analysis.

2.2. Molecular Characterization of miR482 Family Members

The mature and precursor sequences of the miR482 family were classified and analyzed, respectively. Based on the mature sequences, the sequences were divided into two groups, with all the 5p sequences forming one group and the rest forming the other. The two groups were analyzed for base preference separately on the MEME website (https://meme-suite.org/meme/tools/meme, accessed on 4 July 2024) [30]. According to the minimum number of mature bases, a minimum of 17 and 19 motifs were set for base preference analysis, and other settings according to system defaults. Meanwhile, the precursors were submitted to the RNAfold website (http://rna.tbi.univie.ac.at//cgi-bin/RNAWebSuite/RNAfold.cgi, accessed on 22 August 2023) for secondary structure prediction.

2.3. Sequence Alignment and Phylogenetic Tree Estimation

The multiple alignment of the miR482 gene family was conducted by MEGA software (version 11) using the ClustalW method. And MEGA software was utilized to align mature and precursor sequences of the plant miR482 family, and then to construct the phylogenetic tree based on the neighbor-joining (NJ) method [31]. Ultrafast bootstrapping was carried out 1000 times, and other parameters were set to default [32].

2.4. Expression Profiling of miR482

To analyze the expression level of miR482 in different tissues of different species, the expression data of miR482 in Citrus sinensis, Nicotiana tabacum, Solanum lycopersicum, Solanum tuberosum, Fragaria vesca, Glycine max, Manihot esculenta, and Prunus persica were obtained from PmiREN (Plant miRNA ENcyclopedia) and used for comparison.

2.5. Target Gene Prediction and GO Enrichment Analysis

The plant small RNA target analysis server tool (psRNATarget, https://www.zhaolab.org/psRNATarget/home, accessed on 24 August 2023) was employed to predict the potential target genes of miR482 in plants [33]. The analysis was conducted under a strict setting where the maximum expected value was set to 1.5. Then, the National Center for Biotechnology Information (NCBI, https://www.ncbi.nlm.nih.gov/, accessed on 28 August 2023) was used for an online comparison to obtain the gene function annotation. The EGGNOG-MAPPER website (http://eggnog-mapper.embl.de, accessed on 4 July 2024) was employed for the GO annotation of predicted target genes. And the results of the GO annotation were visualized using the online omicshare tool (https://www.omicshare.com/tools/Home/Soft/gogseasenior, accessed on 4 July 2024).

2.6. CeRNA Network Construction

To further analyze the function of plant miR482, we systematically summarized related research on the regulation of miR482 by ncRNAs described in the current literature. The ceRNA network of miR482 was based on the following publications: Liu et al. [34], Liu et al. [35], Xu et al. [36], Si et al. [37], and Li et al. [38]. And the ceRNA network of miR482 was constructed using Cytoscape (version 3.9.1).

3. Results

3.1. Statistical Analysis of Members of the miR482 Family in Plants

During the accession into the miRbase database, the members of the miR482 family in plants were downloaded (Table S1). In total 140 mature sequences from 26 species belonging to 14 families were found, as shown in Table 1 and Figure 1. The largest miR482 family was Picea abies with 25 members including 24 different submembers ranging from miR482a to miR482x, followed by C. sinensis with 14, including 7 different submembers, while only 1 member was in Medicago truncatula, Phaseolus vulgaris, Vigna unguiculata, Vitis vinifera, and Zea mays. Among these sequences, 99 mature sequences of miR482 were located at the 3′ end of the precursor, while 41 mature sequences were at the 5′ end. Meanwhile, 91 mature sequences were experimentally verified, 8 mature sequences were obtained by homologous comparison, and 41 mature sequences remained to be verified.

3.2. Molecular Characterization of the Mature and Precursor Sequences of the miR482 Family

The distribution of the length of mature sequences in the miR482 family is shown in Figure 2. The length of miR482 mature sequences was mainly within 20–23 nt. And the sequences of 22 nt in length accounted for 69.29% and were the most abundant, which is basically in line with the typical miRNA maturation length characteristics of plants.
A sequence multiplex comparison of the mature sequences is shown in Figure 3. Most mature sequences were highly conserved. These sequences could be categorized into two groups based on base characterization and conservatism, where the mature sequences produced at the 5′ end were clustered into group II and the rest were clustered into group I.
The 140 mature sequences of the miR482 family in plants were divided into two groups utilized for base preference analysis at each position, and the results are shown in Figure 4. In group I, a preference for uracil (U) was observed in the 1st, 2nd, and 17th positions, whereas the 4th, 5th, 10th, 11th, 13th, 14th, and 15th positions showed a preference for cytosine (C). In group II, a preference for uracil (U) was observed only in the 3rd position, whereas the 4th, 5th, 6th, 8th, 13th, 14th, and 18th positions showed a preference for guanine (G).
Based on the precursor sequences of the miR482 family, the secondary structure of the precursor sequences was predicted via the RNAfold website. All the precursor sequences of the miR482 family could form the stable hairpin structure (Figure S1). The number of loop structures formed by each precursor was not the same, the miR482 mature sequences could be generated from either 5p or 3p, and the mature sequences might also include some of the loop sequences and stem sequences.

3.3. Phylogenetic Analysis of the miR482 Family in Plants

To explore the evolutionary relationships of the miR482 family members in plants, a phylogenetic tree based on the alignments of 140 mature sequences was constructed (Figure 5). All the mature sequences could be divided into six groups, from group I to group VI. Group VI (Figure 5, blue branches) contained the most mature sequences, including 99 mature sequences of the miR482 family. Groups IV and V contained 16 and 22 sequences, respectively. And only 1 sequence was in groups I, II, and III, respectively. All miR482-3p mature sequences were clustered into group VI, while miR482-5p mature sequences were widely distributed in other groups.
Meanwhile, phylogenetic analysis was also performed based on the 106 precursor sequences of the plant miR482 family, and the results are shown in Figure 6. Different from the phylogenetic tree of mature sequences, all the precursors were divided into five groups, from group I to group V. Of these, group V contained the largest number of precursor sequences, including 89 precursor sequences (Figure 6, green branches). Groups III and IV contained eight and seven sequences, respectively. Groups I and II included only one precursor sequence.

3.4. Tissue-Specific Expression Analysis of miR482

Based on the expression information obtained from PmiREN, the tissue-specific expression patterns of miR482 in C. sinensis, N. tabacum, S. lycopersicum, S. tuberosum, F. vesca, G. max, M. esculenta, and P. persica were analyzed (Figure 7). As expected, the expression of miR482 varied significantly between different species and tissues. miR482 was highly expressed in C. sinensis, F. vesca, G. max, M. esculenta, and P. persica. A low expression level of miR482 was found in N. tabacum, S. lycopersicum, and S. tuberosum. And the expression differences among different members of miR482 in the same species were also large. In C. sinensis, the expression level of miR482e in root was 6646 reads per million (RPM), while the expression level of miR482b in root was only 4 RPM. Significant differences also existed in other species. Compared to other members, miR482a was highly expressed in N. tabacum and F. vesca, miR482e was highly expressed in C. sinensis and S. lycopersicum, miR482d was highly expressed in S. tuberosum and G. max, while miR482b was highly expressed in M. esculenta and P. persica.

3.5. Prediction and Function Annotation of the Targets

To further understand the function of miR482, target gene prediction was performed (Table S2). Of these 140 mature sequences, 149 targets from 63 mature sequences were found across 16 species. miR482 inhibits target genes mainly by cleavage. And multiple members of the miR482 family could target the same gene; for example, fve-miR482a and fve-miR482c co-targeted Garcinia cambogia subspecies sphingosine-like kinase 1 (gene30692-v2.0.a2-hybrid). The 3p and 5p of the miR482 produced from the same precursor also targeted the same gene. The sweet orange csi-miR482f-5p and csi-miR482f-3p co-targeted the citrus α-glucan water bisphosphatase 2 gene (DSWO2A02249).
And the functional annotation and GO enrichment analysis of the target genes were conducted using the online website EGGNOG-MAPPER at the same time (Figure 8). The most enriched GO terms of those target genes were in the biological process category, which comprised single-organism process and cellular process. Within the cellular component category, cell and cell part was enriched. And, in the molecular function category, binding was enriched.

3.6. CeRNA Network of miR482

By integrating the existing related reports, the ceRNA network of miR482 was constructed (Figure 9, Table S3). Tomato sllncRNA15492 and sllncRNA23468 acted as the ceRNAs for miR482a and miR482b, respectively, by regulating NBS-LRR in response to Phytophthora infestans infection (Figure 9a). In addition, sllncRNA08489 and sllncRNA39298 could also act as ceRNAs to decoy slmiR482e-3p and slmiR482e-5p, respectively, in response to tomato P. infestans, which further regulated the expression of NBS-LRR (Figure 9a). Furthermore, CC-NBS-LRR was found to be co-regulated by miR482c-3p and miR482e-3p, which were the targets of five lncRNAs (TCONS-00026121, TCONS-00061862, TCONS-00026124, TCONS-00061875, TCONS-00061867) and three lncRNAs (TCONS-00057529, TCONS-00113349, TCONS-00057528), respectively, in responding to Ralstonia solanacearum (Figure 9b). And the specific lncRNA (TCONS_00003967)–miR482g–mRNA regulatory network was found in cotton response to drought stress (Figure 9c). Additionally, when peanut was infected by root-knot nematodes, four lncRNAs (MSTRG.2115, MSTRG.30601, MSTRG.30599, and MSTRG.31962) and a circRNA (circRNA320) acted as the ceRNAs of miR482c (Figure 9d).

4. Discussion

miR482 is widely present in seed plants, and has been reported in a variety of plants including tomato, cotton, and apple [24,39,40]. In this study, 140 mature sequences of the miR482 family from 26 species were obtained and used for analysis. More conserved miRNA families usually have more members [41]. For example, 14 members of the miR482 family were discovered in cotton involved in root-knot nematode Meloidogyne incognita infection [42], and there were six members of miR482 in Chrysanthemums under salt stress [41]. In this research, 14 out of 26 species had more than three members of the miR482 family, and in particular, P. abies contained 25 members, including 24 different submembers ranging from miR482a to miR482x. The expression of miR482 varied significantly between different species and tissues, and the expression of different members of miR482 in the same species also varied widely. For example, miR482e showed a high expression level in C. sinensis, while a low expression level in P. persica. In F. vesca, miR482e had a higher expression level in seed; however, it was not detected in fruit. Many miRNAs are conserved among plants. In our study, 140 mature sequences of the miR482 family were used for analysis. The sequence analysis showed that they were conserved, and the preference for uracil (U) was observed in the start and end bases of the sequences.
According to the phylogenetic tree, all the precursor sequences of the miR482 family were clustered into five groups, and the mature sequences were divided into six groups. The diverse distributional features were independent of the species itself in terms of their affinity and evolutionary proximity (Figure 5 and Figure 6). This is consistent with previous reports of the evolutionary analysis of miR156 [43]. The miR482 precursor sequences of the same species were distributed in different groups. The precursor sequences of miR482 in P. abies (pab-MIR482) were clustered in groups III and V, while the mature sequences of miR482 in P. abies were clustered in groups IV, V, and VI. This indicates the non-synchronization of the evolutionary rate of different miR482 family member sequences within the same species. All the precursor sequences of miR482 in A. officinalis were clustered in group V. Correspondingly, all the mature sequences were gathered in group VI in A. officinalis, as in Pinus taeda, which indicates that the sequences of miR482 in A. officinalis and P. taeda have recently diversified. Besides this, csi-miR482b-5p independently formed group III, and the expression levels of csi-miR482b in C. sinensis were significantly lower than those of other members. It is hypothesized that there may be some correlation between differences in expression and evolutionary branching.
Previous studies have found that the length of an miRNA mature sequence varies among different species. The lengths of miRNA in poplar [44], ginseng [45], and cotton [46] were at most 21 nt, whereas those of blueberry [47], sweet potato [48], and tea tree [49] were at most 24 nt. And it was reported that the length was related to different functions. miRNAs 21–22 nt in length are mainly associated with mRNA cleavage and gene silencing at the post-transcriptional level [50], while miRNAs 24 nt in length are mainly related to RNA-guided DNA methylation [51]. In this study, the length of miR482 members in plants was at most 21–22 nt, and the mode of regulation between miRNA and target genes was mainly cleavage, which was consistent with our above conclusion.
miR482 is involved in the regulation of stress-related biological activities and mainly regulates NBS-LRR genes [52]. For example, apple miR482 targets NBS-LRR in response to Alternaria alternata f. sp. mali infection [40]. In tomato, miR482a, miR482b, miR482e-3p, and miR482e-5p were all involved in the regulation of P. infestans infection by targeting NBS-LRR, while lncRNA15492, lncRNA23468, lncRNA08489, and lncRNA39298 acted as their target ceRNAs [24,34,35,52]. Additionally, infected with R. solanacearum, miR482c-3p and miR482e-3p targeting CC-NBS-LRR were both responsive, and eight lncRNAs acted as the target ceRNAs [37]. Then, a ceRNA network was constructed (Figure 9a,b). Of these, the target genes of miR482, NBS-LRR genes, are the largest family of resistance genes (R genes) and are widely reported in plants. Also, a similar regulation ceRNA network was found in cotton in response to drought stress, and in peanut after infection with root-knot nematodes (Figure 9c,d). Apart from NBS-LRR, the predicted target genes of miR482 were also found to participate in response to stimulus, binding, and cellular anatomical entity. Lipoxygenases (LOXs) are dioxygenases without heme and iron and are involved in the development and adaptation of many plants to the environment. Strawberry LOX genes are involved in low-temperature, drought, salt, and other abiotic stress responses [53]. In our research, LOXs were targeted by fve-miR482b and fve-miR482c. α-glucan water dikinases are cytosolic enzymes that are important for proper A. thaliana seed development and involved in cold tolerance [54,55], which was the predicted target gene of csi-miR482b-5p and csi-miR482b-3p. Taken together, the miR482 gene family plays a critical role in various biological processes. However, the role of the miR482 gene family remains to be further confirmed by experiments.

5. Conclusions

In this study, the miR482 gene family members, corresponding predicted target genes, and the ceRNA network in plants were comprehensively analyzed and summarized. This research offers an extensive characterization of miR482 in plants, and also provides valuable references for the study of miR482-mediated plant responses to adversity stress.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/genes15081043/s1. Figure S1: Secondary structure of miR482 precursor sequences in plants. Table S1: The mature and precursor sequences of the miR482 family in plants. Table S2: Target genes of miR482 and function annotation. Table S3: ceRNA network of miR482.

Author Contributions

X.C. and M.X. conceived and designed the research plans; W.K., D.Q. and M.X. analyzed the data; W.K. wrote the manuscript; Y.H., Y.L. and X.C. revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (31801868; 32302621, 32102405) and Natural Science Foundation of Shandong Province, China (ZR2018PC023, ZR2021QC100, ZR2020QC156).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Materials. Further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Bartel, D.P. MicroRNAs: Genomics, biogenesis, mechanism, and function. Cell 2004, 116, 281–297. [Google Scholar] [CrossRef] [PubMed]
  2. Lee, R.C.; Feinbaum, R.L.; Ambros, V. The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell 1993, 75, 843–854. [Google Scholar] [CrossRef] [PubMed]
  3. Llave, C.; Xie, Z.; Kasschau, K.D.; Carrington, J.C. Cleavage of Scarecrow-like mRNA targets directed by a class of Arabidopsis miRNA. Science 2002, 297, 2053–2056. [Google Scholar] [CrossRef]
  4. Xu, L.; Wang, Y.; Xu, Y.Y.; Wang, L.J.; Zhai, L.L.; Zhu, X.W.; Gong, Y.Q.; Ye, S.; Liu, L.W. Identification and characterization of novel and conserved microRNAs in radish (Raphanus sativus L.) using high-throughput sequencing. Plant Sci. 2013, 201–202, 108–114. [Google Scholar] [CrossRef] [PubMed]
  5. Zeng, Q.Y.; Yang, C.Y.; Ma, Q.B.; Li, X.P.; Dong, W.W.; Nian, H. Identification of wild soybean miRNAs and their target genes responsive to aluminum stress. BMC Plant Biol. 2012, 12, 182. [Google Scholar] [CrossRef] [PubMed]
  6. Cao, X.; Wu, Z.; Jiang, F.L.; Zhou, R.; Yang, Z.N. Identification of chilling stress-responsive tomato microRNAs and their target genes by high-throughput sequencing and degradome analysis. BMC Genom. 2014, 15, 1130. [Google Scholar] [CrossRef] [PubMed]
  7. Fabian, M.R.; Sonenberg, N.; Filipowicz, W. Regulation of mRNA translation and stability by microRNAs. Annu. Rev. Biochem. 2010, 79, 351–379. [Google Scholar] [CrossRef] [PubMed]
  8. Khan, Z.; Suthanthiran, M.; Muthukumar, T. MicroRNAs and transplantation. Clin. Lab. Med. 2019, 39, 125–143. [Google Scholar] [CrossRef] [PubMed]
  9. Fei, Q.; Xia, R.; Meyers, B.C. Phased, secondary, small interfering RNAs in posttranscriptional regulatory networks. Plant Cell 2013, 25, 2400–2415. [Google Scholar] [CrossRef] [PubMed]
  10. Xu, C.Z.; Tao, Y.X.; Fu, X.K.; Guo, L.; Xing, H.T.; Li, C.F.; Yang, Z.W.; Su, H.L.; Wang, X.Q.; Hu, J.; et al. The microRNA476a-RFL module regulates adventitious root formation through a mitochondria-dependent pathway in Populus. New Phytol. 2021, 230, 2011–2028. [Google Scholar] [PubMed]
  11. Liu, Q.G.; Wang, Z.C.; Yu, S.; Li, W.L.; Zhang, M.Q.; Yang, J.; Li, D.D.; Yang, J.L.; Li, C.H. Pu-miR172d regulates stomatal density and water-use efficiency via targeting PuGTL1 in poplar. J. Exp. Bot. 2021, 72, 1370–1383. [Google Scholar] [CrossRef] [PubMed]
  12. Shen, X.X.; He, J.Q.; Ping, Y.K.; Guo, J.X.; Hou, N.; Cao, F.G.; Li, X.W.; Geng, D.L.; Wang, S.C.; Chen, P.X.; et al. The positive feedback regulatory loop of miR160-Auxin Response Factor 17-HYPONASTIC LEAVES 1 mediates drought tolerance in apple trees. Plant Physiol. 2022, 188, 1686–1708. [Google Scholar] [CrossRef] [PubMed]
  13. Zhao, C.; Ma, J.J.; Zhang, Y.H.; Yang, S.X.; Feng, X.Z.; Yan, J. The miR166 mediated regulatory module controls plant height by regulating gibberellic acid biosynthesis and catabolism in soybean. J. Integr. Plant Biol. 2022, 64, 995–1006. [Google Scholar] [CrossRef] [PubMed]
  14. Li, J.; Song, Q.Q.; Zuo, Z.F.; Liu, L. MicroRNA398: A master regulator of plant development and stress responses. Int. J. Mol. Sci. 2022, 23, 10803. [Google Scholar] [CrossRef] [PubMed]
  15. Xiao, Y.; Luan, Y.S. MiR482 mediates the biotic stress response of plant. J. Plant Physiol. 2014, 50, 741–748. [Google Scholar]
  16. An, Y.; Furber, K.L.; Ji, S.P. Pseudogenes regulate parental gene expression via ceRNA network. J. Cell. Mol. Med. 2017, 21, 185–192. [Google Scholar] [CrossRef] [PubMed]
  17. Bosson, A.D.; Zamudio, J.R.; Sharp, P.A. Endogenous miRNA and target concentrations determine susceptibility to potential ceRNA competition. Mol. Cell 2014, 56, 347–359. [Google Scholar] [CrossRef] [PubMed]
  18. Yamamura, S.; Imai-Sumida, M.; Tanaka, Y.; Dahiya, R. Interaction and cross-talk between non-coding RNAs. Cell. Mol. Life Sci. 2018, 75, 467–484. [Google Scholar] [CrossRef] [PubMed]
  19. Lu, S.F.; Sun, Y.H.; Shi, R.; Clark, C.; Li, L.G.; Chiang, V.L. Novel and mechanical stress-responsive MicroRNAs in Populus trichocarpa that are absent from Arabidopsis. Plant Cell 2005, 17, 2186–2203. [Google Scholar] [CrossRef] [PubMed]
  20. Zhang, Y.Q.; Waseem, M.; Zeng, Z.H.; Xu, J.; Chen, C.J.; Liu, Y.L.; Zhai, J.X.; Xia, R. MicroRNA482/2118, a miRNA superfamily essential for both disease resistance and plant development. New Phytol. 2022, 233, 2047–2057. [Google Scholar] [CrossRef] [PubMed]
  21. Subramanian, S.; Fu, Y.; Sunkar, R.; Barbazuk, W.B.; Zhu, J.K.; Yu, O. Novel and nodulation-regulated microRNAs in soybean roots. BMC Genom. 2008, 9, 160. [Google Scholar] [CrossRef] [PubMed]
  22. Arenas-Huertero, C.; Pérez, B.; Rabanal, F.; Blanco-Melo, D.; De la Rosa, C.; Estrada-Navarrete, G.; Sanchez, F.; Covarrubias, A.A.; Reyes, J.L. Conserved and novel miRNAs in the legume Phaseolus vulgaris in response to stress. Plant Mol. Biol. 2009, 70, 385–401. [Google Scholar] [CrossRef] [PubMed]
  23. Mohorianu, I.; Schwach, F.; Jing, R.C.; Lopez-Gomollon, S.; Moxon, S.; Szittya, G.; Sorefan, K.; Moulton, V.; Dalmay, T. Profiling of short RNAs during fleshy fruit development reveals stage-specific sRNAome expression patterns. Plant J. 2011, 67, 232–246. [Google Scholar] [CrossRef] [PubMed]
  24. Jiang, N.; Meng, J.; Cui, J.; Sun, G.X.; Luan, Y.S. Function identification of miR482b, a negative regulator during tomato resistance to Phytophthora infestans. Hortic. Res. 2018, 5, 9. [Google Scholar] [CrossRef] [PubMed]
  25. Zhu, Q.H.; Jin, S.X.; Yuan, Y.M.; Liu, Q.; Zhang, X.L.; Wilson, I. CRISPR/Cas9-mediated saturated mutagenesis of the cotton MIR482 family for dissecting the functionality of individual members in disease response. Plant Direct 2022, 6, e410. [Google Scholar] [CrossRef] [PubMed]
  26. Zhang, Y.Q.; Zeng, Z.H.; Hu, H.M.; Zhao, M.L.; Chen, C.J.; Ma, X.S.; Li, G.L.; Li, J.G.; Liu, Y.L.; Hao, Y.W.; et al. MicroRNA482/2118 is lineage-specifically involved in gibberellin signalling via the regulation of GID1 expression by targeting noncoding PHAS genes and subsequently instigated phasiRNAs. Plant Biotechnol. J. 2024, 22, 819–832. [Google Scholar] [CrossRef] [PubMed]
  27. Qi, J.Y. Research on Molecular Mechanism of Sly-miR482a Negatively Regulating Plant Fruit Development and Resistance to Botrytis cinerea. Master’s Thesis, Zhejiang Sci-Tech University, Hangzhou, China, 2021. [Google Scholar]
  28. Chen, K. Functional Characterization of the miR482/2118 Family and Its Target Genes in Cotton Resistance to Verticillium wilt. Master’s Thesis, Huazhong Agricultural University, Wuhan, China, 2023. [Google Scholar]
  29. Kozomara, A.; Birgaoanu, M.; Griffiths-Jones, S. miRBase: From microRNA sequences to function. Nucleic Acids Res. 2019, 47, D155–D162. [Google Scholar] [CrossRef] [PubMed]
  30. Bailey, T.L.; Johnson, J.; Grant, C.E.; Noble, W.S. The MEME suite. Nucleic Acids Res. 2015, 43, W39–W49. [Google Scholar] [CrossRef] [PubMed]
  31. Tamura, K.; Stecher, G.; Kumar, S. MEGA11: Molecular evolutionary genetics analysis Version 11. Mol. Biol. Evol. 2021, 38, 3022–3027. [Google Scholar] [CrossRef] [PubMed]
  32. Kumar, S.; Nei, M.; Dudley, J.; Tamura, K. MEGA: A biologist-centric software for evolutionary analysis of DNA and protein sequences. Brief. Bioinform. 2008, 9, 299–306. [Google Scholar] [CrossRef] [PubMed]
  33. Dai, X.; Zhuang, Z.; Zhao, P.X. psRNATarget: A plant small RNA target analysis server (2017 release). Nucleic Acids Res. 2018, 46, W49–W54. [Google Scholar] [CrossRef] [PubMed]
  34. Liu, W.W.; Cui, J.; Luan, Y.S. Overexpression of lncRNA08489 enhances tomato immunity against Phytophthora infestans by decoying miR482e-3p. Biochem. Biophys. Res. Commun. 2022, 587, 36–41. [Google Scholar] [CrossRef] [PubMed]
  35. Liu, W.W.; Cui, J.; Luan, Y.S. SlmiR482e-5p regulates tomato resistance to Phytophthora infestans infection along with slmiR482e-3p via sllncRNA39298-mediated inhibition. Physiol. Mol. Plant Pathol. 2022, 121, 101875. [Google Scholar] [CrossRef]
  36. Xu, P.; Li, H.; Wang, X.H.; Zhao, G.; Lu, X.F.; Dai, S.J.; Cui, X.Y.; Yuan, M.; Liu, Z.N. Integrated analysis of the lncRNA/circRNA-miRNA-mRNA expression profiles reveals novel insights into potential mechanisms in response to root-knot nematodes in peanut. BMC Genom. 2022, 23, 239. [Google Scholar] [CrossRef] [PubMed]
  37. Si, X.Y.; Liu, H.Y.; Cheng, X.; Xu, C.C.; Han, Z.H.; Dai, Z.R.; Wang, R.Q.; Pan, C.T.; Lu, G. Integrative transcriptomic analysis unveils lncRNA-miRNA-mRNA interplay in tomato plants responding to Ralstonia solanacearum. Int. J. Biol. Macromol. 2023, 253, 126891. [Google Scholar] [CrossRef] [PubMed]
  38. Li, B.Q.; Feng, C.; Zhang, W.H.; Sun, S.M.; Yue, D.D.; Zhang, X.L.; Yang, X.Y. Comprehensive non-coding RNA analysis reveals specific lncRNA/circRNA-miRNA-mRNA regulatory networks in the cotton response to drought stress. Int. J. Biol. Macromol. 2023, 253, 126558. [Google Scholar] [CrossRef] [PubMed]
  39. Wu, P.; Lu, C.Z.; Wang, B.T.; Zhang, F.Y.; Shi, L.F.; Xu, Y.J.; Chen, A.M.; Si, H.J.; Su, J.J.; Wu, J.H. Cotton RSG2 mediates plant resistance against Verticillium dahliae by miR482b regulation. Biology 2023, 12, 898. [Google Scholar] [CrossRef] [PubMed]
  40. Liu, F.Y.; Tang, J.Q.; Li, T.Z.; Zhang, Q.L. The microRNA miR482 regulates NBS-LRR genes in response to ALT1 infection in apple. Plant Sci. 2024, 343, 112078. [Google Scholar] [CrossRef] [PubMed]
  41. Nai, J.F.; Ma, T.M.; Liu, Y.J.; Zhou, Y.W. Identification of small RNAs associated with salt stress in Chrysanthemums through high-throughput sequencing and bioinformatics analysis. Genes 2023, 14, 561. [Google Scholar] [CrossRef] [PubMed]
  42. Cai, C.P.; Li, C.; Sun, R.R.; Zhang, B.H.; Nichols, R.L.; Hake, K.D.; Pan, X.P. Small RNA and degradome deep sequencing reveals important roles of microRNAs in cotton (Gossypium hirsutum L.) response to root-knot nematode Meloidogyne incognita infection. Genomics 2021, 113, 1146–1156. [Google Scholar] [CrossRef]
  43. Li, X.K.; Liu, H.F.; Liu, J.L.; Zheng, M.; Hua, W. Sequence identification and bioinformatics of bna-miR156 gene family and target genes in rapeseed (Brassica napus L.). Chin. J. Oil Crop Sci. 2018, 40, 162. [Google Scholar]
  44. Li, H.; Dai, X.R.; Huang, X.; Xu, M.X.; Wang, Q.; Yan, X.J.; Sederoff, R.R.; Li, Q.Z. Single-cell RNA sequencing reveals a high-resolution cell atlas of xylem in Populus. J. Integr. Plant Biol. 2021, 63, 1906–1921. [Google Scholar] [CrossRef] [PubMed]
  45. Wang, Y.F.; Chen, Y.L.; Peng, M.Y.; Yang, C.; Yang, Z.M.; Gong, M.J.; Yin, Y.Q.; Zeng, Y. Identification of microRNA and analysis of target genes in Panax ginseng. Chin. Herb. Med. 2022, 15, 69–75. [Google Scholar] [CrossRef] [PubMed]
  46. Lu, X.K.; Yin, Z.J.; Wang, J.J.; Chen, X.G.; Wang, D.L.; Wang, S.; Guo, L.X.; Fan, W.L.; Chen, C.; Wang, X.G.; et al. Identification and function analysis of drought-specific small RNAs in Gossypium hirsutum L. Plant Sci. 2019, 280, 187–196. [Google Scholar] [CrossRef] [PubMed]
  47. Li, G.P.; Wang, Y.; Lou, X.M.; Li, H.L.; Zhang, C.Q. Identification of Blueberry miRNAs and their targets based on high-throughput sequencing and degradome analyses. Int. J. Mol. Sci. 2018, 19, 983. [Google Scholar] [CrossRef] [PubMed]
  48. Yang, Z.M.; Zhu, P.P.; Kang, H.; Liu, L.; Cao, Q.H.; Sun, J.; Dong, T.T.; Zhu, M.K.; Li, Z.Y.; Xu, T. High-throughput deep sequencing reveals the important role that microRNAs play in the salt response in sweet potato (Ipomoea batatas L.). BMC Genom. 2020, 21, 164. [Google Scholar] [CrossRef] [PubMed]
  49. Jeyaraj, A.; Liu, S.R.; Zhang, X.; Zhang, R.; Shangguan, M.Z.; Wei, C.L. Genome-wide identification of microRNAs responsive to Ectropis oblique feeding in tea plant (Camellia sinensis L.). Sci. Rep. 2017, 7, 13634. [Google Scholar] [CrossRef] [PubMed]
  50. Brodersen, P.; Voinnet, O. The diversity of RNA silencing pathways in plants. Trends Genet. 2006, 22, 268–280. [Google Scholar] [CrossRef] [PubMed]
  51. Wu, L.; Zhou, H.Y.; Zhang, Q.Q.; Zhang, J.G.; Ni, F.R.; Liu, C.; Qi, Y.J. DNA methylation mediated by a microRNA pathway. Mol. Cell. 2010, 38, 465–475. [Google Scholar] [PubMed]
  52. Jiang, N.; Cui, J.; Hou, X.X.; Yang, G.L.; Xiao, Y.; Han, L.; Meng, J.; Luan, Y.S. Sl-lncRNA15492 interacts with Sl-miR482a and affects Solanum lycopersicum immunity against Phytophthora infestans. Plant J. 2020, 103, 1561–1574. [Google Scholar] [CrossRef] [PubMed]
  53. Li, Z.Q.; Xie, Q.; Yan, J.H.; Chen, J.Q.; Chen, Q.X. Genome-wide identification and characterization of the abiotic-stress-responsive GRF gene family in diploid woodland strawberry (Fragaria vesca). Plants 2021, 10, 1916. [Google Scholar] [CrossRef] [PubMed]
  54. Pirone, C.; Gurrieri, L.; Gaiba, I.; Adamiano, A.; Valle, F.; Trost, P.; Sparla, F. The analysis of the different functions of starch-phosphorylating enzymes during the development of Arabidopsis thaliana plants discloses an unexpected role for the cytosolic isoform GWD2. Physiol. Plant 2017, 160, 447–457. [Google Scholar] [CrossRef] [PubMed]
  55. Yano, R.; Nakamura, M.; Yoneyama, T.; Nishida, I. Starch-related α-glucan/water dikinase is involved in the cold-induced development of freezing tolerance in Arabidopsis. Plant Physiol. 2005, 138, 837–846. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Presence and distribution of the miR482 family in plants.
Figure 1. Presence and distribution of the miR482 family in plants.
Genes 15 01043 g001
Figure 2. Length distribution of the miR482 family in plants.
Figure 2. Length distribution of the miR482 family in plants.
Genes 15 01043 g002
Figure 3. Multiple-sequence comparison of mature sequences in the miR482 family: (I) comparison of the 3p and non 5p or 3p miR482 sequence; (II) comparison of the 5p sequences.
Figure 3. Multiple-sequence comparison of mature sequences in the miR482 family: (I) comparison of the 3p and non 5p or 3p miR482 sequence; (II) comparison of the 5p sequences.
Genes 15 01043 g003
Figure 4. Base preference analysis of the mature sequences in the miR482 family. (I) 3p sequences and non 5p or 3p miR482 sequence were analyzed for base preference as a group; (II) 5p sequences were analyzed for base preference as a group.
Figure 4. Base preference analysis of the mature sequences in the miR482 family. (I) 3p sequences and non 5p or 3p miR482 sequence were analyzed for base preference as a group; (II) 5p sequences were analyzed for base preference as a group.
Genes 15 01043 g004
Figure 5. Phylogenetic relationship and classification of the mature sequences from the plant miR482 family.
Figure 5. Phylogenetic relationship and classification of the mature sequences from the plant miR482 family.
Genes 15 01043 g005
Figure 6. Phylogenetic relationship and classification of the precursor sequences from the plant miR482 family.
Figure 6. Phylogenetic relationship and classification of the precursor sequences from the plant miR482 family.
Genes 15 01043 g006
Figure 7. Expression patterns of miR482 in C. sinensis (a), N. tabacum (b), S. lycopersicum (c), S. tuberosum (d), F. vesca (e), G. max (f), M. esculenta (g), and P. persica (h). miR482a, miR482b, miR482c, miR482d, miR482e, and miR482f are the members of the miR482 family.
Figure 7. Expression patterns of miR482 in C. sinensis (a), N. tabacum (b), S. lycopersicum (c), S. tuberosum (d), F. vesca (e), G. max (f), M. esculenta (g), and P. persica (h). miR482a, miR482b, miR482c, miR482d, miR482e, and miR482f are the members of the miR482 family.
Genes 15 01043 g007
Figure 8. GO enrichment of predicted miR482-targeted genes.
Figure 8. GO enrichment of predicted miR482-targeted genes.
Genes 15 01043 g008
Figure 9. The ceRNAs network of miR482 in plants: (a) tomato in response to P. infestans; (b) tomato in response to R. solanacearum; (c) cotton in response to drought stress; (d) peanut in response to root-knot nematodes.
Figure 9. The ceRNAs network of miR482 in plants: (a) tomato in response to P. infestans; (b) tomato in response to R. solanacearum; (c) cotton in response to drought stress; (d) peanut in response to root-knot nematodes.
Genes 15 01043 g009
Table 1. Statistics of the miR482 family in plants.
Table 1. Statistics of the miR482 family in plants.
FamiliesSpeciesNameNumber of 5′ EndNumber of 3′ EndEvidence
LiliaceaeAsparagus officinalisaof-miR482033 (experimental)
MyrtaceaeAquilegia caeruleaaqc-miR482030
RutaceaeC. sinensiscsi-miR4827714 (experimental)
MyrtaceaeEugenia unifloraeun-miR482336 (experimental)
RosaceaeF. vescafve-miR482044 (experimental)
MalvaceaeGossypium hirsutumghr-miR482022 (experimental)
FabaceaeG. maxgma-miR482459 (experimental)
MalvaceaeGossypium raimondiigra-miR482033 (experimental)
FabaceaeGlycine sojagso-miR482022 (experimental)
EuphorbiaceaeHevea brasiliensishbr-miR482202 (experimental)
RosaceaeMalus domesticamdm-miR482144 (experimental)
EuphorbiaceaeM. esculentames-miR482054 (experimental)
FabaceaeM. truncatulamtr-miR482112 (experimental)
SolanaceaeN. tabacumnta-miR482143 (experimental)
PinaceaeP. abiespab-miR48212425 (experimental)
PinaceaePinus densatapde-miR482044 (experimental)
AraliaceaePanax ginsengpgi-miR482112 (experimental)
RosaceaeP. persicappe-miR4824610 (experimental)
PinaceaePinus taedapta-miR482044 (by similarity)
SalicaceaeP. trichocarpaptc-miR482358 (experimental)
FabaceaeP. vulgarispvu-miR482112 (experimental)
SolanaceaeS. lycopersicumsly-miR482255 (experimental)
SolanaceaeSolanum tuberosumstu-miR482458 (experimental)
FabaceaeV. unguiculatavun-miR482010
VitaceaeV. viniferavvi-miR482011 (experimental)
PoaceaeZ. mayszma-miR482110
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kuang, W.; Qin, D.; Huang, Y.; Liu, Y.; Cao, X.; Xu, M. Analysis of the miR482 Gene Family in Plants. Genes 2024, 15, 1043. https://doi.org/10.3390/genes15081043

AMA Style

Kuang W, Qin D, Huang Y, Liu Y, Cao X, Xu M. Analysis of the miR482 Gene Family in Plants. Genes. 2024; 15(8):1043. https://doi.org/10.3390/genes15081043

Chicago/Turabian Style

Kuang, Wei, Danfeng Qin, Ying Huang, Yihua Liu, Xue Cao, and Meng Xu. 2024. "Analysis of the miR482 Gene Family in Plants" Genes 15, no. 8: 1043. https://doi.org/10.3390/genes15081043

APA Style

Kuang, W., Qin, D., Huang, Y., Liu, Y., Cao, X., & Xu, M. (2024). Analysis of the miR482 Gene Family in Plants. Genes, 15(8), 1043. https://doi.org/10.3390/genes15081043

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