Comparative genomic analysis using the UCSC genome browser

Dysbindin Domain-Containing 1 in Prostate Cancer: New Insights into Bioinformatic Validation of Molecular and Immunological Features

Prostate cancer (PCa) is one of the most prevalent cancers in men, yet its pathogenic pathways remain poorly understood. Transcriptomics and high-throughput sequencing can help uncover cancer diagnostic targets and understand biological circuits. Using prostate adenocarcinoma (PRAD) datasets of various web-based applications (GEPIA, UALCAN, cBioPortal, SR Plot, hTFtarget, Genome Browser, and MetaCore), we found that upregulated dysbindin domain-containing 1 (DBNDD1) expression in primary prostate tumors was strongly correlated with pathways involving the cell cycle, mitotic in KEGG, WIKI, and REACTOME database, and transcription factor-binding sites with the DBNDD1 gene in prostate samples. DBNDD1 gene expression was influenced by sample type, cancer stage, and promoter methylation levels of different cancers, such as PRAD, liver hepatocellular carcinoma (LIHC), and lung adenocarcinoma (LUAD). Regulation of glycogen synthase kinase (GSK)-3β in bipolar disorder and ATP/ITP/GTP/XTP/TTP/CTP/UTP metabolic pathways was closely correlated with the DBNDD1 gene and its co-expressed genes in PCa. DBNDD1 gene expression was positively associated with immune infiltration of B cells, Myeloid-derived suppressor cell (MDSC), M2 macrophages, andneutrophil, whereas negatively correlated with CD8+ T cells, T follicular helper cells, M1 macrophages, and NK cells in PCa. These findings suggest that DBNDD1 may serve as a viable prognostic marker not only for early-stage PCa but also for immunotherapies.

Prenatal ultrasonographic features in Blomstrand osteochondrodysplasia: Antenatal case series confirmed by postmortem radiology and molecular diagnosis

OBJECTIVE

Blomstrand osteochondrodysplasia (BOCD, MIM #215045) is an ultrarare lethal skeletal dysplasia (LSD) perinatally, characterized by extremely advanced bone maturation, generalized osteosclerosis, and severe tetramicromelia caused by biallelic loss-of-function mutations in the parathyroid hormone receptor-1 gene (PTHR1). We aim to describe prenatal ultrasonographic features in a retrospective fetal case series of BOCD and emphasize the importance of multidisciplinary antenatal evaluation of LSDs to improve the differential diagnosis.

METHOD

Prenatal ultrasound findings of five fetal cases diagnosed with BOCD between 2000 and 2019 in the Prenatal Diagnosis Unit and Medical Genetics were reviewed, along with postmortem examination results and confirmatory molecular results.

RESULTS

All fetuses presented with severe sonographic findings of LSDs comprising tetramicromelia, thoracic hypoplasia, and retro-micrognathia. Marked cervical hyperextension was present in three fetuses. Flared metaphyses were prenatally identified in only one fetus. X-rays of four fetuses evaluated postmortem showed advanced bone maturation, generalized osteosclerosis, and dumbbell-like appearance of long bones due to metaphyseal enlargement.

CONCLUSION

The presence of retro-micrognathia along with a protruding tongue and severe metaphyseal flaring can suggest a diagnosis of BOCD, when prenatal ultrasound findings are indicative for LSD. The diagnosis can be ascertained through postmortem clinical and radiological evaluation and/or molecular testing.

Polyploidy among salmonid aldehyde dehydrogenase genes and proteins

Bioinformatic analyses of salmon (Salmo salar) ALDH amino acid sequences supported the presence of at least 30 ALDH genes, which is more than for any other higher vertebrate and is greater than the 19 human ALDH genes currently reported. These included 8 polyploid ALDH genes and proteins: ALDH1A2 (chromosomes 11 and 26); ALDH1L2 (chromosomes 7 and 17); ALDH2, encoding mitochondrial ALDH2 (chromosomes 2 and 5); ALDH3A2 (chromosomes 4, 9 and 20), for which evidence for 5 genes was obtained; ALDH3B1 (chromosomes 3, 6 and 24); ALDH4A1 (chromosomes 12 and 22); ALDH6A1 (chromosomes 1, 6 and 15); and ALDH18A1 (chromosomes 19 and 28). In contrast, 7 salmon ALDH gene families (ALDH1A1, ALDH1A3, ALDH5, ALDH7, ALDH8, ALDH9 and ALDH16) possessed only one gene family member. Phylogenetic studies of salmon and rainbow trout ALDH3A2 genes and proteins suggested that salmonid gene tetraploidy has occurred in at least 2 distinct stages of ALDH3A2 gene evolution.

Comparative and evolutionary studies of ALDH18A1 genes and proteins

Vertebrate ALDH18A1 genes encode a bifunctional mitochondrial enzyme, catalyzing a 2-step conversion of glutamate to glutamyl semialdehyde, subsequently converted into proline, ornithine and arginine. Bioinformatic analyses of vertebrate and invertebrate genomes were undertaken using known ALDH18A1 amino acid sequences. G5K (glutamyl kinase) and GPR (glutamyl phosphate reductase) domain sequences were identified for all vertebrate and invertebrate genomes examined, whereas bacterial sequences encoded separate enzymes. Vertebrate ALDH18A1 (also called P5CS) sequences were highly conserved throughout vertebrate evolution. A mechanism for generating two major vertebrate ALDH18A1 isoforms is proposed with 'a' isoform containing Asn239-Val240 with wide tissue expression, whereas the 'b' isoform lacking the dipeptide has been reported in gut tissues. Phylogenetic analyses describe the relationships and potential origins of the ALDH18A1 gene during vertebrate and invertebrate evolution and a proposal for generating the bifunctional vertebrate and invertebrate ALDH18A1 gene from a bacterial operon (proBA) encoding G5K and GPR. A more recent Aldh18a1 gene duplication event has apparently occurred with a primordial rat genome.

Vanno: a visualization-aided variant annotation tool

Next-generation sequencing (NGS) technologies have revolutionized the field of genetics and are trending toward clinical diagnostics. Exome and targeted sequencing in a disease context represent a major NGS clinical application, considering its utility and cost-effectiveness. With the ongoing discovery of disease-associated genes, various gene panels have been launched for both basic research and diagnostic tests. However, the fundamental inconsistencies among the diverse annotation sources, software packages, and data formats have complicated the subsequent analysis. To manage disease-associated NGS data, we developed Vanno, a Web-based application for in-depth analysis and rapid evaluation of disease-causative genome sequence alterations. Vanno integrates information from biomedical databases, functional predictions from available evaluation models, and mutation landscapes from TCGA cancer types. A highly integrated framework that incorporates filtering, sorting, clustering, and visual analytic modules is provided to facilitate exploration of oncogenomics datasets at different levels, such as gene, variant, protein domain, or three-dimensional structure. Such design is crucial for the extraction of knowledge from sequence alterations and translating biological insights into clinical applications. Taken together, Vanno supports almost all disease-associated gene tests and exome sequencing panels designed for NGS, providing a complete solution for targeted and exome sequencing analysis. Vanno is freely available at http://cgts.cgu.edu.tw/vanno.

An alignment-free method to find and visualise rearrangements between pairs of DNA sequences

Species evolution is indirectly registered in their genomic structure. The emergence and advances in sequencing technology provided a way to access genome information, namely to identify and study evolutionary macro-events, as well as chromosome alterations for clinical purposes. This paper describes a completely alignment-free computational method, based on a blind unsupervised approach, to detect large-scale and small-scale genomic rearrangements between pairs of DNA sequences. To illustrate the power and usefulness of the method we give complete chromosomal information maps for the pairs human-chimpanzee and human-orangutan. The tool by means of which these results were obtained has been made publicly available and is described in detail.

Aldehyde dehydrogenase homologous folate enzymes: Evolutionary switch between cytoplasmic and mitochondrial localization

Cytosolic and mitochondrial 10-formyltetrahydrofolate dehydrogenases are products of separate genes in vertebrates but only one such gene is present in invertebrates. There is a significant degree of sequence similarity between the two enzymes due to an apparent origin of the gene for the mitochondrial enzyme (ALDH1L2) from the duplication of the gene for the cytosolic enzyme (ALDH1L1). The primordial ALDH1L gene originated from a natural fusion of three unrelated genes, one of which was an aldehyde dehydrogenase. Such structural organization defined the catalytic mechanism of these enzymes, which is similar to that of aldehyde dehydrogenases. Here we report the analysis of ALDH1L1 and ALDH1L2 genes from different species and their phylogeny and evolution. We also performed sequence and structure comparison of ALDH1L enzymes possessing aldehyde dehydrogenase catalysis to those lacking this feature in an attempt to explain mechanistic differences between cytoplasmic ALDH1L1 and mitochondrial ALDH1L2 enzymes and to better understand their functional roles.

Comparative and evolutionary studies of vertebrate ALDH1A-like genes and proteins

Vertebrate ALDH1A-like genes encode cytosolic enzymes capable of metabolizing all-trans-retinaldehyde to retinoic acid which is a molecular 'signal' guiding vertebrate development and adipogenesis. Bioinformatic analyses of vertebrate and invertebrate genomes were undertaken using known ALDH1A1, ALDH1A2 and ALDH1A3 amino acid sequences. Comparative analyses of the corresponding human genes provided evidence for distinct modes of gene regulation and expression with putative transcription factor binding sites (TFBS), CpG islands and micro-RNA binding sites identified for the human genes. ALDH1A-like sequences were identified for all mammalian, bird, lizard and frog genomes examined, whereas fish genomes displayed a more restricted distribution pattern for ALDH1A1 and ALDH1A3 genes. The ALDH1A1 gene was absent in many bony fish genomes examined, with the ALDH1A3 gene also absent in the medaka and tilapia genomes. Multiple ALDH1A1-like genes were identified in mouse, rat and marsupial genomes. Vertebrate ALDH1A1, ALDH1A2 and ALDH1A3 subunit sequences were highly conserved throughout vertebrate evolution. Comparative amino acid substitution rates showed that mammalian ALDH1A2 sequences were more highly conserved than for the ALDH1A1 and ALDH1A3 sequences. Phylogenetic studies supported an hypothesis for ALDH1A2 as a likely primordial gene originating in invertebrate genomes and undergoing sequential gene duplication to generate two additional genes, ALDH1A1 and ALDH1A3, in most vertebrate genomes.

SP transcription factor paralogs and DNA-binding sites coevolve and adaptively converge in mammals and birds

Functional modification of regulatory proteins can affect hundreds of genes throughout the genome, and is therefore thought to be almost universally deleterious. This belief, however, has recently been challenged. A potential example comes from transcription factor SP1, for which statistical evidence indicates that motif preferences were altered in eutherian mammals. Here, we set out to discover possible structural and theoretical explanations, evaluate the role of selection in SP1 evolution, and discover effects on coregulatory proteins. We show that SP1 motif preferences were convergently altered in birds as well as mammals, inducing coevolutionary changes in over 800 regulatory regions. Structural and phylogenic evidence implicates a single causative amino acid replacement at the same SP1 position along both lineages. Furthermore, paralogs SP3 and SP4, which coregulate SP1 target genes through competitive binding to the same sites, have accumulated convergent replacements at the homologous position multiple times during eutherian and bird evolution, presumably to preserve competitive binding. To determine plausibility, we developed and implemented a simple model of transcription factor and binding site coevolution. This model predicts that, in contrast to prevailing beliefs, even small selective benefits per locus can drive concurrent fixation of transcription factor and binding site mutants under a broad range of conditions. Novel binding sites tend to arise de novo, rather than by mutation from ancestral sites, a prediction substantiated by SP1-binding site alignments. Thus, multiple lines of evidence indicate that selection has driven convergent evolution of transcription factors along with their binding sites and coregulatory proteins.

Comparative genomics, molecular evolution and computational modeling of ALDH1B1 and ALDH2

Vertebrate ALDH2 genes encode mitochondrial enzymes capable of metabolizing acetaldehyde and other biological aldehydes in the body. Mammalian ALDH1B1, another mitochondrial enzyme sharing 72% identity with ALDH2, is also capable of metabolizing acetaldehyde but has a tissue distribution and pattern of activity distinct from that of ALDH2. Bioinformatic analyses of several vertebrate genomes were undertaken using known ALDH2 and ALDH1B1 amino acid sequences. Phylogenetic analysis of many representative vertebrate species (including fish, amphibians, birds and mammals) indicated the presence of ALDH1B1 in many mammalian species and in frogs (Xenopus tropicalis); no evidence was found for ALDH1B1 in the genomes of birds, reptiles or fish. Predicted vertebrate ALDH2 and ALDH1B1 subunit sequences and structures were highly conserved, including residues previously shown to be involved in catalysis and coenzyme binding for human ALDH2. Studies of ALDH1B1 sequences supported the hypothesis that the ALDH1B1 gene originated in early vertebrates from a retrotransposition of the vertebrate ALDH2 gene. Given the high degree of similarity between ALDH2 and ALDH1B1, it is surprising that individuals with an inactivating mutation in ALDH2 (ALDH2*2) do not exhibit a compensatory increase in ALDH1B1 activity. We hypothesized that the similarity between the two ALDHs would allow for dominant negative heterotetramerization between the inactive ALDH2 mutants and ALDH1B1. Computational-based molecular modeling studies examining predicted protein-protein interactions indicated that heterotetramerization between ALDH2 and ALDH1B1 subunits was highly probable and may partially explain a lack of compensation by ALDH1B1 in ALDH2(∗)2 individuals.

CSF1R mutations identified in three families with autosomal dominantly inherited leukoencephalopathy

Genetic and phenotypic heterogeneities are considerably high in adult-onset leukoencephalopathy, in which comprehensive mutational analyses of the candidate genes by conventional methods are too laborious. We applied exome sequencing to conduct a comprehensive mutational analysis of genes for autosomal dominant leukoencephalopathies. Genomic DNA samples from four patients of three families with autosomal dominantly inherited adult-onset leukodystrophy were subjected to exome sequencing. On the basis of the results, 21 patients with adult-onset sporadic leukodystrophy and one patient with pathologically proven HDLS were additionally screened for CSF1R mutations. Exome sequencing identified heterozygous CSF1R mutations (p.I794T and p.R777W) in two families. I794T has recently been reported as a causative mutation for hereditary diffuse leukoencephalopathy with spheroids (HDLS), and R777W is a novel mutation. Although mutational analysis of CSF1R in 21 sporadic cases revealed no mutations, another novel CSF1R mutation, p.C653Y, was identified in one patient with autopsy-proven HDSL. These variants were located in the PTK domain where the causative mutations cluster. Functional prediction of the mutant CSF1R as well as cross-species conservation of the affected amino acids supports the notion that these variants are pathogenic for HDLS. Exome sequencing is useful for a comprehensive mutational analysis of causative genes for hereditary leukoencephalopathies, and CSF1R should be considered a candidate gene for patients with autosomal dominant leukoencephalopathies.

Ank3-dependent SVZ niche assembly is required for the continued production of new neurons

The rodent subventricular/subependymal zone (SVZ/SEZ) houses neural stem cells (NSCs) that generate olfactory bulb interneurons. It is unclear how the SVZ environment sustains neuronal production into adulthood. We discovered that the adapter molecule Ankyrin-3 (Ank3) is specifically upregulated in ventricular progenitors destined to become ependymal cells, but not in NSCs, and is required for SVZ niche assembly through progenitor lateral adhesion. Furthermore, we found that Ank3 expression is controlled by Foxj1, a transcriptional regulator of multicilia formation, and genetic deletion of this pathway led to complete loss of SVZ niche structure. Interestingly, radial glia continued to transition into postnatal NSCs without this niche. However, inducible deletion of Foxj1-Ank3 from mature SVZ ependyma resulted in dramatic depletion of neurogenesis. Targeting a pathway regulating ependymal organization/assembly and showing its requirement for new neuron production, our results have important implications for environmental control of adult neurogenesis and harvesting NSCs for replacement therapy.

Use of a Drosophila genome-wide conserved sequence database to identify functionally related cis-regulatory enhancers

Background: Phylogenetic footprinting has revealed that cis-regulatory enhancers consist of conserved DNA sequence clusters (CSCs). Currently, there is no systematic approach for enhancer discovery and analysis that takes full-advantage of the sequence information within enhancer CSCs. Results: We have generated a Drosophila genome-wide database of conserved DNA consisting of >100,000 CSCs derived from EvoPrints spanning over 90% of the genome. cis-Decoder database search and alignment algorithms enable the discovery of functionally related enhancers. The program first identifies conserved repeat elements within an input enhancer and then searches the database for CSCs that score highly against the input CSC. Scoring is based on shared repeats as well as uniquely shared matches, and includes measures of the balance of shared elements, a diagnostic that has proven to be useful in predicting cis-regulatory function. To demonstrate the utility of these tools, a temporally-restricted CNS neuroblast enhancer was used to identify other functionally related enhancers and analyze their structural organization. Conclusions:cis-Decoder reveals that co-regulating enhancers consist of combinations of overlapping shared sequence elements, providing insights into the mode of integration of multiple regulating transcription factors. The database and accompanying algorithms should prove useful in the discovery and analysis of enhancers involved in any developmental process. Developmental Dynamics 241:169–189, 2012. © 2011 Wiley Periodicals, Inc.

Dynamic changes in genome-wide histone H3 lysine 4 methylation patterns in response to dehydration stress in Arabidopsis thaliana

BACKGROUND

The molecular mechanisms of genome reprogramming during transcriptional responses to stress are associated with specific chromatin modifications. Available data, however, describe histone modifications only at individual plant genes induced by stress. We have no knowledge of chromatin modifications taking place at genes whose transcription has been down-regulated or on the genome-wide chromatin modification patterns that occur during the plant's response to dehydration stress.

RESULTS

Using chromatin immunoprecipitation and deep sequencing (ChIP-Seq) we established the whole-genome distribution patterns of histone H3 lysine 4 mono-, di-, and tri-methylation (H3K4me1, H3K4me2, and H3K4me3, respectively) in Arabidopsis thaliana during watered and dehydration stress conditions. In contrast to the relatively even distribution of H3 throughout the genome, the H3K4me1, H3K4me2, and H3K4me3 marks are predominantly located on genes. About 90% of annotated genes carry one or more of the H3K4 methylation marks. The H3K4me1 and H3K4me2 marks are more widely distributed (80% and 84%, respectively) than the H3K4me3 marks (62%), but the H3K4me2 and H3K4me1 levels changed only modestly during dehydration stress. By contrast, the H3K4me3 abundance changed robustly when transcripts levels from responding genes increased or decreased. In contrast to the prominent H3K4me3 peaks present at the 5'-ends of most transcribed genes, genes inducible by dehydration and ABA displayed atypically broader H3K4me3 distribution profiles that were present before and after the stress.

CONCLUSIONS

A higher number (90%) of annotated Arabidopsis genes carry one or more types of H3K4me marks than previously reported. During the response to dehydration stress the changes in H3K4me1, H3K4me2, and H3K4me3 patterns show different dynamics and specific patterns at up-regulated, down-regulated, and unaffected genes. The different behavior of each methylation mark during the response process illustrates that they have distinct roles in the transcriptional response of implicated genes. The broad H3K4me3 distribution profiles on nucleosomes of stress-induced genes uncovered a specific chromatin pattern associated with many of the genes involved in the dehydration stress response.

Visualizing genomes: techniques and challenges

As our ability to generate sequencing data continues to increase, data analysis is replacing data generation as the rate-limiting step in genomics studies. Here we provide a guide to genomic data visualization tools that facilitate analysis tasks by enabling researchers to explore, interpret and manipulate their data, and in some cases perform on-the-fly computations. We will discuss graphical methods designed for the analysis of de novo sequencing assemblies and read alignments, genome browsing, and comparative genomics, highlighting the strengths and limitations of these approaches and the challenges ahead.

Identifying and characterizing a novel protein kinase STK35L1 and deciphering its orthologs and close-homologs in vertebrates

Background

The human kinome containing 478 eukaryotic protein kinases has over 100 uncharacterized kinases with unknown substrates and biological functions. The Ser/Thr kinase 35 (STK35, Clik1) is a member of the NKF 4 (New Kinase Family 4) in the kinome with unknown substrates and biological functions. Various high throughput studies indicate that STK35 could be involved in various human diseases such as colorectal cancer and malaria.

Methodology/Principal Findings

In this study, we found that the previously published coding sequence of the STK35 gene is incomplete. The newly identified sequence of the STK35 gene codes for a protein of 534 amino acids with a N-terminal elongation of 133 amino acids. It has been designated as STK35L (STK35 long). Since it is the first of further homologous kinases we termed it as STK35L1. The STK35L1 protein (58 kDa on SDS-PAGE), but not STK35 (44 kDa), was found to be expressed in all human cells studied (endothelial cells, HeLa, and HEK cells) and was down-regulated after silencing with specific siRNA. EGFP-STK35L1 was localized in the nucleus and the nucleolus. By combining syntenic and gene structure pattern data and homology searches, two further STK35L1 homologs, STK35L2 (previously known as PDIK1L) and STK35L3, were found. All these protein kinase homologs were conserved throughout the vertebrates. The STK35L3 gene was specifically lost during placental mammalian evolution. Using comparative genomics, we have identified orthologous sets of these three protein kinases genes and their possible ancestor gene in two sea squirt genomes.

Conclusions/Significance

We found the full-length coding sequence of the STK35 gene and termed it as STK35L1. We identified a new third STK35-like gene, STK35L3, in vertebrates and a possible ancestor gene in sea squirt genome. This study will provide a comprehensive platform to explore the role of STK35L kinases in cell functions and human diseases.

In vivo selection of hematopoietic progenitor cells and temozolomide dose intensification in rhesus macaques through lentiviral transduction with a drug resistance gene

Major limitations to gene therapy using HSCs are low gene transfer efficiency and the inability of most therapeutic genes to confer a selective advantage on the gene-corrected cells. One approach to enrich for gene-modified cells in vivo is to include in the retroviral vector a drug resistance gene, such as the P140K mutant of the DNA repair enzyme O6-methylguanine-DNA methyltransferase (MGMT*). We transplanted 5 rhesus macaques with CD34+ cells transduced with lentiviral vectors encoding MGMT* and a fluorescent marker, with or without homeobox B4 (HOXB4), a potent stem cell self-renewal gene. Transgene expression and common integration sites in lymphoid and myeloid lineages several months after transplantation confirmed transduction of long-term repopulating HSCs. However, all animals showed only a transient increase in gene-marked lymphoid and myeloid cells after O6-benzylguanine (BG) and temozolomide (TMZ) administration. In 1 animal, cells transduced with MGMT* lentiviral vectors were protected and expanded after multiple courses of BG/TMZ, providing a substantial increase in the maximum tolerated dose of TMZ. Additional cycles of chemotherapy using 1,3-bis-(2-chloroethyl)-1-nitrosourea (BCNU) resulted in similar increases in gene marking levels, but caused high levels of nonhematopoietic toxicity. Inclusion of HOXB4 in the MGMT* vectors resulted in no substantial increase in gene marking or HSC amplification after chemotherapy treatment. Our data therefore suggest that lentivirally mediated gene transfer in transplanted HSCs can provide in vivo chemoprotection of progenitor cells, although selection of long-term repopulating HSCs was not seen.

Diversification of TOLLIP isoforms in mouse and man

The Toll-interacting protein TOLLIP is an ubiquitin-binding protein that interacts with several components of the Toll-like receptor signaling cascade. The canonical protein consists of three annotated domains: an N-terminal TBD-loop-coil domain that mediates protein-protein interactions, a C2 domain that targets TOLLIP to the endosome, and a CUE domain at the C-terminus that binds monoubiquitin. TOLLIP has been described primarily in trafficking of the interleukin-1 receptor (IL1R) and turnover of the interleukin-1 receptor-associated kinase (IRAK), so it is an essential regulator of inflammatory signaling. Here we describe the expression of numerous alternate transcripts from mouse and human TOLLIP, which are predicted to generate at least five variant proteins between the two species. Most of the variant proteins are predicted to have altered N-terminal domains, altered TBD-loop-coil domains, or a truncated C2 domain. A mouse-specific variant arises from an alternate termination exon, and the resulting protein lacks the CUE domain. Two transcripts arising from alternate initiating exons are highly conserved between mouse and human but exhibit different patterns of expression. The consequent protein isoforms retain (TOLLIP.A) or lack (TOLLIP.D) the protein-binding TBD, so are predicted to traffic monoubiquitinated proteins to alternate protein complexes within the endosomal compartment. In summary, the widespread and inducible expression of Tollip isoforms predicts diversification of its function in rodent and human immune systems. Alternate splicing of critical signaling molecules such as Tollip may provide one mechanism behind the broad repertoire of responses generated by cells of the innate immune system in response to infection.

Rapid evolution of mammalian X-linked testis microRNAs

BACKGROUND

MicroRNAs (miRNAs), which are small, non-coding RNAs approximately 21-nucleotides in length, have become a major focus of research in molecular biology. Mammalian miRNAs are proposed to regulate approximately 30% of all protein-coding genes. Previous studies have focused on highly conserved miRNAs, but nonconserved miRNAs represent a potentially important source of novel functionalities during evolution.

RESULTS

An analysis of the chromosome distribution of miRNAs showed higher densities of miRNAs on the X chromosome compared to the average densities on autosomes in all eight mammalian species analyzed. The distribution pattern did not, however, apply well to species beyond mammals. In addition, by comparing orthologous human and mouse miRNAs, we found that X-linked miRNAs had higher substitution rates than autosomal miRNAs. Since the highest proportion of X-linked miRNAs were found in mouse testis, we tested the hypothesis that testis miRNAs are evolving faster on the X chromosome than on autosomes. Mature X-linked testis miRNAs had an average substitution rate between mouse and human that was almost 25-fold higher than mature testis miRNAs on autosomes. In contrast, for mature miRNAs with precursors not expressed in testis, no significant difference in the substitution rate between the X chromosome and autosomes was found. Among mammals, the rapid evolution of X-linked testis miRNAs was also observed in rodents and primates.

CONCLUSION

The rapid evolution of X-linked testis miRNAs implies possible important male reproductive functions and may contribute to speciation in mammals.