Targeted association mapping demonstrating the complex molecular genetics of fatty acid formation in soybean

Identification and Genetic Dissection of Resistance to Red Crown Rot Disease in a Diverse Soybean Germplasm Population

Red crown rot (RCR) disease caused by Calonectria ilicicola negatively impacts soybean yield and quality. Unfortunately, the knowledge of the genetic architecture of RCR resistance in soybeans is limited. In this study, 299 diverse soybean accessions were used to explore their genetic diversity and resistance to RCR, and to mine for candidate genes via emergence rate (ER), survival rate (SR), and disease severity (DS) by a multi-locus random-SNP-effect mixed linear model of GWAS. All accessions had brown necrotic lesions on the primary root, with five genotypes identified as resistant. Nine single-nucleotide polymorphism (SNP) markers were detected to underlie RCR response (ER, SR, and DS). Two SNPs colocalized with at least two traits to form a haplotype block which possessed nine genes. Based on their annotation and the qRT-PCR, three genes, namely Glyma.08G074600, Glyma.08G074700, and Glyma.12G043600, are suggested to modulate soybean resistance to RCR. The findings from this study could serve as the foundation for breeding RCR-tolerant soybean varieties, and the candidate genes could be validated to deepen our understanding of soybean response to RCR.

Construction of Genetic Map and QTL Mapping for Seed Size and Quality Traits in Soybean (Glycine max L.)

Soybean (Glycine max L.) is the main source of vegetable protein and edible oil for humans, with an average content of about 40% crude protein and 20% crude fat. Soybean yield and quality traits are mostly quantitative traits controlled by multiple genes. The quantitative trait loci (QTL) mapping for yield and quality traits, as well as for the identification of mining-related candidate genes, is of great significance for the molecular breeding and understanding the genetic mechanism. In this study, 186 individual plants of the F2 generation derived from crosses between Changjiangchun 2 and Yushuxian 2 were selected as the mapping population to construct a molecular genetic linkage map. A genetic map containing 445 SSR markers with an average distance of 5.3 cM and a total length of 2375.6 cM was obtained. Based on constructed genetic map, 11 traits including hundred-seed weight (HSW), seed length (SL), seed width (SW), seed length-to-width ratio (SLW), oil content (OIL), protein content (PRO), oleic acid (OA), linoleic acid (LA), linolenic acid (LNA), palmitic acid (PA), stearic acid (SA) of yield and quality were detected by the multiple- d size traits and 113 QTLs related to quality were detected by the multiple QTL model (MQM) mapping method across generations F2, F2:3, F2:4, and F2:5. A total of 71 QTLs related to seed size traits and 113 QTLs related to quality traits were obtained in four generations. With those QTLs, 19 clusters for seed size traits and 20 QTL clusters for quality traits were summarized. Two promising clusters, one related to seed size traits and the other to quality traits, have been identified. The cluster associated with seed size traits spans from position 27876712 to 29009783 on Chromosome 16, while the cluster linked to quality traits spans from position 12575403 to 13875138 on Chromosome 6. Within these intervals, a reference genome of William82 was used for gene searching. A total of 36 candidate genes that may be involved in the regulation of soybean seed size and quality were screened by gene functional annotation and GO enrichment analysis. The results will lay the theoretical and technical foundation for molecularly assisted breeding in soybean.

Genome-Wide Detection of Quantitative Trait Loci and Prediction of Candidate Genes for Seed Sugar Composition in Early Mature Soybean

Seed sugar composition, mainly including fructose, glucose, sucrose, raffinose, and stachyose, is an important indicator of soybean [Glycine max (L.) Merr.] seed quality. However, research on soybean sugar composition is limited. To better understand the genetic architecture underlying the sugar composition in soybean seeds, we conducted a genome-wide association study (GWAS) using a population of 323 soybean germplasm accessions which were grown and evaluated under three different environments. A total of 31,245 single-nucleotide polymorphisms (SNPs) with minor allele frequencies (MAFs) ≥ 5% and missing data ≤ 10% were selected and used in the GWAS. The analysis identified 72 quantitative trait loci (QTLs) associated with individual sugars and 14 with total sugar. Ten candidate genes within the 100 Kb flanking regions of the lead SNPs across six chromosomes were significantly associated with sugar contents. According to GO and KEGG classification, eight genes were involved in the sugar metabolism in soybean and showed similar functions in Arabidopsis. The other two, located in known QTL regions associated with sugar composition, may play a role in sugar metabolism in soybean. This study advances our understanding of the genetic basis of soybean sugar composition and facilitates the identification of genes controlling this trait. The identified candidate genes will help improve seed sugar composition in soybean.

MOTHER-OF-FT-AND-TFL1 regulates the seed oil and protein content in soybean

Soybean is a major crop that produces valuable seed oil and protein for global consumption. Seed oil and protein are regulated by complex quantitative trait loci (QTLs) and have undergone intensive selections during the domestication of soybean. It is essential to identify the major genetic components and understand their mechanism behind seed oil and protein in soybean. We report that MOTHER-OF-FT-AND-TFL1 (GmMFT) is the gene of a classical QTL that has been reported to regulate seed oil and protein content in many studies. Mutation of MFT decreased seeds oil content and weight in both Arabidopsis and soybean, whereas increased expression of GmMFT enhanced seeds oil content and weight. Haplotype analysis showed that GmMFT has undergone selection, which resulted in the extended haplotype homozygosity in the cultivated soybean and the enriching of the oil-favorable allele in modern soybean cultivars. This work unraveled the GmMFT-mediated mechanism regulating seed oil and protein content and seed weight, and revealed a previously unknown function of MFT that provides new insights into targeted soybean improvement and breeding.

Genome-wide association analysis of sucrose and alanine contents in edamame beans

The sucrose and Alanine (Ala) content in edamame beans significantly impacts the sweetness flavor of edamame-derived products as an important attribute to consumers' acceptance. Unlike grain-type soybeans, edamame beans are harvested as fresh beans at the R6 to R7 growth stages when beans are filled 80-90% of the pod capacity. The genetic basis of sucrose and Ala contents in fresh edamame beans may differ from those in dry seeds. To date, there is no report on the genetic basis of sucrose and Ala contents in the edamame beans. In this study, a genome-wide association study was conducted to identify single nucleotide polymorphisms (SNPs) related to sucrose and Ala levels in edamame beans using an association mapping panel of 189 edamame accessions genotyped with a SoySNP50K BeadChip. A total of 43 and 25 SNPs was associated with sucrose content and Ala content in the edamame beans, respectively. Four genes (Glyma.10g270800, Glyma.08g137500, Glyma.10g268500, and Glyma.18g193600) with known effects on the process of sucrose biosynthesis and 37 novel sucrose-related genes were characterized. Three genes (Gm17g070500, Glyma.14g201100 and Glyma.18g269600) with likely relevant effects in regulating Ala content and 22 novel Ala-related genes were identified. In addition, by summarizing the phenotypic data of edamame beans from three locations in two years, three PI accessions (PI 532469, PI 243551, and PI 407748) were selected as the high sucrose and high Ala parental lines for the perspective breeding of sweet edamame varieties. Thus, the beneficial alleles, candidate genes, and selected PI accessions identified in this study will be fundamental to develop edamame varieties with improved consumers' acceptance, and eventually promote edamame production as a specialty crop in the United States.

Identification of QTN-by-environment interactions and their candidate genes for soybean seed oil-related traits using 3VmrMLM

INTRODUCTION

Although seed oil content and its fatty acid compositions in soybean were affected by environment, QTN-by-environment (QEIs) and gene-by-environment interactions (GEIs) were rarely reported in genome-wide association studies.

METHODS

The 3VmrMLM method was used to associate the trait phenotypes, measured in five to seven environments, of 286 soybean accessions with 106,013 SNPs for detecting QTNs and QEIs.

RESULTS

Seven oil metabolism genes (GmSACPD-A, GmSACPD-B, GmbZIP123, GmSWEET39, GmFATB1A, GmDGAT2D, and GmDGAT1B) around 598 QTNs and one oil metabolism gene GmFATB2B around 54 QEIs were verified in previous studies; 76 candidate genes and 66 candidate GEIs were predicted to be associated with these traits, in which 5 genes around QEIs were verified in other species to participate in oil metabolism, and had differential expression across environments. These genes were found to be related to soybean seed oil content in haplotype analysis. In addition, most candidate GEIs were co-expressed with drought response genes in co-expression network, and three KEGG pathways which respond to drought were enriched under drought stress rather than control condition; six candidate genes were hub genes in the co-expression networks under drought stress.

DISCUSSION

The above results indicated that GEIs, together with drought response genes in co-expression network, may respond to drought, and play important roles in regulating seed oil-related traits together with oil metabolism genes. These results provide important information for genetic basis, molecular mechanisms, and soybean breeding for seed oil-related traits.

QTL Mapping of Palmitic Acid Content Using Specific-Locus Amplified Fragment Sequencing (SLAF-Seq) Genotyping in Soybeans (Glycine max L.)

Soybeans are essential crops that supply protein and oil. The composition and contents of soybean fatty acids are relevant to human health and have a significant relationship with soybean oil processing and applications. Identifying quantitative trait locus (QTL) genes related to palmitic acid could facilitate the development of a range of nutritive soybean cultivars using molecular marker-assisted selection. In this study, we used a cultivar with higher palmitic acid content, ‘Dongnong42’, and a lower palmitic acid content cultivar, ‘Hobbit’, to establish F2:6 recombinant inbred lines. A high-density genetic map containing 9980 SLAF markers was constructed and distributed across 20 soybean chromosomes. The genetic map contained a total genetic distance of 2602.58 cM and an average genetic distance of 0.39 cM between adjacent markers. Two QTLs related to palmitic acid content were mapped using inclusive composite interval mapping, explaining 4.2–10.1% of the phenotypic variance in three different years and environments, including the QTL included in seed palmitic 7-3, which was validated by developing SSR markers. Based on the SNP/Indel and significant differential expression analyses of Dongnong42 and Hobbit, two genes, Glyma.15g119700 and Glyma.15g119800, were selected as candidate genes. The high-density genetic map, QTLs, and molecular markers will be helpful for the map-based cloning of palmitic acid content genes. These could be used to accelerate breeding for high nutritive value cultivars via molecular marker-assisted breeding.

Identification of Loci Governing Agronomic Traits and Mutation Hotspots via a GBS-Based Genome-Wide Association Study in a Soybean Mutant Diversity Pool

In this study, we performed a genotyping-by-sequencing analysis and a genome-wide association study of a soybean mutant diversity pool previously constructed by gamma irradiation. A GWAS was conducted to detect significant associations between 37,249 SNPs, 11 agronomic traits, and 6 phytochemical traits. In the merged data set, 66 SNPs on 13 chromosomes were highly associated (FDR p < 0.05) with the following 4 agronomic traits: days of flowering (33 SNPs), flower color (16 SNPs), node number (6 SNPs), and seed coat color (11 SNPs). These results are consistent with the findings of earlier studies on other genetic features (e.g., natural accessions and recombinant inbred lines). Therefore, our observations suggest that the genomic changes in the mutants generated by gamma irradiation occurred at the same loci as the mutations in the natural soybean population. These findings are indicative of the existence of mutation hotspots, or the acceleration of genome evolution in response to high doses of radiation. Moreover, this study demonstrated that the integration of GBS and GWAS to investigate a mutant population derived from gamma irradiation is suitable for dissecting the molecular basis of complex traits in soybeans.

Genome-Wide Identification of Candidate Genes Underlying Soluble Sugar Content in Vegetable Soybean (Glycine max L.) via Association and Expression Analysis

Soluble sugar is a major indicator of the intrinsic quality of vegetable soybean [Glycine max (L.) Merr. ]. The improvement of soluble sugar content in soybean is very important due to its healthcare functions for humans. The genetic mechanism of soluble sugar in soybean is unclear. In this study, 278 diverse soybean accessions were utilized to identify the quantitative trait nucleotides (QTNs) for total soluble sugar content in soybean seeds based on a genome-wide association study (GWAS). A total of 25,921 single-nucleotide polymorphisms (SNPs) with minor allele frequencies (MAFs) ≥ 5% and missing data ≤ 10% were selected for GWAS. Totally, thirteen QTNs associated with total soluble sugar content were identified, which were distributed on ten chromosomes. One hundred and fifteen genes near the 200-kb flanking region of these identified QTNs were considered candidate genes associated with total soluble sugar content in soybean seed. Gene-based association analysis and haplotype analysis were utilized to further identify the effect of candidate genes on total soluble sugar content. Totally, 84 SNPs from seventeen genes across four chromosomes were significantly associated with the total soluble sugar content. Among them, three SNPs from Glyma.02G292900 were identified at two locations, and other eighty-one SNPs from sixteen genes were detected at three locations. Furthermore, expression level analysis of candidate genes revealed that Glyma.02G293200 and Glyma.02G294900 were significantly positively associated with soluble sugar content and Glyma.02G294000 was significantly negatively associated with soluble sugar content. Six genes (i.e., Glyma.02G292600, Glyma.02G292700, Glyma.02G294000, Glyma.02G294300, Glyma.02G294400, and Glyma.15G264200) identified by GWAS were also detected by the analysis of differential expression genes based on soybean germplasms with higher and lower soluble sugar content. Among them, Glyma.02G294000 is the only gene that was identified by gene-based association analysis with total soluble sugar content and was considered an important candidate gene for soluble sugar content. These candidate genes and beneficial alleles would be useful for improving the soluble sugar content of soybean.

Identification of quantitative trait loci associated with seed quality traits between Canadian and Ukrainian mega-environments using genome-wide association study

KEY MESSAGE

Identifying QTL associated with soybean seed quality traits from a diverse GWAS panel cultivated in Canadian and Ukrainian mega-environments may facilitate future cultivar development for foreign markets. Understanding the complex genetic basis of seed quality traits for soybean in the mega-environments (MEs) is critical for developing a marker-assisted selection program that will lead to breeding superior cultivars adapted to specific regions. This study aimed to analyze the accumulation of 14 soybean seed quality traits in Canadian ME and two seed quality traits in Ukrainian ME and identify associated ME specific quantitative trait loci (QTLSP) and ME universal QTL (QTLU) for protein and oil using a genome-wide association study (GWAS) panel consisting of 184 soybean genotypes. The panel was planted in three locations in Canada and two locations in Ukraine in 2018 and 2019. Genotype plus genotype-by-environment biplot analysis was conducted to assess the accumulation of individual seed compounds across different locations. The protein accumulation was high in the Canadian ME and low in the Ukrainian ME, whereas the oil concentration showed the opposite trends between the two MEs. No QTLU were identified across the MEs for protein and oil concentrations. In contrast, nine Canadian QTLSP for protein were identified on various chromosomes, which were co-located with QTL controlling other traits identified in the Canadian ME. The lack of common QTLU for protein and oil suggests that it may be necessary to use QTLSP associated with these traits separately for the Canadian and Ukrainian ME. Additional Ukrainian data for seed compounds other than oil and protein are required to identify novel QTLSP and QTLU for such traits for the individual or combined Canadian and Ukrainian MEs.

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.

Integrating omics approaches to discover and prioritize candidate genes involved in oil biosynthesis in soybean

Soybean is a major source of edible protein and oil. Oil content is a quantitative trait that is significantly determined by genetic and environmental factors. Over the past 30 years, a large volume of soybean genetic, genomic, and transcriptomic data have been accumulated. Nevertheless, integrative analyses of such data remain scarce, in spite of their importance for crop improvement. We hypothesized that the co-occurrence of genomic regions for oil-related traits in different studies may reveal more stable regions encompassing important genetic determinants of oil content and quality in soybean. We integrated publicly available data, obtained with distinct techniques, to discover and prioritize candidate genes involved in oil biosynthesis and regulation in soybean. We detected key fatty acid biosynthesis genes (e.g., BCCP2 and ACCase, FADs, KAS family proteins) and several transcription factors, which are likely regulators of oil biosynthesis. In addition, we identified new candidates for seed oil accumulation and quality, such as Glyma.03G213300 and Glyma.19G160700, which encode a translocator protein homolog and a histone acetyltransferase, respectively. Further, oil and protein genomic hotspots are strongly associated with breeding and not with domestication, suggesting that soybean domestication prioritized other traits. The genes identified here are promising targets for breeding programs and for the development of soybean lines with increased oil content and quality.

Genome-Wide Association Studies of Soybean Yield-Related Hyperspectral Reflectance Bands Using Machine Learning-Mediated Data Integration Methods

In conjunction with big data analysis methods, plant omics technologies have provided scientists with cost-effective and promising tools for discovering genetic architectures of complex agronomic traits using large breeding populations. In recent years, there has been significant progress in plant phenomics and genomics approaches for generating reliable large datasets. However, selecting an appropriate data integration and analysis method to improve the efficiency of phenome-phenome and phenome-genome association studies is still a bottleneck. This study proposes a hyperspectral wide association study (HypWAS) approach as a phenome-phenome association analysis through a hierarchical data integration strategy to estimate the prediction power of hyperspectral reflectance bands in predicting soybean seed yield. Using HypWAS, five important hyperspectral reflectance bands in visible, red-edge, and near-infrared regions were identified significantly associated with seed yield. The phenome-genome association analysis of each tested hyperspectral reflectance band was performed using two conventional genome-wide association studies (GWAS) methods and a machine learning mediated GWAS based on the support vector regression (SVR) method. Using SVR-mediated GWAS, more relevant QTL with the physiological background of the tested hyperspectral reflectance bands were detected, supported by the functional annotation of candidate gene analyses. The results of this study have indicated the advantages of using hierarchical data integration strategy and advanced mathematical methods coupled with phenome-phenome and phenome-genome association analyses for a better understanding of the biology and genetic backgrounds of hyperspectral reflectance bands affecting soybean yield formation. The identified yield-related hyperspectral reflectance bands using HypWAS can be used as indirect selection criteria for selecting superior genotypes with improved yield genetic gains in large breeding populations.

Identification of soybean phosphorous efficiency QTLs and genes using chlorophyll fluorescence parameters through GWAS and RNA-seq

MAIN CONCLUSION

Soybean phosphorous efficiency QTLs were identified and candidate genes were predicted using chlorophyll fluorescence parameters through GWAS and RNA-seq. Phosphorus (P) is an essential nutrient element for crop growth and development, lack of P uptake seriously affects yield in various crops. Photosynthesis is the basis of crop production, while it is very sensitive to P deficiency. It is of great importance to study the genetic relationship between photosynthesis and P efficiency to provide genetic insight for soybean improvement. In this study, a genome-wide association study (GWAS) was performed using 292,035 SNPs and the ratios of four main chlorophyll fluorescence parameters of 219 diverse soybean accessions under P deficiency and normal P across three experiments. In total, 52 SNPs in 12 genomic regions were detected in association with the four main chlorophyll fluorescence parameters under sufficient or deficient P levels. Combined it with RNA-seq analysis, we predicted three candidate genes for the significant genomic regions. For example, the expression level of the candidate gene (Glyma.18g092900) in P deficiency tolerant accession was three times higher than that of P deficiency sensitive one under phosphorous deficiency condition. This study provides insight into genetic links between photosynthetic and phosphorous efficiency and further functional analysis will provide valuable information for understanding the underlying genetic mechanism to facilitate marker-assisted breeding in soybean.

Identification of glutathione transferase gene associated with partial resistance to Sclerotinia stem rot of soybean using genome-wide association and linkage mapping

KEY MESSAGE

Association and linkage mapping techniques were used to identify and verify single nucleotide polymorphisms (SNPs) associated with Sclerotinia sclerotiorum resistance. A novel resistant gene, GmGST , was cloned and shown to be involved in soybean resistance to SSR. Sclerotinia stem rot (SSR), caused by the fungus Sclerotinia sclerotiorum, is one of the most devastating diseases in soybean (Glycine max (Linn.) Merr.) However, the genetic architecture underlying soybean resistance to SSR is poorly understood, despite several mapping and gene mining studies. In the present study, the identification of quantitative trait loci (QTLs) involved in the resistance to S. sclerotiorum was conducted in two segregating populations: an association population that consisted of 261 diverse soybean germplasms, and the MH population, derived from a cross between a partially resistant cultivar (Maple arrow) and a susceptible cultivar (Hefeng25). Three and five genomic regions affecting resistance were detected by genome-wide association study to control the lesion length of stems (LLS) and the death rate of seedling (DRS), respectively. Four QTLs were detected to underlie LLS, and one QTL controlled DRS after SSR infection. A major locus on chromosome (Chr.) 13 (qDRS13-1), which affected both DRS and LLS, was detected in both the natural population and the MH population. GmGST, encoding a glutathione S-transferase, was cloned as a candidate gene in qDRS13-1. GmGST was upregulated by the induction of the partially resistant cultivar Maple arrow. Transgenic experiments showed that the overexpression of GmGST in soybean increased resistance to S. sclerotiorum and the content of soluble pigment in stems of soybean. The results increase our understanding of the genetic architecture of soybean resistance to SSR and provide a framework for the future marker-assisted breeding of resistant soybean cultivars.

Novel Genomic Regions of Fusarium Wilt Resistance in Bottle Gourd [Lagenaria siceraria (Mol.) Standl.] Discovered in Genome-Wide Association Study

Fusarium wilt (FW) is a typical soil-borne disease that seriously affects the yield and fruit quality of bottle gourd. Thus, to improve resistance to FW in bottle gourd, the genetic mechanism underlying FW resistance needs to be explored. In this study, we performed a genome-wide association study (GWAS) based on 5,330 single-nucleotide polymorphisms (SNPs) and 89 bottle gourd accessions. The GWAS results revealed a total of 10 SNPs (P ≤ 0.01, -log₁₀ P ≥ 2.0) significantly associated with FW resistance that were detected in at least two environments (2019DI, 2020DI, and the average across the 2 years); these SNPs were located on chromosomes 1, 2, 3, 4, 8, and 9. Linkage disequilibrium (LD) block structure analysis predicted three potential candidate genes for FW resistance. Genes HG_GLEAN_10001030 and HG_GLEAN_10001042 were within the range of the mean LD block of the marker BGReSe_14202; gene HG_GLEAN_10011803 was 280 kb upstream of the marker BGReSe_00818. Real-time quantitative PCR (qRT-PCR) analysis showed that HG_GLEAN_10011803 was significantly up-regulated in FW-infected plants of YD-4, Yin-10, and Hanbi; HG_GLEAN_10001030 and HG_GLEAN_10001042 were specifically up-regulated in FW-infected plants of YD-4. Therefore, gene HG_GLEAN_10011803 is likely the major effect candidate gene for resistance against FW in bottle gourd. This work provides scientific evidence for the exploration of candidate gene and development of functional markers in FW-resistant bottle gourd breeding programs.

Genetic dissection of soybean partial resistance to sclerotinia stem rot through genome wide association study and high throughout single nucleotide polymorphisms

Sclerotinia stem rot (SSR) is a disease of soybean [Glycine max (L.) Merr] that causes severe yield losses. We studied 185 representative soybean accessions to evaluate partial SSR resistance and sequenced these by the specific-locus amplified fragment sequencing method. In total, 22,048 single-nucleotide polymorphisms (SNPs), with minor allele frequencies (MAF) ≥5% and missing data <3%, were developed and applied to genome-wide association study of SSR responsiveness and assess linkage disequilibrium (LD) level for candidate gene selection. We identified 18 association signals related to SSR partial resistance. Among them, six overlapped the regions of previous quantitative trait loci, and twelve were novel. We identified 243 candidate genes located in the 200 kb genomic region of these peak SNPs. Based on quantitative real-time polymerase chain reaction and haplotype analysis, Glyma.03G196000 and Glyma.20G095100, encoding pentatricopeptide repeat proteins, might be important factors in the resistance response of soybean to SSR.

Identification of SNP markers associated with soybean fatty acids contents by genome-wide association analyses

Composition of fatty acids (FAs) in soybean seed is important for the quality and uses of soybean oil. Using gas chromatography, we have measured soybean FAs profiles of 621 soybean accessions (maturity groups I through IV) grown in five different environments; Columbus, OH (2015), Wooster, OH (2014 and 2015), Plymouth, NC (2015), and Urbana, IL (2015). Using publicly available SoySNP50K genotypic data and the FA profiles from this study, a genome-wide association analysis was completed with a compressed mixed linear model to identify 43 genomic regions significantly associated with a fatty acid at a genome wide significance threshold of 5%. Among these regions, one and three novel genomic regions associated with palmitic acid and stearic acid, respectively, were identified across all five environments. Additionally, nine novel environment-specific FA-related genomic regions were discovered providing new insights into the genetics of soybean FAs. Previously reported FA-related loci, such as FATB1a, SACPD-C, and KASIII, were also confirmed in this study. Our results will be useful for future functional studies and marker-assisted breeding for soybean FAs.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01216-1.

Proteome-wide Systems Genetics to Identify Functional Regulators of Complex Traits

Proteomic technologies now enable the rapid quantification of thousands of proteins across genetically diverse samples. Integration of these data with systems-genetics analyses is a powerful approach to identify new regulators of economically important or disease-relevant phenotypes in various populations. In this review, we summarize the latest proteomic technologies and discuss technical challenges for their use in population studies. We demonstrate how the analysis of correlation structure and loci mapping can be used to identify genetic factors regulating functional protein networks and complex traits. Finally, we provide an extensive summary of the use of proteome-wide systems genetics throughout fungi, plant, and animal kingdoms and discuss the power of this approach to identify candidate regulators and drug targets in large human consortium studies.

Genome-wide association analysis of sucrose concentration in soybean (Glycine max L.) seed based on high-throughput sequencing

The sucrose concentration in soybean seed significantly affects the flavor of soybean-derived products. In this study, an association panel of 178 elite accessions and 33,149 single-nucleotide polymorphisms (SNPs) was utilized to identify quantitative trait nucleotides (QTNs) of sucrose concentration in soybean seeds by genome-wide association study (GWAS). Five QTNs (rs2688589, rs29026218, rs5926884, rs6886889, and rs10299216) distributed across five genomic regions in five chromosomes were identified in two or more locations by GWAS. A total of 60 candidate genes near the 200-kb flanking region of these five identified loci were identified. Three of these genes (Glyma.04G032600, Glyma.04G034600, and Glyma.11G092100) have been reported to be involved in the process of sugar biosynthesis. Based on gene-based association and haplotype analyses, a total of 35 SNPs from 10 genes associated with sucrose concentration were identified. Of them, Glyma.04G032600 was the only gene that has been reported to be related to sucrose content; the other nine genes were novel and may be associated with sucrose content. These beneficial alleles and candidate genes may be of great value in improving sucrose content in soybean seeds.

Towards model-driven characterization and manipulation of plant lipid metabolism

Plant lipids have versatile applications and provide essential fatty acids in human diet. Therefore, there has been a growing interest to better characterize the genetic basis, regulatory networks, and metabolic pathways that shape lipid quantity and composition. Addressing these issues is challenging due to context-specificity of lipid metabolism integrating environmental, developmental, and tissue-specific cues. Here we systematically review the known metabolic pathways and regulatory interactions that modulate the levels of storage lipids in oilseeds. We argue that the current understanding of lipid metabolism provides the basis for its study in the context of genome-wide plant metabolic networks with the help of approaches from constraint-based modeling and metabolic flux analysis. The focus is on providing a comprehensive summary of the state-of-the-art of modeling plant lipid metabolic pathways, which we then contrast with the existing modeling efforts in yeast and microalgae. We then point out the gaps in knowledge of lipid metabolism, and enumerate the recent advances of using genome-wide association and quantitative trait loci mapping studies to unravel the genetic regulations of lipid metabolism. Finally, we offer a perspective on how advances in the constraint-based modeling framework can propel further characterization of plant lipid metabolism and its rational manipulation.

Genome-wide association study of partial resistance to sclerotinia stem rot of cultivated soybean based on the detached leaf method

Sclerotinia stem rot (SSR) is a devastating fungal disease that causes severe yield losses of soybean worldwide. In the present study, a representative population of 185 soybean accessions was selected and utilized to identify the quantitative trait nucleotide (QTN) of partial resistance to soybean SSR via a genome-wide association study (GWAS). A total of 22,048 single-nucleotide polymorphisms (SNPs) with minor allele frequencies (MAF) > 5% and missing data < 3% were used to assess linkage disequilibrium (LD) levels. Association signals associated with SSR partial resistance were identified by two models, including compressed mixed linear model (CMLM) and multi-locus random-SNP-effect mixed linear model (mrMLM). Finally, seven QTNs with major effects (a known locus and six novel loci) via CMLM and nine novel QTNs with minor effects via mrMLM were detected in relation to partial resistance to SSR, respectively. One of all the novel loci (Gm05:14834789 on Chr.05), which was co-located by these two methods, might be a stable one that showed high significance in SSR partial resistance. Additionally, a total of 71 major and 85 minor candidate genes located in the 200-kb genomic region of each peak SNP detected by CMLM and mrMLM were found, respectively. By using a gene-based association, a total of six SNPs from three major effects genes and eight SNPs from four minor effects genes were identified. Of them, Glyma.18G012200 has been characterized as a significant element in controlling fungal disease in plants.

Quantitative trait loci analysis of seed oil content and composition of wild and cultivated soybean

BACKGROUND

Soybean oil is a major source of edible oil, and the domestication of wild soybean has resulted in significant changes in oil content and composition. Extensive efforts have been made to identify genetic loci that are related to soybean oil traits. The objective of this study was to identify quantitative trait loci (QTLs) related to soybean seed oil and compare the fatty acid composition between wild and cultivated soybean.

RESULTS

Using the specific-locus amplified fragment sequencing (SLAF-seq) method, a total of 181 recombinant inbred lines (RILs) derived from a cross between wild soybean ZYD00463 (Glycine soja) and cultivated soybean WDD01514 (Glycine max) were genotyped. Finally, a high-density genetic linkage map comprising 11,398 single-nucleotide polymorphism (SNP) markers on 20 linkage groups (LGs) was constructed. Twenty-four stable QTLs for seed oil content and composition were identified by model-based composite interval mapping (CIM) across multiple environments. Among these QTLs, 23 overlapped with or were adjacent to previously reported QTLs. One QTL, qPA10_1 (5.94-9.98 Mb) on Chr. Ten is a novel locus for palmitic acid. In the intervals of stable QTLs, some interesting genes involved in lipid metabolism were detected.

CONCLUSIONS

We developed 181 RILs from a cross between wild soybean ZYD00463 and cultivated soybean WDD01514 and constructed a high-density genetic map using the SLAF-seq method. We identified 24 stable QTLs for seed oil content and compositions, which includes qPA10_1 on Chr. 10, a novel locus for palmitic acid. Some interesting genes in the QTL regions were also detected. Our study will provide useful information for scientists to learn about genetic variations in lipid metabolism between wild and cultivated soybean.

Identifying Vitamin E Biosynthesis Genes in Elaeis guineensis by Genome-Wide Association Study

Elaeis guineensis is a tropical oil crop and has the highest oil yield per unit area. Palm oil has high palmitic acid content and is also rich in vitamins, including vitamin E. We conducted genome-wide association studies in a diversity panel of 161 E. guineensis accessions to identify single-nucleotide polymorphisms (SNPs) linked with vitamin E and validated candidate genes in these marker-associated intervals. Based on the SNPs reported in our previous research, 47 SNP markers were detected to be significantly associated with the variation of tocopherol and tocotrienol content at a cutoff P value of 6.3 × 10^-7. A total of 656 candidate genes in the flanking regions of the 47 SNPs were identified, followed by pathway enrichment analysis. Of these candidate genes, EgHGGT (homogentisate geranylgeranyl transferase) involved in the biosynthesis of tocotrienols had a higher expression level in the mesocarp compared to other tissues. Expression of the EgHGGT gene was positively correlated with the variation in α-tocotrienol content. Induced overexpression of the gene in Arabidopsis caused a significant increase in vitamin E content and production of α-tocotrienols compared to wild Arabidopsis.

Natural Variation in Fatty Acid Composition of Diverse World Soybean Germplasms Grown in China

Soybean (Glycine max L. Merr.) is one of the most important crops in the world. Its major content of vegetable oil made it widely used for human consumption and several food industries. To investigate the variation in seed fatty acid composition of soybeans from different origins, a set of 633 soybean accessions originated from four diverse germplasm collections—including China, United States of America (USA), Japan, and Russia—were grown in three locations, Beijing, Anhui, and Hainan for two years. The results showed significant differences (P < 0.001) among the four germplasm origins for all fatty acid contents investigated. Higher levels, on average, of palmitic acid (PA) and linolenic acid (LNA) were observed in Russian germplasm (12.31% and 8.15%, respectively), whereas higher levels of stearic acid (SA) and oleic acid (OA) were observed in Chinese germplasm (3.95% and 21.95%, respectively). The highest level of linoleic acid (LA) was noticed in the USA germplasm accessions (56.34%). The largest variation in fatty acid composition was found in LNA, while a large variation was observed between Chinese and USA germplasms for LA level. Maturity group (MG) significantly (P < 0.0001) affected all fatty acids and higher levels of PA, SA, and OA were observed in early maturing accessions, while higher levels of LA and LNA were observed in late maturing accessions. The trends of fatty acids concentrations with different MG in this study further provide an evidence of the importance of MG in breeding for such soybean seed components. Collectively, the unique accessions identified in this study can be used to strengthen the soybean breeding programs for meeting various human nutrition patterns around the globe.

Genome-wide association study of inflorescence length of cultivated soybean based on the high-throughout single-nucleotide markers

As an important and complex trait, inflorescence length (IL) of soybean [Glycine max (L.) Merr.] significantly affected seed yields. Therefore, elucidating molecular basis of inflorescence architecture, especially for IL, was important for improving soybean yield potentials. Longer IL meaned to have more pod and seed in soybean. Hence, increasing IL and improving yield are targets for soybean breeding. In this study, a association panel, comprising 283 diverse samples, was used to dissect the genetic basis of IL based on genome-wide association analysis (GWAS) and haplotype analysis. GWAS and haplotype analysis were conducted through high-throughout single-nucleotide polymorphisms (SNP) developed by SLAF-seq methodology. A total of 39, 057 SNPs (minor allele frequency ≥ 0.2 and missing data ≤ 10%) were utilized to evaluate linkage disequilibrium (LD) level in the tested association panel. A total of 30 association signals were identified to be associated with IL via GWAS. Among them, 13 SNPs were novel, and another 17 SNPs were overlapped or located near the linked regions of known quantitative trait nucleotide (QTN) with soybean seed yield or yield component. The functional genes, located in the 200-kb genomic region of each peak SNP, were considered as candidate genes, such as the cell division/ elongation, specific enzymes, and signaling or transport of specific proteins. These genes have been reported to participant in the regulation of IL. Ten typical long-IL lines and ten typical short-IL lines were re-sequencing, and then, six SNPs from five genes were obtained based on candidate gene-based association. In addition, 42 haplotypes were defined based on haplotype analysis. Of them, 11 haplotypes were found to regulate long IL (> 14 mm) in soybean. The identified 30 QTN with beneficial alleles and their candidate genes might be valuable for dissecting the molecular mechanisms of IL and further improving the yield potential of soybean.

Genome-wide association and transcriptional studies reveal novel genes for unsaturated fatty acid synthesis in a panel of soybean accessions

BACKGROUND

The nutritional value of soybean oil is largely influenced by the proportions of unsaturated fatty acids (FAs), including oleic acid (OA, 18:1), linoleic acid (LLA, 18:2), and linolenic acid (LNA, 18:3). Genome-wide association (GWAS) studies along with gene expression studies in soybean [Glycine max (L.) Merr.] were leveraged to dissect the genetics of unsaturated FAs.

RESULTS

A association panel of 194 diverse soybean accessions were phenotyped in 2013, 2014 and 2015 to identify Single Nucleotide Polymorphisms (SNPs) associated with OA, LLA, and LNA content, and determine putative candidate genes responsible for regulating unsaturated FAs composition. 149 SNPs that represented 73 genomic regions were found to be associated with the unsaturated FA contents in soybean seeds according to the results of GWAS. Twelve novel genes were predicted to be involved in unsaturated FA synthesis in soybean. The relationship between expression pattern of the candidate genes and the accumulation of unsaturated FAs revealed that multiple genes might be involved in unsaturated FAs regulation simultaneously but work in very different ways: Glyma.07G046200 and Glyma.20G245500 promote the OA accumulation in soybean seed in all the tested accessions; Glyma.13G68600 and Glyma.16G200200 promote the OA accumulation only in high OA germplasms; Glyma.07G151300 promotes OA accumulation in higher OA germplasms and suppresses that in lower OA germplasms; Glyma.16G003500 has the effect of increasing LLA accumulation in higher LA germplasms; Glyma.07G254500 suppresses the accumulation of LNA in lower OA germplasms; Glyma.14G194300 might be involved in the accumulation of LNA content in lower LNA germplasms.

CONCLUSIONS

The beneficial alleles and candidate genes identified might be valuable for improving marker-assisted breeding efficiency and exploring the molecular mechanisms underlying unsaturated fatty acid of soybean.

Identification of the Genomic Region Underlying Seed Weight per Plant in Soybean (Glycine max L. Merr.) via High-Throughput Single-Nucleotide Polymorphisms and a Genome-Wide Association Study

Seed weight per plant (SWPP) of soybean (Glycine max (L.) Merr.), a complicated quantitative trait controlled by multiple genes, was positively associated with soybean seed yields. In the present study, a natural soybean population containing 185 diverse accessions primarily from China was used to analyze the genetic basis of SWPP via genome-wide association analysis (GWAS) based on high-throughput single-nucleotide polymorphisms (SNPs) generated by the Specific Locus Amplified Fragment Sequencing (SLAF-seq) method. A total of 33,149 SNPs were finally identified with minor allele frequencies (MAF) > 5% which were present in 97% of all the genotypes. Twenty association signals associated with SWPP were detected via GWAS. Among these signals, eight SNPs were novel loci, and the other twelve SNPs were overlapped or located in the linked genomic regions of the reported QTL from SoyBase database. Several genes belonging to the categories of hormone pathways, RNA regulation of transcription in plant development, ubiquitin, transporting systems, and other metabolisms were considered as candidate genes associated with SWPP. Furthermore, nine genes from the flanking region of Gm07:19488264, Gm08:15768591, Gm08:15768603, or Gm18:23052511 were significantly associated with SWPP and were stable among multiple environments. Nine out of 18 haplotypes from nine genes showed the effect of increasing SWPP. The identified loci along with the beneficial alleles and candidate genes could be of great value for studying the molecular mechanisms underlying SWPP and for improving the potential seed yield of soybean in the future.

Detecting the QTL-Allele System of Seed Oil Traits Using Multi-Locus Genome-Wide Association Analysis for Population Characterization and Optimal Cross Prediction in Soybean

Soybean is one of the world's major vegetative oil sources, while oleic acid and linolenic acid content are the major quality traits of soybean oil. The restricted two-stage multi-locus genome-wide association analysis (RTM-GWAS), characterized with error and false-positive control, has provided a potential approach for a relatively thorough detection of whole-genome QTL-alleles. The Chinese soybean landrace population (CSLRP) composed of 366 accessions was tested under four environments to identify the QTL-allele constitution of seed oil, oleic acid and linolenic acid content (SOC, OAC, and LAC). Using RTM-GWAS with 29,119 SNPLDBs (SNP linkage disequilibrium blocks) as genomic markers, 50, 98, and 50 QTLs with 136, 283, and 154 alleles (2-9 per locus) were detected, with their contribution 82.52, 90.31, and 83.86% to phenotypic variance, corresponding to their heritability 91.29, 90.97, and 90.24% for SOC, OAC, and LAC, respectively. The RTM-GWAS was shown to be more powerful and efficient than previous single-locus model GWAS procedures. For each trait, the detected QTL-alleles were organized into a QTL-allele matrix as the population genetic constitution. From which the genetic differentiation among 6 eco-populations was characterized as significant allele frequency differentiation on 28, 56, and 30 loci for the three traits, respectively. The QTL-allele matrices were also used for genomic selection for optimal crosses, which predicted transgressive potential up to 24.76, 40.30, and 2.37% for the respective traits, respectively. From the detected major QTLs, 38, 27, and 25 candidate genes were annotated for the respective traits, and two common QTL covering eight genes were identified for further study.

Genome-wide association mapping of QTL underlying seed oil and protein contents of a diverse panel of soybean accessions

To investigate the genetic basis of variation in oil and protein contents in soybean seeds, a diverse collection of 421 mainly Chinese soybean cultivars was genotyped using 1536 SNPs, mostly from candidate genes related to acyl-lipid metabolism and from regions harboring known QTL. Six significant associations were identified for each of seed oil and protein contents which individually explained 2.7-5.9% of the phenotypic variance. Six associations occurred in or near known QTL and the remaining are putative novel QTL. Ten significant associations influenced the oil content without decreasing protein content, and vice versa. One SNP was pleiotropic, with opposite effects on oil and protein contents. The genetic region covering Map-6076 and-6077 was shown to be involved in controlling oil content in soybean by integrating the results of association mapping with information on known QTL and tissue-specific expression data. This region was subject to strong selection during the genetic improvement of soybean. Our results not only confirm and refine the map positions of known QTL but also contribute to a further elucidation of the genetic architecture of protein and oil contents in soybean seeds by identifying new associations exhibiting pleiotropic effects on seed protein and oil contents.

Genome-wide association studies dissect the genetic networks underlying agronomical traits in soybean

BACKGROUND

Soybean (Glycine max [L.] Merr.) is one of the most important oil and protein crops. Ever-increasing soybean consumption necessitates the improvement of varieties for more efficient production. However, both correlations among different traits and genetic interactions among genes that affect a single trait pose a challenge to soybean breeding.

RESULTS

To understand the genetic networks underlying phenotypic correlations, we collected 809 soybean accessions worldwide and phenotyped them for two years at three locations for 84 agronomic traits. Genome-wide association studies identified 245 significant genetic loci, among which 95 genetically interacted with other loci. We determined that 14 oil synthesis-related genes are responsible for fatty acid accumulation in soybean and function in line with an additive model. Network analyses demonstrated that 51 traits could be linked through the linkage disequilibrium of 115 associated loci and these links reflect phenotypic correlations. We revealed that 23 loci, including the known Dt1, E2, E1, Ln, Dt2, Fan, and Fap loci, as well as 16 undefined associated loci, have pleiotropic effects on different traits.

CONCLUSIONS

This study provides insights into the genetic correlation among complex traits and will facilitate future soybean functional studies and breeding through molecular design.

High-resolution mapping of QTL for fatty acid composition in soybean using specific-locus amplified fragment sequencing

KEY MESSAGE

We constructed a high-density linkage map comprising 3541 markers developed by specific-locus amplified fragment sequencing, and identified 26 stable QTL including nine novel loci, for fatty acid composition in soybean. Soybean oil quality and stability are mainly determined by the fatty acid composition of the seed. In the present study, we constructed a high-density genetic linkage map using 200 recombinant inbred lines derived from a cross between cultivated soybean varieties Luheidou2 and Nanhuizao, and SNP markers developed by specific-locus amplified fragment sequencing (SLAF-seq). This map comprises 3541 markers on 20 linkage groups and spans a genetic distance of 2534.42 cM, with an average distance of 0.72 cM between adjacent markers. Inclusive composite interval mapping revealed 26 stable QTL for five fatty acids, explaining 0.4-37.0% of the phenotypic variance for individual fatty acids across environments. Of these QTL, nine are novel loci (qLA1, qLNA2_1, qPA4_1, qLA4_1, qPA6_1, qSA12_1, qPA16_1, qOA18_1, and qFA19_1). These stable QTL harbor three fatty acid biosynthesis genes (GmFabG, GmACP, and GmFAD8), and 66 genes encoding lipid-related transcription factors. These stable QTL and tightly linked SNP markers can be used for marker-assisted selection in soybean breeding programs.

Identification of genetic loci associated with crude protein and mineral concentrations in alfalfa (Medicago sativa) using association mapping

BACKGROUND

Alfalfa (Medicago sativa) is one of the most important legume forage species in China and many other countries of the world. It provides a quality source of proteins and minerals to animals. Genetic underpinnings for these important traits, however, are elusive. An alfalfa (M. sativa) association mapping study for six traits, namely crude protein (CP), rumen undegraded protein (RUP), and four mineral elements (Ca, K, Mg and P), was conducted in three consecutive years using a large collection encompassing 336 genotypes genotyped with 85 simple sequence repeat (SSR) markers.

RESULTS

All the traits were significantly influenced by genotype, environment, and genotype × environment interaction. Eight-five significant associations (P < 0.005) were identified. Among these, five associations with Ca were repeatedly observed and six co-localized associations were identified.

CONCLUSIONS

The identified marker alleles significantly associated with the traits provided important information for understanding genetic controls of alfalfa quality. The markers could be used in assisting selection for the individual traits in breeding populations for developing new alfalfa cultivars.

Short communication: Discrimination between retail bovine milks with different fat contents using chemometrics and fatty acid profiling

We used a multivariate chemometric approach to differentiate or associate retail bovine milks with different fat contents and non-dairy beverages, using fatty acid profiles and statistical analysis. We collected samples of bovine milk (whole, semi-skim, and skim; n = 62) and non-dairy beverages (n = 27), and we analyzed them using gas-liquid chromatography. Principal component analysis of the fatty acid data yielded 3 significant principal components, which accounted for 72% of the total variance in the data set. Principal component 1 was related to saturated fatty acids (C4:0, C6:0, C8:0, C12:0, C14:0, C17:0, and C18:0) and monounsaturated fatty acids (C14:1 cis-9, C16:1 cis-9, C17:1 cis-9, and C18:1 trans-11); whole milk samples were clearly differentiated from the rest using this principal component. Principal component 2 differentiated semi-skim milk samples by n-3 fatty acid content (C20:3n-3, C20:5n-3, and C22:6n-3). Principal component 3 was related to C18:2 trans-9,trans-12 and C20:4n-6, and its lower scores were observed in skim milk and non-dairy beverages. A cluster analysis yielded 3 groups: group 1 consisted of only whole milk samples, group 2 was represented mainly by semi-skim milks, and group 3 included skim milk and non-dairy beverages. Overall, the present study showed that a multivariate chemometric approach is a useful tool for differentiating or associating retail bovine milks and non-dairy beverages using their fatty acid profile.

A genome-wide association study of seed composition traits in wild soybean (Glycine soja)

BACKGROUND

Cultivated soybean (Glycine max) is a major agricultural crop that provides a crucial source of edible protein and oil. Decreased amounts of saturated palmitic acid and increased amounts of unsaturated oleic acid in soybean oil are considered optimal for human cardiovascular health and therefore there has considerable interest by breeders in discovering genes affecting the relative concentrations of these fatty acids. Using a genome-wide association (GWA) approach with nearly 30,000 single nucleotide polymorphisms (SNPs), we investigated the genetic basis of protein, oil and all five fatty acid levels in seeds from a sample of 570 wild soybeans (Glycine soja), the progenitor of domesticated soybean, to identify quantitative trait loci (QTLs) affecting these seed composition traits.

RESULTS

We discovered 29 SNPs located on ten different chromosomes that are significantly associated with the seven seed composition traits in our wild soybean sample. Eight SNPs co-localized with QTLs previously uncovered in linkage or association mapping studies conducted with cultivated soybean samples, while the remaining SNPs appeared to be in novel locations. Twenty-four of the SNPs significantly associated with fatty acid variation, with the majority located on chromosomes 14 (6 SNPs) and seven (8 SNPs). Two SNPs were common for two or more fatty acids, suggesting loci with pleiotropic effects. We also identified some candidate genes that are involved in fatty acid metabolism and regulation. For each of the seven traits, most of the SNPs produced differences between the average phenotypic values of the two homozygotes of about one-half standard deviation and contributed over 3% of their total variability.

CONCLUSIONS

This is the first GWA study conducted on seed composition traits solely in wild soybean populations, and a number of QTLs were found that have not been previously discovered. Some of these may be useful to breeders who select for increased protein/oil content or altered fatty acid ratios in the seeds. The results also provide additional insight into the genetic architecture of these traits in a large sample of wild soybean, and suggest some new candidate genes whose molecular effects on these traits need to be further studied.