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.     

Ru2 and Ru encode mouse orthologs of the genes mutated in human Hermansky-Pudlak syndrome types 5 and 6

Hermansky-Pudlak syndrome (HPS) is a genetically heterogeneous disease involving abnormalities of melanosomes, platelet dense granules and lysosomes. Here we have used positional candidate and transgenic rescue approaches to identify the genes mutated in ruby-eye 2 and ruby-eye mice (ru2 and ru, respectively), two 'mimic' mouse models of HPS. We also show that these genes are orthologs of the genes mutated in individuals with HPS types 5 and 6, respectively, and that their protein products directly interact. Both genes are previously unknown and are found only in higher eukaryotes, and together represent a new class of genes that have evolved in higher organisms to govern the synthesis of highly specialized lysosome-related organelles.     

Estimating coverage and power for genetic association studies using near-complete variation data

Although studies suggest that SNPs derived from HapMap provide promising coverage and power for association studies, the lack of alternative variation datasets limits independent analysis. Using near-complete variation data for 76 genes resequenced in HapMap samples, we find that coverage of common variation by commercial genotyping arrays is substantially lower compared to the HapMap-based estimates. We quantify the power offered by these arrays for a range of disease models.     

Deep resequencing of GWAS loci identifies independent rare variants associated with inflammatory bowel disease

More than 1,000 susceptibility loci have been identified through genome-wide association studies (GWAS) of common variants; however, the specific genes and full allelic spectrum of causal variants underlying these findings have not yet been defined. Here we used pooled next-generation sequencing to study 56 genes from regions associated with Crohn's disease in 350 cases and 350 controls. Through follow-up genotyping of 70 rare and low-frequency protein-altering variants in nine independent case-control series (16,054 Crohn's disease cases, 12,153 ulcerative colitis cases and 17,575 healthy controls), we identified four additional independent risk factors in NOD2, two additional protective variants in IL23R, a highly significant association with a protective splice variant in CARD9 (P < 1 × 10(-16), odds ratio ≈ 0.29) and additional associations with coding variants in IL18RAP, CUL2, C1orf106, PTPN22 and MUC19. We extend the results of successful GWAS by identifying new, rare and probably functional variants that could aid functional experiments and predictive models.     

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 genetically isolated populations of Finland and sardinia may not be a panacea for linkage disequilibrium mapping of common disease genes

The choice of which population to study in the mapping of common disease genes may be critical. Isolated founder populations, such as that found in Finland, have already proved extremely useful for mapping the genes for specific rare monogenic disorders and are being used in attempts to map the genes underlying common, complex diseases. But simulation results suggest that, under the common disease-common variant hypothesis, most isolated populations will prove no more useful for linkage disequilibrium (LD) mapping of common disease genes than large outbred populations. There is very little empirical data to either support or refute this conclusion at present. Therefore, we evaluated LD between 21 common microsatellite polymorphisms on chromosome 18q21 in 2 genetic isolates (Finland and Sardinia) and compared the results with those observed in two mixed populations (United Kingdom and United States of America). Mean levels of LD were similar across all four populations. Our results provide empirical support for the expectation that genetic isolates like Finland and Sardinia will not prove significantly more valuable than general populations for LD mapping of common variants underlying complex disease.     

Leukaemia disease genes: large-scale cloning and pathway predictions

Retroviral insertional mutagenesis in BXH2 and AKXD recombinant inbred mice induces a high incidence of myeloid or B- and T-cell leukaemia and the proviral integration sites in the leukaemias provide powerful genetic tags for disease gene identification. Some of the disease genes identified by proviral tagging are also associated with human disease, validating this approach for human disease gene identification. Although many leukaemia disease genes have been identified over the years, many more remain to be cloned. Here we describe an inverse PCR (IPCR) method for proviral tagging that makes use of automated DNA sequencing and the genetic tools provided by the Mouse Genome Project, which increases the throughput for disease gene identification. We also use this IPCR method to clone and analyse more than 400 proviral integration sites from AKXD and BXH2 leukaemias and, in the process, identify more than 90 candidate disease genes. Some of these genes function in pathways already implicated in leukaemia, whereas others are likely to define new disease pathways. Our studies underscore the power of the mouse as a tool for gene discovery and functional genomics.     

Mutations of PKD1 in ADPKD2 cysts suggest a pathogenic effect of trans-heterozygous mutations

Autosomal dominant polycystic kidney disease (ADPKD) is caused by mutations in PKD1 and PKD2. The products of these genes associate to form heteromeric complexes. Several models have been proposed to explain the mechanism of cyst formation. Here we find somatic mutations of PKD2 in 71% of ADPKD2 cysts analysed. Clonal somatic mutations of PKD1 were identified in a subset of cysts that lacked PKD2 mutations.     

Linkage of high-density lipoprotein-cholesterol concentrations to a locus on chromosome 9p in Mexican Americans.

High-density lipoproteins (HDLs) are anti-atherogenic lipoproteins that have a major role in transporting cholesterol from peripheral tissues to the liver, where it is removed. Epidemiologic studies have shown that low levels of high-density lipoprotein-cholesterol (HDL-C) are associated with an increased incidence of coronary heart disease and an increased mortality rate, indicating a protective role of high concentrations of HDL-C against atherogenesis and the development of coronary heart disease. HDL-C level is influenced by several genetic and nongenetic factors. Nongenetic factors include smoking, which has been shown to decrease the HDL-C level. Exercise and alcohol have been shown to increase HDL-C levels. Decreased HDL-C is often associated with other coronary heart disease risk factors such as obesity, hyperinsulinemia and insulin resistance, hypertriglyceridemia and hypertension. Although several genes have been identified for rare forms of dyslipidemia, the genes accounting for major variation in HDL-C levels have yet to be identified. Using a multipoint variance components linkage approach, we found strong evidence of linkage (lod score=3.4; P=0.00004) of a quantitative trait locus (QTL) for HDL-C level to a genetic location between markers D9S925 and D9S741 on chromosome 9p in Mexican Americans. A replication study in an independent set of Mexican American families confirmed the existence of a QTL on chromosome 9p.     

Detection of regulatory variation in mouse genes

Functional polymorphism in genes can be classified as coding variation, altering the amino-acid sequence of the encoded protein, or regulatory variation, affecting the level or pattern of expression of the gene. Coding variation can be recognized directly from DNA sequence, and consequently its frequency and characteristics have been extensively described. By contrast, virtually nothing is known about the extent to which gene regulation varies in populations. Yet it is likely that regulatory variants are important in modulating gene function: alterations in gene regulation have been proposed to influence disease susceptibility and to have been the primary substrate for the evolution of species. Here, we report a systematic study to assess the extent of cis-acting regulatory variation in 69 genes across four inbred mouse strains. We find that at least four of these genes show allelic differences in expression level of 1.5-fold or greater, and that some of these differences are tissue specific. The results show that the impact of regulatory variants can be detected at a significant frequency in a genomic survey and suggest that such variation may have important consequences for organismal phenotype and evolution. The results indicate that larger-scale surveys in both mouse and human could identify a substantial number of genes with common regulatory variation.     

Genetic dissection of autoimmune type I diabetes in the BB rat.

The BB rat is among the best models of insulin-dependent diabetes mellitus--with onset and pathogenesis closely resembling the human disease. One unusual feature is a severe T-cell lymphopenia, which appears to be inherited as a recessive trait controlled by a single gene, Lyp. Based on genetic analysis of several crosses, we show that development of diabetes involves at least three genes: Lyp, which is tightly linked to the neuropeptide Y (Npy) gene on chromosome 4, a gene linked to the major histocompatibility complex (MHC) on chromosome 20, and a third unmapped gene for which the Fischer rat strain carries an allele conferring resistance.     

The Collaborative Cross, a community resource for the genetic analysis of complex traits

The goal of the Complex Trait Consortium is to promote the development of resources that can be used to understand, treat and ultimately prevent pervasive human diseases. Existing and proposed mouse resources that are optimized to study the actions of isolated genetic loci on a fixed background are less effective for studying intact polygenic networks and interactions among genes, environments, pathogens and other factors. The Collaborative Cross will provide a common reference panel specifically designed for the integrative analysis of complex systems and will change the way we approach human health and disease.     

Linkage disequilibrium mapping in isolated founder populations: diastrophic dysplasia in Finland.

Linkage disequilibrium mapping in isolated populations provides a powerful tool for fine structure localization of disease genes. Here, Luria and Delbrück's classical methods for analysing bacterial cultures are adapted to the study of human isolated founder populations in order to estimate (i) the recombination fraction between a disease locus and a marker; (ii) the expected degree of allelic homogeneity in a population; and (iii) the mutation rate of marker loci. Using these methods, we report striking linkage disequilibrium for diastrophic dysplasia (DTD) in Finland indicating that the DTD gene should lie within 0.06 centimorgans (or about 60 kilobases) of the CSF1R gene. Predictions about allelic homogeneity in Finland and mutation rates in simple sequence repeats are confirmed by independent observations.     

Myosin IXB variant increases the risk of celiac disease and points toward a primary intestinal barrier defect

Celiac disease is probably the best-understood immune-related disorder. The disease presents in the small intestine and results from the interplay between multiple genes and gluten, the triggering environmental factor. Although HLA class II genes explain 40% of the heritable risk, non-HLA genes accounting for most of the familial clustering have not yet been identified. Here we report significant and replicable association (P = 2.1 x 10(-6)) to a common variant located in intron 28 of the gene myosin IXB (MYO9B), which encodes an unconventional myosin molecule that has a role in actin remodeling of epithelial enterocytes. Individuals homozygous with respect to the at-risk allele have a 2.3-times higher risk of celiac disease (P = 1.55 x 10(-5)). This result is suggestive of a primary impairment of the intestinal barrier in the etiology of celiac disease, which may explain why immunogenic gluten peptides are able to pass through the epithelial barrier.     

Systematic mapping of genetic interactions in Caenorhabditis elegans identifies common modifiers of diverse signaling pathways

Most heritable traits, including disease susceptibility, are affected by interactions between multiple genes. However, we understand little about how genes interact because very few possible genetic interactions have been explored experimentally. We have used RNA interference in Caenorhabditis elegans to systematically test approximately 65,000 pairs of genes for their ability to interact genetically. We identify approximately 350 genetic interactions between genes functioning in signaling pathways that are mutated in human diseases, including components of the EGF/Ras, Notch and Wnt pathways. Most notably, we identify a class of highly connected 'hub' genes: inactivation of these genes can enhance the phenotypic consequences of mutation of many different genes. These hub genes all encode chromatin regulators, and their activity as genetic hubs seems to be conserved across animals. We propose that these genes function as general buffers of genetic variation and that these hub genes may act as modifier genes in multiple, mechanistically unrelated genetic diseases in humans.     

Common variation in three genes, including a noncoding variant in CFH, strongly influences risk of age-related macular degeneration

Age-related macular degeneration (AMD) is a common, late-onset disease with seemingly typical complexity: recurrence ratios for siblings of an affected individual are three- to sixfold higher than in the general population, and family-based analysis has resulted in only modestly significant evidence for linkage. In a case-control study drawn from a US-based population of European descent, we have identified a previously unrecognized common, noncoding variant in CFH, the gene encoding complement factor H, that substantially increases the influence of this locus on AMD, and we have strongly replicated the associations of four other previously reported common alleles in three genes (P values ranging from 10(-6) to 10(-70)). Despite excellent power to detect epistasis, we observed purely additive accumulation of risk from alleles at these genes. We found no differences in association of these loci with major phenotypic categories of advanced AMD. Genotypes at these five common SNPs define a broad spectrum of interindividual disease risk and explain about half of the classical sibling risk of AMD in our study population.     

In vivo RNA interference demonstrates a role for Nramp1 in modifying susceptibility to type 1 diabetes

Type 1 diabetes is an autoimmune disease influenced by multiple genetic loci. Although more than 20 insulin-dependent diabetes (Idd) loci have been implicated in the nonobese diabetic (NOD) mouse model, few causal gene variants have been identified. Here we show that RNA interference (RNAi) can be used to probe candidate genes in this disease model. Slc11a1 encodes a phagosomal ion transporter, Nramp1, that affects resistance to intracellular pathogens and influences antigen presentation. This gene is the strongest candidate among the 42 genes in the Idd5.2 region; a naturally occurring mutation in the protective Idd5.2 haplotype results in loss of function of the Nramp1 protein. Using lentiviral transgenesis, we generated NOD mice in which Slc11a1 is silenced by RNAi. Silencing reduced the frequency of type 1 diabetes, mimicking the protective Idd5.2 region. Our results demonstrate a role for Slc11a1 in modifying susceptibility to type 1 diabetes and illustrate that RNAi can be used to study causal genes in a mammalian model organism.     

Biased biological functions of horizontally transferred genes in prokaryotic genomes

Horizontal gene transfer is one of the main mechanisms contributing to microbial genome diversification. To clarify the overall picture of interspecific gene flow among prokaryotes, we developed a new method for detecting horizontally transferred genes and their possible donors by Bayesian inference with training models for nucleotide composition. Our method gives the average posterior probability (horizontal transfer index) for each gene sequence, with a low horizontal transfer index indicating recent horizontal transfer. We found that 14% of open reading frames in 116 prokaryotic complete genomes were subjected to recent horizontal transfer. Based on this data set, we quantitatively determined that the biological functions of horizontally transferred genes, except mobile element genes, are biased to three categories: cell surface, DNA binding and pathogenicity-related functions. Thus, the transferability of genes seems to depend heavily on their functions.     

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.