Mouse genome-wide association mapping needs linkage analysis to avoid false-positive Loci

Detection of genomic regions underlying milk production traits in Valle del Belice dairy sheep using regional heritability mapping

The aim of this study was to identify genomic regions underlying milk production traits in the Valle del Belice dairy sheep using regional heritability mapping (RHM). Repeated measurements for milk yield (MY), fat percentage and yield (F% and FY) and protein percentage and yield (P% and PY), collected over a period of 6 years (2006–2012) on 481 Valle del Belice ewes, were used for the analysis. Animals were genotyped with the Illumina 50k SNP chip. Variance components, heritabilities and repeatabilities within and across lactations were estimated, fitting parity, litter size, season of lambing and fortnights in milk, as fixed; and additive genetic, permanent environment within and across lactations, flock by test‐day interaction and residual as random effects. For the RHM analysis, the model included the same fixed and random effects as before, plus an additional regional genomic additive effect (specific for the region being tested) as random. While the whole genomic additive effect was estimated using the genomic relationship matrix (GRM) constructed from all SNPs, the regional genomic additive effect was estimated from a GRM matrix constructed from the SNPs within each region. Heritability estimates ranged between 0.06 and 0.15, with repeatabilities being between 0.14 and 0.24 across lactations and between 0.23 and 0.39 within lactation for all milk production traits. A substantial effect of flock‐test‐day on milk production traits was also estimated. Significant genomic regions at either genome‐wide (p < .05) or suggestive (i.e., one false positive per genome scan) level were identified on chromosome (OAR) 2, 3 and 20 for F% and on OAR3 for P%, with the regions on OAR3 in common between the two traits. Our results confirmed the role of LALBA and AQP genes, on OAR3, as candidate genes for milk production traits in sheep.

Genetic Variation in the Magnitude and Longevity of the IgG Subclass Response to a Diphtheria-Tetanus-Acellular Pertussis (DTaP) Vaccine in Mice

The type of IgG subclasses induced by vaccination is an important determinant of vaccine efficacy because the IgG subclasses vary in their biological function. The goal of this study was to determine the influence of the genetic background on the production and duration of vaccine-induced IgG subclasses. IgG1, IgG2b, and IgG3 titers against diphtheria toxoid (DT), pertussis toxin (PT), filamentous hemagglutinin (FHA), and pertactin (Prn) were measured in mice from 28 different inbred and wild-derived strains vaccinated with an aluminum hydroxide-adjuvanted DTaP vaccine. The titers and duration of vaccine-specific IgG subclass responses were different among mouse strains, indicating that genetic factors contribute to this variation. Statistical associations were used to identify potential mechanisms that contribute to antibody production and longevity. This analysis showed that the mechanisms guiding the magnitude of antibody production were antigen-dependent for IgG1 but antigen-independent for IgG2b and IgG3. However, the mechanisms driving the longevity of antibody titers were antigen-independent for IgG1, IgG2b, and IgG3. The ratio of IgG1 and IgG3 titers identified Th1 and Th2-prone mouse strains. TLR4-deficient C3H/HeJ mice had an enhanced IgG1 response compared with C3H/HeOuJ mice with intact TLR4. This work demonstrates that the genetic background contributes significantly to the magnitude and longevity of vaccine-induced IgG1, IgG2b, and IgG3 titers in mice.

A Solanum neorickii introgression population providing a powerful complement to the extensively characterized Solanum pennellii population

We present a complementary resource for trait fine-mapping in tomato to those based on the intra-specific cross between cultivated tomato and the wild tomato species Solanum pennellii, which have been extensively used for quantitative genetics in tomato over the last 20 years. The current population of backcross inbred lines (BILs) is composed of 107 lines derived after three backcrosses of progeny of the wild species Solanum neorickii (LA2133) and cultivated tomato (cultivar TA209) and is freely available to the scientific community. These S. neorickii BILs were genotyped using the 10K SolCAP single nucleotide polymorphism chip, and 3111 polymorphic markers were used to map recombination break points relative to the physical map of Solanum lycopersicum. The BILs harbor on average 4.3 introgressions per line, with a mean introgression length of 34.7 Mbp, allowing partitioning of the genome into 340 bins and thereby facilitating rapid trait mapping. We demonstrate the power of using this resource in comparison with archival data from the S. pennellii resources by carrying out metabolic quantitative trait locus analysis following gas chromatography-mass spectrometry on fruits harvested from the S. neorickii BILs. The metabolic candidate genes phenylalanine ammonia-lyase and cystathionine gamma-lyase were then tested and validated in F2 populations and via agroinfiltration-based overexpression in order to exemplify the fidelity of this method in identifying the genes that drive tomato metabolic phenotypes.

Variable Susceptibility to Cigarette Smoke-Induced Emphysema in 34 Inbred Strains of Mice Implicates Abi3bp in Emphysema Susceptibility

Chronic obstructive pulmonary disease (COPD) is caused by a complex interaction of environmental exposures, most commonly cigarette smoke, and genetic factors. Chronic cigarette smoke exposure in the mouse is a commonly used animal model of COPD. We aimed to expand our knowledge about the variable susceptibility of inbred strains to this model and test for genetic variants associated with this trait. To that end, we sought to measure differential susceptibility to cigarette smoke-induced emphysema in the mouse, identify genetic loci associated with this quantitative trait, and find homologous human genes associated with COPD. Alveolar chord length (CL) in 34 inbred strains of mice was measured after 6 months of exposure to cigarette smoke. After testing for association, we connected a murine candidate locus to a published meta-analysis of moderate-to-severe COPD. We identified deleterious mutations in a candidate gene in silico and measured gene expression in extreme strains. A/J was the most susceptible strain in our survey (Δ CL 7.0 ± 2.2 μm) and CBA/J was the least susceptible (Δ CL -0.3 ± 1.2 μm). By integrating mouse and human genome-wide scans, we identified the candidate gene Abi3bp. CBA/J mice harbor predicted deleterious variants in Abi3bp, and expression of the gene differs significantly between CBA/J and A/J mice. This is the first report of susceptibility to cigarette smoke-induced emphysema in 34 inbred strains of mice, and Abi3bp is identified as a potential contributor to this phenotype.

Complex genetic control of lung tumorigenesis in resistant mice strains

The SM/J mouse strain is resistant to chemically‐induced lung tumorigenesis despite having a haplotype, in the pulmonary adenoma susceptibility locus (Pas1) locus, that confers tumor susceptibility in other strains. To clarify this inconsistent genotype‐phenotype correlation, we crossed SM/J mice with another resistant strain and conducted genome‐wide linkage analysis in the (C57BL/6J × SM/J)F2 progeny exposed to urethane to induce lung tumors. Overall, >80% of F2 mice of both sexes developed from 1 to 20 lung tumors. Genotyping of 372 F2 mice for 744 informative non‐redundant SNPs dispersed over all autosomal chromosomes revealed four quantitative trait loci (QTLs) affecting lung tumor multiplicity, on chromosomes 3 (near rs13477379), 15 (rs6285067), 17 (rs33373629) and 18 (rs3706601), all with logarithm of the odds (LOD) scores >5. Four QTLs modulated total lung tumor volume, on chromosome 3 (rs13477379), 10 (rs13480702), 15 (rs6285067) and 17 (rs3682923), all with LOD scores >4. No QTL modulating lung tumor multiplicity or total volume was detected in Pas1 on chromosome 6. The present study demonstrates that the SM/J strain carries, at the Pas1 locus, the resistance allele: a finding that will facilitate identification of the Pas1 causal element. More generally, it demonstrates that lung tumorigenesis is under complex polygenic control even in a pedigree with low susceptibility to this neoplasia, suggesting that the genetics of lung tumorigenesis is much more complex than evidenced by the pulmonary adenoma susceptibility and resistance loci that have, so far, been mapped in a small number of crosses between a few inbred strains.

Identification of Padi2 as a novel angiogenesis-regulating gene by genome association studies in mice

Recent findings indicate that growth factor-driven angiogenesis is markedly influenced by genetic variation. This variation in angiogenic responsiveness may alter the susceptibility to a number of angiogenesis-dependent diseases. Here, we utilized the genetic diversity available in common inbred mouse strains to identify the loci and candidate genes responsible for differences in angiogenic response. The corneal micropocket neovascularization assay was performed on 42 different inbred mouse strains using basic fibroblast growth factor (bFGF) pellets. We performed a genome-wide association study utilizing efficient mixed-model association (EMMA) mapping using the induced vessel area from all strains. Our analysis yielded five loci with genome-wide significance on chromosomes 4, 8, 11, 15 and 16. We further refined the mapping on chromosome 4 within a haplotype block containing multiple candidate genes. These genes were evaluated by expression analysis in corneas of various inbred strains and in vitro functional assays in human microvascular endothelial cells (HMVECs). Of these, we found the expression of peptidyl arginine deiminase type II (Padi2), known to be involved in metabolic pathways, to have a strong correlation with a haplotype shared by multiple high angiogenic strains. In addition, inhibition of Padi2 demonstrated a dosage-dependent effect in HMVECs. To investigate its role in vivo, we knocked down Padi2 in transgenic kdrl:zsGreen zebrafish embryos using morpholinos. These embryos had disrupted vessel formation compared to control siblings. The impaired vascular pattern was partially rescued by human PADI2 mRNA, providing evidence for the specificity of the morphant phenotype. Taken together, our study is the first to indicate the potential role of Padi2 as an angiogenesis-regulating gene. The characterization of Padi2 and other genes in associated pathways may provide new understanding of angiogenesis regulation and novel targets for diagnosis and treatment of a wide variety of angiogenesis-dependent diseases.

A small secreted protein in Zymoseptoria tritici is responsible for avirulence on wheat cultivars carrying the Stb6 resistance gene

Zymoseptoria tritici is the causal agent of Septoria tritici blotch, a major pathogen of wheat globally and the most damaging pathogen of wheat in Europe. A gene-for-gene (GFG) interaction between Z. tritici and wheat cultivars carrying the Stb6 resistance gene has been postulated for many years, but the genes have not been identified. We identified AvrStb6 by combining quantitative trait locus mapping in a cross between two Swiss strains with a genome-wide association study using a natural population of c. 100 strains from France. We functionally validated AvrStb6 using ectopic transformations. AvrStb6 encodes a small, cysteine-rich, secreted protein that produces an avirulence phenotype on wheat cultivars carrying the Stb6 resistance gene. We found 16 nonsynonymous single nucleotide polymorphisms among the tested strains, indicating that AvrStb6 is evolving very rapidly. AvrStb6 is located in a highly polymorphic subtelomeric region and is surrounded by transposable elements, which may facilitate its rapid evolution to overcome Stb6 resistance. AvrStb6 is the first avirulence gene to be functionally validated in Z. tritici, contributing to our understanding of avirulence in apoplastic pathogens and the mechanisms underlying GFG interactions between Z. tritici and wheat.

QTL mapping identifies a major locus for resistance in wheat to Sunn pest (Eurygaster integriceps) feeding at the vegetative growth stage

This research provides the first report of a major locus controlling wheat resistance to Sunn pest. It developed and validated SNP markers that will be useful for marker-assisted selection. Sunn pest (Eurygaster integriceps Puton) is the most destructive insect pest of bread wheat and durum wheat in West and Central Asia and East Europe. Breeding for resistance at the vegetative stage of growth is vital in reducing the damage caused by overwintered adult populations that feed on shoot and leaves of seedlings, and in reducing the next generation of pest populations (nymphs and adults), which can cause damage to grain quality by feeding on spikes. In the present study, two doubled haploid (DH) populations involving resistant landraces from Afghanistan were genotyped with the 90k SNP iSelect assay and candidate gene-based KASP markers. The DH lines and parents were phenotyped for resistance to Sunn pest feeding, using artificial infestation cages at Terbol station, in Lebanon, over three years. Quantitative trait locus (QTL) analysis identified a single major locus on chromosome 4BS in the two populations, with the resistance allele derived from the landrace accessions, IG139431 and IG139883. The QTL explained a maximum of 42 % of the phenotypic variation in the Cham6 × IG139431 and 56 % in the Cham6 × IG139883 populations. SNP markers closest to the QTL showed high similarity to rice genes that putatively encode proteins for defense response to herbivory and wounding. The markers were validated in a large, unrelated population of parental wheat genotypes. All wheat lines carrying the 'C-G' haplotype at the identified SNPs were resistant, suggesting that selection based on a haplotype of favourable alleles would be effective in predicting resistance status of unknown genotypes.

Modeling diseases in multiple mouse strains for precision medicine studies

The genetic basis of the phenotypic variability observed in patients can be studied in mice by generating disease models through genetic or chemical interventions in many genetic backgrounds where the clinical phenotypes can be assessed and used for genome-wide association studies (GWAS). This is particularly relevant for rare disorders, where patients sharing identical mutations can present with a wide variety of symptoms, but there are not enough number of patients to ensure statistical power of GWAS. Inbred strains are homozygous for each loci, and their single nucleotide polymorphisms catalogs are known and freely available, facilitating the bioinformatics and reducing the costs of the study, since it is not required to genotype every mouse. This kind of approach can be applied to pharmacogenomics studies as well.

Regional heritability mapping method helps explain missing heritability of blood lipid traits in isolated populations

Single single-nucleotide polymorphism (SNP) genome-wide association studies (SSGWAS) may fail to identify loci with modest effects on a trait. The recently developed regional heritability mapping (RHM) method can potentially identify such loci. In this study, RHM was compared with the SSGWAS for blood lipid traits (high-density lipoprotein (HDL), low-density lipoprotein (LDL), plasma concentrations of total cholesterol (TC) and triglycerides (TG)). Data comprised 2246 adults from isolated populations genotyped using ∼300 000 SNP arrays. The results were compared with large meta-analyses of these traits for validation. Using RHM, two significant regions affecting HDL on chromosomes 15 and 16 and one affecting LDL on chromosome 19 were identified. These regions covered the most significant SNPs associated with HDL and LDL from the meta-analysis. The chromosome 19 region was identified in our data despite the fact that the most significant SNP in the meta-analysis (or any SNP tagging it) was not genotyped in our SNP array. The SSGWAS identified one SNP associated with HDL on chromosome 16 (the top meta-analysis SNP) and one on chromosome 10 (not reported by RHM or in the meta-analysis and hence possibly a false positive association). The results further confirm that RHM can have better power than SSGWAS in detecting causal regions including regions containing crucial ungenotyped variants. This study suggests that RHM can be a useful tool to explain some of the ‘missing heritability' of complex trait variation.

A haplotype block associated with thousand-kernel weight on chromosome 5DS in common wheat (Triticum aestivum L.)

Spike number per unit area, number of grains per spike, and thousand-kernel weight (TKW) are important yield components for wheat (Triticum aestivum L.). TKW has the highest heritability among the three components. We validated 27 simple sequence repeat (SSR) loci associated with TKW in an F(2:5) breeding population grown in four environments. A cfd7(8265bp) marker on chromosome 5DS showed the strongest association with TKW and had a significantly positive effect on TKW compared to allele cfd7(8259bp), with mean increases of 5.17, 3.63, 4.11, and 5.16 g in the four environments. Markers cfd67 and cfd40 flanking cfd78 also showed significantly positive associations with TKW with increases of 5.11, 3.29, 4.31, and 4.50 g for cfd67(205), and 4.98, 3.49, 4.06, and 4.84 g for cfd40(187) compared with cfd67(203) and cfd40(190) in the four environments, respectively. A major quantitative trait locus for TKW spanning 2.94 cM on chromosome 5DS was detected by association mapping. Strong linkage disequilibrium (LD) (r(2)  > 0.2) was detected among the three linked markers, which formed three haplotype blocks in the F(2:5) breeding population. Mean TKW differences between HapB-I and HapB-II were 5.80, 4.41, 4.02, and 5.06 g in the four environments, respectively. Moreover, significant LD was detected only between cfd78 and cfd67 and between cfd67 and cfd40 in a germplasm collection. This study provides a base for cloning genes related to TKW on chromosome 5DS.

Limitation of number of strains and persistence of false positive loci in QTL mapping using recombinant inbred strains

While the identification of causal genes of quantitative trait loci (QTL) remains a difficult problem in the post-genome era, the number of QTL continues to accumulate, mainly identified using the recombinant inbred (RI) strains. Over the last decade, hundreds of publications have reported nearly a thousand QTL identified from RI strains. We hypothesized that the inaccuracy of most of these QTL makes it difficult to identify causal genes. Using data from RI strains derived from C57BL/6J (B6) X DBA/2J (D2), we tested the possibility of detection of reliable QTL with different numbers of strains in the same trait in five different traits. Our results indicated that studies using RI strains of less than 30 in general have a higher probability of failing to detect reliable QTL. Errors in many studies could include false positive loci, switches between QTL with small and major effects, and missing the real major loci. The similar data was obtained from a RI strain population derived from a different pair of parents and a RI strain population of rat. Thus, thousands of reported QTL from studies of RI strains may need to be double-checked for accuracy before proceeding to causal gene identification.

Genome-wide mapping of virulence in brown planthopper identifies loci that break down host plant resistance

Insects and plants have coexisted for over 350 million years and their interactions have affected ecosystems and agricultural practices worldwide. Variation in herbivorous insects' virulence to circumvent host resistance has been extensively documented. However, despite decades of investigation, the genetic foundations of virulence are currently unknown. The brown planthopper (Nilaparvata lugens) is the most destructive rice (Oryza sativa) pest in the world. The identification of the resistance gene Bph1 and its introduction in commercial rice varieties prompted the emergence of a new virulent brown planthopper biotype that was able to break the resistance conferred by Bph1. In this study, we aimed to construct a high density linkage map for the brown planthopper and identify the loci responsible for its virulence in order to determine their genetic architecture. Based on genotyping data for hundreds of molecular markers in three mapping populations, we constructed the most comprehensive linkage map available for this species, covering 96.6% of its genome. Fifteen chromosomes were anchored with 124 gene-specific markers. Using genome-wide scanning and interval mapping, the Qhp7 locus that governs preference for Bph1 plants was mapped to a 0.1 cM region of chromosome 7. In addition, two major QTLs that govern the rate of insect growth on resistant rice plants were identified on chromosomes 5 (Qgr5) and 14 (Qgr14). This is the first study to successfully locate virulence in the genome of this important agricultural insect by marker-based genetic mapping. Our results show that the virulence which overcomes the resistance conferred by Bph1 is controlled by a few major genes and that the components of virulence originate from independent genetic characters. The isolation of these loci will enable the elucidation of the molecular mechanisms underpinning the rice-brown planthopper interaction and facilitate the development of durable approaches for controlling this most destructive agricultural insect.

Allelic variation in a cellulose synthase gene (PtoCesA4) associated with growth and wood properties in Populus tomentosa

Lignocellulosic biomass from trees provides a renewable feedstock for biofuels, lumber, pulp, paper, and other uses. Dissecting the mechanism underlying natural variation of the complex traits controlling growth and lignocellulose biosynthesis in trees can enable marker-assisted breeding to improve wood quality and yield. Here, we combined linkage disequilibrium (LD)-based association analysis with traditional linkage analysis to detect the genetic effect of a Populus tomentosa cellulose synthase gene, PtoCesA4. PtoCesA4 is strongly expressed in developing xylem and leaves. Nucleotide diversity and LD in PtoCesA4, sampled from the P. tomentosa natural distribution, revealed that PtoCesA4 harbors high single nucleotide polymorphism (SNP) diversity (πT = 0.0080 and θw = 0.0098) and low LD (r(2) ≥ 0.1, within 1400 bp), demonstrating that the potential of a candidate-gene-based LD approach in understanding the molecular basis underlying quantitative variation in this species. By combining single SNP, multi-SNP, and haplotype-based associations in an association population of 460 individuals with single SNP linkage analysis in a family-based linkage populations (1200 individuals), we identified three strong associations (false discovery rate Q < 0.05) in both populations. These include two nonsynonymous markers (SNP49 associated with α-cellulose content and SNP59 associated with fiber width) and a noncoding marker (SNP18 associated with α-cellulose content). Variation in RNA transcript abundance among genotypic classes of SNP49 was confirmed in these two populations. Therefore, combining different methods allowed us to examine functional PtoCesA4 allelic variation underlying natural variation in complex quantitative traits related to growth and lignocellulosic biosynthesis.

Genome-wide association studies in mice

Genome-wide association studies (GWASs) have transformed the field of human genetics and have led to the discovery of hundreds of genes that are implicated in human disease. The technological advances that drove this revolution are now poised to transform genetic studies in model organisms, including mice. However, the design of GWASs in mouse strains is fundamentally different from the design of human GWASs, creating new challenges and opportunities. This Review gives an overview of the novel study designs for mouse GWASs, which dramatically improve both the statistical power and resolution compared to classical gene-mapping approaches.

Early gene expression differences in inbred mouse strains with susceptibility to pulmonary adenomas

Lung cancer is the most common cause of cancer-related deaths in both men and women, and effective preventatives are rare due to the difficulty of early detection. Specific gene expression signatures have been identified in individuals that already developed lung cancer. To identify if gene expression differences could be detected in individuals before the onset of the disease, we obtained lung tissues for microarray analysis from young, healthy mice of 9 inbred strains with known differences in their susceptibility to spontaneous pulmonary adenomas when aged. We found that the most common differentially expressed genes among all possible 36 strain comparisons showed significant associations with cancer- and inflammation-related processes. Significant expression differences between susceptible and resistant strains were detected for Aldh3a1, Cxcr1 and 7, Dpt, and Nptx1-genes with known cancer-related functions, and Cd209, Cxcr1 and 7, and Plag2g1b-genes with known inflammatory-related functions. Whereas Aldh3a1, Cd209, Dpt, and Pla2g1b had increased expression, Cxcr1 and 7, and Nptx1 had decreased expression in strains susceptible to pulmonary adenomas. Thus, our study shows that expression differences between susceptible and resistant strains can be detected in young and healthy mice without manifestation of pulmonary adenomas and, thus, may provide an opportunity of early detection. Finally, the identified genes have previously been reported for human non-small cell lung cancer suggesting that molecular pathways may be shared between these two cancer types.

Identification of exercise capacity QTL using association mapping in inbred mice

There are large interindividual differences in exercise capacity. It is well established that there is a genetic basis for these differences. However, the genetic factors underlying this variation are undefined. Therefore, the purpose of this study was to identify novel putative quantitative trait loci (QTL) for exercise capacity by measuring exercise capacity in inbred mice and performing genome-wide association mapping. Exercise capacity, defined as run time and work, was assessed in male mice (n = 6) from 34 strains of classical and wild-derived inbred mice performing a graded treadmill test. Genome-wide association mapping was performed with an efficient mixed-model association (EMMA) algorithm to identify QTL. Exercise capacity was significantly different across strains. Run time varied by 2.7-fold between the highest running strain (C58/J) and the lowest running strain (A/J). These same strains showed a 16.5-fold difference in work. Significant associations were identified for exercise time on chromosomes 1, 2, 7, 11, and 13. The QTL interval on chromosome 2 (~168 Mb) contains one gene, Nfatc2, and overlaps with a suggestive QTL for training responsiveness in humans. These results provide phenotype data on the widest range of inbred strains tested thus far and indicate that genetic background significantly influences exercise capacity. Furthermore, the novel QTLs identified in the current study provide new targets for investigating the underlying mechanisms for variation in exercise capacity.

The genome architecture of the Collaborative Cross mouse genetic reference population

The Collaborative Cross Consortium reports here on the development of a unique genetic resource population. The Collaborative Cross (CC) is a multiparental recombinant inbred panel derived from eight laboratory mouse inbred strains. Breeding of the CC lines was initiated at multiple international sites using mice from The Jackson Laboratory. Currently, this innovative project is breeding independent CC lines at the University of North Carolina (UNC), at Tel Aviv University (TAU), and at Geniad in Western Australia (GND). These institutions aim to make publicly available the completed CC lines and their genotypes and sequence information. We genotyped, and report here, results from 458 extant lines from UNC, TAU, and GND using a custom genotyping array with 7500 SNPs designed to be maximally informative in the CC and used a novel algorithm to infer inherited haplotypes directly from hybridization intensity patterns. We identified lines with breeding errors and cousin lines generated by splitting incipient lines into two or more cousin lines at early generations of inbreeding. We then characterized the genome architecture of 350 genetically independent CC lines. Results showed that founder haplotypes are inherited at the expected frequency, although we also consistently observed highly significant transmission ratio distortion at specific loci across all three populations. On chromosome 2, there is significant overrepresentation of WSB/EiJ alleles, and on chromosome X, there is a large deficit of CC lines with CAST/EiJ alleles. Linkage disequilibrium decays as expected and we saw no evidence of gametic disequilibrium in the CC population as a whole or in random subsets of the population. Gametic equilibrium in the CC population is in marked contrast to the gametic disequilibrium present in a large panel of classical inbred strains. Finally, we discuss access to the CC population and to the associated raw data describing the genetic structure of individual lines. Integration of rich phenotypic and genomic data over time and across a wide variety of fields will be vital to delivering on one of the key attributes of the CC, a common genetic reference platform for identifying causative variants and genetic networks determining traits in mammals.

Investigating natural variation in Drosophila courtship song by the evolve and resequence approach

A primary goal of population genetics is to determine the genetic basis of natural trait variation. We could significantly advance this goal by developing comprehensive genome-wide approaches to link genotype and phenotype in model organisms. Here we combine artificial selection with population-based resequencing to investigate the genetic basis of variation in the interpulse interval (IPI) of Drosophila melanogaster courtship song. We performed divergent selection on replicate populations for only 14 generations, but had considerable power to differentiate alleles that evolved due to selection from those that evolved stochastically. We identified a large number of variants that changed frequency in response to selection for this simple behavior, and they are highly underrepresented on the X chromosome. Though our power was adequate using this experimental technique, the ability to differentiate causal variants from those affected by linked selection requires further development.

Increasing association mapping power and resolution in mouse genetic studies through the use of meta-analysis for structured populations

Genetic studies in mouse models have played an integral role in the discovery of the mechanisms underlying many human diseases. The primary mode of discovery has been the application of linkage analysis to mouse crosses. This approach results in high power to identify regions that affect traits, but in low resolution, making it difficult to identify the precise genomic location harboring the causal variant. Recently, a panel of mice referred to as the hybrid mouse diversity panel (HMDP) has been developed to overcome this problem. However, power in this panel is limited by the availability of inbred strains. Previous studies have suggested combining results across multiple panels as a means to increase power, but the methods employed may not be well suited to structured populations, such as the HMDP. In this article, we introduce a meta-analysis-based method that may be used to combine HMDP studies with F2 cross studies to gain power, while increasing resolution. Due to the drastically different genetic structure of F2s and the HMDP, the best way to combine two studies for a given SNP depends on the strain distribution pattern in each study. We show that combining results, while accounting for these patterns, leads to increased power and resolution. Using our method to map bone mineral density, we find that two previously implicated loci are replicated with increased significance and that the size of the associated is decreased. We also map HDL cholesterol and show a dramatic increase in the significance of a previously identified result.

Flowering time in maize: linkage and epistasis at a major effect locus

In a previous study, we identified a candidate fragment length polymorphism associated with flowering time variation after seven generations of selection for flowering time, starting from the maize inbred line F252. Here, we characterized the candidate region and identified underlying polymorphisms. Then, we combined QTL mapping, association mapping, and developmental characterization to dissect the genetic mechanisms responsible for the phenotypic variation. The candidate region contained the Eukaryotic Initiation Factor (eIF-4A) and revealed a high level of sequence and structural variation beyond the 3'-UTR of eIF-4A, including several insertions of truncated transposable elements. Using a biallelic single-nucleotide polymorphism (SNP) (C/T) in the candidate region, we confirmed its association with flowering time variation in a panel of 317 maize inbred lines. However, while the T allele was correlated with late flowering time within the F252 genetic background, it was correlated with early flowering time in the association panel with pervasive interactions between allelic variation and the genetic background, pointing to underlying epistasis. We also detected pleiotropic effects of the candidate polymorphism on various traits including flowering time, plant height, and leaf number. Finally, we were able to break down the correlation between flowering time and leaf number in the progeny of a heterozygote (C/T) within the F252 background consistent with causal loci in linkage disequilibrium. We therefore propose that both a cluster of tightly linked genes and epistasis contribute to the phenotypic variation for flowering time.

A Cyp2a polymorphism predicts susceptibility to NNK-induced lung tumorigenesis in mice

Lung tumors from smokers as well as lung tumors from mice exposed to tobacco carcinogens such as 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK), often carry mutations in K-ras, which activates downstream-signaling pathways such as PI3K/AKT/mTOR pathway. Mice with genetic deletion of one of three isoforms of AKT were used to investigate the role of AKT in mutant K-ras-induced lung tumorigenesis in mice. Although deletion of Akt1 or Akt2 decreased NNK-induced lung tumor formation by 90%, deletion of Akt2 failed to decrease lung tumorigenesis in two other mouse models driven by mutant K-ras. Genetic mapping showed that Akt2 was tightly linked to the cytochrome P450 Cyp2a locus on chromosome 7. Consequently, targeted deletion of Akt2 created linkage to a strain-specific Cyp2a5 polymorphism that decreased activation of NNK in vitro. Mice with this Cyp2a5 polymorphism had decreased NNK-induced DNA adduct formation in vivo and decreased NNK-induced lung tumorigenesis. These studies support human epidemiological studies linking CYP2A polymorphisms with lung cancer risk in humans and highlight the need to confirm phenotypes of genetically engineered mice in multiple mouse strains.

Identification of fat4 and tsc22d1 as novel candidate genes for spontaneous pulmonary adenomas

Genetic influences that underlie spontaneous lung oncogenesis are poorly understood. The objective of this study was to determine the genetic influences on spontaneous pulmonary adenoma frequency and severity in 28 strains of mice as part of a large-scale aging study conducted at the Jackson Aging Center (http://agingmice.jax.org/). Genome-wide association studies were conducted in these strains with both low-density (132,000) and high-density (4,000,000) panel of single-nucleotide polymorphisms (SNP). Our analysis revealed that adenomas were relatively less frequent and less severe in females than males, and that loci implicated in frequency and severity were often different between male and female mice. While some of the significant loci identified mapped to genomic locations known to be responsible for carcinogen-induced cancers (e.g., Pas1), others were unique to our study. In particular, Fat4 was influential in males and Tsc22d1 was influential in females. SNPs implicated were predicted to alter amino acid sequence and change protein function. In summary, our results suggested that genetic influences that underlie pulmonary adenoma frequency are dependent on gender, and that Fat4 and Tsc22d1 are likely candidate genes to influence formation of spontaneous pulmonary adenoma in aging male and female mice, respectively.

A dual role for the immune response in a mouse model of inflammation-associated lung cancer

Lung cancer is the leading cause of cancer death worldwide. Both principal factors known to cause lung cancer, cigarette smoke and asbestos, induce pulmonary inflammation, and pulmonary inflammation has recently been implicated in several murine models of lung cancer. To further investigate the role of inflammation in the development of lung cancer, we generated mice with combined loss of IFN-γ and the β-common cytokines GM-CSF and IL-3. These immunodeficient mice develop chronic pulmonary inflammation and lung tumors at a high frequency. Examination of the relationship between these tumors and their inflammatory microenvironment revealed a dual role for the immune system in tumor development. The inflammatory cytokine IL-6 promoted optimal tumor growth, yet wild-type mice rejected transplanted tumors through the induction of adaptive immunity. These findings suggest a model whereby cytokine deficiency leads to oncogenic inflammation that combines with defective antitumor immunity to promote lung tumor formation, representing a unique system for studying the role of the immune system in lung tumor development.

Genetic control of susceptibility to infection with Candida albicans in mice

Candida albicans is an opportunistic pathogen that causes acute disseminated infections in immunocompromised hosts, representing an important cause of morbidity and mortality in these patients. To study the genetic control of susceptibility to disseminated C. albicans in mice, we phenotyped a group of 23 phylogenetically distant inbred strains for susceptibility to infection as measured by extent of fungal replication in the kidney 48 hours following infection. Susceptibility was strongly associated with the loss-of-function mutant complement component 5 (C5/Hc) allele, which is known to be inherited by approximately 40% of inbred strains. Our survey identified 2 discordant strains, AKR/J (C5-deficient, resistant) and SM/J (C5-sufficient, susceptible), suggesting that additional genetic effects may control response to systemic candidiasis in these strains. Haplotype association mapping in the 23 strains using high density SNP maps revealed several putative loci regulating the extent of C. albicans replication, amongst which the most significant were C5 (P value = 2.43×10(-11)) and a novel effect on distal chromosome 11 (P value = 7.63×10(-9)). Compared to other C5-deficient strains, infected AKR/J strain displays a reduced fungal burden in the brain, heart and kidney, and increased survival, concomitant with uniquely high levels of serum IFNγ. C5-independent genetic effects were further investigated by linkage analysis in an [A/JxAKR/J]F2 cross (n = 158) where the mutant Hc allele is fixed. These studies identified a chromosome 11 locus (Carg4, Candida albicans resistance gene 4; LOD = 4.59), and a chromosome 8 locus (Carg3; LOD = 3.95), both initially detected by haplotype association mapping. Alleles at both loci were inherited in a co-dominant manner. Our results verify the important effect of C5-deficiency in inbred mouse strains, and further identify two novel loci, Carg3 and Carg4, which regulate resistance to C. albicans infection in a C5-independent manner.

Towards identifying genes underlying ecologically relevant traits in Arabidopsis thaliana

A major challenge in evolutionary biology and plant breeding is to identify the genetic basis of complex quantitative traits, including those that contribute to adaptive variation. Here we review the development of new methods and resources to fine-map intraspecific genetic variation that underlies natural phenotypic variation in plants. In particular, the analysis of 107 quantitative traits reported in the first genome-wide association mapping study in Arabidopsis thaliana sets the stage for an exciting time in our understanding of plant adaptation. We also argue for the need to place phenotype-genotype association studies in an ecological context if one is to predict the evolutionary trajectories of plant species.

In silico QTL mapping of basal liver iron levels in inbred mouse strains

Both iron deficiency and iron excess are detrimental in many organisms, and previous studies in both mice and humans suggest that genetic variation may influence iron status in mammals. However, these genetic factors are not well defined. To address this issue, we measured basal liver iron levels in 18 inbred strains of mice of both sexes on a defined iron diet and found ∼4-fold variation in liver iron in males (lowest 153 μg/g, highest 661 μg/g) and ∼3-fold variation in females (lowest 222 μg/g, highest 658 μg/g). We carried out a genome-wide association mapping to identify haplotypes underlying differences in liver iron and three other related traits (copper and zinc liver levels, and plasma diferric transferrin levels) in a subset of 14 inbred strains for which genotype information was available. We identified two putative quantitative trait loci (QTL) that contain genes with a known role in iron metabolism: Eif2ak1 and Igf2r. We also identified four putative QTL that reside in previously identified iron-related QTL and 22 novel putative QTL. The most promising putative QTL include a 0.22 Mb region on Chromosome 7 and a 0.32 Mb region on Chromosome 11 that both contain only one candidate gene, Adam12 and Gria1, respectively. Identified putative QTL are good candidates for further refinement and subsequent functional studies.

Multiple linked quantitative trait loci within the Tmevd2/Eae3 interval control the severity of experimental allergic encephalomyelitis in DBA/2J mice

Theiler’s murine encephalomyelitis virus-induced demyelination (TMEVD) and experimental allergic encephalomyelitis (EAE) are the principal animal models of multiple sclerosis (MS). Previously we provided evidence that Tmevd2 and Eae3 may represent either a shared susceptibility locus or members of a gene complex controlling susceptibility to CNS inflammatory demyelinating disease. To explore the genetic relationship between Tmevd2 and Eae3, we generated a D2.C-Tmevd2 interval-specific congenic (ISC) line and three overlapping interval-specific recombinant congenic (ISRC) lines in which the Tmevd2 resistant allele from BALB/cByJ was introgressed onto the TMEVD-susceptible DBA/2J background. These mice, all H2^d, were studied for susceptibility to EAE elicited by immunization with PLP_180–199. Compared to DBA/2J mice, D2.C-Tmevd2 mice developed a significantly less severe clinical disease course and EAE pathology in the spinal cord, confirming the existence of Eae3 and its linkage to Tmevd2 in this strain combination. Compare to DBA/2J and D2.C-Tmevd2, all three ISRC lines exhibited clinical disease courses of intermediate severity. Neither differences in ex vivo antigen-specific cytokine nor proliferative responses uniquely cosegregated with differences in disease severity. These results indicate that multiple QTL within the Tmevd2/Eae3 interval influence EAE severity, one of which includes a homology region for a QTL found in MS by admixture mapping.

Sequence variation at multiple loci influences red cell hemoglobin concentration

A substantial genetic contribution underlies variation in baseline peripheral blood counts. We performed quantitative trait locus/loci analyses to identify chromosome regions harboring genes influencing red cell hemoglobin concentration using the cell hemoglobin concentration mean (CHCM), a directly measured parameter analogous to the mean cell hemoglobin concentration. Fourteen significant loci (gene symbols Chcmq1-Chcmq14) were detected. Seven of these influenced CHCM in a sex-specific fashion, and 2 showed significant interactive effects (epistasis). For quantitative trait locus/loci detected in multiple crosses, confidence intervals were narrowed using statistical and bioinformatic approaches. Two strong candidate genes emerged and were further analyzed: adult β-globin (Hbb) for Chcmq3 on Chr 7, and transferrin (Trf) for Chcmq2 on Chr 9. High and low allele parental strains in crosses detecting Chcmq3 segregate 100% with the known ancestral haplotype blocks, hemoglobin (Hb) diffuse (Hbb(d)) and Hb single (Hbb(s)), respectively. Hbb(d) consists of nonidentical major and minor polypeptides and exhibits an increased positive charge relative to Hbb(s) due to the net loss of 2 negative residues in the Hbb(dminor) polypeptide, resulting in a pI of 7.85 versus 7.13. Thus, as shown in human erythrocytes, positively charged Hbs are associated with cell dehydration and increased CHCM in mouse erythrocytes.

Comparison of unrestrained plethysmography and forced oscillation for identifying genetic variability of airway responsiveness in inbred mice

Lung function detection in mice is currently most accurately measured by invasive techniques, which are costly, labor intensive, and terminal. This limits their use for large-scale or longitudinal studies. Noninvasive assays are often used instead, but their accuracy for measuring lung function parameters such as resistance and elastance has been questioned in studies involving small numbers of mouse strains. Here we compared parameters detected by two different methods using 29 inbred mouse strains: enhanced pause (Penh), detected by unrestrained plethysmography, and central airway resistance and lung elastance, detected by a forced oscillation technique. We further tested whether the phenotypic variations were determined by the same genomic location in genome-wide association studies using a linear mixed model algorithm. Penh, resistance, and elastance were measured in nonexposed mice or mice exposed to saline and increasing doses of aerosolized methacholine. Because Penh differed from airway resistance in several strains and because the peak genetic associations found for Penh, resistance, or elastance were located at different genomic regions, we conclude that using Penh as an indicator for lung function changes in high-throughput genetic studies (i.e., genome-wide association studies or quantitative trait locus studies) measures something fundamentally different than airway resistance and lung elastance.

Mapping quantitative traits and strategies to find quantitative trait genes

In 1999 a meeting took place at the Jackson Laboratory, a large mouse research centre in Bar Harbor, Maine, to consider the value of systematically collecting phenotypes on inbred strains of mice (Paigen and Eppig (2000) [1]). The group concluded that cataloguing the extensive phenotypic diversity present among laboratory mice, and in particular providing the research community with data from cohorts of animals, phenotyped according to standardized protocols, was essential if we were to take advantage of the possibilities of mouse genetics. Beginning with the collection of basic physiological, biochemical and behavioral data on nine commonly used inbred strains, the project has expanded so that by the beginning of 2010 data for 178 strains had been collected, with 105 phenotype projects yielding over 2000 different measurements (Bogue et al. (2007) [2].

Overview of techniques to account for confounding due to population stratification and cryptic relatedness in genomic data association analyses

Population-based genomic association analyses are more powerful than within-family analyses. However, population stratification (unknown or ignored origin of individuals from multiple source populations) and cryptic relatedness (unknown or ignored covariance between individuals because of their relatedness) are confounding factors in population-based genomic association analyses, which inflate the false-positive rate. As a consequence, false association signals may arise in genomic data association analyses for reasons other than true association between the tested genomic factor (marker genotype, gene or protein expression) and the study phenotype. It is therefore important to correct or account for these confounders in population-based genomic data association analyses. The common correction techniques for population stratification and cryptic relatedness problems are presented here in the phenotype-marker association analysis context, and comments on their suitability for other types of genomic association analyses (for example, phenotype-expression association) are also provided. Even though many of these techniques have originally been developed in the context of human genetics, most of them are also applicable to model organisms and breeding populations.

Linkage and association mapping of Arabidopsis thaliana flowering time in nature

Flowering time is a key life-history trait in the plant life cycle. Most studies to unravel the genetics of flowering time in Arabidopsis thaliana have been performed under greenhouse conditions. Here, we describe a study about the genetics of flowering time that differs from previous studies in two important ways: first, we measure flowering time in a more complex and ecologically realistic environment; and, second, we combine the advantages of genome-wide association (GWA) and traditional linkage (QTL) mapping. Our experiments involved phenotyping nearly 20,000 plants over 2 winters under field conditions, including 184 worldwide natural accessions genotyped for 216,509 SNPs and 4,366 RILs derived from 13 independent crosses chosen to maximize genetic and phenotypic diversity. Based on a photothermal time model, the flowering time variation scored in our field experiment was poorly correlated with the flowering time variation previously obtained under greenhouse conditions, reinforcing previous demonstrations of the importance of genotype by environment interactions in A. thaliana and the need to study adaptive variation under natural conditions. The use of 4,366 RILs provides great power for dissecting the genetic architecture of flowering time in A. thaliana under our specific field conditions. We describe more than 60 additive QTLs, all with relatively small to medium effects and organized in 5 major clusters. We show that QTL mapping increases our power to distinguish true from false associations in GWA mapping. QTL mapping also permits the identification of false negatives, that is, causative SNPs that are lost when applying GWA methods that control for population structure. Major genes underpinning flowering time in the greenhouse were not associated with flowering time in this study. Instead, we found a prevalence of genes involved in the regulation of the plant circadian clock. Furthermore, we identified new genomic regions lacking obvious candidate genes.

Fine mapping in 94 inbred mouse strains using a high-density haplotype resource

The genetics of phenotypic variation in inbred mice has for nearly a century provided a primary weapon in the medical research arsenal. A catalog of the genetic variation among inbred mouse strains, however, is required to enable powerful positional cloning and association techniques. A recent whole-genome resequencing study of 15 inbred mouse strains captured a significant fraction of the genetic variation among a limited number of strains, yet the common use of hundreds of inbred strains in medical research motivates the need for a high-density variation map of a larger set of strains. Here we report a dense set of genotypes from 94 inbred mouse strains containing 10.77 million genotypes over 121,433 single nucleotide polymorphisms (SNPs), dispersed at 20-kb intervals on average across the genome, with an average concordance of 99.94% with previous SNP sets. Through pairwise comparisons of the strains, we identified an average of 4.70 distinct segments over 73 classical inbred strains in each region of the genome, suggesting limited genetic diversity between the strains. Combining these data with genotypes of 7570 gap-filling SNPs, we further imputed the untyped or missing genotypes of 94 strains over 8.27 million Perlegen SNPs. The imputation accuracy among classical inbred strains is estimated at 99.7% for the genotypes imputed with high confidence. We demonstrated the utility of these data in high-resolution linkage mapping through power simulations and statistical power analysis and provide guidelines for developing such studies. We also provide a resource of in silico association mapping between the complex traits deposited in the Mouse Phenome Database with our genotypes. We expect that these resources will facilitate effective designs of both human and mouse studies for dissecting the genetic basis of complex traits.

Assessing the prospects of genome-wide association studies performed in inbred mice

The remarkable success in mapping genes linked to a number of disease traits using genome-wide association studies (GWAS) in human cohorts has renewed interest in applying this same technique in model organisms such as inbred laboratory mice. Unlike humans, however, the limited genetic diversity in the ancestry of laboratory mice combined with selection pressure over the past decades have yielded an intricate population genetic structure that can complicate the results obtained from association studies. This problem is further exacerbated by the small number of strains typically used in such studies where multiple spurious associations arise as a result of random chance. We sought to empirically assess the viability of GWAS in inbred mice using hundreds of expression traits for which the true location of the expression quantitative trait locus was known a priori. We then measured transcript abundance levels for these expression traits in 16 classical and 3 wild-derived inbred strains and carried out a genome-wide association scan, demonstrating the low statistical power of such studies and empirically estimating the large extent to which allelic association of transcripts gives rise to spurious associations. We provide evidence illustrating that in a large fraction of cases, the marker with the most significant p values fails to map to the location of the true eQTL. Finally, we provide experimental support for hundreds of traits, and that combining linkage analysis with association mapping provides significant increases in statistical power over a stand-alone GWAS as well as significantly higher mapping resolution than either study alone.

A high-resolution association mapping panel for the dissection of complex traits in mice

Systems genetics relies on common genetic variants to elucidate biologic networks contributing to complex disease-related phenotypes. Mice are ideal model organisms for such approaches, but linkage analysis has been only modestly successful due to low mapping resolution. Association analysis in mice has the potential of much better resolution, but it is confounded by population structure and inadequate power to map traits that explain less than 10% of the variance, typical of mouse quantitative trait loci (QTL). We report a novel strategy for association mapping that combines classic inbred strains for mapping resolution and recombinant inbred strains for mapping power. Using a mixed model algorithm to correct for population structure, we validate the approach by mapping over 2500 cis-expression QTL with a resolution an order of magnitude narrower than traditional QTL analysis. We also report the fine mapping of metabolic traits such as plasma lipids. This resource, termed the Hybrid Mouse Diversity Panel, makes possible the integration of multiple data sets and should prove useful for systems-based approaches to complex traits and studies of gene-by-environment interactions.

An experimental assessment of in silico haplotype association mapping in laboratory mice

Background

To assess the utility of haplotype association mapping (HAM) as a quantitative trait locus (QTL) discovery tool, we conducted HAM analyses for red blood cell count (RBC) and high density lipoprotein cholesterol (HDL) in mice. We then experimentally tested each HAM QTL using published crosses or new F2 intercrosses guided by the haplotype at the HAM peaks.

Results

The HAM for RBC, using 33 classic inbred lines, revealed 8 QTLs; 2 of these were true positives as shown by published crosses. A HAM-guided (C57BL/6J × CBA/J)F2 intercross we carried out verified 2 more as true positives and 4 as false positives. The HAM for HDL, using 81 strains including recombinant inbred lines and chromosome substitution strains, detected 46 QTLs. Of these, 36 were true positives as shown by published crosses. A HAM-guided (C57BL/6J × A/J)F2 intercross that we carried out verified 2 more as true positives and 8 as false positives. By testing each HAM QTL for RBC and HDL, we demonstrated that 78% of the 54 HAM peaks were true positives and 22% were false positives. Interestingly, all false positives were in significant allelic association with one or more real QTL.

Conclusion

Because type I errors (false positives) can be detected experimentally, we conclude that HAM is useful for QTL detection and narrowing. We advocate the powerful and economical combined approach demonstrated here: the use of HAM for QTL discovery, followed by mitigation of the false positive problem by testing the HAM-predicted QTLs with small HAM-guided experimental crosses.

Genetic loci that regulate healing and regeneration in LG/J and SM/J mice

MRL mice display unusual healing properties. When MRL ear pinnae are hole punched, the holes close completely without scarring, with regrowth of cartilage and reappearance of both hair follicles and sebaceous glands. Studies using (MRL/lpr x C57BL/6)F(2) and backcross mice first showed that this phenomenon was genetically determined and that multiple loci contributed to this quantitative trait. The lpr mutation itself, however, was not one of them. In the present study we examined the genetic basis of healing in the Large (LG/J) mouse strain, a parent of the MRL mouse and a strain that shows the same healing phenotype. LG/J mice were crossed with Small (SM/J) mice and the F(2) population was scored for healing and their genotypes determined at more than 200 polymorphic markers. As we previously observed for MRL and (MRL x B6)F(2) mice, the wound-healing phenotype was sexually dimorphic, with female mice healing more quickly and more completely than male mice. We found quantitative trait loci (QTLs) on chromosomes (Chrs) 9, 10, 11, and 15. The heal QTLs on Chrs 11 and 15 were linked to differential healing primarily in male animals, whereas QTLs on Chrs 9 and 10 were not sexually dimorphic. A comparison of loci identified in previous crosses with those in the present report using LG/J x SM/J showed that loci on Chrs 9, 11, and 15 colocalized with those seen in previous MRL crosses, whereas the locus on Chr 10 was not seen before and is contributed by SM/J.

Association mapping: critical considerations shift from genotyping to experimental design

The goal of many plant scientists' research is to explain natural phenotypic variation in terms of simple changes in DNA sequence. Traditionally, linkage mapping has been the most commonly employed method to reach this goal: experimental crosses are made to generate a family with known relatedness, and attempts are made to identify cosegregation of genetic markers and phenotypes within this family. In vertebrate systems, association mapping (also known as linkage disequilibrium mapping) is increasingly being adopted as the mapping method of choice. Association mapping involves searching for genotype-phenotype correlations in unrelated individuals and often is more rapid and cost-effective than traditional linkage mapping. We emphasize here that linkage and association mapping are complementary approaches and are more similar than is often assumed. Unlike in vertebrates, where controlled crosses can be expensive or impossible (e.g., in humans), the plant scientific community can exploit the advantages of both controlled crosses and association mapping to increase statistical power and mapping resolution. While the time and money required for the collection of genotype data were critical considerations in the past, the increasing availability of inexpensive DNA sequencing and genotyping methods should prompt researchers to shift their attention to experimental design. This review provides thoughts on finding the optimal experimental mix of association mapping using unrelated individuals and controlled crosses to identify the genes underlying phenotypic variation.

Replication and narrowing of gene expression quantitative trait loci using inbred mice

Gene expression quantitative trait locus (eQTL) mapping has become a powerful tool in systems biology. While many authors have made important discoveries using this approach, one persistent challenge in eQTL studies is the selection of loci and genes that should receive further biological investigation. In this study we compared eQTL generated from gene expression profiling in the livers of two panels of mouse strains: 41 BXD recombinant inbred and 36 Mouse Diversity Panel (MDP) strains. Cis-eQTL, loci in which the transcript and its maximum QTL are colocated, have been shown to be more reproducible than trans-eQTL, which are not colocated with the transcript. We observed that between 9.9 and 12.1% of cis-eQTL and between 2.0 and 12.6% of trans-eQTL replicated between the two panels depending on the degree of statistical stringency. Notably, a significant eQTL hotspot on distal chromosome 12 observed in the BXD panel was reproduced in the MDP. Furthermore, the shorter linkage disequilibrium in the MDP strains allowed us to considerably narrow the locus and limit the number of candidate genes to a cluster of Serpin genes, which code for extracellular proteases. We conclude that this strategy has some utility in increasing confidence and resolution in eQTL mapping studies; however, due to the high false-positive rate in the MDP, eQTL mapping in inbred strains is best carried out in combination with an eQTL linkage study.