An Arabidopsis example of association mapping in structured samples

Association mapping in Salix viminalis L. (Salicaceae) - identification of candidate genes associated with growth and phenology

Willow species (Salix) are important as short-rotation biomass crops for bioenergy, which creates a demand for faster genetic improvement and breeding through deployment of molecular marker-assisted selection (MAS). To find markers associated with important adaptive traits, such as growth and phenology, for use in MAS, we genetically dissected the trait variation of a Salix viminalis (L.) population of 323 accessions. The accessions were sampled throughout northern Europe and were established at two field sites in Pustnäs, Sweden, and at Woburn, UK, offering the opportunity to assess the impact of genotype-by-environment interactions (G × E) on trait-marker associations. Field measurements were recorded for growth and phenology traits. The accessions were genotyped using 1536 SNP markers developed from phenology candidate genes and from genes previously observed to be differentially expressed in contrasting environments. Association mapping between 1233 of these SNPs and the measured traits was performed taking into account population structure and threshold selection bias. At a false discovery rate (FDR) of 0.2, 29 SNPs were associated with bud burst, leaf senescence, number of shoots or shoot diameter. The percentage of accession variation (Radj2) explained by these associations ranged from 0.3% to 4.4%, suggesting that the studied traits are controlled by many loci of limited individual impact. Despite this, a SNP in the EARLY FLOWERING 3 gene was repeatedly associated (FDR < 0.2) with bud burst. The rare homozygous genotype exhibited 0.4-1.0 lower bud burst scores than the other genotype classes on a five-grade scale. Consequently, this marker could be promising for use in MAS and the gene deserves further study. Otherwise, associations were less consistent across sites, likely due to their small Radj2 estimates and to considerable G × E interactions indicated by multivariate association analyses and modest trait accession correlations across sites (0.32-0.61).

VT, Devi SJSR, Mangrauthia SK, Voleti SR, Sarla N. Field level evaluation of rice introgression lines for heat tolerance and validation of markers linked to spikelet fertility

Rice lines derived from wild species and mutants can serve as a good resource for favorable alleles for heat tolerance. In all, 48 stable lines including 17 KMR3/O. rufipogon introgression lines (KMR3 ILs), 15 Swarna/O. nivara ILs (Swarna ILs) along with their parents, Nagina 22 (N22) and its 4 EMS induced mutants and 7 varieties were evaluated for heat tolerance under irrigated conditions under field in two seasons, wet season 2012 using poly cover house method and dry season 2013 using late sown method. Spikelet fertility (SF), yield per plant (YP) and heat susceptibility index (HSI) for these two traits were considered as criteria to assess heat tolerance compared to control. Four KMR3 ILs and eight Swarna ILs were identified as heat tolerant based on SF and YP and their HSIs in both wet and dry seasons. S-65 and S-70 showed low SF and high YP consistently in response to heat in both seasons. We provide evidence that SF alone may not be the best criterion to assess heat tolerance and including YP is important as lines with low SF but high YP and vice versa were identified under heat stress. Out of 49 SSR markers linked to spikelet fertility, 18 were validated for five traits. RM229 in wet season and RM430 and RM210 in dry season were significantly associated with both SF and its HSI under heat stress. RM430 was also significantly associated with both YP and its HSI in dry season. Thirty two candidate genes were identified close to nine markers associated with traits under heat stress.

Association analysis for disease resistance to Fusarium oxysporum in cape gooseberry (Physalis peruviana L)

BACKGROUND

Vascular wilt caused by Fusarium oxysporum is the most important disease in cape gooseberry (Physalis peruviana L.) in Colombia. The development of resistant cultivars is considered one of the most cost-effective means to reduce the impact of this disease. In order to do so, it is necessary to provide breeders with molecular markers and promising germplasm for introgression of different resistance loci as part of breeding schemes. Here we described an association mapping study in cape gooseberry with the goal to: (i) select promising materials for use in plant breeding and (ii) identify SNPs associated with the cape gooseberry resistance response to the F. oxysporum pathogen under greenhouse conditions, as potential markers for cape gooseberry breeding.

RESULTS

We found a total of 21 accessions with different resistance responses within a diversity panel of 100 cape gooseberry accessions. A total of 60,663 SNPs were also identified within the same panel by means of GBS (Genotyping By Sequencing). Model-based population structure and neighbor-joining analyses showed three populations comprising the cape gooseberry panel. After correction for population structure and kinship, we identified SNPs markers associated with the resistance response against F. oxysporum. The identification of markers was based on common tags using the reference genomes of tomato and potato as well as the root/stem transcriptome of cape gooseberry. By comparing their location with the tomato genome, 16 SNPs were found in genes involved in defense/resistance response to pathogens, likewise when compared with the genome of potato, 12 markers were related.

CONCLUSIONS

The work presented herein provides the first association mapping study in cape gooseberry showing both the identification of promising accessions with resistance response phenotypes and the identification of a set of SNP markers mapped to defense/resistance response genes of reference genomes. Thus, the work also provides new knowledge on candidate genes involved in the P. peruviana - F. oxysporum pathosystem as a foundation for further validation in marker-assisted selection. The results have important implications for conservation and breeding strategies in cape gooseberry.

Accounting for Population Structure in Gene-by-Environment Interactions in Genome-Wide Association Studies Using Mixed Models

Although genome-wide association studies (GWASs) have discovered numerous novel genetic variants associated with many complex traits and diseases, those genetic variants typically explain only a small fraction of phenotypic variance. Factors that account for phenotypic variance include environmental factors and gene-by-environment interactions (GEIs). Recently, several studies have conducted genome-wide gene-by-environment association analyses and demonstrated important roles of GEIs in complex traits. One of the main challenges in these association studies is to control effects of population structure that may cause spurious associations. Many studies have analyzed how population structure influences statistics of genetic variants and developed several statistical approaches to correct for population structure. However, the impact of population structure on GEI statistics in GWASs has not been extensively studied and nor have there been methods designed to correct for population structure on GEI statistics. In this paper, we show both analytically and empirically that population structure may cause spurious GEIs and use both simulation and two GWAS datasets to support our finding. We propose a statistical approach based on mixed models to account for population structure on GEI statistics. We find that our approach effectively controls population structure on statistics for GEIs as well as for genetic variants.

Genome-Wide Association Mapping of Fusarium Head Blight Resistance in Wheat using Genotyping-by-Sequencing

Fusarium head blight (FHB) is one of the most important wheat ( L.) diseases worldwide, and host resistance displays complex genetic control. A genome-wide association study (GWAS) was performed on 273 winter wheat breeding lines from the midwestern and eastern regions of the United States to identify chromosomal regions associated with FHB resistance. Genotyping-by-sequencing (GBS) was used to identify 19,992 single-nucleotide polymorphisms (SNPs) covering all 21 wheat chromosomes. Marker-trait associations were performed with different statistical models, the most appropriate being a compressed mixed linear model (cMLM) controlling for relatedness and population structure. Ten significant SNP-trait associations were detected on chromosomes 4A, 6A, 7A, 1D, 4D, and 7D, and multiple SNPs were associated with on chromosome 3B. Although combination of favorable alleles of these SNPs resulted in lower levels of severity (SEV), incidence (INC), and deoxynivalenol concentration (DON), lines carrying multiple beneficial alleles were in very low frequency for most traits. These SNPs can now be used for creating new breeding lines with different combinations of favorable alleles. This is one of the first GWAS using genomic resources from the International Wheat Genome Sequencing Consortium (IWGSC).

50 years of Arabidopsis research: highlights and future directions

922 I. 922 II. 922 III. 925 IV. 925 V. 926 VI. 927 VII. 928 VIII. 929 IX. 930 X. 931 XI. 932 XII. 933 XIII. Natural variation and genome-wide association studies 934 XIV. 934 XV. 935 XVI. 936 XVII. 937 937 References 937 SUMMARY: The year 2014 marked the 25(th) International Conference on Arabidopsis Research. In the 50 yr since the first International Conference on Arabidopsis Research, held in 1965 in Göttingen, Germany, > 54 000 papers that mention Arabidopsis thaliana in the title, abstract or keywords have been published. We present herein a citational network analysis of these papers, and touch on some of the important discoveries in plant biology that have been made in this powerful model system, and highlight how these discoveries have then had an impact in crop species. We also look to the future, highlighting some outstanding questions that can be readily addressed in Arabidopsis. Topics that are discussed include Arabidopsis reverse genetic resources, stock centers, databases and online tools, cell biology, development, hormones, plant immunity, signaling in response to abiotic stress, transporters, biosynthesis of cells walls and macromolecules such as starch and lipids, epigenetics and epigenomics, genome-wide association studies and natural variation, gene regulatory networks, modeling and systems biology, and synthetic biology.

Iterative Usage of Fixed and Random Effect Models for Powerful and Efficient Genome-Wide Association Studies

False positives in a Genome-Wide Association Study (GWAS) can be effectively controlled by a fixed effect and random effect Mixed Linear Model (MLM) that incorporates population structure and kinship among individuals to adjust association tests on markers; however, the adjustment also compromises true positives. The modified MLM method, Multiple Loci Linear Mixed Model (MLMM), incorporates multiple markers simultaneously as covariates in a stepwise MLM to partially remove the confounding between testing markers and kinship. To completely eliminate the confounding, we divided MLMM into two parts: Fixed Effect Model (FEM) and a Random Effect Model (REM) and use them iteratively. FEM contains testing markers, one at a time, and multiple associated markers as covariates to control false positives. To avoid model over-fitting problem in FEM, the associated markers are estimated in REM by using them to define kinship. The P values of testing markers and the associated markers are unified at each iteration. We named the new method as Fixed and random model Circulating Probability Unification (FarmCPU). Both real and simulated data analyses demonstrated that FarmCPU improves statistical power compared to current methods. Additional benefits include an efficient computing time that is linear to both number of individuals and number of markers. Now, a dataset with half million individuals and half million markers can be analyzed within three days.

Identification of new heading date determinants in wheat 5B chromosome

BACKGROUND

Variability of heading date may assist in wheat adaptation to local environments. Thereafter, discovery of new heading date determinants is important for cereal improvement. In this study we used common wheat cultivar Chinese Spring (CS) and the substitution line of CS with 5B chromosome from T. dicoccoides (CS-5Bdic), different in their heading date by two weeks, to detect determinants of heading date on 5B chromosome.

RESULTS

The possible influence of the VRN-B1 gene, the most powerful regulator of flowering, located on 5B chromosome, to differences in heading time between CS and CS-5Bdic was studied. The sequencing of this gene from CS-5Bdic showed that an insertion of a nucleotide triplet produced an additional amino acid in the corresponding protein. No changes in the transcription levels of each homoeologous VRN-1 loci were found in CS-5Bdic by comparison with CS. To ascertain the loci determining heading date difference, a set of 116 recombinant inbred 5В chromosomal lines as a result of hybridization of CS with CS-5Bdic were developed and their heading dates were estimated. Using the Illumina Infinium 15 k Wheat platform, 379 5B-specific polymorphic markers were detected and a genetic map with 82 skeletal markers was constructed. Phenotype (heading date) - genotype association analysis revealed seventy eight markers in pericentromeric region of 5B chromosome significantly associated with heading date variation. Based on this estimation and synteny with model crop genomes we identified the three best candidate genes: WRKY, ERF/AP2 and FHY3/FAR1.

CONCLUSIONS

We supposed that the difference in activity of WRKY, ERF/AP2 and/or FHY3/FAR1 transcription factors between CS and CS-5Bdic to be a probable reason for the observed difference in heading dates. Data obtained in this study provide a good basis for the subsequent investigation of heading time pathways in wheat.

Characterizing Variation of Branch Angle and Genome-Wide Association Mapping in Rapeseed (Brassica napus L.)

Changes in the rapeseed branch angle alter plant architecture, allowing more efficient light capture as planting density increases. In this study, a natural population of rapeseed was grown in three environments and evaluated for branch angle trait to characterize their phenotypic patterns and genotype with a 60K Brassica Infinium SNP array. Significant phenotypic variation was observed from 20 to 70°. As a result, 25 significant quantitative trait loci (QTL) associated with branch angle were identified on chromosomes A2, A3, A7, C3, C5, and C7 by the MLM model in TASSEL 4.0. Orthologs of the functional candidate genes involved in branch angle were identified. Among the key QTL, the peak SNPs were close to the key orthologous genes BnaA.Lazy1 and BnaC.Lazy1 on A3 and C3 homologous genome blocks. With the exception of Lazy (LA) orthologous genes, SQUMOSA PROMOTER BINDING PROTEIN LIKE 14 (SPL14) and an auxin-responsive GRETCHEN HAGEN 3 (GH3) genes from Arabidopsis thaliana were identified close to two clusters of SNPs on the A7 and C7 chromosomes. These findings on multiple novel loci and candidate genes of branch angle will be useful for further understanding and genetic improvement of plant architecture in rapeseed.

Detection of Favorable QTL Alleles and Candidate Genes for Lint Percentage by GWAS in Chinese Upland Cotton

Improving cotton yield is a major breeding goal for Chinese upland cotton. Lint percentage is an important yield component and a critical economic index for cotton cultivars, and raising the lint percentage has a close relationship to improving cotton lint yield. To investigate the genetic architecture of lint percentage, a diversity panel consisting of 355 upland cotton accessions was grown, and the lint percentage was measured in four different environments. Genotyping was performed with specific-locus amplified fragment sequencing (SLAF-seq). Twelve single-nucleotide polymorphisms (SNPs) associated with lint percentage were detected via a genome-wide association study (GWAS), in which five SNP loci distributed on chromosomes A_t3 (A02) and A_t4 (A08) and contained two major-effect QTLs, which were detected in the best linear unbiased predictions (BLUPs) and in more than three environments simultaneously. Furthermore, favorable haplotypes (FHs) of two major-effect QTLs and 47 putative candidate genes in the two linkage disequilibrium (LD) blocks of these associated loci were identified. The expression levels of these putative candidate genes were estimated using RNA-seq data from ten upland cotton tissues. We found that Gh_A02G1268 was very highly expressed during the early fiber development stage, whereas the gene was poorly expressed in the seed. These results implied that Gh_A02G1268 may determine the lint percentage by regulating seed and fiber development. The favorable QTL alleles and candidate genes for lint percentage identified in this study will have high potential for improving lint yield in future Chinese cotton breeding programs.

Genome Wide Association Mapping in Arabidopsis thaliana Identifies Novel Genes Involved in Linking Allyl Glucosinolate to Altered Biomass and Defense

A key limitation in modern biology is the ability to rapidly identify genes underlying newly identified complex phenotypes. Genome wide association studies (GWAS) have become an increasingly important approach for dissecting natural variation by associating phenotypes with genotypes at a genome wide level. Recent work is showing that the Arabidopsis thaliana defense metabolite, allyl glucosinolate (GSL), may provide direct feedback regulation, linking defense metabolism outputs to the growth, and defense responses of the plant. However, there is still a need to identify genes that underlie this process. To start developing a deeper understanding of the mechanism(s) that modulate the ability of exogenous allyl GSL to alter growth and defense, we measured changes in plant biomass and defense metabolites in a collection of natural 96 A. thaliana accessions fed with 50 μM of allyl GSL. Exogenous allyl GSL was introduced exclusively to the roots and the compound transported to the leaf leading to a wide range of heritable effects upon plant biomass and endogenous GSL accumulation. Using natural variation we conducted GWAS to identify a number of new genes which potentially control allyl responses in various plant processes. This is one of the first instances in which this approach has been successfully utilized to begin dissecting a novel phenotype to the underlying molecular/polygenic basis.

Genetic Mapping of Anaerobic Germination-Associated QTLs Controlling Coleoptile Elongation in Rice

BACKGROUND

Increasing numbers of rice farmers are adopting direct-seeding methods to save on costs associated with labor and transplanting. Successful seedling establishment in flooded conditions requires rapid coleoptile growth to ensure access to oxygen near the water surface. It is important that the natural variations in coleoptile growth of submerged rice plants are identified.

RESULTS

Coleoptile responses of submerged plants at the germination stage were analyzed in diverse rice accessions and recombinant inbred lines. Several genomic regions identified from a genome-wide association analysis were significantly associated with anaerobic germination, with many that corresponded to published quantitative trait locus (QTL) intervals. In the recombinant inbred line population derived from a cross between japonica and indica varieties, only one unique and strong signal explaining about 27 % of the phenotypic variation was detected. Distinct haplotypes associated with variations in coleoptile length were identified in diverse germplasm.

CONCLUSIONS

We demonstrated the value of combining genome-wide association analysis and biparental QTL mapping approaches to identify chromosomal regions regulating coleoptile elongation in submerged rice plants. The significant genomic regions detected in this study are potential candidates for incorporation into elite cultivars to improve seedling survival during anaerobic germination. Future studies that map the QTLs and investigate the effects and functions of candidate genes may lead to new rice varieties that can be used in direct-seeding systems.

Genome-wide association study of phosphorus-deficiency-tolerance traits in Aegilops tauschii

Using GWAS, 13 significant SNPs distributed on six of the seven Aegilops tauschii chromosomes (all but 5D) were identified, and several candidate P-deficiency-responsive genes were proposed from searches of public databases. Aegilops tauschii, the wheat (Triticum aestivum) D-genome progenitor, possesses numerous genes for stress resistance, including genes for tolerance of phosphorus (P) deficiency. Investigation of the genetic architecture of A. tauschii will help in developing P-deficiency-tolerant varieties of wheat. We evaluated nine traits in a population of 380 A. tauschii specimens under conditions with and without P application, and we performed genome-wide association studies for these traits using single nucleotide polymorphism (SNP) chips containing 7185 markers. Using a general linear model, we identified 119 SNPs that were significantly associated with all nine traits, and a mixed linear model revealed 18 SNPs associated with all traits. Both models detected 13 significant markers distributed on six of the seven A. tauschii chromosomes (all but 5D). Searches of public databases revealed several candidate/flanking genes related to P-deficiency tolerance. These genes were grouped in five categories by the types of proteins they encoded: defense response proteins, enzymes, promoters and transcription factors, storage proteins, or proteins triggered by P deficiency. The identified SNPs and genes contain essential information for cloning genes related to P-deficiency tolerance in A. tauschii and wheat, and they provide a foundation for breeding P-deficiency tolerant wheat cultivars.

Population Genomics for Understanding Adaptation in Wild Plant Species

Darwin's theory of evolution by natural selection is the foundation of modern biology. However, it has proven remarkably difficult to demonstrate at the genetic, genomic, and population level exactly how wild species adapt to their natural environments. We discuss how one can use large sets of multiple genome sequences from wild populations to understand adaptation, with an emphasis on the small herbaceous plant Arabidopsis thaliana. We present motivation for such studies; summarize progress in describing whole-genome, species-wide sequence variation; and then discuss what insights have emerged from these resources, either based on sequence information alone or in combination with phenotypic data. We conclude with thoughts on opportunities with other plant species and the impact of expected progress in sequencing technology and genome engineering for studying adaptation in nature.

"Missing" G x E Variation Controls Flowering Time in Arabidopsis thaliana

Understanding how genetic variation interacts with the environment is essential for understanding adaptation. In particular, the life cycle of plants is tightly coordinated with local environmental signals through complex interactions with the genetic variation (G x E). The mechanistic basis for G x E is almost completely unknown. We collected flowering time data for 173 natural inbred lines of Arabidopsis thaliana from Sweden under two growth temperatures (10°C and 16°C), and observed massive G x E variation. To identify the genetic polymorphisms underlying this variation, we conducted genome-wide scans using both SNPs and local variance components. The SNP-based scan identified several variants that had common effects in both environments, but found no trace of G x E effects, whereas the scan using local variance components found both. Furthermore, the G x E effects appears to be concentrated in a small fraction of the genome (0.5%). Our conclusion is that G x E effects in this study are mostly due to large numbers of allele or haplotypes at a small number of loci, many of which correspond to previously identified flowering time genes.

Integrating Image-Based Phenomics and Association Analysis to Dissect the Genetic Architecture of Temporal Salinity Responses in Rice

Salinity affects a significant portion of arable land and is particularly detrimental for irrigated agriculture, which provides one-third of the global food supply. Rice (Oryza sativa), the most important food crop, is salt sensitive. The genetic resources for salt tolerance in rice germplasm exist but are underutilized due to the difficulty in capturing the dynamic nature of physiological responses to salt stress. The genetic basis of these physiological responses is predicted to be polygenic. In an effort to address this challenge, we generated temporal imaging data from 378 diverse rice genotypes across 14 d of 90 mm NaCl stress and developed a statistical model to assess the genetic architecture of dynamic salinity-induced growth responses in rice germplasm. A genomic region on chromosome 3 was strongly associated with the early growth response and was captured using visible range imaging. Fluorescence imaging identified four genomic regions linked to salinity-induced fluorescence responses. A region on chromosome 1 regulates both the fluorescence shift indicative of the longer term ionic stress and the early growth rate decline during salinity stress. We present, to our knowledge, a new approach to capture the dynamic plant responses to its environment and elucidate the genetic basis of these responses using a longitudinal genome-wide association model.

Integrating Image-Based Phenomics and Association Analysis to Dissect the Genetic Architecture of Temporal Salinity Responses in Rice

Salinity affects a significant portion of arable land and is particularly detrimental for irrigated agriculture, which provides one-third of the global food supply. Rice (Oryza sativa), the most important food crop, is salt sensitive. The genetic resources for salt tolerance in rice germplasm exist but are underutilized due to the difficulty in capturing the dynamic nature of physiological responses to salt stress. The genetic basis of these physiological responses is predicted to be polygenic. In an effort to address this challenge, we generated temporal imaging data from 378 diverse rice genotypes across 14 d of 90 mm NaCl stress and developed a statistical model to assess the genetic architecture of dynamic salinity-induced growth responses in rice germplasm. A genomic region on chromosome 3 was strongly associated with the early growth response and was captured using visible range imaging. Fluorescence imaging identified four genomic regions linked to salinity-induced fluorescence responses. A region on chromosome 1 regulates both the fluorescence shift indicative of the longer term ionic stress and the early growth rate decline during salinity stress. We present, to our knowledge, a new approach to capture the dynamic plant responses to its environment and elucidate the genetic basis of these responses using a longitudinal genome-wide association model.

Genetic diversity and association mapping of iron and zinc concentrations in chickpea (Cicer arietinum L.)

Chickpea (Cicer arietinum L.) is the world's second most important pulse crop after common bean. Chickpea has historically been an important daily staple in the diet of millions of people, especially in the developing countries. Current chickpea breeding programs have mainly been directed toward high yield, biotic and abiotic stress resilience that has increased global production, but less attention has been directed toward improving micronutrient concentrations in seeds. In an effort to develop micronutrient-dense chickpea lines, a study to examine the variability and to identify SNP alleles associated with seed iron and zinc concentrations was conducted using 94 diverse accessions of chickpea. The results indicated that there is substantial variability present in chickpea germplasm for seed iron and zinc concentrations. In the current set of germplasm, zinc is negatively correlated with grain yield across all locations and years; whereas the negative correlation between iron and grain yield was only significant at the Elrose locality. Eight SNP loci associated with iron and (or) zinc concentrations in chickpea seeds were identified. One SNP located on chromosome 1 (chr1) is associated with both iron and zinc concentrations. On chr4, three SNPs associated with zinc concentration and two SNPs for iron concentration were identified. Two additional SNP loci, one on chr6 and the other on chr7, were also found to be associated with iron and zinc concentrations, respectively. The results show potential opportunity for molecular breeding for improvement of seed iron and zinc concentrations in chickpea.

Genome-Wide Association Mapping for Leaf Tip Necrosis and Pseudo-black Chaff in Relation to Durable Rust Resistance in Wheat

The partial rust resistance genes Lr34 and Sr2 have been used extensively in wheat (Triticum aestivum L.) improvement, as they confer exceptional durability. Interestingly, the resistance of Lr34 is associated with the expression of leaf tip necrosis (LTN) and Sr2 with pseudo-black chaff (PBC). Genome-wide association mapping using CIMMYT's stem rust resistance screening nursery (SRRSN) wheat lines was done to identify genotyping-by-sequencing (GBS) markers linked to LTN and PBC. Phenotyping for these traits was done in Ithaca, New York (fall 2011); Njoro, Kenya (main and off-seasons, 2012), and Wellington, India (winter, 2013). Using the mixed linear model (MLM), 18 GBS markers were significantly associated with LTN. While some markers were linked to loci where the durable leaf rust resistance genes Lr34 (7DS), Lr46 (1BL), and Lr68 (7BL) were mapped, significant associations were also detected with other loci on 2BL, 5B, 3BS, 4BS, and 7BS. Twelve GBS markers linked to the Sr2 locus (3BS) and loci on 2DS, 4AL, and 7DS were significantly associated with PBC. This study provides insight into the complex genetic control of LTN and PBC. Further efforts to validate and study these loci might aid in determining the nature of their association with durable resistance.

Genome-wide DArT and SNP scan for QTL associated with resistance to stripe rust (Puccinia striiformis f. sp. tritici) in elite ICARDA wheat (Triticum aestivum L.) germplasm

Identified DArT and SNP markers including a first reported QTL on 3AS, validated large effect APR on 3BS. The different genes can be used to incorporate stripe resistance in cultivated varieties. Stripe rust [yellow rust, caused by Puccinia striiformis f. sp. tritici (Pst)] is a serious disease in wheat (Triticum aestivum). This study employed genome-wide association mapping (GWAM) to identify markers linked to stripe rust resistance genes using Diversity Arrays Technology (DArT(®)) and single-nucleotide polymorphism (SNP) Infinium 9K assays in 200 ICARDA wheat genotypes, phenotyped for seedling and adult plant resistance in two sites over two growing seasons in Syria. Only 25.8 % of the genotypes showed resistance at seedling stage while about 33 and 44 % showed moderate resistance and resistance response, respectively. Mixed-linear model adjusted for false discovery rate at p < 0.05 identified 12 DArT and 29 SNP markers on chromosome arms 3AS, 3AL, 1AL, 2AL, 2BS, 2BL, 3BS, 3BL, 5BL, 6AL, and 7DS significantly linked to Pst resistance genes. Of these, the locus on 3AS has not been previously reported to confer resistance to stripe rust in wheat. The QTL on 3AS, 3AL, 1AL, 2AL, and 2BS were effective at seedling and adult plant growth stages while those on 3BS, 3BL, 5BL, 6AL and 7DS were effective at adult plant stage. The 3BS QTL was validated in Cham-6 × Cham-8 recombinant inbred line population; composite interval analysis identified a stripe resistance QTL flanked by the DArT marker, wPt-798970, contributed by Cham-6 parent which accounted for 31.2 % of the phenotypic variation. The DArT marker "wPt-798970" lies 1.6 cM away from the 3BS QTL detected within GWAM. Epistatic interactions were also investigated; only the QTL on 1AL, 3AS and 6AL exhibited interactions with other loci. These results suggest that GWAM can be an effective approach for identifying and improving resistance to stripe rust in wheat.

Genome-Wide Association Mapping for Resistance to Leaf and Stripe Rust in Winter-Habit Hexaploid Wheat Landraces

Leaf rust, caused by Puccinia triticina (Pt), and stripe rust, caused by P. striiformis f. sp. tritici (Pst), are destructive foliar diseases of wheat worldwide. Breeding for disease resistance is the preferred strategy of managing both diseases. The continued emergence of new races of Pt and Pst requires a constant search for new sources of resistance. Here we report a genome-wide association analysis of 567 winter wheat (Triticum aestivum) landrace accessions using the Infinium iSelect 9K wheat SNP array to identify loci associated with seedling resistance to five races of Pt (MDCL, MFPS, THBL, TDBG, and TBDJ) and one race of Pst (PSTv-37) frequently found in the Northern Great Plains of the United States. Mixed linear models identified 65 and eight significant markers associated with leaf rust and stripe rust, respectively. Further, we identified 31 and three QTL associated with resistance to Pt and Pst, respectively. Eleven QTL, identified on chromosomes 3A, 4A, 5A, and 6D, are previously unknown for leaf rust resistance in T. aestivum.

Genetic Architecture of Natural Variation in Rice Chlorophyll Content Revealed by a Genome-Wide Association Study

Chlorophyll content is one of the most important physiological traits as it is closely related to leaf photosynthesis and crop yield potential. So far, few genes have been reported to be involved in natural variation of chlorophyll content in rice (Oryza sativa) and the extent of variations explored is very limited. We conducted a genome-wide association study (GWAS) using a diverse worldwide collection of 529 O. sativa accessions. A total of 46 significant association loci were identified. Three F2 mapping populations with parents selected from the association panel were tested for validation of GWAS signals. We clearly demonstrated that Grain number, plant height, and heading date7 (Ghd7) was a major locus for natural variation of chlorophyll content at the heading stage by combining evidence from near-isogenic lines and transgenic plants. The enhanced expression of Ghd7 decreased the chlorophyll content, mainly through down-regulating the expression of genes involved in the biosynthesis of chlorophyll and chloroplast. In addition, Narrow leaf1 (NAL1) corresponded to one significant association region repeatedly detected over two years. We revealed a high degree of polymorphism in the 5' UTR and four non-synonymous SNPs in the coding region of NAL1, and observed diverse effects of the major haplotypes. The loci or candidate genes identified would help to fine-tune and optimize the antenna size of canopies in rice breeding.

Genome-wide association mapping in plants

We present new association mapping methods which address the unique challenges of analyzing genome-wide data from multi-environment plant studies. Association studies on a genome-wide scale are being performed in plants. Unlike human studies, plant studies contain replicates whose data may be recorded across different environments. Plant studies also often employ elaborate experimental designs for controlling extraneous phenotypic variation. As a result, the genome-wide analysis of data from plant studies can be challenging. In this paper, we present QK-based association mapping for the analysis of data from plant association studies. In doing so, we have developed: (a) a general multivariate QK framework for association mapping in plant studies of arbitrary complexity; (b) a new weighted two-stage analysis approach for QK-based association mapping; (c) a heuristic procedure for determining when two-stage analysis is appropriate; and (d) a Monte Carlo sampling procedure for controlling the genome-wide type I error rate. We conduct a simulation study to evaluate the performance of our genome-wide mapping technique. We also analyze data from a multi-environment association study in wheat.

Mixed model with correction for case-control ascertainment increases association power

We introduce a liability-threshold mixed linear model (LTMLM) association statistic for case-control studies and show that it has a well-controlled false-positive rate and more power than existing mixed-model methods for diseases with low prevalence. Existing mixed-model methods suffer a loss in power under case-control ascertainment, but no solution has been proposed. Here, we solve this problem by using a χ(2) score statistic computed from posterior mean liabilities (PMLs) under the liability-threshold model. Each individual's PML is conditional not only on that individual's case-control status but also on every individual's case-control status and the genetic relationship matrix (GRM) obtained from the data. The PMLs are estimated with a multivariate Gibbs sampler; the liability-scale phenotypic covariance matrix is based on the GRM, and a heritability parameter is estimated via Haseman-Elston regression on case-control phenotypes and then transformed to the liability scale. In simulations of unrelated individuals, the LTMLM statistic was correctly calibrated and achieved higher power than existing mixed-model methods for diseases with low prevalence, and the magnitude of the improvement depended on sample size and severity of case-control ascertainment. In a Wellcome Trust Case Control Consortium 2 multiple sclerosis dataset with >10,000 samples, LTMLM was correctly calibrated and attained a 4.3% improvement (p = 0.005) in χ(2) statistics over existing mixed-model methods at 75 known associated SNPs, consistent with simulations. Larger increases in power are expected at larger sample sizes. In conclusion, case-control studies of diseases with low prevalence can achieve power higher than that in existing mixed-model methods.

A genome-wide association study of malting quality across eight U.S. barley breeding programs

We report malt quality QTLs relevant to breeding with greater precision than previous mapping studies. The distribution of favorable alleles suggests strategies for marker-assisted breeding and germplasm exchange. This study leverages the breeding data of 1,862 barley breeding lines evaluated in 97 field trials for genome-wide association study of malting quality traits in barley. The mapping panel consisted of six-row and two-row advanced breeding lines from eight breeding populations established at six public breeding programs across the United States. A total of 4,976 grain samples were subjected to micro-malting analysis and mapping of nine quality traits was conducted with 3,072 SNP markers distributed throughout the genome. Association mapping was performed for individual breeding populations and for combined six-row and two-row populations. Only 16% of the QTL we report here had been detected in prior bi-parental mapping studies. Comparison of the analyses of the combined two-row and six-row panels identified only two QTL regions that were common to both. In total, 108 and 107 significant marker-trait associations were identified in all six-row and all two-row breeding programs, respectively. A total of 102 and 65 marker-trait associations were specific to individual six-row and two-row breeding programs, respectively indicating that most marker-trait associations were breeding population specific. Combining datasets from different breeding program resulted in both the loss of some QTL that were apparent in the analyses of individual programs and the discovery of new QTL not identified in individual programs. This suggests that simply increasing sample size by pooling samples with different breeding history does not necessarily increase the power to detect associations. The genetic architecture of malting quality and the distribution of favorable alleles suggest strategies for marker-assisted selection and germplasm exchange.

Association mapping of seedling resistance to spot form net blotch in a worldwide collection of barley

Spot form net blotch (SFNB), caused by the necrotrophic fungal pathogen Pyrenophora teres f. maculata, is an important foliar disease of barley in major production regions around the world. Deployment of adequate host resistance is challenging because the virulence of P. teres f. maculata is highly variable and characterized minor-effect resistances are typically ineffective against the diverse pathogen populations. A world barley core collection consisting of 2,062 barley accessions of diverse origin and genotype were phenotyped at the seedling stage with four P. teres f. maculata isolates collected from the United States (FGO), New Zealand (NZKF2), Australia (SG1), and Denmark (DEN 2.6). Of the 2,062 barley accessions phenotyped, 1,480 were genotyped with the Illumina barley iSelect chip and passed the quality controls with 5,954 polymorphic markers used for further association mapping analysis. Genome-wide association mapping was utilized to identify and map resistance loci from the seedling disease response data and the single nucleotide polymorphism (SNP) marker data. The best among six different regression models was identified for each isolate and association analysis was performed separately for each. A total of 138 significant (-log10P value>3.0) marker-trait associations (MTA) were detected. Using a 5 cM cutoff, a total of 10, 8, 13, and 10 quantitative trait loci (QTL) associated with SFNB resistance were identified for the FGO, SG1, NZKF2, and DEN 2.6 isolates, respectively. Loci containing from 1 to 34 MTA were identified on all seven barley chromosomes with one locus at 66 to 69 cM on chromosome 2H common to all four isolates. Six distinct loci were identified by the association mapping (AM) analysis that corresponded to previously characterized SFNB resistance QTL identified by biparental population analysis (QRpt4, QRpt6, Rpt4, Rpt6, Rpt7, and a QTL on 4H that was not given a provisional gene or QTL nomenclature). The 21 putative novel loci identified may represent a broad spectrum of resistance and or susceptibility loci. This is the first comprehensive AM study to characterize SFNB resistance loci underlying broad populations of the barley host and P. teres f. maculata pathogen.

Genome-wide association study for flowering time, maturity dates and plant height in early maturing soybean (Glycine max) germplasm

BACKGROUND

Soybean (Glycine max) is a photoperiod-sensitive and self-pollinated species. Days to flowering (DTF) and maturity (DTM), duration of flowering-to-maturity (DFTM) and plant height (PH) are crucial for soybean adaptability and yield. To dissect the genetic architecture of these agronomically important traits, a population consisting of 309 early maturity soybean germplasm accessions was genotyped with the Illumina Infinium SoySNP50K BeadChip and phenotyped in multiple environments. A genome-wide association study (GWAS) was conducted using a mixed linear model that involves both relative kinship and population structure.

RESULTS

The linkage disequilibrium (LD) decayed slowly in soybean, and a substantial difference in LD pattern was observed between euchromatic and heterochromatic regions. A total of 27, 6, 18 and 27 loci for DTF, DTM, DFTM and PH were detected via GWAS, respectively. The Dt1 gene was identified in the locus strongly associated with both DTM and PH. Ten candidate genes homologous to Arabidopsis flowering genes were identified near the peak single nucleotide polymorphisms (SNPs) associated with DTF. Four of them encode MADS-domain containing proteins. Additionally, a pectin lyase-like gene was also identified in a major-effect locus for PH where LD decayed rapidly.

CONCLUSIONS

This study identified multiple new loci and refined chromosomal regions of known loci associated with DTF, DTM, DFTM and/or PH in soybean. It demonstrates that GWAS is powerful in dissecting complex traits and identifying candidate genes although LD decayed slowly in soybean. The loci and trait-associated SNPs identified in this study can be used for soybean genetic improvement, especially the major-effect loci associated with PH could be used to improve soybean yield potential. The candidate genes may serve as promising targets for studies of molecular mechanisms underlying the related traits in soybean.

A genome-wide association study in large white and landrace pig populations for number piglets born alive

The number of piglets born alive (NBA) per litter is one of the most important traits in pig breeding due to its influence on production efficiency. It is difficult to improve NBA because the heritability of the trait is low and it is governed by a high number of loci with low to moderate effects. To clarify the biological and genetic background of NBA, genome-wide association studies (GWAS) were performed using 4,012 Large White and Landrace pigs from herdbook and commercial breeding companies in Germany (3), Austria (1) and Switzerland (1). The animals were genotyped with the Illumina PorcineSNP60 BeadChip. Because of population stratifications within and between breeds, clusters were formed using the genetic distances between the populations. Five clusters for each breed were formed and analysed by GWAS approaches. In total, 17 different significant markers affecting NBA were found in regions with known effects on female reproduction. No overlapping significant chromosome areas or QTL between Large White and Landrace breed were detected.

Genetic diversity and association mapping of bacterial blight and other horticulturally important traits with microsatellite markers in pomegranate from India

This genetic diversity study aimed to estimate the population structure and explore the use of association mapping strategies to identify linked markers for bacterial resistance, growth and fruit quality in pomegranate collections from India. In total, 88 accessions including 37 cultivated types were investigated. A total of 112 alleles were amplified by use of 44 publicly available microsatellites for estimating molecular genetic diversity and population structure. Neighbor-joining analysis, model-based population structure and principal component analysis corroborated the genetic relationships among wild-type and cultivated pomegranate collections from India. Our study placed all 88 germplasm into four clusters. We identified a cultivated clade of pomegranates in close proximity to Daru types of wild-type pomegranates that grow naturally near the foothills of the Himalayas. Admixture analysis sorted various lineages of cultivated pomegranates to their respective ancestral forms. We identified four linked markers for fruit weight, titratable acidity and bacterial blight severity. PGCT001 was found associated with both fruit weight and bacterial blight, and the association with fruit weight during both seasons analyzed was significant after Bonferroni correction. This research demonstrates effectiveness of microsatellites to resolve population structure among the wild and cultivar collection of pomegranates and future use for association mapping studies.

Identification of loci governing eight agronomic traits using a GBS-GWAS approach and validation by QTL mapping in soya bean

Soya bean is a major source of edible oil and protein for human consumption as well as animal feed. Understanding the genetic basis of different traits in soya bean will provide important insights for improving breeding strategies for this crop. A genome-wide association study (GWAS) was conducted to accelerate molecular breeding for the improvement of agronomic traits in soya bean. A genotyping-by-sequencing (GBS) approach was used to provide dense genome-wide marker coverage (>47,000 SNPs) for a panel of 304 short-season soya bean lines. A subset of 139 lines, representative of the diversity among these, was characterized phenotypically for eight traits under six environments (3 sites × 2 years). Marker coverage proved sufficient to ensure highly significant associations between the genes known to control simple traits (flower, hilum and pubescence colour) and flanking SNPs. Between one and eight genomic loci associated with more complex traits (maturity, plant height, seed weight, seed oil and protein) were also identified. Importantly, most of these GWAS loci were located within genomic regions identified by previously reported quantitative trait locus (QTL) for these traits. In some cases, the reported QTLs were also successfully validated by additional QTL mapping in a biparental population. This study demonstrates that integrating GBS and GWAS can be used as a powerful complementary approach to classical biparental mapping for dissecting complex traits in soya bean.

Genome-wide association mapping of salinity tolerance in rice (Oryza sativa)

Salinity tolerance in rice is highly desirable to sustain production in areas rendered saline due to various reasons. It is a complex quantitative trait having different components, which can be dissected effectively by genome-wide association study (GWAS). Here, we implemented GWAS to identify loci controlling salinity tolerance in rice. A custom-designed array based on 6,000 single nucleotide polymorphisms (SNPs) in as many stress-responsive genes, distributed at an average physical interval of <100 kb on 12 rice chromosomes, was used to genotype 220 rice accessions using Infinium high-throughput assay. Genetic association was analysed with 12 different traits recorded on these accessions under field conditions at reproductive stage. We identified 20 SNPs (loci) significantly associated with Na⁺/K⁺ ratio, and 44 SNPs with other traits observed under stress condition. The loci identified for various salinity indices through GWAS explained 5–18% of the phenotypic variance. The region harbouring Saltol, a major quantitative trait loci (QTLs) on chromosome 1 in rice, which is known to control salinity tolerance at seedling stage, was detected as a major association with Na⁺/K⁺ ratio measured at reproductive stage in our study. In addition to Saltol, we also found GWAS peaks representing new QTLs on chromosomes 4, 6 and 7. The current association mapping panel contained mostly indica accessions that can serve as source of novel salt tolerance genes and alleles. The gene-based SNP array used in this study was found cost-effective and efficient in unveiling genomic regions/candidate genes regulating salinity stress tolerance in rice.

Genome-Wide Association Mapping for Tomato Volatiles Positively Contributing to Tomato Flavor

Tomato volatiles, mainly derived from essential nutrients and health-promoting precursors, affect tomato flavor. Taste volatiles present a major challenge for flavor improvement and quality breeding. In this study, we performed genome-wide association studies (GWAS) to investigate potential chromosome regions associated with the tomato flavor volatiles. We observed significant variation (1200x) among the selected 28 most important volatiles in tomato based on their concentration and odor threshold importance across our sampled accessions. Using 174 tomato accessions, GWAS identified 125 significant associations (P < 0.005) among 182 SSR markers and 28 volatiles (27 volatiles with at least one significant association). Several significant associations were co-localized in previously identified quantitative trait loci (QTL). This result provides new potential candidate loci affecting the metabolism of several volatiles.

Principal component regression and linear mixed model in association analysis of structured samples: competitors or complements?

Genome-wide association studies (GWAS) have been established as a major tool to identify genetic variants associated with complex traits, such as common diseases. However, GWAS may suffer from false positives and false negatives due to confounding population structures, including known or unknown relatedness. Another important issue is unmeasured environmental risk factors. Among many methods for adjusting for population structures, two approaches stand out: one is principal component regression (PCR) based on principal component analysis, which is perhaps the most popular due to its early appearance, simplicity, and general effectiveness; the other is based on a linear mixed model (LMM) that has emerged recently as perhaps the most flexible and effective, especially for samples with complex structures as in model organisms. As shown previously, the PCR approach can be regarded as an approximation to an LMM; such an approximation depends on the number of the top principal components (PCs) used, the choice of which is often difficult in practice. Hence, in the presence of population structure, the LMM appears to outperform the PCR method. However, due to the different treatments of fixed vs. random effects in the two approaches, we show an advantage of PCR over LMM: in the presence of an unknown but spatially confined environmental confounder (e.g., environmental pollution or lifestyle), the PCs may be able to implicitly and effectively adjust for the confounder whereas the LMM cannot. Accordingly, to adjust for both population structures and nongenetic confounders, we propose a hybrid method combining the use and, thus, strengths of PCR and LMM. We use real genotype data and simulated phenotypes to confirm the above points, and establish the superior performance of the hybrid method across all scenarios.

Detection of genomic loci associated with environmental variables using generalized linear mixed models

We tested the use of Generalized Linear Mixed Models to detect associations between genetic loci and environmental variables, taking into account the population structure of sampled individuals. We used a simulation approach to generate datasets under demographically and selectively explicit models. These datasets were used to analyze and optimize GLMM capacity to detect the association between markers and selective coefficients as environmental data in terms of false and true positive rates. Different sampling strategies were tested, maximizing the number of populations sampled, sites sampled per population, or individuals sampled per site, and the effect of different selective intensities on the efficiency of the method was determined. Finally, we apply these models to an Arabidopsis thaliana SNP dataset from different accessions, looking for loci associated with spring minimal temperature. We identified 25 regions that exhibit unusual correlations with the climatic variable and contain genes with functions related to temperature stress.

Applying association mapping and genomic selection to the dissection of key traits in elite European wheat

We show the application of association mapping and genomic selection for key breeding targets using a large panel of elite winter wheat varieties and a large volume of agronomic data. The heightening urgency to increase wheat production in line with the needs of a growing population, and in the face of climatic uncertainty, mean new approaches, including association mapping (AM) and genomic selection (GS) need to be validated and applied in wheat breeding. Key adaptive responses are the cornerstone of regional breeding. There is evidence that new ideotypes for long-standing traits such as flowering time may be required. In order to detect targets for future marker-assisted improvement and validate the practical application of GS for wheat breeding we genotyped 376 elite wheat varieties with 3,046 DArT, single nucleotide polymorphism and gene markers and measured seven traits in replicated yield trials over 2 years in France, Germany and the UK. The scale of the phenotyping exceeds the breadth of previous AM and GS studies in these key economic wheat production regions of Northern Europe. Mixed-linear modelling (MLM) detected significant marker-trait associations across and within regions. Genomic prediction using elastic net gave low to high prediction accuracies depending on the trait, and could be experimentally increased by modifying the constituents of the training population (TP). We also tested the use of differentially penalised regression to integrate candidate gene and genome-wide markers to predict traits, demonstrating the validity and simplicity of this approach. Overall, our results suggest that whilst AM offers potential for application in both research and breeding, GS represents an exciting opportunity to select key traits, and that optimisation of the TP is crucial to its successful implementation.

Association mapping for five agronomic traits in the common bean (Phaseolus vulgaris L.)

BACKGROUND

The common bean is the most important grain legume and a major source of protein in many developing countries. We analysed the following traits: pod fibre (PF), seeds per pod (SPP), plant type (PT), growth habit (GH), and days to flowering (DF) for a set of diverse common bean accessions and determined whether such traits were associated with amplified fragment length polymorphism (AFLP), simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers.

RESULTS

In this study, 66 common bean genotypes were used and genotyped with 233 AFLP, 105 SNP and 80 SSR markers. The association analysis between markers and five traits was performed using a General Linear Model (GLM) in Trait Analysis by aSSociation, Evolution and Linkage (TASSEL). The population structure was determined using the STRUCTURE software, and seven groups (K = 7) were identified among genotypes. The associations for such traits were identified and quantified; 62 markers were associated with the five traits.

CONCLUSION

This study demonstrated that association mapping using a reasonable number of markers, distributed across the genome and with the appropriate number of individuals harboured to detect DNA markers linked to the traits of PF, SPP, PT, GH and DF in common bean.

A genome-wide scan of selective sweeps and association mapping of fruit traits using microsatellite markers in watermelon

Our genetic diversity study uses microsatellites of known map position to estimate genome level population structure and linkage disequilibrium, and to identify genomic regions that have undergone selection during watermelon domestication and improvement. Thirty regions that showed evidence of selective sweep were scanned for the presence of candidate genes using the watermelon genome browser (www.icugi.org). We localized selective sweeps in intergenic regions, close to the promoters, and within the exons and introns of various genes. This study provided an evidence of convergent evolution for the presence of diverse ecotypes with special reference to American and European ecotypes. Our search for location of linked markers in the whole-genome draft sequence revealed that BVWS00358, a GA repeat microsatellite, is the GAGA type transcription factor located in the 5' untranslated regions of a structure and insertion element that expresses a Cys2His2 Zinc finger motif, with presumed biological processes related to chitin response and transcriptional regulation. In addition, BVWS01708, an ATT repeat microsatellite, located in the promoter of a DTW domain-containing protein (Cla002761); and 2 other simple sequence repeats that association mapping link to fruit length and rind thickness.

A genome-wide association study of limb bone length using a Large White × Minzhu intercross population

Background

In pig, limb bone length influences ham yield and body height to a great extent and has important economic implications for pig industry. In this study, an intercross population was constructed between the indigenous Chinese Minzhu pig breed and the western commercial Large White pig breed to examine the genetic basis for variation in limb bone length. The aim of this study was to detect potential genetic variants associated with porcine limb bone length.

Methods

A total of 571 F2 individuals from a Large White and Minzhu intercross population were genotyped using the Illumina PorcineSNP60K Beadchip, and phenotyped for femur length (FL), humerus length (HL), hipbone length (HIPL), scapula length (SL), tibia length (TL), and ulna length (UL). A genome-wide association study was performed by applying the previously reported approach of genome-wide rapid association using mixed model and regression. Statistical significance of the associations was based on Bonferroni-corrected P-values.

Results

A total of 39 significant SNPs were mapped to a 11.93 Mb long region on pig chromosome 7 (SSC7). Linkage analysis of these significant SNPs revealed three haplotype blocks of 495 kb, 376 kb and 492 kb, respectively, in the 11.93 Mb region. Annotation based on the pig reference genome identified 15 genes that were located near or contained the significant SNPs in these linkage disequilibrium intervals. Conditioned analysis revealed that four SNPs, one on SSC2 and three on SSC4, showed significant associations with SL and HL, respectively.

Conclusions

Analysis of the 15 annotated genes that were identified in these three haplotype blocks indicated that HMGA1 and PPARD, which are expressed in limbs and influence chondrocyte cell growth and differentiation, could be considered as relevant biological candidates for limb bone length in pig, with potential applications in breeding programs. Our results may also be useful for the study of the mechanisms that underlie human limb length and body height.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-014-0056-6) contains supplementary material, which is available to authorized users.

Genome-wide prediction methods in highly diverse and heterozygous species: proof-of-concept through simulation in grapevine

Nowadays, genome-wide association studies (GWAS) and genomic selection (GS) methods which use genome-wide marker data for phenotype prediction are of much potential interest in plant breeding. However, to our knowledge, no studies have been performed yet on the predictive ability of these methods for structured traits when using training populations with high levels of genetic diversity. Such an example of a highly heterozygous, perennial species is grapevine. The present study compares the accuracy of models based on GWAS or GS alone, or in combination, for predicting simple or complex traits, linked or not with population structure. In order to explore the relevance of these methods in this context, we performed simulations using approx 90,000 SNPs on a population of 3,000 individuals structured into three groups and corresponding to published diversity grapevine data. To estimate the parameters of the prediction models, we defined four training populations of 1,000 individuals, corresponding to these three groups and a core collection. Finally, to estimate the accuracy of the models, we also simulated four breeding populations of 200 individuals. Although prediction accuracy was low when breeding populations were too distant from the training populations, high accuracy levels were obtained using the sole core-collection as training population. The highest prediction accuracy was obtained (up to 0.9) using the combined GWAS-GS model. We thus recommend using the combined prediction model and a core-collection as training population for grapevine breeding or for other important economic crops with the same characteristics.

Enrichment of statistical power for genome-wide association studies

BACKGROUND

The inheritance of most human diseases and agriculturally important traits is controlled by many genes with small effects. Identifying these genes, while simultaneously controlling false positives, is challenging. Among available statistical methods, the mixed linear model (MLM) has been the most flexible and powerful for controlling population structure and individual unequal relatedness (kinship), the two common causes of spurious associations. The introduction of the compressed MLM (CMLM) method provided additional opportunities for optimization by adding two new model parameters: grouping algorithms and number of groups.

RESULTS

This study introduces another model parameter to develop an enriched CMLM (ECMLM). The parameter involves algorithms to define kinship between groups (that is, kinship algorithms). The ECMLM calculates kinship using several different algorithms and then chooses the best combination between kinship algorithms and grouping algorithms.

CONCLUSION

Simulations show that the ECMLM increases statistical power. In some cases, the magnitude of power gained by using ECMLM instead of CMLM is larger than the improvement found by using CMLM instead of MLM.

Combining high-throughput phenotyping and genome-wide association studies to reveal natural genetic variation in rice

Even as the study of plant genomics rapidly develops through the use of high-throughput sequencing techniques, traditional plant phenotyping lags far behind. Here we develop a high-throughput rice phenotyping facility (HRPF) to monitor 13 traditional agronomic traits and 2 newly defined traits during the rice growth period. Using genome-wide association studies (GWAS) of the 15 traits, we identify 141 associated loci, 25 of which contain known genes such as the Green Revolution semi-dwarf gene, SD1. Based on a performance evaluation of the HRPF and GWAS results, we demonstrate that high-throughput phenotyping has the potential to replace traditional phenotyping techniques and can provide valuable gene identification information. The combination of the multifunctional phenotyping tools HRPF and GWAS provides deep insights into the genetic architecture of important traits.

Combining high-throughput phenotyping and genome-wide association studies to reveal natural genetic variation in rice

Even as the study of plant genomics rapidly develops through the use of high-throughput sequencing techniques, traditional plant phenotyping lags far behind. Here we develop a high-throughput rice phenotyping facility (HRPF) to monitor 13 traditional agronomic traits and 2 newly defined traits during the rice growth period. Using genome-wide association studies (GWAS) of the 15 traits, we identify 141 associated loci, 25 of which contain known genes such as the Green Revolution semi-dwarf gene, SD1. Based on a performance evaluation of the HRPF and GWAS results, we demonstrate that high-throughput phenotyping has the potential to replace traditional phenotyping techniques and can provide valuable gene identification information. The combination of the multifunctional phenotyping tools HRPF and GWAS provides deep insights into the genetic architecture of important traits.

A SUPER powerful method for genome wide association study

Genome-Wide Association Studies shed light on the identification of genes underlying human diseases and agriculturally important traits. This potential has been shadowed by false positive findings. The Mixed Linear Model (MLM) method is flexible enough to simultaneously incorporate population structure and cryptic relationships to reduce false positives. However, its intensive computational burden is prohibitive in practice, especially for large samples. The newly developed algorithm, FaST-LMM, solved the computational problem, but requires that the number of SNPs be less than the number of individuals to derive a rank-reduced relationship. This restriction potentially leads to less statistical power when compared to using all SNPs. We developed a method to extract a small subset of SNPs and use them in FaST-LMM. This method not only retains the computational advantage of FaST-LMM, but also remarkably increases statistical power even when compared to using the entire set of SNPs. We named the method SUPER (Settlement of MLM Under Progressively Exclusive Relationship) and made it available within an implementation of the GAPIT software package.

Identification and validation of functional markers in a global rice collection by association mapping

Recent results indicate that marker-assisted selection is an effective approach to reduce the cost and to improve the efficacy and accuracy of selection in plant breeding. This study was conducted to identify and validate molecular markers linked to important breeding traits by association mapping. The association was evaluated between 81 molecular markers (STS, SSR, Indel, CAPS, and PCR-based SNP) and 15 morphological traits in a global panel of 100 rice (Oryza sativa) accessions. The population structure analysis identified three main subpopulations. Obvious kinship relationships were also detected between the rice accessions. Association analysis was performed based on the mixed linear model by considering population structure and family relatedness. In addition, the false discovery rate method was used to correct the multiple testing. A total of 47 marker-trait associations were identified, including 22 markers for 14 traits. Among all, the polymorphism at the loci DDR-GL was highly associated with grain characters (grain length, grain width, and length/width ratio). In addition, marker RM3148 was responsible for five important traits simultaneously. Results demonstrated that such informative markers can be very useful for rice breeding programs using marker-assisted selection. Moreover, the diverse populations of rice accessions are a valuable resource for association mapping of morphological traits.

Genome-wide association studies identifies seven major regions responsible for iron deficiency chlorosis in soybean (Glycine max)

Iron deficiency chlorosis (IDC) is a yield limiting problem in soybean (Glycine max (L.) Merr) production regions with calcareous soils. Genome-wide association study (GWAS) was performed using a high density SNP map to discover significant markers, QTL and candidate genes associated with IDC trait variation. A stepwise regression model included eight markers after considering LD between markers, and identified seven major effect QTL on seven chromosomes. Twelve candidate genes known to be associated with iron metabolism mapped near these QTL supporting the polygenic nature of IDC. A non-synonymous substitution with the highest significance in a major QTL region suggests soybean orthologs of FRE1 on Gm03 is a major gene responsible for trait variation. NAS3, a gene that encodes the enzyme nicotianamine synthase which synthesizes the iron chelator nicotianamine also maps to the same QTL region. Disease resistant genes also map to the major QTL, supporting the hypothesis that pathogens compete with the plant for Fe and increase iron deficiency. The markers and the allelic combinations identified here can be further used for marker assisted selection.

Genome wide association studies using a new nonparametric model reveal the genetic architecture of 17 agronomic traits in an enlarged maize association panel

Association mapping is a powerful approach for dissecting the genetic architecture of complex quantitative traits using high-density SNP markers in maize. Here, we expanded our association panel size from 368 to 513 inbred lines with 0.5 million high quality SNPs using a two-step data-imputation method which combines identity by descent (IBD) based projection and k-nearest neighbor (KNN) algorithm. Genome-wide association studies (GWAS) were carried out for 17 agronomic traits with a panel of 513 inbred lines applying both mixed linear model (MLM) and a new method, the Anderson-Darling (A-D) test. Ten loci for five traits were identified using the MLM method at the Bonferroni-corrected threshold −log₁₀ (P) >5.74 (α = 1). Many loci ranging from one to 34 loci (107 loci for plant height) were identified for 17 traits using the A-D test at the Bonferroni-corrected threshold −log₁₀ (P) >7.05 (α = 0.05) using 556809 SNPs. Many known loci and new candidate loci were only observed by the A-D test, a few of which were also detected in independent linkage analysis. This study indicates that combining IBD based projection and KNN algorithm is an efficient imputation method for inferring large missing genotype segments. In addition, we showed that the A-D test is a useful complement for GWAS analysis of complex quantitative traits. Especially for traits with abnormal phenotype distribution, controlled by moderate effect loci or rare variations, the A-D test balances false positives and statistical power. The candidate SNPs and associated genes also provide a rich resource for maize genetics and breeding.

Genotype imputation has been used widely in the analysis of genome-wide association studies (GWAS) to boost power and fine-map associations. We developed a two-step data imputation method to meet the challenge of large proportion missing genotypes. GWAS have uncovered an extensive genetic architecture of complex quantitative traits using high-density SNP markers in maize in the past few years. Here, GWAS were carried out for 17 agronomic traits with a panel of 513 inbred lines applying both mixed linear model and a new method, the Anderson-Darling (A-D) test. We intend to show that the A-D test is a complement to current GWAS methods, especially for complex quantitative traits controlled by moderate effect loci or rare variations and with abnormal phenotype distribution. In addition, the traits associated QTL identified here provide a rich resource for maize genetics and breeding.

Single nucleotide polymorphisms generated by genotyping by sequencing to characterize genome-wide diversity, linkage disequilibrium, and selective sweeps in cultivated watermelon

Background

A large single nucleotide polymorphism (SNP) dataset was used to analyze genome-wide diversity in a diverse collection of watermelon cultivars representing globally cultivated, watermelon genetic diversity. The marker density required for conducting successful association mapping depends on the extent of linkage disequilibrium (LD) within a population. Use of genotyping by sequencing reveals large numbers of SNPs that in turn generate opportunities in genome-wide association mapping and marker-assisted selection, even in crops such as watermelon for which few genomic resources are available. In this paper, we used genome-wide genetic diversity to study LD, selective sweeps, and pairwise F_ST distributions among worldwide cultivated watermelons to track signals of domestication.

Results

We examined 183 Citrullus lanatus var. lanatus accessions representing domesticated watermelon and generated a set of 11,485 SNP markers using genotyping by sequencing. With a diverse panel of worldwide cultivated watermelons, we identified a set of 5,254 SNPs with a minor allele frequency of ≥ 0.05, distributed across the genome. All ancestries were traced to Africa and an admixture of various ancestries constituted secondary gene pools across various continents. A sliding window analysis using pairwise F_ST values was used to resolve selective sweeps. We identified strong selection on chromosomes 3 and 9 that might have contributed to the domestication process. Pairwise analysis of adjacent SNPs within a chromosome as well as within a haplotype allowed us to estimate genome-wide LD decay. LD was also detected within individual genes on various chromosomes. Principal component and ancestry analyses were used to account for population structure in a genome-wide association study. We further mapped important genes for soluble solid content using a mixed linear model.

Conclusions

Information concerning the SNP resources, population structure, and LD developed in this study will help in identifying agronomically important candidate genes from the genomic regions underlying selection and for mapping quantitative trait loci using a genome-wide association study in sweet watermelon.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-767) contains supplementary material, which is available to authorized users.