Genetic dissection of yield-related traits via genome-wide association analysis across multiple environments in wild soybean (Glycine soja Sieb. and Zucc.)

Genome-wide association studies reveal novel QTLs for agronomic traits in soybean

Introduction

Soybean, as a globally significant crop, has garnered substantial attention due to its agricultural importance. The utilization of molecular approaches to enhance grain yield in soybean has gained popularity.

Methods

In this study, we conducted a genome-wide association study (GWAS) using 156 Chinese soybean accessions over a two-year period. We employed the general linear model (GLM) and the mixed linear model (MLM) to analyze three agronomic traits: pod number, grain number, and grain weight.

Results

Our findings revealed significant associations between qgPNpP-98, qgGNpP-89 and qgHGW-85 QTLs and pod number, grain number, and grain weight, respectively. These QTLs were identified on chromosome 16, a region spanning 413171bp exhibited associations with all three traits.

Discussion

These QTL markers identified in this study hold potential for improving yield and agronomic traits through marker-assisted selection and genomic selection in breeding programs.

Genetic dissection of ten photosynthesis-related traits based on InDel- and SNP-GWAS in soybean

KEY MESSAGE

A total of 416 InDels and 112 SNPs were significantly associated with soybean photosynthesis-related traits. GmIWS1 and GmCDC48 might be related to chlorophyll fluorescence and gas-exchange parameters, respectively. Photosynthesis is one of the main factors determining crop yield. A better understanding of the genetic architecture for photosynthesis is of great significance for soybean yield improvement. Our previous studies identified 5,410,112 single nucleotide polymorphisms (SNPs) from the resequencing data of 219 natural soybean accessions. Here, we identified 634,106 insertions and deletions (InDels) from these 219 accessions and used these InDel variations to perform principal component and linkage disequilibrium analysis of this population. The genome-wide association study (GWAS) were conducted on six chlorophyll fluorescence parameters (chlorophyll content, light energy absorbed per reaction center, quantum yield for electron transport, probability that a trapped exciton moves an electron into the electron transport chain beyond primary quinone acceptor, maximum quantum yield of photosystem II primary photochemistry in the dark-adapted state, performance index on absorption basis) and four gas-exchange parameters (intercellular carbon dioxide concentration, stomatal conductance, net photosynthesis rate, transpiration rate) and revealed 416 significant InDels and 112 significant SNPs. Based on GWAS results, GmIWS1 (encoding a transcription elongation factor) and GmCDC48 (encoding a cell division cycle protein) with the highest expression in the mapping region were determined as the candidate genes responsible for chlorophyll fluorescence and gas-exchange parameters, respectively. Further identification of favorable haplotypes with higher photosynthesis, seed weight and seed yield were carried out for GmIWS1 and GmCDC48. Overall, this study revealed the natural variations and candidate genes underlying the photosynthesis-related traits based on abundant phenotypic and genetic data, providing valuable insights into the genetic mechanisms controlling photosynthesis and yield in soybean.

Molecular, genetic, and genomic basis of seed size and yield characteristics in soybean

Soybean (Glycine max L. Merr.) is a crucial oilseed cash crop grown worldwide and consumed as oil, protein, and food by humans and feed by animals. Comparatively, soybean seed yield is lower than cereal crops, such as maize, rice, and wheat, and the demand for soybean production does not keep up with the increasing consumption level. Therefore, increasing soybean yield per unit area is the most crucial breeding objective and is challenging for the scientific community. Moreover, yield and associated traits are extensively researched in cereal crops, but little is known about soybeans' genetics, genomics, and molecular regulation of yield traits. Soybean seed yield is a complex quantitative trait governed by multiple genes. Understanding the genetic and molecular processes governing closely related attributes to seed yield is crucial to increasing soybean yield. Advances in sequencing technologies have made it possible to conduct functional genomic research to understand yield traits' genetic and molecular underpinnings. Here, we provide an overview of recent progress in the genetic regulation of seed size in soybean, molecular, genetics, and genomic bases of yield, and related key seed yield traits. In addition, phytohormones, such as auxin, gibberellins, cytokinins, and abscisic acid, regulate seed size and yield. Hence, we also highlight the implications of these factors, challenges in soybean yield, and seed trait improvement. The information reviewed in this study will help expand the knowledge base and may provide the way forward for developing high-yielding soybean cultivars for future food demands.

Identification of genetic loci for early maturity in spring bread wheat using the association analysis and gene dissection

Background

Early maturity in spring bread wheat is highly desirable in the regions where it enables the plants to evade high temperatures and plant pathogens at the end of the growing season.

Methods

To reveal the genetic loci responsible for the maturity time association analysis was carried out based on phenotyping for an 11-year period and high-throughput SNP genotyping of a panel of the varieties contrasting for this trait. The expression of candidate genes was verified using qPCR. The association between the SNP markers and the trait was validated using the biparental F2:3 population.

Results

Our data showed that under long-day conditions, the period from seedling to maturity is mostly influenced by the time from heading to maturity, rather than the heading time. The QTLs associated with the trait were located on 2A, 3B, 4A, 5B, 7A and 7B chromosomes with the 7BL locus being the most significant and promising for its SNPs accelerated the maturity time by about 9 days. Gene dissection in this locus detected a number of candidates, the best being TraesCS7B02G391800 (bZIP9) and TraesCS7B02G412200 (photosystem II reaction center). The two genes are predominantly expressed in the flag leaf while flowering. The effect of the SNPs was verified in F2:3 population and confirmed the association of the 4A, 5B and 7BL loci with the maturity time.

The Genotypic Variability among Short-Season Soybean Cultivars for Nitrogen Fixation under Drought Stress

Soybean fixes atmospheric nitrogen through the symbiotic rhizobia bacteria that inhabit root nodules. Drought stress negatively affect symbiotic nitrogen fixation (SNF) in soybean. The main objective of this study was to identify allelic variations associated with SNF in short-season Canadian soybean varieties under drought stress. A diversity panel of 103 early-maturity Canadian soybean varieties was evaluated under greenhouse conditions to determine SNF-related traits under drought stress. Drought was imposed after three weeks of plant growth, where plants were maintained at 30% field capacity (FC) (drought) and 80% FC (well-watered) until seed maturity. Under drought stress, soybean plants had lower seed yield, yield components, seed nitrogen content, % nitrogen derived from the atmosphere (%Ndfa), and total seed nitrogen fixed compared to those under well-watered conditions. Significant genotypic variability among soybean varieties was found for yield, yield parameters, and nitrogen fixation traits. A genome-wide association study (GWAS) was conducted using 2.16 M single nucleotide single nucleotide polymorphisms (SNPs) for different yield and nitrogen fixation related parameters for 30% FC and their relative performance (30% FC/80% FC). In total, five quantitative trait locus (QTL) regions, including candidate genes, were detected as significantly associated with %Ndfa under drought stress and relative performance. These genes can potentially aid in future breeding efforts to develop drought-resistant soybean varieties.

Genome-wide association study identifies candidate genes and favorable haplotypes for seed yield in Brassica napus

Oilseed rape (Brassica napus L.) is one of the most essential oil crops. Genetic improvement of seed yield (SY) is a major aim of B. napus breeding. Several studies have been reported on the genetic mechanisms of SY of B. napus. Here, a genome-wide association study (GWAS) of SY was conducted using a panel of 403 natural accessions of B. napus, with more than five million high-quality single-nucleotide polymorphisms (SNPs). A total of 1773 significant SNPs were detected associated with SY, and 783 significant SNPs were co-located with previously reported QTLs. The lead SNPs chrA01__8920351 and chrA02__4555979 were jointly detected in Trial 2_2 and Trial 2_mean value, and in Trial 1_2 and Trial 1_mean value, respectively. Subsequently, two candidate genes of BnaA01g17200D and BnaA02g08680D were identified through combining transcriptome, candidate gene association analysis, and haplotype analysis. BnaA09g10430D detected through lead SNP chrA09__5160639 was associated with SY of B. napus. Our results provide valuable information for studying the genetic control of seed yield in B. napus and valuable genes, haplotypes, and cultivars resources for the breeding of high seed yield B. napus cultivars.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01332-6.

CALCIUM-DEPENDENT PROTEIN KINASE38 regulates flowering time and common cutworm resistance in soybean

Photoperiod-sensitive plants such as soybean (Glycine max) often face threats from herbivorous insects throughout their whole growth period and especially during flowering; however, little is known about the relationship between plant flowering and insect resistance. Here, we used gene editing, multiple omics, genetic diversity and evolutionary analyses to confirm that the calcium-dependent protein kinase GmCDPK38 plays a dual role in coordinating flowering time regulation and insect resistance of soybean. Haplotype 2 (Hap2)-containing soybeans flowered later and were more resistant to the common cutworm (Spodoptera litura Fabricius) than those of Hap3. gmcdpk38 mutants with Hap3 knocked out exhibited similar flowering and resistance phenotypes as Hap2. Knocking out GmCDPK38 altered numerous flowering- and resistance-related phosphorylated proteins, genes, and metabolites. For example, the S-adenosylmethionine synthase GmSAMS1 was post-translationally upregulated in the gmcdpk38 mutants. GmCDPK38 has abundant genetic diversity in wild soybeans and was likely selected during soybean domestication. We found that Hap2 was mostly distributed at low latitudes and had a higher frequency in cultivars than in wild soybeans, while Hap3 was widely selected at high latitudes. Overall, our results elucidated that the two distinct traits (flowering time and insect resistance) are mediated by GmCDPK38.

Genome-wide association study, haplotype analysis, and genomic prediction reveal the genetic basis of yield-related traits in soybean (Glycine max L.)

Identifying the genetic components underlying yield-related traits in soybean is crucial for improving its production and productivity. Here, 211 soybean genotypes were evaluated across six environments for four yield-related traits, including seed yield per plant (SYP), number of pods per plant number of seeds per plant and 100-seed weight (HSW). Genome-wide association study (GWAS) and genomic prediction (GP) analyses were performed using 12,617 single nucleotide polymorphism markers from NJAU 355K SoySNP Array. A total of 57 SNPs were significantly associated with four traits across six environments and a combined environment using five Genome-wide association study models. Out of these, six significant SNPs were consistently identified in more than three environments using multiple GWAS models. The genomic regions (±670 kb) flanking these six consistent SNPs were considered stable QTL regions. Gene annotation and in silico expression analysis revealed 15 putative genes underlying the stable QTLs that might regulate soybean yield. Haplotype analysis using six significant SNPs revealed various allelic combinations regulating diverse phenotypes for the studied traits. Furthermore, the GP analysis revealed that accurate breeding values for the studied soybean traits is attainable at an earlier generation. Our study paved the way for increasing soybean yield performance within a short breeding cycle.

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.

Genotyping Platforms for Genome-Wide Association Studies: Options and Practical Considerations

Genome-wide association studies (GWAS) in crops requires genotyping platforms that are capable of producing accurate high density genotyping data on hundreds of plants in a cost-effective manner. Currently there are multiple commercial platforms available that are being effectively used across crops. These platforms include genotyping arrays such as the Illumina Infinium arrays and the Applied Biosystems Axiom Arrays along with a variety of resequencing methods. These methods are being used to genotype tens of thousands of markers up to millions of markers on GWAS panels. They are being used on crops with simple genomes to crops with very complex, large, polyploid genomes. Depending on the crop and the goal of the GWAS, there are several options and practical considerations to take into account when selecting a genotyping technology to ensure that the right coverage, accuracy, and cost for the study is achieved.

GmFULa improves soybean yield by enhancing carbon assimilation without altering flowering time or maturity

KEY MESSAGE

GmFULa improved soybean yield by enhancing carbon assimilation. Meanwhile, different from known yield-related genes, it did not alter flowering time or maturity. Soybean (Glycine max (L.) Merr.) is highly demanded by a continuously growing human population. However, increasing soybean yield is a major challenge. FRUITFULL (FUL), a MADS-box transcription factor, plays important roles in multiple developmental processes, especially fruit and pod development, which are crucial for soybean yield formation. However, the functions of its homologs in soybean are not clear. Here, through haplotype analysis, we found that one haplotype of the soybean homolog GmFULa (GmFULa-H02) is dominant in cultivated soybeans, suggesting that GmFULa-H02 was highly selected during domestication and varietal improvement of soybean. Interestingly, transgenic overexpression of GmFULa enhanced vegetative growth with more biomass accumulated and ultimately increased the yield but without affecting the plant height or changing the flowering time and maturity, indicating that it enhances the efficiency of dry matter accumulation. It also promoted the yield factors like branch number, pod number and 100-seed weight, which ultimately increased the yield. It increased the palisade tissue cell number and the chlorophyll content to promote photosynthesis and increase the soluble sugar content in leaves and fresh seeds. Furthermore, GmFULa were found to be sublocalized in the nucleus and positively regulate sucrose synthases (SUSs) and sucrose transporters (SUTs) by binding with the conserved CArG boxes in their promoters. Overall, these results showed GmFULa promotes the capacity of assimilation and the transport of the resultant assimilates to increase yield, and provided insights into the link between GmFULa and sucrose synthesis with transport-related molecular pathways that control seed yield.

Comparative selective signature analysis and high-resolution GWAS reveal a new candidate gene controlling seed weight in soybean

KEY MESSAGE

We detected a QTL qHSW-16 undergone strong selection associated with seed weight and identified a novel candidate gene controlling seed weight candidate gene for this major QTL by qRT-PCT. Soybean [Glycine max (L.) Merr.] provides more than half of the world's oilseed production. To expand its germplasm resources useful for breeding increased yield and oil quality cultivars, it is necessary to resolve the diversity and evolutionary history of this crop. In this work, we resequenced 283 soybean accessions from China and obtained a large number of high-quality SNPs for investigation of the population genetics that underpin variation in seed weight and other agronomic traits. Selective signature analysis detected 78 (~ 25.0 Mb) and 39 (~ 22.60 Mb) novel putative selective signals that were selected during soybean domestication and improvement, respectively. Genome-wide association study (GWAS) identified five loci associated with seed weight. Among these QTLs, qHSW-16, overlapped with the improvement-selective region on chromosome 16, suggesting that this QTL may be underwent strong selection during soybean improvement. Of the 18 candidate genes in qHSW-16, only SoyZH13_16G122400 showed higher expression levels in a large seed variety compared to a small seed variety during seed development. These results identify SoyZH13_16G122400 as a novel candidate gene controlling seed weight and provide foundational insights into the molecular targets for breeding improvement of seed weight and potential seed yield in soybean.

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.