Genome-wide association study of bacterial blight resistance in soybean.

Bacterial blight caused by Pseudomonas syringae pv. glycines (Psg) is a widespread foliar disease. Although four Resistance to Pseudomonas syringae pv. glycinea (Rpg) 1 ~ 4 (Rpg1~4) genes that have been observed to segregate in a Mendelian pattern have been reported to confer resistance to Psg in soybean, the genetic basis of quantitative resistance to bacterial blight in soybean remains unclear. In the present study, the Psg resistance of two soybean association panels consisting of 573 and 213 lines, respectively, were phenotyped in multiple environments in 2014 - 2016. Genome-wide association study (GWAS) were performed using 2 models FarmCPU and BLINK to identify Psg resistance loci. A total of 40 soybean varieties with high level of Psg resistance were identified, and 14 quantitative trait loci (QTLs) were detected on 12 soybean chromosomes. These QTLs were identified for the first time. The majority of the QTLs were only detected in one or the other association panels, while qRPG-18-1 was detected in both association panels for at least one growing season. A total of 46 candidate Psg resistance genes were identified from the qRpg_13_1, qRPG-15-1, and qRPG-18-1 loci based on gene function annotation. In addition, we found the genomic region covering rpg1-b and rpg1-r harbored the synteny with a genomic region on chromosome 15, and identified 16 nucleotide binding site - leucine-rich repeat (NBS-LRR) genes as the candidate Psg resistance genes from the synteny blocks. This study provides new information for dissecting the genetic control of Psg resistance in soybean.

Identification of Novel Genomic Regions for Bacterial Leaf Pustule (BLP) Resistance in Soybean (Glycine max L.) via Integrating Linkage Mapping and Association Analysis.

Bacterial leaf pustule (BLP), caused by Xanthornonas axonopodis pv. glycines (Xag), is a worldwide disease of soybean, particularly in warm and humid regions. To date, little is known about the underlying molecular mechanisms of BLP resistance. The only single recessive resistance gene rxp has not been functionally identified yet, even though the genotypes carrying the gene have been widely used for BLP resistance breeding. Using a linkage mapping in a recombinant inbred line (RIL) population against the Xag strain Chinese C5, we identified that quantitative trait locus (QTL) qrxp-17-2 accounted for 74.33% of the total phenotypic variations. We also identified two minor QTLs, qrxp-05-1 and qrxp-17-1, that accounted for 7.26% and 22.26% of the total phenotypic variations, respectively, for the first time. Using a genome-wide association study (GWAS) in 476 cultivars of a soybean breeding germplasm population, we identified a total of 38 quantitative trait nucleotides (QTNs) on chromosomes (Chr) 5, 7, 8, 9,15, 17, 19, and 20 under artificial infection with C5, and 34 QTNs on Chr 4, 5, 6, 9, 13, 16, 17, 18, and 20 under natural morbidity condition. Taken together, three QTLs and 11 stable QTNs were detected in both linkage mapping and GWAS analysis, and located in three genomic regions with the major genomic region containing qrxp_17_2. Real-time RT-PCR analysis of the relative expression levels of five potential candidate genes in the resistant soybean cultivar W82 following Xag treatment showed that of Glyma.17G086300, which is located in qrxp-17-2, significantly increased in W82 at 24 and 72 h post-inoculation (hpi) when compared to that in the susceptible cultivar Jack. These results indicate that Glyma.17G086300 is a potential candidate gene for rxp and the QTLs and QTNs identified in this study will be useful for marker development for the breeding of Xag-resistant soybean cultivars.

Hydroperoxide lyase modulates defense response and confers lesion-mimic leaf phenotype in soybean (Glycine max (L.) Merr.).

Allene oxide synthase (AOS) and hydroperoxide lyase (HPL) are two important members of P450 enzymes metabolizing hydroperoxy fatty acid to produce jasmonates and aldehydes respectively, which function in response to diverse environmental and developmental stimuli. However, their exact roles in soybean have not been clarified. In present study, we identified a lesion-mimic mutant in soybean named NT302, which exhibits etiolated phenotype together with chlorotic and spontaneous lesions on leaves at R3 podding stage. The underlying gene was identified as GmHPL encoding hydroperoxide lyase by map-based cloning strategy. Sequence analysis demonstrated that a single nucleotide mutation created a premature termination codon (Gln20-Ter), which resulted in a truncated GmHPL protein in NT302. GmHPL RNA was significantly reduced in NT302 mutant, while genes in AOS branch of the 13-LOX pathway were up-regulated in NT302. The mutant exhibited higher susceptibility to bacterial leaf pustule (BLP) disease, but increased resistance against common cutworm (CCW) pest. GmHPL was significantly induced in response to MeJA, wounding, and CCW in wild type soybean. Virus induced gene silencing (VIGS) of GhHPL genes gave rise to similar lesion-mimic leaf phenotypes in upland cotton, coupled with upregulation of the expression of JA biosynthesis and JA-induced genes. Our study provides evidence that competition exist between HPL and AOS branches in 13-LOX pathway of the oxylipin metabolism in soybean, thereby plays essential roles in modulation of plant development and defense.

Double mutation of two homologous genes YL1 and YL2 results in a leaf yellowing phenotype in soybean [Glycine max (L.) Merr].

KEY MESSAGE: Two homologous, chloroplast located CAAX proteases were identified to be functional redundancy in determining soybean leaf color, and they probably play essential roles in regulating light harvesting and absorption during photosynthesis process. Leaf color mutants are ideal materials for studying photosynthesis and chlorophyll metabolism. The soybean [Glycine max (L.) Merr.] yellowing leaf (yl) variation is a recombinant mutant characterized by yellow foliage, which derived from the specific cross between green seed-coated and yellow seed-coated soybean varieties. Molecular cloning and subsequent gene silencing revealed that the yellow leaf trait of yl was controlled by two recessive nuclear genes, glyma11g04660 and glyma01g40650, named as YL1 and YL2 respectively, and the latter was confirmed to be same as the earlier reported green seed-coat gene G. Both YL1 and YL2 belonged to chloroplast-located proteases possessing Abi domain, and these genes were expressed in various tissues, especially in young leaves. In yl, the expression of YL1 and YL2 were suppressed in most tissues, and the young leaf of yl presented an increased maximal photochemical efficiency (Fv/Fm) as well as enhanced net photosynthesis activity (Pn), indicating that YL1 and YL2 are involved in light absorption regulation during photosynthesis process. Collectively, the identification and description of YL1 and YL2 in our study provides insights for the regulatory mechanism of photosynthesis process, and these findings will further assist to clarify the close relationship between photosynthesis and chlorophyll metabolism.

Whole-genome mapping identified novel "QTL hotspots regions" for seed storability in soybean (Glycine max L.).

BACKGROUND: Seed aging in soybean is a serious challenge for agronomic production and germplasm preservation. However, its genetic basis remains largely unclear in soybean. Unraveling the genetic mechanism involved in seed aging, and enhancing seed storability is an imperative goal for soybean breeding. The aim of this study is to identify quantitative trait loci (QTLs) using high-density genetic linkage maps of soybean for seed storability. In this regard, two recombinant inbred line (RIL) populations derived from Zhengyanghuangdou × Meng 8206 (ZM6) and Linhefenqingdou × Meng 8206 (LM6) crosses were evaluated for three seed-germination related traits viz., germination rate (GR), normal seedling length (SL) and normal seedling fresh weight (FW) under natural and artificial aging conditions to map QTLs for seed storability. RESULTS: A total of 34 QTLs, including 13 QTLs for GR, 11 QTLs for SL and 10 QTLs for FW, were identified on 11 chromosomes with the phenotypic variation ranged from 7.30 to 23.16% under both aging conditions. All these QTLs were novel, and 21 of these QTLs were clustered in five QTL-rich regions on four different chromosomes viz., Chr3, Chr5, Chr17 &Chr18, among them the highest concentration of seven and six QTLs were found in "QTL hotspot A" (Chr17) and "QTL hotspot B" (Chr5), respectively. Furthermore, QTLs within all the five QTL clusters are linked to at least two studied traits, which is also supported by highly significant correlation between the three germination-related traits. QTLs for seed-germination related traits in "QTL hotspot B" were found in both RIL populations and aging conditions, and also QTLs underlying "QTL hotspot A" are identified in both RIL populations under artificial aging condition. These are the stable genomic regions governing the inheritance of seed storability in soybean, and will be the main focus for soybean breeders. CONCLUSION: This study uncovers the genetic basis of seed storability in soybean. The newly identified QTLs provides valuable information, and will be main targets for fine mapping, candidate gene identification and marker-assisted breeding. Hence, the present study is the first report for the comprehensive and detailed investigation of genetic architecture of seed storability in soybean.

Genome-Wide Detection of Major and Epistatic Effect QTLs for Seed Protein and Oil Content in Soybean Under Multiple Environments Using High-Density Bin Map.

Seed protein and oil content are the two important traits determining the quality and value of soybean. Development of improved cultivars requires detailed understanding of the genetic basis underlying the trait of interest. However, it is prerequisite to have a high-density linkage map for precisely mapping genomic regions, and therefore the present study used high-density genetic map containing 2267 recombination bin markers distributed on 20 chromosomes and spanned 2453.79 cM with an average distance of 1.08 cM between markers using restriction-site-associated DNA sequencing (RAD-seq) approach. A recombinant inbred line (RIL) population of 104 lines derived from a cross between Linhefenqingdou and Meng 8206 cultivars was evaluated in six different environments to identify main- and epistatic-effect quantitative trait loci (QTLs)as well as their interaction with environments. A total of 44 main-effect QTLs for protein and oil content were found to be distributed on 17 chromosomes, and 15 novel QTL were identified for the first time. Out of these QTLs, four were major and stable QTLs, viz., qPro-7-1, qOil-8-3, qOil-10-2 and qOil-10-4, detected in at least two environments plus combined environment with R² values >10%. Within the physical intervals of these four QTLs, 111 candidate genes were screened for their direct or indirect involvement in seed protein and oil biosynthesis/metabolism processes based on gene ontology and annotation information. Based on RNA sequencing (RNA-seq) data analysis, 15 of the 111 genes were highly expressed during seed development stage and root nodules that might be considered as the potential candidate genes. Seven QTLs associated with protein and oil content exhibited significant additive and additive × environment interaction effects, and environment-independent QTLs revealed higher additive effects. Moreover, three digenic epistatic QTLs pairs were identified, and no main-effect QTLs showed epistasis. In conclusion, the use of a high-density map identified closely linked flanking markers, provided better understanding of genetic architecture and candidate gene information, and revealed the scope available for improvement of soybean quality through marker assisted selection (MAS).

Characterization and Rapid Gene-Mapping of Leaf Lesion Mimic Phenotype of spl-1 Mutant in Soybean (Glycine max (L.) Merr.).

In plants, lesion mimic mutants (LMMs) reveal spontaneous disease-like lesions in the absence of pathogen that constitutes powerful genetic material to unravel genes underlying programmed cell death (PCD), particularly the hypersensitive response (HR). However, only a few LMMs are reported in soybean, and no related gene has been cloned until now. In the present study, we isolated a new LMM named spotted leaf-1 (spl-1) from NN1138-2 cultivar through ethyl methanesulfonate (EMS) treatment. The present study revealed that lesion formation might result from PCD and excessive reactive oxygen species (ROS) accumulation. The chlorophyll content was significantly reduced but antioxidant activities, viz., superoxide dismutase (SOD), peroxidase (POD) and catalase (CAT), as well as the malondialdehyde (MDA) contents, were detected higher in spl-1 than in the wild-type. According to segregation analysis of mutant phenotype in two genetic populations, viz., W82×spl-1 and PI378692×spl-1, the spotted leaf phenotype of spl-1 is controlled by a single recessive gene named lm1. The lm1 locus governing mutant phenotype of spl-1 was first identified in 3.15 Mb genomic region on chromosome 04 through MutMap analysis, which was further verified and fine mapped by simple sequence repeat (SSR) marker-based genetic mapping. Genetic linkage analysis narrowed the genomic region (lm1 locus) for mutant phenotype to a physical distance of ~76.23 kb. By searching against the Phytozome database, eight annotated candidate genes were found within the lm1 region. qRT-PCR expression analysis revealed that, among these eight genes, only Glyma.04g242300 showed highly significant expression levels in wild-type relative to the spl-1 mutant. However, sequencing data of the CDS region showed no nucleotide difference between spl-1 and its wild type within the coding regions of these genes but might be in the non-coding regions such as 5' or 3' UTR. Hence, the data of the present study are in favor of Glyma.04g242300 being the possible candidate genes regulating the mutant phenotype of spl-1. However, further validation is needed to prove this function of the gene as well as its role in PCD, which in turn would be helpful to understand the mechanism and pathways involved in HR disease resistance of soybean.

Genetic Dissection of Extreme Seed-Flooding Tolerance in a Wild Soybean PI342618B by Linkage Mapping and Candidate Gene Analysis.

Seed-flooding stress is one of major abiotic constraints that adversely affects soybean production worldwide. Identifying tolerant germplasms and revealing the genetic basis of seed-flooding tolerance are imperative goals for soybean breeding. In the present study, high-density linkage maps of two inter-specific recombinant inbred line (RIL) populations, named NJIRNP and NJIR4P, were utilized to identify major quantitative trait loci (QTLs) for seed-flooding tolerance using three parameters viz., germination rate (GR), normal seedling rate (NSR), and electrical conductivity (EC). A total of 25 and 18 QTLs were detected by composite interval mapping (CIM) and mixed-model-based composite interval mapping (MCIM), respectively, and 12 common QTLs were identified through both methods. All favorable alleles for the tolerance are notably from the wild soybean parent. Moreover, four digenic epistatic QTL pairs were identified, and three of them showed no main effects. In addition, the pigmented soybean genotypes exhibited high seed-flooding tolerance compared with yellow seed coat genotypes in both populations. Moreover, out of five identified QTLs, one major region containing multiple QTLs associated with all three traits was identified on Chromosome 8, and most of the QTLs within this hotspot were major loci (R2 > 10) and detectable in both populations and multiple environments. Based on the gene expression and functional annotation information, 10 candidate genes from QTL "hotspot 8-2" were screened for further analysis. Furthermore, the results of qRT-PCR and sequence analysis revealed that only one gene, GmDREB2 (Glyma.08G137600), was significantly induced under flooding stress and displayed a TTC tribasic insertion mutation of the nucleotide sequence in the tolerant wild parent (PI342618B). GmDREB2 encodes an ERF transcription factor, and the subcellular localization analysis using green fluorescent protein (GFP) revealed that GmDREB2 protein was localized in the nucleus and plasma membrane. Furthermore, overexpression of GmDREB2 significantly promoted the growth of soybean hairy roots, which might indicate its critical role in seed-flooding stress. Thus, GmDREB2 was considered as the most possible candidate gene for seed-flooding tolerance.

Adsorption of Cd (II) by a novel living and non-living Cupriavidus necator GX_5: optimization, equilibrium and kinetic studies.

Biosorbents have been extensively studied for heavy metal adsorption due to their advantages of low cost and high efficiency. In the study, the living and non-living biomass of Cupriavidus necator GX_5 previously isolated were evaluated for their adsorption capacity and/or removal efficiency for Cd (II) through batch experiments, SEM and FT-IR investigations. The maximum removal efficiency rates for the live and dead biomass were 60.51% and 78.53%, respectively, at an optimum pH of 6, a dosage of 1 g/L and an initial Cd (II) concentration of 5 mg/L. The pseudo-second-order kinetic model was more suitable for fitting the experimental data, indicating that the rate-limiting step might be chemisorption. The Freundlich isotherm model fit better than the Langmuir isotherm model, implying that the adsorption process of both biosorbents was heterogeneous. FT-IR observation reflected that various functional groups were involved in Cd (II) adsorption: -OH, -NH, C=O, C-O and C-C groups for the living biomass and -OH, -NH, C-H, C = O, C-N and N-H groups for the dead biomass. Our results imply that non-living biosorbents have a higher capacity and stronger strength for absorbing Cd (II) than living biomass. Therefore, we suggest that dead GX_5 is a promising adsorbent and can be used in Cd (II)-contaminated environments.

GmNAC8 acts as a positive regulator in soybean drought stress.

NAC proteins represent one of the largest transcription factor (TF) families involved in the regulation of plant development and the response to abiotic stress. In the present study, we elucidated the detailed role of GmNAC8 in the regulation of drought stress tolerance in soybean. The GmNAC8 protein was localized in the nucleus, and expression of the GmNAC8 gene was significantly induced in response to drought, abscisic acid (ABA), ethylene (ETH) and salicylic acid (SA) treatments. Thus, we generated GmNAC8 overexpression (OE1 and OE2) and GmNAC8 knockout (KO1 and KO2) lines to determine the role of GmNAC8 in drought stress tolerance. Our results revealed that, compared with the wild type (WT) plant, GmNAC8 overexpression and GmNAC8 knockout lines exhibited significantly higher and lower drought tolerance, respectively. Furthermore, the SOD activity and proline content were significantly higher in the GmNAC8 overexpression lines and significantly lower in the GmNAC8 knockout lines than in the WT plants under drought stress. In addition, GmNAC8 protein was found to physically interact with the drought-induced protein GmDi19-3 in the nucleus. Moreover, the GmDi19-3 expression pattern showed the same trend as the GmNAC8 gene did under drought and hormone (ABA, ETH and SA) treatments, and GmDi19-3 overexpression lines (GmDi19-3-OE9, GmDi19-3-OE10 and GmDi19-3-OE31) showed enhanced drought tolerance compared to that of the WT plants. Hence, the above results indicated that GmNAC8 acts as a positive regulator of drought tolerance in soybean and inferred that GmNAC8 probably functions by interacting with another positive regulatory protein, GmDi19-3.

QTL architecture of vine growth habit and gibberellin oxidase gene diversity in wild soybean (Glycine soja).

Vine growth habit (VGH) is a beneficial phenotype in many wild plants, and is considered an important domesticated-related trait in soybean. However, its genetic basis remains largely unclear. Hence, in the present study we used an integrated strategy combining linkage mapping and population genome diversity analyses to reveal the genetics of VGH in soybean. In this regard, two recombinant inbred line (RIL) populations derived by crossing a common wild soybean genotype (PI342618B) with two cultivated lines viz., NN 86-4 and NN 493-1 were used to map quantitative trait loci (QTL) for VGH. Here, we identified seven and five QTLs at flowering stage (R1) and maturity stage (R8), respectively, and among them qVGH-18-1, qVGH-18-2, qVGH-19-3, qVGH-19-4 were identified as major loci (R2 > 10% and detection time ≥2). However, qVGH-18-2 was considered as a main QTL for VGH being consistently identified in both RIL populations as well as all growth stages and cropping years. Out of all the annotated genes within qVGH-18-2, Glyma18g06870 was identified as the candidate gene and named as VGH1, which was a gibberellin oxidase (GAox) belongs to 2-oxoglutarate-dependent dioxygenase (2- ODD). Interestingly, there was one member of 2-ODD/GAox in qVGH-18-1 and qVGH-19-4 named as VGH2 and VGH3, respectively. Moreover, from sequencing data analysis VGH1 and three other GAox genes were found significantly divergent between vine and erect soybean with FST value larger than 0.25. Hence, GAox was assumed to play a major role in governing inheritance of VGH in soybean. Therefore, elucidating the genetic mechanism of GAox is very useful for exploring VGH and other stem traits, as well as genetic improvement of plant type in soybean.

Utilization of Interspecific High-Density Genetic Map of RIL Population for the QTL Detection and Candidate Gene Mining for 100-Seed Weight in Soybean.

Seed-weight is one of the most important traits determining soybean yield. Hence, it is prerequisite to have detailed understanding of the genetic basis regulating seed-weight for the development of improved cultivars. In this regard, the present study used high-density interspecific linkage map of NJIR4P recombinant inbred population evaluated in four different environments to detect stable Quantitative trait loci (QTLs) as well as mine candidate genes for 100-seed weight. In total, 19 QTLs distributed on 12 chromosomes were identified in all individual environments plus combined environment, out of which seven were novel and eight are stable identified in more than one environment. However, all the novel QTLs were minor (R 2 < 10%). The remaining 12 QTLs detected in this study were co-localized with the earlier reported QTLs with narrow genomic regions, and out of these only 2 QTLs were major (R 2 > 10%) viz., qSW-17-1 and qSW-17-4. Beneficial alleles of all identified QTLs were derived from cultivated soybean parent (Nannong493-1). Based on Protein ANalysis THrough Evolutionary Relationships, gene annotation information, and literature search, 29 genes within 5 stable QTLs were predicted to be possible candidate genes that might regulate seed-weight/size in soybean. However, it needs further validation to confirm their role in seed development. In conclusion, the present study provides better understanding of trait genetics and candidate gene information through the use high-density inter-specific bin map, and also revealed considerable scope for genetic improvement of 100-seed weight in soybean using marker-assisted breeding.

Deciphering the Genetic Architecture of Plant Height in Soybean Using Two RIL Populations Sharing a Common M8206 Parent.

Plant height (PH) is an important agronomic trait that is closely related to soybean yield and quality. However, it is a complex quantitative trait governed by multiple genes and is influenced by environment. Unraveling the genetic mechanism involved in PH, and developing soybean cultivars with desirable PH is an imperative goal for soybean breeding. In this regard, the present study used high-density linkage maps of two related recombinant inbred line (RIL) populations viz., MT and ZM evaluated in three different environments to detect additive and epistatic effect quantitative trait loci (QTLs) as well as their interaction with environments for PH in Chinese summer planting soybean. A total of eight and 12 QTLs were detected by combining the composite interval mapping (CIM) and mixed-model based composite interval mapping (MCIM) methods in MT and ZM populations, respectively. Among these QTLs, nine QTLs viz., QPH-2, qPH-6-2MT, QPH-6, qPH-9-1ZM, qPH-10-1ZM, qPH-13-1ZM, qPH-16-1MT, QPH-17 and QPH-19 were consistently identified in multiple environments or populations, hence were regarded as stable QTLs. Furthermore, Out of these QTLs, three QTLs viz., qPH-4-2ZM, qPH-15-1MT and QPH-17 were novel. In particular, QPH-17 could detect in both populations, which was also considered as a stable and major QTL in Chinese summer planting soybean. Moreover, eleven QTLs revealed significant additive effects in both populations, and out of them only six showed additive by environment interaction effects, and the environment-independent QTLs showed higher additive effects. Finally, six digenic epistatic QTLs pairs were identified and only four additive effect QTLs viz., qPH-6-2MT, qPH-19-1MT/QPH-19, qPH-5-1ZM and qPH-17-1ZM showed epistatic effects. These results indicate that environment and epistatic interaction effects have significant influence in determining genetic basis of PH in soybean. These results would not only increase our understanding of the genetic control of plant height in summer planting soybean but also provide support for implementing marker assisted selection (MAS) in developing cultivars with ideal plant height as well as gene cloning to elucidate the mechanisms of plant height.

Identification of QTN and Candidate Gene for Seed-flooding Tolerance in Soybean [Glycine max (L.) Merr.] using Genome-Wide Association Study (GWAS).

Seed-flooding stress is one of the major abiotic constraints severely affecting soybean yield and quality. Understanding the molecular mechanism and genetic basis underlying seed-flooding tolerance will be of greatly importance in soybean breeding. However, very limited information is available about the genetic basis of seed-flooding tolerance in soybean. The present study performed Genome-Wide Association Study (GWAS) to identify the quantitative trait nucleotides (QTNs) associated with three seed-flooding tolerance related traits, viz., germination rate (GR), normal seedling rate (NSR) and electric conductivity (EC), using a panel of 347 soybean lines and the genotypic data of 60,109 SNPs with MAF > 0.05. A total of 25 and 21 QTNs associated with all three traits were identified via mixed linear model (MLM) and multi-locus random-SNP-effect mixed linear model (mrMLM) in three different environments (JP14, HY15, and Combined). Among these QTNs, three major QTNs, viz., QTN13, qNSR-10 and qEC-7-2, were identified through both methods MLM and mrMLM. Interestingly, QTN13 located on Chr.13 has been consistently identified to be associated with all three studied traits in both methods and multiple environments. Within the 1.0 Mb physical interval surrounding the QTN13, nine candidate genes were screened for their involvement in seed-flooding tolerance based on gene annotation information and available literature. Based on the qRT-PCR and sequence analysis, only one gene designated as GmSFT (Glyma.13g248000) displayed significantly higher expression level in all tolerant genotypes compared to sensitive ones under flooding treatment, as well as revealed nonsynonymous mutation in tolerant genotypes, leading to amino acid change in the protein. Additionally, subcellular localization showed that GmSFT was localized in the nucleus and cell membrane. Hence, GmSFT was considered as the most likely candidate gene for seed-flooding tolerance in soybean. In conclusion, the findings of the present study not only increase our knowledge of the genetic control of seed-flooding tolerance in soybean, but will also be of great utility in marker-assisted selection and gene cloning to elucidate the mechanisms of seed-flooding tolerance.

Genome-Wide Association Studies for Dynamic Plant Height and Number of Nodes on the Main Stem in Summer Sowing Soybeans.

Plant height (PH) and the number of nodes on the main stem (NN) serve as major plant architecture traits affecting soybean seed yield. Although many quantitative trait loci for the two traits have been reported, their genetic controls at different developmental stages in soybeans remain unclear. Here, 368 soybean breeding lines were genotyped using 62,423 single nucleotide polymorphism (SNP) markers and phenotyped for the two traits at three different developmental stages over two locations in order to identify their quantitative trait nucleotides (QTNs) using compressed mixed linear model (CMLM) and multi-locus random-SNP-effect mixed linear model (mrMLM) approaches. As a result, 11 and 13 QTNs were found by CMLM to be associated with PH and NN, respectively. Among these QTNs, 8, 3, and 4 for PH and 6, 6, and 8 for NN were found at the three stages, and 3 and 6 were repeatedly detected for PH and NN. In addition, 34 and 30 QTNs were found by mrMLM to be associated with PH and NN, respectively. Among these QTNs, 11, 13, and 16 for PH and 11, 15, and 8 for NN were found at the three stages. A majority of these QTNs overlapped with the previously reported loci. Moreover, one QTN within the known E2 locus for flowering time was detected for the two traits at all three stages, and another that overlapped with the Dt1 locus for stem growth habit was also identified for the two traits at the mature stage. This may explain the highly significant correlation between the two traits. Our findings provide evidence for mixed major plus polygenes inheritance for dynamic traits and an extended understanding of their genetic architecture for molecular dissection and breeding utilization in soybeans.

Identification of Major Quantitative Trait Loci for Seed Oil Content in Soybeans by Combining Linkage and Genome-Wide Association Mapping.

Soybean oil is the most widely produced vegetable oil in the world and its content in soybean seed is an important quality trait in breeding programs. More than 100 quantitative trait loci (QTLs) for soybean oil content have been identified. However, most of them are genotype specific and/or environment sensitive. Here, we used both a linkage and association mapping methodology to dissect the genetic basis of seed oil content of Chinese soybean cultivars in various environments in the Jiang-Huai River Valley. One recombinant inbred line (RIL) population (NJMN-RIL), with 104 lines developed from a cross between M8108 and NN1138-2, was planted in five environments to investigate phenotypic data, and a new genetic map with 2,062 specific-locus amplified fragment markers was constructed to map oil content QTLs. A derived F2 population between MN-5 (a line of NJMN-RIL) and NN1138-2 was also developed to confirm one major QTL. A soybean breeding germplasm population (279 lines) was established to perform a genome-wide association study (GWAS) using 59,845 high-quality single nucleotide polymorphism markers. In the NJMN-RIL population, 8 QTLs were found that explained a range of phenotypic variance from 6.3 to 26.3% in certain planting environments. Among them, qOil-5-1, qOil-10-1, and qOil-14-1 were detected in different environments, and qOil-5-1 was further confirmed using the secondary F2 population. Three loci located on chromosomes 5 and 20 were detected in a 2-year long GWAS, and one locus that overlapped with qOil-5-1 was found repeatedly and treated as the same locus. qOil-5-1 was further localized to a linkage disequilibrium block region of approximately 440 kb. These results will not only increase our understanding of the genetic control of seed oil content in soybean, but will also be helpful in marker-assisted selection for breeding high seed oil content soybean and gene cloning to elucidate the mechanisms of seed oil content.

Characterizing Two Inter-specific Bin Maps for the Exploration of the QTLs/Genes that Confer Three Soybean Evolutionary Traits.

Annual wild soybean (Glycine soja Sieb. and Zucc.), the wild progenitor of the cultivated soybean [Glycine max (L.) Merr.], is valuable for improving the later. The construction of a linkage map is crucial for studying the genetic differentiation between these species, but marker density is the main factor limiting the accuracy of such a map. Recent advances in next-generation sequencing technologies allow for the generation of high-density linkage maps. Here, two sets of inter-specific recombinant inbred line populations, named NJIRNP and NJIR4P, composed of 284 and 161 lines, respectively, were generated from the same wild male parent, PI 342618B, and genotyped by restriction-site-associated DNA sequencing. Two linkage maps containing 5,728 and 4,354 bins were constructed based on 89,680 and 80,995 single nucleotide polymorphisms, spanning a total genetic distance of 2204.6 and 2136.7 cM, with an average distance of 0.4 and 0.5 cM between neighboring bins in NJRINP and NJRI4P, respectively. With the two maps, seven well-studied loci, B1 for seed bloom; G and I for seed coat color; E2, E3, qDTF16.1 and two linked FLOWERING LOCUS T for days to flowering, were detected. In addition, two SB and two DTF loci were newly identified in wild soybean. Using two high-density maps, the mapping resolution was enhanced, e.g., G was narrowed to a region of 0.4 Mb on chromosome 1, encompassing 54 gene models, among which only Glyma01g40590 was predicted to be involved in anthocyanin accumulation, and its interaction with I was verified in both populations. In addition, five genes, Glyma16g03030, orthologous to Arabidopsis Phytochrome A (PHYA); Glyma13g28810, Glyma13g29920, and Glyma13g30710 predicted to encode the APETALA 2 (AP2) domain; and Glyma02g00300, involved in response to red or far red light, might be candidate DTF genes. Our results demonstrate that RAD-seq is a cost-effective approach for constructing high-density and high-quality bin maps that can be used to map QTLs/genes into such small enough regions that their candidate genes can be predicted.