Joint GWAS and WGCNA Identify Genes Regulating the Isoflavone Content in Soybean Seeds.

Isoflavone is a secondary metabolite of the soybean phenylpropyl biosynthesis pathway with physiological activity and is beneficial to human health. In this study, the isoflavone content of 205 soybean germplasm resources from 3 locations in 2020 showed wide phenotypic variation. A joint genome-wide association study (GWAS) and weighted gene coexpression network analysis (WGCNA) identified 33 single-nucleotide polymorphisms and 11 key genes associated with soybean isoflavone content. Gene ontology enrichment analysis, gene coexpression, and haplotype analysis revealed natural variations in the Glyma.12G109800 (GmOMT7) gene and promoter region, with Hap1 being the elite haplotype. Transient overexpression and knockout of GmOMT7 increased and decreased the isoflavone content, respectively, in hairy roots. The combination of GWAS and WGCNA effectively revealed the genetic basis of soybean isoflavone and identified potential genes affecting isoflavone synthesis and accumulation in soybean, providing a valuable basis for the functional study of soybean isoflavone.

GWAS and WGCNA Analysis Uncover Candidate Genes Associated with Oil Content in Soybean.

Soybean vegetable oil is an important source of the human diet. However, the analysis of the genetic mechanism leading to changes in soybean oil content is still incomplete. In this study, a total of 227 soybean materials were applied and analyzed by a genome-wide association study (GWAS). There are 44 quantitative trait nucleotides (QTNs) that were identified as associated with oil content. A total of six, four, and 34 significant QTN loci were identified in Xiangyang, Hulan, and Acheng, respectively. Of those, 26 QTNs overlapped with or were near the known oil content quantitative trait locus (QTL), and 18 new QTNs related to oil content were identified. A total of 594 genes were located near the peak single nucleotide polymorphism (SNP) from three tested environments. These candidate genes exhibited significant enrichment in tropane, piperidine, and pyridine alkaloid biosynthesiss (ko00960), ABC transporters (ko02010), photosynthesis-antenna proteins (ko00196), and betalain biosynthesis (ko00965). Combined with the GWAS and weighted gene co-expression network analysis (WGCNA), four candidate genes (Glyma.18G300100, Glyma.11G221100, Glyma.13G343300, and Glyma.02G166100) that may regulate oil content were identified. In addition, Glyma.18G300100 was divided into two main haplotypes in the studied accessions. The oil content of haplotype 1 is significantly lower than that of haplotype 2. Our research findings provide a theoretical basis for improving the regulatory mechanism of soybean oil content.

Integrating Genome-Wide Association Study, Transcriptome and Metabolome Reveal Novel QTL and Candidate Genes That Control Protein Content in Soybean.

Protein content (PC) is crucial to the nutritional quality of soybean [Glycine max (L.) Merrill]. In this study, a total of 266 accessions were used to perform a genome-wide association study (GWAS) in three tested environments. A total of 23,131 high-quality SNP markers (MAF ≥ 0.02, missing data ≤ 10%) were identified. A total of 40 association signals were significantly associated with PC. Among them, five novel quantitative trait nucleotides (QTNs) were discovered, and another 32 QTNs were found to be overlapping with the genomic regions of known quantitative trait loci (QTL) related to soybean PC. Combined with GWAS, metabolome and transcriptome sequencing, 59 differentially expressed genes (DEGs) that might control the change in protein content were identified. Meantime, four commonly upregulated differentially abundant metabolites (DAMs) and 29 commonly downregulated DAMs were found. Remarkably, the soybean gene Glyma.08G136900, which is homologous with Arabidopsis hydroxyproline-rich glycoproteins (HRGPs), may play an important role in improving the PC. Additionally, Glyma.08G136900 was divided into two main haplotype in the tested accessions. The PC of haplotype 1 was significantly lower than that of haplotype 2. The results of this study provided insights into the genetic mechanisms regulating protein content in soybean.

GmPLP1 negatively regulates soybean resistance to high light stress by modulating photosynthetic capacity and reactive oxygen species accumulation in a blue light-dependent manner.

High light stress is an important factor limiting crop yield. Light receptors play an important role in the response to high light stress, but their mechanisms are still poorly understood. Here, we found that the abundance of GmPLP1, a positive blue light receptor protein, was significantly inhibited by high light stress and mainly responded to high blue light. GmPLP1 RNA-interference soybean lines exhibited higher light energy utilization ability and less light damage and reactive oxygen species (ROS) accumulation in leaves under high light stress, while the phenotype of GmPLP1:GmPLP1-Flag overexpression soybean showed the opposite characteristics. Then, we identified a protein-protein interaction between GmPLP1 and GmVTC2, and the intensity of this interaction was primarily affected by sensing the intensity of blue light. More importantly, overexpression of GmVTC2b improved soybean tolerance to high light stress by enhancing the ROS scavenging capability through increasing the biosynthesis of ascorbic acid. This regulation was significantly enhanced after interfering with a GmPLP1-interference fragment in GmVTC2b-ox soybean leaves, but was weakened when GmPLP1 was transiently overexpressed. These findings demonstrate that GmPLP1 regulates the photosynthetic capacity and ROS accumulation of soybean to adapt to changes in light intensity by sensing blue light. In summary, this study discovered a new mechanism through which GmPLP1 participates in high light stress in soybean, which has great significance for improving soybean yield and the adaptability of soybean to high light.

Genome-wide identification and analysis of glyceraldehyde-3-phosphate dehydrogenase family reveals the role of GmGAPDH14 to improve salt tolerance in soybean (Glycine max L.).

Introduction: Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) is an essential key enzyme in the glycolytic pathway and plays an important role in stress responses. Although GAPDH family genes have been found in different plant species, the determination of their gene family analysis and their functional roles in soybean are still unknown. Methods: In this study, gene sequence and expression data were obtained using online tools, and systematic evolution, expression profile analysis, and qRT-PCR analysis were conducted. Results and Discussion: Here a total of 16 GmGAPDH genes were identified on nine chromosomes, which were classified into three clusters. Additionally, all GmGAPDH genes harbor two highly conserved domains, including Gp_dh_N (PF00044) and Gp_dh_C (PF02800). The qRTPCR analysis also showed that most GmGAPDH genes significantly responded to multiple abiotic stresses, including NaHCO3, polyethylene glycol, cold, and salt. Among them, GmGAPDH14 was extraordinarily induced by salt stress. The GmGAPDH14 gene was cloned and overexpressed through soybean hair roots. The overexpressed transgenic soybean plants of the GmGAPDH14 gene have also shown better growth than that of control plants. Moreover, the overexpressed transgenic plants of GmGAPDH14 gene had higher activities of superoxide dismutase but lower malonaldehyde (MDA) content than those of control plants under salt stress. Meanwhile, a total of four haplotypes were found for the GmGAPDH14 gene, and haplotypes 2, 3, and 4 were beneficial for the tolerance of soybean to salt stress. These results suggest that the GmGAPDH14 gene might be involved in the process of soybean tolerance to salt stress. The results of this study will be valuable in understanding the role of GAPDH genes in the abiotic stress response of soybean.

Genome-wide association analysis of resistance to frogeye leaf spot China race 7 in soybean based on high-throughput sequencing.

KEY MESSAGE: FLS is a disease that causes severe yield reduction in soybean. In this study, four genes (Glyma.16G176800, Glyma.16G177300, Glyma.16G177400 and Glyma.16G182300) were tentatively confirmed to play an important role in the resistance of soybean to FLS race 7. Frogeye leaf spot (FLS) causes severe yield loss in soybean and has been found in several countries worldwide. Therefore, it is necessary to select and utilize FLS-resistant varieties for the management of FLS. In the present study, 335 representative soybean materials were assessed for partial resistance to FLS race 7. Quantitative trait nucleotide (QTN) and FLS race 7 candidate genes were identified using genome-wide association analysis (GWAS) based on a site-specific amplified fragment sequencing (SLAF-seq) approach. A total of 23,156 single-nucleotide polymorphisms (SNPs) were used to evaluate the level of linkage disequilibrium with a minor allele frequency ≥ 5 and deletion data < 3%. These SNPs covered about 947.01 MBP, nearly 86.09% of the entire soybean genome. In addition, a compressed mixed linear model was utilized to identify association signals for partial resistance to FLS race 7. A total of 15 QTNs associated with resistance were found to be novel for FLS race 7 resistance. A total of 217 candidate genes located in the 200-kb genomic region of these peak SNPs were identified. Based on gene association analysis, qRT-PCR, haplotype analysis and virus-induced gene silencing (VIGS) systems were used to further verify candidate genes Glyma.16G176800, Glyma.16G177300, Glyma.16G177400 and Glyma.16G182300. This indicates that these four candidate genes may participate in FLS race 7 resistance responses.

Hierarchical Poly(3,4-ethylenedioxythiophene):Poly(styrenesulfonate)/Reduced graphene oxide/Polypyrrole hybrid electrode with excellent rate capability and cycling stability for fiber-shaped supercapacitor.

Fiber-shaped supercapacitor (FSSC) is considered as a promising energy storage device for wearable electronics due to its high power density and outstanding safety. However, it is still a great challenge to simultaneously achieve high specific capacitance especially at rapid charging/discharging rate and long-term cycling stability of fiber electrode in FSSC for practical application. Here, a ternary poly(3,4-ethylenedioxythiophene):poly(styrenesulfonate)/reduced graphene oxide/polypyrrole (PEDOT:PSS/rGO/PPy) fiber electrode was constructed by in situ chemical polymerization of pyrrole on hydrothermally-assembled and acid-treated PEDOT:PSS/rGO (PG) hybrid hydrogel fiber. In this case, the porous PG hybrid fiber framework possesses combined advantages of highly-conductive PEDOT and flexible two-dimensional (2D) small-sized rGO sheets, which provides large surface area for the deposition of high-pseudocapacitance PPy, multiscale electrons/ions transport channels for the efficient utilization of active sites, and buffering layers to accommodate the structure change during electrochemical process. Attributed to the synergy, as-obtained ternary fiber electrode presents ultrahigh volumetric/areal specific capacitance (389 F cm-3 at 1 A cm-3 or 983 mF cm-2 at 2.5 mA cm-2) and outstanding rate performance (56 %, 1-20 A cm-3). In addition, 80 % preservation of initial capacitance after 8000 cycles for the corresponding FSSC also illustrates its greatly improved cycle stability compared with 64 % of binary PEDOT:PSS/PPy based counterpart. Accordingly, here proposed strategy promises a new opportunity to develop high-activity and durable electrode materials with potential applications in supercapacitor and beyond.

Genome-Wide Association Analysis and Gene Mining of Resistance to China Race 1 of Frogeye Leaf Spot in Soybean.

Soybean frogeye leaf spot (FLS) is a worldwide fungal disease. Its higher occurrence frequency and wider distribution range always led to severe yield losses of soybean, therefore, breeding new cultivars with FLS resistance has been an important breeding goal for soybean breeders. In this study, an association panel of 183 representative soybean accessions was used to evaluate their resistance to FLS race 1, and to identify quantitative trait nucleotides (QTNs) and candidate genes based on genome-wide association study (GWAS) and high-throughput single-nucleotide polymorphisms (SNPs). A total of 23,156 high-quality SNPs were developed using the specific locus-amplified fragment sequencing (SLAF-seq) approach. Finally, 13 novel association signals associated with FLS race 1 resistance were identified by the compressed mixed linear model (CMLM). In addition, 119 candidate genes were found within the 200-kb flanking genomic region of these 13 peak SNPs. Based on the gene-based association analysis, haplotype analysis, expression pattern analysis, and virus-induced gene silencing (VIGS) systems, four genes (Glyma.05G121100, Glyma.17G228300, Glyma.19G006900, and Glyma.19G008700) were preliminarily proved to play an important role in the soybean resistance to FLS race 1.

Fine-Mapping and Functional Analyses of a Candidate Gene Controlling Isoflavone Content in Soybeans Seed.

Isoflavones, one of the most important secondary metabolites produced by soybeans (Glycine max (L.) Merr.), are important for a variety of biological processes, and are beneficial for human health. To identify genetic loci underlying soybean isoflavone content, a mapping population containing 119 F5:18 recombinant inbred lines, derived by crossing soybean cultivar "Zhongdou27" with "Dongong8004," was used. We identified 15 QTLs associated with isoflavone contents. A novel loci, qISO19-1, was mapped onto soybean chromosome 19 and was fine-mapped to a 62.8 kb region using a BC2F2 population. We considered GmMT1 as a candidate gene for the qISO19-1 locus due to the significant positive correlation recovered between its expression level and isoflavone content in the seeds of 43 soybean germplasms. Overexpression of GmMT1 in Arabidopsis and soybean cultivars increased isoflavone contents. Transgenic soybeans overexpressing GmMT1 also exhibited improved resistance to pathogenic infection, while transgenic Arabidopsis resisted salt and drought stress.

Comparison of optical coherence tomography-guided and intravascular ultrasound-guided rotational atherectomy for calcified coronary lesions.

BACKGROUND: To compare the effect and outcomes of optical coherence tomography (OCT)-guided rotational atherectomy (RA) with intravascular ultrasound (IVUS)-guided RA in the treatment of calcified coronary lesions. METHODS: Data of calcified coronary lesions treated with RA that underwent OCT-guided or IVUS-guided from January 2016 to December 2019 at a single-center registry were retrospectively analyzed. The effect and outcomes between underwent OCT-guided RA and IVUS-guided RA were compared. RESULTS: A total of 33 lesions in 32 patients received OCT-guided RA and 51 lesions in 47 patients received IVUS-guided RA. There was no significant difference between OCT-guided RA group and IVUS-guided RA group in clinical baselines characteristics. Comparing the procedural and lesions characteristics of the two groups, the contrast volume was larger [(348.8 ± 110.6) ml vs. (275.2 ± 76.8) ml, P = 0.002] and the scoring balloon was more frequently performed (33.3% vs. 3.9%, P = 0.001) after RA and before stenting in the OCT-guided RA group. Comparing the intravascular imaging findings of the two groups, stent expansion was significantly larger in the OCT-guided RA group ([82 ± 8]% vs. [75 ± 9]%, P = 0.001). Both groups achieved procedural success immediately. There were no significantly differences in the incidence of complications. Although there was no statistical difference in the occurrence of MACE at 1 year between OCT-guided RA group and IVUS-guided RA group (3.1% vs. 6.4%, P = 0.517), no cardiovascular death, TVR and stent thrombosis occurred in OCT-guided RA group. CONCLUSIONS: OCT-guided RA compared to IVUS-guided RA for treating calcified coronary lesions resulted in better stent expansion and may have improved prognosis.

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.

GmST1, which encodes a sulfotransferase, confers resistance to soybean mosaic virus strains G2 and G3.

Soybean mosaic virus (SMV) is one of the most widespread and devastating viral diseases worldwide. The genetic architecture of qualitative resistance to SMV in soybean remains unclear. Here, the Rsvg2 locus was identified as underlying soybean resistance to SMV by genome-wide association and linkage analyses. Fine mapping results showed that soybean resistance to SMV strains G2 and G3 was controlled by a single dominant gene, GmST1, on chromosome 13, encoding a sulfotransferase (SOT). A key variation at position 506 in the coding region of GmST1 associated with the structure of the encoded SOT and changed SOT activity levels between RSVG2-S and RSVG2-R alleles. In RSVG2-S allele carrier "Hefeng25", the overexpression of GmST1 carrying the RSVG2-R allele from the SMV-resistant line "Dongnong93-046" conferred resistance to SMV strains G2 and G3. Compared to Hefeng25, the accumulation of SMV was decreased in transgenic plants carrying the RSVG2-R allele. SMV infection differentiated both the accumulation of jasmonates and expression patterns of genes involved in jasmonic acid (JA) signalling, biosynthesis and catabolism in RSVG2-R and RSVG2-S allele carriers. This characterization of GmST1 suggests a new scenario explaining soybean resistance to SMV.

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.

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.

Identification of Loci and Candidate Genes Analyses for Tocopherol Concentration of Soybean Seed.

Tocopherol (Toc) occurs in soybean seeds and is extracted together with the soybean oil. Toc is utilized as an antioxidant in food and an additive in animal feed. A total of 180 representative accessions and 144 recombinant inbred lines (RILs) from the cross of 'Hefeng 25' and 'OAC Bayfield' were selected to evaluate individuals and total Toc concentrations in soybean seeds. The 180 soybean samples were sequenced by the approach of Specific Locus Amplified Fragment Sequencing (SLAF-seq). A total of 22,611 single nucleotide polymorphisms (SNPs) were developed. Nineteen quantitative trait nucleotides (QTNs) were identified associated with individual or total-Toc based on genome-wide association analysis (GWAS). Among them, three QTNs located near known QTLs, and 16 were novel. Eighteen QTLs and nine eQTLs were also detected by linkage mapping. The QTN rs9337368 on Chr.02 was colocalized according to the linkage mapping of the RILs and genome-wide association analysis and regarded as a stable locus for mining the candidate genes in association with Toc. A total of 42 candidate genes near the 200 kbp flanking region of this identified locus were found. Upon a gene-based association, 11 SNPs from five genes out of the 42 candidates were detected. Expression level analysis of five candidate genes revealed that two genes were significantly related to Toc content. The identified loci, along with the candidate genes, might be valuable for increasing the Toc concentration in soybean seeds and improving the nutritional value of soybean oil.

Identification of a candidate gene associated with isoflavone content in soybean seeds using genome-wide association and linkage mapping.

Isoflavone, a secondary metabolite produced by Glycine max (L.) Merr. (soybean), is valuable for human and plant health. The genetic architecture of soybean isoflavone content remains unclear, however, despite several mapping studies. We generated genomic data for 200 soybean cultivars and 150 recombinant inbred lines (RILs) to localize putative loci associated with soybean seed isoflavone content. Using a genome-wide association study (GWAS), we identified 87 single-nucleotide polymorphisms (SNPs) that were significantly associated with isoflavone concentration. Using linkage mapping, we identified 37 quantitative trait loci (QTLs) underlying the content of four isoflavones found in the RILs. A major locus on chromosome 8 (qISO8-1) was co-located by both the GWAS and linkage mapping. qISO8-1 was fine mapped to a 99.5-kb region, flanked by SSR_08_1651 and SSR_08_1656, in a BC2 F5 population. GmMPK1, encoding a mitogen-activated protein kinase, was identified as the causal gene in qISO8-1, and two natural GmMPK1 polymorphisms were significantly associated with isoflavone content. Overexpression of GmMPK1 in soybean hairy roots resulted in increased isoflavone concentrations. Overexpressing GmMPK1 in transgenic soybeans had greater resistance to Phytophthora root rot, suggesting that GmMPK1 might increase soybean resistance to biotic stress by influencing isoflavone content. Our results not only increase our understanding of the genetic architecture of soybean seed isoflavone content, but also provide a framework for the future marker-assisted breeding of high isoflavone content in soybean cultivars.

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.

Genome wide association mapping and candidate gene analysis for hundred seed weight in soybean [Glycine max (L.) Merrill].

BACKGROUND: The hundred seed weight (HSW) is one of the yield components of soybean [Glycine max (L.) Merrill] and is especially critical for various soybean food types. In this study, a representative sample consisting of 185 accessions was selected from Northeast China and analysed in three tested environments to determine the quantitative trait nucleotide (QTN) of HSW through a genome-wide association study (GWAS). RESULT: A total of 24,180 single nucleotide polymorphisms (SNPs) with minor allele frequencies greater than 0.2 and missing data less than 3% were utilized to estimate linkage disequilibrium (LD) levels in the tested association panel. Thirty-four association signals were identified as associated with HSW via GWAS. Among them, nineteen QTNs were novel, and another fifteen QTNs were overlapped or located near the genomic regions of known HSW QTL. A total of 237 genes, derived from 31 QTNs and located near peak SNPs from the three tested environments in 2015 and 2016, were considered candidate genes, were related to plant growth regulation, hormone metabolism, cell, RNA, protein metabolism, development, starch accumulation, secondary metabolism, signalling, and the TCA cycle, some of which have been found to participate in the regulation of HSW. A total of 106 SNPs from 16 candidate genes were significantly associated with HSW in soybean. CONCLUSIONS: The identified loci with beneficial alleles and candidate genes might be valuable for the molecular network and MAS of HSW.

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.