Single-locus Genome-wide Association Studies Research Articles

Unveiling the predominance of Saccharum spontaneum alleles for resistance to orange rust in sugarcane using genome-wide association.

Six QTLs of resistance to sugarcane orange rust were identified in modern interspecific hybrids by GWAS. For five of them, the resistance alleles originated from S. spontaneum. Altogether, they efficiently predict disease resistance. Sugarcane orange rust (SOR) is a threatening emerging disease in many sugarcane industries worldwide. Improving the genetic resistance of commercial cultivars remains the most promising solution to control this disease. In this study, an association panel of 568 modern interspecific sugarcane hybrids (Saccharum officinarum x S. spontaneum) from Réunion's breeding program was evaluated for its resistance to SOR under natural conditions of infection. Two genome-wide association studies (GWAS) were conducted between disease reactions and 183,842 single nucleotide polymorphism (SNP) markers obtained by targeted genotyping-by-sequencing. Five resistance quantitative trait loci (QTLs), named Oru1, Oru2, Oru3, Oru4 and Oru5, were identified using a single-locus GWAS (SL-GWAS). These five QTLs all originated from the species S. spontaneum. A multi-locus GWAS (ML-GWAS) uncovered an additional but less significant resistance QTL named Oru6, which originated from S. officinarum. All six QTLs had a moderate to major phenotypic effect on disease resistance. Prediction accuracy estimated with linear regression models based on each of the five QTLs identified by SL-GWAS was between 0.16-0.41. Altogether, these five QTLs provided a relatively high prediction accuracy of 0.60. In comparison, accuracies obtained with six genome-wide prediction models (i.e., GBLUP, Bayes-A, Bayes-B, Bayes-C, Bayesian Lasso and RKHS) reached only 0.65. The good prediction accuracy of disease resistance provided by the QTLs and the predominant S. spontaneum origin of their resistance alleles pave the way for effective marker-assisted breeding strategies.

Single-Locus and Multi-Locus Genome-Wide Association Studies Identify Genes Associated with Liver Cu Concentration in Merinoland Sheep.

Economic losses due to copper intoxication or deficiency is a problem encountered by sheep farmers. The aim of this study was to investigate the ovine genome for genomic regions and candidate genes responsible for variability in liver copper concentration. Liver samples were collected from slaughtered lambs of the Merinoland breed from two farms, and used for measurement of copper concentration and genome-wide association study (GWAS). A total of 45,511 SNPs and 130 samples were finally used for analysis, in which single-locus and several multi-locus GWAS (SL-GWAS; ML-GWAS) methods were employed. Gene enrichment analysis was performed for identified candidate genes to detect gene ontology (GO) terms significantly associated with hepatic copper levels. The SL-GWAS and a minimum of two ML-GWAS identified two and thirteen significant SNPs, respectively. Within genomic regions surrounding identified SNPs, we observed nine promising candidate genes such as DYNC1I2, VPS35, SLC38A9 and CHMP1A. GO terms such as lysosomal membrane, mitochondrial inner membrane and sodium:proton antiporter activity were significantly enriched. Genes involved in these identified GO terms mediate multivesicular body (MVB) fusion with lysosome for degradation and control mitochondrial membrane permeability. This reveals the polygenic status of this trait and candidate genes for further studies on breeding for copper tolerance in sheep.

Identification of candidate genes associated with carcass component weights in commercial crossbred pigs through a combined GWAS approach.

In the pork industry chain, carcass cutting is crucial for enhancing the commercial value of pork carcasses. However, the genetic mechanisms underlying carcass component weights remain poorly understood. Here, we used a combined genome-wide association study (GWAS) approach that integrated single- and multi-locus models to map genetic markers and genes associated with the weights of seven carcass components in Duroc × Landrace × Yorkshire (DLY) pigs. As multi-locus GWAS captures more single nucleotide polymorphisms (SNPs) with large effects than single-locus GWAS, the combined GWAS approach detected more SNPs than using the single-locus model alone. We identified 177 nonredundant SNPs associated with these traits in 526 DLY pigs, including boneless butt shoulder (BBS), boneless picnic shoulder (BPS), boneless leg (BL), belly (BELLY), front fat (FF), rear fat (RF), and skin-on whole loin (SLOIN). Using single-locus GWAS, we identified a quantitative trait locus (QTL) for SLOIN on Sus scrofa chromosome 15 (SSC15). Notably, a single SNP (ASGA0069883) in the proximity of this QTL was consistently detected by all GWAS models (one single-locus and four multi-locus models) and explained more than 4% of the phenotypic variance. Our findings suggest that the involved gene, MYO3B, is proposed to be a strong candidate for SLOIN. Further analysis also identified several candidate genes related to BBS (PPP3CA and CPEB4), BPS (ECH1), FF (CACNB2 and ZNF217), BELLY (FGFRL1), BL (CHST11), and RF (LRRK2). The identified SNPs can be used as molecular markers for the genetic improvement of pork carcasses in the molecular-guided breeding of modern commercial pigs.

Identification of Genomic Regions and Candidate Genes for Litter Traits in French Large White Pigs Using Genome-Wide Association Studies.

Simple SummarySow reproductive performance is crucial in the entire pig industry. However, the heritability of reproductive traits is slowthat it is difficult to make great genetic progress by traditional breeding methods. Therefore, in this study, we performed genome-wide association analyses of litter traits in French White sows to explore genes related to sow reproductive traits and we finally detected a notable region (10.72–11.06 Mb on chromosome 7) and four promising genes (JARID2, PDIA6, FLRT2 and DICER1). The results of this study provide a theoretical basis for subsequent molecular breeding.The reproductive traits of sows are one of the important economic traits in pig production, and their performance directly affects the economic benefits of the entire pig industry. In this study, a total of 895 French Large White pigs were genotyped by GeneSeek Porcine 50K SNP Beadchip and four phenotypic traits of 1407 pigs were recorded, including total number born (TNB), number born alive (NBA), number healthy piglets (NHP) and litter weight born alive (LWB). To identify genomic regions and genes for these traits, we used two approaches: a single-locus genome-wide association study (GWAS) and a single-step GWAS (ssGWAS). Overall, a total of five SNPs and 36 genomic regions were identified by single-locus GWAS and ssGWAS, respectively. Notably, fourof all five significant SNPs were located in 10.72–11.06 Mb on chromosome 7, were also identified by ssGWAS. These regions explained the highest or second highest genetic variance in the TNB, NBA and NHP traits and harbor the protein coding gene ENSSSCG00000042180. In addition, several candidate genes associated with litter traits were identified, including JARID2, PDIA6, FLRT2 and DICER1. Overall, these novel results reflect the polygenic genetic architecture of the litter traits and provide a theoretical reference for the following implementation of molecular breeding.

Identifying loci controlling total starch content of leaf in Nicotiana tabacum through genome-wide association study.

Starch is an important primary metabolite in plants, which can provide bioenergy for fuel ethanol production. There are many studies focusing on starch metabolism in Arabidopsis, maize, and rice, but few reports have been made on the starch content of tobacco leaves. Hence, to identify the marker-trait associations and isolate the candidate genes related to starch content of tobacco leaf, the genome-wide association study (GWAS) was performed using a multiparent advanced generation intercross (MAGIC) population consisting of 276 accessions genotyped by a 430K SNP array. In this study, we detected the leaf starch content of tobacco plants cultivated in two places (Zhucheng and Chenzhou), which showed a wide variation of starch content in the population. A total of 28 and 45 significant single-nucleotide polymorphism (SNP) loci associated with leaf starch content were identified by single-locus and multi-locus GWAS models, respectively, and the phenotypic variance explained by these loci varied from 1.80 to - 14.73%. Furthermore, among these quantitative trait loci (QTLs), one SNP, AX-106011713 located on chromosome 19, was detected repeatedly in multiple models and two environments, which was selected for linkage disequilibrium (LD) analysis to obtain the target candidate region. Through gene annotation, haplotype, and gene expression analysis, two candidate genes encoding E3 ubiquitin-protein ligase (Ntab0823160) and fructose-bisphosphate aldolase (Ntab0375050) were obtained. Results showed that the variety carrying the beneficial alleles of the two candidate genes had higher gene expression level and leaf starch content, suggesting the potential role of candidate genes in enhancing the level of tobacco leaf starch content. Furthermore, silencing of Ntab0823160 in tobacco leaves reduced the content of total starch to 39.41-69.75% of that in the wide type plants. Taken together, our results provide useful resources for further investigation of the starch metabolic pathway and are also beneficial for the creation of eco-friendly cultivars with increased accumulation of leaf starch content.

Genome-wide association studies dissect the G\u2009\xd7\u2009E interaction for agronomic traits in a worldwide collection of safflowers (Carthamus tinctorius L.)

Genome-wide association studies were conducted using a globally diverse safflower (Carthamus tinctorius L.) Genebank collection for grain yield (YP), days to flowering (DF), plant height (PH), 500 seed weight (SW), seed oil content (OL), and crude protein content (PR) in four environments (sites) that differed in water availability. Phenotypic variation was observed for all traits. YP exhibited low overall genetic correlations (rGoverall) across sites, while SW and OL had high rGoverall and high pairwise genetic correlations (rGij) across all pairwise sites. In total, 92 marker-trait associations (MTAs) were identified using three methods, single locus genome-wide association studies (GWAS) using a mixed linear model (MLM), the Bayesian multi-locus method (BayesR), and meta-GWAS. MTAs with large effects across all sites were detected for OL, SW, and PR, and MTAs specific for the different water stress sites were identified for all traits. Five MTAs were associated with multiple traits; 4 of 5 MTAs were variously associated with the three traits of SW, OL, and PR. This study provided insights into the phenotypic variability and genetic architecture of important safflower agronomic traits under different environments.

Multi-locus genome-wide association studies for five yield-related traits in rice

BackgroundImproving the overall production of rice with high quality is a major target of breeders. Mining potential yield-related loci have been geared towards developing efficient rice breeding strategies. In this study, one single-locus genome-wide association studies (SL-GWAS) method (MLM) in conjunction with five multi-locus genome-wide association studies (ML-GWAS) approaches (mrMLM, FASTmrMLM, pLARmEB, pKWmEB, and ISIS EM-BLASSO) were conducted in a panel consisting of 529 rice core varieties with 607,201 SNPs.ResultsA total of 152, 106, 12, 111, and 64 SNPs were detected by the MLM model associated with the five yield-related traits, namely grain length (GL), grain width (GW), grain thickness (GT), thousand-grain weight (TGW), and yield per plant (YPP), respectively. Furthermore, 74 significant quantitative trait nucleotides (QTNs) were presented across at least two ML-GWAS methods to be associated with the above five traits successively. Finally, 20 common QTNs were simultaneously discovered by both SL-GWAS and ML-GWAS methods. Based on genome annotation, gene expression analysis, and previous studies, two candidate key genes (LOC_Os09g02830 and LOC_Os07g31450) were characterized to affect GW and TGW, separately.ConclusionsThese outcomes will provide an indication for breeding high-yielding rice varieties in the immediate future.

Genome-Wide Association Study for Body Length, Body Height, and Total Teat Number in Large White Pigs

Body length, body height, and total teat number are economically important traits in pig breeding, as these traits are usually associated with the growth, reproductivity, and longevity potential of piglets. Here, we report a genetic analysis of these traits using a population comprising 2,068 Large White pigs. A genotyping-by-sequencing (GBS) approach was used to provide high-density genome-wide SNP discovery and genotyping. Univariate and bivariate animal models were used to estimate heritability and genetic correlations. The results showed that heritability estimates for body length, body height, and total teat number were 0.25 ± 0.04, 0.11 ± 0.03, and 0.22 ± 0.04, respectively. The genetic correlation between body length and body height exhibited a strongly positive correlation (0.63 ± 0.15), while a positive but low genetic correlation was observed between total teat number and body length. Furthermore, we used two different genome-wide association study (GWAS) approaches: single-locus GWAS and weighted single-step GWAS (WssGWAS), to identify candidate genes for these traits. Single-locus GWAS detected 76, 13, and 29 significant single-nucleotide polymorphisms (SNPs) associated with body length, body height, and total teat number. Notably, the most significant SNP (S17_15781294), which is located 20 kb downstream of the BMP2 gene, explained 9.09% of the genetic variance for body length traits, and it also explained 9.57% of the genetic variance for body height traits. In addition, another significant SNP (S7_97595973), which is located in the ABCD4 gene, explained 8.92% of the genetic variance for total teat number traits. GWAS results for these traits identified some candidate genomic regions, such as SSC6: 14.96–15.02 Mb, SSC7: 97.18–98.18 Mb, SSC14: 128.29–131.15 Mb, SSC17: 15.39–17.27 Mb, and SSC17: 22.04–24.15 Mb, providing a starting point for further inheritance research. Most quantitative trait loci were detected by single-locus GWAS and WssGWAS. These findings reveal the complexity of the genetic mechanism of the three traits and provide guidance for subsequent genetic improvement through genome selection.

Uncovering genomic regions controlling plant architectural traits in hexaploid wheat using different GWAS models

Wheat is a major food crop worldwide. The plant architecture is a complex trait mostly influenced by plant height, tiller number, and leaf morphology. Plant height plays a crucial role in lodging and thus affects yield and grain quality. In this study, a wheat population was genotyped by using Illumina iSelect 90K single nucleotide polymorphism (SNP) assay and finally 22,905 high-quality SNPs were used to perform a genome-wide association study (GWAS) for plant architectural traits employing four multi-locus GWAS (ML-GWAS) and three single-locus GWAS (SL-GWAS) models. As a result, 174 and 97 significant SNPs controlling plant architectural traits were detected by ML-GWAS and SL-GWAS methods, respectively. Among these SNP makers, 43 SNPs were consistently detected, including seven across multiple environments and 36 across multiple methods. Interestingly, five SNPs (Kukri_c34553_89, RAC875_c8121_1490, wsnp_Ex_rep_c66315_64480362, Ku_c5191_340, and tplb0049a09_1302) consistently detected across multiple environments and methods, played a role in modulating both plant height and flag leaf length. Furthermore, candidate SNPs (BS00068592_51, Kukri_c4750_452 and BS00022127_51) constantly repeated in different years and methods associated with flag leaf width and number of tillers. We also detected several SNPs (Jagger_c6772_80, RAC875_c8121_1490, BS00089954_51, Excalibur_01167_1207, and Ku_c5191_340) having common associations with more than one trait across multiple environments. By further appraising these GWAS methods, the pLARmEB and FarmCPU models outperformed in SNP detection compared to the other ML-GWAS and SL-GWAS methods, respectively. Totally, 152 candidate genes were found to be likely involved in plant growth and development. These finding will be helpful for better understanding of the genetic mechanism of architectural traits in wheat.

Genome-Wide Association Study on Resistance to Rice Black-Streaked Dwarf Disease Caused by Rice black-streaked dwarf virus.

Rice black-streaked dwarf disease caused by Rice black-streaked dwarf virus (RBSDV) is one of the most destructive viral diseases of rice. Thus, it is imperative that resistant rice germplasms are screened for novel RBSDV-resistant genes. RBSDV resistance of a diverse global collection comprising 1,953 rice accessions was evaluated under natural conditions across 3 years. The average disease incidences of the Xian/indica (XI) subgroup were significantly lower than those of the Geng/japonica (GJ) subgroup. Interestingly, most XI-1A accessions in the Xian subgroup were significantly more susceptible than XI-1B accessions, even though XI-1A and XI-1B have a close phylogenetic relationship. Four Xian accessions stably and highly resistant to RBSDV were consistently identified in 2 years. Ten genomic regions (GRs) with 147 single nucleotide polymorphisms associated with RBSDV resistance were detected by a single-locus genome-wide association study (GWAS), of which five were repeatedly identified in a multilocus GWAS. Two previously reported GRs, grRBSDV-6.1 and grRBSDV-6.3, which were repeatedly detected as stably and highly associated with RBSDV resistance, contained 17 and seven genes, respectively, with significant differences of resistance among haplotypes. Haplotype analyses of the candidate genes LOC_Os06g03150 in grRBSDV-6.1 and LOC_Os06g31190 in grRBSDV-6.3 suggested that the former gene is mainly associated with the differentiation of resistance within the Xian subgroup and the latter gene mainly explains the difference in the resistance between Xian and Geng. Another three novel resistance GRs (grRBSDV-1.1, grRBSDV-7.1, and grRBSDV-9.1) were identified. Our findings may enhance the application of disease-resistant rice germplasms for breeding RBSDV-resistant varieties.

Loci Associated With Antibody Response in Feral Swine (Sus scrofa) Infected With Brucella suis.

Feral swine (Sus scrofa) are a destructive invasive species widespread throughout the United States that disrupt ecosystems, damage crops, and carry pathogens of concern for the health of domestic stock and humans including Brucella suis—the causative organism for swine brucellosis. In domestic swine, brucellosis results in reproductive failure due to abortions and infertility. Contact with infected feral swine poses spillover risks to domestic pigs as well as humans, companion animals, wildlife, and other livestock. Genetic factors influence the outcome of infectious diseases; therefore, genome wide association studies (GWAS) of differential immune responses among feral swine can provide an understanding of disease dynamics and inform management to prevent the spillover of brucellosis from feral swine to domestic pigs. We sought to identify loci associated with differential antibody responses among feral swine naturally infected with B. suis using a case-control GWAS. Tissue, serum, and genotype data (68,516 bi-allelic single nucleotide polymorphisms) collected from 47 feral swine were analyzed in this study. The 47 feral swine were culture positive for Brucella spp. Of these 47, 16 were antibody positive (cases) whereas 31 were antibody negative (controls). Single-locus GWAS were performed using efficient mixed-model association eXpedited (EMMAX) methodology with three genetic models: additive, dominant, and recessive. Eight loci associated with seroconversion were identified on chromosome 4, 8, 9, 10, 12, and 18. Subsequent bioinformatic analyses revealed nine putative candidate genes related to immune function, most notably phagocytosis and induction of an inflammatory response. Identified loci and putative candidate genes may play an important role in host immune responses to B. suis infection, characterized by a detectable bacterial presence yet a differential antibody response. Given that antibody tests are used to evaluate brucellosis infection in domestic pigs and for disease surveillance in invasive feral swine, additional studies are needed to fully understand the genetic component of the response to B. suis infection and to more effectively translate estimates of Brucella spp. antibody prevalence among feral swine to disease control management action.

Increased power from conditional bacterial genome-wide association identifies macrolide resistance mutations in Neisseria gonorrhoeae

The emergence of resistance to azithromycin complicates treatment of Neisseria gonorrhoeae, the etiologic agent of gonorrhea. Substantial azithromycin resistance remains unexplained after accounting for known resistance mutations. Bacterial genome-wide association studies (GWAS) can identify novel resistance genes but must control for genetic confounders while maintaining power. Here, we show that compared to single-locus GWAS, conducting GWAS conditioned on known resistance mutations reduces the number of false positives and identifies a G70D mutation in the RplD 50S ribosomal protein L4 as significantly associated with increased azithromycin resistance (p-value = 1.08 × 10−11). We experimentally confirm our GWAS results and demonstrate that RplD G70D and other macrolide binding site mutations are prevalent (present in 5.42% of 4850 isolates) and widespread (identified in 21/65 countries across two decades). Overall, our findings demonstrate the utility of conditional associations for improving the performance of microbial GWAS and advance our understanding of the genetic basis of macrolide resistance.

MrMLM v4.0.2: An R Platform for Multi-locus Genome-wide Association Studies

Previous studies have reported that some important loci are missed in single-locus genome-wide association studies (GWAS), especially because of the large phenotypic error in field experiments. To solve this issue, multi-locus GWAS methods have been recommended. However, only a few software packages for multi-locus GWAS are available. Therefore, we developed an R software named mrMLM v4.0.2. This software integrates mrMLM, FASTmrMLM, FASTmrEMMA, pLARmEB, pKWmEB, and ISIS EM-BLASSO methods developed by our lab. There are four components in mrMLM v4.0.2, including dataset input, parameter setting, software running, and result output. The fread function in data.table is used to quickly read datasets, especially big datasets, and the doParallel package is used to conduct parallel computation using multiple CPUs. In addition, the graphical user interface software mrMLM.GUI v4.0.2, built upon Shiny, is also available. To confirm the correctness of the aforementioned programs, all the methods in mrMLM v4.0.2 and three widely-used methods were used to analyze real and simulated datasets. The results confirm the superior performance of mrMLM v4.0.2 to other methods currently available. False positive rates are effectively controlled, albeit with a less stringent significance threshold. mrMLM v4.0.2 is publicly available at BioCode (https://bigd.big.ac.cn/biocode/tools/BT007077) or R (https://cran.r-project.org/web/packages/mrMLM.GUI/index.html) as an open-source software.

Identification of QTNs and Their Candidate Genes for 100-Seed Weight in Soybean (Glycine max L.) Using Multi-Locus Genome-Wide Association Studies.

100-seed weight (100-SW) in soybeans is a yield component trait and controlled by multiple genes with different effects, but limited information is available for its quantitative trait nucleotides (QTNs) and candidate genes. To better understand the genetic architecture underlying the trait and improve the precision of marker-assisted selection, a total of 43,834 single nucleotide polymorphisms (SNPs) in 250 soybean accessions were used to identify significant QTNs for 100-SW in four environments and their BLUP values using six multi-locus and one single-locus genome-wide association study methods. As a result, a total of 218 significant QTNs were detected using multi-locus methods, whereas eight QTNs were identified by a single-locus method. Among 43 QTNs or QTN clusters identified repeatedly across various environments and/or approaches, all of them exhibited significant trait differences between their corresponding alleles, 33 were found in the genomic region of previously reported QTLs, 10 were identified as new QTNs, and three (qHSW-4-1, qcHSW-7-3, and qcHSW-10-4) were detected in all the four environments. The number of seed weight (SW) increasing alleles for each accession ranged from 8 (18.6%) to 36 (83.72%), and three accessions (Yixingwuhuangdou, Nannong 95C-5, and Yafanzaodou) had more than 35 SW increasing alleles. Among 36 homologous seed-weight genes in Arabidopsis underlying the above 43 stable QTNs, more importantly, Glyma05g34120, GmCRY1, and GmCPK11 had known seed-size/weight-related genes in soybean, and Glyma07g07850, Glyma10g03440, and Glyma10g36070 were candidate genes identified in this study. These results provide useful information for genetic foundation, marker-assisted selection, genomic prediction, and functional genomics of 100-SW.

Genome-wide association analyses identify known and novel loci for teat number in Duroc pigs using single-locus and multi-locus models

BackgroundMore teats are necessary for sows to nurse larger litters to provide immunity and nutrient for piglets prior to weaning. Previous studies have reported the strong effect of an insertion mutation in the Vertebrae Development Associated (VRTN) gene on Sus scrofa chromosome 7 (SSC7) that increased the number of thoracic vertebrae and teat number in pigs. We used genome-wide association studies (GWAS) to map genetic markers and genes associated with teat number in two Duroc pig populations with different genetic backgrounds. A single marker method and several multi-locus methods were utilized. A meta-analysis that combined the effects and P-values of 34,681 single nucleotide polymorphisms (SNPs) that were common in the results of single marker GWAS of American and Canadian Duroc pigs was conducted. We also performed association tests between the VRTN insertion and teat number in the same populations.ResultsA total of 97 SNPs were found to be associated with teat number. Among these, six, eight and seven SNPs were consistently detected with two, three and four multi-locus methods, respectively. Seven SNPs were concordantly identified between single marker and multi-locus methods. Moreover, 26 SNPs were newly found by multi-locus methods to be associated with teat number. Notably, we detected one consistent quantitative trait locus (QTL) on SSC7 for teat number using single-locus and meta-analysis of GWAS and the top SNP (rs692640845) explained 8.68% phenotypic variance of teat number in the Canadian Duroc pigs. The associations between the VRTN insertion and teat number in two Duroc pig populations were substantially weaker. Further analysis revealed that the effect of VRTN on teat number may be mediated by its LD with the true causal mutation.ConclusionsOur study suggested that VRTN insertion may not be a strong or the only candidate causal mutation for the QTL on SSC7 for teat number in the analyzed Duroc pig populations. The combination of single-locus and multi-locus GWAS detected additional SNPs that were absent using only one model. The identified SNPs will be useful for the genetic improvement of teat number in pigs by assigning higher weights to associated SNPs in genomic selection.

Benchmarking bacterial genome-wide association study methods using simulated genomes and phenotypes.

Genome-wide association studies (GWASs) have the potential to reveal the genetics of microbial phenotypes such as antibiotic resistance and virulence. Capitalizing on the growing wealth of bacterial sequence data, microbial GWAS methods aim to identify causal genetic variants while ignoring spurious associations. Bacteria reproduce clonally, leading to strong population structure and genome-wide linkage, making it challenging to separate true ‘hits’ (i.e. mutations that cause a phenotype) from non-causal linked mutations. GWAS methods attempt to correct for population structure in different ways, but their performance has not yet been systematically and comprehensively evaluated under a range of evolutionary scenarios. Here, we developed a bacterial GWAS simulator (BacGWASim) to generate bacterial genomes with varying rates of mutation, recombination and other evolutionary parameters, along with a subset of causal mutations underlying a phenotype of interest. We assessed the performance (recall and precision) of three widely used single-locus GWAS approaches (cluster-based, dimensionality-reduction and linear mixed models, implemented in plink, pyseer and gemma) and one relatively new multi-locus model implemented in pyseer, across a range of simulated sample sizes, recombination rates and causal mutation effect sizes. As expected, all methods performed better with larger sample sizes and effect sizes. The performance of clustering and dimensionality reduction approaches to correct for population structure were considerably variable according to the choice of parameters. Notably, the multi-locus elastic net (lasso) approach was consistently amongst the highest-performing methods, and had the highest power in detecting causal variants with both low and high effect sizes. Most methods reached the level of good performance (recall >0.75) for identifying causal mutations of strong effect size [log odds ratio (OR) ≥2] with a sample size of 2000 genomes. However, only elastic nets reached the level of reasonable performance (recall=0.35) for detecting markers with weaker effects (log OR ~1) in smaller samples. Elastic nets also showed superior precision and recall in controlling for genome-wide linkage, relative to single-locus models. However, all methods performed relatively poorly on highly clonal (low-recombining) genomes, suggesting room for improvement in method development. These findings show the potential for multi-locus models to improve bacterial GWAS performance. BacGWASim code and simulated data are publicly available to enable further comparisons and benchmarking of new methods.

Genome-Wide Association Study and QTL Meta-Analysis Identified Novel Genomic Loci Controlling Potassium Use Efficiency and Agronomic Traits in Bread Wheat.

Potassium use efficiency, a complex trait, directly impacts the yield potential of crop plants. Low potassium efficiency leads to a high use of fertilizers, which is not only farmer unfriendly but also deteriorates the environment. Genome-wide association studies (GWAS) are widely used to dissect complex traits. However, most studies use single-locus one-dimensional GWAS models which do not provide true information about complex traits that are controlled by multiple loci. Here, both single-locus GWAS (MLM) and multi-locus GWAS (pLARmEB, FASTmrMLM, mrMLM, FASTmrEMMA) models were used with genotyping from 90 K Infinium SNP array and phenotype derived from four normal and potassium-stress environments, which identified 534 significant marker-trait associations (MTA) for agronomic and potassium related traits: pLARmEB = 279, FASTmrMLM = 213, mrMLM = 35, MLM = 6, FASTmrEMMA = 1. Further screening of these MTA led to the detection of eleven stable loci: q1A, q1D, q2B-1, q2B-2, q2D, q4D, q5B-1, q5B-2, q5B-3, q6D, and q7A. Moreover, Meta-QTL (MQTL) analysis of four independent QTL studies for potassium deficiency in bread wheat located 16 MQTL on 13 chromosomes. One locus identified in this study (q5B-1) colocalized with an MQTL (MQTL_11), while the other ten loci were novel associations. Gene ontology of these loci identified 20 putative candidate genes encoding functional proteins involved in key pathways related to stress tolerance, sugar metabolism, and nutrient transport. These findings provide potential targets for breeding potassium stress resistant wheat cultivars and advocate the advantages of multi-locus GWAS models for studying complex traits.

Genome-Wide Association Analysis Dissects the Genetic Basis of the Grain Carbon and Nitrogen Contents in Milled Rice

BackgroundCarbon (C) and nitrogen (N) are two fundamental components of starch and protein, which are important determinants of grain yield and quality. The food preferences of consumers and the expected end-use of grains in different rice-growing regions require diverse varieties that differ in terms of the grain N content (GNC) and grain C content (GCC) of milled rice. Thus, it is important that quantitative trait loci (QTLs)/genes with large effects on the variation of GNC and GCC are identified in breeding programs.ResultsTo dissect the genetic basis of the variation of GNC and GCC in rice, the Dumas combustion method was used to analyze 751 diverse accessions regarding the GNC, GCC, and C/N ratio of the milled grains. The GCC and GNC differed significantly among the rice subgroups, especially between Xian/Indica (XI) and Geng/Japonica (GJ). Interestingly, in the GJ subgroup, the GNC was significantly lower in modern varieties (MV) than in landraces (LAN). In the XI subgroup, the GCC was significantly higher in MV than in LAN. One, six, and nine QTLs, with 55 suggestively associated single nucleotide polymorphisms, were detected for the GNC, GCC, and C/N ratio in three panels during a single-locus genome-wide association study (GWAS). Three of these QTLs were also identified in a multi-locus GWAS. We screened 113 candidate genes in the 16 QTLs in gene-based haplotype analyses. Among these candidate genes, LOC_Os01g06240 at qNC-1.1, LOC_Os05g33300 at qCC-5.1, LOC_Os01g04360 at qCN-1.1, and LOC_Os05g43880 at qCN-5.2 may partially explain the significant differences between the LAN and MV. These candidate genes should be cloned and may be useful for molecular breeding to rapidly improve the GNC, GCC, and C/N ratio of rice.ConclusionsOur findings represent valuable information regarding the genetic basis of the GNC and GCC and may be relevant for enhancing the application of favorable haplotypes of candidate genes for the molecular breeding of new rice varieties with specific grain N and C contents.

Comparing Single-SNP, Multi-SNP, and Haplotype-Based Approaches in Association Studies for Major Traits in Barley.

The multiple single nucleotide polymorphism (multi-SNP) and haplotype-based approaches that jointly consider multiple markers unveiled a larger number of associations, some of which were shared with the single-SNP approach. A larger overlap of quantitative trait loci (QTLs) between the single-SNP and haplotype-based approaches was obtained than with the multi-SNP approach. Despite a limited overlap between the QTLs detected by these approaches, each uncovered QTLs reported previously, suggesting that each approach is capable of uncovering a different subset of QTLs. We demonstrated the efficiency of an integrated genome-wide association study (GWAS) procedure, combining single-locus and multilocus approaches to improve the capacity and reliability of association analysis to detect key QTLs. The efficiency of barley breeding programs may be improved by the practical use of QTLs identified in this study. Genome-wide association studies (GWAS) have been widely used to identify quantitative trait loci (QTLs) underlying complex agronomic traits. The conventional GWAS model is based on a single-locus model, which may prove inaccurate if a trait is controlled by multiple loci, which is the case for most agronomic traits in barley (Hordeum vulgare L.). Additionally, an individual single nucleotide polymorphism (SNP) will prove incapable of capturing underlying allelic diversity. A multilocus model could potentially represent a better alternative for QTL identification. This study aimed to explore different GWAS approaches (single-SNP, multi-SNP, and haplotype-based) to establish SNP-trait associations and to potentially describe the complex genetic architecture of seven key traits in spring barley. The multi-SNP and haplotype-based approaches unveiled a larger number of significant associations, some of which were shared with the single-SNP approach. Globally, the multi-SNP approach explained more of the phenotypic variance (cumulative R2 ) and provided the best fit with the genetic model [Bayesian information criterion (BIC)]. Compared with the multi-SNP approach, the single-SNP and haplotype-based approaches were relatively similar in terms of cumulative R2 and BIC, with an improvement with the haplotype-based approach. Despite limited overlap between detected QTLs, each approach discovered QTLs that had been validated previously, suggesting that each approach can uncover a different subset of QTLs. An integrated GWAS procedure, considering single-locus and multilocus GWAS approaches jointly, may improve the capacity of association studies to detect key QTLs and to provide a more complete picture of the genetic architecture of complex traits in barley.

Genome-wide association studies in the genetic dissection of ovule number, seed number, and seed weight in Brassica napus L.

Ovule number (ON), seed per silique (SS), and thousand seed weight (TSW) are the most important but complex traits affecting yield. This study undertook genome-wide association studies (GWAS) on 521 accessions of rapeseed genotyped with the Brassica 60 K SNP array by six multi-locus GWAS (ML-GWAS) and four single-locus GWAS (SL-GWAS) methods. The findings of our study showed that 280 and 31 significant quantitative trait nucleotides/loci (QTNs/QTLs) were detected above six multi-locus and four single-locus models, respectively. Among these sequences, 74 common significant QTNs were repeatedly detected by more than three ML-GWAS models and in multiple environments. Among the QTNs, 26 were detected via multiple environments, while 13 were detected via multiple methods and environments. Interestingly, 119 QTNs were detected by a single model (pLARmEB), demonstrating that this model is largely reliable and stable. However, in SL-GWAS, the GLM model detected the highest number of QTNs. The distribution of the superior allele results showed that, among 74 common significant QTNs, 28 QTNs were >50%, while 45 QTNs were <50%. The highest percentage of superior alleles indicates their probable involvement in yield-determining traits (YDTs) improvement which may facilitate marker-assisted selection. Furthermore, on the basis of common significant QTNs, 42 candidate genes were detected. We strongly believe that the genetic manipulation of these putative genes may further improve rapeseed molecular breeding for the creation of eco-friendly cultivars with improved yield. The results obtained by this strategy may lead to a breakthrough in rapeseed production at the industrial level.