Genome-wide association study and RNA-seq identifies GmWRI1-like transcription factor related to the seed weight in soybean.

The cultivated soybean (Glycine max (L.) Merrill) is domesticated from wild soybean (Glycine soja) and has heavier seeds with a higher oil content than the wild soybean. In this study, we identified a novel candidate gene associated with SW using a genome-wide association study (GWAS). The candidate gene GmWRI14-like was detected by GWAS analysis in three consecutive years. By constructing transgenic soybeans overexpressing the GmWRI14-like gene and gmwri14-like soybean mutants, we found that overexpression of GmWRI14-like increased the SW and increased total fatty acid content. We then used RNA-seq and qRT-PCR to identify the target genes directly or indirectly regulated by GmWRI14-like. Transgenic soyabeans overexpressing GmWRI14-like showed increased accumulation of GmCYP78A50 and GmCYP78A69 than non-transgenic soybean lines. Interestingly, we also found that GmWRI14-like proteins could interact with GmCYP78A69/GmCYP78A50 using yeast two-hybrid and bimolecular fluorescence complementation. Our results not only shed light on the genetic architecture of cultivated soybean SW, but also lays a theoretical foundation for improving the SW and oil content of soybeans.

Candidate loci for breeding compact plant-type soybean varieties.

Plant height and node number are important agronomic traits that influence yield in soybean (Glycine max L.). Here, to better understand the genetic basis of the traits, we used two recombinant inbred line (RIL) populations to detect quantitative trait loci (QTLs) associated with plant height and node number in different environments. This analysis detected 9 and 21 QTLs that control plant height and node number, respectively. Among them, we identified two genomic regions that overlap with Determinate stem 1 (Dt1) and Dt2, which are known to influence both plant height and node number. Furthermore, different combinations of Dt1 and Dt2 alleles were enriched in distinct latitudes. In addition, we determined that the QTLs qPH-13-SE and qPH-13-DW in the two RIL populations overlap with genomic intervals associated with plant height and the QTL qNN-04-DW overlaps with an interval associated with node number. Combining the dwarf allele of qPH-13-SE/qPH-13-DW and the multiple-node allele of qNN-04-DW produced plants with ideal plant architecture, i.e., shorter main stems with more nodes. This plant type may help increase yield at high planting density. This study thus provides candidate loci for breeding elite soybean cultivars for plant height and node number. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-022-01352-2.

The genetic basis of high-latitude adaptation in wild soybean.

In many plants, flowering time is influenced by daylength as an adaptive response. In soybean (Glycine max) cultivars, however, photoperiodic flowering reduces crop yield and quality in high-latitude regions. Understanding the genetic basis of wild soybean (Glycine soja) adaptation to high latitudes could aid breeding of improved cultivars. Here, we identify the Tof4 (Time of flowering 4) locus, which encodes by an E1-like protein, E1La, that represses flowering and enhances adaptation to high latitudes in wild soybean. Moreover, we found that Tof4 physically associates with the promoters of two important FLOWERING LOCUS T (FT2a and FT5a) and with Tof5 to inhibit their transcription under long photoperiods. The effect of Tof4 on flowering and maturity is mediated by FT2a and FT5a proteins. Intriguingly, Tof4 and the key flowering repressor E1 independently but additively regulate flowering time, maturity, and grain yield in soybean. We determined that weak alleles of Tof4 have undergone natural selection, facilitating adaptation to high latitudes in wild soybean. Notably, over 71.5% of wild soybean accessions harbor the mutated alleles of Tof4 or a previously reported gain-of-function allele Tof5H2, suggesting that these two loci are the genetic basis of wild soybean adaptation to high latitudes. Almost no cultivated soybean carries the mutated tof4 allele. Introgression of the tof4-1 and Tof5H2 alleles into modern soybean or editing E1 family genes thus represents promising avenues to obtain early-maturity soybean, thereby improving productivity in high latitudes.

Diverse flowering responses subjecting to ambient high temperature in soybean under short-day conditions.

Flowering time is one of important agronomic traits determining the crop yield and affected by high temperature. When facing high ambient temperature, plants often initiate early flowering as an adaptive strategy to escape the stress and ensure successful reproduction. However, here we find opposing ways in the short-day crop soybean to respond to different levels of high temperatures, in which flowering accelerates when temperature changes from 25 to 30 °C, but delays when temperature reaches 35 °C under short day. phyA-E1, possibly photoperiodic pathway, is crucial for 35 °C-mediated late flowering, however, does not contribute to promoting flowering at 30 °C. 30 °C-induced up-regulation of FT2a and FT5a leads to early flowering, independent of E1. Therefore, distinct responsive mechanisms are adopted by soybean when facing different levels of high temperatures for successful flowering and reproduction.

GIGANTEA orthologs, E2 members, redundantly determine photoperiodic flowering and yield in soybean.

Soybean (Glycine max L.) is a typical photoperiod-sensitive crop, such that photoperiod determines its flowering time, maturity, grain yield, and phenological adaptability. During evolution, the soybean genome has undergone two duplication events, resulting in about 75% of all genes being represented by multiple copies, which is associated with rampant gene redundancy. Among duplicated genes, the important soybean maturity gene E2 has two homologs, E2-Like a (E2La) and E2-Like b (E2Lb), which encode orthologs of Arabidopsis GIGANTEA (GI). Although E2 was cloned a decade ago, we still know very little about its contribution to flowering time and even less about the function of its homologs. Here, we generated single and double mutants in E2, E2La, and E2Lb by genome editing and determined that E2 plays major roles in the regulation of flowering time and yield, with the two E2 homologs depending on E2 function. At high latitude regions, e2 single mutants showed earlier flowering and high grain yield. Remarkably, in terms of genetic relationship, genes from the legume-specific transcription factor family E1 were epistatic to E2. We established that E2 and E2-like proteins form homodimers or heterodimers to regulate the transcription of E1 family genes, with the homodimer exerting a greater function than the heterodimers. In addition, we established that the H3 haplotype of E2 is the ancestral allele and is mainly restricted to low latitude regions, from which the loss-of-function alleles of the H1 and H2 haplotypes were derived. Furthermore, we demonstrated that the function of the H3 allele is stronger than that of the H1 haplotype in the regulation of flowering time, which has not been shown before. Our findings provide excellent allelic combinations for classical breeding and targeted gene disruption or editing.

Novel and multifaceted regulations of photoperiodic flowering by phytochrome A in soybean.

Photoperiod is an important environmental cue. Plants can distinguish the seasons and flower at the right time through sensing the photoperiod. Soybean is a sensitive short-day crop, and the timing of flowering varies greatly at different latitudes, thus affecting yields. Soybean cultivars in high latitudes adapt to the long day by the impairment of two phytochrome genes, PHYA3 and PHYA2, and the legume-specific flowering suppressor, E1. However, the regulating mechanism underlying phyA and E1 in soybean remains largely unknown. Here, we classified the regulation of the E1 family by phyA2 and phyA3 at the transcriptional and posttranscriptional levels, revealing that phyA2 and phyA3 regulate E1 by directly binding to LUX proteins, the critical component of the evening complex, to regulate the stability of LUX proteins. In addition, phyA2 and phyA3 can also directly associate with E1 and its homologs to stabilize the E1 proteins. Therefore, phyA homologs control the core flowering suppressor E1 at both the transcriptional and posttranscriptional levels, to double ensure the E1 activity. Thus, our results disclose a photoperiod flowering mechanism in plants by which the phytochrome A regulates LUX and E1 activity.

Agronomical selection on loss-of-function of GIGANTEA simultaneously facilitates soybean salt tolerance and early maturity.

Salt stress and flowering time are major factors limiting geographic adaptation and yield productivity in soybean (Glycine max). Although improving crop salt tolerance and latitude adaptation are essential for efficient agricultural production, whether and how these two traits are integrated remains largely unknown. Here, we used a genome-wide association study to identify a major salt-tolerance locus controlled by E2, an ortholog of Arabidopsis thaliana GIGANTEA (GI). Loss of E2 function not only shortened flowering time and maturity, but also enhanced salt-tolerance in soybean. E2 delayed soybean flowering by enhancing the transcription of the core flowering suppressor gene E1, thereby repressing Flowering Locus T (FT) expression. An E2 knockout mutant e2CR displayed reduced accumulation of reactive oxygen species (ROS) during the response to salt stress by releasing peroxidase, which functions in ROS scavenging to avoid cytotoxicity. Evolutionary and population genetic analyses also suggested that loss-of-function e2 alleles have been artificially selected during breeding for soybean adaptation to high-latitude regions with greater salt stress. Our findings provide insights into the coupled selection for adaptation to both latitude and salt stress in soybean; and offer an ideal target for molecular breeding of early-maturing and salt-tolerant cultivars.

A functionally divergent SOC1 homolog improves soybean yield and latitudinal adaptation.

Soybean (Glycine max) grows in a wide range of latitudes, but it is extremely sensitive to photoperiod, which reduces its yield and ability to adapt to different environments. Therefore, understanding of the genetic basis of soybean adaptation is of great significance for breeding and improvement. Here, we characterized Tof18 (SOC1a) that conditions early flowering and growth habit under both short-day and long-day conditions. Molecular analysis confirmed that the two SOC1 homologs present in soybeans (SOC1a and SOC1b) underwent evolutionary functional divergence, with SOC1a having stronger effects on flowering time and stem node number than SOC1b due to transcriptional differences. soc1a soc1b double mutants showed stronger functional effects than either of the single mutants, perhaps due to the formation of SOC1a and SOC1b homodimers or heterodimers. Additionally, Tof18/SOC1a improves the latitudinal adaptation of cultivated soybeans, highlighting the functional importance of SOC1a. The Tof18G allele facilitates adaptation to high latitudes, whereas Tof18A facilitates adaptation to low latitudes. We demonstrated that SOC1s contribute to floral induction in both leaves and shoot apex through inter-regulation with FTs. The SOC1a-SOC1b-Dt2 complex plays essential roles in stem growth habit by directly binding to the regulatory sequence of Dt1, making the genes encoding these proteins potential targets for genome editing to improve soybean yield via molecular breeding. Since the natural Tof18A allele increases node number, introgressing this allele into modern cultivars could improve yields, which would help optimize land use for food production in the face of population growth and global warming.

Identifying Quantitative Trait Loci and Candidate Genes Conferring Resistance to Soybean Mosaic Virus SC7 by Quantitative Trait Loci-Sequencing in Soybean.

Soybean mosaic virus (SMV) is detrimental to soybean (Glycine max) breeding, seed quality, and yield worldwide. Improving the basic resistance of host plants is the most effective and economical method to reduce damage from SMV. Therefore, it is necessary to identify and clone novel SMV resistance genes. Here, we report the characterization of two soybean cultivars, DN50 and XQD, with different levels of resistance to SMV. Compared with XQD, DN50 exhibits enhanced resistance to the SMV strain SC7. By combining bulked-segregant analysis (BSA)-seq and fine-mapping, we identified a novel resistance locus, R SMV -11, spanning an approximately 207-kb region on chromosome 11 and containing 25 annotated genes in the reference Williams 82 genome. Of these genes, we identified eleven with non-synonymous single-nucleotide polymorphisms (SNPs) or insertion-deletion mutations (InDels) in their coding regions between two parents. One gene, GmMATE68 (Glyma.11G028900), harbored a frameshift mutation. GmMATE68 encodes a multidrug and toxic compound extrusion (MATE) transporter that is expressed in all soybean tissues and is induced by SC7. Given that MATE transporter families have been reported to be linked with plant disease resistance, we suggest that GmMATE68 is responsible for SC7 resistance in DN50. Our results reveal a novel SMV-resistance locus, improving understanding of the genetics of soybean disease resistance and providing a potential new tool for marker-assisted selection breeding in soybean.

Genome-Wide Analysis of DREB Genes Identifies a Novel Salt Tolerance Gene in Wild Soybean (Glycine soja).

Salt stress is a major factor limiting the growth and yield of soybean (Glycine max). Wild soybeans (Glycine soja) contain high allelic diversity and beneficial alleles that can be re-introduced into domesticated soybeans to improve adaption to the environment. However, very few beneficial alleles have been identified from wild soybean. Here, we demonstrate that wild soybean is more salt tolerant than cultivated soybean and examine dehydration responsive element-binding (DREB) family transcription factor genes to look for advantageous alleles that might improve drought tolerance in cultivated soybean. Our genome-wide analysis identified 103 DREB genes from the Glycine max genome. By combined RNA-sequencing and population genetics of wild, landrace, and cultivated soybean accessions, we show that the natural variation in DREB3a and DREB3b is related to differences in salt tolerance in soybean accessions. Interestingly, DREB3b, but not DREB3a, appears to have undergone artificial selection. Soybean plants carrying the wild soybean DREB3b allele (DREB3b39Del ) are more salt tolerant than those containing the reference genome allele (DREB3bRef ). Together, our results suggest that the loss of the DREB3b39Del allele through domestication of cultivated soybean may be associated with a reduction in salt tolerance. Our findings provide crucial information for improving salt tolerance in soybean through molecular breeding.

Parallel selection of distinct Tof5 alleles drove the adaptation of cultivated and wild soybean to high latitudes.

Photoperiod responsiveness is a key factor limiting the geographic distribution of cultivated soybean and its wild ancestor. In particular, the genetic basis of the adaptation in wild soybean remains poorly understood. In this study, by combining whole-genome resequencing and genome-wide association studies we identified a novel locus, Time of Flowering 5 (Tof5), which promotes flowering and enhances adaptation to high latitudes in both wild and cultivated soybean. By genomic, genetic and transgenic analyses we showed that Tof5 encodes a homolog of Arabidopsis thaliana FRUITFULL (FUL). Importantly, further analyses suggested that different alleles of Tof5 have undergone parallel selection. The Tof5H1 allele was strongly selected by humans after the early domestication of cultivated soybean, while Tof5H2 allele was naturally selected in wild soybean, and in each case facilitating adaptation to high latitudes. Moreover, we found that the key flowering repressor E1 suppresses the transcription of Tof5 by binding to its promoter. In turn, Tof5 physically associates with the promoters of two important FLOWERING LOCUS T (FT), FT2a and FT5a, to upregulate their transcription and promote flowering under long photoperiods. Collectively, our findings provide insights into how wild soybean adapted to high latitudes through natural selection and indicate that cultivated soybean underwent changes in the same gene but evolved a distinct allele that was artificially selected after domestication.

The legume-specific transcription factor E1 controls leaf morphology in soybean.

BACKGROUND: The leaf is a determinate organ essential for photosynthesis, whose size and shape determine plant architecture and strongly affect agronomic traits. In soybean, the molecular mechanism of leaf development is not well understood. The flowering repressor gene E1, which encodes a legume-specific B3-like protein, is known to be the gene with the largest influence on soybean flowering and maturity. However, knowledge of its potential other functions remains poor. RESULTS: Here, we identified a novel function of E1 protein in leaf development. Unifoliolate leaves of E1-overexpression (E1-OE) lines were smaller and curlier than those of wild type DongNong 50 (DN50) and Williams 82 (W82). Transverse histological sections showed disorganized cells and significantly elevated palisade tissue number, spongy tissue number, and bulliform cell number in E1-OE lines. Our results indicate that E1 binds to the promoters of the leaf- development-related CINCINNATA (CIN)-like TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTOR (TCP) transcription factor genes to negatively regulate their expression. CONCLUSIONS: Our findings identify E1 as an important new factor in soybean leaf development.

GmRAV confers ecological adaptation through photoperiod control of flowering time and maturity in soybean.

Photoperiod strictly controls vegetative and reproductive growth stages in soybean (Glycine max). A soybean GmRAV (Related to ABI3/VP1) transcription factor containing both AP2 and B3 domains was shown to be a key component of this process. We identified six polymorphisms in the GmRAV promoter that showed significant association with flowering time and maturity of soybean in one or multiple environments. Soybean varieties with minor polymorphism exhibited a longer growth period contributing to soybean adaptation to lower latitudes. The cis-acting element GT1CONSENSUS motif of the GmRAV promoter controlled the growth period, and the major allele in this motif shortened duration of late reproductive stages by reducing GmRAV expression levels. Three GmRAV-overexpressing (GmRAV-ox) transgenic lines displayed later flowering time and maturity, shorter height and fewer numbers of leaves compared with control plants, whereas transgenic inhibition of GmRAV expression resulted in earlier flowering time and maturity and increased plant height. Combining DNA affinity purification sequencing and RNA sequencing analyses revealed 154 putative target genes directly bound and transcriptionally regulated by GmRAV. Two GmRAV binding motifs [C(A/G)AACAA(G/T)A(C/T)A(G/T)] and [C(T/A)A(C)C(T/G)CTG] were identified, and acting downstream of E3E4, GmRAV repressed GmFT5a transcriptional activity through binding a CAACA motif, thereby delaying soybean growth and extending both vegetative and reproductive phases.

Genetic basis and adaptation trajectory of soybean from its temperate origin to tropics.

Soybean (Glycine max) serves as a major source of protein and edible oils worldwide. The genetic and genomic bases of the adaptation of soybean to tropical regions remain largely unclear. Here, we identify the novel locus Time of Flowering 16 (Tof16), which confers delay flowering and improve yield at low latitudes and determines that it harbors the soybean homolog of LATE ELONGATED HYPOCOTYL (LHY). Tof16 and the previously identified J locus genetically additively but independently control yield under short-day conditions. More than 80% accessions in low latitude harbor the mutations of tof16 and j, which suggests that loss of functions of Tof16 and J are the major genetic basis of soybean adaptation into tropics. We suggest that maturity and yield traits can be quantitatively improved by modulating the genetic complexity of various alleles of the LHY homologs, J and E1. Our findings uncover the adaptation trajectory of soybean from its temperate origin to the tropics.

Overcoming the genetic compensation response of soybean florigens to improve adaptation and yield at low latitudes.

The classical soybean (Glycine max) trait long juvenile (LJ) is essentially a reduction in sensitivity to short-day (SD) conditions for induction and completion of flowering, and has been introduced into soybean cultivars to improve yield in tropical environments. However, only one locus, J, is known to confer LJ in low-latitude varieties. Here, we defined two quantitative trait loci contributing to the LJ trait, LJ16.1 and LJ16.2, and identified them as the florigen (FT) homologs FT2a and FT5a, respectively. The two selected florigen variations both delay flowering time under SD conditions by repressing the floral meristem identity gene GmAPETALA1. Single mutants have a relatively subtle effect on flowering time and displayed a substantial genetic compensation response, but this was absent in ft2a ft5a double mutants, which showed an enhanced LJ phenotype that translated to higher yields under SD conditions. A survey of sequence diversity suggests that FT2a and FT5a variants have diverse origins and have played distinct roles as soybean spread to lower latitudes. Our results show that integration of variants in the florigen genes offers a strategy for customizing flowering time to adjust adaptation and improve crop productivity in tropical regions.

A critical role of the soybean evening complex in the control of photoperiod sensitivity and adaptation.

Photoperiod sensitivity is a key factor in plant adaptation and crop production. In the short-day plant soybean, adaptation to low latitude environments is provided by mutations at the J locus, which confer extended flowering phase and thereby improve yield. The identity of J as an ortholog of Arabidopsis ELF3, a component of the circadian evening complex (EC), implies that orthologs of other EC components may have similar roles. Here we show that the two soybean homeologs of LUX ARRYTHMO interact with J to form a soybean EC. Characterization of mutants reveals that these genes are highly redundant in function but together are critical for flowering under short day, where the lux1 lux2 double mutant shows extremely late flowering and a massively extended flowering phase. This phenotype exceeds that of any soybean flowering mutant reported to date, and is strongly reminiscent of the "Maryland Mammoth" tobacco mutant that featured in the seminal 1920 study of plant photoperiodism by Garner and Allard [W. W. Garner, H. A. Allard, J. Agric. Res. 18, 553-606 (1920)]. We further demonstrate that the J-LUX complex suppresses transcription of the key flowering repressor E1 and its two homologs via LUX binding sites in their promoters. These results indicate that the EC-E1 interaction has a central role in soybean photoperiod sensitivity, a phenomenon also first described by Garner and Allard. EC and E1 family genes may therefore constitute key targets for customized breeding of soybean varieties with precise flowering time adaptation, either by introgression of natural variation or generation of new mutants by gene editing.

A flowering time locus dependent on E2 in soybean.

Soybean [Glycine max (L.) Merrill] is very sensitive to changes in photoperiod as a typical short-day plant. Photoperiodic flowering influences soybean latitudinal adaptability and yield to a considerable degree. Identifying new quantitative trait loci (QTLs) controlling flowering time is a powerful initial approach for elucidating the mechanisms underlying flowering time and adaptation to different latitudes in soybean. In this study, we developed a Recombinant Inbred Lines (RILs) population and recorded flowering time under natural long-day conditions. We also constructed a high-density genetic map by genotyping-by-sequencing and used it for QTL mapping. In total, we detected twelve QTLs, four of which are stable and named by qR1-2, qR1-4, qR1-6.1, and qR1-10, respectively. Among these four QTLs, qR1-4 and qR1-6.1 are novel. QTL mapping in two sub-populations classified by the genotype of the maturity locus E2, genetic interaction evaluation between E2 and qR1-2, and qRT-PCR indicated that E2 has an epistatic effect on qR1-2, and that causal gene of qR1-2 acts upstream of E2. We presumed the most likely candidate genes according to the resequencing data and briefly analyzed the geographic distributions of these genes. These findings will be beneficial for our understanding of the mechanisms underlying photoperiodic flowering in soybean, contribute to further investigate of E2, and provide genetic resources for molecular breeding of soybean. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-021-01224-1.

A recent retrotransposon insertion of J caused E6 locus facilitating soybean adaptation into low latitude.

Soybean (Glycine max) is an important legume crop that was domesticated in temperate regions. Soybean varieties from these regions generally mature early and exhibit extremely low yield when grown under inductive short-day (SD) conditions at low latitudes. The long-juvenile (LJ) trait, which is characterized by delayed flowering and maturity, and improved yield under SD conditions, allowed the cultivation of soybean to expand to lower latitudes. Two major loci control the LJ trait: J and E6. In the current study, positional cloning, sequence analysis, and transgenic complementation confirmed that E6 is a novel allele of J, the ortholog of Arabidopsis thaliana EARLY FLOWERING 3 (ELF3). The mutant allele e6PG , which carries a Ty1/Copia-like retrotransposon insertion, does not suppress the legume-specific flowering repressor E1, allowing E1 to inhibit Flowering Locus T (FT) expression and thus delaying flowering and increasing yields under SD conditions. The e6PG allele is a rare allele that has not been incorporated into modern breeding programs. The dysfunction of J might have greatly facilitated the adaptation of soybean to low latitudes. Our findings increase our understanding of the molecular mechanisms underlying the LJ trait and provide valuable resources for soybean breeding.

Two homologous LHY pairs negatively control soybean drought tolerance by repressing the abscisic acid responses.

The circadian clock plays essential roles in diverse plant biological processes, such as flowering, phytohormone biosynthesis and abiotic stress responses. The manner in which circadian clock genes regulate drought stress responses in model plants has been well established, but comparatively little is known in crop species, such as soybean, a major global crop. This paper reports that the core clock components GmLHYs, the orthologues of CCA1/LHY in Arabidopsis, negatively control drought tolerance in soybean. The expressions of four GmLHYs were all induced by drought, and the quadruple mutants of GmLHYs demonstrated significantly improved drought tolerance. Transcriptome profiling suggested that the abscisic acid (ABA) signaling pathway is regulated by GmLHYs to respond to drought tolerance. Genetic dissections showed that two homologous pairs of LHY1a and LHY1b redundantly control the drought response. Functional characterization of LHY1a and LHY1b in Arabidopsis and soybean further supported the notion that GmLHYs can maintain cellular homeostasis through the ABA signaling pathway under drought stress. This study improves our understanding of the underlying molecular mechanisms on soybean drought tolerance. Furthermore, the two homologues of LHY1a and LHY1b provide alternative targets for genome editing to rapidly generate mutant alleles in elite soybean cultivars to enhance their drought tolerance.

High-density QTL mapping of leaf-related traits and chlorophyll content in three soybean RIL populations.

BACKGROUND: Leaf size and shape, which affect light capture, and chlorophyll content are important factors affecting photosynthetic efficiency. Genetic variation of these components significantly affects yield potential and seed quality. Identification of the genetic basis for these traits and the relationship between them is of great practical significance for achieving ideal plant architecture and high photosynthetic efficiency for improved yield. RESULTS: Here, we undertook a large-scale linkage mapping study using three mapping populations to determine the genetic interplay between soybean leaf-related traits and chlorophyll content across two environments. Correlation analysis revealed a significant negative correlation between leaf size and shape, while both traits were positively correlated with chlorophyll content. This phenotypic relationship was verified across the three mapping populations as determined by principal component analysis, suggesting that these traits are under the control of complex and interrelated genetic components. The QTLs for leaf-related traits and chlorophyll are partly shared, which further supports the close genetic relationship between the two traits. The largest-effect major loci, q20, was stably identified across all population and environments and harbored the narrow leaflet gene Gm-JAG1 (Ln/ln), which is a key regulator of leaflet shape in soybean. CONCLUSION: Our results uncover several major QTLs (q4-1, q4-2, q11, q13, q18 and q20) and its candidate genes specific or common to leaf-related traits and chlorophyll, and also show a complex epistatic interaction between the two traits. The SNP markers closely linked to these valuable QTLs could be used for molecular design breeding with improved plant architecture, photosynthetic capacity and even yield.