Age-Related Alterations in Macrophage Distribution and Function Are Associated With Delayed Cutaneous Wound Healing.

Ageing-related delays and dysregulated inflammation in wound healing are well-documented in both human and animal models. However, cellular and molecular changes underlying this impairment in healing progression are not fully understood. In this study, we characterised ageing-associated changes to macrophages in wounds of young and aged mice and investigated transcriptomic differences that may impact the progression of wound healing. Full-thickness wounds created on the dorsum of C57BL/6J young and aged mice were excised on Days 3 and 7 post-wounding for analysis by immunohistochemistry, flow cytometry, and RNA sequencing. Our data revealed that macrophages were significantly reduced in aged wounds in comparison to young. Functional transcriptomic analyses showed that macrophages from aged wounds exhibited significantly reduced expression of cell cycle, DNA replication, and repair pathway genes. Furthermore, we uncovered an elevated pro-inflammatory gene expression program in the aged macrophages correlated with poor inflammation resolution and excessive tissue damage observed in aged wounds. Altogether, our work provides insights into how poorly healing aged wounds are phenotypically defined by the presence of macrophages with reduced proliferative capacity and an exacerbated inflammatory response, both of which are pathways that can be targeted to improve healing in the elderly.

miR-9a regulates levels of both rhomboid mRNA and protein in the early Drosophila melanogaster embryo.

MicroRNAs can have subtle and combinatorial effects on the levels of the targets and pathways they act on. Studying the consequences of a single microRNA knockout often proves difficult as many such knockouts exhibit phenotypes only under stress conditions. This has often led to the hypothesis that microRNAs buffer the effects of intrinsic and environmental stochasticity on gene expression. Observing and understanding this buffering effect entails quantitative analysis of microRNA and target expression in single cells. To this end, we have employed single-molecule fluorescence in situ hybridization, immunofluorescence, and high-resolution confocal microscopy to investigate the effects of miR-9a loss on the expression of the serine-protease Rhomboid in Drosophila melanogaster early embryos. Our single-cell quantitative approach shows that spatially, the rhomboid mRNA pattern is identical in WT and miR-9a knockout embryos. However, we find that the number of mRNA molecules per cell is higher when miR-9a is absent, and the level and temporal accumulation of rhomboid protein shows a more dramatic increase in the miR-9a knockout. Specifically, we see accumulation of rhomboid protein in miR-9a mutants by stage 5, much earlier than in WT. The data, therefore, show that miR-9a functions in the regulation of rhomboid mRNA and protein levels. While further work is required to establish whether rhomboid is a direct target of miR-9 in Drosophila, our results further establish the miR-9 family microRNAs as conserved regulators of timing in neurogenic processes. This study shows the power of single-cell quantification as an experimental tool to study phenotypic consequences of microRNA mis-regulation.

The embryonic transcriptome of Parhyale hawaiensis reveals different dynamics of microRNAs and mRNAs during the maternal-zygotic transition.

Parhyale hawaiensis has emerged as the crustacean model of choice due to its tractability, ease of imaging, sequenced genome, and development of CRISPR/Cas9 genome editing tools. However, transcriptomic datasets spanning embryonic development are lacking, and there is almost no annotation of non-protein-coding RNAs, including microRNAs. We have sequenced microRNAs, together with mRNAs and long non-coding RNAs, in Parhyale using paired size-selected RNA-seq libraries at seven time-points covering important transitions in embryonic development. Focussing on microRNAs, we annotate 175 loci in Parhyale, 88 of which have no known homologs. We use these data to annotate the microRNAome of 37 crustacean genomes, and suggest a core crustacean microRNA set of around 61 sequence families. We examine the dynamic expression of microRNAs and mRNAs during the maternal-zygotic transition. Our data suggest that zygotic genome activation occurs in two waves in Parhyale with microRNAs transcribed almost exclusively in the second wave. Contrary to findings in other arthropods, we do not predict a general role for microRNAs in clearing maternal transcripts. These data significantly expand the available transcriptomics resources for Parhyale, and facilitate its use as a model organism for the study of small RNAs in processes ranging from embryonic development to regeneration.

smiFISH and embryo segmentation for single-cell multi-gene RNA quantification in arthropods.

Recently, advances in fluorescent in-situ hybridization techniques and in imaging technology have enabled visualization and counting of individual RNA molecules in single cells. This has greatly enhanced the resolution in our understanding of transcriptional processes. Here, we adapt a recently published smiFISH protocol (single-molecule inexpensive fluorescent in-situ hybridization) to whole embryos across a range of arthropod model species, and also to non-embryonic tissues. Using multiple fluorophores with distinct spectra and white light laser confocal imaging, we simultaneously detect and separate single RNAs from up to eight different genes in a whole embryo. We also combine smiFISH with cell membrane immunofluorescence, and present an imaging and analysis pipeline for 3D cell segmentation and single-cell RNA counting in whole blastoderm embryos. Finally, using whole embryo single-cell RNA count data, we propose two alternative single-cell variability measures to the commonly used Fano factor, and compare the capacity of these three measures to address different aspects of single-cell expression variability.

Single-cell visualization of mir-9a and Senseless co-expression during Drosophila melanogaster embryonic and larval peripheral nervous system development.

The Drosophila melanogaster peripheral nervous system (PNS) comprises the sensory organs that allow the fly to detect environmental factors such as temperature and pressure. PNS development is a highly specified process where each sensilla originates from a single sensory organ precursor (SOP) cell. One of the major genetic orchestrators of PNS development is Senseless, which encodes a zinc finger transcription factor (Sens). Sens is both necessary and sufficient for SOP differentiation. Senseless expression and SOP number are regulated by the microRNA miR-9a. However, the reciprocal dynamics of Senseless and miR-9a are still obscure. By coupling single-molecule FISH with immunofluorescence, we are able to visualize transcription of the mir-9a locus and expression of Sens simultaneously. During embryogenesis, we show that the expression of mir-9a in SOP cells is rapidly lost as Senseless expression increases. However, this mutually exclusive expression pattern is not observed in the third instar imaginal wing disk, where some Senseless-expressing cells show active sites of mir-9a transcription. These data challenge and extend previous models of Senseless regulation and show complex co-expression dynamics between mir-9a and Senseless. The differences in this dynamic relationship between embryonic and larval PNS development suggest a possible switch in miR-9a function. Our work brings single-cell resolution to the understanding of dynamic regulation of PNS development by Senseless and miR-9a.

The Transcription Factor-microRNA Regulatory Network during hESC-chondrogenesis.

Human embryonic stem cells (ESCs) offer a promising therapeutic approach for osteoarthritis (OA). The unlimited source of cells capable of differentiating to chondrocytes has potential for repairing damaged cartilage or to generate disease models via gene editing. However their use is limited by the efficiency of chondrogenic differentiation. An improved understanding of the transcriptional and post-transcriptional regulation of chondrogenesis will enable us to improve hESC chondrogenic differentiation protocols. Small RNA-seq and whole transcriptome sequencing was performed on distinct stages of hESC-directed chondrogenesis. This revealed significant changes in the expression of several microRNAs including upregulation of known cartilage associated microRNAs and those transcribed from the Hox complexes, and the downregulation of pluripotency associated microRNAs. Integration of miRomes and transcriptomes generated during hESC-directed chondrogenesis identified key functionally related clusters of co-expressed microRNAs and protein coding genes, associated with pluripotency, primitive streak, limb development and extracellular matrix. Analysis identified regulators of hESC-directed chondrogenesis such as miR-29c-3p with 10 of its established targets identified as co-regulated 'ECM organisation' genes and miR-22-3p which is highly co-expressed with ECM genes and may regulate these genes indirectly by targeting the chondrogenic regulators SP1 and HDAC4. We identified several upregulated transcription factors including HOXA9/A10/D13 involved in limb patterning and RELA, JUN and NFAT5, which have targets enriched with ECM associated genes. We have developed an unbiased approach for integrating transcriptome and miRome using protein-protein interactions, transcription factor regulation and miRNA target interactions and identified key regulatory networks prominent in hESC chondrogenesis.

Enhanced genome assembly and a new official gene set for Tribolium castaneum.

BACKGROUND: The red flour beetle Tribolium castaneum has emerged as an important model organism for the study of gene function in development and physiology, for ecological and evolutionary genomics, for pest control and a plethora of other topics. RNA interference (RNAi), transgenesis and genome editing are well established and the resources for genome-wide RNAi screening have become available in this model. All these techniques depend on a high quality genome assembly and precise gene models. However, the first version of the genome assembly was generated by Sanger sequencing, and with a small set of RNA sequence data limiting annotation quality. RESULTS: Here, we present an improved genome assembly (Tcas5.2) and an enhanced genome annotation resulting in a new official gene set (OGS3) for Tribolium castaneum, which significantly increase the quality of the genomic resources. By adding large-distance jumping library DNA sequencing to join scaffolds and fill small gaps, the gaps in the genome assembly were reduced and the N50 increased to 4753kbp. The precision of the gene models was enhanced by the use of a large body of RNA-Seq reads of different life history stages and tissue types, leading to the discovery of 1452 novel gene sequences. We also added new features such as alternative splicing, well defined UTRs and microRNA target predictions. For quality control, 399 gene models were evaluated by manual inspection. The current gene set was submitted to Genbank and accepted as a RefSeq genome by NCBI. CONCLUSIONS: The new genome assembly (Tcas5.2) and the official gene set (OGS3) provide enhanced genomic resources for genetic work in Tribolium castaneum. The much improved information on transcription start sites supports transgenic and gene editing approaches. Further, novel types of information such as splice variants and microRNA target genes open additional possibilities for analysis.

Visualizing gene expression during zebrafish pronephros development and regeneration.

The vertebrate kidney is comprised of functional units known as nephrons. Defects in nephron development or activity are a common feature of kidney disease. Current medical treatments are unable to ameliorate the dire consequences of nephron deficit or injury. Although there have been tremendous advancements in our understanding of nephron ontogeny and the response to damage, many significant knowledge gaps still remain. The zebrafish embryo kidney, or pronephros, is an ideal model for many renal development and regeneration studies because it is comprised of nephrons that share conserved features with the nephron units that comprise the mammalian metanephric kidney. In this chapter, we provide an overview about the benefits of using the zebrafish pronephros to study the mechanisms underlying nephrogenesis as well as epithelial repair and regeneration. We subsequently detail methods for the spatiotemporal assessment of gene and protein expression in zebrafish embryos that can be used to extend the understanding of nephron development and disease, and thereby create new opportunities to identify therapeutic strategies for regenerative medicine.

Prostaglandin signaling regulates renal multiciliated cell specification and maturation.

Multiciliated cells (MCCs) are specialized epithelia with apical bundles of motile cilia that direct fluid flow. MCC dysfunction is associated with human diseases of the respiratory, reproductive, and central nervous systems. Further, the appearance of renal MCCs has been cataloged in several kidney conditions, where their function is unknown. Despite their pivotal health importance, many aspects of MCC development remain poorly understood. Here, we utilized a chemical screen to identify molecules that affect MCC ontogeny in the zebrafish embryo kidney, and found prostaglandin signaling is essential both for renal MCC progenitor formation and terminal differentiation. Moreover, we show that prostaglandin activity is required downstream of the transcription factor ets variant 5a (etv5a) during MCC fate choice, where modulating prostaglandin E2 (PGE2) levels rescued MCC number. The discovery that prostaglandin signaling mediates renal MCC development has broad implications for other tissues, and could provide insight into a multitude of pathological states.

Abundant expression of somatic transposon-derived piRNAs throughout Tribolium castaneum embryogenesis.

BACKGROUND: Piwi-interacting RNAs (piRNAs) are a class of short (~26-31-nucleotide) non-protein-coding RNAs expressed in the metazoan germline. The piRNA pathway in arthropods is best understood in the ovary of Drosophila melanogaster, where it acts to silence active transposable elements (TEs). Maternal loading of piRNAs in oocytes is further required for the inheritance of piRNA-mediated transposon defence. However, our understanding of the diversity, evolution and function of the piRNA complement beyond drosophilids is limited. The red flour beetle, Tribolium castaneum, is an emerging model organism separated from Drosophila by ~ 350 million years of evolution that displays a number of features ancestral to arthropods, including short germ embryogenesis. Here, we characterize the maternally deposited and zygotically expressed small RNA and mRNA complements throughout T. castaneum embryogenesis. RESULTS: We find that beetle oocytes and embryos of all stages are abundant in heterogeneous ~ 28-nucleotide RNAs. These small RNAs originate from discrete genomic loci enriched in TE sequences and display the molecular signatures of transposon-derived piRNAs. In addition to the maternally loaded primary piRNAs, Tribolium embryos produce secondary piRNAs by the cleavage of zygotically activated TE transcripts via the ping-pong mechanism. The two Tribolium piRNA pathway effector proteins, Tc-Piwi/Aub and Tc-Ago3, are also expressed throughout the soma of early embryos. CONCLUSIONS: Our results show that the piRNA pathway in Tribolium is not restricted to the germline, but also operates in the embryo and may act to antagonize zygotically activated transposons. Taken together, these data highlight a functional divergence of the piRNA pathway between insects.

The house spider genome reveals an ancient whole-genome duplication during arachnid evolution.

BACKGROUND: The duplication of genes can occur through various mechanisms and is thought to make a major contribution to the evolutionary diversification of organisms. There is increasing evidence for a large-scale duplication of genes in some chelicerate lineages including two rounds of whole genome duplication (WGD) in horseshoe crabs. To investigate this further, we sequenced and analyzed the genome of the common house spider Parasteatoda tepidariorum. RESULTS: We found pervasive duplication of both coding and non-coding genes in this spider, including two clusters of Hox genes. Analysis of synteny conservation across the P. tepidariorum genome suggests that there has been an ancient WGD in spiders. Comparison with the genomes of other chelicerates, including that of the newly sequenced bark scorpion Centruroides sculpturatus, suggests that this event occurred in the common ancestor of spiders and scorpions, and is probably independent of the WGDs in horseshoe crabs. Furthermore, characterization of the sequence and expression of the Hox paralogs in P. tepidariorum suggests that many have been subject to neo-functionalization and/or sub-functionalization since their duplication. CONCLUSIONS: Our results reveal that spiders and scorpions are likely the descendants of a polyploid ancestor that lived more than 450 MYA. Given the extensive morphological diversity and ecological adaptations found among these animals, rivaling those of vertebrates, our study of the ancient WGD event in Arachnopulmonata provides a new comparative platform to explore common and divergent evolutionary outcomes of polyploidization events across eukaryotes.

Pervasive microRNA Duplication in Chelicerates: Insights from the Embryonic microRNA Repertoire of the Spider Parasteatoda tepidariorum.

MicroRNAs are small (∼22 nt) noncoding RNAs that repress translation and therefore regulate the production of proteins from specific target mRNAs. microRNAs have been found to function in diverse aspects of gene regulation within animal development and many other processes. Among invertebrates, both conserved and novel, lineage specific, microRNAs have been extensively studied predominantly in holometabolous insects such as Drosophila melanogaster However little is known about microRNA repertoires in other arthropod lineages such as the chelicerates. To understand the evolution of microRNAs in this poorly sampled subphylum, we characterized the microRNA repertoire expressed during embryogenesis of the common house spider Parasteatoda tepidariorum We identified a total of 148 microRNAs in P. tepidariorum representing 66 families. Approximately half of these microRNA families are conserved in other metazoans, while the remainder are specific to this spider. Of the 35 conserved microRNAs families 15 had at least two copies in the P. tepidariorum genome. A BLAST-based approach revealed a similar pattern of duplication in other spiders and a scorpion, but not among other chelicerates and arthropods, with the exception of a horseshoe crab. Among the duplicated microRNAs we found examples of lineage-specific tandem duplications, and the duplication of entire microRNA clusters in three spiders, a scorpion, and in a horseshoe crab. Furthermore, we found that paralogs of many P. tepidariorum microRNA families exhibit arm switching, which suggests that duplication was often followed by sub- or neofunctionalization. Our work shows that understanding the evolution of microRNAs in the chelicerates has great potential to provide insights into the process of microRNA duplication and divergence and the evolution of animal development.

Functional and Genetic Analysis of Spectraplakins in Drosophila.

The cytoskeleton is a dynamic network of filamentous protein polymers required for virtually all cellular processes. It consists of three major classes, filamentous actin (F-actin), intermediate filaments, and microtubules, all displaying characteristic structural properties, functions, cellular distributions, and sets of interacting regulatory proteins. One unique class of proteins, the spectraplakins, bind, regulate, and integrate the functions of all three classes of cytoskeleton proteins. Spectraplakins are giant, evolutionary conserved multidomain proteins (spanning up to 9000 aa) that are true members of the plakin, spectrin, and Gas2-like protein families. They have OMIM-listed disease links to epidermolysis bullosa and hereditary sensory and autonomic neuropathy. Their role in disease is likely underrepresented since studies in model animal systems have revealed critical roles in polarity, morphogenesis, differentiation and maintenance, migration, signaling, and intracellular trafficking in a variety of tissues. This enormous diversity of spectraplakin function is consistent with the numerous isoforms produced from single genomic loci that combine different sets of functional domains in distinct cellular contexts. To study the broad range of functions and complexity of these proteins, Drosophila is a powerful model. Thus, the fly spectraplakin Short stop (Shot) acts as an actin-microtubule linker and plays important roles in many developmental processes, which provide experimentally amenable and relevant contexts in which to study spectraplakin functions. For these studies, a versatile range of relevant experimental resources that facilitate genetics and transgenic approaches, highly refined genomics tools, and an impressive set of spectraplakin-specific genetic and molecular tools are readily available. Here, we use the example of Shot to illustrate how the various tools and strategies available for Drosophila can be employed to decipher and dissect cellular roles and molecular mechanisms of spectraplakins.

MicroRNA evolution, expression, and function during short germband development in Tribolium castaneum.

MicroRNAs are well-established players in the development of multicellular animals. Most of our understanding of microRNA function in arthropod development comes from studies in Drosophila. Despite their advantages as model systems, the long germband embryogenesis of fruit flies is an evolutionary derived state restricted to several holometabolous insect lineages. MicroRNA evolution and expression across development in animals exhibiting the ancestral and more widespread short germband mode of embryogenesis has not been characterized. We sequenced small RNA libraries of oocytes and successive intervals covering the embryonic development of the short germband model organism, Tribolium castaneum. We analyzed the evolution and temporal expression of the microRNA complement and sequenced libraries of total RNA to investigate the relationships with microRNA target expression. We show microRNA maternal loading and sequence-specific 3' end nontemplate oligoadenylation of maternally deposited microRNAs that is conserved between Tribolium and Drosophila. We further uncover large clusters encoding multiple paralogs from several Tribolium-specific microRNA families expressed during a narrow interval of time immediately after the activation of zygotic transcription. These novel microRNAs, together with several early expressed conserved microRNAs, target a significant number of maternally deposited transcripts. Comparison with Drosophila shows that microRNA-mediated maternal transcript targeting is a conserved process in insects, but the number and sequences of microRNAs involved have diverged. The expression of fast-evolving and species-specific microRNAs in the early blastoderm of T. castaneum is consistent with previous findings in Drosophila and shows that the unique permissiveness for microRNA innovation at this stage is a conserved phenomenon.

Conserved temporal patterns of microRNA expression in Drosophila support a developmental hourglass model.

The spatiotemporal control of gene expression is crucial for the successful completion of animal development. The evolutionary constraints on development are particularly strong for the mid-embryonic stage when body segments are specified, as evidenced by a high degree of morphological and protein-coding gene conservation during this period--a phenomenon known as the developmental hourglass. The discovery of microRNA-mediated gene control revealed an entirely new layer of complexity of the molecular networks that orchestrate development. However, the constraints on microRNA developmental expression and evolution, and the implications for animal evolution are less well understood. To systematically explore the conservation of microRNAs during development, we carried out a genome-wide comparative study of microRNA expression levels throughout the ontogenesis of two divergent fruit flies, Drosophila melanogaster and D. virilis. We show that orthologous microRNAs display highly similar temporal profiles regardless of their mutation rates, suggesting that the timely expression of microRNA genes can be more constrained than their sequence. Furthermore, transitions between key developmental events in the different species are accompanied by conserved shifts in microRNA expression profiles, with the mid-embryonic period between gastrulation and segmentation characterized by the highest similarity of microRNA expression. The conservation of microRNA expression therefore displays an hourglass pattern similar to that observed for protein-coding genes.

Target repression induced by endogenous microRNAs: large differences, small effects.

MicroRNAs are small RNAs that regulate protein levels. It is commonly assumed that the expression level of a microRNA is directly correlated with its repressive activity - that is, highly expressed microRNAs will repress their target mRNAs more. Here we investigate the quantitative relationship between endogenous microRNA expression and repression for 32 mature microRNAs in Drosophila melanogaster S2 cells. In general, we find that more abundant microRNAs repress their targets to a greater degree. However, the relationship between expression and repression is nonlinear, such that a 10-fold greater microRNA concentration produces only a 10% increase in target repression. The expression/repression relationship is the same for both dominant guide microRNAs and minor mature products (so-called passenger strands/microRNA* sequences). However, we find examples of microRNAs whose cellular concentrations differ by several orders of magnitude, yet induce similar repression of target mRNAs. Likewise, microRNAs with similar expression can have very different repressive abilities. We show that the association of microRNAs with Argonaute proteins does not explain this variation in repression. The observed relationship is consistent with the limiting step in target repression being the association of the microRNA/RISC complex with the target site. These findings argue that modest changes in cellular microRNA concentration will have minor effects on repression of targets.

Fast-evolving microRNAs are highly expressed in the early embryo of Drosophila virilis.

MicroRNAs are short non-protein-coding RNAs that regulate gene expression at the post-transcriptional level and are essential for the embryonic development of multicellular animals. Comparative genome-scale analyses have revealed that metazoan evolution is accompanied by the continuous acquisition of novel microRNA genes. This suggests that novel microRNAs may promote innovation and diversity in development. We determined the evolutionary origins of extant Drosophila microRNAs and estimated the sequence divergence between the 130 orthologous microRNAs in Drosophila melanogaster and Drosophila virilis, separated by 63 million years of evolution. We then generated small RNA sequencing data sets covering D. virilis development and explored the relationship between microRNA conservation and expression in a developmental context. We find that late embryonic, larval, and adult stages are dominated by conserved microRNAs. This pattern, however, does not hold for the early embryo, where rapidly evolving microRNAs are uniquely present at high levels in both species. The group of fast-evolving microRNAs that are highly expressed in the early embryo belong to two Drosophilid lineage-specific clusters: mir-310 ∼ 313 and mir-309 ∼ 6. These clusters have particularly complex evolutionary histories of duplication, gain, and loss. Our analyses suggest that the early embryo is a more permissive environment for microRNA changes and innovations. Fast-evolving microRNAs, therefore, have the opportunity to become preferentially integrated in early developmental processes, and may impact the evolution of development. The relationship between microRNA conservation and expression throughout the development of Drosophila differs from that previously observed for protein-coding genes.

Sex-biased expression of microRNAs in Schistosoma mansoni.

Schistosomiasis is an important neglected tropical disease caused by digenean helminth parasites of the genus Schistosoma. Schistosomes are unusual in that they are dioecious and the adult worms live in the blood system. MicroRNAs play crucial roles during gene regulation and are likely to be important in sex differentiation in dioecious species. Here we characterize 112 microRNAs from adult Schistosoma mansoni individuals, including 84 novel microRNA families, and investigate the expression pattern in different sexes. By deep sequencing, we measured the relative expression levels of conserved and newly identified microRNAs between male and female samples. We observed that 13 microRNAs exhibited sex-biased expression, 10 of which are more abundant in females than in males. Sex chromosomes showed a paucity of female-biased genes, as predicted by theoretical evolutionary models. We propose that the recent emergence of separate sexes in Schistosoma had an effect on the chromosomal distribution and evolution of microRNAs, and that microRNAs are likely to participate in the sex differentiation/maintenance process.

Clusters of microRNAs emerge by new hairpins in existing transcripts.

Genetic linkage may result in the expression of multiple products from a polycistronic transcript, under the control of a single promoter. In animals, protein-coding polycistronic transcripts are rare. However, microRNAs are frequently clustered in the genomes of animals, and these clusters are often transcribed as a single unit. The evolution of microRNA clusters has been the subject of much speculation, and a selective advantage of clusters of functionally related microRNAs is often proposed. However, the origin of microRNA clusters has not been so far explored. Here, we study the evolution of microRNA clusters in Drosophila melanogaster. We observed that the majority of microRNA clusters arose by the de novo formation of new microRNA-like hairpins in existing microRNA transcripts. Some clusters also emerged by tandem duplication of a single microRNA. Comparative genomics show that these clusters are unlikely to split or undergo rearrangements. We did not find any instances of clusters appearing by rearrangement of pre-existing microRNA genes. We propose a model for microRNA cluster evolution in which selection over one of the microRNAs in the cluster interferes with the evolution of the other linked microRNAs. Our analysis suggests that the study of microRNAs and small RNAs must consider linkage associations.