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.     

Class II cytokine receptors and their ligands: key antiviral and inflammatory modulators

Class II cytokine receptors were originally defined on the basis of sequence homologies in the extracellular domains of receptors for interferons (IFNs) and interleukin-10 (IL-10), and the ligands, known as class II cytokines, also have a common structure. More recently, a series of new receptors and cytokines that belong to this family have been discovered. The therapeutic potential of the 'old' members of this family, IFNs and IL-1, is recognized in the clinic, and the existence of structurally related molecules is raising expectations for additional clinical applications. In this review, I discuss both structural and biological data that are emerging about this family of receptors and ligands, to highlight the potential applications of modulating the activity of these cytokines.     

The gamma-chain cytokine/receptor system in fish: more ligands and receptors

The mammalian gamma-chain (γC) cytokine family consists of interleukin (IL)-2, IL-4, IL-7, IL-9, IL-15 and IL-21. They signal through a receptor complex containing the common γC and a private alpha chain, and in the case of IL-2 and IL-15 an additional common IL-2/15Rβ chain. Deficiency of γC signalling in mammals prevents CD4⁺ T cells from developing effector functions and CD8⁺ T cells from developing immunological memory. Thus γC cytokines are critical for the generation and peripheral homeostasis of naïve and memory T cells. This review will give an update on the γC ligands and receptor subunits in fish, and also present some new data on the cloning and expression of a second γC and two IL-2Rβ chains in rainbow trout Oncorhynchus mykiss. In recent years, aided by the availability of sequenced fish genomes and expressed sequence tag databases, five of the six mammalian γC cytokines and their cognate receptors have been discovered in fish, with only the IL-9/IL-9R homologues apparently absent. Paralogues have been discovered in diploid fish and all the receptors described in the tetraploid rainbow trout, including γC itself, IL-2Rβ, IL-4Rα, IL-13Rα1, IL-13Rα2 and IL-2/15Rα, have duplicates. As a consequence of the teleost and salmonid whole genome duplications, even more paralogues may yet be discovered. Some of the paralogues have changes in domain structures and show differential expression and modulation, suggesting the potential for a change in function. Functional characterisation of fish γC cytokines is beginning but made more difficult by the co-existence of so many paralogues of the ligands and their receptors. Initial functional studies have shown that fish γC cytokines can modulate the expression of key cytokines (e.g. interferon-γ, IL-10 and IL-22) of the adaptive immune response, and may thus have promise as adjuvants to improve vaccination efficiency in fish.     

Evolution of trace amine associated receptor (TAAR) gene family in vertebrates: lineage-specific expansions and degradations of a second class of vertebrate chemosensory receptors expressed in the olfactory epithelium

The trace amine-associated receptors (TAARs) form a specific family of G protein-coupled receptors in vertebrates. TAARs were initially considered neurotransmitter receptors, but recent study showed that mouse TAARs function as chemosensory receptors in the olfactory epithelium. To clarify the evolutionary dynamics of the TAAR gene family in vertebrates, near-complete repertoires of TAAR genes and pseudogenes were identified from the genomic assemblies of 4 teleost fishes (zebrafish, fugu, stickleback, and medaka), western clawed frogs, chickens, 3 mammals (humans, mice, and opossum), and sea lampreys. Database searches revealed that fishes had many putatively functional TAAR genes (13-109 genes), whereas relatively small numbers of TAAR genes (3-22 genes) were identified in tetrapods. Phylogenetic analysis of these genes indicated that the TAAR gene family was subdivided into 5 subfamilies that diverged before the divergence of ray-finned fishes and tetrapods. In tetrapods, virtually all TAAR genes were located in 1 specific region of their genomes as a gene cluster; however, in fishes, TAAR genes were scattered throughout more than 2 genomic locations. This possibly reflects a whole-genome duplication that occurred in the common ancestor of ray-finned fishes. Expression analysis of zebrafish and stickleback TAAR genes revealed that many TAARs in these fishes were expressed in the olfactory organ, suggesting the relatively high importance of TAARs as chemosensory receptors in fishes. A possible evolutionary history of the vertebrate TAAR gene family was inferred from the phylogenetic and comparative genomic analyses.     

Toll-like receptors: a family of pattern-recognition receptors in mammals

The innate immune system uses a variety of germline-encoded pattern-recognition receptors that recognize conserved microbial structures or pathogen-associated molecular patterns, such as those that occur in the bacterial cell-wall components peptidoglycan and lipopolysaccharide. Recent studies have highlighted the importance of Toll-like receptors (TLRs) as a family of pattern-recognition receptors in mammals that can discriminate between chemically diverse classes of microbial products. First identified on the basis of sequence similarity with the Drosophila protein Toll, TLRs are members of an ancient superfamily of proteins, which includes related proteins in invertebrates and plants. TLRs activate innate immune defense reactions, such as the release of inflammatory cytokines, but increasing evidence supports an additional critical role for TLRs in orchestrating the development of adaptive immune responses. The sequence similarity between the intracellular domains of the TLRs and the mammalian interleukin-1 and interleukin-18 cytokine receptors reflects the use of a common intracellular signal-transduction cascade triggered by these receptor classes. But more recent findings have demonstrated that there are in fact TLR-specific signaling pathways and cellular responses. Thus, TLRs function as sentinels of the mammalian immune system that can discriminate between diverse pathogen-associated molecular patterns and then elicit pathogen-specific cellular immune responses.     

Comparative analysis of Phytophthora genes encoding secreted proteins reveals conserved synteny and lineage-specific gene duplications and deletions

Comparative analysis of two Phytophthora genomes revealed overall colinearity in four genomic regions consisting of a 1.5-Mb sequence of Phytophthora sojae and a 0.9-Mb sequence of P. ramorum. In these regions with conserved synteny, the gene order is largely similar; however, genome rearrangements also have occurred. Deletions and duplications often were found in association with genes encoding secreted proteins, including effectors that are important for interaction with host plants. Among secreted protein genes, different evolutionary patterns were found. Elicitin genes that code for a complex family of highly conserved Phytophthora-specific elicitors show conservation in gene number and order, and often are clustered. In contrast, the race-specific elicitor gene Avrlb-1 appeared to be missing from the region with conserved synteny, as were its five homologs that are scattered over the four genomic regions. Some gene families encoding secreted proteins were found to be expanded in one species compared with the other. This could be the result of either repeated gene duplications in one species or specific deletions in the other. These different evolutionary patterns may shed light on the functions of these secreted proteins in the biology and pathology of the two Phytophthora spp.     

Comparative genomic analysis of mitogen activated protein kinase gene family in grapevine

Mitogen activated protein kinases (MAPKs) are important proteins involved in the signal transduction of extracellular information to intracellular targets, and play a crucial role in the response to biotic and abiotic stresses. Although Arabidopsis MAPKs are used for identification of the putative MAPKs in higher plants, no grapevine MAPK gene nomenclature has yet been appeared in the literature. In this study, we have identified 12 members of grapevine MAPK gene (VvMPK) family via In-silico analysis of current grapevine genome database. The structural comparison of 12 VvMPKs through the analysis of chromosome locations, sequence annotation and paralogous gene pair indicated that VvMPKs have evolved by segmental duplication, rather than by tandem amplification. Although further functional analysis of VvMPKs through in vivo and in vivo experiments will be required, our study provides the basis for future research on the diverse signaling pathways medicated by MAPKs in grapevine.     

Comparative genomic and expression analysis of the conserved NTPDase gene family in Xenopus

The purines, ATP and adenosine, are important signaling molecules in the nervous system. ATP is sequentially degraded to adenosine by the ectonucleotidase proteins. The NTPDase (or CD39) family is a subfamily of these enzymes, which consists of nine members in mammals. In Xenopus embryos, we have shown that ATP, and its antagonist adenosine, regulate the rundown of swimming and we therefore proposed that ectonucleotidase proteins are key regulators of locomotor activity. Here, we report the cloning of all nine members of the NTPDase family in Xenopus laevis and Xenopus tropicalis. Our phylogenetic analysis shows that this family is highly conserved between the frog species and also during vertebrate evolution. In the adult frog, NTPDase genes are broadly expressed. During development, all NTPDase genes, except for NTPDase8, are expressed and display a distinct specific expression pattern, suggesting potentially different functions of these proteins during embryogenesis of X. laevis.     

Large-scale, lineage-specific expansion of a bric-a-brac/tramtrack/broad complex ubiquitin-ligase gene family in rice

Selective ubiquitination of proteins is directed by diverse families of ubiquitin-protein ligases (or E3s) in plants. One important type uses Cullin-3 as a scaffold to assemble multisubunit E3 complexes containing one of a multitude of bric-a-brac/tramtrack/broad complex (BTB) proteins that function as substrate recognition factors. We previously described the 80-member BTB gene superfamily in Arabidopsis thaliana. Here, we describe the complete BTB superfamily in rice (Oryza sativa spp japonica cv Nipponbare) that contains 149 BTB domain-encoding genes and 43 putative pseudogenes. Amino acid sequence comparisons of the rice and Arabidopsis superfamilies revealed a near equal repertoire of putative substrate recognition module types. However, phylogenetic comparisons detected numerous gene duplication and/or loss events since the rice and Arabidopsis BTB lineages split, suggesting possible functional specialization within individual BTB families. In particular, a major expansion and diversification of a subset of BTB proteins containing Meprin and TRAF homology (MATH) substrate recognition sites was evident in rice and other monocots that likely occurred following the monocot/dicot split. The MATH domain of a subset appears to have evolved significantly faster than those in a smaller core subset that predates flowering plants, suggesting that the substrate recognition module in many monocot MATH-BTB E3s are diversifying to ubiquitinate a set of substrates that are themselves rapidly changing. Intriguing possibilities include pathogen proteins attempting to avoid inactivation by the monocot host.     

Comparative genomic analysis of the Sm gene family in rice and maize

Sm proteins are a group of ubiquitous ring-shaped oligomers that function in multiple aspects of RNA metabolism. However, until this study, no comprehensive study incorporating phylogeny, chromosomal location, gene organization, adaptive evolution, expression profiling and functional networks has been reported for rice and maize. In this study, twenty-five and thirty-three Sm genes have been identified in rice and maize, respectively. Phylogenetic analyses identified eighteen gene groups. Results by gene locations indicated that segmental duplication contributes to the expansion of this gene family in rice and maize. Gene organization and motif compositions of the Sm members are highly conserved in each group, indicative of their functional conservation. Expression profiles have provided insights into the possible functional divergence among members of the Sm gene family. Adaptive evolution analyses suggested that purifying selection was the main force driving Sm evolution, but some critical sites might be responsible for functional divergence. In addition, four hundred and seventy-nine interactions were identified by functional network analyses, and most of which were associated with binding, cellular macromolecule biosynthesis, pre-mRNA processing and transferase activity. Overall, the data contribute to a better understanding of the complexity of Sm gene family in rice and maize and will provide a solid foundation for future functional studies.     

Comparative analysis of genomic sequence and expression of a lipid transfer protein gene family in winter barley

Three genomic clones (gb/*4.2, gb/f4.6, gb/t4.9) have been isolatedfrom a barley [Hordeum vulgare L.) cv. Igri genomic libraryusing a cDNA clone of the low temperature responsive gene b/r4.1as the probe. The genomic clones are closely related and fromsequence homology are believed to be three further members ofa barley gene family encoding lipid transfer proteins. One clone,gb/t4.9, corresponds to a second cDNA clone (Jb/f4.9). Northernblot analysis revealed that the three genes reported here areinduced by low temperature and drought; however, the relativelevel of response to these stimuli was shown to differ betweenthe genes. Further, the developmental response of gb/?4.2 isdifferent from that found for other members of the gene family.The degree of low temperature induction of the blt4 gene familyappears to be cultivar-dependent. Comparative analysis of the putative promoter regions of thegenomic clones indicates three regions that show over 80 sequenceconservation between gb/f4.9 and gb/f4.6. Also identified areputative cis acting elements previously identified in ABA responsivegenes and low temperature inducible genes. The possible functionalsignificance of the conserved regions and elements is discussedin relation to the expression of the blt4 gene family. Key words: Barley, gene expression, genomic sequence, lipid transfer protein, low temperature     

The type 1 (EGFR-related) family of growth factor receptors and their ligands

This review considers the biology of the type 1 growth factor receptor family which is increasingly recognised as important in the control of normal cell proliferation and in the pathogenesis of human cancer. The family currently comprises three closely related members: the epidermal growth factor (EGF) receptor, c-erbB-2 and c-erbB-3, all of which show abnormalities of expression in various human tumours. The family of factors related to EGF has also expanded recently and now includes transforming growth factor alpha, heparin-binding EGF, amphiregulin, cripto and heregulin, as well as several other potential ligands for the c-erbB2-2 receptor. The involvement of these receptors and growth factors in human cancer has implications for the design of novel forms of therapy for cancer, and we review recent advances and future avenues for investigation.     

Divergent evolution of the myosin heavy chain gene family in fish and tetrapods: evidence from comparative genomic analysis.

Myosin heavy chain genes (MYHs) are the most important functional domains of myosins, which are highly conserved throughout evolution. The human genome contains 15 MYHs, whereas the corresponding number in teleost appears to be much higher. Although teleosts comprise more than one-half of all vertebrate species, our knowledge of MYHs in teleosts is rather limited. A comprehensive analysis of the torafugu (Takifugu rubripes) genome database enabled us to detect at least 28 MYHs, almost twice as many as in humans. RT-PCR revealed that at least 16 torafugu MYH representatives (5 fast skeletal, 3 cardiac, 2 slow skeletal, 1 superfast, 2 smooth, and 3 nonmuscle types) are actually transcribed. Among these, MYH(M743-2) and MYH(M5) of fast and slow skeletal types, respectively, are expressed during development of torafugu embryos. Syntenic analysis reveals that torafugu fast skeletal MYHs are distributed across five genomic regions, three of which form clusters. Interestingly, while human fast skeletal MYHs form one cluster, its syntenic region in torafugu is duplicated, although each locus contains just a single MYH in torafugu. The results of the syntenic analysis were further confirmed by corresponding analysis of MYHs based on databases from Tetraodon, zebrafish, and medaka genomes. Phylogenetic analysis suggests that fast skeletal MYHs evolved independently in teleosts and tetrapods after fast skeletal MYHs had diverged from four ancestral MYHs.     

Interactions of Ly49 family receptors with MHC class I ligands in trans and cis

The Ly49A NK cell receptor interacts with MHC class I (MHC-I) molecules on target cells and negatively regulates NK cell-mediated target cell lysis. We have recently shown that the MHC-I ligand-binding capacity of the Ly49A NK cell receptor is controlled by the NK cells' own MHC-I. To see whether this property was unique to Ly49A, we have investigated the binding of soluble MHC-I multimers to the Ly49 family receptors expressed in MHC-I-deficient and -sufficient C57BL/6 mice. In this study, we confirm the binding of classical MHC-I to the inhibitory Ly49A, C and I receptors, and demonstrate that detectable MHC-I binding to MHC-I-deficient NK cells is exclusively mediated by these three receptors. We did not detect significant multimer binding to stably transfected or NK cell-expressed Ly49D, E, F, G, and H receptors. Yet, we identified the more distantly related Ly49B and Ly49Q, which are not expressed by NK cells, as two novel MHC-I receptors in mice. Furthermore, we show using MHC-I-sufficient mice that the NK cells' own MHC-I significantly masks the Ly49A and Ly49C, but not the Ly49I receptor. Nevertheless, Ly49I was partly masked on transfected tumor cells, suggesting that the structure of Ly49I is compatible in principal with cis binding of MHC-I. Finally, masking of Ly49Q by cis MHC-I was minor, whereas masking of Ly49B was not detected. These data significantly extend the MHC-I specificity of Ly49 family receptors and show that the accessibility of most, but not all, MHC-I-binding Ly49 receptors is modulated by the expression of MHC-I in cis.     

Molecular analysis of an MHC class II deficiency patient reveals a novel mutation in the RFX5 gene

 Patients suffering from major histocompatibility complex (MHC) class II deficiency, a rare primary immunodeficiency, are characterized by a lack of MHC class II expression which is the result of defects in trans-acting factors. At least four complementation groups, A, B, C, and D, can be discerned. The gene affected in group C patients is known to be RFX5 and encodes one of the subunits of the multimeric phosphoprotein complex, RFX. In the present study we fused fibroblasts of a recently identified MHC class II deficiency patient, OSE, with fibroblasts derived from patients representative of each of the four complementation groups. Transient heterokaryon analysis indicated that OSE belonged to complementation group C. Furthermore, transfection of wild-type RFX5 cDNA into OSE fibroblasts resulted in restoration of the defect. Mutation analysis revealed that the RFX5 mRNA lacked four nucleotides and that this deletion was the consequence of a G to A transition in a splice acceptor site. Genomic oligotyping demonstrated that OSE was homozygous for the splice site mutation. Received: 15 July 1998 / Revised: 3 September 1998     

Comparative genomic analysis reveals occurrence of genetic recombination in virulent Cryptosporidium hominis subtypes and telomeric gene duplications in Cryptosporidium parvum

Background

Cryptosporidium hominis is a dominant species for human cryptosporidiosis. Within the species, IbA10G2 is the most virulent subtype responsible for all C. hominis–associated outbreaks in Europe and Australia, and is a dominant outbreak subtype in the United States. In recent yearsIaA28R4 is becoming a major new subtype in the United States. In this study, we sequenced the genomes of two field specimens from each of the two subtypes and conducted a comparative genomic analysis of the obtained sequences with those from the only fully sequenced Cryptosporidium parvum genome.

Results

Altogether, 8.59-9.05 Mb of Cryptosporidium sequences in 45–767 assembled contigs were obtained from the four specimens, representing 94.36-99.47% coverage of the expected genome. These genomes had complete synteny in gene organization and 96.86-97.0% and 99.72-99.83% nucleotide sequence similarities to the published genomes of C. parvum and C. hominis, respectively. Several major insertions and deletions were seen between C. hominis and C. parvum genomes, involving mostly members of multicopy gene families near telomeres. The four C. hominis genomes were highly similar to each other and divergent from the reference IaA25R3 genome in some highly polymorphic regions. Major sequence differences among the four specimens sequenced in this study were in the 5′ and 3′ ends of chromosome 6 and the gp60 region, largely the result of genetic recombination.

Conclusions

The sequence similarity among specimens of the two dominant outbreak subtypes and genetic recombination in chromosome 6, especially around the putative virulence determinant gp60 region, suggest that genetic recombination plays a potential role in the emergence of hyper-transmissible C. hominis subtypes. The high sequence conservation between C. parvum and C. hominis genomes and significant differences in copy numbers of MEDLE family secreted proteins and insulinase-like proteases indicate that telomeric gene duplications could potentially contribute to host expansion in C. parvum.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1517-1) contains supplementary material, which is available to authorized users.