Role of genetic modifiers in an orthologous rat model of ARPKD

Human data and animal models of autosomal recessive polycystic kidney disease (ARPKD) suggest that genetic factors modulate the onset and severity of the disease. We report here for the first time that ARPKD susceptibility is attenuated by introgressing the mutated Pkhd1 disease allele from the polycystic kidney (PCK) rat onto the FHH (Fawn-Hooded Hypertensive) genetic background. Compared with PCK, the FHH.Pkhd1 strain had significantly decreased renal cyst formation that coincided with a threefold reduction in mean kidney weights. Further analysis revealed that the FHH. Pkhd1 is protected from increased blood pressure as well as elevated plasma creatinine and blood urea nitrogen levels. On the other hand, liver weight and biliary cystogenesis revealed no differences between PCK and FHH.Pkdh1, indicating that genes within the FHH genetic background prevent the development of renal, but not hepatic, manifestations of ARPKD. Microarray expression analysis of kidneys from 30-day-old PCK rats revealed increased expression of genes previously identified in PKD renal expression profiles, such as inflammatory response, extracellular matrix synthesis, and cell proliferation genes among others, whereas the FHH.Pkhd1 did not show activation of these common markers of disease. This newly developed strain can serve as a tool to map modifier genes for renal disease in ARPKD and provides further insight into disease variability and pathophysiology.     

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.     

Structural model analysis of multiple quantitative traits

We introduce a method for the analysis of multilocus, multitrait genetic data that provides an intuitive and precise characterization of genetic architecture. We show that it is possible to infer the magnitude and direction of causal relationships among multiple correlated phenotypes and illustrate the technique using body composition and bone density data from mouse intercross populations. Using these techniques we are able to distinguish genetic loci that affect adiposity from those that affect overall body size and thus reveal a shortcoming of standardized measures such as body mass index that are widely used in obesity research. The identification of causal networks sheds light on the nature of genetic heterogeneity and pleiotropy in complex genetic systems.     

Quantitative trait loci for urinary albumin in crosses between C57BL/6J and A/J inbred mice in the presence and absence of Apoe

We investigated the effect of apolipoprotein E (Apoe) on albuminuria in the males of two independent F2 intercrosses between C57BL/6J and A/J mice, using wild-type inbred strains in the first cross and B6-Apoe(-/-) animals in the second cross. In the first cross, we identified three quantitative trait loci (QTL): chromosome (Chr) 2 [LOD 3.5, peak at 70 cM, confidence interval (C.I.) 28-88 cM]; Chr 9 (LOD 2.0, peak 5 cM, C.I. 5-25 cM); and Chr 19 (LOD 1.9, peak 49 cM, C.I. 23-54 cM). The Chr 2 and Chr 19 QTL were concordant with previously found QTL for renal damage in rat and human. The Chr 9 QTL was concordant with a locus found in rat. The second cross, testing only Apoe(-/-) progeny, did not identify any of these loci, but detected two other loci on Chr 4 (LOD 3.2, peak 54 cM, C.I. 29-73 cM) and Chr 6 (LOD 2.6, peak 33 cM, C.I. 11-61 cM), one of which was concordant with a QTL found in rat. The dependence of QTL detection on the presence of Apoe and the concordance of these QTL with rat and human kidney disease QTL suggest that Apoe plays a role in renal damage.     

Aging in inbred strains of mice: study design and interim report on median lifespans and circulating IGF1 levels

To better characterize aging in mice, the Jackson Aging Center carried out a lifespan study of 31 genetically-diverse inbred mouse strains housed in a specific pathogen-free facility. Clinical assessments were carried out every 6 months, measuring multiple age-related phenotypes including neuromuscular, kidney and heart function, body composition, bone density, hematology, hormonal levels, and immune system parameters. In a concurrent cross-sectional study of the same 31 strains at 6, 12, and 20 months, more invasive measurements were carried out followed by necropsy to assess apoptosis, DNA repair, chromosome fragility, and histopathology. In this report, which is the initial paper of a series, the study design, median lifespans, and circulating insulin-like growth factor 1 (IGF1) levels at 6, 12, and 18 months are described for the first cohort of 32 females and 32 males of each strain. Survival curves varied dramatically among strains with the median lifespans ranging from 251 to 964 days. Plasma IGF1 levels, which also varied considerably at each time point, showed an inverse correlation with a median lifespan at 6 months ( R  = −0.33, P  = 0.01). This correlation became stronger if the short-lived strains with a median lifespan < 600 days were removed from the analysis ( R  = −0.53, P  < 0.01). These results support the hypothesis that the IGF1 pathway plays a key role in regulating longevity in mice and indicates that common genetic mechanisms may exist for regulating IGF1 levels and lifespan.     

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.     

Genetic basis of HDL variation in 129/SvImJ and C57BL/6J mice: importance of testing candidate genes in targeted mutant mice

To evaluate the effect of genetic background on high-density lipoprotein cholesterol (HDL) levels in Soat1(-/-) mice, we backcrossed sterol O-acyltransferase 1 (Soat1)(-/-) mice, originally reported to have elevated HDL levels, to C57BL/6 mice and constructed a congenic strain with only a small region (3.3Mb) of 129 alleles, specifically excluding the nearby apolipoprotein A-II (Apoa2) gene from 129. HDL levels in these Soat1(-/-) mice were no different from C57BL/6, indicating that the passenger gene Apoa2 caused the previously reported elevation of HDL in these Soat1(-/-) mice. Because many knockouts are made in strain 129 and then subsequently backcrossed into C57BL/6, it is important to identify quantitative trait loci (QTL) that differ between 129 and C57BL/6 so that one can guard against effects ascribed to a knockout but really caused by a passenger gene from 129. To provide such data, we generated 528 F(2) progeny from an intercross of 129S1/SvImJ and C57BL/6 and measured HDL concentrations in F(2) animals first fed chow and then atherogenic diet. A genome wide scan using 508 single-nucleotide polymorphisms (SNPs) identified 19 QTL, 2 of which were male specific and 2 were female specific. Using comparative genomics and haplotype analysis, we narrowed QTL on chromosomes 3, 5, 8, 17, and 18 to 0.5, 6.3, 2.6, 1.1, and 0.6 Mb, respectively. These data will serve as a reference for any effort to test the impact of candidate genes on HDL using a knockout strategy.     

A survey of airway responsiveness in 36 inbred mouse strains facilitates gene mapping studies and identification of quantitative trait loci

Airway hyper-responsiveness (AHR) is a critical phenotype of human asthma and animal models of asthma. Other studies have measured AHR in nine mouse strains, but only six strains have been used to identify genetic loci underlying AHR. Our goals were to increase the genetic diversity of available strains by surveying 27 additional strains, to apply haplotype association mapping to the 36-strain survey, and to identify new genetic determinants for AHR. We derived AHR from the increase in airway resistance in females subjected to increasing levels of methacholine concentrations. We used haplotype association mapping to identify associations between AHR and haplotypes on chromosomes 3, 5, 8, 12, 13, and 14. And we used bioinformatics techniques to narrow the identified region on chromosome 13, reducing the region to 29 candidate genes, with 11 of considerable interest. Our combined use of haplotype association mapping with bioinformatics tools is the first study of its kind for AHR on these 36 strains of mice. Our analyses have narrowed the possible QTL genes and will facilitate the discovery of novel genes that regulate AHR in mice.     

Quantitative trait mapping in a diallel cross of recombinant inbred lines

A recombinant inbred intercross (RIX) is created by generating diallel F₁ progeny from one or more panels of recombinant inbred (RI) strains. This design was originally introduced to extend the power of small RI panels for the confirmation of quantitative trait loci (QTL) provisionally detected in a parental RI set. For example, the set of 13 C × B (C57BL/6ByJ × BALB/cByJ) RI strains can, in principle, be supplemented with 156 isogenic F₁s. We describe and test a method of analysis, based on a linear mixed model, that accounts for the correlation structure of RIX populations. This model suggests a novel permutation algorithm that is needed to obtain appropriate threshold values for genome-wide scans of an RIX population. Despite the combinational multiplication of unique genotypes that can be generated using an RIX design, the effective sample size of the RIX population is limited by the number of progenitor RI genomes that are combined. When using small RI panels such as the C × B there appears to be only modest advantage of the RIX design when compared with the original RI panel for detecting QTLs with additive effects. The RIX, however, does have an inherent ability to detect dominance effects, and, unlike RI strains, the RIX progeny are genetically reproducible but are not fully inbred, providing somewhat more natural genetic context. We suggest a breeding strategy, the balanced partial RIX, that balances the advantage of RI and RIX designs. This involves the use of a partial RIX population derived from a large RI panel in which the available information is maximized by minimizing correlations among RIX progeny.