Genome-wide genetic association of complex traits in heterogeneous stock mice

Difficulties in fine-mapping quantitative trait loci (QTLs) are a major impediment to progress in the molecular dissection of complex traits in mice. Here we show that genome-wide high-resolution mapping of multiple phenotypes can be achieved using a stock of genetically heterogeneous mice. We developed a conservative and robust bootstrap analysis to map 843 QTLs with an average 95% confidence interval of 2.8 Mb. The QTLs contribute to variation in 97 traits, including models of human disease (asthma, type 2 diabetes mellitus, obesity and anxiety) as well as immunological, biochemical and hematological phenotypes. The genetic architecture of almost all phenotypes was complex, with many loci each contributing a small proportion to the total variance. Our data set, freely available at http://gscan.well.ox.ac.uk, provides an entry point to the functional characterization of genes involved in many complex traits.     

Molecular evolution of the human interleukin-8 receptor gene cluster.

Interleukin-8 (IL-8) is the prototype for a family of at least eight neutrophil chemoattractants whose genes map to human chromosome 4q13-q21. Two human IL-8 receptors, IL8RA and IL8RB, are known from cDNA cloning; IL8RA is a promiscuous receptor for at least two other related ligands, GRO alpha and NAP-2. We now report cloning of the genes for IL8RA, IL8RB and a recently inactivated pseudogene of receptor A (IL8RAP). These form a cluster of only three genes in the superfamily of G protein-coupled receptors (GPCRs) and map to 2q34-q35. The coevolutionary diversity displayed by the IL-8 ligand-receptor complex--ligand promiscuity for IL-8, receptor promiscuity for IL8RA, gene duplication for both ligands and receptors and gene extinction in the case of IL8RAP--is unprecedented for the GPCR superfamily.     

Analysing complex genetic traits with chromosome substitution strains

Many valuable animal models of human disease are known and new models are continually being generated in existing inbred strains,. Some disease models are simple mendelian traits, but most have a polygenic basis. The current approach to identifying quantitative trait loci (QTLs) that underlie such traits is to localize them in crosses, construct congenic strains carrying individual QTLs, and finally map and clone the genes. This process is time-consuming and expensive, requiring the genotyping of large crosses and many generations of breeding. Here we describe a different approach in which a panel of chromosome substitution strains (CSSs) is used for QTL mapping. Each of these strains has a single chromosome from the donor strain substituting for the corresponding chromosome in the host strain. We discuss the construction, applications and advantages of CSSs compared with conventional crosses for detecting and analysing QTLs, including those that have weak phenotypic effects.     

Dopa decarboxylase (Ddc) affects variation in Drosophila longevity

Mutational analyses in model organisms have shown that genes affecting metabolism and stress resistance regulate life span, but the genes responsible for variation in longevity in natural populations are largely unidentified. Previously, we mapped quantitative trait loci (QTLs) affecting variation in longevity between two Drosophila melanogaster strains. Here, we show that the longevity QTL in the 36E;38B cytogenetic interval on chromosome 2 contains multiple closely linked QTLs, including the Dopa decarboxylase (Ddc) locus. Complementation tests to mutations show that Ddc is a positional candidate gene for life span in these strains. Linkage disequilibrium (LD) mapping in a sample of 173 alleles from a single population shows that three common molecular polymorphisms in Ddc account for 15.5% of the genetic contribution to variance in life span from chromosome 2. The polymorphisms are in strong LD, and the effects of the haplotypes on longevity suggest that the polymorphisms are maintained by balancing selection. DDC catalyzes the final step in the synthesis of the neurotransmitters, dopamine and serotonin. Thus, these data implicate variation in the synthesis of bioamines as a factor contributing to natural variation in individual life span.     

Variants modulating the expression of a chromosome domain encompassing PLAG1 influence bovine stature

We report mapping of a quantitative trait locus (QTL) with a major effect on bovine stature to a ～780-kb interval using a Hidden Markov Model-based approach that simultaneously exploits linkage and linkage disequilibrium. We re-sequenced the interval in six sires with known QTL genotype and identified 13 clustered candidate quantitative trait nucleotides (QTNs) out of >9,572 discovered variants. We eliminated five candidate QTNs by studying the phenotypic effect of a recombinant haplotype identified in a breed diversity panel. We show that the QTL influences fetal expression of seven of the nine genes mapping to the ～780-kb interval. We further show that two of the eight candidate QTNs, mapping to the PLAG1-CHCHD7 intergenic region, influence bidirectional promoter strength and affect binding of nuclear factors. By performing expression QTL analyses, we identified a splice site variant in CHCHD7 and exploited this naturally occurring null allele to exclude CHCHD7 as single causative gene.     

Fine mapping of regulatory loci for mammalian gene expression using radiation hybrids

We mapped regulatory loci for nearly all protein-coding genes in mammals using comparative genomic hybridization and expression array measurements from a panel of mouse-hamster radiation hybrid cell lines. The large number of breaks in the mouse chromosomes and the dense genotyping of the panel allowed extremely sharp mapping of loci. As the regulatory loci result from extra gene dosage, we call them copy number expression quantitative trait loci, or ceQTLs. The -2log10P support interval for the ceQTLs was <150 kb, containing an average of <2-3 genes. We identified 29,769 trans ceQTLs with -log10P > 4, including 13 hotspots each regulating >100 genes in trans. Further, this work identifies 2,761 trans ceQTLs harboring no known genes, and provides evidence for a mode of gene expression autoregulation specific to the X chromosome.     

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.     

Integrated transcriptional profiling and linkage analysis for identification of genes underlying disease

Integration of genome-wide expression profiling with linkage analysis is a new approach to identifying genes underlying complex traits. We applied this approach to the regulation of gene expression in the BXH/HXB panel of rat recombinant inbred strains, one of the largest available rodent recombinant inbred panels and a leading resource for genetic analysis of the highly prevalent metabolic syndrome. In two tissues important to the pathogenesis of the metabolic syndrome, we mapped cis- and trans-regulatory control elements for expression of thousands of genes across the genome. Many of the most highly linked expression quantitative trait loci are regulated in cis, are inherited essentially as monogenic traits and are good candidate genes for previously mapped physiological quantitative trait loci in the rat. By comparative mapping we generated a data set of 73 candidate genes for hypertension that merit testing in human populations. Mining of this publicly available data set is expected to lead to new insights into the genes and regulatory pathways underlying the extensive range of metabolic and cardiovascular disease phenotypes that segregate in these recombinant inbred strains.     

A systematic, genome-wide, phenotype-driven mutagenesis programme for gene function studies in the mouse

As the human genome project approaches completion, the challenge for mammalian geneticists is to develop approaches for the systematic determination of mammalian gene function. Mouse mutagenesis will be a key element of studies of gene function. Phenotype-driven approaches using the chemical mutagen ethylnitrosourea (ENU) represent a potentially efficient route for the generation of large numbers of mutant mice that can be screened for novel phenotypes. The advantage of this approach is that, in assessing gene function, no a priori assumptions are made about the genes involved in any pathway. Phenotype-driven mutagenesis is thus an effective method for the identification of novel genes and pathways. We have undertaken a genome-wide, phenotype-driven screen for dominant mutations in the mouse. We generated and screened over 26,000 mice, and recovered some 500 new mouse mutants. Our work, along with the programme reported in the accompanying paper, has led to a substantial increase in the mouse mutant resource and represents a first step towards systematic studies of gene function in mammalian genetics.     

A humanized model for multiple sclerosis using HLA-DR2 and a human T-cell receptor

Multiple sclerosis (MS) is a complex chronic neurologic disease with a suspected autoimmune pathogenesis. Although there is evidence that the development of MS is determined by both environmental influences and genes, these factors are largely undefined, except for major histocompatibility (MHC) genes. Linkage analyses and association studies have shown that susceptibility to MS is associated with genes in the human histocompatibility leukocyte antigens (HLA) class II region, but the contribution of these genes to MS disease development less compared with their contribution to disorders such as insulin-dependent diabetes mellitus. Due to the strong linkage disequilibrium in the MHC class II region, it has not been possible to determine which gene(s) is responsible for the genetic predisposition. In transgenic mice, we have expressed three human components involved in T-cell recognition of an MS-relevant autoantigen presented by the HLA-DR2 molecule: DRA*0101/DRB1*1501 (HLA-DR2), an MHC class II candidate MS susceptibility genes found in individuals of European descent; a T-cell receptor (TCR) from an MS-patient-derived T-cell clone specific for the HLA-DR2 bound immunodominant myelin basic protein (MBP) 4102 peptide; and the human CD4 coreceptor. The amino acid sequence of the MBP 84-102 peptide is the same in both human and mouse MBP. Following administration of the MBP peptide, together with adjuvant and pertussis toxin, transgenic mice developed focal CNS inflammation and demyelination that led to clinical manifestations and disease courses resembling those seen in MS. Spontaneous disease was observed in 4% of mice. When DR2 and TCR double-transgenic mice were backcrossed twice to Rag2 (for recombination-activating gene 2)-deficient mice, the incidence of spontaneous disease increased, demonstrating that T cells specific for the HLA-DR2 bound MBP peptide are sufficient and necessary for development of disease. Our study provides evidence that HLA-DR2 can mediate both induced and spontaneous disease resembling MS by presenting an MBP self-peptide to T cells.     

The gene encoding proline dehydrogenase modulates sensorimotor gating in mice

Hemizygous cryptic deletions of the q11 band of human chromosome 22 have been associated with a number of psychiatric and behavioural phenotypes, including schizophrenia. Here we report the isolation and characterization of PRODH, a human homologue of Drosophila melanogaster sluggish-A (slgA), which encodes proline dehydrogenase responsible for the behavioural phenotype of the slgA mutant. PRODH is localized at chromosome 22q11 in a region deleted in some psychiatric patients. We also isolated the mouse homologue of slgA (Prodh), identified a mutation in this gene in the Pro/Re hyperprolinaemic mouse strain and found that these mice have a deficit in sensorimotor gating accompanied by regional neurochemical alterations in the brain. Sensorimotor gating is a neural filtering process that allows attention to be focused on a given stimulus, and is affected in patients with neuropsychiatric disorders. Furthermore, several lines of evidence suggest that proline may serve as a modulator of synaptic transmission in the mammalian brain. Our observations, in conjunction with the chromosomal location of PRODH, suggest a potential involvement of this gene in the 22q11-associated psychiatric and behavioural phenotypes.     

Maternal imprinting of the mouse Snrpn gene and conserved linkage homology with the human Prader-Willi syndrome region.

Prader-Willi syndrome (PWS) is associated with paternal gene deficiencies in human chromosome 15q11-13, suggesting that PWS is caused by a deficiency in one or more maternally imprinted genes. We have now mapped a gene, Snrpn, encoding a brain-enriched small nuclear ribonucleoprotein (snRNP)-associated polypeptide SmN, to mouse chromosome 7 in a region of homology with human chromosome 15q11-13 and demonstrated that Snrpn is a maternally imprinted gene in mouse. These studies, in combination with the accompanying human mapping studies showing that SNRPN maps in the Prader-Willi critical region, identify SNRPN as a candidate gene involved in PWS and suggest that PWS may be caused, in part, by defects in mRNA processing.     

Candidate lung tumor susceptibility genes identified through whole-genome association analyses in inbred mice

We performed a whole-genome association analysis of lung tumor susceptibility using dense SNP maps ( approximately 1 SNP per 20 kb) in inbred mice. We reproduced the pulmonary adenoma susceptibility 1 (Pas1) locus identified in previous linkage studies and further narrowed this quantitative trait locus (QTL) to a region of less than 0.5 Mb in which at least two genes, Kras2 (Kirsten rat sarcoma oncogene 2) and Casc1 (cancer susceptibility candidate 1; also known as Las1), are strong candidates. Casc1 knockout mouse tumor bioassays showed that Casc1-deficient mice were susceptible to chemical induction of lung tumors. We also found three more genetic loci for lung adenoma development. Analysis of one of these candidate loci identified a previously uncharacterized gene Lasc1, bearing a nonsynonymous substitution (D102E). We found that the Lasc1 Glu102 allele preferentially promotes lung tumor cell growth. Our findings demonstrate the prospects for using dense SNP maps in laboratory mice to refine previous QTL regions and identify genetic determinants of complex traits.     

High-throughput retroviral tagging to identify components of specific signaling pathways in cancer

Genetic screens carried out in lower organisms such as yeast, Drosophila melanogaster and Caenorhabditis elegans have revealed many signaling pathways. For example, components of the RAS signaling cascade were identified using a mutant eye phenotype in D. melanogaster as a readout. Screening is usually based on enhancing or suppressing a phenotype by way of a known mutation in a particular signaling pathway. Such in vivo screens have been difficult to carry out in mammals, however, owing to their relatively long generation times and the limited number of animals that can be screened. Here we describe an in vivo mammalian genetic screen used to identify components of pathways contributing to oncogenic transformation. We applied retroviral insertional mutagenesis in Myc transgenic (E mu Myc) mice lacking expression of Pim1 and Pim2 to search for genes that can substitute for Pim1 and Pim2 in lymphomagenesis. We determined the chromosomal positions of 477 retroviral insertion sites (RISs) derived from 38 tumors from E mu Myc Pim1(-/-) Pim2(-/-) mice and 27 tumors from E mu Myc control mice using the Ensembl and Celera annotated mouse genome databases. There were 52 sites occupied by proviruses in more than one tumor. These common insertion sites (CISs) are likely to contain genes contributing to tumorigenesis. Comparison of the RISs in tumors of Pim-null mice with the RISs in tumors of E mu Myc control mice indicated that 10 of the 52 CISs belong to the Pim complementation group. In addition, we found that Pim3 is selectively activated in Pim-null tumor cells, which supports the validity of our approach.     

Positional cloning of a novel gene influencing asthma from chromosome 2q14

Asthma is a common disease in children and young adults. Four separate reports have linked asthma and related phenotypes to an ill-defined interval between 2q14 and 2q32 (refs. 1-4), and two mouse genome screens have linked bronchial hyper-responsiveness to the region homologous to 2q14 (refs. 5,6). We found and replicated association between asthma and the D2S308 microsatellite, 800 kb distal to the IL1 cluster on 2q14. We sequenced the surrounding region and constructed a comprehensive, high-density, single-nucleotide polymorphism (SNP) linkage disequilibrium (LD) map. SNP association was limited to the initial exons of a solitary gene of 3.6 kb (DPP10), which extends over 1 Mb of genomic DNA. DPP10 encodes a homolog of dipeptidyl peptidases (DPPs) that cleave terminal dipeptides from cytokines and chemokines, and it presents a potential new target for asthma therapy.