Genetic basis of individual differences in the response to small-molecule drugs in yeast

Individual response to small-molecule drugs is variable; a drug that provides a cure for some may confer no therapeutic benefit or trigger an adverse reaction in others. To begin to understand such differences systematically, we treated 104 genotyped segregants from a cross between two yeast strains with a collection of 100 diverse small molecules. We used linkage analysis to identify 124 distinct linkages between genetic markers and response to 83 compounds. The linked markers clustered at eight genomic locations, or quantitative-trait locus 'hotspots', that contain one or more polymorphisms that affect response to multiple small molecules. We also experimentally verified that a deficiency in leucine biosynthesis caused by a deletion of LEU2 underlies sensitivity to niguldipine, which is structurally related to therapeutic calcium channel blockers, and that a natural coding-region polymorphism in the inorganic phosphate transporter PHO84 underlies sensitivity to two polychlorinated phenols that uncouple oxidative phosphorylation. Our results provide a step toward a systematic understanding of small-molecule drug action in genetically distinct individuals.     

Genetic analysis of the mouse brain proteome

Proteome analysis is a fundamental step in systematic functional genomics. Here we have resolved 8,767 proteins from the mouse brain proteome by large-gel two-dimensional electrophoresis. We detected 1,324 polymorphic proteins from the European collaborative interspecific backcross. Of these, we mapped 665 proteins genetically and identified 466 proteins by mass spectrometry. Qualitatively polymorphic proteins, to 96%, reflect changes in conformation and/or mass. Quantitatively polymorphic proteins show a high frequency (73%) of allele-specific transmission in codominant heterozygotes. Variations in protein isoforms and protein quantity often mapped to chromosomal positions different from that of the structural gene, indicating that single proteins may act as polygenic traits. Genetic analysis of proteomes may detect the types of polymorphism that are most relevant in disease-association studies.     

Genetic variation in laboratory mice

Characterizing the patterns of genetic variation in an organism provides fundamental insight into the evolutionary history of the organism and defines the scope and nature of studies that must be designed to correlate genotype to phenotype. Given the pre-eminent role of the inbred mouse in biomedical research, considerable effort has been undertaken in recent years to describe more fully the nature and amount of genetic variation among the numerous strains of mice that are in widest use. Here, we discuss recent studies that have contributed to an emerging understanding of the unique variation patterns found in inbred strains of mice and how they have arisen through a combination of natural evolution and human-directed breeding. These preliminary results have ramifications for genetic research into complex biomedical traits and are the basis for the development of future variation resources.     

Predicting phenotypic variation in yeast from individual genome sequences

A central challenge in genetics is to predict phenotypic variation from individual genome sequences. Here we construct and evaluate phenotypic predictions for 19 strains of Saccharomyces cerevisiae. We use conservation-based methods to predict the impact of protein-coding variation within genes on protein function. We then rank strains using a prediction score that measures the total sum of function-altering changes in different sets of genes reported to influence over 100 phenotypes in genome-wide loss-of-function screens. We evaluate our predictions by comparing them with the observed growth rate and efficiency of 15 strains tested across 20 conditions in quantitative experiments. The median predictive performance, as measured by ROC AUC, was 0.76, and predictions were more accurate when the genes reported to influence a trait were highly connected in a functional gene network.     

Genetic variation in selenoprotein S influences inflammatory response

Chronic inflammation has a pathological role in many common diseases and is influenced by both genetic and environmental factors. Here we assess the role of genetic variation in selenoprotein S (SEPS1, also called SELS or SELENOS), a gene involved in stress response in the endoplasmic reticulum and inflammation control. After resequencing SEPS1, we genotyped 13 SNPs in 522 individuals from 92 families. As inflammation biomarkers, we measured plasma levels of IL-6, IL-1beta and TNF-alpha. Bayesian quantitative trait nucleotide analysis identified associations between SEPS1 polymorphisms and all three proinflammatory cytokines. One promoter variant, -105G --> A, showed strong evidence for an association with each cytokine (multivariate P = 0.0000002). Functional analysis of this polymorphism showed that the A variant significantly impaired SEPS1 expression after exposure to endoplasmic reticulum stress agents (P = 0.00006). Furthermore, suppression of SEPS1 by short interfering RNA in macrophage cells increased the release of IL-6 and TNF-alpha. To investigate further the significance of the observed associations, we genotyped -105G --> A in 419 Mexican American individuals from 23 families for replication. This analysis confirmed a significant association with both TNF-alpha (P = 0.0049) and IL-1beta (P = 0.0101). These results provide a direct mechanistic link between SEPS1 and the production of inflammatory cytokines and suggest that SEPS1 has a role in mediating inflammation.     

Genetic basis of total colourblindness among the Pingelapese islanders

Complete achromatopsia is a rare, autosomal recessive disorder characterized by photophobia, low visual acuity, nystagmus and a total inability to distinguish colours. In this disease, cone photoreceptors, the retinal sensory neurons mediating colour vision, seem viable but fail to generate an electrical response to light. Achromatopsia, or rod monochromatism, was first mapped to 2p11-2q12 (MIM 216900; ref. 3), where it is associated with missense mutations in CNGA3 (ref. 4). CNGA3 encodes the alpha-subunit of the cone cyclic nucleotide-gated cation channel, which generates the light-evoked electrical responses of cone photoreceptors. A second locus at 8q21-q22 has been identified among the Pingelapese islanders of Micronesia, who have a high incidence of recessive achromatopsia (MIM 262300). Here we narrow the achromatopsia locus to 1.4 cM and show that Pingelapese achromatopsia segregates with a missense mutation at a highly conserved site in CNGB3, a new gene that encodes the beta-subunit of the cone cyclic nucleotide-gated cation channel. Two independent frameshift deletions establish that achromatopsia is the null phenotype of CNGB3. Combined with earlier findings, our results demonstrate that both alpha- and beta-subunits of the cGMP-gated channel are essential for phototransduction in all three classes of cones.     

Genetic variation in DLG5 is associated with inflammatory bowel disease

Crohn disease and ulcerative colitis are two subphenotypes of inflammatory bowel disease (IBD), a complex disorder resulting from gene-environment interaction. We refined our previously defined linkage region for IBD on chromosome 10q23 and used positional cloning to identify genetic variants in DLG5 associated with IBD. DLG5 encodes a scaffolding protein involved in the maintenance of epithelial integrity. We identified two distinct haplotypes with a replicable distortion in transmission (P = 0.000023 and P = 0.004 for association with IBD, P = 0.00012 and P = 0.04 for association with Crohn disease). One of the risk-associated DLG5 haplotypes is distinguished from the common haplotype by a nonsynonymous single-nucleotide polymorphism 113G-->A, resulting in the amino acid substitution R30Q in the DUF622 domain of DLG5. This mutation probably impedes scaffolding of DLG5. We stratified the study sample according to the presence of risk-associated CARD15 variants to study potential gene-gene interaction. We found a significant difference in association of the 113A DLG5 variant with Crohn disease in affected individuals carrying the risk-associated CARD15 alleles versus those carrying non-risk-associated CARD15 alleles. This is suggestive of a complex pattern of gene-gene interaction between DLG5 and CARD15, reflecting the complex nature of polygenic diseases. Further functional studies will evaluate the biological significance of DLG5 variants.     

Genetic variation in PNPLA3 confers susceptibility to nonalcoholic fatty liver disease

Nonalcoholic fatty liver disease (NAFLD) is a burgeoning health problem of unknown etiology that varies in prevalence among ancestry groups. To identify genetic variants contributing to differences in hepatic fat content, we carried out a genome-wide association scan of nonsynonymous sequence variations (n = 9,229) in a population comprising Hispanic, African American and European American individuals. An allele in PNPLA3 (rs738409[G], encoding I148M) was strongly associated with increased hepatic fat levels (P = 5.9 x 10(-10)) and with hepatic inflammation (P = 3.7 x 10(-4)). The allele was most common in Hispanics, the group most susceptible to NAFLD; hepatic fat content was more than twofold higher in PNPLA3 rs738409[G] homozygotes than in noncarriers. Resequencing revealed another allele of PNPLA3 (rs6006460[T], encoding S453I) that was associated with lower hepatic fat content in African Americans, the group at lowest risk of NAFLD. Thus, variation in PNPLA3 contributes to ancestry-related and inter-individual differences in hepatic fat content and susceptibility to NAFLD.     

Genetic variation in the KIF1B locus influences susceptibility to multiple sclerosis

The few loci associated with multiple sclerosis (MS) are all related to immune function. We report a GWA study identifying a new locus replicated in 2,679 cases and 3,125 controls. An rs10492972[C] variant located in the KIF1B gene was associated with MS with an odds ratio of 1.35 (P = 2.5 x 10(-10)). KIF1B is a neuronally expressed gene plausibly implicated in the irreversible axonal loss characterizing MS in the long term.     

Genetic variation in DPP6 is associated with susceptibility to amyotrophic lateral sclerosis

We identified a SNP in the DPP6 gene that is consistently strongly associated with susceptibility to amyotrophic lateral sclerosis (ALS) in different populations of European ancestry, with an overall P value of 5.04 x 10(-8) in 1,767 cases and 1,916 healthy controls and with an odds ratio of 1.30 (95% confidence interval (CI) of 1.18-1.43). Our finding is the first report of a genome-wide significant association with sporadic ALS and may be a target for future functional studies.     

Genetic variation in Mon1a affects protein trafficking and modifies macrophage iron loading in mice

We undertook a quantitative trait locus (QTL) analysis in mice to identify modifier genes that might influence the severity of human iron disorders. We identified a strong QTL on mouse chromosome 9 that differentially affected macrophage iron burden in C57BL/10J and SWR/J mice. A C57BL/10J missense allele of an evolutionarily conserved gene, Mon1a, cosegregated with the QTL in congenic mouse lines. We present evidence that Mon1a is involved in trafficking of ferroportin, the major mammalian iron exporter, to the surface of iron-recycling macrophages. Differences in amounts of surface ferroportin correlate with differences in cellular iron content. Mon1a is also important for trafficking of cell-surface and secreted molecules unrelated to iron metabolism, suggesting that it has a fundamental role in the mammalian secretory apparatus.     

Genetic variation in the 5q31 cytokine gene cluster confers susceptibility to Crohn disease

Linkage disequilibrium (LD) mapping provides a powerful method for fine-structure localization of rare disease genes, but has not yet been widely applied to common disease. We sought to design a systematic approach for LD mapping and apply it to the localization of a gene (IBD5) conferring susceptibility to Crohn disease. The key issues are: (i) to detect a significant LD signal (ii) to rigorously bound the critical region and (iii) to identify the causal genetic variant within this region. We previously mapped the IBD5 locus to a large region spanning 18 cM of chromosome 5q31 (P<10(-4)). Using dense genetic maps of microsatellite markers and single-nucleotide polymorphisms (SNPs) across the entire region, we found strong evidence of LD. We bound the region to a common haplotype spanning 250 kb that shows strong association with the disease (P< 2 x 10(-7)) and contains the cytokine gene cluster. This finding provides overwhelming evidence that a specific common haplotype of the cytokine region in 5q31 confers susceptibility to Crohn disease. However, genetic evidence alone is not sufficient to identify the causal mutation within this region, as strong LD across the region results in multiple SNPs having equivalent genetic evidence-each consistent with the expected properties of the IBD5 locus. These results have important implications for Crohn disease in particular and LD mapping in general.     

Genetic variation in the gene encoding calpain-10 is associated with type 2 diabetes mellitus

Type 2 or non-insulin-dependent diabetes mellitus (NIDDM) is the most common form of diabetes worldwide, affecting approximately 4% of the world's adult population. It is multifactorial in origin with both genetic and environmental factors contributing to its development. A genome-wide screen for type 2 diabetes genes carried out in Mexican Americans localized a susceptibility gene, designated NIDDM1, to chromosome 2. Here we describe the positional cloning of a gene located in the NIDDM1 region that shows association with type 2 diabetes in Mexican Americans and a Northern European population from the Botnia region of Finland. This putative diabetes-susceptibility gene encodes a ubiquitously expressed member of the calpain-like cysteine protease family, calpain-10 (CAPN10). This finding suggests a novel pathway that may contribute to the development of type 2 diabetes.     

Genetic variation near IRS1 associates with reduced adiposity and an impaired metabolic profile

Genome-wide association studies have identified 32 loci influencing body mass index, but this measure does not distinguish lean from fat mass. To identify adiposity loci, we meta-analyzed associations between ～2.5 million SNPs and body fat percentage from 36,626 individuals and followed up the 14 most significant (P < 10(-6)) independent loci in 39,576 individuals. We confirmed a previously established adiposity locus in FTO (P = 3 × 10(-26)) and identified two new loci associated with body fat percentage, one near IRS1 (P = 4 × 10(-11)) and one near SPRY2 (P = 3 × 10(-8)). Both loci contain genes with potential links to adipocyte physiology. Notably, the body-fat-decreasing allele near IRS1 is associated with decreased IRS1 expression and with an impaired metabolic profile, including an increased visceral to subcutaneous fat ratio, insulin resistance, dyslipidemia, risk of diabetes and coronary artery disease and decreased adiponectin levels. Our findings provide new insights into adiposity and insulin resistance.     

Modular epistasis in yeast metabolism

Epistatic interactions, manifested in the effects of mutations on the phenotypes caused by other mutations, may help uncover the functional organization of complex biological networks. Here, we studied system-level epistatic interactions by computing growth phenotypes of all single and double knockouts of 890 metabolic genes in Saccharomyces cerevisiae, using the framework of flux balance analysis. A new scale for epistasis identified a distinctive trimodal distribution of these epistatic effects, allowing gene pairs to be classified as buffering, aggravating or noninteracting. We found that the ensuing epistatic interaction network could be organized hierarchically into function-enriched modules that interact with each other 'monochromatically' (i.e., with purely aggravating or purely buffering epistatic links). This property extends the concept of epistasis from single genes to functional units and provides a new definition of biological modularity, which emphasizes interactions between, rather than within, functional modules. Our approach can be used to infer functional gene modules from purely phenotypic epistasis measurements.     

Construction of a mouse yeast artificial chromosome library in a recombination-deficient strain of yeast.

We have constructed a new generation yeast artificial chromosome (YAC) library from female C57BL/10 mice in a recombination-deficient strain of Saccharomyces cerevisiae carrying a mutation in the RAD52 gene. The YAC library contains 41,568 clones with an average insert size of 240 kilobases, representing a greater than threefold coverage of the mouse genome. Currently, the library can be screened by polymerase chain reaction and we have isolated positive clones at a number of loci in the mouse genome. This rad52 library should enable a long-term assessment of the effect of one of the yeast recombination pathway genes on both, genome-wide YAC clone stability and the frequency of chimaeric clones.     

Genomic and functional evolution of the Drosophila melanogaster sperm proteome

In addition to delivering a haploid genome to the egg, sperm have additional critical functions, including egg activation, origination of the zygote centrosome and delivery of paternal factors. Despite this, existing knowledge of the molecular basis of sperm form and function is limited. We used whole-sperm mass spectrometry to identify 381 proteins of the Drosophila melanogaster sperm proteome (DmSP). This approach identified mitochondrial, metabolic and cytoskeletal proteins, in addition to several new functional categories. We also observed nonrandom genomic clustering of sperm genes and underrepresentation on the X chromosome. Identification of widespread functional constraint on the proteome indicates that sexual selection has had a limited role in the overall evolution of D. melanogaster sperm. The relevance of the DmSP to the study of mammalian sperm function and fertilization mechanisms is demonstrated by the identification of substantial homology between the DmSP and proteins of the mouse axoneme accessory structure.     

Robustness against mutations in genetic networks of yeast

There are two principal mechanisms that are responsible for the ability of an organism's physiological and developmental processes to compensate for mutations. In the first, genes have overlapping functions, and loss-of-function mutations in one gene will have little phenotypic effect if there are one or more additional genes with similar functions. The second mechanism has its origin in interactions between genes with unrelated functions, and has been documented in metabolic and regulatory gene networks. Here I analyse, on a genome-wide scale, which of these mechanisms of robustness against mutations is more prevalent. I used functional genomics data from the yeast Saccharomyces cerevisiae to test hypotheses related to the following: if gene duplications are mostly responsible for robustness, then a correlation is expected between the similarity of two duplicated genes and the effect of mutations in one of these genes. My results demonstrate that interactions among unrelated genes are the major cause of robustness against mutations. This type of robustness is probably an evolved response of genetic networks to stabilizing selection.