Modularity and evolutionary constraint on proteins

Modularity, which has been found in the functional and physical protein interaction networks of many organisms, has been postulated to affect both the mode and tempo of evolution. Here I show that in the yeast Saccharomyces cerevisiae, protein interaction hubs situated in single modules are highly constrained, whereas those connecting different modules are more plastic. This pattern of change could reflect a tendency for evolutionary innovations to occur by altering the proteins and interactions between rather than within modules, in a manner somewhat similar to the evolution of new proteins through the shuffling of conserved protein domains.     

Multiple knockout analysis of genetic robustness in the yeast metabolic network

Genetic robustness characterizes the constancy of the phenotype in face of heritable perturbations. Previous investigations have used comprehensive single and double gene knockouts to study gene essentiality and pairwise gene interactions in the yeast Saccharomyces cerevisiae. Here we conduct an in silico multiple knockout investigation of a flux balance analysis model of the yeast's metabolic network. Cataloging gene sets that provide mutual functional backup, we identify sets of up to eight interacting genes and characterize the 'k robustness' (the depth of backup interactions) of each gene. We find that 74% (360) of the metabolic genes participate in processes that are essential to growth in a standard laboratory environment, compared with only 13% previously found to be essential using single knockouts. The genes' k robustness is shown to be a solid indicator of their biological buffering capacity and is correlated with both the genes' environmental specificity and their evolutionary retention.     

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.     

A high-resolution HLA and SNP haplotype map for disease association studies in the extended human MHC

The proteins encoded by the classical HLA class I and class II genes in the major histocompatibility complex (MHC) are highly polymorphic and are essential in self versus non-self immune recognition. HLA variation is a crucial determinant of transplant rejection and susceptibility to a large number of infectious and autoimmune diseases. Yet identification of causal variants is problematic owing to linkage disequilibrium that extends across multiple HLA and non-HLA genes in the MHC. We therefore set out to characterize the linkage disequilibrium patterns between the highly polymorphic HLA genes and background variation by typing the classical HLA genes and >7,500 common SNPs and deletion-insertion polymorphisms across four population samples. The analysis provides informative tag SNPs that capture much of the common variation in the MHC region and that could be used in disease association studies, and it provides new insight into the evolutionary dynamics and ancestral origins of the HLA loci and their haplotypes.     

Flowering-time genes modulate meristem determinacy and growth form in Arabidopsis thaliana

Plants have evolved annual and perennial life forms as alternative strategies to adapt reproduction and survival to environmental constraints. In isolated situations, such as islands, woody perennials have evolved repeatedly from annual ancestors. Although the molecular basis of the rapid evolution of insular woodiness is unknown, the molecular difference between perennials and annuals might be rather small, and a change between these life strategies might not require major genetic innovations. Developmental regulators can strongly affect evolutionary variation and genes involved in meristem transitions are good candidates for a switch in growth habit. We found that the MADS box proteins SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) and FRUITFULL (FUL) not only control flowering time, but also affect determinacy of all meristems. In addition, downregulation of both proteins established phenotypes common to the lifestyle of perennial plants, suggesting their involvement in the prevention of secondary growth and longevity in annual life forms.     

Distribution and functional impact of DNA copy number variation in the rat

The abundance and dynamics of copy number variants (CNVs) in mammalian genomes poses new challenges in the identification of their impact on natural and disease phenotypes. We used computational and experimental methods to catalog CNVs in rat and found that they share important functional characteristics with those in human. In addition, 113 one-to-one orthologous genes overlap CNVs in both human and rat, 80 of which are implicated in human disease. CNVs are nonrandomly distributed throughout the genome. Chromosome 18 is a cold spot for CNVs as well as evolutionary rearrangements and segmental duplications, suggesting stringent selective mechanisms underlying CNV genesis or maintenance. By exploiting gene expression data available for rat recombinant inbred lines, we established the functional relationship of CNVs underlying 22 expression quantitative trait loci. These characteristics make the rat an excellent model for studying phenotypic effects of structural variation in relation to human complex traits and disease.     

Evolutionary conservation of motif constituents in the yeast protein interaction network

Understanding why some cellular components are conserved across species but others evolve rapidly is a key question of modern biology. Here we show that in Saccharomyces cerevisiae, proteins organized in cohesive patterns of interactions are conserved to a substantially higher degree than those that do not participate in such motifs. We find that the conservation of proteins in distinct topological motifs correlates with the interconnectedness and function of that motif and also depends on the structure of the overall interactome topology. These findings indicate that motifs may represent evolutionary conserved topological units of cellular networks molded in accordance with the specific biological function in which they participate.     

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.     

Epistasis and the release of genetic variation during long-term selection

It is an enigma how long-term selection in model organisms and agricultural species can lead to marked phenotypic changes without exhausting genetic variation for the selected trait. Here, we show that the genetic architecture of an apparently major locus for growth in chicken dissects into a genetic network of four interacting loci. The interactions in this radial network mediate a considerably larger selection response than predicted by a single-locus model.     

Module networks: identifying regulatory modules and their condition-specific regulators from gene expression data

Much of a cell's activity is organized as a network of interacting modules: sets of genes coregulated to respond to different conditions. We present a probabilistic method for identifying regulatory modules from gene expression data. Our procedure identifies modules of coregulated genes, their regulators and the conditions under which regulation occurs, generating testable hypotheses in the form 'regulator X regulates module Y under conditions W'. We applied the method to a Saccharomyces cerevisiae expression data set, showing its ability to identify functionally coherent modules and their correct regulators. We present microarray experiments supporting three novel predictions, suggesting regulatory roles for previously uncharacterized proteins.     

Control of neurulation by the nucleosome assembly protein-1-like 2

Neurulation is a complex process of histogenesis involving the precise temporal and spatial organization of gene expression. Genes influencing neurulation include proneural genes determining primary cell fate, neurogenic genes involved in lateral inhibition pathways and genes controlling the frequency of mitotic events. This is reflected in the aetiology and genetics of human and mouse neural tube defects, which are of both multifactorial and multigenic origin. The X-linked gene Nap1l2, specifically expressed in neurons, encodes a protein that is highly similar to the nucleosome assembly (NAP) and SET proteins. We inactivated Nap1l2 in mice by gene targeting, leading to embryonic lethality from mid-gestation onwards. Surviving mutant chimaeric embryos showed extensive surface ectoderm defects as well as the presence of open neural tubes and exposed brains similar to those observed in human spina bifida and anencephaly. These defects correlated with an overproduction of neuronal precursor cells. Protein expression studies showed that the Nap1l2 protein binds to condensing chromatin during S phase and in apoptotic cells, but remained cytoplasmic during G1 phase. Nap1l2 therefore likely represents a class of tissue-specific factors interacting with chromatin to regulate neuronal cell proliferation.     

A direct functional link between the multi-PDZ domain protein GRIP1 and the Fraser syndrome protein Fras1

Cell adhesion to extracellular matrix (ECM) proteins is crucial for the structural integrity of tissues and epithelial-mesenchymal interactions mediating organ morphogenesis. Here we describe how the loss of a cytoplasmic multi-PDZ scaffolding protein, glutamate receptor interacting protein 1 (GRIP1), leads to the formation of subepidermal hemorrhagic blisters, renal agenesis, syndactyly or polydactyly and permanent fusion of eyelids (cryptophthalmos). Similar malformations are characteristic of individuals with Fraser syndrome and animal models of this human genetic disorder, such as mice carrying the blebbed mutation (bl) in the gene encoding the Fras1 ECM protein. GRIP1 can physically interact with Fras1 and is required for the localization of Fras1 to the basal side of cells. In one animal model of Fraser syndrome, the eye-blebs (eb) mouse, Grip1 is disrupted by a deletion of two coding exons. Our data indicate that GRIP1 is required for normal cell-matrix interactions during early embryonic development and that inactivation of Grip1 causes Fraser syndrome-like defects in mice.