Multi-locus Genome-wide Association Studies Research Articles

Insights into the genetic architecture of Phytophthora capsici root rot resistance in chile pepper (Capsicum spp.) from multi-locus genome-wide association study

BackgroundPhytophthora root rot, a major constraint in chile pepper production worldwide, is caused by the soil-borne oomycete, Phytophthora capsici. This study aimed to detect significant regions in the Capsicum genome linked to Phytophthora root rot resistance using a panel consisting of 157 Capsicum spp. genotypes. Multi-locus genome wide association study (GWAS) was conducted using single nucleotide polymorphism (SNP) markers derived from genotyping-by-sequencing (GBS). Individual plants were separately inoculated with P. capsici isolates, ‘PWB-185’, ‘PWB-186’, and ‘6347’, at the 4–8 leaf stage and were scored for disease symptoms up to 14-days post-inoculation. Disease scores were used to calculate disease parameters including disease severity index percentage, percent of resistant plants, area under disease progress curve, and estimated marginal means for each genotype.ResultsMost of the genotypes displayed root rot symptoms, whereas five accessions were completely resistant to all the isolates and displayed no symptoms of infection. A total of 55,117 SNP markers derived from GBS were used to perform multi-locus GWAS which identified 330 significant SNP markers associated with disease resistance. Of these, 56 SNP markers distributed across all the 12 chromosomes were common across the isolates, indicating association with more durable resistance. Candidate genes including nucleotide-binding site leucine-rich repeat (NBS-LRR), systemic acquired resistance (SAR8.2), and receptor-like kinase (RLKs), were identified within 0.5 Mb of the associated markers.ConclusionsResults will be used to improve resistance to Phytophthora root rot in chile pepper by the development of Kompetitive allele-specific markers (KASP®) for marker validation, genomewide selection, and marker-assisted breeding.

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.

Multi-locus genome-wide association studies reveal the dynamic genetic architecture of flowering time in chrysanthemum.

The dynamic genetic architecture of flowering time in chrysanthemum was elucidated by GWAS. Thirty-six known genes and 14 candidate genes were identified around the stable QTNs and QEIs, among which ERF-1 was highlighted. Flowering time (FT) adaptation is one of the major breeding goals in chrysanthemum, a multipurpose ornamental plant. In order to reveal the dynamic genetic architecture of FT in chrysanthemum, phenotype investigation of ten FT-related traits was conducted on 169 entries in 2 environments. The broad-sense heritability of five non-conditional FT traits, i.e., budding (FBD), visible coloring (VC), early opening (EO), full-bloom (OF) and decay period (DP), ranged from 56.93 to 84.26%, which were higher than that of the five derived conditional FT traits (38.51-75.13%). The phenotypic variation coefficients of OF_EO and DP_OF were relatively large ranging from 30.59 to 36.17%. Based on 375,865 SNPs, the compressed variance component mixed linear model 3VmrMLM was applied for a multi-locus genome-wide association study (GWAS). As a result, 313 quantitative trait nucleotides (QTNs) were identified for the non-conditional FT traits in single-environment analysis, while 119 QTNs and 67 QTN-by-environment interactions (QEIs) were identified in multi-environment analysis. As for the conditional traits, 343 QTNs were detected in single-environment analysis, and 119 QTNs and 83 QEIs were identified in multi- environment analysis. Among the genes around stable QTNs and QEIs, 36 were orthologs of known FT genes in Arabidopsis and other plants; 14 candidates were mined by combining the transcriptomics data and functional annotation, including ERF-1, ACA10, and FOP1. Furthermore, the haplotype analysis of ERF-1 revealed six elite accessions with extreme FBD. Our findings contribute to the understanding of dynamic genetic architecture of FT and provide valuable resources for future chrysanthemum molecular breeding programs.

Multi-locus genome-wide association study reveal genomic regions underlying root system architecture traits in Ethiopian sorghum germplasm.

The identification of genomic regions underlying the root system architecture (RSA) is vital for improving crop abiotic stress tolerance. To improve sorghum (Sorghum bicolor L. Moench) for environmental stress tolerance, information on genetic variability and genomic regions linked to RSA traits is paramount. The aim of this study was, therefore, to investigate common quantitative trait nucleotides (QTNs) via multiple methodologies and identify genomic regions linked to RSA traits in a panel of 274 Ethiopian sorghum accessions. Multi-locus genome-wide association study was conducted using 265,944 high-quality single nucleotide polymorphism markers. Considering the QTN detected by at least three different methods, a total of 17 reliable QTNs were found to be significantly associated with root angle, number, length, and dry weight. Four QTNs were detected on chromosome SBI-05, followed by SBI-01 and SBI-02 with three QTNs each. Among the 17 QTNs, 11 are colocated with previously identified root traits quantitative trait loci and the remaining six are genome regions with novel genes. A total of 118 genes are colocated with these up- and down-streams of the QTNs. Moreover, five QTNs were found intragenic. These QTNs are S5_8994835 (number of nodal roots), S10_55702393 (number of nodal roots), S1_56872999 (nodal root angle), S9_1212069 (nodal root angle), and S5_5667192 (root dry weight) intragenic regions of Sobic.005G073101, Sobic.010G198000, Sobic.001G273000, Sobic.009G013600, and Sobic.005G054700, respectively. Particularly, Sobic.005G073101, Sobic.010G198000, and Sobic.009G013600 were found responsible for the plant growth hormone-induced RSA. These genes may regulate root development in the seedling stage. Further analysis on these genes might be important to explore the genetic structure of RSA of sorghum.

Multi-model genome-wide association studies for appearance quality in rice.

Improving the quality of the appearance of rice is critical to meet market acceptance. Mining putative quality-related genes has been geared towards the development of effective breeding approaches for rice. In the present study, two SL-GWAS (CMLM and MLM) and three ML-GWAS (FASTmrEMMA, mrMLM, and FASTmrMLM) genome-wide association studies were conducted in a subset of 3K-RGP consisting of 198 rice accessions with 553,831 SNP markers. A total of 594 SNP markers were identified using the mixed linear model method for grain quality traits. Additionally, 70 quantitative trait nucleotides (QTNs) detected by the ML-GWAS models were strongly associated with grain aroma (AR), head rice recovery (HRR, %), and percentage of grains with chalkiness (PGC, %). Finally, 39 QTNs were identified using single- and multi-locus GWAS methods. Among the 39 reliable QTNs, 20 novel QTNs were identified for the above-mentioned three quality-related traits. Based on annotation and previous studies, four functional candidate genes (LOC_Os01g66110, LOC_Os01g66140, LOC_Os07g44910, and LOC_Os02g14120) were found to influence AR, HRR (%), and PGC (%), which could be utilized in rice breeding to improve grain quality traits.

Genome-wide association with transcriptomics reveals a shade-tolerance gene network in soybean

Shade tolerance is essential for soybeans in inter/relay cropping systems. A genome-wide association study (GWAS) integrated with transcriptome sequencing was performed to identify genes and construct a genetic network governing the trait in a set of recombinant inbred lines derived from two soybean parents with contrasting shade tolerance. An improved GWAS procedure, restricted two-stage multi-locus genome-wide association study based on gene/allele sequence markers (GASM-RTM-GWAS), identified 140 genes and their alleles associated with shade-tolerance index (STI), 146 with relative pith cell length (RCL), and nine with both. Annotation of these genes by biological categories allowed the construction of a protein–protein interaction network by 187 genes, of which half were differentially expressed under shading and non-shading conditions as well as at different growth stages. From the identified genes, three ones jointly identified for both traits by both GWAS and transcriptome and two genes with maximum links were chosen as beginners for entrance into the network. Altogether, both STI and RCL gene systems worked for shade-tolerance with genes interacted each other, this confirmed that shade-tolerance is regulated by more than single group of interacted genes, involving multiple biological functions as a gene network.

The Genetic Dissection of Nitrogen Use-Related Traits in Flax (Linum usitatissimum L.) at the Seedling Stage through the Integration of Multi-Locus GWAS, RNA-seq and Genomic Selection.

Nitrogen (N), the most important macro-nutrient for plant growth and development, is a key factor that determines crop yield. Yet its excessive applications pollute the environment and are expensive. Hence, studying nitrogen use efficiency (NUE) in crops is fundamental for sustainable agriculture. Here, an association panel consisting of 123 flax accessions was evaluated for 21 NUE-related traits at the seedling stage under optimum N (N+) and N deficiency (N-) treatments to dissect the genetic architecture of NUE-related traits using a multi-omics approach integrating genome-wide association studies (GWAS), transcriptome analysis and genomic selection (GS). Root traits exhibited significant and positive correlations with NUE under N- conditions (r = 0.33 to 0.43, p < 0.05). A total of 359 QTLs were identified, accounting for 0.11% to 23.1% of the phenotypic variation in NUE-related traits. Transcriptomic analysis identified 1034 differentially expressed genes (DEGs) under contrasting N conditions. DEGs involved in N metabolism, root development, amino acid transport and catabolism and others, were found near the QTLs. GS models to predict NUE stress tolerance index (NUE_STI) trait were tested using a random genome-wide SNP dataset and a GWAS-derived QTLs dataset. The latter produced superior prediction accuracy (r = 0.62 to 0.79) compared to the genome-wide SNP marker dataset (r = 0.11) for NUE_STI. Our results provide insights into the QTL architecture of NUE-related traits, identify candidate genes for further studies, and propose genomic breeding tools to achieve superior NUE in flax under low N input.

Multi-locus genome-wide association study and genomic prediction for flowering time in chrysanthemum.

Multi-locus GWAS detected several known and candidate genes responsible for flowering time in chrysanthemum. The associations could greatly increase the predictive ability of genome selection that accelerates the possible application of GS in chrysanthemum breeding. Timely flowering is critical for successful reproduction and determines the economic value for ornamental plants.To investigate the genetic architecture of flowering time in chrysanthemum, a multi-locus genome-wide association study (GWAS) was performed using a collection of 200 accessions and 330,710 single-nucleotide polymorphisms (SNPs) via 3VmrMLM method. Five flowering time traits including budding (FBD), visible colouring (VC), early opening (EO), full-bloom (OF) and senescing (SF) stages, plus five derived conditional traits were recorded in two environments. Extensive phenotypic variations were observed for these flowering time traits with coefficients of variation ranging from 6.42 to 38.27%, and their broad-sense heritability ranged from 71.47 to 96.78%. GWAS revealed 88 stable quantitative trait nucleotides (QTNs) and 93 QTN-by-environment interactions (QEIs) associated with flowering time traits, accounting for 0.50-8.01% and 0.30-10.42% of the phenotypic variation, respectively. Amongst the genes around these stable QTNs and QEIs, 21 and 10 were homologous to known flowering genes in Arabidopsis; 20 and 11 candidate genes were mined by combining the functional annotation and transcriptomics data, respectively, such as MYB55, FRIGIDA-like, WRKY75 and ANT. Furthermore, genomic selection (GS) was assessed using three models and seven unique marker datasets. We found the prediction accuracy (PA) using significant SNPs identified by GWAS under SVM model exhibited the best performance with PA ranging from 0.90 to 0.95. Our findings provide new insights into the dynamic genetic architecture of flowering time and the identified significant SNPs and candidate genes will accelerate the future molecular improvement of chrysanthemum.

Single Nucleotide Polymorphisms Associated with AA-Amyloidosis in Siamese and Oriental Shorthair Cats.

AA-amyloidosis in Siamese and Oriental shorthair cats is a lethal condition in which amyloid deposits accumulate systemically, especially in the liver and the thyroid gland. The age at death of affected cats varies between one and seven years. A previous study indicated a complex mode of inheritance involving a major locus. In the present study, we performed a multi-locus genome-wide association study (GWAS) using five methods (mrMLM, FASTmrMLM, FASTmrEMMA, pLARmEB and ISIS EM-BLASSO) to identify variants associated with AA-amyloidosis in Siamese/Oriental cats. We genotyped 20 affected mixed Siamese/Oriental cats from a cattery and 48 healthy controls from the same breeds using the Illumina Infinium Feline 63 K iSelect DNA array. The multi-locus GWAS revealed eight significantly associated single nucleotide polymorphisms (SNPs) on FCA A1, D1, D2 and D3. The genomic regions harboring these SNPs contain 55 genes, of which 3 are associated with amyloidosis in humans or mice. One of these genes is SAA1, which encodes for a member of the Serum Amyloid A family, the precursor protein of Amyloid A, and a mutation in the promotor of this gene causes hereditary AA-amyloidosis in humans. These results provide novel knowledge regarding the complex genetic background of hereditary AA-amyloidosis in Siamese/Oriental cats and, therefore, contribute to future genomic studies of this disease in cats.

The Identification of a Quantative Trait Loci-Allele System of Antixenosis against the Common Cutworm (Spodoptera litura Fabricius) at the Seedling Stage in the Chinese Soybean Landrace Population.

Common cutworm (CCW) is an omnivorous insect causing severe yield losses in soybean crops. The seedling-stage mini-tray identification system with the damaged leaf percentage (DLP) as an indicator was used to evaluate antixenosis against CCW in the Chinese soybean landrace population (CSLRP) under three environments. Using the innovative restricted two-stage multi-locus genome-wide association study procedure (RTM-GWAS), 86 DLP QTLs with 243 alleles (2-11/QTL) were identified, including 66 main-effect loci with 203 alleles and 57 QTL-environment interaction loci with 172 alleles. Among the main-effect loci, 12 large-contribution loci (R2 ≥ 1%) explained 25.45% of the phenotypic variation (PV), and 54 small-contribution loci (R2 < 1%) explained 16.55% of the PV. This indicates that the CSLRP can be characterized with a DLP QTL-allele system complex that has not been found before, except for a few individual QTLs without alleles involved. From the DLP QTL-allele matrix, the recombination potentials expressed in the 25th percentile of the DLP of all possible crosses were predicted to be reduced by 41.5% as the maximum improvement and 14.2% as the maximum transgression, indicating great breeding potential in the antixenosis of the CSLRP. From the QTLs, 62 candidate genes were annotated, which were involved in eight biological function categories as a gene network of the DLP. Changing from susceptible to moderate plus resistant varieties in the CSLRP, 26 QTLs had 32 alleles involved, in which 19 genes were annotated from 25 QTL-alleles, including eight increased negative alleles on seven loci and 11 decreased positive alleles on 11 loci, showing the major genetic constitution changes for the antixenosis enhancement at the seedling stage in the CSLRP.

Genetic analysis of genotype-specific parameters in the DSSTA-CROPGRO-soybean phenology simulation model via a multi-GWAS method

ContextGenotype-specific parameters (GSP) are the fundamental elements for phenotypic prediction in the DSSAT-CROPGRO soybean phenology simulation model (DCS-PSM). Yet, it is unclear about the quantitative relationship between GSP and genomics, which limits the crop models used for breeding research. ObjectiveExploring the genetic composition of GSP related to four soybean phenological traits in DCS-PSM for providing a basis for whether GSP are genetic coefficient. Highlighting the breeding value of DCS-PSM from a genetic perspective and strengthening its genetic understanding, which could support a further crop model-assisted molecular breeding. MethodsWe conducted a multi-genome-wide association study (GWAS) method by intrgrating the single-locus with multi-locus GWAS models for GSP associated quantitative trait nucleotides (QTN) detection. And the heritability, functional annotation of QTN and genotype-environment interaction (GEI) effect were further analyzed. Results(i). All GSP had high broad-sense heritability, ranging form 88.20–98.80%. (ii). Multi-locus GWAS significantly increased the QTN detection efficiency, the GSP associated QTN detected by multi-GWAS method contributed 22.11–50.25% phenotypic variation explained (PVE) of GSP. Besides, the GSP associated QTN accounted for a high PVE of four phenological traits, ranging from 54.25% to 74.75%. (iii). Most QTN overlapped with the function-known genes which regulate soybean phenology through photoperiod-temperature pathways. (iv). GSP associated QTN could reveal the GEI effect and contributed 1.45–6.27% PVE of soybean phenology. ConclusionsThis study indicated that the genetic factor was the primary source of GSP’s variation, and the multi-GWAS method could be an alternative way for the genetic analysis of GSP. Phenology-related GSP had genetic potential and could be regarded as genetic coefficients and breeding targets.

Multi-locus genome-wide association studies reveal the genetic architecture of Fusarium head blight resistance in durum wheat.

Durum wheat is more susceptible to Fusarium head blight (FHB) than other types or classes of wheat. The disease is one of the most devastating in wheat; it reduces yield and end-use quality and contaminates the grain with fungal mycotoxins such as deoxynivalenol (DON). A panel of 265 Canadian and European durum wheat cultivars, as well as breeding and experimental lines, were tested in artificially inoculated field environments (2019-2022, inclusive) and two greenhouse trials (2019 and 2020). The trials were assessed for FHB severity and incidence, visual rating index, Fusarium-damaged kernels, DON accumulation, anthesis or heading date, maturity date, and plant height. In addition, yellow pigment and protein content were analyzed for the 2020 field season. To capture loci underlying FHB resistance and related traits, GWAS was performed using single-locus and several multi-locus models, employing 13,504 SNPs. Thirty-one QTL significantly associated with one or more FHB-related traits were identified, of which nine were consistent across environments and associated with multiple FHB-related traits. Although many of the QTL were identified in regions previously reported to affect FHB, the QTL QFhb-3B.2, associated with FHB severity, incidence, and DON accumulation, appears to be novel. We developed KASP markers for six FHB-associated QTL that were consistently detected across multiple environments and validated them on the Global Durum Panel (GDP). Analysis of allelic diversity and the frequencies of these revealed that the lines in the GDP harbor between zero and six resistance alleles. This study provides a comprehensive assessment of the genetic basis of FHB resistance and DON accumulation in durum wheat. Accessions with multiple favorable alleles were identified and will be useful genetic resources to improve FHB resistance in durum breeding programs through marker-assisted recurrent selection and gene stacking.

Genetic potential of grain-related traits in rice landraces: phenomics and multi-locus association analyses.

Grain length, width, weight, and the number of grains per panicle are crucial determinants contributing to yield in cereal crops. Understanding the genetic basis of grain-related traits has been the main research object in crop science. Kerala has a collection of different rice landraces. Characterization of these valuable genetic resources for 39 distinct agro-morphological traits was carried out in two seasons from 2017 to 2019 directly in farmers field. Most characteristics were polymorphic except ligule shape, leaf angle, and panicle axis. The results of principal component analysis implied that leaf length, plant height, culm length, flag leaf length, and grain-related traits were the principal discriminatory characteristics of rice landraces. For identifying the genetic basis of key grain traits of rice, three multi locus GWAS models were performed based on 1,47,994 SNPs in 73 rice accessions. As a result, 48 quantitative trait nucleotides (QTNs) were identified to be associated with these traits. After characterization of their function and expression, 15 significant candidate genes involved in regulating grain width, number of grains per panicle, and yield were identified. The detected QTNs and candidate genes in this study could be further used for marker-assisted high-quality breeding of rice.

Genetic analysis of iron, zinc and grain yield in wheat-Aegilops derivatives using multi-locus GWAS.

Wheat is a major staple crop and helps to reduce worldwide micronutrient deficiency. Investigating the genetics that control the concentrations of iron (Fe) and zinc (Zn) in wheat is crucial. Hence, we undertook a comprehensive study aimed at elucidating the genomic regions linked to the contents of Fe and Zn in the grain. We performed the multi-locus genome-wide association (ML-GWAS) using a panel of 161 wheat-Aegilops substitution and addition lines to dissect the genomic regions controlling grain iron (GFeC), and grain zinc (GZnC) contents. The wheat panel was genotyped using 10,825 high-quality SNPs and phenotyped in three different environments (E1-E3) during 2017-2019. A total of 111 marker-trait associations (MTAs) (at p-value < 0.001) were detected that belong to all three sub-genomes of wheat. The highest number of MTAs were identified for GFeC (58), followed by GZnC (44) and yield (9). Further, six stable MTAs were identified for these three traits and also two pleiotropic MTAs were identified for GFeC and GZnC. A total of 1291 putative candidate genes (CGs) were also identified for all three traits. These CGs encode a diverse set of proteins, including heavy metal-associated (HMA), bZIP family protein, AP2/ERF, and protein previously associated with GFeC, GZnC, and grain yield. The significant MTAs and CGs pinpointed in this current study are poised to play a pivotal role in enhancing both the nutritional quality and yield of wheat, utilizing marker-assisted selection (MAS) techniques.

Identification of southern corn rust resistance QTNs in Chinese summer maize germplasm via multi-locus GWAS and post-GWAS analysis.

Southern corn rust (SCR) caused by Puccinia polysora Underw is a major disease leading to severe yield losses in China Summer Corn Belt. Using six multi-locus GWAS methods, we identified a set of SCR resistance QTNs from a diversity panel of 140 inbred lines collected from China Summer Corn Belt. Thirteen QTNs on chromosomes 1, 2, 4, 5, 6, and 8 were grouped into three types of allele effects and their associations with SCR phenotypes were verified by post-GWAS case-control sampling, allele/haplotype effect analysis. Relative resistance (RRR) and relative susceptibility (RRs) catering to its inbred carrier were estimated from single QTN and QTN-QTN combos and epistatitic effects were estimated for QTN-QTN combos. By transcriptomic annotation, a set of candidate genes were predicted to be involved in transcriptional regulation (S5_145, Zm00001d01613, transcription factor GTE4), phosphorylation (S8_123, Zm00001d010672, Pgk2- phosphoglycerate kinase 2), and temperature stress response (S6_164a/S6_164b, Zm00001d038806, hsp101, and S5_211, Zm00001d017978, cellulase25). The breeding implications of the above findings were discussed.

Exploring GWAS and genomic prediction to improve Septoria tritici blotch resistance in wheat

Septoria tritici blotch (STB) is a destructive foliar diseases threatening wheat grain yield. Wheat breeding for STB disease resistance has been identified as the most sustainable and environment-friendly approach. In this work, a panel of 316 winter wheat breeding lines from a commercial breeding program were evaluated for STB resistance at the seedling stage under controlled conditions followed by genome-wide association study (GWAS) and genomic prediction (GP). The study revealed a significant genotypic variation for STB seedling resistance, while disease severity scores exhibited a normal frequency distribution. Moreover, we calculated a broad-sense heritability of 0.62 for the trait. Nine single- and multi-locus GWAS models identified 24 marker-trait associations grouped into 20 quantitative trait loci (QTLs) for STB seedling-stage resistance. The seven QTLs located on chromosomes 1B, 2A, 2B, 5B (two), 7A, and 7D are reported for the first time and could potentially be novel. The GP cross-validation analysis in the RR-BLUP model estimated the genomic-estimated breeding values (GEBVs) of STB resistance with a prediction accuracy of 0.49. Meanwhile, the GWAS assisted wRR-BLUP model improved the accuracy to 0.58. The identified QTLs can be used for marker-assisted backcrossing against STB in winter wheat. Moreover, the higher prediction accuracy recorded from the GWAS-assisted GP analysis implies its power to successfully select superior candidate lines based on their GEBVs for STB resistance.

Genome-wide association study of soluble solids content, flesh color, and fruit shape in citron watermelon.

Fruit quality traits are crucial determinants of consumers' willingness to purchase watermelon produce, making them major goals for breeding programs. There is limited information on the genetic underpinnings of fruit quality traits in watermelon. A total of 125 citron watermelon (Citrullus amarus) accessions were genotyped using single nucleotide polymorphisms (SNPs) molecular markers generated via whole-genome resequencing. A total of 2,126,759 genome-wide SNP markers were used to uncover marker-trait associations using single and multi-locus GWAS models. High broad-sense heritability for fruit quality traits was detected. Correlation analysis among traits revealed positive relationships, with the exception of fruit diameter and fruit shape index (ratio of fruit length to fruit diameter), which was negative. A total of 37 significant SNP markers associated with soluble solids content, flesh color, fruit length, fruit diameter, and fruit shape index traits were uncovered. These peak SNPs accounted for 2.1%-23.4% of the phenotypic variation explained showing the quantitative inheritance nature of the evaluated traits. Candidate genes relevant to fruit quality traits were uncovered on chromosomes Ca01, Ca03, Ca06, and Ca07. These significant molecular markers and candidate genes will be useful in marker-assisted breeding of fruit quality traits in watermelon.

Multi-locus genome-wide association studies reveal genetic loci and candidate genes for leaf angle traits in cut chrysanthemum

Chrysanthemum (Chrysanthemum morifolium Ramat.), one of the four most popular cut flowers in the world, occupies a large share of global floriculture trade. Leaf angle (LA) is an important factor affecting plant architecture, photosynthesis and yields of cut chrysanthemum. However, the genetic basis and related genes underlying LA remain largely unknown. In this research, phenotypes of five LA-related traits in 200 cut chrysanthemum accessions were investigated across two consecutive years (Y2020 and Y2021), and their broad-sense genetic heritabilities ranged from 68.91% to 84.70%. Leaf orientation value (LOV) exhibited the maximum phenotypic coefficient of variation of 35.66% in Y2020 and 32.34% in Y2021. A new methodological framework named 3VmrMLM was applied for multi-locus genome-wide association study (GWAS) based on a set of 330,710 single nucleotide polymorphisms (SNPs) to explore the genetic loci and candidate genes of LA. As a result, 303 quantitative trait nucleotides (QTNs) in single-environment, 120 QTNs and 75 QTN × environments interactions (QEIs) in multiple-environments were identified. Among the genes around these QTNs and QEIs, 26 known genes were considered to be associated with LA, including ARF16, IBH1, AN3, etc. Furthermore, the haplotype analysis of ARF16 revealed four accessions with elite LA. Our findings provide a fundamental understanding to the genetic architecture of LA and contribute to further gene cloning and molecular improvement of chrysanthemum plant architecture.

Multi-locus genome-wide association study of fusarium head blight in relation to days to anthesis and plant height in a spring wheat association panel.

Fusarium head blight (FHB) is a highly destructive fungal disease of wheat to which host resistance is quantitatively inherited and largely influenced by the environment. Resistance to FHB has been associated with taller height and later maturity; however, a further understanding of these relationships is needed. An association mapping panel (AMP) of 192 predominantly Canadian spring wheat was genotyped with the wheat 90K single-nucleotide polymorphism (SNP) array. The AMP was assessed for FHB incidence (INC), severity (SEV) and index (IND), days to anthesis (DTA), and plant height (PLHT) between 2015 and 2017 at three Canadian FHB-inoculated nurseries. Seven multi-environment trial (MET) datasets were deployed in a genome-wide association study (GWAS) using a single-locus mixed linear model (MLM) and a multi-locus random SNP-effect mixed linear model (mrMLM). MLM detected four quantitative trait nucleotides (QTNs) for INC on chromosomes 2D and 3D and for SEV and IND on chromosome 3B. Further, mrMLM identified 291 QTNs: 50 (INC), 72 (SEV), 90 (IND), 41 (DTA), and 38 (PLHT). At two or more environments, 17 QTNs for FHB, DTA, and PLHT were detected. Of these 17, 12 QTNs were pleiotropic for FHB traits, DTA, and PLHT on chromosomes 1A, 1D, 2D, 3B, 5A, 6B, 7A, and 7B; two QTNs for DTA were detected on chromosomes 1B and 7A; and three PLHT QTNs were located on chromosomes 4B and 6B. The 1B DTA QTN and the three pleiotropic QTNs on chromosomes 1A, 3B, and 6B are potentially identical to corresponding quantitative trait loci (QTLs) in durum wheat. Further, the 3B pleiotropic QTN for FHB INC, SEV, and IND co-locates with TraesCS3B02G024900 within the Fhb1 region on chromosome 3B and is ~3 Mb from a cloned Fhb1 candidate gene TaHRC. While the PLHT QTN on chromosome 6B is putatively novel, the 1B DTA QTN co-locates with a disease resistance protein located ~10 Mb from a Flowering Locus T1-like gene TaFT3-B1, and the 7A DTA QTN is ~5 Mb away from a maturity QTL QMat.dms-7A.3 of another study. GWAS and QTN candidate genes enabled the characterization of FHB resistance in relation to DTA and PLHT. This approach should eventually generate additional and reliable trait-specific markers for breeding selection, in addition to providing useful information for FHB trait discovery.

Genetic dissection of thousand-seed weight in linseed (Linum usitatissimum L.) using multi-locus genome-wide association study.

Flaxseed/linseed is an important oilseed crop having applications in the food, nutraceutical, and paint industry. Seed weight is one of the most crucial determinants of seed yield in linseed. Here, quantitative trait nucleotides (QTNs) associated with thousand-seed weight (TSW) have been identified using multi-locus genome-wide association study (ML-GWAS). Field evaluation was carried out in five environments in multi-year-location trials. SNP genotyping information of the AM panel of 131 accessions comprising 68,925 SNPs was employed for ML-GWAS. From the six ML-GWAS methods employed, five methods helped identify a total of 84 unique significant QTNs for TSW. QTNs identified in ≥ 2 methods/environments were designated as stable QTNs. Accordingly, 30 stable QTNs have been identified for TSW accounting up to 38.65% trait variation. Alleles with positive effect on trait were analyzed for 12 strong QTNs with r 2 ≥ 10.00%, which showed significant association of specific alleles with higher trait value in three or more environments. A total of 23 candidate genes have been identified for TSW, which included B3 domain-containing transcription factor, SUMO-activating enzyme, protein SCARECROW, shaggy-related protein kinase/BIN2, ANTIAUXIN-RESISTANT 3, RING-type E3 ubiquitin transferase E4, auxin response factors, WRKY transcription factor, and CBS domain-containing protein. In silico expression analysis of candidate genes was performed to validate their possible role in different stages of seed development process. The results from this study provide significant insight and elevate our understanding on genetic architecture of TSW trait in linseed.