Integration of Mouse and Human Genome-Wide Association Data Identifies KCNIP4 as an Asthma Gene

Asthma is a common chronic respiratory disease characterized by airway hyperresponsiveness (AHR). The genetics of asthma have been widely studied in mouse and human, and homologous genomic regions have been associated with mouse AHR and human asthma-related phenotypes. Our goal was to identify asthma-related genes by integrating AHR associations in mouse with human genome-wide association study (GWAS) data. We used Efficient Mixed Model Association (EMMA) analysis to conduct a GWAS of baseline AHR measures from males and females of 31 mouse strains. Genes near or containing SNPs with EMMA p-values <0.001 were selected for further study in human GWAS. The results of the previously reported EVE consortium asthma GWAS meta-analysis consisting of 12,958 diverse North American subjects from 9 study centers were used to select a subset of homologous genes with evidence of association with asthma in humans. Following validation attempts in three human asthma GWAS (i.e., Sepracor/LOCCS/LODO/Illumina, GABRIEL, DAG) and two human AHR GWAS (i.e., SHARP, DAG), the Kv channel interacting protein 4 (KCNIP4) gene was identified as nominally associated with both asthma and AHR at a gene- and SNP-level. In EVE, the smallest KCNIP4 association was at rs6833065 (P-value 2.9e-04), while the strongest associations for Sepracor/LOCCS/LODO/Illumina, GABRIEL, DAG were 1.5e-03, 1.0e-03, 3.1e-03 at rs7664617, rs4697177, rs4696975, respectively. At a SNP level, the strongest association across all asthma GWAS was at rs4697177 (P-value 1.1e-04). The smallest P-values for association with AHR were 2.3e-03 at rs11947661 in SHARP and 2.1e-03 at rs402802 in DAG. Functional studies are required to validate the potential involvement of KCNIP4 in modulating asthma susceptibility and/or AHR. Our results suggest that a useful approach to identify genes associated with human asthma is to leverage mouse AHR association data.     

Detection of loci for allergic asthma using SMXA recombinant inbred strains of mice

Asthma is regarded as a multifactorial inflammatory disorder arising as a result of inappropriate immune responses in genetically susceptible individuals to common environmental antigens. However, the precise molecular basis is unknown. To identify genes for susceptibility to three asthma-related traits, airway hyperresponsiveness (AHR), eosinophil infiltration, and allergen-specific serum IgE levels, we conducted a genetic analysis using SMXA recombinant inbred (RI) strains of mice. Quantitative trait locus analysis detected a significant locus for AHR on chromosome 17. For eosinophil infiltration, significant loci were detected on chromosomes 9 and 16. Although we could not detect any significant loci for allergen-specific serum IgE, analysis of consomic strains showed that chromosomes 17 and 19 carried genes that affected this trait. We detected genetic susceptibility loci that separately regulated the three asthma-related phenotypes. Our results suggested that different genetic mechanisms regulate these asthma-related phenotypes. Genetic analyses using murine RI and consomic strains enhance understanding of the molecular mechanisms of asthma in human.     

Host genetics and Chlamydia disease: prediction and validation of disease severity mechanisms

Genetic mapping studies may provide association between sequence variants and disease susceptibility that can, with further experimental and computational analysis, lead to discovery of causal mechanisms and effective intervention. We have previously demonstrated that polymorphisms in immunity-related GTPases (IRG) confer a significant difference in susceptibility to Chlamydia psittaci infection in BXD recombinant mice. Here we combine genetic mapping and network modeling to identify causal pathways underlying this association. We infected a large panel of BXD strains with C. psittaci and assessed host genotype, IRG protein polymorphisms, pathogen load, expression of 32 cytokines, inflammatory cell populations, and weight change. Proinflammatory cytokines correlated with each other and were controlled by a novel genetic locus on chromosome 1, but did not affect disease status, as quantified by weight change 6 days after infection In contrast, weight change correlated strongly with levels of inflammatory cell populations and pathogen load that were controlled by an IRG encoding genetic locus (Ctrq3) on chromosome 11. These data provided content to generate a predictive model of infection using a Bayesian framework incorporating genotypes, immune system parameters, and weight change as a measure of disease severity. Two predictions derived from the model were tested and confirmed in a second round of experiments. First, strains with the susceptible IRG haplotype lost weight as a function of pathogen load whereas strains with the resistant haplotype were almost completely unaffected over a very wide range of pathogen load. Second, we predicted that macrophage activation by Ctrq3 would be central in conferring pathogen tolerance. We demonstrated that macrophage depletion in strains with the resistant haplotype led to neutrophil influx and greater weight loss despite a lower pathogen burden. Our results show that genetic mapping and network modeling can be combined to identify causal pathways underlying chlamydial disease susceptibility.     

Haplotype-based genetics in mice and rats

Haplotype blocks are conceptually defined as genomic segments harbouring sets of coupled polymorphisms that reflect a common ancestral origin. Experimentally, however, haplotype blocks are characterized using computational algorithms based on incomplete inventories of polymorphisms. Haplotype blocks and their deduced strain-distribution patterns are considered to be extremely powerful for use in genetic association and mapping experiments in laboratory mice and rats. However, recent high-density SNP screening in commonly used mouse inbred strains reveals a complex pattern, suggesting that the current expectations for the use of haplotype blocks in genetic mapping will have to be revisited.     

Genetic Analysis of Intracapillary Glomerular Lipoprotein Deposits in Aging Mice

Background

Renal aging is characterized by functional and structural changes like decreased glomerular filtration rate, and glomerular, tubular and interstitial damage. To gain insight in pathways involved in renal aging, we studied aged mouse strains and used genetic analysis to identify genes associated with aging phenotypes.

Methods

Upon morphological screening in kidneys from 20-month-old mice from 26 inbred strains we noted intracapillary PAS-positive deposits. The severity of these deposits was quantified by scoring of a total of 50 glomeruli per section (grade 0–4). Electron microscopy and immunohistochemical staining for apoE, apoB, apoA-IV and perilipin-2 was performed to further characterize the lesions. To identify loci associated with these PAS-positive intracapillary glomerular deposits, we performed haplotype association mapping.

Results

Six out of 26 mouse strains showed glomerular PAS-positive deposits. The severity of these deposits varied: NOD(0.97), NZW(0.41), NON(0.30), B10(0.21), C3 H(0.9) and C57BR(0.7). The intracapillary deposits were strongly positive for apoE and weakly positive for apoB and apoA-IV. Haplotype association mapping showed a strong association with a 30-Kb haplotype block on Chr 1 within the Esrrg gene. We investigated 1 Mb on each site of this region, which includes the genes Spata17, Gpatch2, Esrrg, Ush2a and Kctd3.

Conclusions

By analyzing 26 aged mouse strains we found that some strains developed an intracapillary PAS and apoE-positive lesion and identified a small haplotype block on Chr 1 within the Esrrg gene to be associated with these lipoprotein deposits. The region spanning this haplotype block contains the genes Spata17, Gpatch2, Esrrg, Ush2a and Kctd3, which are all highly expressed in the kidney. Esrrg might be involved in the evolvement of these glomerular deposits by influencing lipid metabolism and possibly immune reponses.     

An Integrative Genomic Analysis of the Superior Fecundity Phenotype in QSi5 Mice

Laboratory inbred mouse models are a valuable resource to identify quantitative trait loci (QTL) for complex reproductive performance traits. Advances in mouse genomics and high density single nucleotide polymorphism mapping has enabled genome-wide association studies to identify genes linked with specific phenotypes. Gene expression profiles of reproductive tissues also provide potentially useful information for identifying genes that play an important role. We have developed a highly fecund inbred strain, QSi5, with accompanying genotyping for comparative analysis of reproductive performance. Here we analyzed the QSi5 phenotype using a comparative analysis with fecundity data derived from 22 inbred strains of mice from the Mouse Phenome Project, and integration with published expression data from mouse ovary development. Using a haplotype association approach, 400 fecundity-associated regions (FDR < 0.05) with 499 underlying genes were identified. The most significant associations were located on Chromosomes 14, 8, and 6, and the genes underlying these regions were extracted. When these genes were analyzed for expression in an ovarian development profile (GSE6916) several distinctive co-expression patterns across each developmental stage were identified. The genetic analysis also refined 21 fecundity associated intervals on Chromosomes 1, 6, 9, 13, and 17 that overlapped with previously reported reproductive performance QTL. The combined use of phenotypic and in silico data with an integrative genomic analysis provides a powerful tool for elucidating the molecular mechanisms underlying fecundity.     

Identification of candidate genes for nicotine withdrawal in C57BL/6J × DBA/2J recombinant inbred mice

Nicotine is the reinforcing ingredient in tobacco. Following chronic exposure, sudden cessation of nicotine use produces negative symptoms of withdrawal that contribute to dependence. The molecular mechanisms underlying nicotine withdrawal behaviors, however, are poorly understood. Using recombinant inbred mice, chronic nicotine was delivered by minipump and withdrawal induced using mecamylamine. Somatic signs of withdrawal, and anxiety-like behavior using elevated plus maze, were then assessed. Interval mapping was used to identify associations between genetic variation and withdrawal behaviors, and with basal gene expression. Differential gene expression following nicotine exposure and withdrawal was also assessed in progenitor mice using microarrays. Quantitative trait loci mapping identified chromosome intervals with significant genetic associations to somatic signs of withdrawal or withdrawal-induced anxiety-like behavior. Using bioinformatics, and association with basal gene expression in nucleus accumbens, we implicated Rb1, Bnip3l, Pnma2, Itm2b, and Kif13b as candidate genes for somatic signs of withdrawal, and Galr1, which showed trans-regulation from a region of chromosome 14 that was associated with somatic signs of withdrawal. Candidate genes within the chromosome 9 region associated with anxiety-like withdrawal behavior included Dixdc1, Ncam1, and Sorl1. Bioinformatics identified six genes that were also significantly associated with nicotine or alcohol traits in recent human genome-wide association studies. Withdrawal-associated somatic signs and anxiety-like behavior had strong non-overlapping genetic associations, respectively, with regions of chromosome 14 and chromosome 9. Genetic, behavioral and gene expression correlations, and bioinformatics analysis identified several candidate genes that may represent novel molecular targets for modulating nicotine withdrawal symptoms.     

In silico QTL mapping of maternal nurturing ability with the mouse diversity panel

Significant variation exists for maternal nurturing ability in inbred mice. Although classical mapping approaches have identified quantitative trait loci (QTL) that may account for this variation, the underlying genes are unknown. In this study, lactation performance data among the mouse diversity panel were used to map genomic regions associated with this variation. Females from each of 32 inbred strains (n = 8-19 dams/strain) were studied during the first 8 days of lactation by allowing them to raise weight- and size-normalized cross-foster litters (10 pups/litter). Average daily weight gain (ADG) of litters served as the primary indicator of milk production. The number of pups successfully reared to 8 days (PNUM8) also served as a related indicator of maternal performance. Initial haplotype association analysis using a Bonferroni-corrected, genome-wide threshold revealed 10 and 15 associations encompassing 11 and 13 genes for ADG and PNUM8, respectively. The most significant of these associated haplotype blocks were found on MMU 8, 11, and 19 and contained the genes Nr3c2, Egfr, Sec61g, and Gnaq. Lastly, two haplotype blocks on MMU9 were detected in association with PNUM8. These overlapped with the previously described maternal performance QTL, Neogq1. These results suggest that the application of in silico QTL mapping is a useful tool in discovering the presence of novel candidate genes involved in determining lactation capacity in mice.     

Linkage disequilibrium and haplotype mapping of a skin cancer susceptibility locus in outbred mice

Car-R (carcinogenesis-resistant) and Car-S (carcinogenesis-susceptible) outbred mice, obtained by phenotypic selection from an initial intercross of eight inbred strains, show a >100-fold difference in their susceptibility to two-stage skin tumorigenesis. We found that the lines carry a high degree of genetic polymorphism, with an average heterozygosity of 0.39. This polymorphism allowed the use of linkage disequilibrium (LD) and haplotype analysis for the mapping of a skin cancer modifier locus on Chr 7, in a short region of 6 cM, around the Tyr gene. Car-S mice inherited the susceptibility allele at this locus from the A/J, BALB/c, SJL/J, and SWR/J strains. Our results demonstrate the feasibility and usefulness of mapping disease genes by LD in phenotypically selected, genetically heterogeneous animals. Received: 16 March 2000 / Accepted: 9 June 2000     

An integrative approach for the identification of quantitative trait loci

The genetic dissection of complex traits is one of the most difficult and most important challenges facing science today. We discuss here an integrative approach to quantitative trait loci (QTL) mapping in mice. This approach makes use of the wealth of genetic tools available in mice, as well as the recent advances in genome sequence data already available for a number of inbred mouse strains. We have developed mapping strategies that allow a stepwise narrowing of a QTL mapping interval, prioritizing candidate genes for further analysis with the potential of identifying the most probable candidate gene for the given trait. This approach integrates traditional mapping tools, fine mapping tools, sequence-based analysis, bioinformatics and gene expression.     

Efficient control of population structure in model organism association mapping

Genomewide association mapping in model organisms such as inbred mouse strains is a promising approach for the identification of risk factors related to human diseases. However, genetic association studies in inbred model organisms are confronted by the problem of complex population structure among strains. This induces inflated false positive rates, which cannot be corrected using standard approaches applied in human association studies such as genomic control or structured association. Recent studies demonstrated that mixed models successfully correct for the genetic relatedness in association mapping in maize and Arabidopsis panel data sets. However, the currently available mixed-model methods suffer from computational inefficiency. In this article, we propose a new method, efficient mixed-model association (EMMA), which corrects for population structure and genetic relatedness in model organism association mapping. Our method takes advantage of the specific nature of the optimization problem in applying mixed models for association mapping, which allows us to substantially increase the computational speed and reliability of the results. We applied EMMA to in silico whole-genome association mapping of inbred mouse strains involving hundreds of thousands of SNPs, in addition to Arabidopsis and maize data sets. We also performed extensive simulation studies to estimate the statistical power of EMMA under various SNP effects, varying degrees of population structure, and differing numbers of multiple measurements per strain. Despite the limited power of inbred mouse association mapping due to the limited number of available inbred strains, we are able to identify significantly associated SNPs, which fall into known QTL or genes identified through previous studies while avoiding an inflation of false positives. An R package implementation and webserver of our EMMA method are publicly available.     

Identification of genetic determinants of IGF‐1 levels and longevity among mouse inbred strains

The IGF‐1 signaling pathway plays an important role in regulating longevity. To identify the genetic loci and genes that regulate plasma IGF‐1 levels, we intercrossed MRL/MpJ and SM/J, inbred mouse strains that differ in IGF‐1 levels. Quantitative trait loci (QTL) analysis of IGF‐1 levels of these F2 mice detected four QTL on chromosomes (Chrs) 9 (48 Mb), 10 (86 Mb), 15 (18 Mb), and 17 (85 Mb). Haplotype association mapping of IGF‐1 levels in 28 domesticated inbred strains identified three suggestive loci in females on Chrs 2 (13 Mb), 10 (88 Mb), and 17 (28 Mb) and in four males on Chrs 1 (159 Mb), 3 (52 and 58 Mb), and 16 (74 Mb). Except for the QTL on Chr 9 and 16, all loci co‐localized with IGF‐1 QTL previously identified in other mouse crosses. The most significant locus was the QTL on Chr 10, which contains the Igf1 gene and which had a LOD score of 31.8. Haplotype analysis among 28 domesticated inbred strains revealed a major QTL on Chr 10 overlapping with the QTL identified in the F2 mice. This locus showed three major haplotypes; strains with haplotype 1 had significantly lower plasma IGF‐1 and extended longevity (P < 0.05) than strains with haplotype 2 or 3. Bioinformatic analysis, combined with sequencing and expression studies, showed that Igf1 is the most likely QTL gene, but that other genes may also play a role in this strong QTL.     

Evidence for association and linkage between atopy, airway hyper-responsiveness, and the β subunit Glu237Gly variant of the high-affinity receptor for immunoglobulin E in the French-Canadian population

Following detection of linkage between atopy and chromosome 11q13 markers, association between this disorder and variants of the beta subunit of the high-affinity receptor for immunoglobulin E (FcepsilonRI-beta, a candidate gene for asthma-related conditions co-localizing within the same region) was reported in Australian, British and Japanese populations. Investigations in several other ethnic groups failed to replicate these observations. Due to the complexity of defining intermediate phenotypes related to asthma, detection of such associations may have been hampered by clinical misclassifications. To assess whether the FcepsilonRI-beta gene was involved in atopy and/or airway hyperresponsiveness (AHR) in the French-Canadian population, we conducted a case-control study in 200 subjects using strict criteria for asthma and related conditions. The Ile181Leu and Glu237Gly FcepsilonRI-beta sequence variants were tested exploiting two amplification refractory mutation systems. No association was detected between atopy or AHR and the Ile181Leu FcepsilonRI-beta variant. However, a strong association was observed between atopy and the Glu237Gly FcepsilonRI-beta variant (odds ratio=12.25). Four large Eastern Québec families (n=106 subjects) were also recruited to perform a genetic linkage study. We observed suggestive evidence of linkage between atopy and the Glu237Gly FcepsilonRI-beta variant (Zmax=2.30). This study is the first to detect the presence of an association between atopy and the Glu237Gly FcepsilonRI-beta variant in French-Canadians. Our data suggest that a susceptibility locus for atopy is located on chromosome 11q13 in this population.     

Fine haplotype structure of a chromosome 17 region in the laboratory and wild mouse

Extensive linkage disequilibrium among classical laboratory strains represents an obstacle in the high-resolution haplotype mapping of mouse quantitative trait loci (QTL). To determine the potential of wild-derived mouse strains for fine QTL mapping, we constructed a haplotype map of a 250-kb region of the t-complex on chromosome 17 containing the Hybrid sterility 1 (Hst1) gene. We resequenced 33 loci from up to 80 chromosomes of five mouse (sub)species. Trans-species single-nucleotide polymorphisms (SNPs) were rare between Mus m. musculus (Mmmu) and Mus m. domesticus (Mmd). The haplotypes in Mmmu and Mmd differed and therefore strains from these subspecies should not be combined for haplotype-associated mapping. The haplotypes of t-chromosomes differed from all non-t Mmmu and Mmd haplotypes. Half of the SNPs and SN indels but only one of seven longer rearrangements found in classical laboratory strains were useful for haplotype mapping in the wild-derived M. m. domesticus. The largest Mmd haplotype block contained three genes of a highly conserved synteny. The lengths of the haplotype blocks deduced from 36 domesticus chromosomes were in tens of kilobases, suggesting that the wild-derived Mmd strains are suitable for fine interval-specific mapping.     

Integration of QTL and bioinformatic tools to identify candidate genes for triglycerides in mice

To identify genetic loci influencing lipid levels, we performed quantitative trait loci (QTL) analysis between inbred mouse strains MRL/MpJ and SM/J, measuring triglyceride levels at 8 weeks of age in F2 mice fed a chow diet. We identified one significant QTL on chromosome (Chr) 15 and three suggestive QTL on Chrs 2, 7, and 17. We also carried out microarray analysis on the livers of parental strains of 282 F2 mice and used these data to find cis-regulated expression QTL. We then narrowed the list of candidate genes under significant QTL using a "toolbox" of bioinformatic resources, including haplotype analysis; parental strain comparison for gene expression differences and nonsynonymous coding single nucleotide polymorphisms (SNP); cis-regulated eQTL in livers of F2 mice; correlation between gene expression and phenotype; and conditioning of expression on the phenotype. We suggest Slc25a7 as a candidate gene for the Chr 7 QTL and, based on expression differences, five genes (Polr3 h, Cyp2d22, Cyp2d26, Tspo, and Ttll12) as candidate genes for Chr 15 QTL. This study shows how bioinformatics can be used effectively to reduce candidate gene lists for QTL related to complex traits.     

The mouse QTL map helps interpret human genome-wide association studies for HDL cholesterol

Genome-wide association (GWA) studies represent a powerful strategy for identifying susceptibility genes for complex diseases in human populations but results must be confirmed and replicated. Because of the close homology between mouse and human genomes, the mouse can be used to add evidence to genes suggested by human studies. We used the mouse quantitative trait loci (QTL) map to interpret results from a GWA study for genes associated with plasma HDL cholesterol levels. We first positioned single nucleotide polymorphisms (SNPs) from a human GWA study on the genomic map for mouse HDL QTL. We then used mouse bioinformatics, sequencing, and expression studies to add evidence for one well-known HDL gene (Abca1) and three newly identified genes (Galnt2, Wwox, and Cdh13), thus supporting the results of the human study. For GWA peaks that occur in human haplotype blocks with multiple genes, we examined the homologous regions in the mouse to prioritize the genes using expression, sequencing, and bioinformatics from the mouse model, showing that some genes were unlikely candidates and adding evidence for candidate genes Mvk and Mmab in one haplotype block and Fads1 and Fads2 in the second haplotype block. Our study highlights the value of mouse genetics for evaluating genes found in human GWA studies.     

Critical role for IL-13 in the development of allergen-induced airway hyperreactivity

Airway hyperresponsiveness to a variety of specific and nonspecific stimuli is a cardinal feature of asthma, which affects nearly 10% of the population in industrialized countries. Eosinophilic pulmonary inflammation, eosinophil-derived products, as well as Th2 cytokines IL-13, IL-4, and IL-5, have been associated with the development of airway hyperreactivity (AHR), but the specific immunological basis underlying the development of AHR remains controversial. Herein we show that mice with targeted deletion of IL-13 failed to develop allergen-induced AHR, despite the presence of vigorous Th2-biased, eosinophilic pulmonary inflammation. However, AHR was restored in IL-13(-/-) mice by the administration of recombinant IL-13. Moreover, adoptive transfer of OVA-specific Th2 cells generated from TCR-transgenic IL-13(-/-) mice failed to induce AHR in recipient SCID mice, although such IL-13(-/-) Th2 cells produced high levels of IL-4 and IL-5 and induced significant airway inflammation. These studies definitively demonstrate that IL-13 is necessary and sufficient for the induction of AHR and that eosinophilic airway inflammation in the absence of IL-13 is inadequate for the induction of AHR. Therefore, treatment of human asthma with antagonists of IL-13 may be very effective.     

Intrinsic and antigen-induced airway hyperresponsiveness are the result of diverse physiological mechanisms.

Airway hyperresponsiveness (AHR) is a defining feature of asthma. We have previously shown, in mice sensitized and challenged with antigen, that AHR is attributable to normal airway smooth muscle contraction with exaggerated airway closure. In the present study we sought to determine if the same was true for mice known to have intrinsic AHR, the genetic strain of mice, A/J. We found that A/J mice have AHR characterized by minimal increase in elastance following aerosolized methacholine challenge compared with mice (BALB/c) that have been antigen sensitized and challenged [concentration that evokes 50% change in elastance (PC(50)): 22.9 +/- 5.7 mg/ml for A/J vs. 3.3 +/- 0.4 mg/ml for antigen-challenged and -sensitized mice; P < 0.004]. Similar results were found when intravenous methacholine was used (PC(30) 0.22 +/- 0.08 mg/ml for A/J vs. 0.03 +/- 0.004 mg/ml for antigen-challenged and -sensitized mice). Computational model analysis revealed that the AHR in A/J mice is dominated by exaggerated airway smooth muscle contraction and that when the route of methacholine administration was changed to intravenous, central airway constriction dominates. Absorption atelectasis was used to provide evidence of the lack of airway closure in A/J mice. Bronchoconstriction during ventilation with 100% oxygen resulted in a mean 9.8% loss of visible lung area in A/J mice compared with 28% in antigen-sensitized and -challenged mice (P < 0.02). We conclude that the physiology of AHR depends on the mouse model used and the route of bronchial agonist administration.