Exciton engineering of 2D Ruddlesden-Popper perovskites by synergistically tuning the intra and interlayer structures.

Designing two-dimensional halide perovskites for high-performance optoelectronic applications requires deep understanding of the structure-property relationship that governs their excitonic behaviors. However, a design framework that considers both intra and interlayer structures modified by the A-site and spacer cations, respectively, has not been developed. Here, we use pressure to synergistically tune the intra and interlayer structures and uncover the structural modulations that result in improved optoelectronic performance. Under applied pressure, (BA)2(GA)Pb2I7 exhibits a 72-fold boost of photoluminescence and 10-fold increase of photoconductivity. Based on the observed structural change, we introduce a structural descriptor χ that describes both the intra and interlayer characteristics and establish a general quantitative relationship between χ and photoluminescence quantum yield: smaller χ correlates with minimized trapped excitons and more efficient emission from free excitons. Building on this principle, we design a perovskite (CMA)2(FA)Pb2I7 that exhibits a small χ and an impressive photoluminescence quantum yield of 59.3%.

Two-dimensional lead halide perovskite lateral homojunctions enabled by phase pinning.

Two-dimensional organic-inorganic hybrid halide perovskites possess diverse structural polymorphs with versatile physical properties, which can be controlled by order-disorder transition of the spacer cation, making them attractive for constructing semiconductor homojunctions. Here, we demonstrate a space-cation-dopant-induced phase stabilization approach to creating a lateral homojunction composed of ordered and disordered phases within a two-dimensional perovskite. By doping a small quantity of pentylammonium into (butylammonium)2PbI4 or vice versa, we effectively suppress the ordering transition of the spacer cation and the associated out-of-plane octahedral tilting in the inorganic framework, resulting in phase pining of the disordered phase when decreasing temperature or increasing pressure. This enables epitaxial growth of a two-dimensional perovskite homojunction with tunable optical properties under temperature and pressure stimuli, as well as directional exciton diffusion across the interface. Our results demonstrate a previously unexplored strategy for constructing two-dimensional perovskite heterostructures by thermodynamic tuning and spacer cation doping.

The domestication-associated L1 gene encodes a eucomic acid synthase pleiotropically modulating pod pigmentation and shattering in soybean.

Pod coloration is a domestication-related trait in soybean, with modern cultivars typically displaying brown or tan pods, while their wild relative, Glycine soja, possesses black pods. However, the factors regulating this color variation remain unknown. In this study, we cloned and characterized L1, the classical locus responsible for black pods in soybean. By using map-based cloning and genetic analyses, we identified the causal gene of L1 and revealed that it encodes a hydroxymethylglutaryl-coenzyme A (CoA) lyase-like (HMGL-like) domain protein. Biochemical assays showed that L1 functions as a eucomic acid synthase and facilitates the synthesis of eucomic acid and piscidic acid, both of which contribute to coloration of pods and seed coats in soybean. Interestingly, we found that L1 plants are more prone to pod shattering under light exposure than l1 null mutants because dark pigmentation increases photothermal efficiency. Hence, pleiotropic effects of L1 on pod color and shattering, as well as seed pigmentation, likely contributed to the preference for l1 alleles during soybean domestication and improvement. Collectively, our study provides new insights into the mechanism of pod coloration and identifies a new target for future de novo domestication of legume crops.

Genome-wide signatures of the geographic expansion and breeding of soybean.

Soybean is a leguminous crop that provides oil and protein. Exploring the genomic signatures of soybean evolution is crucial for breeding varieties with improved adaptability to environmental extremes. We analyzed the genome sequences of 2,214 soybeans and proposed a soybean evolutionary route, i.e., the expansion of annual wild soybean (Glycine soja Sieb. & Zucc.) from southern China and its domestication in central China, followed by the expansion and local breeding selection of its landraces (G. max (L.) Merr.). We observed that the genetic introgression in soybean landraces was mostly derived from sympatric rather than allopatric wild populations during the geographic expansion. Soybean expansion and breeding were accompanied by the positive selection of flowering time genes, including GmSPA3c. Our study sheds light on the evolutionary history of soybean and provides valuable genetic resources for its future breeding.

Genome-wide association studies for soybean epicotyl length in two environments using 3VmrMLM.

Germination of soybean seed is the imminent vital process after sowing. The status of plumular axis and radicle determine whether soybean seed can emerge normally. Epicotyl, an organ between cotyledons and first functional leaves, is essential for soybean seed germination, seedling growth and early morphogenesis. Epicotyl length (EL) is a quantitative trait controlled by multiple genes/QTLs. Here, the present study analyzes the phenotypic diversity and genetic basis of EL using 951 soybean improved cultivars and landraces from Asia, America, Europe and Africa. 3VmrMLM was used to analyze the associations between EL in 2016 and 2020 and 1,639,846 SNPs for the identification of QTNs and QTN-by-environment interactions (QEIs)".A total of 180 QTNs and QEIs associated with EL were detected. Among them, 74 QTNs (ELS_Q) and 16 QEIs (ELS_QE) were identified to be associated with ELS (epicotyl length of single plant emergence), and 60 QTNs (ELT_Q) and 30 QEIs (ELT_QE) were identified to be associated with ELT (epicotyl length of three seedlings). Based on transcript abundance analysis, GO (Gene Ontology) enrichment and haplotype analysis, ten candidate genes were predicted within nine genic SNPs located in introns, upstream or downstream, which were supposed to be directly or indirectly involved in the process of seed germination and seedling development., Of 10 candidate genes, two of them (Glyma.04G122400 and Glyma.18G183600) could possibly affect epicotyl length elongation. These results indicate the genetic basis of EL and provides a valuable basis for specific functional studies of epicotyl traits.

Soybean genetic resources contributing to sustainable protein production.

KEY MESSAGE: Genetic resources contributes to the sustainable protein production in soybean. Soybean is an important crop for food, oil, and forage and is the main source of edible vegetable oil and vegetable protein. It plays an important role in maintaining balanced dietary nutrients for human health. The soybean protein content is a quantitative trait mainly controlled by gene additive effects and is usually negatively correlated with agronomic traits such as the oil content and yield. The selection of soybean varieties with high protein content and high yield to secure sustainable protein production is one of the difficulties in soybean breeding. The abundant genetic variation of soybean germplasm resources is the basis for overcoming the obstacles in breeding for soybean varieties with high yield and high protein content. Soybean has been cultivated for more than 5000 years and has spread from China to other parts of the world. The rich genetic resources play an important role in promoting the sustainable production of soybean protein worldwide. In this paper, the origin and spread of soybean and the current status of soybean production are reviewed; the genetic characteristics of soybean protein and the distribution of resources are expounded based on phenotypes; the discovery of soybean seed protein-related genes as well as transcriptomic, metabolomic, and proteomic studies in soybean are elaborated; the creation and utilization of high-protein germplasm resources are introduced; and the prospect of high-protein soybean breeding is described.

Understanding Electron-Phonon Interactions in 3D Lead Halide Perovskites from the Stereochemical Expression of 6s2 Lone Pairs.

The electron-phonon (e-ph) interaction in lead halide perovskites (LHPs) plays a role in a variety of physical phenomena. Unveiling how the local lattice distortion responds to charge carriers is a critical step toward understanding the e-ph interaction in LHPs. Herein, we advance a fundamental understanding of the e-ph interaction in LHPs from the perspective of stereochemical activity of 6s2 lone-pair electrons on the Pb2+ cation. We demonstrate a model system based on three LHPs with distinctive lone-pair activities for studying the structure-property relationships. By tuning the A-cation chemistry, we synthesized single-crystal CsPbBr3, (MA0.13EA0.87)PbBr3 (MA+ = methylammonium; EA+ = ethylammonium), and (MHy)PbBr3 (MHy+ = methylhydrazinium), which exhibit stereo-inactive, dynamic stereo-active, and static stereo-active lone pairs, respectively. This gives rise to distinctive local lattice distortions and low-frequency vibrational modes. We find that the e-ph interaction leads to a blue shift of the band gap as temperature increases in the structure with the dynamic stereo-active lone pair but to a red shift in the structure with the static stereo-active lone pair. Furthermore, analyses of the temperature-dependent low-energy photoluminescence tails reveal that the strength of the e-ph interaction increases with increasing lone-pair activity, leading to a transition from a large polaron to a small polaron, which has significant influence on the emission spectra and charge carrier dynamics. Our results highlight the role of the lone-pair activity in controlling the band gap, phonon, and polaronic effect in LHPs and provide guidelines for optimizing the optoelectronic properties, especially for tin-based and germanium-based halide perovskites, where stereo-active lone pairs are more prominent than their lead counterparts.

Dissection of the practical soybean breeding pipeline by developing ZDX1, a high-throughput functional array.

KEY MESSAGE: We developed the ZDX1 high-throughput functional soybean array for high accuracy evaluation and selection of both parents and progeny, which can greatly accelerate soybean breeding. Microarray technology facilitates rapid, accurate, and economical genotyping. Here, using resequencing data from 2214 representative soybean accessions, we developed the high-throughput functional array ZDX1, containing 158,959 SNPs, covering 90.92% of soybean genes and sites related to important traits. By application of the array, a total of 817 accessions were genotyped, including three subpopulations of candidate parental lines, parental lines and their progeny from practical breeding. The fixed SNPs were identified in progeny, indicating artificial selection during the breeding process. By identifying functional sites of target traits, novel soybean cyst nematode-resistant progeny and maturity-related novel sources were identified by allele combinations, demonstrating that functional sites provide an efficient method for the rapid screening of desirable traits or gene sources. Notably, we found that the breeding index (BI) was a good indicator for progeny selection. Superior progeny were derived from the combination of distantly related parents, with at least one parent having a higher BI. Furthermore, new combinations based on good performance were proposed for further breeding after excluding redundant and closely related parents. Genomic best linear unbiased prediction (GBLUP) analysis was the best analysis method and achieved the highest accuracy in predicting four traits when comparing SNPs in genic regions rather than whole genomic or intergenic SNPs. The prediction accuracy was improved by 32.1% by using progeny to expand the training population. Collectively, a versatile assay demonstrated that the functional ZDX1 array provided efficient information for the design and optimization of a breeding pipeline for accelerated soybean breeding.

Transcriptome analysis reveals differing response and tolerance mechanism of EPSPS and GAT genes among transgenic soybeans.

BACKGROUND: Glyphosate is a broad-spectrum, non-selective systemic herbicide. Introduction of glyphosate tolerance genes such as EPSPS or detoxification genes such as GAT can confer glyphosate tolerance on plants. Our previous study revealed that co-expression of EPSPS and GAT genes conferred higher glyphosate tolerance without "yellow flashing". However, the plant response to glyphosate at the transcriptional level was not investigated. METHODS AND RESULTS: To investigate the glyphosate tolerance mechanism, RNA-seq was conducted using four soybean genotypes, including two non-transgenic (NT) soybeans, ZH10 and MD12, and two GM soybeans, HJ698 and ZH10-6. Differentially expressed genes (DEGs) were identified in these soybeans before and after glyphosate treatment. Similar response to glyphosate in the two NT soybeans and the different effects of glyphosate on the two GM soybeans were identified. As treatment time was prolonged, the expression level of some DEGs involved in shikimate biosynthetic pathway and herbicide targeted cross-pathways was increased or declined continuously in NT soybeans, and altered slightly in HJ698. However, the expression level of some DEGs was altered in ZH10-6 at 12 hpt, while similar expression level of some DEGs involved in shikimate biosynthetic pathway and herbicide targeted cross-pathways was observed in ZH10-6 at 0 hpt and 72 hpt. These observations likely explain the higher glyphosate tolerance in ZH10-6 than in HJ698 and NT soybeans. CONCLUSIONS: These results suggested that GAT and EPSPS genes together play a crucial role in response to glyphosate, the GAT gene may work at the early stage of glyphosate exposure, whereas the EPSPS gene may be activated after the uptake of glyphosate by plants. These findings will provide valuable insight for the molecular basis underlying glyphosate tolerance or glyphosate detoxication.

GmMs1 encodes a kinesin-like protein essential for male fertility in soybean (Glycine max L.).

The application of heterosis is a promising approach for greatly increasing yield in soybean (Glycine max L.). Nuclear male sterility is essential for hybrid seed production and the utilization of heterosis. Here we report the cloning of the gene underlying the soybean male-sterile mutant ms-1, which has been widely used for recurrent selection in soybean breeding programs. We initially delimited the ms1 locus to a 16.15 kb region on chromosome 13, based on SLAF_BSA sequencing followed by genotyping of an F2 population segregating for the locus. Compared with the same region in fertile plants, the mutant region lacks a sequence of approximately 38.7 kb containing five protein-coding genes, including an ortholog of the kinesin-like protein gene NACK2, named GmMs1. The GmMs1 knockout plants generated via CRISPR/Cas-mediated gene editing displayed a complete male-sterile phenotype. Metabolic profiling showed that fertile anthers accumulated starch and sucrose normally, whereas sterile anthers had higher anthocyanin levels and lower flavonoid levels and lower antioxidant enzyme activities. These results provide insights into the molecular mechanisms governing male sterility and demonstrate that GmMs1 could be used to create male-sterile lines through targeted mutagenesis. These findings pave the way for designing seed production technology and an intelligent male-sterile line system to utilize heterosis in soybean.

Seed isoflavone profiling of 1168 soybean accessions from major growing ecoregions in China.

Soybean (Glycine max L. Merrill) isoflavones are secondary metabolites of great interest because of their beneficial impact on human health. We profiled the seed isoflavone composition of 1168 soybean accessions collected from diverse ecoregions of China in three locations over two years. We observed significant differences in isoflavone content among the accessions, accession types, years of growth and ecoregions of origin. Total isoflavone (TIF) concentration of the soybean accessions ranged from 745 µg g-1 to 5253.98 µg g-1, which represents a 7-fold difference. The highest mean TIF concentration (2689.27 µg g-1) was observed in the Huang Huai Hai Valley Region (HR) accessions, followed by accessions from the Southern Region (SR) and Northern Region (NR) with TIF concentration of 2518.91 µg g-1 and 1942.78 µg g-1, respectively. Thirty-five accessions were identified as elite soybean resources based on their higher TIF concentration (4024.74 µg g-1 to 5253.98 µg g-1). Pairwise correlation analysis showed significant positive correlations between individual isoflavones and TIF concentrations. Malonyldaidzin and malonylgenistin showed the highest correlations with TIF concentration (r = 0.90 and r = 0.92, respectively), whereas acetyldaidzin showed the lowest correlation. The main isoflavone components had significant negative correlations with latitude and longitude, indicating that the geographical origin of the accessions influenced their seed isoflavone composition. Based on principal component analysis, glycosides and malonylglycosides of isoflavones were the major discriminative components for the soybean accessions. The present study demonstrated the geographical distribution of soybean seed isoflavone concentrations across the main ecoregion of China. The identified soybean accessions with both high and low TIF concentrations, which are desirable materials for industrial uses and could also be used as parents to breed soybean lines with improved isoflavone quantity and composition in the seeds.

A Domestication-Associated Gene GmPRR3b Regulates the Circadian Clock and Flowering Time in Soybean.

Improved soybean cultivars have been adapted to grow at a wide range of latitudes, enabling expansion of cultivation worldwide. However, the genetic basis of this broad adaptation is still not clear. Here, we report the identification of GmPRR3b as a major flowering time regulatory gene that has been selected during domestication and genetic improvement for geographic expansion. Through a genome-wide association study of a diverse soybean landrace panel consisting of 279 accessions, we identified 16 candidate quantitative loci associated with flowering time and maturity time. The strongest signal resides in the known flowering gene E2, verifying the effectiveness of our approach. We detected strong signals associated with both flowering and maturity time in a genomic region containing GmPRR3b. Haplotype analysis revealed that GmPRR3bH6 is the major form of GmPRR3b that has been utilized during recent breeding of modern cultivars. mRNA profiling analysis showed that GmPRR3bH6 displays rhythmic and photoperiod-dependent expression and is preferentially induced under long-day conditions. Overexpression of GmPRR3bH6 increased main stem node number and yield, while knockout of GmPRR3bH6 using CRISPR/Cas9 technology delayed growth and the floral transition. GmPRR3bH6 appears to act as a transcriptional repressor of multiple predicted circadian clock genes, including GmCCA1a, which directly upregulates J/GmELF3a to modulate flowering time. The causal SNP (Chr12:5520945) likely endows GmPRR3bH6 a moderate but appropriate level of activity, leading to early flowering and vigorous growth traits preferentially selected during broad adaptation of landraces and improvement of cultivars.

Identification of loci controlling adaptation in Chinese soya bean landraces via a combination of conventional and bioclimatic GWAS.

Landraces often contain genetic diversity that has been lost in modern cultivars, including alleles that confer enhanced local adaptation. To comprehensively identify loci associated with adaptive traits in soya bean landraces, for example flowering time, a population of 1938 diverse landraces and 97 accessions of the wild progenitor of cultivated soya bean, Glycine soja was genotyped using tGBS® . Based on 99 085 high-quality SNPs, landraces were classified into three sub-populations which exhibit geographical genetic differentiation. Clustering was inferred from STRUCTURE, principal component analyses and neighbour-joining tree analyses. Using phenotypic data collected at two locations separated by 10 degrees of latitude, 17 trait-associated SNPs (TASs) for flowering time were identified, including a stable locus Chr12:5914898 and previously undetected candidate QTL/genes for flowering time in the vicinity of the previously cloned flowering genes, E1 and E2. Using passport data associated with the collection sites of the landraces, 27 SNPs associated with adaptation to three bioclimatic variables (temperature, daylength, and precipitation) were identified. A series of candidate flowering genes were detected within linkage disequilibrium (LD) blocks surrounding 12 bioclimatic TASs. Nine of these TASs exhibit significant differences in flowering time between alleles within one or more of the three individual sub-populations. Signals of selection during domestication and/or subsequent landrace diversification and adaptation were detected at 38 of the 44 flowering and bioclimatic TASs. Hence, this study lays the groundwork to begin breeding for novel environments predicted to arise following global climate change.

Identification of the dwarf gene GmDW1 in soybean (Glycine max L.) by combining mapping-by-sequencing and linkage analysis.

KEY MESSAGE: GmDW1 encodes an ent-kaurene synthase (KS) acting at the early step of the biosynthesis pathway for gibberellins (GAs) and regulates the development of plant height in soybean. Plant height is an important component of plant architecture, and significantly affects crop breeding practices and yield. Here, we report the characterization of an EMS-induced dwarf mutant (dw) of the soybean cultivar Zhongpin 661 (ZDD23893). The dw mutant displayed reduced plant height and shortened internodes, both of which were mainly attributed to the longitudinally decreased cell length. The bioactive GA1 (gibberellin A1) and GA4 (gibberellin A4) were not detectable in the stem of dw, and the dwarf phenotype could be rescued by treatment with exogenous GA3. Genetic analysis showed that the dwarf trait of dw was controlled by a recessive nuclear gene. By combining linkage analysis and mapping-by-sequencing, we mapped the GmDW1 gene to an approximately 460-kb region on chromosome (Chr.) 8, containing 36 annotated genes in the reference Willliams 82 genome. Of these genes, we identified two nonsynonymous single nucleotide polymorphisms (SNPs) that are present in the encoding regions of Gmdw1 and Glyma.08G165100 in dw, respectively. However, only the SNP mutation (T>A) at nucleotide 1224 in Gmdw1 cosegregated with the dwarf phenotype. GmDW1 encodes an ent-kaurene synthase, and was expressed in various tissues including root, stem, and leaf. Further phenotypic analysis of the allelic variations in soybean accessions strongly indicated that GmDW1 is responsible for the dwarf phenotype in dw. Our results provide important information for improving our understanding of the genetics of soybean plant height and crop breeding.

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.

Development and utilization of a new chemically-induced soybean library with a high mutation density .

Mutagenized populations have provided important materials for introducing variation and identifying gene function in plants. In this study, an ethyl methanesulfonate (EMS)-induced soybean (Glycine max) population, consisting of 21,600 independent M2 lines, was developed. Over 1,000 M4 (5) families, with diverse abnormal phenotypes for seed composition, seed shape, plant morphology and maturity that are stably expressed across different environments and generations were identified. Phenotypic analysis of the population led to the identification of a yellow pigmentation mutant, gyl, that displayed significantly decreased chlorophyll (Chl) content and abnormal chloroplast development. Sequence analysis showed that gyl is allelic to MinnGold, where a different single nucleotide polymorphism variation in the Mg-chelatase subunit gene (ChlI1a) results in golden yellow leaves. A cleaved amplified polymorphic sequence marker was developed and may be applied to marker-assisted selection for the golden yellow phenotype in soybean breeding. We show that the newly developed soybean EMS mutant population has potential for functional genomics research and genetic improvement in soybean.

GmSALT3, Which Confers Improved Soybean Salt Tolerance in the Field, Increases Leaf Cl- Exclusion Prior to Na+ Exclusion But Does Not Improve Early Vigor under Salinity.

Soil salinity reduces soybean growth and yield. The recently identified GmSALT3 (Glycine max salt Tolerance-associated gene on chromosome 3) has the potential to improve soybean yields in salinized conditions. Here we evaluate the impact of GmSALT3 on soybean performance under saline or non-saline conditions. Three sets of near isogenic lines (NILs), with genetic similarity of 95.6-99.3% between each pair of NIL-T and NIL-S, were generated from a cross between two varieties 85-140 (salt-sensitive, S) and Tiefeng 8 (salt-tolerant, T) by using marker-assisted selection. Each NIL-T; 782-T, 820-T and 860-T, contained a common ~1000 kb fragment on chromosome 3 where GmSALT3 was located. We show that GmSALT3 does not contribute to an improvement in seedling emergence rate or early vigor under salt stress. However, when 12-day-old seedlings were exposed to NaCl stress, the NIL-T lines accumulated significantly less leaf Na+ compared with their corresponding NIL-S, while no significant difference of K+ concentration was observed between NIL-T and NIL-S; the magnitude of Na+ accumulation within each NIL-T set was influenced by the different genetic backgrounds. In addition, NIL-T lines accumulated less Cl- in the leaf and more in the root prior to any difference in Na+; in the field they accumulated less pod wall Cl- than the corresponding NIL-S lines. Under non-saline field conditions, no significant differences were observed for yield related traits within each pair of NIL-T and NIL-S lines, indicating there was no yield penalty for having the GmSALT3 gene. In contrast, under saline field conditions the NIL-T lines had significantly greater plant seed weight and 100-seed weight than the corresponding NIL-S lines, meaning GmSALT3 conferred a yield advantage to soybean plants in salinized fields. Our results indicated that GmSALT3 mediated regulation of both Na+ and Cl- accumulation in soybean, and contributes to improved soybean yield through maintaining a higher seed weight under saline stress.

Identification of Genomic Insertion and Flanking Sequence of G2-EPSPS and GAT Transgenes in Soybean Using Whole Genome Sequencing Method.

Molecular characterization of sequence flanking exogenous fragment insertion is essential for safety assessment and labeling of genetically modified organism (GMO). In this study, the T-DNA insertion sites and flanking sequences were identified in two newly developed transgenic glyphosate-tolerant soybeans GE-J16 and ZH10-6 based on whole genome sequencing (WGS) method. More than 22.4 Gb sequence data (∼21 × coverage) for each line was generated on Illumina HiSeq 2500 platform. The junction reads mapped to boundaries of T-DNA and flanking sequences in these two events were identified by comparing all sequencing reads with soybean reference genome and sequence of transgenic vector. The putative insertion loci and flanking sequences were further confirmed by PCR amplification, Sanger sequencing, and co-segregation analysis. All these analyses supported that exogenous T-DNA fragments were integrated in positions of Chr19: 50543767-50543792 and Chr17: 7980527-7980541 in these two transgenic lines. Identification of genomic insertion sites of G2-EPSPS and GAT transgenes will facilitate the utilization of their glyphosate-tolerant traits in soybean breeding program. These results also demonstrated that WGS was a cost-effective and rapid method for identifying sites of T-DNA insertions and flanking sequences in soybean.

Identification and Validation of Loci Governing Seed Coat Color by Combining Association Mapping and Bulk Segregation Analysis in Soybean.

Soybean seed coat exists in a range of colors from yellow, green, brown, black, to bicolor. Classical genetic analysis suggested that soybean seed color was a moderately complex trait controlled by multi-loci. However, only a couple of loci could be detected using a single biparental segregating population. In this study, a combination of association mapping and bulk segregation analysis was employed to identify genes/loci governing this trait in soybean. A total of 14 loci, including nine novel and five previously reported ones, were identified using 176,065 coding SNPs selected from entire SNP dataset among 56 soybean accessions. Four of these loci were confirmed and further mapped using a biparental population developed from the cross between ZP95-5383 (yellow seed color) and NY279 (brown seed color), in which different seed coat colors were further dissected into simple trait pairs (green/yellow, green/black, green/brown, yellow/black, yellow/brown, and black/brown) by continuously developing residual heterozygous lines. By genotyping entire F2 population using flanking markers located in fine-mapping regions, the genetic basis of seed coat color was fully dissected and these four loci could explain all variations of seed colors in this population. These findings will be useful for map-based cloning of genes as well as marker-assisted breeding in soybean. This work also provides an alternative strategy for systematically isolating genes controlling relative complex trait by association analysis followed by biparental mapping.