Mapping in structured populations by resample model averaging

Genetic analysis of pericarp pigmentation variation in Corn Belt dent maize

The US standard for maize commercially grown for grain specifies that yellow corn can contain at maximum 5% corn of other colors. Inbred parents of commercial hybrids typically have clear pericarp, but transgressive segregants in breeding populations can display variation in pericarp pigmentation. We identified 10 doubled haploid biparental populations segregating for pigmented pericarp and evaluated qualitative genetic models using chi-square tests of observed and expected frequencies. Pigmentation ranged from light to dark brown color, and pigmentation intensity was quantitatively measured across 1,327 inbred lines using hue calculated from RGB pixel values. Genetic mapping was used to identify loci associated with pigmentation intensity. For 9 populations, pigmentation inheritance best fit a hypothesis of a 2- or 3-gene epistatic model. Significant differences in pigment intensity were observed across populations. W606S-derived inbred lines with the darkest pericarp often had clear glumes, suggesting the presence of a novel P1-rw allele, a hypothesis supported by a significant quantitative trait locus peak at P1. A separate quantitative trait locus region on chromosome 2 between 221.64 and 226.66 Mbp was identified in LH82-derived populations, and the peak near p1 was absent. A genome-wide association study using 416 inbred lines from the Wisconsin Diversity panel with full genome resequencing revealed 4 significant associations including the region near P1. This study supports that pericarp pigmentation among dent maize inbreds can arise by transgressive segregation when pigmentation in the parental generation is absent and is partially explained by functional allelic variation at the P1 locus.

Multi-omics analysis identifies drivers of protein phosphorylation

BACKGROUND

Phosphorylation of proteins is a key step in the regulation of many cellular processes including activation of enzymes and signaling cascades. The abundance of a phosphorylated peptide (phosphopeptide) is determined by the abundance of its parent protein and the proportion of target sites that are phosphorylated.

RESULTS

We quantified phosphopeptides, proteins, and transcripts in heart, liver, and kidney tissue samples of mice from 58 strains of the Collaborative Cross strain panel. We mapped ~700 phosphorylation quantitative trait loci (phQTL) across the three tissues and applied genetic mediation analysis to identify causal drivers of phosphorylation. We identified kinases, phosphatases, cytokines, and other factors, including both known and potentially novel interactions between target proteins and genes that regulate site-specific phosphorylation. Our analysis highlights multiple targets of pyruvate dehydrogenase kinase 1 (PDK1), a regulator of mitochondrial function that shows reduced activity in the NZO/HILtJ mouse, a polygenic model of obesity and type 2 diabetes.

CONCLUSIONS

Together, this integrative multi-omics analysis in genetically diverse CC strains provides a powerful tool to identify regulators of protein phosphorylation. The data generated in this study provides a resource for further exploration.

Genetic Mapping of Multiple Traits Identifies Novel Genes for Adiposity, Lipids, and Insulin Secretory Capacity in Outbred Rats

Despite the successes of human genome-wide association studies, the causal genes underlying most metabolic traits remain unclear. We used outbred heterogeneous stock (HS) rats, coupled with expression data and mediation analysis, to identify quantitative trait loci (QTLs) and candidate gene mediators for adiposity, glucose tolerance, serum lipids, and other metabolic traits. Physiological traits were measured in 1,519 male HS rats, with liver and adipose transcriptomes measured in >410 rats. Genotypes were imputed from low-coverage whole-genome sequencing. Linear mixed models were used to detect physiological and expression QTLs (pQTLs and eQTLs, respectively), using both single nucleotide polymorphism (SNP)- and haplotype-based models for pQTL mapping. Genes with cis-eQTLs that overlapped pQTLs were assessed as causal candidates through mediation analysis. We identified 14 SNP-based pQTLs and 19 haplotype-based pQTLs, of which 10 were in common. Using mediation, we identified the following genes as candidate mediators of pQTLs: Grk5 for fat pad weight and serum triglyceride pQTLs on Chr1, Krtcap3 for fat pad weight and serum triglyceride pQTLs on Chr6, Ilrun for a fat pad weight pQTL on Chr20, and Rfx6 for a whole pancreatic insulin content pQTL on Chr20. Furthermore, we verified Grk5 and Ktrcap3 using gene knockdown/out models, thereby shedding light on novel regulators of obesity.

Architecture of the outbred brown fat proteome defines regulators of metabolic physiology

Brown adipose tissue (BAT) regulates metabolic physiology. However, nearly all mechanistic studies of BAT protein function occur in a single inbred mouse strain, which has limited the understanding of generalizable mechanisms of BAT regulation over physiology. Here, we perform deep quantitative proteomics of BAT across a cohort of 163 genetically defined diversity outbred mice, a model that parallels the genetic and phenotypic variation found in humans. We leverage this diversity to define the functional architecture of the outbred BAT proteome, comprising 10,479 proteins. We assign co-operative functions to 2,578 proteins, enabling systematic discovery of regulators of BAT. We also identify 638 proteins that correlate with protection from, or sensitivity to, at least one parameter of metabolic disease. We use these findings to uncover SFXN5, LETMD1, and ATP1A2 as modulators of BAT thermogenesis or adiposity, and provide OPABAT as a resource for understanding the conserved mechanisms of BAT regulation over metabolic physiology.

Association mapping across a multitude of traits collected in diverse environments in maize

Classical genetic studies have identified many cases of pleiotropy where mutations in individual genes alter many different phenotypes. Quantitative genetic studies of natural genetic variants frequently examine one or a few traits, limiting their potential to identify pleiotropic effects of natural genetic variants. Widely adopted community association panels have been employed by plant genetics communities to study the genetic basis of naturally occurring phenotypic variation in a wide range of traits. High-density genetic marker data-18M markers-from 2 partially overlapping maize association panels comprising 1,014 unique genotypes grown in field trials across at least 7 US states and scored for 162 distinct trait data sets enabled the identification of of 2,154 suggestive marker-trait associations and 697 confident associations in the maize genome using a resampling-based genome-wide association strategy. The precision of individual marker-trait associations was estimated to be 3 genes based on a reference set of genes with known phenotypes. Examples were observed of both genetic loci associated with variation in diverse traits (e.g., above-ground and below-ground traits), as well as individual loci associated with the same or similar traits across diverse environments. Many significant signals are located near genes whose functions were previously entirely unknown or estimated purely via functional data on homologs. This study demonstrates the potential of mining community association panel data using new higher-density genetic marker sets combined with resampling-based genome-wide association tests to develop testable hypotheses about gene functions, identify potential pleiotropic effects of natural genetic variants, and study genotype-by-environment interaction.

Genome-wide association mapping of resistance to the sorghum aphid in Sorghum bicolor

Since 2013, the sorghum aphid (SA), Melanaphis sorghi (Theobald), has been a serious pest that hampers all types of sorghum production in the U.S. Known sorghum aphid resistance in sorghum is limited to a few genetic regions on SBI-06. In this study, a subset of the Sorghum Association Panel (SAP) was used along with some additional lines to identify genomic regions that confer sorghum aphid resistance. SAP lines were grown in the field and visually evaluated for SA resistance during the growing seasons of 2019 and 2020 in Tifton, GA. In 2020, the SAP accessions were also evaluated for SA resistance in the field using drone-based high throughput phenotyping (HTP). Flowering time was recorded in the field to confirm that our methods were sufficient for identifying known quantitative trait loci (QTL). This study combined phenotypic data from field-based visual ratings and reflectance data to identify genome-wide associated (GWAS) marker-trait associations (MTA) using genotyping-by-sequencing (GBS) data. Several MTAs were identified for SA-related traits across the genome, with a few common markers that were consistently identified on SBI-08 and SBI-10 for aphid count and plant damage, as well as loci for reflectance-based traits on SBI-02, SBI-03, and SBI-05. Candidate genes encoding leucine-rich repeats (LRR), Avr proteins, lipoxygenases (LOXs), calmodulins (CAM) dependent protein kinase, WRKY transcription factors, flavonoid biosynthesis genes, and 12-oxo-phytodienoic acid reductase were identified near SNPs that had significant associations with different SA traits. In this study, flowering time-related genes were also identified as a positive control for the methods. The total phenotypic variation explained by significant SNPs across SA-scored traits, reflectance data, and flowering time ranged from 6 to 61%, while the heritability value ranged from 4 to 69%. This study identified three new sources of resistant lines to sorghum aphid. These results supported the existing literature, and also revealed several new loci. Markers identified in this study will support marker-assisted breeding for sorghum aphid resistance.

A Bayesian model selection approach to mediation analysis

Genetic studies often seek to establish a causal chain of events originating from genetic variation through to molecular and clinical phenotypes. When multiple phenotypes share a common genetic association, one phenotype may act as an intermediate for the genetic effects on the other. Alternatively, the phenotypes may be causally unrelated but share genetic loci. Mediation analysis represents a class of causal inference approaches used to determine which of these scenarios is most plausible. We have developed a general approach to mediation analysis based on Bayesian model selection and have implemented it in an R package, bmediatR. Bayesian model selection provides a flexible framework that can be tailored to different analyses. Our approach can incorporate prior information about the likelihood of models and the strength of causal effects. It can also accommodate multiple genetic variants or multi-state haplotypes. Our approach reports posterior probabilities that can be useful in interpreting uncertainty among competing models. We compared bmediatR with other popular methods, including the Sobel test, Mendelian randomization, and Bayesian network analysis using simulated data. We found that bmediatR performed as well or better than these alternatives in most scenarios. We applied bmediatR to proteome data from Diversity Outbred (DO) mice, a multi-parent population, and demonstrate the power of mediation with multi-state haplotypes. We also applied bmediatR to data from human cell lines to identify transcripts that are mediated through or are expressed independently from local chromatin accessibility. We demonstrate that Bayesian model selection provides a powerful and versatile approach to identify causal relationships in genetic studies using model organism or human data.

Simultaneous dissection of grain carotenoid levels and kernel color in biparental maize populations with yellow-to-orange grain

Maize enriched in provitamin A carotenoids could be key in combatting vitamin A deficiency in human populations relying on maize as a food staple. Consumer studies indicate that orange maize may be regarded as novel and preferred. This study identifies genes of relevance for grain carotenoid concentrations and kernel color, through simultaneous dissection of these traits in 10 families of the US maize nested association mapping panel that have yellow to orange grain. Quantitative trait loci were identified via joint-linkage analysis, with phenotypic variation explained for individual kernel color quantitative trait loci ranging from 2.4% to 17.5%. These quantitative trait loci were cross-analyzed with significant marker-trait associations in a genome-wide association study that utilized ∼27 million variants. Nine genes were identified: four encoding activities upstream of the core carotenoid pathway, one at the pathway branchpoint, three within the α- or β-pathway branches, and one encoding a carotenoid cleavage dioxygenase. Of these, three exhibited significant pleiotropy between kernel color and one or more carotenoid traits. Kernel color exhibited moderate positive correlations with β-branch and total carotenoids and negligible correlations with α-branch carotenoids. These findings can be leveraged to simultaneously achieve desirable kernel color phenotypes and increase concentrations of provitamin A and other priority carotenoids.

3D reconstruction identifies loci linked to variation in angle of individual sorghum leaves

Selection for yield at high planting density has reshaped the leaf canopy of maize, improving photosynthetic productivity in high density settings. Further optimization of canopy architecture may be possible. However, measuring leaf angles, the widely studied component trait of leaf canopy architecture, by hand is a labor and time intensive process. Here, we use multiple, calibrated, 2D images to reconstruct the 3D geometry of individual sorghum plants using a voxel carving based algorithm. Automatic skeletonization and segmentation of these 3D geometries enable quantification of the angle of each leaf for each plant. The resulting measurements are both heritable and correlated with manually collected leaf angles. This automated and scaleable reconstruction approach was employed to measure leaf-by-leaf angles for a population of 366 sorghum plants at multiple time points, resulting in 971 successful reconstructions and 3,376 leaf angle measurements from individual leaves. A genome wide association study conducted using aggregated leaf angle data identified a known large effect leaf angle gene, several previously identified leaf angle QTL from a sorghum NAM population, and novel signals. Genome wide association studies conducted separately for three individual sorghum leaves identified a number of the same signals, a previously unreported signal shared across multiple leaves, and signals near the sorghum orthologs of two maize genes known to influence leaf angle. Automated measurement of individual leaves and mapping variants associated with leaf angle reduce the barriers to engineering ideal canopy architectures in sorghum and other grain crops.

Regulation of protein abundance in genetically diverse mouse populations

Genetically diverse mouse populations are powerful tools for characterizing the regulation of the proteome and its relationship to whole-organism phenotypes. We used mass spectrometry to profile and quantify the abundance of 6,798 proteins in liver tissue from mice of both sexes across 58 Collaborative Cross (CC) inbred strains. We previously collected liver proteomics data from the related Diversity Outbred (DO) mice and their founder strains. We show concordance across the proteomics datasets despite being generated from separate experiments, allowing comparative analysis. We map protein abundance quantitative trait loci (pQTLs), identifying 1,087 local and 285 distal in the CC mice and 1,706 local and 414 distal in the DO mice. We find that regulatory effects on individual proteins are conserved across the mouse populations, in particular for local genetic variation and sex differences. In comparison, proteins that form complexes are often co-regulated, displaying varying genetic architectures, and overall show lower heritability and map fewer pQTLs. We have made this resource publicly available to enable quantitative analyses of the regulation of the proteome.

Eleven biosynthetic genes explain the majority of natural variation in carotenoid levels in maize grain

Vitamin A deficiency remains prevalent in parts of Asia, Latin America, and sub-Saharan Africa where maize (Zea mays) is a food staple. Extensive natural variation exists for carotenoids in maize grain. Here, to understand its genetic basis, we conducted a joint linkage and genome-wide association study of the US maize nested association mapping panel. Eleven of the 44 detected quantitative trait loci (QTL) were resolved to individual genes. Six of these were correlated expression and effect QTL (ceeQTL), showing strong correlations between RNA-seq expression abundances and QTL allelic effect estimates across six stages of grain development. These six ceeQTL also had the largest percentage of phenotypic variance explained, and in major part comprised the three to five loci capturing the bulk of genetic variation for each trait. Most of these ceeQTL had strongly correlated QTL allelic effect estimates across multiple traits. These findings provide an in-depth genome-level understanding of the genetic and molecular control of carotenoids in plants. In addition, these findings provide a roadmap to accelerate breeding for provitamin A and other priority carotenoid traits in maize grain that should be readily extendable to other cereals.

Sept8/SEPTIN8 involvement in cellular structure and kidney damage is identified by genetic mapping and a novel human tubule hypoxic model

Chronic kidney disease (CKD), which can ultimately progress to kidney failure, is influenced by genetics and the environment. Genes identified in human genome wide association studies (GWAS) explain only a small proportion of the heritable variation and lack functional validation, indicating the need for additional model systems. Outbred heterogeneous stock (HS) rats have been used for genetic fine-mapping of complex traits, but have not previously been used for CKD traits. We performed GWAS for urinary protein excretion (UPE) and CKD related serum biochemistries in 245 male HS rats. Quantitative trait loci (QTL) were identified using a linear mixed effect model that tested for association with imputed genotypes. Candidate genes were identified using bioinformatics tools and targeted RNAseq followed by testing in a novel in vitro model of human tubule, hypoxia-induced damage. We identified two QTL for UPE and five for serum biochemistries. Protein modeling identified a missense variant within Septin 8 (Sept8) as a candidate for UPE. Sept8/SEPTIN8 expression increased in HS rats with elevated UPE and tubulointerstitial injury and in the in vitro hypoxia model. SEPTIN8 is detected within proximal tubule cells in human kidney samples and localizes with acetyl-alpha tubulin in the culture system. After hypoxia, SEPTIN8 staining becomes diffuse and appears to relocalize with actin. These data suggest a role of SEPTIN8 in cellular organization and structure in response to environmental stress. This study demonstrates that integration of a rat genetic model with an environmentally induced tubule damage system identifies Sept8/SEPTIN8 and informs novel aspects of the complex gene by environmental interactions contributing to CKD risk.

Inferring the Allelic Series at QTL in Multiparental Populations

Multiparental populations (MPPs) are experimental populations in which the genome of every individual is a mosaic of known founder haplotypes. These populations are useful for detecting quantitative trait loci (QTL) because tests of association can leverage inferred founder haplotype descent. It is difficult, however, to determine how haplotypes at a locus group into distinct functional alleles, termed the allelic series. The allelic series is important because it provides information about the number of causal variants at a QTL and their combined effects. In this study, we introduce a fully Bayesian model selection framework for inferring the allelic series. This framework accounts for sources of uncertainty found in typical MPPs, including the number and composition of functional alleles. Our prior distribution for the allelic series is based on the Chinese restaurant process, a relative of the Dirichlet process, and we leverage its connection to the coalescent to introduce additional prior information about haplotype relatedness via a phylogenetic tree. We evaluate our approach via simulation and apply it to QTL from two MPPs: the Collaborative Cross (CC) and the Drosophila Synthetic Population Resource (DSPR). We find that, although posterior inference of the exact allelic series is often uncertain, we are able to distinguish biallelic QTL from more complex multiallelic cases. Additionally, our allele-based approach improves haplotype effect estimation when the true number of functional alleles is small. Our method, Tree-Based Inference of Multiallelism via Bayesian Regression (TIMBR), provides new insight into the genetic architecture of QTL in MPPs.

Container-aided integrative QTL and RNA-seq analysis of Collaborative Cross mice supports distinct sex-oriented molecular modes of response in obesity

BACKGROUND

The Collaborative Cross (CC) mouse population is a valuable resource to study the genetic basis of complex traits, such as obesity. Although the development of obesity is influenced by environmental factors, underlying genetic mechanisms play a crucial role in the response to these factors. The interplay between the genetic background and the gene expression pattern can provide further insight into this response, but we lack robust and easily reproducible workflows to integrate genomic and transcriptomic information in the CC mouse population.

RESULTS

We established an automated and reproducible integrative workflow to analyse complex traits in the CC mouse genetic reference panel at the genomic and transcriptomic levels. We implemented the analytical workflow to assess the underlying genetic mechanisms of host susceptibility to diet induced obesity and integrated these results with diet induced changes in the hepatic gene expression of susceptible and resistant mice. Hepatic gene expression differs significantly between obese and non-obese mice, with a significant sex effect, where male and female mice exhibit different responses and coping mechanisms.

CONCLUSION

Integration of the data showed that different genes but similar pathways are involved in the genetic susceptibility and disturbed in diet induced obesity. Genetic mechanisms underlying susceptibility to high-fat diet induced obesity are different in female and male mice. The clear distinction we observed in the systemic response to the high-fat diet challenge and to obesity between male and female mice points to the need for further research into distinct sex-related mechanisms in metabolic disease.

Genomic mapping in outbred mice reveals overlap in genetic susceptibility for HZE ion- and γ-ray-induced tumors

Cancer risk from galactic cosmic radiation exposure is considered a potential "showstopper" for a manned mission to Mars. Calculating the actual risks confronted by spaceflight crews is complicated by our limited understanding of the carcinogenic effects of high-charge, high-energy (HZE) ions, a radiation type for which no human exposure data exist. Using a mouse model of genetic diversity, we find that the histotype spectrum of HZE ion-induced tumors is similar to the spectra of spontaneous and γ-ray-induced tumors and that the genomic loci controlling susceptibilities overlap between groups for some tumor types. Where it occurs, this overlap indicates shared tumorigenesis mechanisms regardless of the type of radiation exposure and supports the use of human epidemiological data from γ-ray exposures to predict cancer risks from galactic cosmic rays.

Genetic architecture of phenotypic means and plasticities of kernel size and weight in maize

Genetic relationships between the phenotypic means and plasticities of kernel size and weight revealed the common genetic control of these traits in maize. Kernel size and weight are crucial components of grain yield in maize, and phenotypic plasticity in these traits facilitates adaptations to changing environments. Elucidating the genetic architecture of the mean phenotypic values and plasticities of kernel size and weight may be essential for breeding climate-robust maize varieties. Here, a maize nested association mapping (CN-NAM) population and association panel were grown in different environments. A joint linkage analysis and genome-wide association mapping were performed for five kernel size and weight phenotypic traits and two phenotypic plasticity measures. The mean phenotypes and plasticities were significantly correlated. The overall results of quantitative trait locus (QTL) and candidate gene analyses indicated moderate and high levels of common genetic control for the two traits. Furthermore, the mean phenotypes or plasticities of the hundred-kernel weight and volume were commonly regulated to a high degree. One pleiotropic locus on chromosome 10 simultaneously controlled the mean phenotypic values and plasticities of kernel size and weight. Therefore, the plasticity of kernel size and weight might be indirectly selected during maize breeding; however, selecting for high or low plasticity in combination with high or low mean phenotypic values of kernel size and weight traits may be difficult.

Using Heterogeneous Stocks for Fine-Mapping Genetically Complex Traits

In this chapter we will review both the rationale and experimental design for using Heterogeneous Stock (HS) populations for fine-mapping of complex traits in mice and rats. We define an HS as an outbred population derived from an intercross between two or more inbred strains. HS have been used to perform genome-wide association studies (GWAS) for multiple behavioral, physiological, and gene expression traits. GWAS using HS require four key steps, which we review: selection of an appropriate HS population, phenotyping, genotyping, and statistical analysis. We provide advice on the selection of an HS, comment on key issues related to phenotyping, discuss genotyping methods relevant to these populations, and describe statistical genetic analyses that are applicable to genetic analyses that use HS.

Genome-wide association analysis for fumonisin content in maize kernels

BACKGROUND

Plant breeding has been proposed as one of the most effective and environmentally safe methods to control fungal infection and to reduce fumonisin accumulation. However, conventional breeding can be hampered by the complex genetic architecture of resistance to fumonisin accumulation and marker-assisted selection is proposed as an efficient alternative. In the current study, GWAS has been performed for the first time for detecting high-resolution QTL for resistance to fumonisin accumulation in maize kernels complementing published GWAS results for Fusarium ear rot.

RESULTS

Thirty-nine SNPs significantly associated with resistance to fumonisin accumulation in maize kernels were found and clustered into 17 QTL. Novel QTLs for fumonisin content would be at bins 3.02, 5.02, 7.05 and 8.07. Genes with annotated functions probably implicated in resistance to pathogens based on previous studies have been highlighted.

CONCLUSIONS

Breeding approaches to fix favorable functional variants for genes implicated in maize immune response signaling may be especially useful to reduce kernel contamination with fumonisins without significantly interfering in mycelia development and growth and, consequently, in the beneficial endophytic behavior of Fusarium verticillioides.

The Collaborative Cross: A Systems Genetics Resource for Studying Host-Pathogen Interactions

Host genetic variation has a major impact on infectious disease susceptibility. The study of pathogen resistance genes, largely aided by mouse models, has significantly advanced our understanding of infectious disease pathogenesis. The Collaborative Cross (CC), a newly developed multi-parental mouse genetic reference population, serves as a tractable model system to study how pathogens interact with genetically diverse populations. In this review, we summarize progress utilizing the CC as a platform to develop improved models of pathogen-induced disease and to map polymorphic host response loci associated with variation in susceptibility to pathogens.

Genome-wide association studies for yield-related traits in soft red winter wheat grown in Virginia

Grain yield is a trait of paramount importance in the breeding of all cereals. In wheat (Triticum aestivum L.), yield has steadily increased since the Green Revolution, though the current rate of increase is not forecasted to keep pace with demand due to growing world population and increasing affluence. While several genome-wide association studies (GWAS) on yield and related component traits have been performed in wheat, the previous lack of a reference genome has made comparisons between studies difficult. In this study, a GWAS for yield and yield-related traits was carried out on a population of 322 soft red winter wheat lines across a total of four rain-fed environments in the state of Virginia using single-nucleotide polymorphism (SNP) marker data generated by a genotyping-by-sequencing (GBS) protocol. Two separate mixed linear models were used to identify significant marker-trait associations (MTAs). The first was a single-locus model utilizing a leave-one-chromosome-out approach to estimating kinship. The second was a sub-setting kinship estimation multi-locus method (FarmCPU). The single-locus model identified nine significant MTAs for various yield-related traits, while the FarmCPU model identified 74 significant MTAs. The availability of the wheat reference genome allowed for the description of MTAs in terms of both genetic and physical positions, and enabled more extensive post-GWAS characterization of significant MTAs. The results indicate a number of promising candidate genes contributing to grain yield, including an ortholog of the rice aberrant panicle organization (APO1) protein and a gibberellin oxidase protein (GA2ox-A1) affecting the trait grains per square meter, an ortholog of the Arabidopsis thaliana mother of flowering time and terminal flowering 1 (MFT) gene affecting the trait seeds per square meter, and a B2 heat stress response protein affecting the trait seeds per head.

A Genome Wide Association Study Reveals Markers and Genes Associated with Resistance to Fusarium verticillioides Infection of Seedlings in a Maize Diversity Panel

Fusarium verticillioides infects maize, causing ear rot, yield loss and contamination by fumonisin mycotoxins. The fungus can be transmitted via kernels and cause systemic infection in maize. Maize resistance to the fungus may occur at different developmental stages, from seedling to maturity. Resistance during kernel germination is part of the plant-pathogen interaction and so far this aspect has not been investigated. In the present study, a genome wide association study (GWAS) of resistance to Fusarium during the seedling developmental stage was conducted in a maize diversity panel using 226,446 SNP markers. Seedling germination and disease phenotypes were scored on artificially inoculated kernels using the rolled towel assay. GWAS identified 164 SNPs significantly associated with the traits examined. Four SNPs were associated with disease severity score after inoculation, 153 were associated with severity in asymptomatic kernels and 7 with the difference between the severity ratings in inoculated and non-inoculated seeds. A set of genes containing or physically near the significant SNPs were identified as candidates for Fusarium resistance at the seedling stage. Functional analysis revealed that many of these genes are directly involved in plant defense against pathogens and stress responses, including transcription factors, chitinase, cytochrome P450, and ubiquitination proteins. In addition, 25 genes were found in high linkage disequilibrium with the associated SNPs identified by GWAS and four of them directly involved in disease resistance. These findings contribute to understanding the complex system of maize-F. verticillioides and may improve genomic selection for Fusarium resistance at the seedling stage.

Integrating Coexpression Networks with GWAS to Prioritize Causal Genes in Maize

Genome-wide association studies (GWAS) have identified loci linked to hundreds of traits in many different species. Yet, because linkage equilibrium implicates a broad region surrounding each identified locus, the causal genes often remain unknown. This problem is especially pronounced in nonhuman, nonmodel species, where functional annotations are sparse and there is frequently little information available for prioritizing candidate genes. We developed a computational approach, Camoco, that integrates loci identified by GWAS with functional information derived from gene coexpression networks. Using Camoco, we prioritized candidate genes from a large-scale GWAS examining the accumulation of 17 different elements in maize (Zea mays) seeds. Strikingly, we observed a strong dependence in the performance of our approach based on the type of coexpression network used: expression variation across genetically diverse individuals in a relevant tissue context (in our case, roots that are the primary elemental uptake and delivery system) outperformed other alternative networks. Two candidate genes identified by our approach were validated using mutants. Our study demonstrates that coexpression networks provide a powerful basis for prioritizing candidate causal genes from GWAS loci but suggests that the success of such strategies can highly depend on the gene expression data context. Both the software and the lessons on integrating GWAS data with coexpression networks generalize to species beyond maize.

ComPaSS-GWAS: A method to reduce type I error in genome-wide association studies when replication data are not available

Results from association studies are traditionally corroborated by replicating the findings in an independent data set. Although replication studies may be comparable for the main trait or phenotype of interest, it is unlikely that secondary phenotypes will be comparable across studies, making replication problematic. Alternatively, there may simply not be a replication sample available because of the nature or frequency of the phenotype. In these situations, an approach based on complementary pairs stability selection for genome-wide association study (ComPaSS-GWAS), is proposed as an ad-hoc alternative to replication. In this method, the sample is randomly split into two conditionally independent halves multiple times (resamples) and a GWAS is performed on each half in each resample. Similar in spirit to testing for association with independent discovery and replication samples, a marker is corroborated if its p-value is significant in both halves of the resample. Simulation experiments were performed for both nongenetic and genetic models. The type I error rate and power of ComPaSS-GWAS were determined and compared to the statistical properties of a traditional GWAS. Simulation results show that the type I error rate decreased as the number of resamples increased with only a small reduction in power and that these results were comparable with those from a traditional GWAS. Blood levels of vitamin pyridoxal 5'-phosphate from the Trinity Student Study (TSS) were used to validate this approach. The results from the validation study were compared to, and were consistent with, those obtained from previously published independent replication data and functional studies.

Genome-wide mapping of quantitative trait loci in admixed populations using mixed linear model and Bayesian multiple regression analysis

Background

Population stratification and cryptic relationships have been the main sources of excessive false-positives and false-negatives in population-based association studies. Many methods have been developed to model these confounding factors and minimize their impact on the results of genome-wide association studies. In most of these methods, a two-stage approach is applied where: (1) methods are used to determine if there is a population structure in the sample dataset and (2) the effects of population structure are corrected either by modeling it or by running a separate analysis within each sub-population. The objective of this study was to evaluate the impact of population structure on the accuracy and power of genome-wide association studies using a Bayesian multiple regression method.

Methods

We conducted a genome-wide association study in a stochastically simulated admixed population. The genome was composed of six chromosomes, each with 1000 markers. Fifteen segregating quantitative trait loci contributed to the genetic variation of a quantitative trait with heritability of 0.30. The impact of genetic relationships and breed composition (BC) on three analysis methods were evaluated: single marker simple regression (SMR), single marker mixed linear model (MLM) and Bayesian multiple-regression analysis (BMR). Each method was fitted with and without BC. Accuracy, power, false-positive rate and the positive predictive value of each method were calculated and used for comparison.

Results

SMR and BMR, both without BC, were ranked as the worst and the best performing approaches, respectively. Our results showed that, while explicit modeling of genetic relationships and BC is essential for models SMR and MLM, BMR can disregard them and yet result in a higher power without compromising its false-positive rate.

Conclusions

This study showed that the Bayesian multiple-regression analysis is robust to population structure and to relationships among study subjects and performs better than a single marker mixed linear model approach.

Heterogeneous Stock Populations for Analysis of Complex Traits

Heterogeneous Stock (HS) populations allow for fine-resolution genetic mapping of a variety of complex traits. HS mice and rats were created from breeding together eight inbred strains, followed by maintaining the colony in a manner that minimizes inbreeding. After 50 or more generations of breeding, the resulting animals' chromosomes represent a genetic mosaic of the founders' haplotypes, with the average distance between recombination events in the centiMorgan range. This allows for genetic mapping to only a few Mb, a much smaller region than what can be identified using traditional F2 intercross or backcross mapping strategies. HS animals have been used to fine-map a variety of complex traits including anxiety and fear behaviors, diabetes, asthma, and heart disease, among others. Once a quantitative trait locus (QTL) has been identified, founder sequence and expression analysis can be used to identify underlying causal genes. In the following review, we provide an overview of how HS rats and mice have been used to identify genetic loci, and in some cases the causal genes, underlying complex traits. We discuss the creation and breeding strategies for both HS rats and mice. We then discuss the statistical analyses used to identify genetic loci, as well as strategies to identify causal genes underlying these loci. We end the chapter by discussing limitations faced when using HS populations, including several statistical challenges that have not been fully resolved.

Genetic Fine-Mapping and Identification of Candidate Genes and Variants for Adiposity Traits in Outbred Rats

OBJECTIVE

Obesity is a major risk factor for multiple diseases and is in part heritable, yet the majority of causative genetic variants that drive excessive adiposity remain unknown. Here, outbred heterogeneous stock (HS) rats were used in controlled environmental conditions to fine-map novel genetic modifiers of adiposity.

METHODS

Body weight and visceral fat pad weights were measured in male HS rats that were also genotyped genome-wide. Quantitative trait loci (QTL) were identified by genome-wide association of imputed single-nucleotide polymorphism (SNP) genotypes using a linear mixed effect model that accounts for unequal relatedness between the HS rats. Candidate genes were assessed by protein modeling and mediation analysis of expression for coding and noncoding variants, respectively.

RESULTS

HS rats exhibited large variation in adiposity traits, which were highly heritable and correlated with metabolic health. Fine-mapping of fat pad weight and body weight revealed three QTL and prioritized five candidate genes. Fat pad weight was associated with missense SNPs in Adcy3 and Prlhr and altered expression of Krtcap3 and Slc30a3, whereas Grid2 was identified as a candidate within the body weight locus.

CONCLUSIONS

These data demonstrate the power of HS rats for identification of known and novel heritable mediators of obesity traits.

Novel Loci Underlie Natural Variation in Vitamin E Levels in Maize Grain

Tocopherols, tocotrienols, and plastochromanols (collectively termed tocochromanols) are lipid-soluble antioxidants synthesized by all plants. Their dietary intake, primarily from seed oils, provides vitamin E and other health benefits. Tocochromanol biosynthesis has been dissected in the dicot Arabidopsis thaliana, which has green, photosynthetic seeds, but our understanding of tocochromanol accumulation in major crops, whose seeds are nonphotosynthetic, remains limited. To understand the genetic control of tocochromanols in grain, we conducted a joint linkage and genome-wide association study in the 5000-line U.S. maize (Zea mays) nested association mapping panel. Fifty-two quantitative trait loci for individual and total tocochromanols were identified, and of the 14 resolved to individual genes, six encode novel activities affecting tocochromanols in plants. These include two chlorophyll biosynthetic enzymes that explain the majority of tocopherol variation, which was not predicted given that, like most major cereal crops, maize grain is nonphotosynthetic. This comprehensive assessment of natural variation in vitamin E levels in maize establishes the foundation for improving tocochromanol and vitamin E content in seeds of maize and other major cereal crops.

Identification of trans Protein QTL for Secreted Airway Mucins in Mice and a Causal Role for Bpifb1

Mucus hyper-secretion is a hallmark feature of asthma and other muco-obstructive airway diseases. The mucin proteins MUC5AC and MUC5B are the major glycoprotein components of mucus and have critical roles in airway defense. Despite the biomedical importance of these two proteins, the loci that regulate them in the context of natural genetic variation have not been studied. To identify genes that underlie variation in airway mucin levels, we performed genetic analyses in founder strains and incipient lines of the Collaborative Cross (CC) in a house dust mite mouse model of asthma. CC founder strains exhibited significant differences in MUC5AC and MUC5B, providing evidence of heritability. Analysis of gene and protein expression of Muc5ac and Muc5b in incipient CC lines (n = 154) suggested that post-transcriptional events were important regulators of mucin protein content in the airways. Quantitative trait locus (QTL) mapping identified distinct, trans protein QTL for MUC5AC (chromosome 13) and MUC5B (chromosome 2). These two QTL explained 18 and 20% of phenotypic variance, respectively. Examination of the MUC5B QTL allele effects and subsequent phylogenetic analysis allowed us to narrow the MUC5B QTL and identify Bpifb1 as a candidate gene. Bpifb1 mRNA and protein expression were upregulated in parallel to MUC5B after allergen challenge, and Bpifb1 knockout mice exhibited higher MUC5B expression. Thus, BPIFB1 is a novel regulator of MUC5B.

Genome-wide dissection of heterosis for yield traits in two-line hybrid rice populations

Heterosis has been widely utilized in agriculture and is important for world food safety. Many genetic models have been proposed as mechanisms underlying heterosis during the past century, yet more evidence is needed to support such models. To investigate heterosis in two-line hybrid rice, we generated a partial diallel crossing scheme, which consisted of approximately 500 F1 hybrids derived from 14 male sterile lines and 39 restorer lines. In this population, increased panicle number played the most important role in yield heterosis of hybrid rice. Genome-wide association studies identified many QTLs related to the yield traits of F1 hybrids, better paternal heterosis and special combining ability. Relevant genes, including Hd3a, qGL3, OsmiR156h, and LAX2, were identified as candidates within these QTLs. Nearly forty percent of the QTLs had only two genotypes in the F1 hybrids, mainly because the maternal lines were under intense selective pressure. Further analysis found male sterile lines and restorer lines made different superior allele contributions to F1 hybrids, and their contributions varied among different traits. These results extend our understanding of the molecular basis of heterosis in two-line hybrid rice.

Maize YABBY Genes drooping leaf1 and drooping leaf2 Regulate Plant Architecture

Leaf architecture directly influences canopy structure, consequentially affecting yield. We discovered a maize (Zea mays) mutant with aberrant leaf architecture, which we named drooping leaf1 (drl1). Pleiotropic mutations in drl1 affect leaf length and width, leaf angle, and internode length and diameter. These phenotypes are enhanced by natural variation at the drl2 enhancer locus, including reduced expression of the drl2-Mo17 allele in the Mo17 inbred. A second drl2 allele, produced by transposon mutagenesis, interacted synergistically with drl1 mutants and reduced drl2 transcript levels. The drl genes are required for proper leaf patterning, development and cell proliferation of leaf support tissues, and for restricting auricle expansion at the midrib. The paralogous loci encode maize CRABS CLAW co-orthologs in the YABBY family of transcriptional regulators. The drl genes are coexpressed in incipient and emergent leaf primordia at the shoot apex, but not in the vegetative meristem or stem. Genome-wide association studies using maize NAM-RIL (nested association mapping-recombinant inbred line) populations indicated that the drl loci reside within quantitative trait locus regions for leaf angle, leaf width, and internode length and identified rare single nucleotide polymorphisms with large phenotypic effects for the latter two traits. This study demonstrates that drl genes control the development of key agronomic traits in maize.

Novel Resampling Improves Statistical Power for Multiple-Trait QTL Mapping

Multiple-trait analysis typically employs models that associate a quantitative trait locus (QTL) with all of the traits. As a result, statistical power for QTL detection may not be optimal if the QTL contributes to the phenotypic variation in only a small proportion of the traits. Excluding QTL effects that contribute little to the test statistic can improve statistical power. In this article, we show that an optimal power can be achieved when the number of QTL effects is best estimated, and that a stringent criterion for QTL effect selection may improve power when the number of QTL effects is small but can reduce power otherwise. We investigate strategies for excluding trivial QTL effects, and propose a method that improves statistical power when the number of QTL effects is relatively small, and fairly maintains the power when the number of QTL effects is large. The proposed method first uses resampling techniques to determine the number of nontrivial QTL effects, and then selects QTL effects by the backward elimination procedure for significance test. We also propose a method for testing QTL-trait associations that are desired for biological interpretation in applications. We validate our methods using simulations and Arabidopsis thaliana transcript data.

Enhancing genomic prediction with genome-wide association studies in multiparental maize populations

Genome-wide association mapping using dense marker sets has identified some nucleotide variants affecting complex traits that have been validated with fine-mapping and functional analysis. However, many sequence variants associated with complex traits in maize have small effects and low repeatability. In contrast to genome-wide association study (GWAS), genomic prediction (GP) is typically based on models incorporating information from all available markers, rather than modeling effects of individual loci. We considered methods to integrate results of GWASs into GP models in the context of multiple interconnected families. We compared association tests based on a biallelic additive model constraining the effect of a single-nucleotide polymorphism (SNP) to be equal across all families in which it segregates to a model in which the effect of a SNP can vary across families. Association SNPs were then included as fixed effects into a GP model that also included the random effects of the whole genome background. Simulation studies revealed that the effectiveness of this joint approach depends on the extent of polygenicity of the traits. Congruent with this finding, cross-validation studies indicated that GP including the fixed effects of the most significantly associated SNPs along with the polygenic background was more accurate than the polygenic background model alone for moderately complex but not highly polygenic traits measured in the maize nested association mapping population. Individual SNPs with strong and robust association signals can effectively improve GP. Our approach provides a new integrative modeling approach for both reliable gene discovery and robust GP.

Imputation-Based Fine-Mapping Suggests That Most QTL in an Outbred Chicken Advanced Intercross Body Weight Line Are Due to Multiple, Linked Loci

The Virginia chicken lines have been divergently selected for juvenile body weight for more than 50 generations. Today, the high- and low-weight lines show a >12-fold difference for the selected trait, 56-d body weight. These lines provide unique opportunities to study the genetic architecture of long-term, single-trait selection. Previously, several quantitative trait loci (QTL) contributing to weight differences between the lines were mapped in an F₂-cross between them, and these were later replicated and fine-mapped in a nine-generation advanced intercross of them. Here, we explore the possibility to further increase the fine-mapping resolution of these QTL via a pedigree-based imputation strategy that aims to better capture the genetic diversity in the divergently selected, but outbred, founder lines. The founders of the intercross were high-density genotyped, and then pedigree-based imputation was used to assign genotypes throughout the pedigree. Imputation increased the marker density 20-fold in the selected QTL, providing 6911 markers for the subsequent analysis. Both single-marker association and multi-marker backward-elimination analyses were used to explore regions associated with 56-d body weight. The approach revealed several statistically and population structure independent associations and increased the mapping resolution. Further, most QTL were also found to contain multiple independent associations to markers that were not fixed in the founder populations, implying a complex underlying architecture due to the combined effects of multiple, linked loci perhaps located on independent haplotypes that still segregate in the selected lines.

Estimating parent-specific QTL effects through cumulating linked identity-by-state SNP effects in multiparental populations

The emergence of multiparental mapping populations enabled plant geneticists to gain deeper insights into the genetic architecture of major agronomic traits and to map quantitative trait loci (QTLs) controlling the expression of these traits. Although the investigated mapping populations are similar, one open question is whether genotype data should be modelled as identical by state (IBS) or identical by descent (IBD). Whereas IBS simply makes use of raw genotype scores to distinguish alleles, IBD data are derived from parental offspring information. We report on comparing IBS and IBD by applying two multiple regression models on four traits studied in the barley nested association mapping (NAM) population HEB-25. We observed that modelling parent-specific IBD genotypes produced a lower number of significant QTLs with increased prediction abilities compared with modelling IBS genotypes. However, at lower trait heritabilities the IBS model produced higher prediction abilities. We developed a method to estimate multiallelic QTL effects in multiparental populations from simple biallelic IBS data. This method is based on cumulating IBS-derived single-nucleotide polymorphism (SNP) effect estimates in a defined genetic region surrounding a QTL. Comparing the resulting parent-specific QTL effects with those obtained from IBD approaches revealed high accordance that could be confirmed through simulations. The method turned out to be also applicable to a barley multiparent advanced generation inter-cross (MAGIC) population. The 'cumulation method' represents a universal approach to differentiate parent-specific QTL effects in multiparental populations, even if no IBD information is available. In future, the method could further benefit from the availability of much denser SNP maps.

A multi-marker association method for genome-wide association studies without the need for population structure correction

All common genome-wide association (GWA) methods rely on population structure correction, to avoid false genotype-to-phenotype associations. However, population structure correction is a stringent penalization, which also impedes identification of real associations. Using recent statistical advances, we developed a new GWA method, called Quantitative Trait Cluster Association Test (QTCAT), enabling simultaneous multi-marker associations while considering correlations between markers. With this, QTCAT overcomes the need for population structure correction and also reflects the polygenic nature of complex traits better than single-marker methods. Using simulated data, we show that QTCAT clearly outperforms linear mixed model approaches. Moreover, using QTCAT to reanalyse public human, mouse and Arabidopsis GWA data revealed nearly all known and some previously undetected associations. Following up on the most significant novel association in the Arabidopsis data allowed us to identify a so far unknown component of root growth.

Numerous genetic loci identified for drought tolerance in the maize nested association mapping populations

Background

Maize requires more water than most other crops; therefore, the water use efficiency of this crop must be improved for maize production under undesirable land and changing environmental conditions.

Results

To elucidate the genetic control of drought in maize, we evaluated approximately 5000 inbred lines from 30 linkage-association joint mapping populations under two contrasting water regimes for seven drought-related traits, including yield and anthesis-silking interval (ASI). The joint linkage analysis was conducted to identify 220 quantitative trait loci (QTLs) under well-watered conditions and 169 QTLs under water-stressed conditions. The genome-wide association analysis identified 365 single nucleotide polymorphisms (SNPs) associated with drought-related traits, and these SNPs were located in 354 candidate genes. Fifty-two of these genes showed significant differential expression in the inbred line B73 under the well-watered and water-stressed conditions. In addition, genomic predictions suggested that the moderate-density SNPs obtained through genotyping-by-sequencing were able to make accurate predictions in the nested association mapping population for drought-related traits with moderate-to-high heritability under the water-stressed conditions.

Conclusions

The results of the present study provide important information that can be used to understand the genetic basis of drought stress responses and facilitate the use of beneficial alleles for the improvement of drought tolerance in maize.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3170-8) contains supplementary material, which is available to authorized users.

The Genetics of Leaf Flecking in Maize and Its Relationship to Plant Defense and Disease Resistance

Physiological leaf spotting, or flecking, is a mild-lesion phenotype observed on the leaves of several commonly used maize (Zea mays) inbred lines and has been anecdotally linked to enhanced broad-spectrum disease resistance. Flecking was assessed in the maize nested association mapping (NAM) population, comprising 4,998 recombinant inbred lines from 25 biparental families, and in an association population, comprising 279 diverse maize inbreds. Joint family linkage analysis was conducted with 7,386 markers in the NAM population. Genome-wide association tests were performed with 26.5 million single-nucleotide polymorphisms (SNPs) in the NAM population and with 246,497 SNPs in the association population, resulting in the identification of 18 and three loci associated with variation in flecking, respectively. Many of the candidate genes colocalizing with associated SNPs are similar to genes that function in plant defense response via cell wall modification, salicylic acid- and jasmonic acid-dependent pathways, redox homeostasis, stress response, and vesicle trafficking/remodeling. Significant positive correlations were found between increased flecking, stronger defense response, increased disease resistance, and increased pest resistance. A nonlinear relationship with total kernel weight also was observed whereby lines with relatively high levels of flecking had, on average, lower total kernel weight. We present evidence suggesting that mild flecking could be used as a selection criterion for breeding programs trying to incorporate broad-spectrum disease resistance.

Genome-wide association analysis of forage quality in maize mature stalk

BACKGROUND

Plant digestibility of silage maize (Zea mays L.) has a large influence on nutrition intake for animal feeding. Improving forage quality will enhance the utilization efficiency and feeding value of forage maize. Dissecting the genetic basis of forage quality will improve our understanding of the complex nature of cell wall biosynthesis and degradation, which is also helpful for breeding good quality silage maize.

RESULTS

Acid detergent fiber (ADF), neutral detergent fiber (NDF) and in vitro dry matter digestibility (IVDMD) of stalk were evaluated in a diverse maize population, which is comprised of 368 inbred lines and planted across seven environments. Using a mixed model accounting for population structure and polygenic background effects, a genome-wide association study was conducted to identify single nucleotide polymorphisms (SNPs) significantly associated with forage quality. Scanning 559,285 SNPs across the whole genome, 73, 41 and 82 SNPs were found to be associated with ADF, NDF, and IVDMD, respectively. Each significant SNP explained 4.2 %-6.2 % of the phenotypic variation. Underlying these associated loci, 56 genes were proposed as candidate genes for forage quality.

CONCLUSIONS

Of all the candidate genes proposed by GWAS, we only found a C3H gene (ZmC3H2) that is directly involved in cell wall component biosynthesis. The candidate genes found in this study are mainly involved in signal transduction, stress resistance, and transcriptional regulation of cell wall biosynthetic gene expression. Adding high digestibility maize into the association panel would be helpful for increasing genetic variability and identifying more genes associated with forage quality traits. Cloning and functional validation of these genes would be helpful for understanding the molecular mechanism of the fiber content and digestibility. These findings provide us new insights into cell wall formation and deposition.

Characterization of Biosynthetic Pathways for the Production of the Volatile Homoterpenes DMNT and TMTT in Zea mays

Plant volatiles not only have multiple defense functions against herbivores, fungi, and bacteria, but also have been implicated in signaling within the plant and toward other organisms. Elucidating the function of individual plant volatiles will require more knowledge of their biosynthesis and regulation in response to external stimuli. By exploiting the variation of herbivore-induced volatiles among 26 maize (Zea mays) inbred lines, we conducted a nested association mapping and genome-wide association study (GWAS) to identify a set of quantitative trait loci (QTLs) for investigating the pathways of volatile terpene production. The most significant identified QTL affects the emission of (E)-nerolidol, linalool, and the two homoterpenes (E)-3,8-dimethyl-1,4,7-nonatriene (DMNT) and (E,E)-4,8,12-trimethyltrideca-1,3,7,11-tetraene (TMTT). GWAS associated a single nucleotide polymorphism in the promoter of the gene encoding the terpene synthase TPS2 with this QTL Biochemical characterization of TPS2 verified that this plastid-localized enzyme forms linalool, (E)-nerolidol, and (E,E)-geranyllinalool. The subsequent conversion of (E)-nerolidol into DMNT maps to a P450 monooxygenase, CYP92C5, which is capable of converting nerolidol into DMNT by oxidative degradation. A QTL influencing TMTT accumulation corresponds to a similar monooxygenase, CYP92C6, which is specific for the conversion of (E,E)-geranyllinalool to TMTT The DMNT biosynthetic pathway and both monooxygenases are distinct from those previously characterized for DMNT and TMTT synthesis in Arabidopsis thaliana, suggesting independent evolution of these enzymatic activities.

A Genome-Wide Association Study for Partial Resistance to Maize Common Rust

Quantitative resistance to maize common rust (causal agent Puccinia sorghi) was assessed in an association mapping population of 274 diverse inbred lines. Resistance to common rust was found to be moderately correlated with resistance to three other diseases and with the severity of the hypersensitive defense response previously assessed in the same population. Using a mixed linear model accounting for the confounding effects of population structure and flowering time, genome-wide association tests were performed based at 246,497 single-nucleotide polymorphism loci. Three loci associated with maize common rust resistance were identified. Candidate genes at each locus had predicted roles, mainly in cell wall modification. Other candidate genes included a resistance gene and a gene with a predicted role in regulating accumulation of reactive oxygen species.

Genome-wide dissection of the maize ear genetic architecture using multiple populations

Improvement of grain yield is an essential long-term goal of maize (Zea mays) breeding to meet continual and increasing food demands worldwide, but the genetic basis remains unclear. We used 10 different recombination inbred line (RIL) populations genotyped with high-density markers and phenotyped in multiple environments to dissect the genetic architecture of maize ear traits. Three methods were used to map the quantitative trait loci (QTLs) affecting ear traits. We found 17-34 minor- or moderate-effect loci that influence ear traits, with little epistasis and environmental interactions, totally accounting for 55.4-82% of the phenotypic variation. Four novel QTLs were validated and fine mapped using candidate gene association analysis, expression QTL analysis and heterogeneous inbred family validation. The combination of multiple different populations is a flexible and manageable way to collaboratively integrate widely available genetic resources, thereby boosting the statistical power of QTL discovery for important traits in agricultural crops, ultimately facilitating breeding programs.

Increasing the power of genome wide association studies in natural populations using repeated measures - evaluation and implementation

Genomewide association studies (GWAS) enable detailed dissections of the genetic basis for organisms' ability to adapt to a changing environment. In long‐term studies of natural populations, individuals are often marked at one point in their life and then repeatedly recaptured. It is therefore essential that a method for GWAS includes the process of repeated sampling. In a GWAS, the effects of thousands of single‐nucleotide polymorphisms (SNPs) need to be fitted and any model development is constrained by the computational requirements. A method is therefore required that can fit a highly hierarchical model and at the same time is computationally fast enough to be useful.

Our method fits fixed SNP effects in a linear mixed model that can include both random polygenic effects and permanent environmental effects. In this way, the model can correct for population structure and model repeated measures. The covariance structure of the linear mixed model is first estimated and subsequently used in a generalized least squares setting to fit the SNP effects. The method was evaluated in a simulation study based on observed genotypes from a long‐term study of collared flycatchers in Sweden.

The method we present here was successful in estimating permanent environmental effects from simulated repeated measures data. Additionally, we found that especially for variable phenotypes having large variation between years, the repeated measurements model has a substantial increase in power compared to a model using average phenotypes as a response.

The method is available in the r package RepeatABEL. It increases the power in GWAS having repeated measures, especially for long‐term studies of natural populations, and the R implementation is expected to facilitate modelling of longitudinal data for studies of both animal and human populations.

Positional Gene Cloning in Experimental Populations

Positional cloning is a technique that identifies a trait-associated gene based on its location in the genome and involves methods such as linkage analysis, association mapping, and bioinformatics. This approach can be used for gene identification even when little is known about the molecular basis of the trait. Vast majority of traits are regulated by multiple genomic loci called quantitative trait loci (QTL). We describe experimental populations and designs that can be used for positional cloning, including backcrosses, intercrosses, and heterogeneous stocks, and advantages and disadvantages of different approaches. Once the phenotype and genotype of each individual in an experimental population have been determined, QTL identification can be accomplished. We describe the statistical tools used to identify the existence, location, and significance of QTLs. These different methods have advantages and disadvantages to consider when selecting the appropriate model to be used, which is briefly discussed.Although the objective of QTL mapping is to identify genomic regions associated with a trait, the ultimate goal is to identify the gene and the genetic variation (which is often quantitative trait nucleotide, QTN) or haplotype that is responsible for the phenotype. By discovering the function of causative variants or haplotypes we can understand the molecular changes that lead to the phenotype. We briefly describe how the genomic sequences can be exploited to identify QTNs and how these can be validated in congenic strains and functionally tested to understand their influence on phenotype expression.

Genome Wide Identification of SARS-CoV Susceptibility Loci Using the Collaborative Cross

New systems genetics approaches are needed to rapidly identify host genes and genetic networks that regulate complex disease outcomes. Using genetically diverse animals from incipient lines of the Collaborative Cross mouse panel, we demonstrate a greatly expanded range of phenotypes relative to classical mouse models of SARS-CoV infection including lung pathology, weight loss and viral titer. Genetic mapping revealed several loci contributing to differential disease responses, including an 8.5Mb locus associated with vascular cuffing on chromosome 3 that contained 23 genes and 13 noncoding RNAs. Integrating phenotypic and genetic data narrowed this region to a single gene, Trim55, an E3 ubiquitin ligase with a role in muscle fiber maintenance. Lung pathology and transcriptomic data from mice genetically deficient in Trim55 were used to validate its role in SARS-CoV-induced vascular cuffing and inflammation. These data establish the Collaborative Cross platform as a powerful genetic resource for uncovering genetic contributions of complex traits in microbial disease severity, inflammation and virus replication in models of outbred populations.

Standing genetic variation as a major contributor to adaptation in the Virginia chicken lines selection experiment

Background

Artificial selection provides a powerful approach to study the genetics of adaptation. Using selective-sweep mapping, it is possible to identify genomic regions where allele-frequencies have diverged during selection. To avoid false positive signatures of selection, it is necessary to show that a sweep affects a selected trait before it can be considered adaptive. Here, we confirm candidate, genome-wide distributed selective sweeps originating from the standing genetic variation in a long-term selection experiment on high and low body weight of chickens.

Results

Using an intercross between the two divergent chicken lines, 16 adaptive selective sweeps were confirmed based on their association with the body weight at 56 days of age. Although individual additive effects were small, the fixation for alternative alleles across the loci contributed at least 40 % of the phenotypic difference for the selected trait between these lines. The sweeps contributed about half of the additive genetic variance present within and between the lines after 40 generations of selection, corresponding to a considerable portion of the additive genetic variance of the base population.

Conclusions

Long-term, single-trait, bi-directional selection in the Virginia chicken lines has resulted in a gradual response to selection for extreme phenotypes without a drastic reduction in the genetic variation. We find that fixation of several standing genetic variants across a highly polygenic genetic architecture made a considerable contribution to long-term selection response. This provides new fundamental insights into the dynamics of standing genetic variation during long-term selection and adaptation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0785-z) contains supplementary material, which is available to authorized users.

Identification of microRNAs associated with allergic airway disease using a genetically diverse mouse population

Background

Allergic airway diseases (AADs) such as asthma are characterized in part by granulocytic airway inflammation. The gene regulatory networks that govern granulocyte recruitment are poorly understood, but evidence is accruing that microRNAs (miRNAs) play an important role. To identify miRNAs that may underlie AADs, we used two complementary approaches that leveraged the genotypic and phenotypic diversity of the Collaborative Cross (CC) mouse population. In the first approach, we sought to identify miRNA expression quantitative trait loci (eQTL) that overlap QTL for AAD-related phenotypes. Specifically, CC founder strains and incipient lines of the CC were sensitized and challenged with house dust mite allergen followed by measurement of granulocyte recruitment to the lung. Total lung RNA was isolated and miRNA was measured using arrays for CC founders and qRT-PCR for incipient CC lines.

Results

Among CC founders, 92 miRNAs were differentially expressed. We measured the expression of 40 of the most highly expressed of these 92 miRNAs in the incipient lines of the CC and identified 18 eQTL corresponding to 14 different miRNAs. Surprisingly, half of these eQTL were distal to the corresponding miRNAs, and even on different chromosomes. One of the largest-effect local miRNA eQTL was for miR-342-3p, for which we identified putative causal variants by bioinformatic analysis of the effects of single nucleotide polymorphisms on RNA structure. None of the miRNA eQTL co-localized with QTL for eosinophil or neutrophil recruitment. In the second approach, we constructed putative miRNA/mRNA regulatory networks and identified three miRNAs (miR-497, miR-351 and miR-31) as candidate master regulators of genes associated with neutrophil recruitment. Analysis of a dataset from human keratinocytes transfected with a miR-31 inhibitor revealed two target genes in common with miR-31 targets correlated with neutrophils, namely Oxsr1 and Nsf.

Conclusions

miRNA expression in the allergically inflamed murine lung is regulated by genetic loci that are smaller in effect size compared to mRNA eQTL and often act in trans. Thus our results indicate that the genetic architecture of miRNA expression is different from mRNA expression. We identified three miRNAs, miR-497, miR-351 and miR-31, that are candidate master regulators of genes associated with neutrophil recruitment. Because miR-31 is expressed in airway epithelia and is predicted to target genes with known links to neutrophilic inflammation, we suggest that miR-31 is a potentially novel regulator of airway inflammation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1732-9) contains supplementary material, which is available to authorized users.

Ensemble Learning of QTL Models Improves Prediction of Complex Traits

Quantitative trait locus (QTL) models can provide useful insights into trait genetic architecture because of their straightforward interpretability but are less useful for genetic prediction because of the difficulty in including the effects of numerous small effect loci without overfitting. Tight linkage between markers introduces near collinearity among marker genotypes, complicating the detection of QTL and estimation of QTL effects in linkage mapping, and this problem is exacerbated by very high density linkage maps. Here we developed a thinning and aggregating (TAGGING) method as a new ensemble learning approach to QTL mapping. TAGGING reduces collinearity problems by thinning dense linkage maps, maintains aspects of marker selection that characterize standard QTL mapping, and by ensembling, incorporates information from many more markers-trait associations than traditional QTL mapping. The objective of TAGGING was to improve prediction power compared with QTL mapping while also providing more specific insights into genetic architecture than genome-wide prediction models. TAGGING was compared with standard QTL mapping using cross validation of empirical data from the maize (Zea mays L.) nested association mapping population. TAGGING-assisted QTL mapping substantially improved prediction ability for both biparental and multifamily populations by reducing both the variance and bias in prediction. Furthermore, an ensemble model combining predictions from TAGGING-assisted QTL and infinitesimal models improved prediction abilities over the component models, indicating some complementarity between model assumptions and suggesting that some trait genetic architectures involve a mixture of a few major QTL and polygenic effects.

Identification of a novel gene for diabetic traits in rats, mice, and humans

The genetic basis of type 2 diabetes remains incompletely defined despite the use of multiple genetic strategies. Multiparental populations such as heterogeneous stocks (HS) facilitate gene discovery by allowing fine mapping to only a few megabases, significantly decreasing the number of potential candidate genes compared to traditional mapping strategies. In the present work, we employed expression and sequence analysis in HS rats (Rattus norvegicus) to identify Tpcn2 as a likely causal gene underlying a 3.1-Mb locus for glucose and insulin levels. Global gene expression analysis on liver identified Tpcn2 as the only gene in the region that is differentially expressed between HS rats with glucose intolerance and those with normal glucose regulation. Tpcn2 also maps as a cis-regulating expression QTL and is negatively correlated with fasting glucose levels. We used founder sequence to identify variants within this region and assessed association between 18 variants and diabetic traits by conducting a mixed-model analysis, accounting for the complex family structure of the HS. We found that two variants were significantly associated with fasting glucose levels, including a nonsynonymous coding variant within Tpcn2. Studies in Tpcn2 knockout mice demonstrated a significant decrease in fasting glucose levels and insulin response to a glucose challenge relative to those in wild-type mice. Finally, we identified variants within Tpcn2 that are associated with fasting insulin in humans. These studies indicate that Tpcn2 is a likely causal gene that may play a role in human diabetes and demonstrate the utility of multiparental populations for positionally cloning genes within complex loci.