Search

Power analyses are often used to determine the number of animals required for a genome-wide association study (GWAS). These analyses are typically intended to estimate the sample size needed for at least 1 locus to exceed a genome-wide significance threshold. A related question that is less commonly considered is the number of significant loci that will be discovered with a given sample size. We used simulations based on a real data set that consisted of 3,173 male and female adult N/NIH heterogeneous stock rats to explore the relationship between sample size and the number of significant loci discovered. Our simulations examined the number of loci identified in subsamples of the full data set. The subsampling analysis was conducted for 4 traits with low (0.15 ± 0.03), medium (0.31 ± 0.03 and 0.36 ± 0.03), and high (0.46 ± 0.03) SNP-based heritabilities. For each trait, we subsampled the data 100 times at different sample sizes (500, 1,000, 1,500, 2,000, and 2,500). We observed an exponential increase in the number of significant loci with larger sample sizes. Our results are consistent with similar observations in human GWAS and imply that future rodent GWAS should use sample sizes that are significantly larger than those needed to obtain a single significant result.

The development of model systems requires a detailed assessment of standing genetic variation across natural populations. The Brachypodium species complex has been promoted as a plant model for grass genomics with translation to small grain and biomass crops. To capture the genetic diversity within this species complex, thousands of Brachypodium accessions from around the globe were collected and genotyped by sequencing. Overall, 1897 samples were classified into two diploid or allopolyploid species, and then further grouped into distinct inbred genotypes. A core set of diverse B. distachyon diploid lines was selected for whole genome sequencing and high resolution phenotyping. Genome-wide association studies across simulated seasonal environments was used to identify candidate genes and pathways tied to key life history and agronomic traits under current and future climatic conditions. A total of 8, 22, and 47 QTL were identified for flowering time, early vigor, and energy traits, respectively. The results highlight the genomic structure of the Brachypodium species complex, and the diploid lines provided a resource that allows complex trait dissection within this grass model species.

A major consideration in multitrait analysis is which traits should be jointly analyzed. As a common strategy, multitrait analysis is performed either on pairs of traits or on all of traits. To fully exploit the power of multitrait analysis, we propose variable selection to choose a subset of informative traits for multitrait quantitative trait locus (QTL) mapping. The proposed method is very useful for achieving optimal statistical power for QTL identification and for disclosing the most relevant traits. It is also a practical strategy to effectively take advantage of multitrait analysis when the number of traits under consideration is too large, making the usual multivariate analysis of all traits challenging. We study the impact of selection bias and the usage of permutation tests in the context of variable selection and develop a powerful implementation procedure of variable selection for genome scanning. We demonstrate the proposed method and selection procedure in a backcross population, using both simulated and real data. The extension to other experimental mapping populations is straightforward.

The JAX Diversity Outbred population is a new mouse resource derived from partially inbred Collaborative Cross strains and maintained by randomized outcrossing. As such, it segregates the same allelic variants as the Collaborative Cross but embeds these in a distinct population architecture in which each animal has a high degree of heterozygosity and carries a unique combination of alleles. Phenotypic diversity is striking and often divergent from phenotypes seen in the founder strains of the Collaborative Cross. Allele frequencies and recombination density in early generations of Diversity Outbred mice are consistent with expectations based on simulations of the mating design. We describe analytical methods for genetic mapping using this resource and demonstrate the power and high mapping resolution achieved with this population by mapping a serum cholesterol trait to a 2-Mb region on chromosome 3 containing only 11 genes. Analysis of the estimated allele effects in conjunction with complete genome sequence data of the founder strains reduced the pool of candidate polymorphisms to seven SNPs, five of which are located in an intergenic region upstream of the Foxo1 gene.

Model organisms offer many advantages for the genetic analysis of complex traits. However, identification of specific genes is often hampered by a lack of recombination between the genomes of inbred progenitors. Recently, genome-wide association studies (GWAS) in humans have offered gene-level mapping resolution that is possible because of the large number of accumulated recombinations among unrelated human subjects. To obtain analogous improvements in mapping resolution in mice, we used a 34th generation advanced intercross line (AIL) derived from two inbred strains (SM/J and LG/J). We used simulations to show that familial relationships among subjects must be accounted for when analyzing these data; we then used a mixed model that included polygenic effects to address this problem in our own analysis. Using a combination of F2 and AIL mice derived from the same inbred progenitors, we identified genome-wide significant, subcentimorgan loci that were associated with methamphetamine sensitivity, (e.g., chromosome 18; LOD = 10.5) and non-drug-induced locomotor activity (e.g., chromosome 8; LOD = 18.9). The 2-LOD support interval for the former locus contains no known genes while the latter contains only one gene (Csmd1). This approach is broadly applicable in terms of phenotypes and model organisms and allows GWAS to be performed in multigenerational crosses between and among inbred strains where familial relatedness is often unavoidable.

Climate change has altered life history events in many plant species; however, little is known about genetic variation underlying seasonal thermal response. In this study, we simulated current and three future warming climates and measured flowering time across a globally diverse set of Arabidopsis thaliana accessions. We found that increased diurnal and seasonal temperature (1°–3°) decreased flowering time in two fall cohorts. The early fall cohort was unique in that both rapid cycling and overwintering life history strategies were revealed; the proportion of rapid cycling plants increased by 3–7% for each 1° temperature increase. We performed genome-wide association studies (GWAS) to identify the underlying genetic basis of thermal sensitivity. GWAS identified five main-effect quantitative trait loci (QTL) controlling flowering time and another five QTL with thermal sensitivity. Candidate genes include known flowering loci; a cochaperone that interacts with heat-shock protein 90; and a flowering hormone, gibberellic acid, a biosynthetic enzyme. The identified genetic architecture allowed accurate prediction of flowering phenotypes (R2 > 0.95) that has application for genomic selection of adaptive genotypes for future environments. This work may serve as a reference for breeding and conservation genetic studies under changing environments.

We used simulations to evaluate methods for assessing statistical significance in association studies. When the statistical model appropriately accounted for relatedness among individuals, unrestricted permutation tests and a few other simulation-based methods effectively controlled type I error rates; otherwise, only gene dropping controlled type I error but at the expense of statistical power.