Mammalian gene evolution: nucleotide sequence divergence between mouse and rat

Functional sensory circuits built from neurons of two species

A central question for regenerative neuroscience is whether synthetic neural circuits, such as those built from two species, can function in an intact brain. Here, we apply blastocyst complementation to selectively build and test interspecies neural circuits. Despite approximately 10-20 million years of evolution, and prominent species differences in brain size, rat pluripotent stem cells injected into mouse blastocysts develop and persist throughout the mouse brain. Unexpectedly, the mouse niche reprograms the birth dates of rat neurons in the cortex and hippocampus, supporting rat-mouse synaptic activity. When mouse olfactory neurons are genetically silenced or killed, rat neurons restore information flow to odor processing circuits. Moreover, they rescue the primal behavior of food seeking, although less well than mouse neurons. By revealing that a mouse can sense the world using neurons from another species, we establish neural blastocyst complementation as a powerful tool to identify conserved mechanisms of brain development, plasticity, and repair.

A Genome-Wide Profiling of Glioma Patients with an IDH1 Mutation Using the Catalogue of Somatic Mutations in Cancer Database

Simple Summary

Glioma patients that present a somatic mutation in the isocitrate dehydrogenase 1 (IDH1) gene have a significantly better prognosis and overall survival than patients with the wild-type genotype. An IDH1 mutation is hypothesized to occur early during cellular transformation and leads to further genetic instability. A genome-wide profiling of glioma patients in the Catalogue of Somatic Mutations in Cancer (COSMIC) database was performed to classify the genetic differences in IDH1-mutant versus IDH1-wildtype patients. This classification will aid in a better understanding of how this specific mutation influences the genetic make-up of glioma and the resulting prognosis. Key differences in co-mutation and gene expression levels were identified that correlate with an improved prognosis.

Abstract

Gliomas are differentiated into two major disease subtypes, astrocytoma or oligodendroglioma, which are then characterized as either IDH (isocitrate dehydrogenase)-wild type or IDH-mutant due to the dramatic differences in prognosis and overall survival. Here, we investigated the genetic background of IDH1-mutant gliomas using the Catalogue of Somatic Mutations in Cancer (COSMIC) database. In astrocytoma patients, we found that IDH1 is often co-mutated with TP53, ATRX, AMBRA1, PREX1, and NOTCH1, but not CHEK2, EGFR, PTEN, or the zinc finger transcription factor ZNF429. The majority of the mutations observed in these genes were further confirmed to be either drivers or pathogenic by the Cancer-Related Analysis of Variants Toolkit (CRAVAT). Gene expression analysis showed down-regulation of DRG2 and MSN expression, both of which promote cell proliferation and invasion. There was also significant over-expression of genes such as NDRG3 and KCNB1 in IDH1-mutant astrocytoma patients. We conclude that IDH1-mutant glioma is characterized by significant genetic changes that could contribute to a better prognosis in glioma patients.

Extra base hits: Widespread empirical support for instantaneous multiple-nucleotide changes

Despite many attempts to introduce evolutionary models that permit substitutions to instantly alter more than one nucleotide in a codon, the prevailing wisdom remains that such changes are rare and generally negligible or are reflective of non-biological artifacts, such as alignment errors. Codon models continue to posit that only single nucleotide change have non-zero rates. Here, we develop and test a simple hierarchy of codon-substitution models with non-zero evolutionary rates for only one-nucleotide (1H), one- and two-nucleotide (2H), or any (3H) codon substitutions. Using over 42, 000 empirical alignments, we find widespread statistical support for multiple hits: 61% of alignments prefer models with 2H allowed, and 23%-with 3H allowed. Analyses of simulated data suggest that these results are not likely to be due to simple artifacts such as model misspecification or alignment errors. Further modeling reveals that synonymous codon island jumping among codons encoding serine, especially along short branches, contributes significantly to this 3H signal. While serine codons were prominently involved in multiple-hit substitutions, there were other common exchanges contributing to better model fit. It appears that a small subset of sites in most alignments have unusual evolutionary dynamics not well explained by existing model formalisms, and that commonly estimated quantities, such as dN/dS ratios may be biased by model misspecification. Our findings highlight the need for continued evaluation of assumptions underlying workhorse evolutionary models and subsequent evolutionary inference techniques. We provide a software implementation for evolutionary biologists to assess the potential impact of extra base hits in their data in the HyPhy package and in the Datamonkey.org server.

Deciphering Ancestral Sex Chromosome Turnovers Based on Analysis of Male Mutation Bias

The age of sex chromosomes is commonly obtained by comparing the substitution rates of XY gametologs. Coupled with phylogenetic reconstructions, one can refine the origin of a sex chromosome system relative to specific speciation events. However, these approaches are insufficient to determine the presence and duration of ancestral sex chromosome systems that were lost in some species. In this study, we worked with genomic and transcriptomic data from mammals and squamates and analyzed the effect of male mutation bias on X-linked sequences in these groups. We searched for signatures indicating whether monotremes shared the same sex chromosomes with placental mammals or whether pleurodonts and acrodonts had a common ancestral sex chromosome system. Our analyses indicate that platypus did not share the XY chromosomes with placental mammals, in agreement with previous work. In contrast, analyses of agamids showed that this lineage maintained the pleurodont XY chromosomes for several million years. We performed multiple simulations using different strengths of male mutation bias to confirm the results. Overall, our work shows that variations in substitution rates due to male mutation bias could be applied to uncover signatures of ancestral sex chromosome systems.

Mutations in disordered proteins as early indicators of nucleic acid changes triggering speciation

In this study, we analyze the role of different structural variants of proteins in the speciation processes. We separate human and mouse proteomes (taken as a reference) into three previously defined variants of disorder: ordered proteins (ORDPs), structured proteins with intrinsically disordered protein regions (IDPRs), and intrinsically disordered proteins (IDPs). Then, using the representation we call here Forsdyke plot, we study the correlation of DNA divergence with the corresponding protein (phenotypic) divergence in the three variants, comparing human and mouse coding sequences with their homologs from 26 eukaryotes. The parameters of the correlation are related to the speciation process. We find that the three variants of disordered proteins are differently related to the speciation process. Specifically, IDPs phenotypically diverge earlier than ORDPs and IDPRs. ORDPs diverge later but are phenotypically more reactive to nucleotide mutations than IDPRs and IDPs. Finally, IDPRs appear to diverge phenotypically later than IDPs, like ORDPs, but they are prone to accept mutations with rates that are similar to those of IDPs. We conclude that IDPs are involved in the early stages of the speciation process, whereas mutations in ORDPs, once speciation is initiated, accelerate phenotypic divergence.

Male Mutation Bias Is the Main Force Shaping Chromosomal Substitution Rates in Monotreme Mammals

In many species, spermatogenesis involves more cell divisions than oogenesis, and the male germline, therefore, accumulates more DNA replication errors, a phenomenon known as male mutation bias. The extent of male mutation bias (α) is estimated by comparing substitution rates of the X, Y, and autosomal chromosomes, as these chromosomes spend different proportions of their time in the germlines of the two sexes. Male mutation bias has been characterized in placental and marsupial mammals as well as birds, but analyses in monotremes failed to detect any such bias. Monotremes are an ancient lineage of egg-laying mammals with distinct biological properties, which include unique germline features. Here, we sought to assess the presence and potential characteristics of male mutation bias in platypus and the short-beaked echidna based on substitution rate analyses of X, Y, and autosomes. We established the presence of moderate male mutation bias in monotremes, corresponding to an α value of 2.12–3.69. Given that it has been unclear what proportion of the variation in substitution rates on the different chromosomal classes is really due to differential number of replications, we analyzed the influence of other confounding forces (selection, replication-timing, etc.) and found that male mutation bias is the main force explaining the between-chromosome classes differences in substitution rates. Finally, we estimated the proportion of variation at the gene level in substitution rates that is owing to replication effects and found that this phenomenon can explain >68% of these variations in monotremes, and in control species, rodents, and primates.

Rat astrocytes are more supportive for mouse OPC self-renewal than mouse astrocytes in culture

Mouse primary oligodendrocyte precursor cells (OPCs) are increasingly used to study the molecular mechanisms underlying the phenotype changes in oligodendrocyte differentiation and axonal myelination observed in transgenic or mutant mouse models. However, mouse OPCs are much more difficult to be isolated by the simple dissociation culture of brain tissues than their rat counterparts. To date, the mechanisms underlying the species difference in OPC preparation remain obscure. In this study, we showed that astrocytes from rats have a stronger effect than those from mouse in promoting OPC proliferation and survival in vitro. Mouse astrocytes displayed significantly weaker viability in culture and reduced potential in maintaining OPC self-renewal, as confirmed by culturing OPCs with conditioned media from rat or mouse astrocytes. These results explained the reason for why stratified cultures of OPCs and astrocytes are difficult to be achieved in mouse CNS tissues. Based on these findings, we adopted inactivated rat astrocytes as feeder cells to support the self-renewal of mouse cortical OPCs and preparation of high-purity mouse OPCs. © 2016 Wiley Periodicals, Inc. Develop Neurobiol 77: 907-916, 2017.

Genomes of Ellobius species provide insight into the evolutionary dynamics of mammalian sex chromosomes

The X and Y sex chromosomes of placental mammals show hallmarks of a tumultuous evolutionary past. The X Chromosome has a rich and conserved gene content, while the Y Chromosome has lost most of its genes. In the Transcaucasian mole vole Ellobius lutescens, the Y Chromosome including Sry has been lost, and both females and males have a 17,X diploid karyotype. Similarly, the closely related Ellobius talpinus, has a 54,XX karyotype in both females and males. Here, we report the sequencing and assembly of the E. lutescens and E. talpinus genomes. The results indicate that the loss of the Y Chromosome in E. lutescens and E. talpinus occurred in two independent events. Four functional homologs of mouse Y-Chromosomal genes were detected in both female and male E. lutescens, of which three were also detected in the E. talpinus genome. One of these is Eif2s3y, known as the only Y-derived gene that is crucial for successful male meiosis. Female and male E. lutescens can carry one and the same X Chromosome with a largely conserved gene content, including all genes known to function in X Chromosome inactivation. The availability of the genomes of these mole vole species provides unique models to study the dynamics of sex chromosome evolution.

Phylogenetic relationships and protein modelling revealed two distinct subfamilies of group II HKT genes between crop and model grasses

Molecular evolution of large protein families in closely related species can provide useful insights on structural functional relationships. Phylogenetic analysis of the grass-specific group II HKT genes identified two distinct subfamilies, I and II. Subfamily II was represented in all species, whereas subfamily I was identified only in the small grain cereals and possibly originated from an ancestral gene duplication post divergence from the coarse grain cereal lineage. The core protein structures were highly analogous despite there being no more than 58% amino acid identity between members of the two subfamilies. Distinctly variable regions in known functional domains, however, indicated functional divergence of the two subfamilies. The subsets of codons residing external to known functional domains predicted signatures of positive Darwinian selection potentially identifying new domains of functional divergence and providing new insights on the structural function and relationships between protein members of the two subfamilies.

The utility of Apc-mutant rats in modeling human colon cancer

Prior to the advent of genetic engineering in the mouse, the rat was the model of choice for investigating the etiology of cancer. Now, recent advances in the manipulation of the rat genome, combined with a growing recognition of the physiological differences between mice and rats, have reignited interest in the rat as a model of human cancer. Two recently developed rat models, the polyposis in the rat colon (Pirc) and Kyoto Apc Delta (KAD) strains, each carry mutations in the intestinal-cancer-associated adenomatous polyposis coli (Apc) gene. In contrast to mouse models carrying Apc mutations, in which cancers develop mainly in the small intestine rather than in the colon and there is no gender bias, these rat models exhibit colonic predisposition and gender-specific susceptibility, as seen in human colon cancer. The rat also provides other experimental resources as a model organism that are not provided by the mouse: the structure of its chromosomes facilitates the analysis of genomic events, the size of its colon permits longitudinal analysis of tumor growth, and the size of biological samples from the animal facilitates multiplexed molecular analyses of the tumor and its host. Thus, the underlying biology and experimental resources of these rat models provide important avenues for investigation. We anticipate that advances in disease modeling in the rat will synergize with resources that are being developed in the mouse to provide a deeper understanding of human colon cancer.

Mitochondrial-nuclear interactions: compensatory evolution or variable functional constraint among vertebrate oxidative phosphorylation genes?

Oxidative phosphorylation (OXPHOS), the major energy-producing pathway in aerobic organisms, includes protein subunits encoded by both mitochondrial (mt) and nuclear (nu) genomes. How these independent genomes have coevolved is a long-standing question in evolutionary biology. Although mt genes evolve faster than most nu genes, maintenance of OXPHOS structural stability and functional efficiency may involve correlated evolution of mt and nu OXPHOS genes. The nu OXPHOS genes might be predicted to exhibit accelerated evolutionary rates to accommodate the elevated substitution rates of mt OXPHOS subunits with which they interact. Evolutionary rates of nu OXPHOS genes should, therefore, be higher than that of nu genes that are not involved in OXPHOS (nu non-OXPHOS). We tested the compensatory evolution hypothesis by comparing the evolutionary rates (synonymous substitution rate dS and nonsynonymous substitution rate dN) among 13 mt OXPHOS genes, 60 nu OXPHOS genes, and 77 nu non-OXPHOS genes in vertebrates (7 fish and 40 mammal species). The results from a combined analysis of all OXPHOS subunits fit the predictions of the hypothesis. However, results from two OXPHOS complexes did not fit this pattern when analyzed separately. We found that the d(N) of nu OXPHOS genes for "core" subunits (those involved in the major catalytic activity) was lower than that of "noncore" subunits, whereas there was no significant difference in d(N) between genes for nu non-OXPHOS and core subunits. This latter finding suggests that compensatory changes play a minor role in the evolution of OXPHOS genes and that the observed accelerated nu substitution rates are due largely to reduced functional constraint on noncore subunits.

Engineering Xenopus embryos for phenotypic drug discovery screening

Many rare human inherited diseases remain untreatable despite the fact that the disease causing genes are known and adequate mouse disease models have been developed. In vivo phenotypic drug screening relies on isolating drug candidates by their ability to produce a desired therapeutic phenotype in whole organisms. Embryos of zebrafish and Xenopus frogs are abundant, small and free-living. They can be easily arrayed in multi-well dishes and treated with small organic molecules. With the development of novel genome modification tools, such a zinc-finger nucleases (ZFNs), transcription activator-like effector nucleases (TALENs), and CRISPR/Cas, it is now possible to efficiently engineer non-mammalian models of inherited human diseases. Here, we will review the rapid progress made in adapting these novel genome editing tools to Xenopus. The advantages of Xenopus embryos as in vivo models to study human inherited diseases will be presented and their utility for drug discovery screening will be discussed. Being a tetrapod, Xenopus complements zebrafish as an indispensable non-mammalian animal model for the study of human disease pathologies and the discovery of novel therapeutics for inherited diseases.

Directional Darwinian Selection in proteins

Background

Molecular evolution is a very active field of research, with several complementary approaches, including dN/dS, HON90, MM01, and others. Each has documented strengths and weaknesses, and no one approach provides a clear picture of how natural selection works at the molecular level. The purpose of this work is to present a simple new method that uses quantitative amino acid properties to identify and characterize directional selection in proteins.

Methods

Inferred amino acid replacements are viewed through the prism of a single physicochemical property to determine the amount and direction of change caused by each replacement. This allows the calculation of the probability that the mean change in the single property associated with the amino acid replacements is equal to zero (H₀: μ = 0; i.e., no net change) using a simple two-tailed t-test.

Results

Example data from calanoid and cyclopoid copepod cytochrome oxidase subunit I sequence pairs are presented to demonstrate how directional selection may be linked to major shifts in adaptive zones, and that convergent evolution at the whole organism level may be the result of convergent protein adaptations.

Conclusions

Rather than replace previous methods, this new method further complements existing methods to provide a holistic glimpse of how natural selection shapes protein structure and function over evolutionary time.

Rate heterogeneity in six protein-coding genes from the holoparasite Balanophora (Balanophoraceae) and other taxa of Santalales

BACKGROUND AND AIMS

The holoparasitic flowering plant Balanophora displays extreme floral reduction and was previously found to have enormous rate acceleration in the nuclear 18S rDNA region. So far, it remains unclear whether non-ribosomal, protein-coding genes of Balanophora also evolve in an accelerated fashion and whether the genes with high substitution rates retain their functionality. To tackle these issues, six different genes were sequenced from two Balanophora species and their rate variation and expression patterns were examined.

METHODS

Sequences including nuclear PI, euAP3, TM6, LFY and RPB2 and mitochondrial matR were determined from two Balanophora spp. and compared with selected hemiparasitic species of Santalales and autotrophic core eudicots. Gene expression was detected for the six protein-coding genes and the expression patterns of the three B-class genes (PI, AP3 and TM6) were further examined across different organs of B. laxiflora using RT-PCR.

KEY RESULTS

Balanophora mitochondrial matR is highly accelerated in both nonsynonymous (d(N)) and synonymous (d(S)) substitution rates, whereas the rate variation of nuclear genes LFY, PI, euAP3, TM6 and RPB2 are less dramatic. Significant d(S) increases were detected in Balanophora PI, TM6, RPB2 and d(N) accelerations in euAP3. All of the protein-coding genes are expressed in inflorescences, indicative of their functionality. PI is restrictively expressed in tepals, synandria and floral bracts, whereas AP3 and TM6 are widely expressed in both male and female inflorescences.

CONCLUSIONS

Despite the observation that rates of sequence evolution are generally higher in Balanophora than in hemiparasitic species of Santalales and autotrophic core eudicots, the five nuclear protein-coding genes are functional and are evolving at a much slower rate than 18S rDNA. The mechanism or mechanisms responsible for rapid sequence evolution and concomitant rate acceleration for 18S rDNA and matR are currently not well understood and require further study in Balanophora and other holoparasites.

BlastR--fast and accurate database searches for non-coding RNAs

We present and validate BlastR, a method for efficiently and accurately searching non-coding RNAs. Our approach relies on the comparison of di-nucleotides using BlosumR, a new log-odd substitution matrix. In order to use BlosumR for comparison, we recoded RNA sequences into protein-like sequences. We then showed that BlosumR can be used along with the BlastP algorithm in order to search non-coding RNA sequences. Using Rfam as a gold standard, we benchmarked this approach and show BlastR to be more sensitive than BlastN. We also show that BlastR is both faster and more sensitive than BlastP used with a single nucleotide log-odd substitution matrix. BlastR, when used in combination with WU-BlastP, is about 5% more accurate than WU-BlastN and about 50 times slower. The approach shown here is equally effective when combined with the NCBI-Blast package. The software is an open source freeware available from www.tcoffee.org/blastr.html.

The positive correlation between dN/dS and dS in mammals is due to runs of adjacent substitutions

A positive correlation between ω, the ratio of the nonsynonymous and synonymous substitution rates, and dS, the synonymous substitution rate has recently been reported. This correlation is unexpected under simple evolutionary models. Here, we investigate two explanations for this correlation: first, whether it is a consequence of a statistical bias in the estimation of ω and second, whether it is due to substitutions at adjacent sites. Using simulations, we show that estimates of ω are biased when levels of divergence are low. This is true using the methods of Yang and Nielsen, Nei and Gojobori, and Muse and Gaut. Although the bias could generate a positive correlation between ω and dS, we show that it is unlikely to be the main determinant. Instead we show that the correlation is reduced when genes that are high quality in sequence, annotation, and alignment are used. The remaining--likely genuine--positive correlation appears to be due to adjacent tandem substitutions; single substitutions, though far more numerous, do not contribute to the correlation. Genuine adjacent substitutions may be due to mutation or selection.

Modeling compositional dynamics based on GC and purine contents of protein-coding sequences

Background

Understanding the compositional dynamics of genomes and their coding sequences is of great significance in gaining clues into molecular evolution and a large number of publically-available genome sequences have allowed us to quantitatively predict deviations of empirical data from their theoretical counterparts. However, the quantification of theoretical compositional variations for a wide diversity of genomes remains a major challenge.

Results

To model the compositional dynamics of protein-coding sequences, we propose two simple models that take into account both mutation and selection effects, which act differently at the three codon positions, and use both GC and purine contents as compositional parameters. The two models concern the theoretical composition of nucleotides, codons, and amino acids, with no prerequisite of homologous sequences or their alignments. We evaluated the two models by quantifying theoretical compositions of a large collection of protein-coding sequences (including 46 of Archaea, 686 of Bacteria, and 826 of Eukarya), yielding consistent theoretical compositions across all the collected sequences.

Conclusions

We show that the compositions of nucleotides, codons, and amino acids are largely determined by both GC and purine contents and suggest that deviations of the observed from the expected compositions may reflect compositional signatures that arise from a complex interplay between mutation and selection via DNA replication and repair mechanisms.

Reviewers

This article was reviewed by Zhaolei Zhang (nominated by Mark Gerstein), Guruprasad Ananda (nominated by Kateryna Makova), and Daniel Haft.

Translational selection on SHH genes

Codon usage bias has been observed in various organisms. In this study, the correlation between SHH genes expression in some tissues and codon usage features was analyzed by bioinformatics. We found that translational selection may act on compositional features of this set of genes.

Measuring and detecting molecular adaptation in codon usage against nonsense errors during protein translation

Codon usage bias (CUB) has been documented across a wide range of taxa and is the subject of numerous studies. While most explanations of CUB invoke some type of natural selection, most measures of CUB adaptation are heuristically defined. In contrast, we present a novel and mechanistic method for defining and contextualizing CUB adaptation to reduce the cost of nonsense errors during protein translation. Using a model of protein translation, we develop a general approach for measuring the protein production cost in the face of nonsense errors of a given allele as well as the mean and variance of these costs across its coding synonyms. We then use these results to define the nonsense error adaptation index (NAI) of the allele or a contiguous subset thereof. Conceptually, the NAI value of an allele is a relative measure of its elevation on a specific and well-defined adaptive landscape. To illustrate its utility, we calculate NAI values for the entire coding sequence and across a set of nonoverlapping windows for each gene in the Saccharomyces cerevisiae S288c genome. Our results provide clear evidence of adaptation to reduce the cost of nonsense errors and increasing adaptation with codon position and expression. The magnitude and nature of this adaptation are also largely consistent with simulation results in which nonsense errors are the only selective force driving CUB evolution. Because NAI is derived from mechanistic models, it is both easier to interpret and more amenable to future refinement than other commonly used measures of codon bias. Further, our approach can also be used as a starting point for developing other mechanistically derived measures of adaptation such as for translational accuracy.

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.

The genomic distribution of intraspecific and interspecific sequence divergence of human segmental duplications relative to human/chimpanzee chromosomal rearrangements

Background

It has been suggested that chromosomal rearrangements harbor the molecular footprint of the biological phenomena which they induce, in the form, for instance, of changes in the sequence divergence rates of linked genes. So far, all the studies of these potential associations have focused on the relationship between structural changes and the rates of evolution of single-copy DNA and have tried to exclude segmental duplications (SDs). This is paradoxical, since SDs are one of the primary forces driving the evolution of structure and function in our genomes and have been linked not only with novel genes acquiring new functions, but also with overall higher DNA sequence divergence and major chromosomal rearrangements.

Results

Here we take the opposite view and focus on SDs. We analyze several of the features of SDs, including the rates of intraspecific divergence between paralogous copies of human SDs and of interspecific divergence between human SDs and chimpanzee DNA. We study how divergence measures relate to chromosomal rearrangements, while considering other factors that affect evolutionary rates in single copy DNA.

Conclusion

We find that interspecific SD divergence behaves similarly to divergence of single-copy DNA. In contrast, old and recent paralogous copies of SDs do present different patterns of intraspecific divergence. Also, we show that some relatively recent SDs accumulate in regions that carry inversions in sister lineages.

Weak correlation between sequence conservation in promoter regions and in protein-coding regions of human-mouse orthologous gene pairs

Background

Interspecies sequence comparison is a powerful tool to extract functional or evolutionary information from the genomes of organisms. A number of studies have compared protein sequences or promoter sequences between mammals, which provided many insights into genomics. However, the correlation between protein conservation and promoter conservation remains controversial.

Results

We examined promoter conservation as well as protein conservation for 6,901 human and mouse orthologous genes, and observed a very weak correlation between them. We further investigated their relationship by decomposing it based on functional categories, and identified categories with significant tendencies. Remarkably, the 'ribosome' category showed significantly low promoter conservation, despite its high protein conservation, and the 'extracellular matrix' category showed significantly high promoter conservation, in spite of its low protein conservation.

Conclusion

Our results show the relation of gene function to protein conservation and promoter conservation, and revealed that there seem to be nonparallel components between protein and promoter sequence evolution.

The pattern of evolution of smaller-scale gene duplicates in mammalian genomes is more consistent with neo- than subfunctionalisation

Gene duplication and the accompanying release of negative selective pressure on the duplicate pair is thought to be the key process that makes functional change in the coding and regulatory regions of genomes possible. However, the nature of these changes remains unresolved. There are a number of models for the fate of gene duplicates, the two most prominent of which are neofunctionalisation and subfunctionalisation, but it is still unclear which is the dominant fate. Using a dataset consisting of smaller-scale (tandem and segmental) duplications identified from the genomes of four fully sequenced mammalian genomes, we characterise two key features of smaller-scale duplicate evolution: the rate of pseudogenisation and the rate of accumulation of replacement substitutions in the coding sequence. We show that the best fitting model for gene duplicate survival is a Weibull function with a downward sloping convex hazard function which implies that the rate of pseudogenisation of a gene declines rapidly with time since duplication. Our analysis of the accumulation of replacement substitutions per replacement site shows that they accumulate on average at 64% of the neutral expectation immediately following duplication and as high as 73% in the human lineage. Although this rate declines with time since duplication, it takes several tens of millions of years before it has declined to half its initial value. We show that the properties of the gene death rate and of the accumulation of replacement substitutions are more consistent with neofunctionalisation (or subfunctionalisation followed by neofunctionalisation) than they are with subfunctionalisation alone or any of the other alternative modes of evolution of smaller-scale duplicates.

Sequence variation in the Tbx4 gene in marine mammals

The amino-acid sequences of the T-domain region of the Tbx4 gene, which is required for hindlimb development, are 100% identical in humans and mice. Cetaceans have lost most of their hindlimb structure, although hindlimb buds are present in very early cetacean embryos. To examine whether the Tbx4 gene has the same function in cetaceans as in other mammals, we analyzed Tbx4 sequences from cetaceans, dugong, artiodactyls and marine carnivores. A total of 39 primers were designed using human and dog Tbx4 nucleotide sequences. Exons 3, 4, 5, 6, 7, and 8 of the Tbx4 genes from cetaceans, artiodactyls, and marine carnivores were sequenced. Non-synonymous substitution sites were detected in the T-domain regions from some cetacean species, but were not detected in those from artiodactyls, the dugong, or the carnivores. The C-terminal regions contained a number of non-synonymous substitutions. Although some indels were present, they were in groups of three nucleotides and therefore did not cause frame shifts. The dN/dS values for the T-domain and C-terminal regions of the cetacean and artiodactylous Tbx4 genes were much lower than 1, indicating that the Tbx4 gene maintains it function in cetaceans, although full expression leading to hindlimb development is suppressed.

A novel olfactory receptor gene family in teleost fish

While for two of three mammalian olfactory receptor families (OR and V2R) ortholog teleost families have been identified, the third family (V1R) has been thought to be represented by a single, closely linked gene pair. We identified four further V1R-like genes in every teleost species analyzed (Danio rerio, Gasterosteus aculeatus, Oryzias latipes, Tetraodon nigroviridis, Takifugu rubripes). In the phylogenetic analysis these ora genes (olfactory receptor class A-related) form a single clade, which includes the entire mammalian V1R superfamily. Homologies are much lower in paralogs than in orthologs, indicating that all six family members are evolutionarily much older than the speciation events in the teleost lineage analyzed here. These ora genes are under strong negative selection, as evidenced by very small d(N)/d(S) values in comparisons between orthologs. A pairwise configuration in the phylogenetic tree suggests the existence of three ancestral Ora subclades, one of which has been lost in amphibia, and a further one in mammals. Unexpectedly, two ora genes exhibit a highly conserved multi-exonic structure and four ora genes are organized in closely linked gene pairs across all fish species studied. All ora genes are expressed specifically in the olfactory epithelium of zebrafish, in sparse cells within the sensory surface, consistent with the expectation for olfactory receptors. The ora gene repertoire is highly conserved across teleosts, in striking contrast to the frequent species-specific expansions observed in tetrapod, especially mammalian V1Rs, possibly reflecting a major shift in gene regulation as well as gene function upon the transition to tetrapods.

Human oncogene tissue-specific expression level significantly correlates with sequence compositional features

To explore whether there exist correlations between human gene expressions and corresponding sequence features, the expression levels for 81 oncogenes in 24 human tissues were collected and correlated with 159 sequence features. It was found that there do exist significant correlations between them, some of which are of significance to understanding translational selection on sequence features of human genes and some have important implications for diagnosing cancers.

Developmental stage related patterns of codon usage and genomic GC content: searching for evolutionary fingerprints with models of stem cell differentiation

Developmental-stage-related patterns of gene expression correlate with codon usage and genomic GC content in stem cell hierarchies.

Background

The usage of synonymous codons shows considerable variation among mammalian genes. How and why this usage is non-random are fundamental biological questions and remain controversial. It is also important to explore whether mammalian genes that are selectively expressed at different developmental stages bear different molecular features.

Results

In two models of mouse stem cell differentiation, we established correlations between codon usage and the patterns of gene expression. We found that the optimal codons exhibited variation (AT- or GC-ending codons) in different cell types within the developmental hierarchy. We also found that genes that were enriched (developmental-pivotal genes) or specifically expressed (developmental-specific genes) at different developmental stages had different patterns of codon usage and local genomic GC (GCg) content. Moreover, at the same developmental stage, developmental-specific genes generally used more GC-ending codons and had higher GCg content compared with developmental-pivotal genes. Further analyses suggest that the model of translational selection might be consistent with the developmental stage-related patterns of codon usage, especially for the AT-ending optimal codons. In addition, our data show that after human-mouse divergence, the influence of selective constraints is still detectable.

Conclusion

Our findings suggest that developmental stage-related patterns of gene expression are correlated with codon usage (GC3) and GCg content in stem cell hierarchies. Moreover, this paper provides evidence for the influence of natural selection at synonymous sites in the mouse genome and novel clues for linking the molecular features of genes to their patterns of expression during mammalian ontogenesis.

Widespread positive selection in synonymous sites of mammalian genes

Evolution of protein sequences is largely governed by purifying selection, with a small fraction of proteins evolving under positive selection. The evolution at synonymous positions in protein-coding genes is not nearly as well understood, with the extent and types of selection remaining, largely, unclear. A statistical test to identify purifying and positive selection at synonymous sites in protein-coding genes was developed. The method compares the rate of evolution at synonymous sites (Ks) to that in intron sequences of the same gene after sampling the aligned intron sequences to mimic the statistical properties of coding sequences. We detected purifying selection at synonymous sites in approximately 28% of the 1,562 analyzed orthologous genes from mouse and rat, and positive selection in approximately 12% of the genes. Thus, the fraction of genes with readily detectable positive selection at synonymous sites is much greater than the fraction of genes with comparable positive selection at nonsynonymous sites, i.e., at the level of the protein sequence. Unlike other genes, the genes with positive selection at synonymous sites showed no correlation between Ks and the rate of evolution in nonsynonymous sites (Ka), indicating that evolution of synonymous sites under positive selection is decoupled from protein evolution. The genes with purifying selection at synonymous sites showed significant anticorrelation between Ks and expression level and breadth, indicating that highly expressed genes evolve slowly. The genes with positive selection at synonymous sites showed the opposite trend, i.e., highly expressed genes had, on average, higher Ks. For the genes with positive selection at synonymous sites, a significantly lower mRNA stability is predicted compared to the genes with negative selection. Thus, mRNA destabilization could be an important factor driving positive selection in nonsynonymous sites, probably, through regulation of expression at the level of mRNA degradation and, possibly, also translation rate. So, unexpectedly, we found that positive selection at synonymous sites of mammalian genes is substantially more common than positive selection at the level of protein sequences. Positive selection at synonymous sites might act through mRNA destabilization affecting mRNA levels and translation.

On the association between chromosomal rearrangements and genic evolution in humans and chimpanzees

BACKGROUND

The role that chromosomal rearrangements might have played in the speciation processes that have separated the lineages of humans and chimpanzees has recently come into the spotlight. To date, however, results are contradictory. Here we revisit this issue by making use of the available human and chimpanzee genome sequence to study the relationship between chromosomal rearrangements and rates of DNA sequence evolution.

RESULTS

Contrary to previous findings for this pair of species, we show that genes located in the rearranged chromosomes that differentiate the genomes of humans and chimpanzees, especially genes within rearrangements themselves, present lower divergence than genes elsewhere in the genome. Still, there are considerable differences between individual chromosomes. Chromosome 4, in particular, presents higher divergence in genes located within its rearrangement.

CONCLUSION

A first conclusion of our analysis is that divergence is lower for genes located in rearranged chromosomes than for those in colinear chromosomes. We also report that non-coding regions within rearranged regions tend to have lower divergence than non-coding regions outside them. These results suggest an association between chromosomal rearrangements and lower non-coding divergence that has not been reported before, even if some chromosomes do not follow this trend and could be potentially associated with a speciation episode. In summary, without excluding it, our results suggest that chromosomal speciation has not been common along the human and chimpanzee lineage.

Isolation of the feline alpha1,3-galactosyltransferase gene, expression in transfected human cells and its phylogenetic analysis

The enzyme alpha 1,3-galactosyltransferase (alpha1,3-GT), which catalyzes synthesis of terminal alpha-galactosyl epitopes (Gal alpha1,3Gal beta1-4GlcNAc-R), is produced in non-primate mammals, prosimians and new-world monkeys, but not in old-world monkeys, apes and humans. We cloned and sequenced a cDNA that contains the coding sequence of the feline alpha1,3-GT gene. Flow cytometric analysis demonstrated that the alpha-galactosyl epitope was expressed on the surface of a human cell line transduced with an expression vector containing this cDNA, and this alpha-galactosyl epitope expression subsided by alpha-galactosidase treatment. The open reading frame of the feline alpha1,3-GT cDNA is 1,113 base pairs in length and encodes 371 amino acids. The nucleotide sequence and its deduced amino acid sequence of the feline alpha1,3-GT gene are 88-90% and 85-87%, respectively, similar to the reported sequences of the bovine, porcine, marmoset and cebus monkey alpha1,3-GT genes, while they are 88% and 82-83%, respectively, similar to those of the orangutan and human alpha1,3-GT pseudogenes, and 81% and 77%, respectively, similar to the murine alpha1,3-GT gene. Thus, the alpha1,3-GT genes and pseudogenes of mammals are highly similar. Ratios of non-synonymous nucleotide changes among the primate pseudogenes as well as the primate genes are still higher than the ratios of non-primates, suggesting that the primate alpha1,3-GT genes tend to be divergent.

Accelerated evolutionary rate may be responsible for the emergence of lineage-specific genes in ascomycota

The evolutionary origin of "orphan" genes, genes that lack sequence similarity to any known gene, remains a mystery. One suggestion has been that most orphan genes evolve rapidly so that similarity to other genes cannot be traced after a certain evolutionary distance. This can be tested by examining the divergence rates of genes with different degrees of lineage specificity. Here the lineage specificity (LS) of a gene describes the phylogenetic distribution of that gene's orthologues in related species. Highly lineage-specific genes will be distributed in fewer species in a phylogeny. In this study, we have used the complete genomes of seven ascomycotan fungi and two animals to define several levels of LS, such as Eukaryotes-core, Ascomycota-core, Euascomycetes-specific, Hemiascomycetes-specific, Aspergillus-specific, and Saccharomyces-specific. We compare the rates of gene evolution in groups of higher LS to those in groups with lower LS. Molecular evolutionary analyses indicate an increase in nonsynonymous nucleotide substitution rates in genes with higher LS. Several analyses suggest that LS is correlated with the evolutionary rate of the gene. This correlation is stronger than those of a number of other factors that have been proposed as predictors of a gene's evolutionary rate, including the expression level of genes, gene essentiality or dispensability, and the number of protein-protein interactions. The accelerated evolutionary rates of genes with higher LS may reflect the influence of selection and adaptive divergence during the emergence of orphan genes. These analyses suggest that accelerated rates of gene evolution may be responsible for the emergence of apparently orphan genes.

High guanine and cytosine content increases mRNA levels in mammalian cells

Mammalian genes are highly heterogeneous with respect to their nucleotide composition, but the functional consequences of this heterogeneity are not clear. In the previous studies, weak positive or negative correlations have been found between the silent-site guanine and cytosine (GC) content and expression of mammalian genes. However, previous studies disregarded differences in the genomic context of genes, which could potentially obscure any correlation between GC content and expression. In the present work, we directly compared the expression of GC-rich and GC-poor genes placed in the context of identical promoters and UTR sequences. We performed transient and stable transfections of mammalian cells with GC-rich and GC-poor versions of Hsp70, green fluorescent protein, and IL2 genes. The GC-rich genes were expressed several-fold to over a 100-fold more efficiently than their GC-poor counterparts. This effect was not due to different translation rates of GC-rich and GC-poor mRNA. On the contrary, the efficient expression of GC-rich genes resulted from their increased steady-state mRNA levels. mRNA degradation rates were not correlated with GC content, suggesting that efficient transcription or mRNA processing is responsible for the high expression of GC-rich genes. We conclude that silent-site GC content correlates with gene expression efficiency in mammalian cells.

The effect of nucleotide composition on gene transcription is investigated for Hsp70, GFP, and IL-2, which all show increased expression to correlate with increased GC content at codon position 3.

Large-scale analyses of synonymous substitution rates can be sensitive to assumptions about the process of mutation

A popular approach to examine the roles of mutation and selection in the evolution of genomes has been to consider the relationship between codon bias and synonymous rates of molecular evolution. A significant relationship between these two quantities is taken to indicate the action of weak selection on substitutions among synonymous codons. The neutral theory predicts that the rate of evolution is inversely related to the level of functional constraint. Therefore, selection against the use of non-preferred codons among those coding for the same amino acid should result in lower rates of synonymous substitution as compared with sites not subject to such selection pressures. However, reliably measuring the extent of such a relationship is problematic, as estimates of synonymous rates are sensitive to our assumptions about the process of molecular evolution. Previous studies showed the importance of accounting for unequal codon frequencies, in particular when synonymous codon usage is highly biased. Yet, unequal codon frequencies can be modeled in different ways, making different assumptions about the mutation process. Here we conduct a simulation study to evaluate two different ways of modeling uneven codon frequencies and show that both model parameterizations can have a dramatic impact on rate estimates and affect biological conclusions about genome evolution. We reanalyze three large data sets to demonstrate the relevance of our results to empirical data analysis.

Enhanced adaptive evolution of sperm-expressed genes on the mammalian X chromosome

Genes on the mammalian X chromosome may be under unique evolutionary pressure due to their hemizygous expression in males. Since any recessive deleterious mutation would immediately be expressed in males and, therefore, efficiently removed from the population, selective constraint could be more pronounced in X-linked genes. Conversely, if a recessive mutation were beneficial, its immediate exposure to selection would be advantageous, and would facilitate adaptive evolution. We tested for positive selection in a total of 86 genes using a maximum likelihood approach, including 40 sperm-expressed and 46 non-sperm, tissue-specific genes. We find evidence to suggest that X-linkage enhances the effects of positive selection in sperm-expressed genes in terms of the number of codons affected, and report a general trend for positively selected genes to reside on the X chromosome rather than on the autosomes. Our data suggest that hemizygous expression in males makes the X chromosome a preferred location for positively selected sperm genes that do not require postmeiotic transcription.

Further understanding human disease genes by comparing with housekeeping genes and other genes

Background

Several studies have compared various features of heritable disease genes with other so called non-disease genes, but they have yielded some conflicting results. A potential problem in those studies is that the non-disease genes contained a large number of essential genes – genes which are indispensable for humans to survive and reproduce. Since a functional disruption of an essential gene has fatal consequences, it's more reasonable to regard essential genes as extremely severe "disease" genes. Here we perform a comparative study on the features of human essential, disease, and other genes.

Results

In the absence of a set of well defined human essential genes, we consider a set of 1,789 ubiquitously expressed human genes (UEHGs), also known as housekeeping genes, as an approximation. We demonstrate that UEHGs are very likely to contain a large proportion of essential genes. We show that the UEHGs, disease genes and other genes are different in their evolutionary conservation rates, DNA coding lengths, gene functions, etc. Our findings systematically confirm that disease genes have an intermediate essentiality which is less than housekeeping genes but greater than other human genes.

Conclusion

The human genome may contain thousands of essential genes having features which differ significantly from disease and other genes. We propose to classify them as a unique group for comparisons of disease genes with non-disease genes. This new way of classification and comparison enables us to have a clearer understanding of disease genes.

Identification of a tandem duplicated array in the Rhox alpha locus on mouse chromosome X

The Rhox family of homeobox-containing genes maps to three gene dense clusters, alpha, beta, and gamma, on the mouse X chromosome. Rhox genes are expressed primarily in reproductive and extra-embryonic tissues and may be regulated by colinear mechanisms. Uniquely among family members, Rhox4 is also expressed in the developing and adult thymus, where its striking restriction to the third pharyngeal pouch has suggested a role in lineage specification. In this article we describe a previously unreported duplication within the Rhox alpha locus that contains seven copies of Rhox4 and eight copies of Rhox2 and Rhox3 in a tandem array. We further show that all seven Rhox4 copies are expressed, although preferential expression occurs and differs between tissues. In contrast to reproductive tissues, we found no evidence of colinear expression of the Rhox alpha cluster during thymus development. All subspecies of mice examined contained multiple copies of Rhox4. However, the only predicted ortholog of Rhox4, rat Rhox4, is present in a single copy, suggesting that the duplications arose at the time of rat and mouse lineage divergence. Finally, no changes in Rhox4 expression were detected in mice with defects in thymus organogenesis, placing it upstream or outside of established transcriptional pathways.

Evolutionary rates and expression level in Chlamydomonas

In many biological systems, especially bacteria and unicellular eukaryotes, rates of synonymous and nonsynonymous nucleotide divergence are negatively correlated with the level of gene expression, a phenomenon that has been attributed to natural selection. Surprisingly, this relationship has not been examined in many important groups, including the unicellular model organism Chlamydomonas reinhardtii. Prior to this study, comparative data on protein-coding sequences from C. reinhardtii and its close noninterfertile relative C. incerta were very limited. We compiled and analyzed protein-coding sequences for 67 nuclear genes from these taxa; the sequences were mostly obtained from the C. reinhardtii EST database and our C. incerta EST data. Compositional and synonymous codon usage biases varied among genes within each species but were highly correlated between the orthologous genes of the two species. Relative rates of synonymous and nonsynonymous substitution across genes varied widely and showed a strong negative correlation with the level of gene expression estimated by the codon adaptation index. Our comparative analysis of substitution rates in introns of lowly and highly expressed genes suggests that natural selection has a larger contribution than mutation to the observed correlation between evolutionary rates and gene expression level in Chlamydomonas.

Chromosomal rearrangements are associated with higher rates of molecular evolution in mammals

Evolutionary rates are not uniformly distributed across the genome. Knowledge about the biological causes of this observation is still incomplete, but its exploration has provided valuable insight into the genomical, historical and demographical variables that influence rates of genetic divergence. Recent studies suggest a possible association between chromosomal rearrangements and regions of greater divergence, but evidence is limited and contradictory. Here, we test the hypothesis of a relationship between chromosomal rearrangements and higher rates of molecular evolution by studying the genomic distribution of divergence between 12,000 human-mouse orthologous genes. Our results clearly show that genes located in genomic regions that have been highly rearranged between the two species present higher rates of synonymous (0.7686 vs. 0.7076) and non-synonymous substitution (0.1014 vs. 0.0871), and that synonymous substitution rates are higher in genes close to the breakpoints of individual rearrangements. The many potential causes of such striking are discussed, particularly in the light of speciation models suggesting that chromosomal rearrangements may have contributed to some of the speciation processes along the human and mouse lineages. Still, there are other possible causes and further research is needed to properly explore them.

Enhanced Synonymous Site Divergence in Positively Selected VertebrateAntimicrobial Peptide Genes

Nonrandom patterns associated with adaptively evolving genes can shed light on how selection and mutation produce rapid changes in sequences. I examine such patterns in two independent families of antimicrobial peptide genes: those in frogs, which are known to have evolved under positive selection, and those in flatfishes, which I show have also evolved under positive selection. I address two recently proposed hypotheses about the molecular evolution of antimicrobial peptide genes. The first is that the mature peptide region is replicated by an error-prone polymerase that increases the mutation rate and the transversion/transition ratio compared to the signal sequence of the same genes. The second is that mature peptides evolve in a coordinated fashion with their propieces, such that a change in net charge in one molecular region prompts an opposite change in charge in the other region. I test these hypotheses using alternative methods that minimize alignment errors, correct for phylogenetic nonindependence, reduce sequence saturation, and account for differing selection pressures on different regions of the gene. In both gene families I show that divergence at both synonymous and nonsynonymous sites within the mature peptide region is enhanced. However, in neither gene family is there evidence of an increased mutational transversion/transition ratio or coordinated evolution. My observations are consistent with either an elevated mutation rate in an adaptively evolving gene region or widespread selection on "silent" sites. These hypotheses challenge the assumption that mutations are random and can be measured by the synonymous substitution rate.

Ancient divergence in bathypelagic lake tanganyika deepwater cichlids: mitochondrial phylogeny of the tribe bathybatini

The cichlid species flock of Lake Tanganyika represents a polyphyletic assemblage of eight ancestral lineages, which colonized the emerging lake independently. Our study is focused on one of these lineages, the Bathybatini, a tribe of specialized piscivorous cichlids of the deep pelagic zone. By analyzing three mtDNA gene segments of all eight species of the tribe and two species of the closely related Trematocarini, we propose on the basis of a linearized tree analysis that the Bathybatini comprise two distinct lineages, the genera Hemibates and Bathybates, that seeded the primary lacustrine Tanganyika radiation independently. The genus Hemibates is likely to represent a distinct lineage that emerged simultaneously with the tribe Trematocarini and the genus Bathybates and should be therefore treated as a distinct tribe. Within the genus Bathybates, B. minor clearly represents the most ancestral split and is likely to have diverged from the remaining species in the course of the "primary lacustrine Tanganyika radiation" during which also the radiations of the Lamprologini and the H-lineage took place. The remaining "large" Bathybates species also diversified almost simultaneously and in step with the diversification of other Tanganyikan lineages-the Limnochromini and Cyprichromini-with B. graueri occupying the most ancestral branch, suggesting that these were induced by the same environmental changes. The lack of geographic color morphs suggests that competition and resource partitioning, rather than allopatric speciation, promoted speciation within the genus Bathybates.

A highly unexpected strong correlation between fixation probability of nonsynonymous mutations and mutation rate

Under prevailing theories, the nonsynonymous-to-synonymous substitution ratio (i.e. K(a)/K(s)), which measures the fixation probability of nonsynonymous mutations, is correlated with the strength of selection. In this article, we report that K(a)/K(s) is also strongly correlated with the mutation rate as measured by K(s), and that this correlation appears to have a similar magnitude as the correlation between K(a)/K(s) and selective strength. This finding cannot be reconciled with current theories. It suggests that we should re-evaluate the current paradigms of coding-sequence evolution, and that the wide use of K(a)/K(s) as a measure of selective strength needs reassessment.

Size of the protein-coding genome and rate of molecular evolution

In diploid populations of size N, there will be 2 Nmu mutations per nucleotide (nt) site (or per locus) per generation (mu stands for mutation rate). If either the population or the coding genome double in size, one expects 4 Nmu mutations. What is important is not the population size per se but the number of genes (coding sites), the two being often interconverted. Here we compared the total physical length of protein-coding genomes (n) with the corresponding absolute rates of synonymous substitution (K(S)), an empirical neutral reference. In the classical occupancy problem and in the coupons collector (CC) problem, n was expressed as the mean rate of change (K(CC)). Despite inherently very low power of the approaches involving averaging of rates, the mode of molecular evolution of the total size phenotype of the coding genome could be evidenced through differences between the genomic estimates of K(CC) [K(CC)=1/(ln n + 0.57721) n] and rate of molecular evolution, K(S). We found that (1) the estimates of n and K(S) are reciprocally correlated across taxa (r=0.812; p<< 0.001); (2) the gamete-cell division hypothesis (Chang et al. Proc Natl Acad Sci USA 91:827-831, 1994) can be confirmed independently in terms of K(CC)/K(S) ratios; (3) the time scale of molecular evolution changes with change in mutation rate, as previously shown by Takahata (Proc Natl Acad Sci USA 87:2419-2423, 1990), Takahata et al. (Genetics 130:925-938, 1992), and Vekemans and Slatkin (Genetics 137:1157-1165, 1994); (4) the generation time and population size (Lynch and Conery, Science 302:1401-1404, 2003) effects left their "signatures" at the level of the size phenotype of the protein-coding genome.

The influence of neighboring-nucleotide composition on single nucleotide polymorphisms (SNPs) in the mouse genome and its comparison with human SNPs

We analyzed the neighboring-nucleotide composition of 433,192 biallelic substitutions, representing the largest public collection of SNPs across the mouse genome. Large neighboring-nucleotide biases relative to the genome- or chromosome-specific average were observed at the immediate adjacent sites and small biases extended farther from the substitution site. For all substitutions, the biases for A, C, G, and T were 0.21, 2.63, 0.71, and -3.55%, respectively, on the immediate adjacent 5' site and -3.67, 0.75, 2.69, and 0.23%, respectively, on the immediate adjacent 3' side. Further examination of the six categories of substitution revealed that the neighboring-nucleotide patterns for transitions were strongly influenced by the hypermutability of dinucleotide CpG and the neighboring effects on transversions were complex. Probability of a transversion increased with increasing A + T content of the two immediate adjacent sites, which was similarly observed in the human and Arabidopsis genomes. Overall, the bias patterns for the neighboring nucleotides in the mouse and human genomes were essentially the same; however, the extent of the biases was notably less in mice. Our results provide the first comprehensive view of the neighboring-nucleotide effects in the mouse genome and are important for understanding the mutational mechanisms and sequence evolution in the mammalian genomes.

Fur Seal Adaptations to Lactation: Insights into Mammary Gland Function

The fur seal (Arctocephalus spp. and Callorhinus spp., members of the pinniped family) is a mammal with the unusual capability to modulate its lactation cycle by turning milk production on and off without the typical mammalian regression and involution of the mammary gland. Lactation has evolved from constraints arising from the spatial and temporal separation of infant nursing and maternal foraging as the mother gives birth and feeds the pup on land while acquisition of nutrients for milk production occurs at sea. The lactation cycle begins with the female fur seal undergoing a perinatal fast of approximately 1 wk, after which time she departs the breeding colony to forage at sea. For the remainder of the long lactation period (116-540 days), the mother alternates between short periods ashore suckling the young with longer periods of up to 4 wk of foraging at sea. Milk production continues while foraging at sea, but at less than 20% the rate of production on land. Fur seals produce one of the richest milk reported, with a very high lipid content contributing up to 85% of total energy. This feature serves as an adaptation to the young's need to produce an insulating blubber layer against heat loss and to serve as an energy store when the mother is away foraging at sea. This atypical pattern of lactation means mothers have long periods with no suckling stimulus and can transfer high-energy milk rapidly while on land to minimize time away from foraging grounds. The absence of suckling stimulus and milk removal during foraging does not result in the onset of involution with associated apoptosis of mammary secretory cells and a subsequent progressive breakdown of the cellular structure of the mammary gland. The mechanisms controlling lactation in the fur seal mammary gland have been investigated using molecular and cellular techniques. These findings have shed light on the processes by which the unique features of lactation in the fur seal are regulated.

Inverse relationship between evolutionary rate and age of mammalian genes

A large number of genes is shared by all living organisms, whereas many others are unique to some specific lineages, indicating their different times of origin. The availability of a growing number of eukaryotic genomes allows us to estimate which mammalian genes are novel genes and, approximately, when they arose. In this article, we classify human genes into four different age groups and estimate evolutionary rates in human and mouse orthologs. We show that older genes tend to evolve more slowly than newer ones; that is, proteins that arose earlier in evolution currently have a larger proportion of sites subjected to negative selection. Interestingly, this property is maintained when a fraction of the fastest-evolving genes is excluded or when only genes belonging to a given functional class are considered. One way to explain this relationship is by assuming that genes maintain their functional constraints along all their evolutionary history, but the nature of more recent evolutionary innovations is such that the functional constraints operating on them are increasingly weaker. Alternatively, our results would also be consistent with a scenario in which the functional constraints acting on a gene would not need to be constant through evolution. Instead, starting from weak functional constraints near the time of origin of a gene-as supported by mechanisms proposed for the origin of orphan genes-there would be a gradual increase in selective pressures with time, resulting in fewer accepted mutations in older versus more novel genes.

Insertions and deletions are male biased too: a whole-genome analysis in rodents

It is presently accepted that, in mammals, due to the greater number of cell divisions in the male germline than in the female germline, nucleotide substitutions occur more frequently in males. The data on mutation bias in insertions and deletions (indels) are contradictory, with some studies indicating no sex bias and others indicating either female or male bias. The sequenced rat and mouse genomes provide a unique opportunity to investigate a potential sex bias for different types of mutations. Indeed, mutation rates can be accurately estimated from a large number of orthologous loci in organisms similar in generation time and in the number of germline cell divisions. Here we compare the mutation rates between chromosome X and autosomes for likely neutral sites in eutherian ancestral interspersed repetitive elements present at orthologous locations in the rat and mouse genomes. We find that small indels are male biased: The male-to-female mutation rate ratio (alpha) for indels in rodents is approximately 2. Similarly, our whole-genome analysis in rodents indicates an approximately twofold excess of nucleotide substitutions originating in males over that in females. This is the same as the male-to-female ratio of the number of germline cell divisions in rat and mouse. Thus, this is consistent with nucleotide substitutions and small indels occurring primarily during DNA replication.

Genome sequence of the Brown Norway rat yields insights into mammalian evolution

The laboratory rat (Rattus norvegicus) is an indispensable tool in experimental medicine and drug development, having made inestimable contributions to human health. We report here the genome sequence of the Brown Norway (BN) rat strain. The sequence represents a high-quality 'draft' covering over 90% of the genome. The BN rat sequence is the third complete mammalian genome to be deciphered, and three-way comparisons with the human and mouse genomes resolve details of mammalian evolution. This first comprehensive analysis includes genes and proteins and their relation to human disease, repeated sequences, comparative genome-wide studies of mammalian orthologous chromosomal regions and rearrangement breakpoints, reconstruction of ancestral karyotypes and the events leading to existing species, rates of variation, and lineage-specific and lineage-independent evolutionary events such as expansion of gene families, orthology relations and protein evolution.