The germ cell-specific TAP-like protein NXF-2 forms a novel granular structure and is required for tra-2 3'UTR-dependent mRNA export in Caenorhabditis elegans

TAP is a general mRNA export receptor and is highly conserved among eukaryotes. The nematode Caenorhabditis elegans has another TAP-like protein, NXF-2, but little is known about its function. In this study, we show that NXF-2 is specifically expressed in germ cells and forms a novel granular structure that is different from that of P granules and that NXF-2 granules are anchored to the nuclear periphery in the mitotic region of the hermaphrodite gonad. In contrast, NXF-2 granules are released within the whole cytoplasm in the meiotic region, where the feminization gene tra-2 starts to function. Both inhibition of XPO-1 (an ortholog of the export receptor CRM1) and mutation of the nuclear export signal of NXF-2 caused the release of NXF-2 granules from the nuclear periphery, indicating that anchoring of NXF-2 granules depends on XPO-1 function. Moreover, inhibition of NXF-2 resulted in a substantial nuclear accumulation of the reporter mRNA carrying the tra-2 3'UTR. These results suggest that, together with XPO-1, NXF-2 exports and anchors tra-2 mRNA to the nuclear periphery to avoid precocious translation until the germ cells reach the meiotic region, thereby contributing to the regulation of tra-2 mRNA expression. This article is protected by copyright. All rights reserved.

Bicc1 and Dicer regulate left-right patterning through post-transcriptional control of the Nodal inhibitor Dand5

Rotating cilia at the vertebrate left-right organizer (LRO) generate an asymmetric leftward flow, which is sensed by cells at the left LRO margin. Ciliary activity of the calcium channel Pkd2 is crucial for flow sensing. How this flow signal is further processed and relayed to the laterality-determining Nodal cascade in the left lateral plate mesoderm (LPM) is largely unknown. We previously showed that flow down-regulates mRNA expression of the Nodal inhibitor Dand5 in left sensory cells. De-repression of the co-expressed Nodal, complexed with the TGFß growth factor Gdf3, drives LPM Nodal cascade induction. Here, we show that post-transcriptional repression of dand5 is a central process in symmetry breaking of Xenopus, zebrafish and mouse. The RNA binding protein Bicc1 was identified as a post-transcriptional regulator of dand5 and gdf3 via their 3'-UTRs. Two distinct Bicc1 functions on dand5 mRNA were observed at pre- and post-flow stages, affecting mRNA stability or flow induced translational inhibition, respectively. To repress dand5, Bicc1 co-operates with Dicer1, placing both proteins in the process of flow sensing. Intriguingly, Bicc1 mediated translational repression of a dand5 3'-UTR mRNA reporter was responsive to pkd2, suggesting that a flow induced Pkd2 signal triggers Bicc1 mediated dand5 inhibition during symmetry breakage.

Evolution after Whole-Genome Duplication: Teleost MicroRNAs

MicroRNAs (miRNAs) are important gene expression regulators implicated in many biological processes, but we lack a global understanding of how miRNA genes evolve and contribute to developmental canalization and phenotypic diversification. Whole-genome duplication events likely provide a substrate for species divergence and phenotypic change by increasing gene numbers and relaxing evolutionary pressures. To understand the consequences of genome duplication on miRNA evolution, we studied miRNA genes following the teleost genome duplication (TGD). Analysis of miRNA genes in four teleosts and in spotted gar, whose lineage diverged before the TGD, revealed that miRNA genes were retained in ohnologous pairs more frequently than protein-coding genes, and that gene losses occurred rapidly after the TGD. Genomic context influenced retention rates, with clustered miRNA genes retained more often than nonclustered miRNA genes and intergenic miRNA genes retained more frequently than intragenic miRNA genes, which often shared the evolutionary fate of their protein-coding host. Expression analyses revealed both conserved and divergent expression patterns across species in line with miRNA functions in phenotypic canalization and diversification, respectively. Finally, major strands of miRNA genes experienced stronger purifying selection, especially in their seeds and 3'-complementary regions, compared with minor strands, which nonetheless also displayed evolutionary features compatible with constrained function. This study provides the first genome-wide, multispecies analysis of the mechanisms influencing metazoan miRNA evolution after whole-genome duplication.

Evolutionary conservation in noncoding genomic regions

Humans may share more genomic commonalities with other species than previously thought. According to current estimates, ~5% of the human genome is functionally constrained, which is a much larger fraction than the ~1.5% occupied by annotated protein-coding genes. Hence, ~3.5% of the human genome comprises likely functional conserved noncoding elements (CNEs) preserved among organisms, whose common ancestors existed throughout hundreds of millions of years of evolution. As whole-genome sequencing emerges as a standard procedure in genetic analyses, interpretation of variations in CNEs, including the elucidation of mechanistic and functional roles, becomes a necessity. Here, we discuss the phenomenon of noncoding conservation via four dimensions (sequence, regulatory conservation, spatiotemporal expression, and structure) and the potential significance of CNEs in phenotype variation and disease.

Genome-Wide Identification and Transcriptomic Analysis of MicroRNAs Across Various Amphioxus Organs Using Deep Sequencing

Amphioxus is the closest living invertebrate proxy of the vertebrate ancestor. Systematic gene identification and expression profile analysis of amphioxus organs are thus important for clarifying the molecular mechanisms of organ function formation and further understanding the evolutionary origin of organs and genes in vertebrates. The precise regulation of microRNAs (miRNAs) is crucial for the functional specification and differentiation of organs. In particular, those miRNAs that are expressed specifically in organs (OSMs) play key roles in organ identity, differentiation, and function. In this study, the genome-wide miRNA transcriptome was analyzed in eight organs of adult amphioxus Branchiostoma belcheri using deep sequencing. A total of 167 known miRNAs and 23 novel miRNAs (named novel_mir), including 139 conserved miRNAs, were discovered, and 79 of these were identified as OSMs. Additionally, analyses of the expression patterns of eight randomly selected known miRNAs demonstrated the accuracy of the miRNA deep sequencing that was used in this study. Furthermore, potentially OSM-regulated genes were predicted for each organ type. Functional enrichment of these predicted targets, as well as further functional analyses of known OSMs, was conducted. We found that the OSMs were potentially to be involved in organ-specific functions, such as epidermis development, gonad development, muscle cell development, proteolysis, lipid metabolism, and generation of neurons. Moreover, OSMs with non-organ-specific functions were detected and primarily include those related to innate immunity and response to stimuli. These findings provide insights into the regulatory roles of OSMs in various amphioxus organs.

The expression and regulation of miR-1 in goat skeletal muscle and satellite cell during muscle growth and development

MicroRNA-1 (miR-1) has been shown to play an important role in muscle growth and development, however, it was mainly discovered in model animals. To explore the function and mechanism of miR-1 in goat, we firstly explored the expression profile of miR-1 in goat tissues and cells. Furthermore, the target gene of miR-1 was predicted, and the relationship between miR-1 and one of its target genes, histone deacetylase 4 (HDAC4), was analyzed through double luciferase reporter assay, real-time PCR, and western blot. It was found that the miR-1 is most abundantly expressed in goat heart and skeletal muscle tissue. Meanwhile, the expression of miR-1 showed an increasing tendency from new-born goats to the 7-month-old goats, and then its expression decreases as the goats mature further. In addition, the expression levels of miR-1 decreased in goat skeletal muscle satellite cells with the algebraic increasing of cells. At last, the results showed that HDAC4 is a target gene of miR-1 in goat, and miR-1 can inhibit the post-transcriptional expression of HDAC4, but had no significant influence on the mRNA level of HDAC4. It was hypothesized that miR-1 promotes muscle development by inhibiting the post-transcriptional expression of HDAC4 in goat.

The Ciona myogenic regulatory factor functions as a typical MRF but possesses a novel N-terminus that is essential for activity

Electroporation-based assays were used to test whether the myogenic regulatory factor (MRF) of Ciona intestinalis (CiMRF) interferes with endogenous developmental programs, and to evaluate the importance of its unusual N-terminus for muscle development. We found that CiMRF suppresses both notochord and endoderm development when it is expressed in these tissues by a mechanism that may involve activation of muscle-specific microRNAs. Because these results add to a large body of evidence demonstrating the exceptionally high degree of functional conservation among MRFs, we were surprised to discover that non-ascidian MRFs were not myogenic in Ciona unless they formed part of a chimeric protein containing the CiMRF N-terminus. Equally surprising, we found that despite their widely differing primary sequences, the N-termini of MRFs of other ascidian species could form chimeric MRFs that were also myogenic in Ciona. This domain did not rescue the activity of a Brachyury protein whose transcriptional activation domain had been deleted, and so does not appear to constitute such a domain. Our results indicate that ascidians have previously unrecognized and potentially novel requirements for MRF-directed myogenesis. Moreover, they provide the first example of a domain that is essential to the core function of an important family of gene regulatory proteins, one that, to date, has been found in only a single branch of the family.

Developmental mechanisms of migratory muscle precursors in medaka pectoral fin formation

Limb muscles are formed from migratory muscle precursor cells (MMPs) that delaminate from the ventral region of dermomyotomes and migrate into the limb bud. MMPs remain undifferentiated during migration, commencing differentiation into skeletal muscle after arrival in the limb. However, it is still unclear whether the developmental mechanisms of MMPs are conserved in teleost fishes. Here, we investigate the development of pectoral fin muscles in the teleost medaka Oryzias latipes. Expression of the MMP marker lbx1 is first observed in several somites prior to the appearance of fin buds. lbx1-positive cells subsequently move anteriorly and localize in the prospective fin bud region to differentiate into skeletal muscle cells. To address the developmental mechanisms underlying fin muscle formation, we knocked down tbx5, a gene that is required for fin bud formation. tbx5 morphants showed loss of fin buds, whereas lbx1 expression initiated normally in anterior somites. Unlike in normal embryos, expression of lbx1 was not maintained in migrating fin MMPs or within the fin buds. We suggest that fin MMPs appear to undergo two phases in their development, with an initial specification of MMPs occurring independent of fin buds and a second fin bud-dependent phase of MMP migration and proliferation. Our results showed that medaka fin muscle is composed of MMPs. It is suggested that the developmental mechanism of fin muscle formation is conserved in teleost fishes including medaka. Through this study, we also propose new insights into the developmental mechanisms of MMPs in fin bud formation.

Exosomes as Mediators of the Systemic Adaptations to Endurance Exercise

Habitual endurance exercise training is associated with multisystemic metabolic adaptations that lower the risk of inactivity-associated disorders such as obesity and type 2 diabetes mellitus (T2DM). Identification of complex systemic signaling networks responsible for these benefits are of great interest because of their therapeutic potential in metabolic diseases; however, specific signals that modulate the multisystemic benefits of exercise in multiple tissues and organs are only recently being discovered. Accumulated evidence suggests that muscle and other tissues have an endocrine function and release peptides and nucleic acids into the circulation in response to acute endurance exercise to mediate the multisystemic adaptations. Factors released from skeletal muscle have been termed myokines and we propose that the total of all factors released in response to endurance exercise (including peptides, nucleic acids, and metabolites) be termed, "exerkines." We propose that many of the exerkines are released within extracellular vesicles called exosomes, which regulate peripheral organ cross talk. Exosomes (30-140 nm) and larger microvesicles [MVs] (100-1000 nm) are subcategories of extracellular vesicles that are released into the circulation. Exosomes contain peptides and several nucleic acids (microRNA [miRNA], messenger RNA [mRNA], mitochondrial DNA [mtDNA]) and are involved in intercellular/tissue exchange of their contents. An acute bout of endurance exercise increases circulating exosomes that are hypothesized to mediate organ cross talk to promote systemic adaptation to endurance exercise. Further support for the role of exosomes (and possibly MVs) in mediating the systemic benefits of exercise comes from the fact that the majority of the previously reported myokines/exerkines are found in extracellular vesicles databases (Vesiclepedia and ExoCarta). We propose that exosomes isolated from athletes following exercise or exosomes bioengineered to incorporate one or many of known exerkines will be therapeutically useful in the treatment of obesity, T2DM, and other aging-associated metabolic disorders.

The potential of endurance exercise-derived exosomes to treat metabolic diseases

Endurance exercise-mediated multisystemic adaptations are known to mitigate metabolism-related disorders such as obesity and type 2 diabetes mellitus (T2DM). However, the underlying molecular mechanisms that promote crosstalk between organs and orchestrate the pro-metabolic effects of endurance exercise remain unclear. Exercise-induced release of peptides and nucleic acids from skeletal muscle and other organs (collectively termed 'exerkines') has been implicated in mediating these systemic adaptations. Given that the extracellular milieu is probably not a hospitable environment for labile exerkines, a lipid vehicle-based mode of delivery has originated over the course of evolution. Two types of extracellular vesicles, exosomes and microvesicles, have been shown to contain proteins and nucleic acids that participate in a variety of physiological and pathological processes. Exosomes, in particular, have been shown to facilitate the exchange of peptides, microRNA, mRNA and mitochondrial DNA between cells and tissues. Intriguingly, circulatory extracellular vesicle content increases in an intensity-dependant manner in response to endurance exercise. We propose that the systemic benefits of exercise are modulated by exosomes and/or microvesicles functioning in an autocrine, paracrine and/or endocrine manner. Furthermore, we posit that native or modified exosomes, and/or microvesicles enriched with exerkines will have therapeutic utility in the treatment of obesity and T2DM.

Automated detection of ncRNAs in the draft genome sequence of a colonial tunicate: the carpet sea squirt Didemnum vexillum

BACKGROUND

The colonial ascidian Didemnum vexillum, sea carpet squirt, is not only a key marine organism to study morphological ancestral patterns of chordates evolution but it is also of great ecological importance due to its status as a major invasive species. Non-coding RNAs, in particular microRNAs (miRNAs), are important regulatory genes that impact development and environmental adaptation. Beyond miRNAs, not much in known about tunicate ncRNAs.

RESULTS

We provide here a comprehensive homology-based annotation of non-coding RNAs in the recently sequenced genome of D. vexillum. To this end we employed a combination of several computational approaches, including blast searches with a wide range of parameters, and secondary structured centered survey with infernal. The resulting candidate set was curated extensively to produce a high-quality ncRNA annotation of the first draft of the D. vexillum genome. It comprises 57 miRNA families, 4 families of ribosomal RNAs, 22 isoacceptor classes of tRNAs (of which more than 72 % of loci are pseudogenes), 13 snRNAs, 12 snoRNAs, and 1 other RNA family. Additionally, 21 families of mitochondrial tRNAs and 2 of mitochondrial ribosomal RNAs and 1 long non-coding RNA.

CONCLUSIONS

The comprehensive annotation of the D. vexillum non-coding RNAs provides a starting point towards a better understanding of the restructuring of the small RNA system in ascidians. Furthermore it provides a valuable research for efforts to establish detailed non-coding RNA annotations for other recently published and recently sequences in tunicate genomes.

Cardiac-specific miRNA in cardiogenesis, heart function, and cardiac pathology (with focus on myocardial infarction)

Cardiac miRNAs (miR-1, miR133a, miR-208a/b, and miR-499) are abundantly expressed in the myocardium. They play a central role in cardiogenesis, heart function and pathology. While miR-1 and miR-133a predominantly control early stages of cardiogenesis supporting commitment of cardiac-specific muscle lineage from embryonic stem cells and mesodermal precursors, miR-208 and miR-499 are involved in the late cardiogenic stages mediating differentiation of cardioblasts to cardiomyocytes and fast/slow muscle fiber specification. In the heart, miR-1/133a control cardiac conductance and automaticity by regulating all phases of the cardiac action potential. miR-208/499 located in introns of the heavy chain myosin genes regulate expression of sarcomeric contractile proteins. In cardiac pathology including myocardial infarction (MI), expression of cardiac miRNAs is markedly altered that leads to deleterious effects associated with heart wounding, arrhythmia, increased apoptosis, fibrosis, hypertrophy, and tissue remodeling. In acute MI, circulating levels of cardiac miRNAs are significantly elevated making them to be a promising diagnostic marker for early diagnosis of acute MI. Great cardiospecific capacity of these miRNAs is very helpful for enhancing regenerative properties and survival of stem cell and cardiac progenitor transplants and for reprogramming of mature non-cardiac cells to cardiomyocytes.

Restricted distribution of mrg-1 mRNA in C. elegans primordial germ cells through germ granule-independent regulation

The chromodomain protein MRG-1 is an essential maternal factor for proper germline development that protects germ cells from cell death in C. elegans. Unlike germ granules, which are exclusively segregated to the germline blastomeres at each cell division from the first cleavage of the embryo, MRG-1 is abundant in all cells in early embryos and is then gradually restricted to the primordial germ cells (PGCs) by the morphogenesis stage. Here, we show that this characteristic spatiotemporal expression pattern is dictated by the mrg-1 3'UTR and is differentially regulated at the RNA level between germline and somatic cells. Asymmetric segregation of germ granules is not necessary to localize MRG-1 to the PGCs. We found that MES-4, an essential chromatin regulator in germ cells, also accumulates in the PGCs in a germ granule-independent manner. We propose that C.elegans PGCs have a novel mechanism to accumulate at least some chromatin-associated proteins that are essential for germline immortality.

Developmental regulation and evolution of muscle-specific microRNAs

MicroRNAs (miRs) are a group of small RNAs that play a major role in post-transcriptional regulation of gene expression. In animals, many of the miRs are expressed in a conserved spatiotemporal manner. Muscle tissues, the major cellular systems involved in the locomotion and physiological functions of animals, have been one of the main sites for verification of miR targets and analysis of their developmental functions. During the determination and differentiation of muscle cells, numerous miRs bind to and repress target mRNAs in a highly specific but redundant manner. Interspecific comparisons of the sequences and expression of miRs have suggested that miR regulation became increasingly important during the course of vertebrate evolution. However, the detailed molecular interactions that have led to the highly complex morphological structures still await investigation. In this review, we will summarize the recent findings on the functional and developmental characteristics of miRs that have played major roles in vertebrate myogenesis, and discuss how the evolution of miRs is related to the morphological complexity of the vertebrates.

MicroRNA-499 expression distinctively correlates to target genes sox6 and rod1 profiles to resolve the skeletal muscle phenotype in Nile tilapia

A class of small non-coding RNAs, the microRNAs (miRNAs), has been shown to be essential for the regulation of specific cell pathways, including skeletal muscle development, maintenance and homeostasis in vertebrates. However, the relative contribution of miRNAs for determining the red and white muscle cell phenotypes is far from being fully comprehended. To better characterize the role of miRNA in skeletal muscle cell biology, we investigated muscle-specific miRNA (myomiR) signatures in Nile tilapia fish. Quantitative (RT-qPCR) and spatial (FISH) expression analyses revealed a highly differential expression (forty-four-fold) of miR-499 in red skeletal muscle compared to white skeletal muscle, whereas the remaining known myomiRs were equally expressed in both muscle cell types. Detailed examination of the miR-499 targets through bioinformatics led us to the sox6 and rod1 genes, which had low expression in red muscle cells according to RT-qPCR, FISH, and protein immunofluorescence profiling experiments. Interestingly, we verified that the high expression of miR-499 perfectly correlates with a low expression of sox6 and rod1 target genes, as verified by a distinctive predominance of mRNA destabilization and protein translational decay to these genes, respectively. Through a genome-wide comparative analysis of SOX6 and ROD1 protein domains and through an in silico gene regulatory network, we also demonstrate that both proteins are essentially similar in vertebrate genomes, suggesting their gene regulatory network may also be widely conserved. Overall, our data shed light on the potential regulation of targets by miR-499 associated with the slow-twitch muscle fiber type phenotype. Additionally the results provide novel insights into the evolutionary dynamics of miRNA and target genes enrolled in a putative constrained molecular pathway in the skeletal muscle cells of vertebrates.

Regulation of Cardiac Cell Fate by microRNAs: Implications for Heart Regeneration

microRNAs are post-transcriptional regulators of gene expression that have been shown to be central players in the establishment of cellular programs, often acting as switches that control the choice between proliferation and differentiation during development and in adult tissues. The heart develops from two small patches of cells in the mesoderm, the heart fields, which originate the different cardiac cell types, including cardiomyocytes, vascular smooth muscle and endothelial cells. These progenitors proliferate and differentiate to establish a highly connected three-dimensional structure, involving a robust succession of gene expression programs strongly influenced by microRNAs. Although the mammalian heart has conventionally been viewed as a post-mitotic organ, cardiac cells have recently been shown to display some regenerative potential, which is nonetheless insufficient to regenerate heart lesions, in contrast with other vertebrates like the zebrafish. Both the proliferation of adult cardiac stem cells and the ability of cardiomyocytes to re-enter the cell cycle have been proposed to sustain these regenerative processes. Here we review the role of microRNAs in the control of stem cell and cardiomyocyte dependent cardiac regeneration processes, and discuss potential applications for the treatment of cardiac injury.

Evolution and genomic organization of muscle microRNAs in fish genomes

Background

MicroRNAs (miRNAs) are small non-coding RNA molecules with an important role upon post-transcriptional regulation. These molecules have been shown essential for several cellular processes in vertebrates, including muscle biology. Many miRNAs were described as exclusively or highly expressed in skeletal and/or cardiac muscle. However, knowledge on the genomic organization and evolution of muscle miRNAs has been unveiled in a reduced number of vertebrates and mostly only reflects their organization in mammals, whereas fish genomes remain largely uncharted. The main goal of this study was to elucidate particular features in the genomic organization and the putative evolutionary history of muscle miRNAs through a genome-wide comparative analysis of cartilaginous and bony fish genomes.

Results

As major outcomes we show that (1) miR-208 was unexpectedly absent in cartilaginous and ray-finned fish genomes whereas it still exist in other vertebrate groups; (2) miR-499 was intergenic in medaka and stickleback conversely to other vertebrates where this miRNA is intronic; (3) the zebrafish genome is the unique harboring two extra paralogous copies of miR-499 and their host gene (Myh7b); (4) a rare deletion event of the intergenic and bicistronic cluster miR-1-1/133a-2 took place only into Tetraodontiformes genomes (pufferfish and spotted green puffer); (5) the zebrafish genome experienced a duplication event of miR-206/-133b; and (6) miR-214 was specifically duplicated in species belonging to superorder Acanthopterygii.

Conclusions

Despite of the aforementioned singularities in fish genomes, large syntenic blocks containing muscle-enriched miRNAs were found to persist, denoting colligated functionality between miRNAs and neighboring genes. Based on the genomic data here obtained, we envisioned a feasible scenario for explaining muscle miRNAs evolution in vertebrates.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0196-x) contains supplementary material, which is available to authorized users.

Functional effects of dopamine transporter gene genotypes on in vivo dopamine transporter functioning: a meta-analysis

Much psychiatric genetic research has focused on a 40-base pair variable number of tandem repeats (VNTR) polymorphism located in the 3'-untranslated region (3'UTR) of the dopamine active transporter (DAT) gene (SLC6A3). This variant produces two common alleles with 9- and 10-repeats (9R and 10R). Studies associating this variant with in vivo DAT activity in humans have had mixed results. We searched for studies using positron emission tomography (PET) or single-photon emission computed tomography (SPECT) to evaluate this association. Random effects meta-analyses assessed the association of the 3'UTR variant with DAT activity. We also evaluated heterogeneity among studies and evidence for publication bias. We found twelve studies comprising 511 subjects, 125 from PET studies and 386 from SPECT studies. The PET studies provided highly significant evidence that the 9R allele was associated with increased DAT activity in human adults. The SPECT studies were highly heterogeneous. As a group, they suggested no association between the 3'UTR polymorphism and DAT activity. When the analysis was limited to the most commonly used ligand, [123I]β-CIT, stratification by affection status dramatically reduced heterogeneity and revealed a significant association of the 9R allele with increased DAT activity for healthy subjects. In humans, the 9R allele of the 3'UTR polymorphism of SLC6A3 regulates dopamine activity in the striatal brain regions independent of the presence of neuropsychiatric illness. Differences in study methodology account for the heterogeneous results across individual studies.

A tri-component conservation strategy reveals highly confident microRNA-mRNA interactions and evolution of microRNA regulatory networks

MicroRNAs are small non-coding RNAs that can regulate expressions of their target genes at the post-transcriptional level. In this study, we propose a tri-component strategy that combines the conservation of microRNAs, homology of mRNA coding regions, and conserved microRNA binding sites in the 3′ untranslated regions to discover conserved microRNA-mRNA interactions. To validate the performance of our conservation strategy, we collected the experimentally validated microRNA-mRNA interactions from three databases as the golden standard. We found that the proposed strategy can improve the performance of existing target prediction algorithms by approximately 2–4 fold. In addition, we demonstrated that the proposed strategy could efficiently retain highly confident interactions from the intersection results of the existing algorithms and filter out the possible false positive predictions in the union one. Furthermore, this strategy can facilitate our ability to trace the homologues in different species that are targeted by the same miRNA family because it combines these three features to identify the conserved miRNA-mRNA interactions during evolution. Through an extensive application of the proposed conservation strategy to a study of the miR-1/206 regulatory network, we demonstrate that the target mRNA recruiting process could be associated with expansion of miRNA family during its evolution. We also uncovered the functional evolution of the miR-1/206 regulatory network. In this network, the early targeted genes tend to participate in more general and development-related functions. In summary, the conservation strategy is capable of helping to highlight the highly confident miRNA-mRNA interactions and can be further applied to reveal the evolutionary features of miRNA regulatory network and functions.