Evolutionary conservation and selection of human disease gene orthologs in the rat and mouse genomes

Background

Model organisms have contributed substantially to our understanding of the etiology of human disease as well as having assisted with the development of new treatment modalities. The availability of the human, mouse and, most recently, the rat genome sequences now permit the comprehensive investigation of the rodent orthologs of genes associated with human disease. Here, we investigate whether human disease genes differ significantly from their rodent orthologs with respect to their overall levels of conservation and their rates of evolutionary change.

Results

Human disease genes are unevenly distributed among human chromosomes and are highly represented (99.5%) among human-rodent ortholog sets. Differences are revealed in evolutionary conservation and selection between different categories of human disease genes. Although selection appears not to have greatly discriminated between disease and non-disease genes, synonymous substitution rates are significantly higher for disease genes. In neurological and malformation syndrome disease systems, associated genes have evolved slowly whereas genes of the immune, hematological and pulmonary disease systems have changed more rapidly. Amino-acid substitutions associated with human inherited disease occur at sites that are more highly conserved than the average; nevertheless, 15 substituting amino acids associated with human disease were identified as wild-type amino acids in the rat. Rodent orthologs of human trinucleotide repeat-expansion disease genes were found to contain substantially fewer of such repeats. Six human genes that share the same characteristics as triplet repeat-expansion disease-associated genes were identified; although four of these genes are expressed in the brain, none is currently known to be associated with disease.

Conclusions

Most human disease genes have been retained in rodent genomes. Synonymous nucleotide substitutions occur at a higher rate in disease genes, a finding that may reflect increased mutation rates in the chromosomal regions in which disease genes are found. Rodent orthologs associated with neurological function exhibit the greatest evolutionary conservation; this suggests that rodent models of human neurological disease are likely to most faithfully represent human disease processes. However, with regard to neurological triplet repeat expansion-associated human disease genes, the contraction, relative to human, of rodent trinucleotide repeats suggests that rodent loci may not achieve a 'critical repeat threshold' necessary to undergo spontaneous pathological repeat expansions. The identification of six genes in this study that have multiple characteristics associated with repeat expansion-disease genes raises the possibility that not all human loci capable of facilitating neurological disease by repeat expansion have as yet been identified.     

VH gene segments in the mouse and human genomes

We have examined the mouse genome sequence to determine its VH gene segment repertoire. In all, 141 segments are mapped to a 3 Mb region of chromosome 12. There is evidence that 92 of these are functional in the mouse strain used for the genome sequence, C57BL/6J; 12 are functional in other mouse strains, and 37 are pseudogenes. The mouse VH gene segment repertoire is therefore twice the size of that in humans. The mouse and human loci bear no large-scale similarity to each other. The 104 functional segments belong to one of the 15 known sequence subgroups, which have been further clustered into eight sets here. Seven of these sets, comprising 101 sequences, are related to five of the human VH families and have the same canonical structures in their hypervariable regions. Duplication of members of one set in the distal half of the locus is mainly responsible for the larger size of the mouse repertoire. Phylogenetic analysis of the VH segments indicates that most of the sequences in the human and mouse VH loci have arisen subsequent to the divergence of the two organisms from their common ancestor.     

A cluster of proton/amino acid transporter genes in the human and mouse genomes

We recently cloned and functionally characterized two novel proton/amino acid transporters (PAT1 and PAT2) from mouse. Here we report the isolation of the corresponding cDNAs of the human orthologues and one additional mouse and human PAT-like transporter cDNA, designated PAT3. The PAT proteins comprise 470 to 483 amino acids. The mouse PAT3 mRNA is expressed in testis of adult mice. In the human and mouse genomes the genes of the PAT transporters (designated SLC36A1-3 and Slc36a1-3, respectively) are clustered on human chromosome 5q33.1 and in the syntenic region of mouse chromosome 11B1.3. PAT-like transporter genes are present as well in the genomes of other eukaryotic organisms such as Drosophila melanogaster and Caenorhabditis elegans. For the PAT3 subtype transporter, we could not yet identify its function. The human PAT1 and PAT2 transporters when functionally expressed in Xenopus laevis oocytes show characteristics similar to those of their mouse counterparts.     

Comparative analysis of processed pseudogenes in the mouse and human genomes

Pseudogenes are important resources in evolutionary and comparative genomics because they provide molecular records of the ancient genes that existed in the genome millions of years ago. We have systematically identified approximately 5000 processed pseudogenes in the mouse genome, and estimated that approximately 60% are lineage specific, created after the mouse and human diverged. In both mouse and human genomes, similar types of genes give rise to many processed pseudogenes. These tend to be housekeeping genes, which are highly expressed in the germ line. Ribosomal-protein genes, in particular, form the largest sub-group. The processed pseudogenes in the mouse occur with a distinctly different chromosomal distribution than LINEs or SINEs - preferentially in GC-poor regions. Finally, the age distribution of mouse-processed pseudogenes closely resembles that of LINEs, in contrast to human, where the age distribution closely follows Alus (SINEs).     

Substitution pattern of the CArG element in human and mouse genomes

cis-Elements CArG bound by serum response factor (SRF) are presently being intensively studied, but little is known about the substitution pattern of functional CArG elements. Here, we have performed the first evolutionary analysis of CArGome in the human and mouse genome through bioinformatic methods and statistical tests. We calculated the substitution rate at each site of the functional CArG elements. The results showed that the core sites of the functional CArG elements evolved faster than did the background DNA, indicating that these sites were likely to evolve under positive selection. Moreover, a strong TATA "motif" was evident in the core region within the functional CArG elements in both human and mouse promoters. This motif could probably be a major contribution to the formation of the spatial structure, which was important for CArG-SRF recognition. Thus, the study further revealed the sequence character and substitution pattern of CArG elements and provided useful information for the study of the SRF-binding efficiencies of CArG promoters in functional assays.     

Correlation between sequence conservation and structural thermodynamics of microRNA precursors from human,mouse, and chicken genomes

Background

Previous studies have shown that microRNA precursors (pre-miRNAs) have considerably more stable secondary structures than other native RNAs (tRNA, rRNA, and mRNA) and artificial RNA sequences. However, pre-miRNAs with ultra stable secondary structures have not been investigated. It is not known if there is a tendency in pre-miRNA sequences towards or against ultra stable structures? Furthermore, the relationship between the structural thermodynamic stability of pre-miRNA and their evolution remains unclear.

Results

We investigated the correlation between pre-miRNA sequence conservation and structural stability as measured by adjusted minimum folding free energies in pre-miRNAs isolated from human, mouse, and chicken. The analysis revealed that conserved and non-conserved pre-miRNA sequences had structures with similar average stabilities. However, the relatively ultra stable and unstable pre-miRNAs were more likely to be non-conserved than pre-miRNAs with moderate stability. Non-conserved pre-miRNAs had more G+C than A+U nucleotides, while conserved pre-miRNAs contained more A+U nucleotides. Notably, the U content of conserved pre-miRNAs was especially higher than that of non-conserved pre-miRNAs. Further investigations showed that conserved and non-conserved pre-miRNAs exhibited different structural element features, even though they had comparable levels of stability.

Conclusions

We proposed that there is a correlation between structural thermodynamic stability and sequence conservation for pre-miRNAs from human, mouse, and chicken genomes. Our analyses suggested that pre-miRNAs with relatively ultra stable or unstable structures were less favoured by natural selection than those with moderately stable structures. Comparison of nucleotide compositions between non-conserved and conserved pre-miRNAs indicated the importance of U nucleotides in the pre-miRNA evolutionary process. Several characteristic structural elements were also detected in conserved pre-miRNAs.

     

The six hyaluronidase-like genes in the human and mouse genomes.

The human genome contains six hyaluronidase-like genes. Three genes (HYAL1, HYAL2 and HYAL3) are clustered on chromosome 3p21.3, and another two genes (HYAL4 and PH-20/SPAM1) and one expressed pseudogene (HYALP1) are similarly clustered on chromosome 7q31.3. The extensive homology between the different hyaluronidase genes suggests ancient gene duplication, followed by en masse block duplication, events that occurred before the emergence of modern mammals. Very recently we have found that the mouse genome also has six hyaluronidase-like genes that are also grouped into two clusters of three, in regions syntenic with the human genome. Surprisingly, the mouse ortholog of HYALP1 does not contain any mutations, and unlike its human counterpart may actually encode an active enzyme. Hyal-1 is the only hyaluronidase in mammalian plasma and urine, and is also found at high levels in major organs such as liver, kidney, spleen, and heart. A model is proposed suggesting that Hyal-2 and Hyal-1 are the major mammalian hyaluronidases in somatic tissues, and that they act in concert to degrade high molecular weight hyaluronan to the tetrasaccharide. Twenty-kDa hyaluronan fragments are generated at the cell surface in unique endocytic vesicles resulting from digestion by the glycosylphosphatidyl-inositol-anchored Hyal-2, transported intracellularly by an unknown process, and then further digested by Hyal-1. The two beta-exoglycosidases, beta-glucuronidase and beta-N-acetyl glucosaminidase, remove sugars from reducing termini of hyaluronan oligomers, and supplement the hyaluronidases in the catabolism of hyaluronan.     

Density peaks of paralog pairs in human and mouse genomes

Paralog gene trees, which reflect the increase of genomic complexity in the evolution, can be complicated and ambiguous. A simpler complementary approach is analysis of density distribution of paralog pairs. It can reveal general features of genome evolution, which may be hidden in the forest of gene trees. It is known that distribution of human paralog pairs along the axis of protein divergence between pair members forms two main peaks. Here I show that there are three main peaks in the mouse genome. Thus, the multimodality of paralog pair distribution seems to be a fundamental feature of mammalian genomes. Despite the great diversity of domains presented in small amounts or in multidomain architectures with a few predominant domains, both in human and mouse the first peak consists mostly of gene pairs with zinc finger domains or olfactory receptor domain. In the mouse the olfactory receptor predominates, which stipulates the three-peak distribution (since in the olfactory receptors the second peak is closer to the first peak than in other genes). The mammalian-wide zinc finger orthologs are biased towards the second peak. Thus, the marsupial orthologs are nearly absent in the first peak of human and mouse. The gene pairs in the first peak show a lower ratio of nonsynonymous to synonymous substitutions, which suggests that their evolution is more constrained. The plausible explanation is that they are in subfunctionalization state (partition of initial function of ancestral gene), whereas the second peak contains gene pairs that are already in neofunctionalization state (acquiring of novel functions). These data suggest that the adaptive radiation of mammals was accompanied by a burst of duplication of zinc finger genes, which are located in the first (most recent) peak of paralog pairs.     

An analysis on gene architecture in human and mouse genomes

A comparative genome analysis on exon-intron distribution profiles is performed for human and mouse genomes to deduce similarities and differences between them. Interestingly, both in human and mouse genomes, the total length in introns and intergenic DNA on each chromosome is significantly correlated to the chromosome size. The results presented provide a framework for understanding the nature and patterns of exon-intron length distributions, the constraints on them and their role in genome design and evolution.     

Selective constraint in intergenic regions of human and mouse genomes.

We aligned and analyzed 100 pairs of complete, orthologous intergenic regions from the human and mouse genomes (average length approximately 12 000 nucleotides). The alignments alternate between highly similar segments and dissimilar segments, indicating a wide variation of selective constraint. The average number of selectively constrained nucleotides within a mammalian intergenic region is at least 2000. This is threefold higher than within a nematode intergenic region and at least twofold higher than the number of selectively constrained nucleotides coding for an average protein. Because mammals possess only two- to threefold more proteins than Caenorhabditis elegans, the higher complexity of mammals might be primarily because of the functioning of intergenic DNA.     

Organization and stability of endogenous xenotropic murine leukemia virus proviral DNA in mouse genomes.

In a series of blot hybridization experiments, using a xenotropic envelope probe and restriction enzymes known to cut xenotropic proviral DNA a single time (EcoRI) or not at all (HindIII), we have studied the organization and relationship of endogenous xenotropic env-related sequences in various mouse strains. Multiple copies (18 to 28) of xenotropic env-reactive fragments were found in all mouse DNAs after digestion with either HindIII or EcoRI, and the majority of fragments were of sizes compatible with their origin from full-length proviral DNA. Five HindIII and five EcoRI restriction fragments were common to all inbred mouse DNAs tested. In addition, each strain exhibited unique characteristic xenotropic env-reactive bands; these bands were remarkably stable during many years of inbreeding. The cleavage patterns characteristic of each strain were also useful for showing genealogical relatedness among the various inbred mice.     

The scattered distribution of actin genes in the mouse and human genomes.

A hamster actin cDNA probe was used to localize actin genes on the major components of mouse and human DNAs, namely on the four families of fragments forming the bulk of these genomes. Over 20 EcoRI fragments hybridizing the probe could be detected; a different subset of these fragments was found in each component. Since the fragment families forming the major components of the mouse and human DNAs derive from very long chromosomal segments, the isochores , the presence of actin genes on all components provides evidence for their dispersion in both genomes. In situ hybridization of 125I-labeled probe to metaphase chromosomes in the presence of dextran sulfate confirmed this dispersion by showing that the 29-30 actin gene sites so identified are distributed on almost all chromosomes. Moreover, some human actin genes could be mapped on specific chromosomal segments; in particular, one gene was localized on the long arm of the X chromosome. Finally, three different mouse actin genes were isolated from a recombinant DNA library and previously investigated interspersed repeated sequences were identified in the vicinity of these genes.     

Computational comparison of two mouse draft genomes and the human golden path

Background  

The availability of both mouse and human draft genomes has marked the beginning of a new era of comparative mammalian genomics. The two available mouse genome assemblies, from the public mouse genome sequencing consortium and Celera Genomics, were obtained using different clone libraries and different assembly methods.