Related Gene Family Research Articles

No evidence for association of autism with rare heterozygous point mutations in Contactin-Associated Protein-Like 2 (CNTNAP2), or in Other Contactin-Associated Proteins or Contactins.

Contactins and Contactin-Associated Proteins, and Contactin-Associated Protein-Like 2 (CNTNAP2) in particular, have been widely cited as autism risk genes based on findings from homozygosity mapping, molecular cytogenetics, copy number variation analyses, and both common and rare single nucleotide association studies. However, data specifically with regard to the contribution of heterozygous single nucleotide variants (SNVs) have been inconsistent. In an effort to clarify the role of rare point mutations in CNTNAP2 and related gene families, we have conducted targeted next-generation sequencing and evaluated existing sequence data in cohorts totaling 2704 cases and 2747 controls. We find no evidence for statistically significant association of rare heterozygous mutations in any of the CNTN or CNTNAP genes, including CNTNAP2, placing marked limits on the scale of their plausible contribution to risk.

CRISPR-Cas9-mediated single-gene and gene family disruption in Trypanosoma cruzi.

ABSTRACTTrypanosoma cruzi is a protozoan parasite of humans and animals, affecting 10 to 20 million people and innumerable animals, primarily in the Americas. Despite being the largest cause of infection-induced heart disease worldwide, even among the neglected tropical diseases (NTDs) T. cruzi is considered one of the least well understood and understudied. The genetic complexity of T. cruzi as well as the limited set of efficient techniques for genome engineering contribute significantly to the relative lack of progress in and understanding of this pathogen. Here, we adapted the CRISPR-Cas9 system for the genetic engineering of T. cruzi, demonstrating rapid and efficient knockout of multiple endogenous genes, including essential genes. We observed that in the absence of a template, repair of the Cas9-induced double-stranded breaks (DSBs) in T. cruzi occurs exclusively by microhomology-mediated end joining (MMEJ) with various-sized deletions. When a template for DNA repair is provided, DSB repair by homologous recombination is achieved at an efficiency several orders of magnitude higher than that in the absence of CRISPR-Cas9-induced DSBs. We also demonstrate the high multiplexing capacity of CRISPR-Cas9 in T. cruzi by knocking down expression of an enzyme gene family consisting of 65 members, resulting in a significant reduction of enzymatic product with no apparent off-target mutations. Lastly, we show that Cas9 can mediate disruption of its own coding sequence, rescuing a growth defect in stable Cas9-expressing parasites. These results establish a powerful new tool for the analysis of gene functions in T. cruzi, enabling the study of essential genes and their functions and analysis of the many large families of related genes that occupy a substantial portion of the T. cruzi genome.

Hypothetical Structure of Prostate-Specific Antigen

Prostate specific antigen (PSA) is the premier biomarker for diagnosis and monitoring of prostate cancer. PSA is a trypsin-like serine protease and belongs to a family of related genes which are known as the kallikrein-related peptidases. The …

Functional evolution of Erg potassium channel gating reveals an ancient origin for I Kr

Mammalian Ether-a-go-go related gene (Erg) family voltage-gated K(+) channels possess an unusual gating phenotype that specializes them for a role in delayed repolarization. Mammalian Erg currents rectify during depolarization due to rapid, voltage-dependent inactivation, but rebound during repolarization due to a combination of rapid recovery from inactivation and slow deactivation. This is exemplified by the mammalian Erg1 channel, which is responsible for IKr, a current that repolarizes cardiac action potential plateaus. The Drosophila Erg channel does not inactivate and closes rapidly upon repolarization. The dramatically different properties observed in mammalian and Drosophila Erg homologs bring into question the evolutionary origins of distinct Erg K(+) channel functions. Erg channels are highly conserved in eumetazoans and first evolved in a common ancestor of the placozoans, cnidarians, and bilaterians. To address the ancestral function of Erg channels, we identified and characterized Erg channel paralogs in the sea anemone Nematostella vectensis. N. vectensis Erg1 (NvErg1) is highly conserved with respect to bilaterian homologs and shares the IKr-like gating phenotype with mammalian Erg channels. Thus, the IKr phenotype predates the divergence of cnidarians and bilaterians. NvErg4 and Caenorhabditis elegans Erg (unc-103) share the divergent Drosophila Erg gating phenotype. Phylogenetic and sequence analysis surprisingly indicates that this alternate gating phenotype arose independently in protosomes and cnidarians. Conversion from an ancestral IKr-like gating phenotype to a Drosophila Erg-like phenotype correlates with loss of the cytoplasmic Ether-a-go-go domain. This domain is required for slow deactivation in mammalian Erg1 channels, and thus its loss may partially explain the change in gating phenotype.

Library Generation by Gene Shuffling

This unit describes the process of gene shuffling, also known as sexual PCR. Gene shuffling is a facile method for the generation of sequence libraries containing the information from a family of related genes. Essentially, related genes are fragmented by DNase I digestion and reassembled by primer-less PCR. The resulting chimeric genes can then be screened or selected for a desired function.

EggNOG v4.0: nested orthology inference across 3686 organisms

With the increasing availability of various ‘omics data, high-quality orthology assignment is crucial for evolutionary and functional genomics studies. We here present the fourth version of the eggNOG database (available at http://eggnog.embl.de) that derives nonsupervised orthologous groups (NOGs) from complete genomes, and then applies a comprehensive characterization and analysis pipeline to the resulting gene families. Compared with the previous version, we have more than tripled the underlying species set to cover 3686 organisms, keeping track with genome project completions while prioritizing the inclusion of high-quality genomes to minimize error propagation from incomplete proteome sets. Major technological advances include (i) a robust and scalable procedure for the identification and inclusion of high-quality genomes, (ii) provision of orthologous groups for 107 different taxonomic levels compared with 41 in eggNOGv3, (iii) identification and annotation of particularly closely related orthologous groups, facilitating analysis of related gene families, (iv) improvements of the clustering and functional annotation approach, (v) adoption of a revised tree building procedure based on the multiple alignments generated during the process and (vi) implementation of quality control procedures throughout the entire pipeline. As in previous versions, eggNOGv4 provides multiple sequence alignments and maximum-likelihood trees, as well as broad functional annotation. Users can access the complete database of orthologous groups via a web interface, as well as through bulk download.

Caste-Specific Expression Patterns of Immune Response and Chemosensory Related Genes in the Leaf-Cutting Ant, Atta vollenweideri

Leaf-cutting ants are evolutionary derived social insects with elaborated division of labor and tremendous colony sizes with millions of workers. Their social organization is mainly based on olfactory communication using different pheromones and is promoted by a pronounced size-polymorphism of workers that perform different tasks within the colony. The size polymorphism and associated behaviors are correlated to distinct antennal lobe (AL) phenotypes. Two worker phenotypes differ in number of olfactory glomeruli in the AL and the presence or absence of an extremely large glomerulus (macroglomerulus), involved in trail-pheromone reception. The males' AL contains three macroglomeruli which are presumably involved in detection of sex-pheromone components. We investigated the antennal transcriptome data of all major castes (males, queens and workers) and two worker subcastes (large and tiny workers). In order to identify putative odorant receptor genes involved in pheromone detection, we identified differentially expressed odorant receptor genes (OR-genes) using custom microarrays. In total, we found 185 OR-gene fragments that are clearly related to ORs and we identified orthologs for 70 OR-genes. Among them one OR-gene differs in relative expression between the two worker subcastes by a factor of >3 and thus is a very promising candidate gene for the trail-pheromone receptor. Using the relative expression of OR-genes in males versus queens, we identified 2 candidates for sex-pheromone receptor genes in males. In addition, we identified genes from all other chemosensory related gene families (13 chemosensory protein genes, 8 odorant binding protein genes, 2 sensory-neuron membrane protein genes, 7 ionotropic receptor genes, 2 gustatory receptor genes), and we found ant-specific expansions in the chemosensory protein gene family. In addition, a large number of genes involved in immune defense exhibited differential expression across the three different castes, and some genes even between the two worker subcastes.

English

Designed degenerate primers unlike conventional primers are superior in matching and amplification of large number of genes, from related gene families. DPPrimer tool was designed to predict primers for PCR amplification of homologous gene from related or diverse plant species. The key features of this tool include platform independence and user friendliness in primer design. Embedded features such as search for functional domains, similarity score selection and phylogebetic tree further enhance the user friendliness of DPPrimer tool. Performance of DPPrimer tool was evaluated by successful PCR amplification of ADP-glucose phosphorylase genes from wheat, barley and rice.AvailabilityDPPrimer is freely accessible at http://202.141.12.147/DGEN_tool/index.html

Transcriptome Analysis of Androgenic Gland for Discovery of Novel Genes from the Oriental River Prawn, Macrobrachium nipponense, Using Illumina Hiseq 2000

BackgroundThe oriental river prawn, Macrobrachium nipponense, is an important aquaculture species in China, even in whole of Asia. The androgenic gland produces hormones that play crucial roles in sexual differentiation to maleness. This study is the first de novo M. nipponense transcriptome analysis using cDNA prepared from mRNA isolated from the androgenic gland. Illumina/Solexa was used for sequencing.Methodology and Principal FindingThe total volume of RNA sample was more than 5 ug. We generated 70,853,361 high quality reads after eliminating adapter sequences and filtering out low-quality reads. A total of 78,408 isosequences were obtained by clustering and assembly of the clean reads, producing 57,619 non-redundant transcripts with an average length of 1244.19 bp. In total 70,702 isosequences were matched to the Nr database, additional analyses were performed by GO (33,203), KEGG (17,868), and COG analyses (13,817), identifying the potential genes and their functions. A total of 47 sex-determination related gene families were identified from the M. nipponense androgenic gland transcriptome based on the functional annotation of non-redundant transcripts and comparisons with the published literature. Furthermore, a total of 40 candidate novel genes were found, that may contribute to sex-determination based on their extremely high expression levels in the androgenic compared to other sex glands,. Further, 437 SSRs and 65,535 high-confidence SNPs were identified in this EST dataset from which 14 EST-SSR markers have been isolated.ConclusionOur study provides new sequence information for M. nipponense, which will be the basis for further genetic studies on decapods crustaceans. More importantly, this study dramatically improves understanding of sex-determination mechanisms, and advances sex-determination research in all crustacean species. The huge number of potential SSR and SNP markers isolated from the transcriptome may shed the lights on research in many fields, including the evolution and molecular ecology of Macrobrachium species.

Integrated evolutionary analysis of human miRNA gene clusters and families implicates evolutionary relationships

Many microRNAs (miRNAs) are clustered on chromosomes and co-transcribed as polycistronic transcripts. Here, an integrated evolutionary analysis of human miRNA gene clusters and families was performed. Generally, miRNA gene clusters include 2–8 members, but some larger clusters have been found to have more members (over 40 miRNAs). 62.22% of them have been shown to be involved in homologous miRNA genes, including multicopy pre-miRNAs and sense/antisense homologous miRNAs. Multicopy pre-miRNAs can enrich the distribution and relationship between miRNA clusters and families. An miRNA family may be located in one or more clusters, and a cluster may be involved in one or more families. Members of different families have been shown to be prone to appear in clusters, and vice versa. Reconstructed phylogenetic trees and networks may indicate potential evolutionary relationships, which also indicate duplication history in specific related gene clusters and families. Related miRNA families are always found to share common target mRNAs and biological pathways. Some clusters containing non-homologous miRNAs also tend to be clustered together as well as homologous miRNAs. In the present work, it is shown that homologous miRNAs are prone to appear in clusters based on functional and evolutionary pressures. The phenomenon of miRNA clusters containing homologous or genetic relationships is quite common. The integrative evolutionary analysis will provide more potential evolutionary and functional relationships between homologous and clustered miRNAs.

CorSig: A General Framework for Estimating Statistical Significance of Correlation and Its Application to Gene Co-Expression Analysis

With the rapid increase of omics data, correlation analysis has become an indispensable tool for inferring meaningful associations from a large number of observations. Pearson correlation coefficient (PCC) and its variants are widely used for such purposes. However, it remains challenging to test whether an observed association is reliable both statistically and biologically. We present here a new method, CorSig, for statistical inference of correlation significance. CorSig is based on a biology-informed null hypothesis, i.e., testing whether the true PCC (ρ) between two variables is statistically larger than a user-specified PCC cutoff (τ), as opposed to the simple null hypothesis of ρ = 0 in existing methods, i.e., testing whether an association can be declared without a threshold. CorSig incorporates Fisher's Z transformation of the observed PCC (r), which facilitates use of standard techniques for p-value computation and multiple testing corrections. We compared CorSig against two methods: one uses a minimum PCC cutoff while the other (Zhu's procedure) controls correlation strength and statistical significance in two discrete steps. CorSig consistently outperformed these methods in various simulation data scenarios by balancing between false positives and false negatives. When tested on real-world Populus microarray data, CorSig effectively identified co-expressed genes in the flavonoid pathway, and discriminated between closely related gene family members for their differential association with flavonoid and lignin pathways. The p-values obtained by CorSig can be used as a stand-alone parameter for stratification of co-expressed genes according to their correlation strength in lieu of an arbitrary cutoff. CorSig requires one single tunable parameter, and can be readily extended to other correlation measures. Thus, CorSig should be useful for a wide range of applications, particularly for network analysis of high-dimensional genomic data.Software AvailabilityA web server for CorSig is provided at http://202.127.200.1:8080/probeWeb. R code for CorSig is freely available for non-commercial use at http://aspendb.uga.edu/downloads.

The functional significance of cattle major histocompatibility complex class I genetic diversity.

Current concerns about food security highlight the importance of maintaining productive and disease-resistant livestock populations. Major histocompatibility complex (MHC) class I genes have a central role in immunity. A high level of diversity in these genes allows populations to survive despite exposure to rapidly evolving pathogens. This review aims to describe the key features of MHC class I genetic diversity in cattle and to discuss their role in disease resistance. Discussion centers on data derived from the cattle genome sequence and studies addressing MHC class I gene expression and function. The impact of intensive selection on MHC diversity is also considered. A high level of complexity in MHC class I genes and functionally related gene families is revealed. This highlights the need for increased efforts to determine key genetic components that govern cattle immune responses to disease, which is increasingly important in the face of changing human and environmental demands.

Illumina-based de novo transcriptome sequencing and analysis of Amanita exitialis basidiocarps

Amanita exitialis is a lethal mushroom that was first discovered in Guangdong Province, China. The high content of amanitin in its basidiocarps makes it lethal to humans. To comprehensively characterize the A. exitialis transcriptome and analyze the Amanita toxins as well as their related gene family, transcriptome sequencing of A. exitialis was performed using Illumina HiSeq 2000 technology. A total of 25,563,688 clean reads were collected and assembled into 62,137 cDNA contigs with an average length of 481bp and N50 length of 788bp. A total of 27,826 proteins and 39,661 unigenes were identified among the assembled contigs. All of the unigenes were classified into 166 functional categories for understanding the gene functions and regulation pathways. The genes contributing to toxic peptide biosynthesis were analyzed. From this set, eleven gene sequences encoding the toxins or related cyclic peptides were discovered in the transcriptome. Three of these sequences matched the peptide toxins α-amanitin, β-amanitin, and phallacidin, while others matched amanexitide and seven matched unknown peptides. All of the genes encoding peptide toxins were confirmed by polymerase chain reaction (PCR) in A. exitialis, and the phylogenetic relationships among these proprotein sequences were discussed. The gene polymorphism and degeneracy of the toxin encoding sequences were found and analyzed. This study provides the first primary transcriptome of A. exitialis, which provided comprehensive gene expression information on the lethal amanitas at the transcriptional level, and could lay a strong foundation for functional genomics studies in those fungi.

Cloning and tissue-specific expression of predicted Pisum sativum actin isoform PEAc14-1.

Actin is an essential component of the cytoskeleton in eukaryotes and participates in many important subcellular processes, making it critical in plant development and morphogenesis (Pei et al. 2012). Under the control of actin-binding proteins, the length, polarity, stability, and three-dimensional structure of the actin microfilament are dynamic (Henty et al. 2011; Fan et al. 2013). The actin microfilament is linked to other organelles in cells to play important roles in a variety of subcellular processes, such as cell division, cell elongation, endocytosis, cell signaling, gravity sensing, apical growth, organelle movement, and programmed cell death (Pei et al. 2012; Li et al. 2012; Smertenko and Franklin-Tong 2011). First isolated from muscle cells, actin constitutes approximately 20% of total cell protein. The existence of plant actin was proved for the first time by Yan and Shi (1963). In plants, actin plays an essential role in many processes; many genes of actin isoforms have been cloned from plants and characterized (Zhang and Liu 2006; Zhao et al. 2012). In general, plants express several isoforms of actin encoded by a family of related genes. In addition, the actin isoforms are often differentially expressed in many tissues and organs, and their functions are also different (Jiang and Zhao 2002; Diaz-Camino et al. 2005); for example, the Arabidopsis thaliana actin gene family contains 10 isovariants and can be divided into reproductive and vegetative classes (Pei et al. 2012). Cao et al. (1994) constructed a cDNA library of Pisum sativum tendrils and obtained 18 highly conservative actin genes divided into three classes: PEAcI, PEAcII, and PEAcIII. These isoforms tend to be expressed in rapidly developing organs, and their activity

Derepression of Cancer/Testis Antigens in cancer is associated with distinct patterns of DNA Hypomethylation

BackgroundThe Cancer/Testis Antigens (CTAs) are a heterogeneous group of proteins whose expression is typically restricted to the testis. However, they are aberrantly expressed in most cancers that have been examined to date. Broadly speaking, the CTAs can be divided into two groups: the CTX antigens that are encoded by the X-linked genes and the non-X CT antigens that are encoded by the autosomes. Unlike the non-X CTAs, the CTX antigens form clusters of closely related gene families and their expression is frequently associated with advanced disease with poorer prognosis. Regardless however, the mechanism(s) underlying their selective derepression and stage-specific expression in cancer remain poorly understood, although promoter DNA demethylation is believed to be the major driver.MethodsHere, we report a systematic analysis of DNA methylation profiling data from various tissue types to elucidate the mechanism underlying the derepression of the CTAs in cancer. We analyzed the methylation profiles of 501 samples including sperm, several cancer types, and their corresponding normal somatic tissue types.ResultsWe found strong evidence for specific DNA hypomethylation of CTA promoters in the testis and cancer cells but not in their normal somatic counterparts. We also found that hypomethylation was clustered on the genome into domains that coincided with nuclear lamina-associated domains (LADs) and that these regions appeared to be insulated by CTCF sites. Interestingly, we did not observe any significant differences in the hypomethylation pattern between the CTAs without CpG islands and the CTAs with CpG islands in the proximal promoter.ConclusionOur results corroborate that widespread DNA hypomethylation appears to be the driver in the derepression of CTA expression in cancer and furthermore, demonstrate that these hypomethylated domains are associated with the nuclear lamina-associated domains (LADS). Taken together, our results suggest that wide-spread methylation changes in cancer are linked to derepression of germ-line-specific genes that is orchestrated by the three dimensional organization of the cancer genome.

The B-hordein prolamin family of barley

The spectrum of B-hordein prolamins and genes in the single barley cultivar Barke is described from an in silico analysis of 1452 B-hordein ESTs and available genomic DNA. Eleven unique B-hordein proteins are derived from EST contigs. Ten contigs encode apparent full-length B-hordeins and the eleventh contains a premature stop codon that will lead to a truncated B-hordein. The 11 sequences are placed within the two previously described classes, i.e., the B1- and B3-type B-hordeins. The number of ESTs assigned to each sequence is used as an estimate of relative gene transcription and expression. Three of the sequences account for 79% of the total ESTs, with one sequence comprises 32% of the total ESTs and has a variant C-terminus caused by an undefined sequence change history near the 3' coding terminus. The 70× difference in EST distribution among sequences points to the importance of understanding differential rates of expression within closely related gene families. Analysis of available genomic sequences confirms the EST assembly and reveals one full-length and two partial sequences of pseudogenes as evidenced by no matching ESTs for the sequences and premature stop codons and frame shifts.

Exploiting a fast neutron mutant genetic resource in Pisum sativum (pea) for functional genomics.

A fast neutron (FN)-mutagenised population was generated in Pisum sativum L. (pea) to enable the identification and isolation of genes underlying traits and processes. Studies of several phenotypic traits have clearly demonstrated the utility of the resource by associating gene deletions with phenotype followed by functional tests exploiting additional mutant sources, from both induced and natural variant germplasm. For forward genetic screens, next generation sequencing methodologies provide an opportunity for identifying genes associated with deletions rapidly and systematically. The application of rapid reverse genetic screens of the fast neutron mutant pea population supports conclusions on the frequency of deletions based on phenotype alone. These studies also suggest that large deletions affecting one or more loci can be non-deleterious to the pea genome, yielding mutants that could not be obtained by other means. Deletion mutants affecting genes associated with seed metabolism and storage are providing unique opportunities to identify the products of complex and related gene families, and to study the downstream consequences of such deletions.

Expression of MS4A and TMEM176 Genes in Human B Lymphocytes

The MS4A gene family in humans includes CD20 and at least 15 other genes. CD20 exists as homo-oligomers in the plasma membrane, however different MS4A proteins expressed in the same cell may hetero-oligomerize. Given the importance of CD20 in B-cell function and as a therapeutic target, we sought to explore the potential for CD20 hetero-oligomerization with other MS4A proteins. We investigated expression in primary human B-cells of the four MS4A genes previously shown to be expressed in human B-cell lines (MS4A4A, MS4A6A, MS4A7, MS4A8B), as well as two genes comprising the closely related TMEM176 gene family, with a view to identifying candidates for future investigation at the protein level. TMEM176A and TMEM176B transcripts were either not detected, or were detected at relatively low levels in a minority of donor B-cell samples. MS4A4A and MS4A8B transcripts were not detected in any normal B-cell sample. MS4A6A and MS4A7 transcripts were detected at low levels in most samples, however the corresponding proteins were not at the plasma membrane when expressed as GFP conjugates in BJAB cells. We also examined expression of these genes in chronic lymphocytic leukemia (CLL), and found that it was similar to normal B-cells with two exceptions. First, whereas MS4A4A expression was undetected in normal B-cells, it was expressed in 1/14 CLL samples. Second, compared to expression levels in normal B-cells, MS4A6A transcripts were elevated in 4/14 CLL samples. In summary, none of the MS4A/TMEM176 genes tested was expressed at high levels in normal or in most CLL B-cells. MS4A6A and MS4A7 were expressed at low levels in most B-cell samples, however the corresponding proteins may not be positioned at the plasma membrane. Altogether, these data suggest that CD20 normally does not form hetero-oligomers with other MS4A proteins and that there are unlikely to be other MS4A proteins in CLL that might provide useful alternate therapeutic targets.

RISC in PD: the impact of microRNAs in Parkinson's disease cellular and molecular pathogenesis

Parkinson's disease (PD) is a debilitating neurodegenerative disease characterized primarily by the selective death of dopaminergic (DA) neurons in the substantia nigra pars compacta of the midbrain. Although several genetic forms of PD have been identified, the precise molecular mechanisms underlying DA neuron loss in PD remain elusive. In recent years, microRNAs (miRNAs) have been recognized as potent post-transcriptional regulators of gene expression with fundamental roles in numerous biological processes. Although their role in PD pathogenesis is still a very active area of investigation, several seminal studies have contributed significantly to our understanding of the roles these small non-coding RNAs play in the disease process. Among these are studies which have demonstrated specific miRNAs that target and down-regulate the expression of PD-related genes as well as those demonstrating a reciprocal relationship in which PD-related genes act to regulate miRNA processing machinery. Concurrently, a wealth of knowledge has become available regarding the molecular mechanisms that unify the underlying etiology of genetic and sporadic PD pathogenesis, including dysregulated protein quality control by the ubiquitin-proteasome system and autophagy pathway, activation of programmed cell death, mitochondrial damage and aberrant DA neurodevelopment and maintenance. Following a discussion of the interactions between PD-related genes and miRNAs, this review highlights those studies which have elucidated the roles of these pathways in PD pathogenesis. We highlight the potential of miRNAs to serve a critical regulatory role in the implicated disease pathways, given their capacity to modulate the expression of entire families of related genes. Although few studies have directly linked miRNA regulation of these pathways to PD, a strong foundation for investigation has been laid and this area holds promise to reveal novel therapeutic targets for PD.

Dynamic evolution of mir-17–92 gene cluster and related miRNA gene families in vertebrates

mir-17-92 gene cluster is widely distributed in vertebrates and plays an important role in regulating multiple biological processes. Its dysregulation may be associated with risk of some human cancers. The microRNA (miRNA) members are identified in the three gene families: mir-17, mir-19 and mir-25. Herein we attempted to understand the evolutionary processes and patterns in vertebrates. The three miRNA gene families showed difference in distribution, number of miRNA genes and clustered miRNA genes in the five animal species. Compared to other related gene clusters, mir-17-92 cluster was well-conserved and had more abundant roles in multiple biological processes. These clustered miRNAs showed inconsistent nucleotide divergence patterns across different animal species, even between homologous miRNA genes. Simultaneously, they also indicated inconsistent expression patterns although they were co-transcribed as a polycistronic transcript. Phylogenetic tree based on human pre-miRNA sequences showed that mir-19 gene family was an older miRNA species, while tree based on miRNA gene cluster indicated evolutionary positions of different animal species. The study shows dynamic evolution of the mir-17-92 gene cluster and related miRNA gene families across vertebrates, which may be derived from potential functional implication. miRNA gene cluster should be a better phylogenetic marker than a single miRNA gene to reveal functional and evolutionary relationships.