Identification and characterization of novel QTL conferring internal detoxification of aluminium in soybean.

Aluminium (Al) toxicity inhibits soybean root growth, leading to insufficient water and nutrient uptake. Two soybean lines ('Magellan' and PI 567731) were identified differing in Al tolerance, as determined by primary root length ratio, total root length ratio, and root tip number ratio under Al stress. Serious root necrosis was observed in PI 567731, but not in Magellan under Al stress. An F8 recombinant inbred line population derived from a cross between Magellan and PI 567731 was used to map the quantitative trait loci (QTL) for Al tolerance. Three QTL on chromosomes 3, 13, and 20, with tolerant alleles from Magellan, were identified. qAl_Gm13 and qAl_Gm20 explained large phenotypic variations (13-27%) and helped maintain root elongation and initiation under Al stress. In addition, qAl_Gm13 and qAl_Gm20 were confirmed in near-isogenic backgrounds and were identified to epistatically regulate Al tolerance via internal detoxification instead of Al3+ exclusion. Phylogenetic and pedigree analysis identified the tolerant alleles of both loci derived from the US ancestral line, A.K.[FC30761], originally from China. Our results provide novel genetic resources for breeding Al-tolerant soybean and suggest that internal detoxification contributes to soybean tolerance to excessive soil Al.

Mapping Quantitative Trait Loci for Soybean Seedling Shoot and Root Architecture Traits in an Inter-Specific Genetic Population.

Wild soybean species (Glycine soja Siebold & Zucc.) comprise a unique resource to widen the genetic base of cultivated soybean [Glycine max (L.) Merr.] for various agronomic traits. An inter-specific mapping population derived from a cross of cultivar Williams 82 and PI 483460B, a wild soybean accession, was utilized for genetic characterization of root architecture traits. The objectives of this study were to identify and characterize quantitative trait loci (QTL) for seedling shoot and root architecture traits, as well as to determine additive/epistatic interaction effects of identified QTLs. A total of 16,469 single nucleotide polymorphisms (SNPs) developed for the Illumina beadchip genotyping platform were used to construct a high resolution genetic linkage map. Among the 11 putative QTLs identified, two significant QTLs on chromosome 7 were determined to be associated with total root length (RL) and root surface area (RSA) with favorable alleles from the wild soybean parent. These seedling root traits, RL (BARC_020495_04641 ~ BARC_023101_03769) and RSA (SNP02285 ~ SNP18129_Magellan), could be potential targets for introgression into cultivated soybean background to improve both tap and lateral roots. The RL QTL region harbors four candidate genes with higher expression in root tissues: Phosphofructokinase (Glyma.07g126400), Snf7 protein (Glyma.07g127300), unknown functional gene (Glyma.07g127900), and Leucine Rich-Repeat protein (Glyma.07g127100). The novel alleles inherited from the wild soybean accession could be used as molecular markers to improve root system architecture and productivity in elite soybean lines.

The importance of slow canopy wilting in drought tolerance in soybean.

Slow canopy wilting (SW) is a water conservation trait controlled by quantitative trait loci (QTLs) in late maturity group soybeans [Glycine max (L.) Merr.]. Recently, two exotic (landraces) plant introductions (PI 567690 and PI 567731) were identified as new SW lines in early maturity groups. Here, we show that the two PIs share the same water conservation strategy of limited maximum transpiration rates as PI 416937. However, in contrast to PI 416937, the transpiration rates of these PIs were sensitive to an aquaporin inhibitor, indicating an independence between limited maximum transpiration and the lack of silver-sensitive aquaporins. Yield tests of selected recombinant inbred lines from two elite/exotic crosses provide direct evidence to support the benefit of SW in drought tolerance. Four SW QTLs mapped in a Pana×PI 567690 cross at multiple environments were found to be co-located with previous reports. Moreover, two new SW QTLs were mapped on chromosomes 6 and 10 from a Magellan×PI 567731 cross. These two QTLs explain the observed relatively large contributions of 20-30% and were confirmed in a near-isogenic background. These findings demonstrate the importance of SW in yield protection under drought and provide genetic resources for improving drought tolerance in early maturity group soybeans.

Identification of new loci for salt tolerance in soybean by high-resolution genome-wide association mapping.

BACKGROUND: Salinity is an abiotic stress that negatively affects soybean [Glycine max (L.) Merr.] seed yield. Although a major gene for salt tolerance was identified and consistently mapped to chromosome (Chr.) 3 by linkage mapping studies, it does not fully explain genetic variability for tolerance in soybean germplasm. In this study, a genome-wide association study (GWAS) was performed to map genomic regions for salt tolerance in a diverse panel of 305 soybean accessions using a single nucleotide polymorphism (SNP) dataset derived from the SoySNP50K iSelect BeadChip. A second GWAS was also conducted in a subset of 234 accessions using another 3.7 M SNP dataset derived from a whole-genome resequencing (WGRS) study. In addition, three gene-based markers (GBM) of the known gene, Glyma03g32900, on Chr. 3 were also integrated into the two datasets. Salt tolerance among soybean lines was evaluated by leaf scorch score (LSS), chlorophyll content ratio (CCR), leaf sodium content (LSC), and leaf chloride content (LCC). RESULTS: For both association studies, a major locus for salt tolerance on Chr. 3 was confirmed by a number of significant SNPs, of which three gene-based SNP markers, Salt-20, Salt14056 and Salt11655, had the highest association with all four traits studied. Also, additional genomic regions on Chrs. 1, 8, and 18 were found to be associated with various traits measured in the second GWAS using the WGRS-derived SNP dataset. CONCLUSIONS: A region identified on Chr. 8 was identified to be associated with all four traits and predicted as a new minor locus for salt tolerance in soybean. The candidate genes harbored in this minor locus may help reveal the molecular mechanism involved in salt tolerance and to improve tolerance in soybean cultivars. The significant SNPs will be useful for marker-assisted selection for salt tolerance in soybean breeding programs.

Harnessing the Potential of Forage Legumes, Alfalfa, Soybean, and Cowpea for Sustainable Agriculture and Global Food Security.

Substantial improvements in access to food and increased purchasing power are driving many people toward consuming nutrition-rich foods causing an unprecedented demand for protein food worldwide, which is expected to rise further. Forage legumes form an important source of feed for livestock and have potential to provide a sustainable solution for food and protein security. Currently, alfalfa is a commercially grown source of forage and feed in many countries. However, soybean and cowpea also have the potential to provide quality forage and fodder for animal use. The cultivation of forage legumes is under threat from changing climatic conditions, indicating the need for breeding cultivars that can sustain and acclimatize to the negative effects of climate change. Recent progress in genetic and genomic tools have facilitated the identification of quantitative trait loci and genes/alleles that can aid in developing forage cultivars through genomics-assisted breeding. Furthermore, transgenic technology can be utilized to manipulate the genetic makeup of plants to improve forage digestibility for better animal performance. In this article, we assess the genetic potential of three important legume crops, alfalfa, soybean, and cowpea in supplying quality fodder and feed for livestock. In addition, we examine the impact of climate change on forage quality and discuss efforts made in enhancing the adaptation of the plant to the abiotic stress conditions. Subsequently, we suggest the application of integrative approaches to achieve adequate forage production amid the unpredictable climatic conditions.

Dissecting genomic hotspots underlying seed protein, oil, and sucrose content in an interspecific mapping population of soybean using high-density linkage mapping.

The cultivated [Glycine max (L) Merr.] and wild [Glycine soja Siebold & Zucc.] soybean species comprise wide variation in seed composition traits. Compared to wild soybean, cultivated soybean contains low protein, high oil, and high sucrose. In this study, an interspecific population was derived from a cross between G. max (Williams 82) and G. soja (PI 483460B). This recombinant inbred line (RIL) population of 188 lines was sequenced at 0.3× depth. Based on 91 342 single nucleotide polymorphisms (SNPs), recombination events in RILs were defined, and a high-resolution bin map was developed (4070 bins). In addition to bin mapping, quantitative trait loci (QTL) analysis for protein, oil, and sucrose was performed using 3343 polymorphic SNPs (3K-SNP), derived from Illumina Infinium BeadChip sequencing platform. The QTL regions from both platforms were compared, and a significant concordance was observed between bin and 3K-SNP markers. Importantly, the bin map derived from next-generation sequencing technology enhanced mapping resolution (from 1325 to 50 Kb). A total of five, nine, and four QTLs were identified for protein, oil, and sucrose content, respectively, and some of the QTLs coincided with soybean domestication-related genomic loci. The major QTL for protein and oil were mapped on Chr. 20 (qPro_20) and suggested negative correlation between oil and protein. In terms of sucrose content, a novel and major QTL were identified on Chr. 8 (qSuc_08) and harbours putative genes involved in sugar transport. In addition, genome-wide association using 91 342 SNPs confirmed the genomic loci derived from QTL mapping. A QTL-based haplotype using whole-genome resequencing of 106 diverse soybean lines identified unique allelic variation in wild soybean that could be utilized to widen the genetic base in cultivated soybean.

A major natural genetic variation associated with root system architecture and plasticity improves waterlogging tolerance and yield in soybean.

Natural genetic variations in waterlogging tolerance are controlled by multiple genes mapped as quantitative trait loci (QTLs) in major crops, including soybean (Glycine max L.). In this research, 2 novel QTLs associated with waterlogging tolerance were mapped from an elite/exotic soybean cross. The subsequent research was focused on a major QTL (qWT_Gm03) with the tolerant allele from the exotic parent. This QTL was isolated into near-isogenic backgrounds, and its effects on waterlogging tolerance were validated in multiple environments. Fine mapping narrowed qWT_Gm03 into a genomic region of <380 Kbp excluding Rps1 gene for Phytophthora sojae resistance. The tolerant allele of qWT_Gm03 promotes root growth under nonstress conditions and favourable root plasticity under waterlogging, resulting in improved waterlogging tolerance, yield, and drought tolerance-related traits, possibly through more efficient water/nutrient uptakes. Meanwhile, involvement of auxin pathways was also identified in the regulation of waterlogging tolerance, as the genotypic differences of qWT_Gm03 in waterlogging tolerance and formation of adventitious/aerial roots can be complemented by an exogenous auxin-biosynthesis inhibitor. These findings provided genetic resources to address the urgent demand of improving waterlogging tolerance in soybean and revealed the determinant roles of root architecture and plasticity in the plant adaptation to waterlogging.

Comparative genome analysis to identify SNPs associated with high oleic acid and elevated protein content in soybean.

The objective of this study was to determine the genetic relationship between the oleic acid and protein content. The genotypes having high oleic acid and elevated protein (HOEP) content were crossed with five elite lines having normal oleic acid and average protein (NOAP) content. The selected accessions were grown at six environments in three different locations and phenotyped for protein, oil, and fatty acid components. The mean protein content of parents, HOEP, and NOAP lines was 34.6%, 38%, and 34.9%, respectively. The oleic acid concentration of parents, HOEP, and NOAP lines was 21.7%, 80.5%, and 20.8%, respectively. The HOEP plants carried both FAD2-1A (S117N) and FAD2-1B (P137R) mutant alleles contributing to the high oleic acid phenotype. Comparative genome analysis using whole-genome resequencing data identified six genes having single nucleotide polymorphism (SNP) significantly associated with the traits analyzed. A single SNP in the putative gene Glyma.10G275800 was associated with the elevated protein content, and palmitic, oleic, and linoleic acids. The genes from the marker intervals of previously identified QTL did not carry SNPs associated with protein content and fatty acid composition in the lines used in this study, indicating that all the genes except Glyma.10G278000 may be the new genes associated with the respective traits.

Genome-Wide Analysis of Grain Yield Stability and Environmental Interactions in a Multiparental Soybean Population.

Genetic improvement toward optimized and stable agronomic performance of soybean genotypes is desirable for food security. Understanding how genotypes perform in different environmental conditions helps breeders develop sustainable cultivars adapted to target regions. Complex traits of importance are known to be controlled by a large number of genomic regions with small effects whose magnitude and direction are modulated by environmental factors. Knowledge of the constraints and undesirable effects resulting from genotype by environmental interactions is a key objective in improving selection procedures in soybean breeding programs. In this study, the genetic basis of soybean grain yield responsiveness to environmental factors was examined in a large soybean nested association population. For this, a genome-wide association to performance stability estimates generated from a Finlay-Wilkinson analysis and the inclusion of the interaction between marker genotypes and environmental factors was implemented. Genomic footprints were investigated by analysis and meta-analysis using a recently published multiparent model. Results indicated that specific soybean genomic regions were associated with stability, and that multiplicative interactions were present between environments and genetic background. Seven genomic regions in six chromosomes were identified as being associated with genotype-by-environment interactions. This study provides insight into genomic assisted breeding aimed at achieving a more stable agronomic performance of soybean, and documented opportunities to exploit genomic regions that were specifically associated with interactions involving environments and subpopulations.

Genetic mapping and haplotype analysis of a locus for quantitative resistance to Fusarium graminearum in soybean accession PI 567516C.

KEY MESSAGE: A major novel quantitative disease resistance locus, qRfg_Gm06, for Fusarium graminearum was genetically mapped to chromosome 6. Genomic-assisted haplotype analysis within this region identified three putative candidate genes. Fusarium graminearum causes seed, root rot, and seedling damping-off in soybean which contributes to reduced stands and yield. A cultivar Magellan and PI 567516C were identified with low and high levels of partial resistance to F. graminearum, respectively. Quantitative disease resistance loci (QDRL) were mapped with 241 F7:8 recombinant inbred lines (RILs) derived from a cross of Magellan × PI 567516C. Phenotypic evaluation for resistance to F. graminearum used the rolled towel assay in a randomized incomplete block design. The genetic map was constructed from 927 polymorphic single nucleotide polymorphism (SNP) and simple sequence repeat (SSR) markers. One major QDRL qRfg_Gm06 was detected and mapped to chromosome 6 with a LOD score of 20.3 explaining 40.2% of the total phenotypic variation. This QDRL was mapped to a ~400 kb genomic region of the Williams 82 reference genome. Genome mining of this region identified 14 putative candidate disease resistance genes. Haplotype analysis of this locus using whole genome re-sequencing (WGRS) of 106 diverse soybean lines narrowed the list to three genes. A SNP genotyping Kompetitive allele-specific PCR (KASP) assay was designed for one of the genes and was validated in a subset of the RILs and all 106 diverse lines.

Root xylem plasticity to improve water use and yield in water-stressed soybean.

We tested the hypothesis that increasing the number of metaxylem vessels would enhance the efficiency of water uptake in soybean (Glycine max) and decrease the yield gap in water-limited environments. A panel of 41 soybean accessions was evaluated in greenhouse, rainout shelter, and rain-fed field environments. The metaxylem number influenced the internal capture of CO2 and improved stomatal conductance, enhancing water uptake/use in soybeans exposed to stress during the reproductive stage. We determined that other root anatomical features, such as cortex cell area and the percentage of stele that comprised cortical cells, also affected seed yield under similar growth parameters. Seed yield was also impacted by pod retention rates under drought stress (24-80 pods/plant). We surmise that effective biomass allocation, that is, the transport of available photosynthates to floral structures at late reproductive growth stages (R6-R7), enables yield protection under drought stress. A mesocosm study of contrasting lines for yield under drought stress and root anatomical features revealed that increases in metaxylem number as an adaptation to drought in the high-yielding lines improved root hydraulic conductivity, which reduced the metabolic cost of exploring water in deeper soil strata and enhanced water transport. This allowed the maintenance of shoot physiological processes under water-limited conditions.

Genetic diversity and genomic strategies for improving drought and waterlogging tolerance in soybeans.

Drought and its interaction with high temperature are the major abiotic stress factors affecting soybean yield and production stability. Ongoing climate changes are anticipated to intensify drought events, which will further impact crop production and food security. However, excessive water also limits soybean production. The success of soybean breeding programmes for crop improvement is dependent on the extent of genetic variation present in the germplasm base. Screening for natural genetic variation in drought- and flooding tolerance-related traits, including root system architecture, water and nitrogen-fixation efficiency, and yield performance indices, has helped to identify the best resources for genetic studies in soybean. Genomic resources, including whole-genome sequences of diverse germplasms, millions of single-nucleotide polymorphisms, and high-throughput marker genotyping platforms, have expedited gene and marker discovery for translational genomics in soybean. This review highlights the current knowledge of the genetic diversity and quantitative trait loci associated with root system architecture, canopy wilting, nitrogen-fixation ability, and flooding tolerance that contributes to the understanding of drought- and flooding-tolerance mechanisms in soybean. Next-generation mapping approaches and high-throughput phenotyping will facilitate a better understanding of phenotype-genotype associations and help to formulate genomic-assisted breeding strategies, including genomic selection, in soybean for tolerance to drought and flooding stress.

Landscape of genomic diversity and trait discovery in soybean.

Cultivated soybean [Glycine max (L.) Merr.] is a primary source of vegetable oil and protein. We report a landscape analysis of genome-wide genetic variation and an association study of major domestication and agronomic traits in soybean. A total of 106 soybean genomes representing wild, landraces, and elite lines were re-sequenced at an average of 17x depth with a 97.5% coverage. Over 10 million high-quality SNPs were discovered, and 35.34% of these have not been previously reported. Additionally, 159 putative domestication sweeps were identified, which includes 54.34 Mbp (4.9%) and 4,414 genes; 146 regions were involved in artificial selection during domestication. A genome-wide association study of major traits including oil and protein content, salinity, and domestication traits resulted in the discovery of novel alleles. Genomic information from this study provides a valuable resource for understanding soybean genome structure and evolution, and can also facilitate trait dissection leading to sequencing-based molecular breeding.

Genomic-assisted haplotype analysis and the development of high-throughput SNP markers for salinity tolerance in soybean.

Soil salinity is a limiting factor of crop yield. The soybean is sensitive to soil salinity, and a dominant gene, Glyma03g32900 is primarily responsible for salt-tolerance. The identification of high throughput and robust markers as well as the deployment of salt-tolerant cultivars are effective approaches to minimize yield loss under saline conditions. We utilized high quality (15x) whole-genome resequencing (WGRS) on 106 diverse soybean lines and identified three major structural variants and allelic variation in the promoter and genic regions of the GmCHX1 gene. The discovery of single nucleotide polymorphisms (SNPs) associated with structural variants facilitated the design of six KASPar assays. Additionally, haplotype analysis and pedigree tracking of 93 U.S. ancestral lines were performed using publically available WGRS datasets. Identified SNP markers were validated, and a strong correlation was observed between the genotype and salt treatment phenotype (leaf scorch, chlorophyll content and Na(+) accumulation) using a panel of 104 soybean lines and, an interspecific bi-parental population (F8) from PI483463 x Hutcheson. These markers precisely identified salt-tolerant/sensitive genotypes (>91%), and different structural-variants (>98%). These SNP assays, supported by accurate phenotyping, haplotype analyses and pedigree tracking information, will accelerate marker-assisted selection programs to enhance the development of salt-tolerant soybean cultivars.

Genomic-assisted phylogenetic analysis and marker development for next generation soybean cyst nematode resistance breeding.

Soybean cyst nematode (SCN, Heterodera glycines Ichinohe) is a serious soybean pest. The use of resistant cultivars is an effective approach for preventing yield loss. In this study, 19,652 publicly available soybean accessions that were previously genotyped with the SoySNP50K iSelect BeadChip were used to evaluate the phylogenetic diversity of SCN resistance genes Rhg1 and Rhg4 in an attempt to identify novel sources of resistance. The sequence information of soybean lines was utilized to develop KASPar (KBioscience Competitive Allele-Specific PCR) assays from single nucleotide polymorphisms (SNPs) of Rhg1, Rhg4, and other novel quantitative trait loci (QTL). These markers were used to genotype a diverse set of 95 soybean germplasm lines and three recombinant inbred line (RIL) populations. SNP markers from the Rhg1 gene were able to differentiate copy number variation (CNV), such as resistant-high copy (PI 88788-type), low copy (Peking-type), and susceptible-single copy (Williams 82) numbers. Similarly, markers for the Rhg4 gene were able to detect Peking-type (resistance) genotypes. The phylogenetic information of SCN resistance loci from a large set of soybean accessions and the gene/QTL specific markers that were developed in this study will accelerate SCN resistance breeding programs.

Comparative analysis of the drought-responsive transcriptome in soybean lines contrasting for canopy wilting.

Drought stress causes significant yield losses in major oil seed crops, such as soybean [Glycine max (L.) Merr]. Few soybean lines have been identified as canopy-wilting tolerant; however, the molecular mechanism conferring tolerance is not fully understood. To understand the biological process, a whole genome transcriptome analysis was performed for leaf tissues of two contrasting soybean lines: drought-susceptible (DS) Pana and drought-tolerant (DT) PI 567690. A pairwise comparison of the DS and DT lines under drought and control conditions detected 1914 and 670 genes with a greater than two-fold change in expression under drought conditions. Pairwise treatment comparison and gene enrichment analysis on the DT line showed the down-regulation of genes associated with protein binding, hydrolase activity, carbohydrate/lipid metabolism, xyloglucan endo-transglycosylases associated with cell-wall, apoplast, and chlorophyll a/b binding proteins. On the other hand, genes that were associated with the biotic stress response, ion binding and transport, the oxido-reductive process and electron carrier activity were up-regulated. Gene enrichment analysis detected UDP glucuronosyl transferase activity-encoding genes to be differentially expressed in PI 567690 under drought stress conditions. We found valuable SNPs variation in aquaporin genes of the DT line that are conserved in known slower canopy-wilting lines, this should facilitate marker-assisted selection in soybeans with improved drought tolerance.

Identification and characterization of a stachyose synthase gene controlling reduced stachyose content in soybean.

KEY MESSAGE: We identified and characterized a mutant of soybean stachyose synthase gene controlling reduced stachyose content which benefit the soybean seed composition breeding program in the future. It has been shown that in soybean, increased sucrose and reduced raffinose family oligosaccharides would have a positive impact on the world's feed industry by improving digestibility and feed efficiency. We searched for new sources of modified oligosaccharide content in a subset of the USDA Soybean Germplasm Collection and then identified plant introduction (PI) 603176A as having ultra-low stachyose content (0.5%). We identified a 33-bp deletion mutant in the putative stachyose synthase gene (STS gene, Glyma19g40550) of PI 603176A. A co-dominate indel marker was successfully developed from this 33-bp deletion area and was genetically mapped into two F 2:3 populations and a F 4:5 population, which associated with low stachyose content in the progeny lines. These observations provided strong evidence that the STS gene is responsible for stachyose biosynthesis in the soybean plant. Expression of the sts gene remained at the normal level, suggesting the loss of function in the gene is due to defective protein function. This gene-based perfect genetic marker for low stachyose content can be useful for marker-assisted selection in soybean molecular breeding programs.

Identification of quantitative trait loci underlying resistance to southern root-knot and reniform nematodes in soybean accession PI 567516C.

Soybean cyst nematode (SCN, Heterodera glycine Ichinohe), southern root-knot nematode [SRKN, Meloidogyne incognita (Kofoid and White) Chitwood] and reniform nematode (RN, Rotylenchulus reniformis Linford and Oliveira) are three important plant-parasitic pests in soybean. Previous study showed that plant introduction (PI) 567516C harbored novel quantitative trait loci (QTL) conferring SCN resistance to soybean. However, QTL underlying resistance to SRKN and RN in PI 567516C remain unknown. The objectives of this study were to identify QTL for resistance to SRKN and RN in PI 567516C. Two hundred and forty-seven F6:9 recombinant inbred lines, derived from a cross between cultivar Magellan and PI 567516C, were evaluated for resistance to SRKN and RN. Two hundred and thirty-eight simple sequence repeats and 687 single nucleotide polymorphism markers were used to construct a genetic linkage map. Three significant QTL associated with resistance to SRKN were mapped on chromosomes (Chrs.) 10, 13 and 17. Two significant QTL associated with resistance to RN were detected on Chrs. 11 and 18. Whole-genome resequencing revealed that there might be Peking-type Rhg1 in PI 567516C. Our study provides useful information to employ PI 567516C in soybean breeding in order to develop new cultivars with resistance to multiple nematodes.

Genetic variants in root architecture-related genes in a Glycine soja accession, a potential resource to improve cultivated soybean.

BACKGROUND: Root system architecture is important for water acquisition and nutrient acquisition for all crops. In soybean breeding programs, wild soybean alleles have been used successfully to enhance yield and seed composition traits, but have never been investigated to improve root system architecture. Therefore, in this study, high-density single-feature polymorphic markers and simple sequence repeats were used to map quantitative trait loci (QTLs) governing root system architecture in an inter-specific soybean mapping population developed from a cross between Glycine max and Glycine soja. RESULTS: Wild and cultivated soybean both contributed alleles towards significant additive large effect QTLs on chromosome 6 and 7 for a longer total root length and root distribution, respectively. Epistatic effect QTLs were also identified for taproot length, average diameter, and root distribution. These root traits will influence the water and nutrient uptake in soybean. Two cell division-related genes (D type cyclin and auxin efflux carrier protein) with insertion/deletion variations might contribute to the shorter root phenotypes observed in G. soja compared with cultivated soybean. Based on the location of the QTLs and sequence information from a second G. soja accession, three genes (slow anion channel associated 1 like, Auxin responsive NEDD8-activating complex and peroxidase), each with a non-synonymous single nucleotide polymorphism mutation were identified, which may also contribute to changes in root architecture in the cultivated soybean. In addition, Apoptosis inhibitor 5-like on chromosome 7 and slow anion channel associated 1-like on chromosome 15 had epistatic interactions for taproot length QTLs in soybean. CONCLUSION: Rare alleles from a G. soja accession are expected to enhance our understanding of the genetic components involved in root architecture traits, and could be combined to improve root system and drought adaptation in soybean.

Identification of novel QTL governing root architectural traits in an interspecific soybean population.

Cultivated soybean (Glycine max L.) cv. Dunbar (PI 552538) and wild G. soja (PI 326582A) exhibited significant differences in root architecture and root-related traits. In this study, phenotypic variability for root traits among 251 BC2F5 backcross inbred lines (BILs) developed from the cross Dunbar/PI 326582A were identified. The root systems of the parents and BILs were evaluated in controlled environmental conditions using a cone system at seedling stage. The G. max parent Dunbar contributed phenotypically favorable alleles at a major quantitative trait locus on chromosome 8 (Satt315-I locus) that governed root traits (tap root length and lateral root number) and shoot length. This QTL accounted for >10% of the phenotypic variation of both tap root and shoot length. This QTL region was found to control various shoot- and root-related traits across soybean genetic backgrounds. Within the confidence interval of this region, eleven transcription factors (TFs) were identified. Based on RNA sequencing and Affymetrix expression data, key TFs including MYB, AP2-EREBP and bZIP TFs were identified in this QTL interval with high expression in roots and nodules. The backcross inbred lines with different parental allelic combination showed different expression pattern for six transcription factors selected based on their expression pattern in root tissues. It appears that the marker interval Satt315-I locus on chromosome 8 contain an essential QTL contributing to early root and shoot growth in soybean.