Comparative genomic analysis reveals evolutionary characteristics and patterns of microRNA clusters in vertebrates

MicroRNAs (miRNAs) are a class of small non-coding RNAs that can play important regulatory roles in many important biological processes. Although clustering patterns of miRNA clusters have been uncovered in animals, the origin and evolution of miRNA clusters in vertebrates are still poorly understood. Here, we performed comparative genomic analyses to construct 51 sets of orthologous miRNA clusters (SOMCs) across seven test vertebrate species, a collection of miRNA clusters from two or more species that are likely to have evolved from a common ancestral miRNA cluster, and used these to systematically examine the evolutionary characteristics and patterns of miRNA clusters in vertebrates. We found that miRNA clusters are continuously generated, and most of them tend to be conserved and maintained in vertebrate genomes, although some adaptive gains and losses of miRNA cluster have occurred during evolution. Furthermore, miRNA clusters appeared relatively early in the evolutionary history might suffer from more complicated adaptive gain-and-loss than those young miRNA clusters. Detailed analysis showed that genomic duplication events of ancestral miRNAs or miRNA clusters are likely to be major driving force and apparently contribute to origin and evolution of miRNA clusters. Comparison of conserved with lineage-specific miRNA clusters revealed that the contribution of duplication events for the formation of miRNA cluster appears to be more important for conserved miRNA clusters than lineage-specific. Our study provides novel insights for further exploring the origins and evolution of miRNA clusters in vertebrates at a genome scale.     

Systematic investigation of Amphioxus (Branchiostoma floridae) microRNAs

MicroRNAs (miRNAs) participate in various biological processes via controlling gene activity. Amphioxus is the best available stand-in as the proximate invertebrate ancestor of the vertebrates. Here, we systematically investigated the miRNAs in amphioxus. First, we identified 245 candidate amphioxus miRNAs, in which 183 miRNAs were firstly reported. Second, we gave evidences to support a birth-and-death process of miRNA genes in some families and gave implications for the functional diversification of miRNA during evolution. Third, we identified 47 development-specific expression miRNAs. We found that only 19 miRNAs were expressed in all developmental stages, 16 miRNAs were neurula-specific and 13 miRNAs were larva-specific. In addition, these potential miRNA-targeting genes were mainly classified into development, muscle formation, cell adhesion, and gene regulation categories. Finally, we found 79 immune related genes targeted by 136 miRNAs in amphioxus. In conclusion, our results take an insight into both the function and evolution of the amphioxus miRNAs.     

Characterization of microRNAs in cephalochordates reveals a correlation between microRNA repertoire homology and morphological similarity in chordate evolution

SUMMARY Cephalochordates, urochordates, and vertebrates comprise the three extant groups of chordates. Although higher morphological and developmental similarity exists between cephalochordates and vertebrates, molecular phylogeny studies have instead suggested that the morphologically simplified urochordates are the closest relatives to vertebrates. MicroRNAs (miRNAs) are regarded as the major factors driving the increase of morphological complexity in early vertebrate evolution, and are extensively characterized in vertebrates and in a few species of urochordates. However, the comprehensive set of miRNAs in the basal chordates, namely the cephalochordates, remains undetermined. Through extensive sequencing of a small RNA library and genomic homology searches, we characterized 100 miRNAs from the cephalochordate amphioxus, Branchiostoma japonicum , and B. floridae . Analysis of the evolutionary history of the cephalochordate miRNAs showed that cephalochordates possess 54 miRNA families homologous to those of vertebrates, which is threefold higher than those shared between urochordates and vertebrates. The miRNA contents demonstrated a clear correlation between the extent of miRNA overlapping and morphological similarity among the three chordate groups, providing a strong evidence of miRNAs being the major genetic factors driving morphological complexity in early chordate evolution.     

Systematic analysis of genomic organization and heterogeneities of miRNA cluster in vertebrates

The clustering propensity of microRNA genes is a common biological phenomenon in various animal and plant species. To gain novel insight into genomic organization and potential functional heterogeneities of miRNA clusters in vertebrates from a genome scale, we used large scale data and presented a comprehensive analysis to examine various features of genomic organization of miRNA clusters across seven vertebrates by a combination of comparative genomics and bioinformatics approaches. The results of pair-wise distance analysis of same-strand consecutive miRNAs suggested that the fractions of the miRNA gene pairs are higher at relatively short pair-wise distances than those of protein-coding genes and other non-coding RNA genes. Especially relatively small number of miRNAs is more clustered at very short pair-wise distances than expected at random. We further observed significant difference between real miRNA clusters and randomly organized clusters for different aspects, including higher overlap of target genes, fewer seed types and significant enrichment in diseases. However, the extent of these features of clustered miRNAs has a different tendency and largely depends on inter-miRNA distances because of diverse clustering propensity of miRNAs in vertebrates, suggesting that this cooperated function or cooperative effects between miRNAs in clusters perhaps be affected by inter-miRNA distances.     

Phylogenetic shadowing and computational identification of human microRNA genes

We sequenced 122 miRNAs in 10 primate species to reveal conservation characteristics of miRNA genes. Strong conservation is observed in stems of miRNA hairpins and increased variation in loop sequences. Interestingly, a striking drop in conservation was found for sequences immediately flanking the miRNA hairpins. This characteristic profile was employed to predict novel miRNAs using cross-species comparisons. Nine hundred and seventy-six candidate miRNAs were identified by scanning whole-genome human/mouse and human/rat alignments. Most of the novel candidates are conserved also in other vertebrates (dog, cow, chicken, opossum, zebrafish). Northern blot analysis confirmed the expression of mature miRNAs for 16 out of 69 representative candidates. Additional support for the expression of 179 novel candidates can be found in public databases, their presence in gene clusters, and literature that appeared after these predictions were made. Taken together, these results suggest the presence of significantly higher numbers of miRNAs in the human genome than previously estimated.     

Identification of novel chicken microRNAs and analysis of their genomic organization

MicroRNAs (miRNAs) represent a family of small noncoding RNAs with important regulatory roles in diverse biological processes ranging from cell differentiation to organism development. In chickens, the full set of miRNAs and the expression patterns of miRNAs during development are still poorly understood when compared to the other vertebrates. In this study, we identified 29 novel miRNAs and 140 potential miRNA loci in the chicken genome by combining the experimental and computational analyses. Detailed expression patterns of 49 miRNAs were first characterized by Northern blotting and indicated the cooperativity of the miRNA expression with their function in embryogenesis and organogenesis. Twenty-seven miRNA clusters were systematically evaluated in the chicken genome and diverse expression patterns for closely linked miRNAs were observed. Our results significantly expand the set of known miRNAs in the chicken and provide the basis for understanding the structural and functional evolution of miRNA genes in vertebrates.     

Characterization of novel GPCR gene coding locus in amphioxus genome: gene structure, expression, and phylogenetic analysis with implications for its involvement in chemoreception

Chemosensation is the primary sensory modality in almost all metazoans. The vertebrate olfactory receptor genes exist as tandem clusters in the genome, so that identifying their evolutionary origin would be useful for understanding the expansion of the sensory world in relation to a large-scale genomic duplication event in a lineage leading to the vertebrates. In this study, I characterized a novel GPCR (G-protein-coupled receptor) gene-coding locus from the amphioxus genome. The genomic DNA contains an intronless ORF whose deduced amino acid sequence encodes a seven-transmembrane protein with some amino acid residues characteristic of vertebrate olfactory receptors (ORs). Surveying counterparts in the Ciona intestinalis (Asidiacea, Urochordata) genome by querying BLAST programs against the Ciona genomic DNA sequence database resulted in the identification of a remotely related gene. In situ hybridization analysis labeled primary sensory neurons in the rostral epithelium of amphioxus adults. Based on these findings, together with comparison of the developmental gene expression between amphioxus and vertebrates, I postulate that chemoreceptive primary sensory neurons in the rostrum are an ancient cell population traceable at least as far back in phylogeny as the common ancestor of amphioxus and vertebrates.     

The amphioxus genome enlightens the evolution of the thyroid hormone signaling pathway

Thyroid hormones (THs) have pleiotropic effects on vertebrate development, with amphibian metamorphosis as the most spectacular example. However, developmental functions of THs in non-vertebrate chordates are largely hypothetical and even TH endogenous production has been poorly investigated. In order to get better insight into the evolution of the thyroid hormone signaling pathway in chordates, we have taken advantage of the recent release of the amphioxus genome. We found amphioxus homologous sequences to most of the genes encoding proteins involved in thyroid hormone signaling in vertebrates, except the fast-evolving thyroglobulin: sodium iodide symporter, thyroid peroxidase, deiodinases, thyroid hormone receptor, TBG, and CTHBP. As only some genes encoding proteins involved in TH synthesis regulation were retrieved (TRH, TSH receptor, and CRH receptor but not their corresponding receptors and ligands), there may be another mode of upstream regulation of TH synthesis in amphioxus. In accord with the notion that two whole genome duplications took place at the base of the vertebrate tree, one amphioxus gene often corresponded to several vertebrate homologs. However, some amphioxus specific duplications occurred, suggesting that several steps of the TH pathway were independently elaborated in the cephalochordate and vertebrate lineages. The present results therefore indicate that amphioxus is capable of producing THs. As several genes of the TH signaling pathway were also found in the sea urchin genome, we propose that the thyroid hormone signaling pathway is of ancestral origin in chordates, if not in deuterostomes, with specific elaborations in each lineage, including amphioxus.     

Identification of microRNAs of the herpesvirus family

Epstein-Barr virus (EBV or HHV4), a member of the human herpesvirus (HHV) family, has recently been shown to encode microRNAs (miRNAs). In contrast to most eukaryotic miRNAs, these viral miRNAs do not have close homologs in other viral genomes or in the genome of the human host. To identify other miRNA genes in pathogenic viruses, we combined a new miRNA gene prediction method with small-RNA cloning from several virus-infected cell types. We cloned ten miRNAs in the Kaposi sarcoma-associated virus (KSHV or HHV8), nine miRNAs in the mouse gammaherpesvirus 68 (MHV68) and nine miRNAs in the human cytomegalovirus (HCMV or HHV5). These miRNA genes are expressed individually or in clusters from either polymerase (pol) II or pol III promoters, and share no substantial sequence homology with one another or with the known human miRNAs. Generally, we predicted miRNAs in several large DNA viruses, and we could neither predict nor experimentally identify miRNAs in the genomes of small RNA viruses or retroviruses.     

Altered miRNA repertoire in the simplified chordate, Oikopleura dioica

Recent studies reveal correlation between microRNA (miRNA) innovation and increased developmental complexity. This is exemplified by dramatic expansion of the miRNA inventory in vertebrates, a lineage where genome duplication has played a significant evolutionary role. Urochordates, the closest extant group to the vertebrates, exhibit an opposite trend to genome and morphological simplification. We show that the urochordate, larvacean, Oikopleura dioica, possesses the requisite miRNA biogenic machinery. The miRNAs isolated by small RNA cloning were expressed throughout the short life cycle, a number of which were stocked as maternal determinants prior to rapid embryonic development. We identify sex-specific miRNAs that appeared as male/female gonad differentiation became apparent and were maintained throughout spermatogenesis. Whereas 80% of mammalian miRNAs are hosted in introns of protein-coding genes, the majority of O. dioica miRNA loci were located in antisense orientations to such genes. Including sister group ascidians in analysis of the urochordate miRNA repertoire, we find that 11 highly conserved bilaterian miRNA families have been lost or derived to the point they are not recognizable in urochordates and a further 4 of these families are absent in larvaceans. Subsequent to this loss/derivation, at least 29 novel miRNA families have been acquired in larvaceans. This suggests a profound reorganization of the miRNA repertoire integral to evolution in the urochordate lineage.     

EST analysis of mRNAs expressed in neurula of Chinese amphioxus

Amphioxus, a cephalochordate, is the closest living relative to the vertebrates. In order to investigate the molecular mechanisms of the early embryogenesis of amphioxus, we constructed a neurula embryo cDNA library of Chinese amphioxus (Branchiostoma belcheri tsingtauense) and generated 5235 expressed sequenced tags in the present study. The initial ESTs consisted of 638 clusters and 1855 singletons, which revealed approximately 2493 unique genes in the data set. Of these sequences, 35.52% ESTs matched to known genes, 12.76% matched to other ESTs, and 51.71% had no match to any known sequences in GenBank. Interestingly we found homologous genes related to neural development and human disease. Bioinformatic analysis showed the direct evidence that the gene homologue found only in vertebrates in previous studies also exists in the amphioxus genome. This study provides a preliminary view of the gene information involved in the development of neurula embryos of Chinese amphioxus and helps our understanding of vertebrate evolution at gene level.