A genome‐wide association study of 23 agronomic traits in Chinese wheat landraces

Uncovering the genetic basis of agronomic traits in wheat landraces is important for ensuring global food security via the development of improved varieties. Here, 723 wheat landraces from 10 Chinese agro‐ecological zones were evaluated for 23 agronomic traits in six environments. All accessions could be clustered into five subgroups based on phenotypic data via discriminant function analysis, which was highly consistent with genotypic classification. A genome‐wide association study was conducted for these traits using 52 303 DArT‐seq markers to identify marker‐trait associations and candidate genes. Using both the general linear model and the mixed linear model, 149 significant markers were identified for 21 agronomic traits based on best linear unbiased prediction values. Considering the linkage disequilibrium decay distance in this study, significant markers within 10 cM were combined as a quantitative trait locus (QTL), with a total of 29 QTL identified for 15 traits. Of these, five QTL for heading date, flag leaf width, peduncle length, and thousand kernel weight had been reported previously. Twenty‐five candidate genes associated with significant markers were identified. These included the known vernalization genes VRN‐B1 and vrn‐B3 and the photoperiod response genes Ppd and PRR. Overall, this study should be helpful in elucidating the underlying genetic mechanisms of complex agronomic traits and performing marker‐assisted selection in wheat.     

Variation in recombination rate may bias human genetic disease mapping studies

The availability of the human genome sequence and variability information (as from the International HapMap project) will enhance our ability to map genetic disorders and choose targets for therapeutic intervention. However, several factors, such as regional variation in recombination rate, can bias conclusions from genetic mapping studies. Here, we examine the impact of regional variation in recombination rate across the human genome. Through computer simulations and literature surveys, we conclude that genetic disorders have been mapped to regions of low recombination more often than expected if such diseases were randomly distributed across the genome. This concentration in low recombination regions may be an artifact, and disorders appearing to be caused by a few genes of large effect may be polygenic. Future genetic mapping studies should be conscious of this potential complication by noting the regional recombination rate of regions implicated in diseases.     

Modelling the genetic architecture of flowering time control in barley through nested association mapping

Background

Barley, globally the fourth most important cereal, provides food and beverages for humans and feed for animal husbandry. Maximizing grain yield under varying climate conditions largely depends on the optimal timing of flowering. Therefore, regulation of flowering time is of extraordinary importance to meet future food and feed demands. We developed the first barley nested association mapping (NAM) population, HEB-25, by crossing 25 wild barleys with one elite barley cultivar, and used it to dissect the genetic architecture of flowering time.

Results

Upon cultivation of 1,420 lines in multi-field trials and applying a genome-wide association study, eight major quantitative trait loci (QTL) were identified as main determinants to control flowering time in barley. These QTL accounted for 64% of the cross-validated proportion of explained genotypic variance (p_G). The strongest single QTL effect corresponded to the known photoperiod response gene Ppd-H1. After sequencing the causative part of Ppd-H1, we differentiated twelve haplotypes in HEB-25, whereof the strongest exotic haplotype accelerated flowering time by 11 days compared to the elite barley haplotype. Applying a whole genome prediction model including main effects and epistatic interactions allowed predicting flowering time with an unmatched accuracy of 77% of cross-validated p_G.

Conclusions

The elaborated causal models represent a fundamental step to explain flowering time in barley. In addition, our study confirms that the exotic biodiversity present in HEB-25 is a valuable toolbox to dissect the genetic architecture of important agronomic traits and to replenish the elite barley breeding pool with favorable, trait-improving exotic alleles.

     

Metabolome‐genome‐wide association study dissects genetic architecture for generating natural variation in rice secondary metabolism

Plants produce structurally diverse secondary (specialized) metabolites to increase their fitness for survival under adverse environments. Several bioactive compounds for new drugs have been identified through screening of plant extracts. In this study, genome‐wide association studies (GWAS) were conducted to investigate the genetic architecture behind the natural variation of rice secondary metabolites. GWAS using the metabolome data of 175 rice accessions successfully identified 323 associations among 143 single nucleotide polymorphisms (SNPs) and 89 metabolites. The data analysis highlighted that levels of many metabolites are tightly associated with a small number of strong quantitative trait loci (QTLs). The tight association may be a mechanism generating strains with distinct metabolic composition through the crossing of two different strains. The results indicate that one plant species produces more diverse phytochemicals than previously expected, and plants still contain many useful compounds for human applications.     

Wild GWAS—association mapping in natural populations

The increasing affordability of sequencing and genotyping technologies has transformed the field of molecular ecology in recent decades. By correlating marker variants with trait variation using association analysis, large‐scale genotyping and phenotyping of individuals from wild populations has enabled the identification of genomic regions that contribute to phenotypic differences among individuals. Such “gene mapping” studies are enabling us to better predict evolutionary potential and the ability of populations to adapt to challenges, such as changing environment. These studies are also allowing us to gain insight into the evolutionary processes maintaining variation in natural populations, to better understand genotype‐by‐environment and epistatic interactions and to track the dynamics of allele frequency change at loci contributing to traits under selection. Gene mapping in the wild using genomewide association scans (GWAS) do, however, come with a number of methodological challenges, not least the population structure in space and time inherent to natural populations. We here provide an overview of these challenges, summarize the exciting methodological advances and applications of association mapping in natural populations reported in this special issue and provide some guidelines for future “wild GWAS” research.     

The relationship between parental genetic or phenotypic divergence and progeny variation in the maize nested association mapping population

Appropriate selection of parents for the development of mapping populations is pivotal to maximizing the power of quantitative trait loci detection. Trait genotypic variation within a family is indicative of the family's informativeness for genetic studies. Accurate prediction of the most useful parental combinations within a species would help guide quantitative genetics studies. We tested the reliability of genotypic and phenotypic distance estimators between pairs of maize inbred lines to predict genotypic variation for quantitative traits within families derived from biparental crosses. We developed 25 families composed of ~200 random recombinant inbred lines each from crosses between a common reference parent inbred, B73, and 25 diverse maize inbreds. Parents and families were evaluated for 19 quantitative traits across up to 11 environments. Genetic distances (GDs) among parents were estimated with 44 simple sequence repeat and 2303 single-nucleotide polymorphism markers. GDs among parents had no predictive value for progeny variation, which is most likely due to the choice of neutral markers. In contrast, we observed for about half of the traits measured a positive correlation between phenotypic parental distances and within-family genetic variance estimates. Consequently, the choice of promising segregating populations can be based on selecting phenotypically diverse parents. These results are congruent with models of genetic architecture that posit numerous genes affecting quantitative traits, each segregating for allelic series, with dispersal of allelic effects across diverse genetic material. This architecture, common to many quantitative traits in maize, limits the predictive value of parental genotypic or phenotypic values on progeny variance.     

A genome‐wide association study reveals candidate genes for the supernumerary nipple phenotype in sheep (Ovis aries)

Genome‐wide association studies (GWASs) have been widely applied in livestock to identify genes associated with traits of economic interest. Here, we conducted the first GWAS of the supernumerary nipple phenotype in Wadi sheep, a native Chinese sheep breed, based on Ovine Infinium HD SNP BeadChip genotypes in a total of 144 ewes (75 cases with four teats, including two normal and two supernumerary teats, and 69 control cases with two teats). We detected 63 significant SNPs at the chromosome‐wise threshold. Additionally, one candidate region (chr1: 170.723–170.734 Mb) was identified by haplotype‐based association tests, with one SNP (rs413490006) surrounding functional genes BBX and CD47 on chromosome 1 being commonly identified as significant by the two mentioned analyses. Moreover, Gene Ontology enrichment for the significant SNPs identified by the GWAS analysis was functionally clustered into the categories of receptor activity and synaptic membrane. In addition, pathway mapping revealed four promising pathways (Wnt, oxytocin, MAPK and axon guidance) involved in the development of the supernumerary nipple phenotype. Our results provide novel and important insights into the genetic mechanisms underlying the phenotype of supernumerary nipples in mammals, including humans. These findings may be useful for future breeding and genetics in sheep and other livestock.     

Metabolite‐based genome‐wide association study enables dissection of the flavonoid decoration pathway of wheat kernels

The marriage of metabolomic approaches with genetic design has proven a powerful tool in dissecting diversity in the metabolome and has additionally enhanced our understanding of complex traits. That said, such studies have rarely been carried out in wheat. In this study, we detected 805 metabolites from wheat kernels and profiled their relative contents among 182 wheat accessions, conducting a metabolite‐based genome‐wide association study (mGWAS) utilizing 14 646 previously described polymorphic SNP markers. A total of 1098 mGWAS associations were detected with large effects, within which 26 candidate genes were tentatively designated for 42 loci. Enzymatic assay of two candidates indicated they could catalyse glucosylation and subsequent malonylation of various flavonoids and thereby the major flavonoid decoration pathway of wheat kernel was dissected. Moreover, numerous high‐confidence genes associated with metabolite contents have been provided, as well as more subdivided metabolite networks which are yet to be explored within our data. These combined efforts presented the first step towards realizing metabolomics‐associated breeding of wheat.     

Nested‐association mapping (NAM)‐based genetic dissection uncovers candidate genes for seed and pod weights in peanut (Arachis hypogaea)

Multiparental genetic mapping populations such as nested‐association mapping (NAM) have great potential for investigating quantitative traits and associated genomic regions leading to rapid discovery of candidate genes and markers. To demonstrate the utility and power of this approach, two NAM populations, NAM_Tifrunner and NAM_Florida‐07, were used for dissecting genetic control of 100‐pod weight (PW) and 100‐seed weight (SW) in peanut. Two high‐density SNP‐based genetic maps were constructed with 3341 loci and 2668 loci for NAM_Tifrunner and NAM_Florida‐07, respectively. The quantitative trait locus (QTL) analysis identified 12 and 8 major effect QTLs for PW and SW, respectively, in NAM_Tifrunner, and 13 and 11 major effect QTLs for PW and SW, respectively, in NAM_Florida‐07. Most of the QTLs associated with PW and SW were mapped on the chromosomes A05, A06, B05 and B06. A genomewide association study (GWAS) analysis identified 19 and 28 highly significant SNP–trait associations (STAs) in NAM_Tifrunner and 11 and 17 STAs in NAM_Florida‐07 for PW and SW, respectively. These significant STAs were co‐localized, suggesting that PW and SW are co‐regulated by several candidate genes identified on chromosomes A05, A06, B05, and B06. This study demonstrates the utility of NAM population for genetic dissection of complex traits and performing high‐resolution trait mapping in peanut.     

Analysis of the genetic architecture of maize kernel size traits by combined linkage and association mapping

Kernel size‐related traits are the most direct traits correlating with grain yield. The genetic basis of three kernel traits of maize, kernel length (KL), kernel width (KW) and kernel thickness (KT), was investigated in an association panel and a biparental population. A total of 21 single nucleotide polymorphisms (SNPs) were detected to be most significantly (P < 2.25 × 10^?6) associated with these three traits in the association panel under four environments. Furthermore, 50 quantitative trait loci (QTL) controlling these traits were detected in seven environments in the intermated B73 × Mo17 (IBM) Syn10 doubled haploid (DH) population, of which eight were repetitively identified in at least three environments. Combining the two mapping populations revealed that 56 SNPs (P < 1 × 10^?3) fell within 18 of the QTL confidence intervals. According to the top significant SNPs, stable‐effect SNPs and the co‐localized SNPs by association analysis and linkage mapping, a total of 73 candidate genes were identified, regulating seed development. Additionally, seven miRNAs were found to situate within the linkage disequilibrium (LD) regions of the co‐localized SNPs, of which zma‐miR164e was demonstrated to cleave the mRNAs of Arabidopsis CUC1, CUC2 and NAC6 in vitro. Overexpression of zma‐miR164e resulted in the down‐regulation of these genes above and the failure of seed formation in Arabidopsis pods, with the increased branch number. These findings provide insights into the mechanism of seed development and the improvement of molecular marker‐assisted selection (MAS) for high‐yield breeding in maize.     

Genome‐wide association study of six quality traits reveals the association of the TaRPP13L1 gene with flour colour in Chinese bread wheat

Flour colour, kernel hardness, grain protein content and wet gluten content are important quality properties that determine end use in bread wheat. Here, a wheat 90K genotyping assay was used for a genome‐wide association study (GWAS) of the six quality‐related traits in Chinese wheat cultivars in eight environments over four years. A total of 846 significant single nucleotide polymorphisms (SNPs) were identified, explaining approximately 30% of the phenotypic variation on average, and 103 multienvironment‐significant SNPs were detected in more than four environments. Quantitative trait loci (QTL) mapping in the biparent population confirmed some important SNP loci. Moreover, it was determined that some important genes were associated with the six quality traits, including some known functional genes and annotated unknown functional genes. Of the annotated unknown functional genes, it was verified that TaRPP13L1 was associated with flour colour. Wheat cultivars or lines with TaRPP13L1‐B1a showed extremely significantly higher flour redness and lower yellowness than those with TaRPP13L1‐B1b in the Chinese wheat natural population and the doubled haploid (DH) population. Two tetraploid wheat lines with premature stop codons of the TaRPP13L1 gene mutagenized by ethyl methanesulfonate (EMS) showed extremely significantly higher flour redness and lower yellowness than wild type. Our data suggest that the TaRPP13L1 gene plays an important role in modulating wheat flour colour. This study provides useful information for further dissection of the genetic basis of flour colour and also provides valuable genes or genetic loci for marker‐assisted selection to improve the process of breeding quality wheat in China.     

Mapping QTL for production traits in segregating Piétrain pig populations using genome‐wide association study results of F2 crosses

In this study, genome‐wide association study (GWAS) results of porcine F2 crosses were used to map QTL in outcross Piétrain populations. For this purpose, two F2 crosses (Piétrain × Meishan, n = 304; Piétrain × Wild Boar, n = 291) were genotyped with the PorcineSNP60v2 BeadChip and phenotyped for the dressing yield, carcass length, daily gain and drip loss traits. GWASs were conducted in the pooled F2 cross applying single marker mixed linear models. For the investigated traits, between two and five (in total 15) QTL core regions, spanning 250 segregating SNPs around a significant trait‐associated peak SNP, were identified. The SNPs within the QTL core regions were subsequently tested for trait association in two outcross Piétrain populations consisting of 771 progeny‐tested boars and 210 sows with their own performance records. In the sow (boar) dataset, five (eight) of the 15 mapped QTL were validated. Hence, many QTL mapped in the F2 crosses (with Piétrain as a common founder breed) are still segregating in the current Piétrain breed. This confirms the usefulness of existing F2 crosses for mapping QTL that are still segregating in the recent founder breed generation. The approach utilizes the high power of an F2 cross to map QTL in a breeding population for which it is not guaranteed that they would be found using a GWAS in this population.     

Population structure and genetic basis of the agronomic traits of upland cotton in China revealed by a genome‐wide association study using high‐density SNPs

Gossypium hirsutum L. represents the largest source of textile fibre, and China is one of the largest cotton‐producing and cotton‐consuming countries in the world. To investigate the genetic architecture of the agronomic traits of upland cotton in China, a diverse and nationwide population containing 503 G. hirsutum accessions was collected for a genome‐wide association study (GWAS) on 16 agronomic traits. The accessions were planted in four places from 2012 to 2013 for phenotyping. The CottonSNP63K array and a published high‐density map based on this array were used for genotyping. The 503 G. hirsutum accessions were divided into three subpopulations based on 11 975 quantified polymorphic single‐nucleotide polymorphisms (SNPs). By comparing the genetic structure and phenotypic variation among three genetic subpopulations, seven geographic distributions and four breeding periods, we found that geographic distribution and breeding period were not the determinants of genetic structure. In addition, no obvious phenotypic differentiations were found among the three subpopulations, even though they had different genetic backgrounds. A total of 324 SNPs and 160 candidate quantitative trait loci (QTL) regions were identified as significantly associated with the 16 agronomic traits. A network was established for multieffects in QTLs and interassociations among traits. Thirty‐eight associated regions had pleiotropic effects controlling more than one trait. One candidate gene, Gh_D08G2376, was speculated to control the lint percentage (LP). This GWAS is the first report using high‐resolution SNPs in upland cotton in China to comprehensively investigate agronomic traits, and it provides a fundamental resource for cotton genetic research and breeding.     

Linkage disequilibrium pattern and genome‐wide association mapping for meat traits in multiple porcine F2 crosses

In the present study, data from four F2 crosses were analysed and used to study the linkage disequilibrium (LD) structure within and across the crosses. Genome‐wide association analyses (GWASes) for conductivity and dressing out meat traits were conducted using single‐marker and Bayesian multi‐marker models using the pooled data from all F2 crosses. Porcine F2 crosses generated from the distantly related founder breeds Wild Boar, Piétrain and Meishan, as well as from a porcine F2 cross from the closely related founder breed Piétrain and an F1 Large White × Landrace cross were pooled. A total of 2572 F2 animals were genotyped using a 62K SNP chip. The positions of the SNPs were based on genome assembly Sscrofa11.1. After post‐alignment and genotype filtering, approximately 50K SNPs were usable for LD studies and GWASes. The main findings of the present study are that the breakdown of LD was faster in crosses from closely related founder breeds compared to crosses from distantly related founders. The fastest breakdown of LD was observed by pooling the data. Based on the single‐marker results and LD structure, clusters and windows were built for 1‐Mb intervals. For conductivity and dressing out, 183 and 191 nominal significant associations respectively and six and five clusters respectively were found. Dominance was important for conductivity, and considering dominance in GWASes improved the mapping signals. Most clear signals were found for conductivity on SSC6, 8 and 15 and for dressing out on SSC2 and 7. Considering dominance might contribute to the accuracy of genomic selection and serve as a guide for choosing mating pairs with good combining abilities. However, further research is needed to investigate if dominance is also important in crossbreed pig breeding schemes.     

Dissection of genetic architecture for glucosinolate accumulations in leaves and seeds of Brassica napus by genome‐wide association study

Glucosinolates (GSLs), whose degradation products have been shown to be increasingly important for human health and plant defence, compose important secondary metabolites found in the order Brassicales. It is highly desired to enhance pest and disease resistance by increasing the leaf GSL content while keeping the content low in seeds of Brassica napus, one of the most important oil crops worldwide. Little is known about the regulation of GSL accumulation in the leaves. We quantified the levels of 9 different GSLs and 15 related traits in the leaves of 366 accessions and found that the seed and leaf GSL content were highly correlated (r = 0.79). A total of 78 loci were associated with GSL traits, and five common and eleven tissue‐specific associated loci were related to total leaf and seed GSL content. Thirty‐six candidate genes were inferred to be involved in GSL biosynthesis. The candidate gene BnaA03g40190D (BnaA3.MYB28) was validated by DNA polymorphisms and gene expression analysis. This gene was responsible for high leaf/low seed GSL content and could explain 30.62% of the total leaf GSL variation in the low seed GSL panel and was not fixed during double‐low rapeseed breeding. Our results provide new insights into the genetic basis of GSL variation in leaves and seeds and may facilitate the metabolic engineering of GSLs and the breeding of high leaf/low seed GSL content in B. napus.     

Quantitative trait locus mapping reveals complex genetic architecture of quantitative virulence in the wheat pathogen Zymoseptoria tritici

We conducted a comprehensive analysis of virulence in the fungal wheat pathogen Zymoseptoria tritici using quantitative trait locus (QTL) mapping. High‐throughput phenotyping based on automated image analysis allowed the measurement of pathogen virulence on a scale and with a precision that was not previously possible. Across two mapping populations encompassing more than 520 progeny, 540 710 pycnidia were counted and their sizes and grey values were measured. A significant correlation was found between pycnidia size and both spore size and number. Precise measurements of percentage leaf area covered by lesions provided a quantitative measure of host damage. Combining these large and accurate phenotypic datasets with a dense panel of restriction site‐associated DNA sequencing (RADseq) genetic markers enabled us to genetically dissect pathogen virulence into components related to host damage and those related to pathogen reproduction. We showed that different components of virulence can be under separate genetic control. Large‐ and small‐effect QTLs were identified for all traits, with some QTLs specific to mapping populations, cultivars and traits and other QTLs shared among traits within the same mapping population. We associated the presence of four accessory chromosomes with small, but significant, increases in several virulence traits, providing the first evidence for a meaningful function associated with accessory chromosomes in this organism. A large‐effect QTL involved in host specialization was identified on chromosome 7, leading to the identification of candidate genes having a large effect on virulence.     

Genetic‐based dissection of arsenic accumulation in maize using a genome‐wide association analysis method

Understanding the mechanism of arsenic (As) accumulation in plants is important in reducing As's toxicity to plants and its potential risks to human health. Here, we performed a genome‐wide association study to dissect the genetic basis of the As contents of different maize tissues in Xixian, which was irrigated with As‐rich surface water, and Changge using an association population consisting of 230 representative maize inbred lines. Phenotypic data revealed a wide normal distribution and high repeatability for the As contents in maize tissues. The As concentrations in maize tissues followed the same trend in the two locations: kernels < axes < stems < bracts < leaves. In total, 15, 16 and 15 non‐redundant quantitative trait loci (QTL s) associated with As concentrations were identified (P  ≤ 2.04 × 10^?6) in five tissues from Xixian, Changge, and the combination of the locations, respectively, explaining 9.70%–24.65% of the phenotypic variation for each QTL , on average. Additionally, four QTL s [involving 15 single nucleotide polymorphisms (SNP s)] were detected in the single and the combined locations, indicating that these loci/SNP s might be stable across different environments. The candidate genes associated with these four loci were predicted. In addition, four non‐redundant QTL s (6 SNP s), including a QTL that was detected in multiple locations according to the genome‐wide association study, were found to co‐localize with four previously reported QTL intervals. These results are valuable to understand the genetic architecture of As mechanism in maize and facilitate the genetic improvement of varieties without As toxicity.     

Genome‐wide association mapping of time‐dependent growth responses to moderate drought stress in Arabidopsis

Large areas of arable land are often confronted with irregular rainfall resulting in limited water availability for part(s) of the growing seasons, which demands research for drought tolerance of plants. Natural variation was observed for biomass accumulation upon controlled moderate drought stress in 324 natural accessions of Arabidopsis. Improved performance under drought stress was correlated with early flowering and lack of vernalization requirement, indicating overlap in the regulatory networks of flowering time and drought response or correlated responses of these traits to natural selection. In addition, plant size was negatively correlated with relative water content (RWC) independent of the absolute water content (WC), indicating a prominent role for soluble compounds. Growth in control and drought conditions was determined over time and was modelled by an exponential function. Genome‐wide association (GWA) mapping of temporal plant size data and of model parameters resulted in the detection of six time‐dependent quantitative trait loci (QTLs) strongly associated with drought. Most QTLs would not have been identified if plant size was determined at a single time point. Analysis of earlier reported gene expression changes upon drought enabled us to identify for each QTL the most likely candidates.