Biochemical networking contributes more to genetic buffering in human and mouse metabolic pathways than does gene duplication

During evolution different genes evolve at unequal rates, reflecting the varying functional constraints on phenotype. An important contributor to this variation is genetic buffering, which reduces the potential detrimental effects of mutations. We studied whether gene duplication and redundant metabolic networks affect genetic buffering by comparing the evolutionary rate of 242 human and mouse orthologous genes involved in metabolic pathways. A gene with a redundant network is defined as one for which the structural layout of metabolic pathways provides an alternative metabolic route that can, in principle, compensate for the loss of a protein function encoded by the gene. We found that genes with redundant networks evolve at similar rates as did genes without redundant networks, [corrected] but no significant difference was detected between single-copy genes and gene families. This implies that redundancy in metabolic networks provides significantly more genetic buffering than do gene families. We also found that genes encoding proteins involved in glycolysis and gluconeogenesis showed as a group a distinct pattern of variation, in contrast to genes involved in other pathways. These results suggest that redundant networks provide genetic buffering and contribute to the functional diversification of metabolic pathways.     

Birth of a metabolic gene cluster in yeast by adaptive gene relocation

Although most eukaryotic genomes lack operons, they contain some physical clusters of genes that are related in function despite being unrelated in sequence. How these clusters are formed during evolution is unknown. The DAL cluster is the largest metabolic gene cluster in yeast and consists of six adjacent genes encoding proteins that enable Saccharomyces cerevisiae to use allantoin as a nitrogen source. We show here that the DAL cluster was assembled, quite recently in evolutionary terms, through a set of genomic rearrangements that happened almost simultaneously. Six of the eight genes involved in allantoin degradation, which were previously scattered around the genome, became relocated to a single subtelomeric site in an ancestor of S. cerevisiae and Saccharomyces castellii. These genomic rearrangements coincided with a biochemical reorganization of the purine degradation pathway, which switched to importing allantoin instead of urate. This change eliminated urate oxidase, one of several oxygen-consuming enzymes that were lost by yeasts that can grow vigorously in anaerobic conditions. The DAL cluster is located in a domain of modified chromatin involving both H2A.Z histone exchange and Hst1-Sum1-mediated histone deacetylation, and it may be a coadapted gene complex formed by epistatic selection.     

Universal trees based on large combined protein sequence data sets.

Universal trees of life based on small-subunit (SSU) ribosomal RNA (rRNA) support the separate mono/holophyly of the domains Archaea (archaebacteria), Bacteria (eubacteria) and Eucarya (eukaryotes) and the placement of extreme thermophiles at the base of the Bacteria. The concept of universal tree reconstruction recently has been upset by protein trees that show intermixing of species from different domains. Such tree topologies have been attributed to either extensive horizontal gene transfer or degradation of phylogenetic signals because of saturation for amino acid substitutions. Here we use large combined alignments of 23 orthologous proteins conserved across 45 species from all domains to construct highly robust universal trees. Although individual protein trees are variable in their support of domain integrity, trees based on combined protein data sets strongly support separate monophyletic domains. Within the Bacteria, we placed spirochaetes as the earliest derived bacterial group. However, elimination from the combined protein alignment of nine protein data sets, which were likely candidates for horizontal gene transfer, resulted in trees showing thermophiles as the earliest evolved bacterial lineage. Thus, combined protein universal trees are highly congruent with SSU rRNA trees in their strong support for the separate monophyly of domains as well as the early evolution of thermophilic Bacteria.     

Natural variation in cardiac metabolism and gene expression in Fundulus heteroclitus

Individual variation in gene expression is important for evolutionary adaptation and susceptibility to diseases and pathologies. In this study, we address the functional importance of this variation by comparing cardiac metabolism to patterns of mRNA expression using microarrays. There is extensive variation in both cardiac metabolism and the expression of metabolic genes among individuals of the teleost fish Fundulus heteroclitus from natural outbred populations raised in a common environment: metabolism differed among individuals by a factor of more than 2, and expression levels of 94% of genes were significantly different (P < 0.01) between individuals in a population. This unexpectedly high variation in metabolic gene expression explains much of the variation in metabolism, suggesting that it is biologically relevant. The patterns of gene expression that are most important in explaining cardiac metabolism differ between groups of individuals. Apparently, the variation in metabolism seems to be related to different patterns of gene expression in the different groups of individuals. The magnitude of differences in gene expression in these groups is not important; large changes in expression have no greater predictive value than small changes. These data suggest that variation in physiological performance is related to the subtle variation in gene expression and that this relationship differs among individuals.     

Heterogeneous mutation processes in human microsatellite DNA sequences

Although microsatellite polymorphisms are one of the most commonly used tools in genetic analyses, it remains to be understood how microsatellite DNA has evolved as a ubiquitous and highly abundant class of repetitive sequences in eukaryotic genomes. On the basis of analyses of spontaneous human microsatellite mutations of germline origin, I show here that different mutation biases underlie the evolution of microsatellite repeats. The within-locus mutation rate increases with allele length, but is not affected by the size difference between an individual's two alleles (allele span). Within loci, long alleles tend to mutate to shorter lengths, thereby acting to prevent infinite growth. Expansions are more common than contractions among dinucleotide repeats, whereas no such trend is evident among tetranucleotide repeats. This observation is consistent with the longer repeat lengths and higher frequency of di- compared with tetranucleotide repeats. An excess of paternally transmitted mutations (male-to-female ratio of 4.9) supports a male-biased mutation rate in the human genome.     

Adaptive evolution of bacterial metabolic networks by horizontal gene transfer

Numerous studies have considered the emergence of metabolic pathways, but the modes of recent evolution of metabolic networks are poorly understood. Here, we integrate comparative genomics with flux balance analysis to examine (i) the contribution of different genetic mechanisms to network growth in bacteria, (ii) the selective forces driving network evolution and (iii) the integration of new nodes into the network. Most changes to the metabolic network of Escherichia coli in the past 100 million years are due to horizontal gene transfer, with little contribution from gene duplicates. Networks grow by acquiring genes involved in the transport and catalysis of external nutrients, driven by adaptations to changing environments. Accordingly, horizontally transferred genes are integrated at the periphery of the network, whereas central parts remain evolutionarily stable. Genes encoding physiologically coupled reactions are often transferred together, frequently in operons. Thus, bacterial metabolic networks evolve by direct uptake of peripheral reactions in response to changed environments.     

Adaptive evolution of a duplicated pancreatic ribonuclease gene in a leaf-eating monkey

Although the complete genome sequences of over 50 representative species have revealed the many duplicated genes in all three domains of life, the roles of gene duplication in organismal adaptation and biodiversity are poorly understood. In addition, the evolutionary forces behind the functional divergence of duplicated genes are often unknown, leading to disagreement on the relative importance of positive Darwinian selection versus relaxation of functional constraints in this process. The methodology of earlier studies relied largely on DNA sequence analysis but lacked functional assays of duplicated genes, frequently generating contentious results. Here we use both computational and experimental approaches to address these questions in a study of the pancreatic ribonuclease gene (RNASE1) and its duplicate gene (RNASE1B) in a leaf-eating colobine monkey, douc langur. We show that RNASE1B has evolved rapidly under positive selection for enhanced ribonucleolytic activity in an altered microenvironment, a response to increased demands for the enzyme for digesting bacterial RNA. At the same time, the ability to degrade double-stranded RNA, a non-digestive activity characteristic of primate RNASE1, has been lost in RNASE1B, indicating functional specialization and relaxation of purifying selection. Our findings demonstrate the contribution of gene duplication to organismal adaptation and show the power of combining sequence analysis and functional assays in delineating the molecular basis of adaptive evolution.     

ASPM is a major determinant of cerebral cortical size

One of the most notable trends in mammalian evolution is the massive increase in size of the cerebral cortex, especially in primates. Humans with autosomal recessive primary microcephaly (MCPH) show a small but otherwise grossly normal cerebral cortex associated with mild to moderate mental retardation. Genes linked to this condition offer potential insights into the development and evolution of the cerebral cortex. Here we show that the most common cause of MCPH is homozygous mutation of ASPM, the human ortholog of the Drosophila melanogaster abnormal spindle gene (asp), which is essential for normal mitotic spindle function in embryonic neuroblasts. The mouse gene Aspm is expressed specifically in the primary sites of prenatal cerebral cortical neurogenesis. Notably, the predicted ASPM proteins encode systematically larger numbers of repeated 'IQ' domains between flies, mice and humans, with the predominant difference between Aspm and ASPM being a single large insertion coding for IQ domains. Our results and evolutionary considerations suggest that brain size is controlled in part through modulation of mitotic spindle activity in neuronal progenitor cells.     

Diversity of microRNAs in human and chimpanzee brain

We used massively parallel sequencing to compare the microRNA (miRNA) content of human and chimpanzee brains, and we identified 447 new miRNA genes. Many of the new miRNAs are not conserved beyond primates, indicating their recent origin, and some miRNAs seem species specific, whereas others are expanded in one species through duplication events. These data suggest that evolution of miRNAs is an ongoing process and that along with ancient, highly conserved miRNAs, there are a number of emerging miRNAs.     

Mutations in the gene encoding 11-cis retinol dehydrogenase cause delayed dark adaptation and fundus albipunctatus.

The metabolic pathways that produce 11-cis retinal are important for vision because this retinoid is the chromophore residing in rhodopsin and the cone opsins. The all-trans retinal that is generated after cone and rod photopigments absorb photons of light is recycled back to 11-cis retinal by the retinal pigment epithelium and Müller cells of the retina. Several of the enzymes involved have recently been purified and molecularly cloned; here we focus on 11-cis retinol dehydrogenase (encoded by the gene RDH5; chromosome 12q13-14; ref. 4), the first cloned enzyme in this pathway. This microsomal enzyme is abundant in the retinal pigment epithelium, where it has been proposed to catalyse the conversion of 11-cis retinol to 11-cis retinal. We evaluated patients with hereditary retinal diseases featuring subretinal spots (retinitis punctata albescens and fundus albipunctatus) and patients with typical dominant or recessive retinitis pigmentosa for mutations in RDH5. Mutations were found only in two unrelated patients, both with fundus albipunctatus; they segregated with disease in the respective families. Recombinant mutant 11-cis retinol dehydrogenases had reduced activity compared with recombinant enzyme with wild-type sequence. Our results suggest that mutant alleles in RDH5 are a cause of fundus albipunctatus, a rare form of stationary night blindness characterized by a delay in the regeneration of cone and rod photopigments.     

The Pristionchus pacificus genome provides a unique perspective on nematode lifestyle and parasitism

Here we present a draft genome sequence of the nematode Pristionchus pacificus, a species that is associated with beetles and is used as a model system in evolutionary biology. With 169 Mb and 23,500 predicted protein-coding genes, the P. pacificus genome is larger than those of Caenorhabditis elegans and the human parasite Brugia malayi. Compared to C. elegans, the P. pacificus genome has more genes encoding cytochrome P450 enzymes, glucosyltransferases, sulfotransferases and ABC transporters, many of which were experimentally validated. The P. pacificus genome contains genes encoding cellulase and diapausin, and cellulase activity is found in P. pacificus secretions, indicating that cellulases can be found in nematodes beyond plant parasites. The relatively higher number of detoxification and degradation enzymes in P. pacificus is consistent with its necromenic lifestyle and might represent a preadaptation for parasitism. Thus, comparative genomics analysis of three ecologically distinct nematodes offers a unique opportunity to investigate the association between genome structure and lifestyle.     

Population genomic analysis of outcrossing and recombination in yeast

The budding yeast Saccharomyces cerevisiae has been used by humans for millennia to make wine, beer and bread. More recently, it became a key model organism for studies of eukaryotic biology and for genomic analysis. However, relatively little is known about the natural lifestyle and population genetics of yeast. One major question is whether genetically diverse yeast strains mate and recombine in the wild. We developed a method to infer the evolutionary history of a species from genome sequences of multiple individuals and applied it to whole-genome sequence data from three strains of Saccharomyces cerevisiae and the sister species Saccharomyces paradoxus. We observed a pattern of sequence variation among yeast strains in which ancestral recombination events lead to a mosaic of segments with shared genealogy. Based on sequence divergence and the inferred median size of shared segments (approximately 2,000 bp), we estimated that although any two strains have undergone approximately 16 million cell divisions since their last common ancestor, only 314 outcrossing events have occurred during this time (roughly one every 50,000 divisions). Local correlations in polymorphism rates indicate that linkage disequilibrium in yeast should extend over kilobases. Our results provide the initial foundation for population studies of association between genotype and phenotype in S. cerevisiae.     

A mutation in OTOF, encoding otoferlin, a FER-1-like protein, causes DFNB9, a nonsyndromic form of deafness

Using a candidate gene approach, we identified a novel human gene, OTOF, underlying an autosomal recessive, nonsyndromic prelingual deafness, DFNB9. The same nonsense mutation was detected in four unrelated affected families of Lebanese origin. OTOF is the second member of a mammalian gene family related to Caenorhabditis elegans fer-1. It encodes a predicted cytosolic protein (of 1,230 aa) with three C2 domains and a single carboxy-terminal transmembrane domain. The sequence homologies and predicted structure of otoferlin, the protein encoded by OTOF, suggest its involvement in vesicle membrane fusion. In the inner ear, the expression of the orthologous mouse gene, mainly in the sensory hair cells, indicates that such a role could apply to synaptic vesicles.     

Distributions of epistasis in microbes fit predictions from a fitness landscape model

How do the fitness effects of several mutations combine? Despite its simplicity, this question is central to the understanding of multilocus evolution. Epistasis (the interaction between alleles at different loci), especially epistasis for fitness traits such as reproduction and survival, influences evolutionary predictions "almost whenever multilocus genetics matters". Yet very few models have sought to predict epistasis, and none has been empirically tested. Here we show that the distribution of epistasis can be predicted from the distribution of single mutation effects, based on a simple fitness landscape model. We show that this prediction closely matches the empirical measures of epistasis that have been obtained for Escherichia coli and the RNA virus vesicular stomatitis virus. Our results suggest that a simple fitness landscape model may be sufficient to quantitatively capture the complex nature of gene interactions. This model may offer a simple and widely applicable alternative to complex metabolic network models, in particular for making evolutionary predictions.     

Duplication-degeneration as a mechanism of gene fission and the origin of new genes in Drosophila species

Gene fission and fusion, the processes by which a single gene is split into two separate genes and two adjacent genes are fused into a single gene, respectively, are among the primary processes that generate new genes. Despite their seeming reversibility, nothing is known about the mechanism of gene fission. Because the nucleotide sequences of fission genes record little about their origination process, conventional analysis of duplicate genes may not be powerful enough to unravel the underlying mechanism. In a survey for young genes in species of the Drosophila melanogaster subgroup using fluorescence in situ hybridization, we identified a young gene family, monkey king, whose genesis sheds light on the evolutionary process of gene fission. Its members originated 1-2 million years ago as retroposed duplicates and evolved into fission genes that separately encode protein domains from a multidomain ancestor. The mechanism underlying this process is gene duplication with subsequent partial degeneration.     

The genome of the extremophile crucifer Thellungiella parvula

Thellungiella parvula is related to Arabidopsis thaliana and is endemic to saline, resource-poor habitats, making it a model for the evolution of plant adaptation to extreme environments. Here we present the draft genome for this extremophile species. Exclusively by next generation sequencing, we obtained the de novo assembled genome in 1,496 gap-free contigs, closely approximating the estimated genome size of 140 Mb. We anchored these contigs to seven pseudo chromosomes without the use of maps. We show that short reads can be assembled to a near-complete chromosome level for a eukaryotic species lacking prior genetic information. The sequence identifies a number of tandem duplications that, by the nature of the duplicated genes, suggest a possible basis for T. parvula's extremophile lifestyle. Our results provide essential background for developing genomically influenced testable hypotheses for the evolution of environmental stress tolerance.