IntroductionMicroRNAs (miRNAs) are an extensive class of noncoding small endogenous RNAs of 18 to 26 nt in length that are derived from self-complementary fold-back structures of longer precursor sequences (pre-miRNAs) and are generated by Dicer-like 1 (DCL1) in plants (Bartel, 2004). Mature miRNAs inhibit gene expression at posttranscriptional levels by either targeting mRNAs for degradation or inhibiting protein translation. Both processes are accomplished by the complementary base pairing of miRNAs to their target mRNA sequences (Ambros, 2004). In plants, for a majority of cases, miRNAs interact with their targets through perfect or near-perfect base pairing and lead to target mRNA degradation (Jones-Rhoades et al., 2006). Increasing evidence has revealed that miRNAs play an important role in a wide range of development processes in plants, including cell proliferation, stress response, metabolism, inflammation, and signal transduction (Ambros, 2004;Jones-Rhoades et al., 2006;Zhang et al., 2007a;Ali et al., 2016) and crosskingdom gene regulation (Barozai and Din, 2017). To date, more than 28,645 miRNAs have been identified from 223 species of plants and animals and deposited in the publicly available database miRBase (Release 21) (Griffiths-Jones et al., 2008). The majority of plant miRNAs have been found in species with fully sequenced genomes including 713 from Oryza sativa, 401 from Populus trichocarpa, 384 from Arabidopsis thaliana, 343 from Solanum tuberosum, 321 from Zea mays, and 241 from Sorghum bicolor (Griffiths-Jones et al., 2008). miRNA-related research is continuously growing and miRNAs, along with their functions, are being identified and elucidated using a wide variety of computational tools and experimental methods including direct cloning, deep sequencing, and other approaches. Comparison of miRNAs across multiple plant species has demonstrated that some miRNAs are highly evolutionary conserved from species to species, such as from mosses to higher flowering eudicots in the plant kingdom (Floyd and Bowman, 2004;Zhang et al., 2006aZhang et al., , 2006b. Conservation of miRNA sequences has provided a powerful strategy for identifying miRNAs in other species. Currently, comparative genome-based homolog searches have been used to identify conserved miRNAs in many plant species, including cotton (