Genome-Wide Association Study Reveals Novel Powdery Mildew Resistance Loci in Bread Wheat

Powdery mildew (PM), caused by the fungal pathogen Blumeria graminis f. sp. tritici (Bgt), significantly threatens global bread wheat production. Although the use of resistant cultivars is an effective strategy for managing PM, currently available wheat cultivars lack sufficient levels of resistance. To tackle this challenge, we conducted a comprehensive genome-wide association study (GWAS) using a diverse panel of 286 bread wheat genotypes. Over three consecutive years (2020-2021, 2021-2022, and 2022-2023), these genotypes were extensively evaluated for PM severity under field conditions following inoculation with virulent Bgt isolates. The panel was previously genotyped using the Illumina 90K Infinium iSelect assay to obtain genome-wide single-nucleotide polymorphism (SNP) marker coverage. By applying FarmCPU, a multilocus mixed model, we identified a total of 113 marker-trait associations (MTAs) located on chromosomes 1A, 1B, 2B, 3A, 3B, 4A, 4B, 5A, 5B, 6B, 7A, and 7B at a significance level of p ≤ 0.001. Notably, four novel MTAs on chromosome 6B were consistently detected in 2020-2021 and 2021-2022. Furthermore, within the confidence intervals of the identified SNPs, we identified 96 candidate genes belonging to different proteins including 12 disease resistance/host-pathogen interaction-related protein families. Among these, protein kinases, leucine-rich repeats, and zinc finger proteins were of particular interest due to their potential roles in PM resistance. These identified loci can serve as targets for breeding programs aimed at developing disease-resistant wheat cultivars.

Genome-Wide Association Study of Body Conformation Traits in a Three-Way Crossbred Commercial Pig Population

Body conformation is the most direct production index, which can fully reflect pig growth status and is closely related to critical economic traits. In this study, we conducted a genome-wide association study (GWAS) on body conformation traits in a population of 1518 Duroc × (Landrace × Yorkshire) commercial pigs. These traits included body length (BL), body height (BH), chest circumference (CC), abdominal circumference (AC), and waist circumference (WC). Both the mixed linear model (MLM) and fixed and random model circulating probability unification (FarmCPU) approaches were employed for the analysis. Our findings revealed 60 significant single nucleotide polymorphisms (SNPs) associated with these body conformation traits in the crossbred pig population. Specifically, sixteen SNPs were significantly associated with BL, three SNPs with BH, thirteen SNPs with CC, twelve SNPs with AC, and sixteen SNPs with WC. Moreover, we identified several promising candidate genes located within the genomic regions associated with body conformation traits. These candidate genes include INTS10, KIRREL3, SOX21, BMP2, MAP4K3, SOD3, FAM160B1, ATL2, SPRED2, SEC16B, and RASAL2. Furthermore, our analysis revealed a novel significant quantitative trait locus (QTL) on SSC7 specifically associated with waist circumference, spanning an 84 kb interval. Overall, the identification of these significant SNPs and potential candidate genes in crossbred commercial pigs enhances our understanding of the genetic basis underlying body conformation traits. Additionally, these findings provide valuable genetic resources for pig breeding programs.

Application of SVR-Mediated GWAS for Identification of Durable Genetic Regions Associated with Soybean Seed Quality Traits

Soybean (Glycine max L.) is an important food-grade strategic crop worldwide because of its high seed protein and oil contents. Due to the negative correlation between seed protein and oil percentage, there is a dire need to detect reliable quantitative trait loci (QTL) underlying these traits in order to be used in marker-assisted selection (MAS) programs. Genome-wide association study (GWAS) is one of the most common genetic approaches that is regularly used for detecting QTL associated with quantitative traits. However, the current approaches are mainly focused on estimating the main effects of QTL, and, therefore, a substantial statistical improvement in GWAS is required to detect associated QTL considering their interactions with other QTL as well. This study aimed to compare the support vector regression (SVR) algorithm as a common machine learning method to fixed and random model circulating probability unification (FarmCPU), a common conventional GWAS method in detecting relevant QTL associated with soybean seed quality traits such as protein, oil, and 100-seed weight using 227 soybean genotypes. The results showed a significant negative correlation between soybean seed protein and oil concentrations, with heritability values of 0.69 and 0.67, respectively. In addition, SVR-mediated GWAS was able to identify more relevant QTL underlying the target traits than the FarmCPU method. Our findings demonstrate the potential use of machine learning algorithms in GWAS to detect durable QTL associated with soybean seed quality traits suitable for genomic-based breeding approaches. This study provides new insights into improving the accuracy and efficiency of GWAS and highlights the significance of using advanced computational methods in crop breeding research.

Identification of Novel QTLs/Defense Genes in Spring Wheat Germplasm Panel for Seedling and Adult Plant Resistance to Stem Rust and Their Validation Through KASP Marker Assays

Stem rust is one of the major diseases threatening wheat production globally. To identify novel resistance quantitative trait loci (QTLs), we performed 35K Axiom Array SNP genotyping assays on an association mapping panel of 400 germplasm accessions, including Indian landraces, in conjunction with phenotyping for stem rust at seedling and adult plant stages. Association analyses using three genome wide association study (GWAS) models (CMLM, MLMM, and FarmCPU) revealed 20 reliable QTLs for seedling and adult plant resistance. Among these 20 QTLs, five QTLs were found consistent with three models, i.e., four QTLs on chromosome 2AL, 2BL, 2DL, and 3BL for seedling resistance and one QTL on chromosome 7DS for adult plant resistance. Further, we identified a total of 21 potential candidate genes underlying QTLs using gene ontology analysis, including a leucine rich repeat receptor (LRR) and P-loop nucleoside triphosphate hydrolase, which have a role in pathogen recognition and disease resistance. Furthermore, four QTLs (Qsr.nbpgr-3B_11, QSr.nbpgr-6AS_11, QSr.nbpgr-2AL_117-6, and QSr.nbpgr-7BS_APR) were validated through KASP located on chromosomes 3B, 6A, 2A, and 7B. Out of these QTLs, QSr.nbpgr-7BS_APR was identified as a novel QTL for stem rust resistance which has been found effective in both seedling as well as the adult plant stages. Identified novel genomic regions and validated QTLs have the potential to be deployed in wheat improvement programs to develop disease resistant varieties for stem rust and can diversify the genetic basis of resistance.

Genome-wide association study to identify genomic loci associated with early vigor in bread wheat under simulated water deficit complemented with quantitative trait loci meta-analysis

A genome-wide association study (GWAS) was used to identify associated loci with early vigor under simulated water deficit and grain yield under field drought in a diverse collection of Iranian bread wheat landraces. In addition, a meta-quantitative trait loci (MQTL) analysis was used to further expand our approach by retrieving already published quantitative trait loci (QTL) from recombinant inbred lines, double haploids, back-crosses, and F2 mapping populations. In the current study, around 16%, 14%, and 16% of SNPs were in significant linkage disequilibrium (LD) in the A, B, and D genomes, respectively, and varied between 5.44% (4A) and 21.85% (6A). Three main subgroups were identified among the landraces with different degrees of admixture, and population structure was further explored through principal component analysis. Our GWAS identified 54 marker-trait associations (MTAs) that were located across the wheat genome but with the highest number found in the B sub-genome. The gene ontology (GO) analysis of MTAs revealed that around 75% were located within or closed to protein-coding genes. In the MQTL analysis, 23 MQTLs, from a total of 215 QTLs, were identified and successfully projected onto the reference map. MQT-YLD4, MQT-YLD9, MQT-YLD13, MQT-YLD17, MQT-YLD18, MQT-YLD19, and MQTL-RL1 contributed to the highest number of projected QTLs and were therefore regarded as the most reliable and stable QTLs under water deficit conditions. These MQTLs greatly facilitate the identification of putative candidate genes underlying at each MQTL interval due to the reduced confidence of intervals associated with MQTLs. These findings provide important information on the genetic basis of early vigor traits and grain yield under water deficit conditions and set the foundation for future investigations into adaptation to water deficit in bread wheat.

Association mapping unravels the genetics controlling seedling drought stress tolerance in winter wheat

Drought is a major constraint in wheat (Triticum aestivum L.) grain yield. The present work aimed to identify quantitative trait nucleotides (QTNs)/ candidate genes influencing drought tolerance-related traits at the seedling stage in 261 accessions of a diverse winter wheat panel. Seeds from three consecutive years were exposed to polyethylene glycol 12% (PEG-6000) and a control treatment (distilled water). The Farm-CPU method was used for the association analysis with 17,093 polymorphic SNPs. PEG treatment reduced shoot length (SL) (-36.3%) and root length (RL) (-11.3%) compared with control treatments, while the coleoptile length (CL) was increased by 11% under drought conditions, suggesting that it might be considered as an indicator of stress-tolerance. Interestingly, we revealed 70 stable QTN across 17 chromosomes. Eight QTNs related to more than one trait were detected on chromosomes 1B, 2A (2), 2B, 2D, 4B, 7A, and 7B and located nearby or inside candidate genes within the linkage disequilibrium (LD) interval. For instance, the QTN on chromosome 2D is located inside the gene TraesCS2D02G133900 that controls the variation of CL_S and SL_C. The allelic variation at the candidate genes showed significant influence on the associated traits, demonstrating their role in controlling the natural variation of multi-traits of drought stress tolerance. The gene expression of these candidate genes under different stress conditions validates their biological role in stress tolerance. Our findings offer insight into understanding the genetic factors and diverse mechanisms in response to water shortage conditions that are important for wheat improvement and adaptation at early developmental stages.

Genome-Wide Association Study (GWAS) of White Mold Resistance in Snap Bean

White mold can result in snap bean yield losses of 90 to 100% when field conditions favor the pathogen. A genome-wide association study (GWAS) was conducted to detect loci significantly associated with white mold resistance in a panel of snap bean (Phaseolus vulgaris L.) cultivars. Two populations of snap bean were used in this study. The first population was the BeanCAP (Coordinated Agriculture Project) Snap Bean Diversity Panel (SBDP) (n = 136), and the second population was the Snap Bean Association Panel (SnAP) (n = 378). SBDP was evaluated for white mold reaction in the field in 2012 and 2013, and SnAP was screened in a greenhouse only using the seedling straw test in 2016. Two reference genomes representing the Andean and Middle American centers of domestication were utilized to align the genotyping-by-sequencing (GBS) data. A GWAS was performed using FarmCPU with one principal component after comparing five models. Thirty-four single-nucleotide polymorphisms (SNPs) significantly associated with white mold resistance were detected. Eleven significant SNPs were identified by the seedling straw test, and 23 significant SNPs were identified by field data. Fifteen SNPs were identified within a 100 kb window containing pentatricopeptide repeat (PPR)-encoding genes, and eleven were close to leucine-rich repeat (LRR)-encoding genes, suggesting that these two classes are of outsized importance for snap bean resistance to white mold.

Genome-Wide Association Study for Udder Conformation Traits in Chinese Holstein Cattle

Udder conformation traits are one of the most economic traits in dairy cows, greatly affecting animal health, milk production, and producer profitability in the dairy industry. Genetic analysis of udder structure and scores have been developed in Holstein cattle. In our research, we conducted a genome-wide association study for five udder traits, including anterior udder attachment (AUA), central suspensory ligament (CSL), posterior udder attachment height (PUAH), posterior udder attachment width (PUAW), and udder depth (UD), in which the fixed and random model circulating probability unification (FarmCPU) model was applied for the association analysis. The heritability and the standard errors of these five udder traits ranged from 0.04 ± 0.00 to 0.49 ± 0.03. Phenotype data were measured from 1000 Holstein cows, and the GeneSeek Genomic Profiler (GGP) Bovine 100 K SNP chip was used to analyze genotypic data in Holstein cattle. For GWAS analysis, 984 individual cows and 84,407 single-nucleotide polymorphisms (SNPs) remained after quality control; a total of 18 SNPs were found at the GW significant threshold (p < 5.90 × 10−7). Many candidate genes were identified within 200kb upstream or downstream of the significant SNPs, which include MGST1, MGST2, MTUS1, PRKN, STXBP6, GRID2, E2F8, CDH11, FOXP1, SLF1, TMEM117, SBF2, GC, ADGRB3, and GCLC. Pathway analysis revealed that 58 Gene Ontology (GO) terms and 18 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were enriched with adjusted p values, and these GO terms and the KEGG pathway analysis were associated with biological information, metabolism, hormonal growth, and development processes. These results could give valuable biological information for the genetic architecture of udder conformation traits in dairy Holstein cattle.

Genome-Wide Association Analysis Reveals Trait-Linked Markers for Grain Nutrient and Agronomic Traits in Diverse Set of Chickpea Germplasm

Chickpea is an inexpensive source of protein, minerals, and vitamins to the poor people living in arid and semi-arid regions of Southern Asia and Sub-Saharan Africa. New chickpea cultivars with enhanced levels of protein, Fe and Zn content are a medium-term strategy for supplying essential nutrients for human health and reducing malnutrition. In the current study, a chickpea reference set of 280 accessions, including landraces, breeding lines, and advanced cultivars, was evaluated for grain protein, Fe, Zn content and agronomic traits over two seasons. Using a mid-density 5k SNP array, 4603 highly informative SNPs distributed across the chickpea genome were used for GWAS analysis. Population structure analysis revealed three subpopulations (K = 3). Linkage disequilibrium (LD) was extensive, and LD decay was relatively low. A total of 20 and 46 marker-trait associations (MTAs) were identified for grain nutrient and agronomic traits, respectively, using FarmCPU and BLINK models. Of which seven SNPs for grain protein, twelve for Fe, and one for Zn content were distributed on chromosomes 1, 4, 6, and 7. The marker S4_4477846 on chr4 was found to be co-associated with grain protein over seasons. The markers S1_11613376 and S1_2772537 co-associated with grain Fe content under NSII and pooled seasons and S7_9379786 marker under NSI and pooled seasons. The markers S4_31996956 co-associated with grain Fe and days to maturity. SNP annotation of associated markers were found to be related to gene functions of metal ion binding, transporters, protein kinases, transcription factors, and many more functions involved in plant metabolism along with Fe and protein homeostasis. The identified significant MTAs has potential use in marker-assisted selection for developing nutrient-rich chickpea cultivars after validation in the breeding populations.

Machine-Learning-Based Genome-Wide Association Studies for Uncovering QTL Underlying Soybean Yield and Its Components

A genome-wide association study (GWAS) is currently one of the most recommended approaches for discovering marker-trait associations (MTAs) for complex traits in plant species. Insufficient statistical power is a limiting factor, especially in narrow genetic basis species, that conventional GWAS methods are suffering from. Using sophisticated mathematical methods such as machine learning (ML) algorithms may address this issue and advance the implication of this valuable genetic method in applied plant-breeding programs. In this study, we evaluated the potential use of two ML algorithms, support-vector machine (SVR) and random forest (RF), in a GWAS and compared them with two conventional methods of mixed linear models (MLM) and fixed and random model circulating probability unification (FarmCPU), for identifying MTAs for soybean-yield components. In this study, important soybean-yield component traits, including the number of reproductive nodes (RNP), non-reproductive nodes (NRNP), total nodes (NP), and total pods (PP) per plant along with yield and maturity, were assessed using a panel of 227 soybean genotypes evaluated at two locations over two years (four environments). Using the SVR-mediated GWAS method, we were able to discover MTAs colocalized with previously reported quantitative trait loci (QTL) with potential causal effects on the target traits, supported by the functional annotation of candidate gene analyses. This study demonstrated the potential benefit of using sophisticated mathematical approaches, such as SVR, in a GWAS to complement conventional GWAS methods for identifying MTAs that can improve the efficiency of genomic-based soybean-breeding programs.

Corrigendum: Multi-Locus Genome-Wide Association Studies Reveal Fruit Quality Hotspots in Peach Genome

[This corrects the article DOI: 10.3389/fpls.2021.644799.].

Performing Genome-Wide Association Studies Using rMVP

Genome wide association study (GWAS), which is a powerful tool to detect the relationship between the traits of interest and high-density markers, has provided unprecedented insights into the genetic basis of quantitative variation for complex traits. Along with the development of high-throughput sequencing technology, both sample sizes and marker sizes are increasing rapidly, which make computations more challenging than ever. Therefore, to efficiently process big data with limited computing resources in a reasonable time and to use state-of-the-art statistical models to reduce false positive and false negative rates have always been hot topics in the domain of GWAS. In this chapter, we describe how to perform GWAS using an R package, rMVP, which includes data preparation, evaluation of population structure, association tests by different models, and high-quality visualization of GWAS results.

Pyramiding of scald resistance genes in four spring barley MAGIC populations

Genome-Wide Association Studies (GWAS) of four Multi-parent Advanced Generation Inter-Cross (MAGIC) populations identified nine regions on chromosomes 1H, 3H, 4H, 5H, 6H and 7H associated with resistance against barley scald disease. Three of these regions are putatively novel resistance Quantitative Trait Loci (QTL). Barley scald is caused by Rhynchosporium commune, one of the most important barley leaf diseases that are prevalent in most barley-growing regions. Up to 40% yield losses can occur in susceptible barley cultivars. Four MAGIC populations were generated in a Nordic Public-Private Pre-breeding of spring barley project (PPP Barley) to introduce resistance to several important diseases. Here, these MAGIC populations consisting of six to eight founders each were tested for scald resistance in field trials in Finland and Iceland. Eight different model covariate combinations were compared for GWAS studies, and the models that deviated the least from the expected p-values were selected. For all QTL, candidate genes were identified that are predicted to be involved in pathogen defence. The MAGIC progenies contained new haplotypes of significant SNP-markers with high resistance levels. The lines with successfully pyramided resistance against scald and mildew and the significant markers are now distributed among Nordic plant breeders and will benefit development of disease-resistant cultivars.

Whole-genome re-sequencing association study on yearling wool traits in Chinese fine-wool sheep

To investigate single nucleotide polymorphism (SNP) loci associated with yearling wool traits of fine-wool sheep for optimizing marker-assisted selection and dissection of the genetic architecture of wool traits, we conducted a genome-wide association study (GWAS) based on the fixed and random model circulating probability unification (FarmCPU) for yearling staple length (YSL), yearling mean fiber diameter (YFD), yearling greasy fleece weight (YGFW), and yearling clean fleece rate (YCFR) by using the whole-genome re-sequenced data (totaling 577 sheep) from the following four fine-wool sheep breeds in China: Alpine Merino sheep (AMS), Chinese Merino sheep (CMS), Qinghai fine-wool sheep (QHS), and Aohan fine-wool sheep (AHS). A total of 16 SNPs were detected above the genome-wise significant threshold (P = 5.45E-09), and 79 SNPs were located above the suggestive significance threshold (P = 5.00E-07) from the GWAS results. For YFD and YGFW traits, 7 and 9 SNPs reached the genome-wise significance thresholds, whereas 10 and 12 SNPs reached the suggestive significance threshold, respectively. For YSL and YCFR traits, none of the SNPs reached the genome-wise significance thresholds, whereas 57 SNPs exceeded the suggestive significance threshold. We recorded 14 genes located at the region of ±50-kb near the genome-wise significant SNPs and 59 genes located at the region of ±50-kb near the suggestive significant SNPs. Meanwhile, we used the Average Information Restricted Maximum likelihood algorithm (AI-REML) in the “HIBLUP” package to estimate the heritability and variance components of the four desired yearling wool traits. The estimated heritability values (h²) of YSL, YFD, YGFW, and YCFR were 0.6208, 0.7460, 0.6758, and 0.5559, respectively. We noted that the genetic parameters in this study can be used for fine-wool sheep breeding. The newly detected significant SNPs and the newly identified candidate genes in this study would enhance our understanding of yearling wool formation, and significant SNPs can be applied to genome selection in fine-wool sheep breeding.

Genome-Wide Association Study Identifies Candidate Genes Associated with Feet and Leg Conformation Traits in Chinese Holstein Cattle

Simple Summary

Feet and leg problems are among the major reasons for dairy cows leaving the herd, as well as having direct association with production and reproduction efficiency, health (e.g., claw disorders and lameness) and welfare. Hence, understanding the genetic architecture underlying feet and conformation traits in dairy cattle offers new opportunities toward the genetic improvement and long-term selection. Through a genome-wide association study on Chinese Holstein cattle, we identified several candidate genes associated with feet and leg conformation traits. These results could provide useful information about the molecular breeding basis of feet and leg traits, thus improving the longevity and productivity of dairy cattle.

Abstract

Feet and leg conformation traits are considered one of the most important economical traits in dairy cattle and have a great impact on the profitability of milk production. Therefore, identifying the single nucleotide polymorphisms (SNPs), genes and pathways analysis associated with these traits might contribute to the genomic selection and long-term plan selection for dairy cattle. We conducted genome-wide association studies (GWASs) using the fixed and random model circulating probability unification (FarmCPU) method to identify SNPs associated with bone quality, heel depth, rear leg side view and rear leg rear view of Chinese Holstein cows. Phenotypic measurements were collected from 1000 individuals of Chinese Holstein cattle and the GeneSeek Genomic Profiler Bovine 100 K SNP chip was utilized for individual genotyping. After quality control, 984 individual cows and 84,906 SNPs remained for GWAS work; as a result, we identified 20 significant SNPs after Bonferroni correction. Several candidate genes were identified within distances of 200 kb upstream or downstream to the significant SNPs, including ADIPOR2, INPP4A, DNMT3A, ALDH1A2, PCDH7, XKR4 and CADPS. Further bioinformatics analyses showed 34 gene ontology terms and two signaling pathways were significantly enriched (p ≤ 0.05). Many terms and pathways are related to biological quality, metabolism and development processes; these identified SNPs and genes could provide useful information about the genetic architecture of feet and leg traits, thus improving the longevity and productivity of Chinese Holstein dairy cattle.

Multi-Locus Genome-Wide Association Studies Reveal Fruit Quality Hotspots in Peach Genome

Peach is one of the most important fruit crops in the world, with the global annual production about 24.6 million tons. The United States is the fourth-largest producer after China, Spain, and Italy. Peach consumption has decreased over the last decade, most likely due to inconsistent quality of the fruit on the market. Thus, marker-assisted selection for fruit quality traits is highly desired in fresh market peach breeding programs and one of the major goals of the RosBREED project. The ability to use DNA information to select for desirable traits would enable peach breeders to efficiently plan crosses and select seedlings with desired quality traits early in the selection process before fruiting. Therefore, we assembled a multi-locus genome wide association study (GWAS) of 620 individuals from three public fresh market peach breeding programs (Arkansas, Texas, and South Carolina). The material was genotyped using 9K SNP array and the traits were phenotyped for three phenological (bloom date, ripening date, and days after bloom) and 11 fruit quality-related traits (blush, fruit diameter, fruit weight, adherence, fruit firmness, redness around pit, fruit texture, pit weight, soluble solid concentration, titratable acidity, and pH) over three seasons (2010, 2011, and 2012). Multi-locus association analyses, carried out using mrMLM 4.0 and FarmCPU R packages, revealed a total of 967 and 180 quantitative trait nucleotides (QTNs), respectively. Among the 88 consistently reliable QTNs detected using multiple multi-locus GWAS methods and/or at least two seasons, 44 were detected for the first time. Fruit quality hotspots were identified on chromosomes 1, 3, 4, 5, 6, and 8. Out of 566 candidate genes detected in the genomic regions harboring the QTN clusters, 435 were functionally annotated. Gene enrichment analyses revealed 68 different gene ontology (GO) terms associated with fruit quality traits. Data reported here advance our understanding of genetic mechanisms underlying important fruit quality traits and further support the development of DNA tools for breeding.

GWAS-Based Identification of New Loci for Milk Yield, Fat, and Protein in Holstein Cattle

Simple Summary

Understanding the genetic architecture underlying milk production traits in cattle is beneficial so that genetic variants can be targeted toward the genetic improvement. In this study, we performed a genome-wide association study for milk production and quality traits in Holstein cattle. In the total of ten significant single-nucleotide polymorphisms (SNPs) associated with milk fat and protein, six are located in previously reported quantitative traits locus (QTL) regions. The study not only identified the effect of DGAT1 gene on milk fat and protein but also found several novel candidate genes. In addition, some pleiotropic SNPs and QTLs were identified that associated with more than two traits, these results could provide some basis for molecular breeding in dairy cattle.

Abstract

High-yield and high-quality of milk are the primary goals of dairy production. Understanding the genetic architecture underlying these milk-related traits is beneficial so that genetic variants can be targeted toward the genetic improvement. In this study, we measured five milk production and quality traits in Holstein cattle population from China. These traits included milk yield, fat, and protein. We used the estimated breeding values as dependent variables to conduct the genome-wide association studies (GWAS). Breeding values were estimated through pedigree relationships by using a linear mixed model. Genotyping was carried out on the individuals with phenotypes by using the Illumina BovineSNP150 BeadChip. The association analyses were conducted by using the fixed and random model Circulating Probability Unification (FarmCPU) method. A total of ten single-nucleotide polymorphisms (SNPs) were detected above the genome-wide significant threshold (p < 4.0 × 10⁻⁷), including six located in previously reported quantitative traits locus (QTL) regions. We found eight candidate genes within distances of 120 kb upstream or downstream to the associated SNPs. The study not only identified the effect of DGAT1 gene on milk fat and protein, but also discovered novel genetic loci and candidate genes related to milk traits. These novel genetic loci would be an important basis for molecular breeding in dairy cattle.

Genome-wide association studies for the concentrations of insulin, triiodothyronine, and thyroxine in Chinese Holstein cattle

To further understand the genetic structure that is associated with insulin (INS) and thyroid hormones (TH), including triiodothyronine (T₃) and thyroxine (T₄), in Chinese Holstein cows, we conducted a genome-wide association study (GWAS) of thyroid hormones and insulin in cows. We conducted GWAS analysis on 1217 Chinese Holstein cows raised in southern China and found 19 significant single nucleotide polymorphisms (SNPs) in this study: 10 SNPs were associated with INS, 5 SNPs were associated with T₃, and 4 SNPs were associated with T₄. In our study, the GWAS method was used for preliminary screening on related genes of traits, and due to insufficient relevant literature, a functional analysis of genes could only be based on human studies. We observed that DGKB from Bos taurus chromosome (BTA)4 is strongly associated with insulin secretion. We found that EXOC4 gene was significantly correlated with T₃ and T₄ traits. Another significant SNP was located in the CYP7A1 gene, which has been confirmed to be affected by thyroid hormones.

Last-Generation Genome-Environment Associations Reveal the Genetic Basis of Heat Tolerance in Common Bean (Phaseolus vulgaris L.)

Genome-environment associations (GEAs) are a powerful strategy for the study of adaptive traits in wild plant populations, yet they still lack behind in the use of modern statistical methods as the ones suggested for genome-wide association studies (GWASs). In order to bridge this gap, we couple GEA with last-generation GWAS algorithms in common bean to identify novel sources of heat tolerance across naturally heterogeneous ecosystems. Common bean (Phaseolus vulgaris L.) is the most important legume for human consumption, and breeding it for resistance to heat stress is key because annual increases in atmospheric temperature are causing decreases in yield of up to 9% for every 1°C. A total of 78 geo-referenced wild accessions, spanning the two gene pools of common bean, were genotyped by sequencing (GBS), leading to the discovery of 23,373 single-nucleotide polymorphism (SNP) markers. Three indices of heat stress were developed for each accession and inputted in last-generation algorithms (i.e. SUPER, FarmCPU, and BLINK) to identify putative associated loci with the environmental heterogeneity in heat stress. Best-fit models revealed 120 significantly associated alleles distributed in all 11 common bean chromosomes. Flanking candidate genes were identified using 1-kb genomic windows centered in each associated SNP marker. Some of these genes were directly linked to heat-responsive pathways, such as the activation of heat shock proteins (MED23, MED25, HSFB1, HSP40, and HSP20). We also found protein domains related to thermostability in plants such as S1 and Zinc finger A20 and AN1. Other genes were related to biological processes that may correlate with plant tolerance to high temperature, such as time to flowering (MED25, MBD9, and PAP), germination and seedling development (Pkinase_Tyr, Ankyrin-B, and Family Glicosil-hydrolase), cell wall stability (GAE6), and signaling pathway of abiotic stress via abscisic acid (histone-like transcription factors NFYB and phospholipase C) and auxin (Auxin response factor and AUX_IAA). This work offers putative associated loci for marker-assisted and genomic selection for heat tolerance in common bean. It also demonstrates that it is feasible to identify genome-wide environmental associations with modest sample sizes by using a combination of various carefully chosen environmental indices and last-generation GWAS algorithms.

Single-plant GWAS coupled with bulk segregant analysis allows rapid identification and corroboration of plant-height candidate SNPs

BACKGROUND

Genome wide association studies (GWAS) are a powerful tool for identifying quantitative trait loci (QTL) and causal single nucleotide polymorphisms (SNPs)/genes associated with various important traits in crop species. Typically, GWAS in crops are performed using a panel of inbred lines, where multiple replicates of the same inbred are measured and the average phenotype is taken as the response variable. Here we describe and evaluate single plant GWAS (sp-GWAS) for performing a GWAS on individual plants, which does not require an association panel of inbreds. Instead sp-GWAS relies on the phenotypes and genotypes from individual plants sampled from a randomly mating population. Importantly, we demonstrate how sp-GWAS can be efficiently combined with a bulk segregant analysis (BSA) experiment to rapidly corroborate evidence for significant SNPs.

RESULTS

In this study we used the Shoepeg maize landrace, collected as an open pollinating variety from a farm in Southern Missouri in the 1960's, to evaluate whether sp-GWAS coupled with BSA can efficiently and powerfully used to detect significant association of SNPs for plant height (PH). Plant were grown in 8 locations across two years and in total 768 individuals were genotyped and phenotyped for sp-GWAS. A total of 306 k polymorphic markers in 768 individuals evaluated via association analysis detected 25 significant SNPs (P ≤ 0.00001) for PH. The results from our single-plant GWAS were further validated by bulk segregant analysis (BSA) for PH. BSA sequencing was performed on the same population by selecting tall and short plants as separate bulks. This approach identified 37 genomic regions for plant height. Of the 25 significant SNPs from GWAS, the three most significant SNPs co-localize with regions identified by BSA.

CONCLUSION

Overall, this study demonstrates that sp-GWAS coupled with BSA can be a useful tool for detecting significant SNPs and identifying candidate genes. This result is particularly useful for species/populations where association panels are not readily available.

Genotyping of Soybean Cultivars With Medium-Density Array Reveals the Population Structure and QTNs Underlying Maturity and Seed Traits

Soybean was domesticated about 5,000 to 6,000 years ago in China. Although genotyping technologies such as genotyping by sequencing (GBS) and high-density array are available, it is convenient and economical to genotype cultivars or populations using medium-density SNP array in genetic study as well as in molecular breeding. In this study, 235 cultivars, collected from China, Japan, USA, Canada and some other countries, were genotyped using SoySNP8k iSelect BeadChip with 7,189 single nucleotide polymorphisms (SNPs). In total, 4,471 polymorphic SNP markers were used to analyze population structure and perform genome-wide association study (GWAS). The most likely K value was 7, indicating this population can be divided into 7 subpopulations, which is well in accordance with the geographic origins of cultivars or accession studied. The LD decay rate was estimated at 184 kb, where r2 dropped to half of its maximum value (0.205). GWAS using FarmCPU detected a stable quantitative trait nucleotide (QTN) for hilum color and seed color, which is consistent with the known loci or genes. Although no universal QTNs for flowering time and maturity were identified across all environments, a total of 30 consistent QTNs were detected for flowering time (R1) or maturity (R7 and R8) on 16 chromosomes, most of them were corresponding to known E1 to E4 genes or QTL region reported in SoyBase (soybase.org). Of 16 consistent QTNs for protein and oil contents, 11 QTNs were detected having antagonistic effects on protein and oil content, while 4 QTNs soly for oil content, and one QTN soly for protein content. The information gained in this study demonstrated that the usefulness of the medium-density SNP array in genotyping for genetic study and molecular breeding.

Genome-wide association mapping of resistance to a Brazilian isolate of Sclerotinia sclerotiorum in soybean genotypes mostly from Brazil

BACKGROUND

Sclerotinia Stem Rot (SSR), caused by the fungal pathogen Sclerotinia sclerotiorum, is ubiquitous in cooler climates where soybean crops are grown. Breeding for resistance to SSR remains challenging in crops like soybean, where no single gene provides strong resistance, but instead, multiple genes work together to provide partial resistance. In this study, a genome-wide association study (GWAS) was performed to dissect the complex genetic architecture of soybean quantitative resistance to SSR and to provide effective molecular markers that could be used in breeding programs. A collection of 420 soybean genotypes were selected based on either reports of resistance, or from one of three different breeding programs in Brazil, two commercial, one public. Plant genotype sensitivity to SSR was evaluated by the cut stem inoculation method, and lesion lengths were measured at 4 days post inoculation.

RESULTS

Genotyping-by-sequencing was conducted to genotype the 420 soybean lines. The TASSEL 5 GBSv2 pipeline was used to call SNPs under optimized parameters, and with the extra step of trimming adapter sequences. After filtering missing data, heterozygosity, and minor allele frequency, a total of 11,811 SNPs and 275 soybean genotypes were obtained for association analyses. Using a threshold of FDR-adjusted p-values <0.1, the Compressed Mixed Linear Model (CMLM) with Genome Association and Prediction Integrated Tool (GAPIT), and the Fixed and Random Model Circulating Probability Unification (FarmCPU) methods, both approaches identified SNPs with significant association to disease response on chromosomes 1, 11, and 18. The CMLM also found significance on chromosome 19, whereas FarmCPU also identified significance on chromosomes 4, 9, and 16.

CONCLUSIONS

These similar and yet different results show that the computational methods used can impact SNP associations in soybean, a plant with a high degree of linkage disequilibrium, and in SSR resistance, a trait that has a complex genetic basis. A total of 125 genes were located within linkage disequilibrium of the three loci shared between the two models. Their annotations and gene expressions in previous studies of soybean infected with S. sclerotiorum were examined to narrow down the candidates.