SUnSeT: spectral unmixing of hyperspectral images for phenotyping soybean seed traits.

KEY MESSAGE: Hyperspectral features enable accurate classification of soybean seeds using linear discriminant analysis and GWAS for novel seed trait genes. Evaluating crop seed traits such as size, shape, and color is crucial for assessing seed quality and improving agricultural productivity. The introduction of the SUnSet toolbox, which employs hyperspectral sensor-derived image analysis, addresses this necessity. In a validation test involving 420 seed accessions from the Korean Soybean Core Collections, the pixel purity index algorithm identified seed- specific hyperspectral endmembers to facilitate segmentation. Various metrics extracted from ventral and lateral side images facilitated the categorization of seeds into three size groups and four shape groups. Additionally, quantitative RGB triplets representing seven seed coat colors, averaged reflectance spectra, and pigment indices were acquired. Machine learning models, trained on a dataset comprising 420 accession seeds and 199 predictors encompassing seed size, shape, and reflectance spectra, achieved accuracy rates of 95.8% for linear discriminant analysis model. Furthermore, a genome-wide association study utilizing hyperspectral features uncovered associations between seed traits and genes governing seed pigmentation and shapes. This comprehensive approach underscores the effectiveness of SUnSet in advancing precision agriculture through meticulous seed trait analysis.

The discovery of novel SNPs associated with group A soyasaponin biosynthesis from Korea soybean core collection.

Soyasaponin is a type of glycoside such as steroids, steroidal alkaloids or triterpenes, which enhance the body immunity. In order to efficiently identify genes and markers related to the soyasaponin, we used a 180K Axiom® SoyaSNP array and whole genome resequencing data from the Korean soybean core collection. As a result of conducting GWAS for group A soyasaponin (Aa and Ab derivatives), 16 significant common markers associated with Aa and Ab derivatives were mapped to chromosome 7, and three candidate genes including Glyma.07g254600 were detected. The functional haplotypes for candidate genes showed that Aa and Ab contents were mainly determined by alleles of AX-90322128, the marker of Glyma.07g254600. In addition, 14 novel SNPs variants closely associated with Aa and Ab derivatives were discovered for Glyma.07g254600. Therefore, the results of this study that identified soyasaponin-associated markers and useful genes utilizing various genomic information could provide insight into functional soybean breeding.

Quantitative Trait Loci (QTL) Analysis of Seed Protein and Oil Content in Wild Soybean (Glycine soja).

Soybean seeds consist of approximately 40% protein and 20% oil, making them one of the world's most important cultivated legumes. However, the levels of these compounds are negatively correlated with each other and regulated by quantitative trait loci (QTL) that are controlled by several genes. In this study, a total of 190 F2 and 90 BC1F2 plants derived from a cross of Daepung (Glycine max) with GWS-1887 (G. soja, a source of high protein), were used for the QTL analysis of protein and oil content. In the F2:3 populations, the average protein and oil content was 45.52% and 11.59%, respectively. A QTL associated with protein levels was detected at Gm20_29512680 on chr. 20 with a likelihood of odds (LOD) of 9.57 and an R2 of 17.2%. A QTL associated with oil levels was also detected at Gm15_3621773 on chr. 15 (LOD: 5.80; R2: 12.2%). In the BC1F2:3 populations, the average protein and oil content was 44.25% and 12.14%, respectively. A QTL associated with both protein and oil content was detected at Gm20_27578013 on chr. 20 (LOD: 3.77 and 3.06; R2 15.8% and 10.7%, respectively). The crossover to the protein content of BC1F3:4 population was identified by SNP marker Gm20_32603292. Based on these results, two genes, Glyma.20g088000 (S-adenosyl-l-methionine-dependent methyltransferases) and Glyma.20g088400 (oxidoreductase, 2-oxoglutarate-Fe(II) oxygenase family protein), in which the amino acid sequence had changed and a stop codon was generated due to an InDel in the exon region, were identified.

Effect of Soybean Volatiles on the Behavior of the Bean Bug, Riptortus pedestris.

The bean bug, Riptortus pedestris, is a polyphagous insect that feeds primarily on leguminous plants, especially soybean (Glycine max). Although the bean bug is an economically important pest of soybean, little is known about how the insect locates soybean fields. In this study, we examined the electroantennogram responses of R. pedestris to soybean volatiles and examined the behavioral responses of the adult bean bugs. R. pedestris adults were attracted more to their host-plant soybean, even when physical contact was absent, than to air or a non-host plant. Accordingly, we hypothesized that R. pedestris can recognize soybean through a plant's volatile organic compounds (VOCs). Five VOCs were identified from intact soybean plants at the vegetative stage: (Z)-3-hexen-1-ol, (Z)-3-hexenyl acetate, 4-ethylbenzaldehyde, α-farnesene, and methyl salicylate. Response spectra of the antennae to these volatiles clearly showed that both male and female R. pedestris can detect soybean volatiles. The adult bean bugs did not show behavioral orientation to any individual compounds but showed significant orientation to a particular blend of synthetic soybean volatiles when tested under laboratory conditions. In the field, this soybean volatile blend did not significantly attract the bean bugs, but it did interact synergistically with the aggregation pheromone to attract the bean bugs. These results highlight the role of host plant volatiles in the sensory ecology of R. pedestris and help explain colonization pattern of the bean bugs in soybean fields.

Chromosomal features revealed by comparison of genetic maps of Glycine max and Glycine soja.

Recombination is a crucial component of evolution and breeding. New combinations of variation on chromosomes are shaped by recombination. Recombination is also involved in chromosomal rearrangements. However, recombination rates vary tremendously among chromosome segments. Genome-wide genetic maps are one of the best tools to study variation of recombination. Here, we describe high density genetic maps of Glycine max and Glycine soja constructed from four segregating populations. The maps were used to identify chromosomal rearrangements and find the highly predictable pattern of cross-overs on the broad scale in soybean. Markers on these genetic maps were used to evaluate assembly quality of the current soybean reference genome sequence. We find a strong inversion candidate larger than 3 Mb based on patterns of cross-overs. We also identify quantitative trait loci (QTL) that control number of cross-overs. This study provides fundamental insights relevant to practical strategy for breeding programs and for pan-genome researches.

High Throughput Phenotyping for Various Traits on Soybean Seeds Using Image Analysis.

Data phenotyping traits on soybean seeds such as shape and color has been obscure because it is difficult to define them clearly. Further, it takes too much time and effort to have sufficient number of samplings especially length and width. These difficulties prevented seed morphology to be incorporated into efficient breeding program. Here, we propose methods for an image acquisition, a data processing, and analysis for the morphology and color of soybean seeds by high-throughput method using images analysis. As results, quantitative values for colors and various types of morphological traits could be screened to create a standard for subsequent evaluation of the genotype. Phenotyping method in the current study could define the morphology and color of soybean seeds in highly accurate and reliable manner. Further, this method enables the measurement and analysis of large amounts of plant seed phenotype data in a short time, which was not possible before. Fast and precise phenotype data obtained here may facilitate Genome Wide Association Study for the gene function analysis as well as for development of the elite varieties having desirable seed traits.

Genetic diversity patterns and domestication origin of soybean.

KEY MESSAGE: Genotyping data of a comprehensive Korean soybean collection obtained using a large SNP array were used to clarify global distribution patterns of soybean and address the evolutionary history of soybean. Understanding diversity and evolution of a crop is an essential step to implement a strategy to expand its germplasm base for crop improvement research. Accessions intensively collected from Korea, which is a small but central region in the distribution geography of soybean, were genotyped to provide sufficient data to underpin population genetic questions. After removing natural hybrids and duplicated or redundant accessions, we obtained a non-redundant set comprising 1957 domesticated and 1079 wild accessions to perform population structure analyses. Our analysis demonstrates that while wild soybean germplasm will require additional sampling from diverse indigenous areas to expand the germplasm base, the current domesticated soybean germplasm is saturated in terms of genetic diversity. We then showed that our genome-wide polymorphism map enabled us to detect genetic loci underlying flower color, seed-coat color, and domestication syndrome. A representative soybean set consisting of 194 accessions was divided into one domesticated subpopulation and four wild subpopulations that could be traced back to their geographic collection areas. Population genomics analyses suggested that the monophyletic group of domesticated soybeans was likely originated at a Japanese region. The results were further substantiated by a phylogenetic tree constructed from domestication-associated single nucleotide polymorphisms identified in this study.

GmBRC1 is a Candidate Gene for Branching in Soybean (Glycine max (L.) Merrill).

Branch number is one of the main factors affecting the yield of soybean (Glycine max (L.)). In this study, we conducted a genome-wide association study combined with linkage analysis for the identification of a candidate gene controlling soybean branching. Five quantitative trait nucleotides (QTNs) were associated with branch numbers in a soybean core collection. Among these QTNs, a linkage disequilibrium (LD) block qtnBR6-1 spanning 20 genes was found to overlap a previously identified major quantitative trait locus qBR6-1. To validate and narrow down qtnBR6-1, we developed a set of near-isogenic lines (NILs) harboring high-branching (HB) and low-branching (LB) alleles of qBR6-1, with 99.96% isogenicity and different branch numbers. A cluster of single nucleotide polymorphisms (SNPs) segregating between NIL-HB and NIL-LB was located within the qtnBR6-1 LD block. Among the five genes showing differential expression between NIL-HB and NIL-LB, BRANCHED1 (BRC1; Glyma.06G210600) was down-regulated in the shoot apex of NIL-HB, and one missense mutation and two SNPs upstream of BRC1 were associated with branch numbers in 59 additional soybean accessions. BRC1 encodes TEOSINTE-BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTORS 1 and 2 transcription factor and functions as a regulatory repressor of branching. On the basis of these results, we propose BRC1 as a candidate gene for branching in soybean.

Development, validation and genetic analysis of a large soybean SNP genotyping array.

Cultivated soybean (Glycine max) suffers from a narrow germplasm relative to other crop species, probably because of under-use of wild soybean (Glycine soja) as a breeding resource. Use of a single nucleotide polymorphism (SNP) genotyping array is a promising method for dissecting cultivated and wild germplasms to identify important adaptive genes through high-density genetic mapping and genome-wide association studies. Here we describe a large soybean SNP array for use in diversity analyses, linkage mapping and genome-wide association analyses. More than four million high-quality SNPs identified from high-depth genome re-sequencing of 16 soybean accessions and low-depth genome re-sequencing of 31 soybean accessions were used to select 180,961 SNPs for creation of the Axiom(®) SoyaSNP array. Validation analysis for a set of 222 diverse soybean lines showed that 170,223 markers were of good quality for genotyping. Phylogenetic and allele frequency analyses of the validation set data indicated that accessions showing an intermediate morphology between cultivated and wild soybeans collected in Korea were natural hybrids. More than 90 unanchored scaffolds in the current soybean reference sequence were assigned to chromosomes using this array. Finally, dense average spacing and preferential distribution of the SNPs in gene-rich chromosomal regions suggest that this array may be suitable for genome-wide association studies of soybean germplasm. Taken together, these results suggest that use of this array may be a powerful method for soybean genetic analyses relating to many aspects of soybean breeding.

Identification of haplotypes at the Rsv4 genomic region in soybean associated with durable resistance to soybean mosaic virus.

KEY MESSAGE: Discovery of new germplasm sources and identification of haplotypes for the durable Soybean mosaic virus resistance gene, Rsv 4, provide novel resources for map-based cloning and genetic improvement efforts in soybean. The Soybean mosaic virus (SMV) resistance locus Rsv4 is of interest because it provides a durable type of resistance in soybean [Glycine max (L.) Merr.]. To better understand its molecular basis, we used a population of 309 BC3F2 individuals to fine-map Rsv4 to a ~120 kb interval and leveraged this genetic information in a second study to identify accessions 'Haman' and 'Ilpumgeomjeong' as new sources of Rsv4. These two accessions along with three other Rsv4 and 14 rsv4 accessions were used to examine the patterns of nucleotide diversity at the Rsv4 region based on high-depth resequencing data. Through a targeted association analysis of these 19 accessions within the ~120 kb interval, a cluster of four intergenic single-nucleotide polymorphisms (SNPs) was found to perfectly associate with SMV resistance. Interestingly, this ~120 kb interval did not contain any genes similar to previously characterized dominant disease resistance genes. Therefore, a haplotype analysis was used to further resolve the association signal to a ~94 kb region, which also resulted in the identification of at least two Rsv4 haplotypes. A haplotype phylogenetic analysis of this region suggests that the Rsv4 locus in G. max is recently introgressed from G. soja. This integrated study provides a strong foundation for efforts focused on the cloning of this durable virus resistance gene and marker-assisted selection of Rsv4-mediated SMV resistance in soybean breeding programs.

Ln is a key regulator of leaflet shape and number of seeds per pod in soybean.

Narrow leaflet soybean (Glycine max) varieties tend to have more seeds per pod than broad leaflet varieties. Narrow leaflet in soybean is conferred by a single recessive gene, ln. Here, we show that the transition from broad (Ln) to narrow leaflet (ln) is associated with an amino acid substitution in the EAR motif encoded by a gene (designated Gm-JAGGED1) homologous to Arabidopsis JAGGED (JAG) that regulates lateral organ development and the variant exerts a pleiotropic effect on fruit patterning. The genomic region that regulates both the traits was mapped to a 12.6-kb region containing only one gene, Gm-JAG1. Introducing the Gm-JAG1 allele into a loss-of-function Arabidopsis jagged mutant partially restored the wild-type JAG phenotypes, including leaf shape, flower opening, and fruit shape, but the Gm-jag1 (ln) and EAR-deleted Gm-JAG1 alleles in the jagged mutant did not result in an apparent phenotypic change. These observations indicate that despite some degree of functional change of Gm-JAG1 due to the divergence from Arabidopsis JAG, Gm-JAG1 complemented the functions of JAG in Arabidopsis thaliana. However, the Gm-JAG1 homoeolog, Gm-JAG2, appears to be sub- or neofunctionalized, as revealed by the differential expression of the two genes in multiple plant tissues, a complementation test, and an allelic analysis at both loci.

Sequence variability in the HC-Pro coding regions of Korean soybean mosaic virus isolates is associated with differences in RNA silencing suppression.

Soybean mosaic virus (SMV), a member of the family Potyviridae, is an important viral pathogen affecting soybean production in Korea. Variations in helper component proteinase (HC-Pro) sequences and the pathogenicity of SMV samples from seven Korean provinces were compared with those of previously characterized SMV isolates from China, Korea and the United States. Phylogenetic analysis separated 16 new Korean SMV isolates into two groups. Fourteen of the new Korean SMV samples belonged to group II and were very similar to U.S. strain SMV G7 and Chinese isolate C14. One isolate in group II, A297-13, differed at three amino acid positions (L54F, N286D, D369N) in the HC-Pro coding sequence from severe isolates and SMV 413, showed very weak silencing suppressor activity, and produced only mild symptoms in soybean. To test the role of each amino acid substitution in RNA silencing and viral RNA accumulation, a series of point mutations was constructed. Substitution of N for D at position 286 in HC-Pro of SMV A297-12 significantly reduced silencing suppression activity. When the mutant HC-Pro of A297-13 was introduced into an infectious clone of SMV 413, accumulation of viral RNA was reduced to less than 3 % of the level of SMV 413 containing HC-Pro of A297-12 at 10 days post-inoculation (dpi) but increased to 40 % of SMV 413(HC-Pro A297-12) at 40 dpi. At 50 dpi RNA accumulation of SMV 413(HC-Pro A297-13) was similar to that of SMV 413(HC-Pro A297-12). However, at 50 dpi, the D at position 286 of HC-Pro in SMV 413(HC-Pro A297-13) was found to have reverted to N. The results showed that 1) a naturally occurring mutation in HC-Pro significantly reduced silencing suppression activity and accumulation of transgene and viral RNAs, and 2) that there was strong selection for revision to wild type when the mutation was introduced into an infectious clone of SMV.

Genome-wide SNP discovery in mungbean by Illumina HiSeq.

Mungbean [Vigna radiata (L.) Wilczek], a self-pollinated diploid plant with 2n = 22 chromosomes, is an important legume crop with a high-quality amino acid profile. Sequence variation at the whole-genome level was examined by comparing two mungbean cultivars, Sunhwanokdu and Gyeonggijaerae 5, using Illumina HiSeq sequencing data. More than 40 billion bp from both mungbean cultivars were sequenced to a depth of 72×. After de novo assembly of Sunhwanokdu contigs by ABySS 1.3.2 (N50 = 9,958 bp), those longer than 10 kb were aligned with Gyeonggijaerae 5 reads using the Burrows-Wheeler Aligner. SAMTools was used for retrieving single nucleotide polymorphisms (SNPs) between Sunhwanokdu and Gyeonggijaerae 5, defining the lowest and highest depths as 5 and 100, respectively, and the sequence quality as 100. Of the 305,504 single-base changes identified, 40,503 SNPs were considered heterozygous in Gyeonggijaerae 5. Among the remaining 265,001 SNPs, 65.9 % (174,579 cases) were transitions and 34.1 % (90,422 cases) were transversions. For SNP validation, a total of 42 SNPs were chosen among Sunhwanokdu contigs longer than 10 kb and sharing at least 80 % sequence identity with common bean expressed sequence tags as determined with est2genome. Using seven mungbean cultivars from various origins in addition to Sunhwanokdu and Gyeonggijaerae 5, most of the SNPs identified by bioinformatics tools were confirmed by Sanger sequencing. These genome-wide SNP markers could enrich the current molecular resources and might be of value for the construction of a mungbean genetic map and the investigation of genetic diversity.

Dynamic genetic features of chromosomes revealed by comparison of soybean genetic and sequence-based physical maps.

Despite the intensive soybean [Glycine max (L.) Merrill] genome studies, the high chromosome number (20) of the soybean plant relative to many other major crops has hindered the development of a high-resolution genomewide genetic map derived from a single population. Here, we report such a map, which was constructed in an F15 population derived from a cross between G. max and G. soja lines using indel polymorphisms detected via a G. soja genome resequencing. By targeting novel indel markers to marker-poor regions, all marker intervals were reduced to under 6 cM on a genome scale. Comparison of the Williams 82 soybean reference genome sequence and our genetic map indicated that marker orders of 26 regions were discrepant with each other. In addition, our comparison showed seven misplaced and two absent markers in the current Williams 82 assembly and six markers placed on the scaffolds that were not incorporated into the pseudomolecules. Then, we showed that, by determining the missing sequences located at the presumed beginning points of the five major discordant segments, these observed discordant regions are mostly errors in the Williams 82 assembly. Distributions of the recombination rates along the chromosomes were similar to those of other organisms. Genotyping of indel markers and genome resequencing of the two parental lines suggested that some marker-poor chromosomal regions may represent introgression regions, which appear to be prevalent in soybean. Given the even and dense distribution of markers, our genetic map can serve as a bridge between genomics research and breeding programs.

Whole-genome sequencing and intensive analysis of the undomesticated soybean (Glycine soja Sieb. and Zucc.) genome.

The genome of soybean (Glycine max), a commercially important crop, has recently been sequenced and is one of six crop species to have been sequenced. Here we report the genome sequence of G. soja, the undomesticated ancestor of G. max (in particular, G. soja var. IT182932). The 48.8-Gb Illumina Genome Analyzer (Illumina-GA) short DNA reads were aligned to the G. max reference genome and a consensus was determined for G. soja. This consensus sequence spanned 915.4 Mb, representing a coverage of 97.65% of the G. max published genome sequence and an average mapping depth of 43-fold. The nucleotide sequence of the G. soja genome, which contains 2.5 Mb of substituted bases and 406 kb of small insertions/deletions relative to G. max, is ∼0.31% different from that of G. max. In addition to the mapped 915.4-Mb consensus sequence, 32.4 Mb of large deletions and 8.3 Mb of novel sequence contigs in the G. soja genome were also detected. Nucleotide variants of G. soja versus G. max confirmed by Roche Genome Sequencer FLX sequencing showed a 99.99% concordance in single-nucleotide polymorphism and a 98.82% agreement in insertion/deletion calls on Illumina-GA reads. Data presented in this study suggest that the G. soja/G. max complex may be at least 0.27 million y old, appearing before the relatively recent event of domestication (6,000∼9,000 y ago). This suggests that soybean domestication is complicated and that more in-depth study of population genetics is needed. In any case, genome comparison of domesticated and undomesticated forms of soybean can facilitate its improvement.

Genome-wide analysis of KIX gene family for organ size regulation in soybean (Glycine max L.).

The KIX domain, conserved among various nuclear and co-activator factors, acts as a binding site that interacts with other transcriptional activators and co-activators, playing a crucial role in gene expression regulation. In plants, the KIX domain is involved in plant hormone signaling, stress response regulation, cell cycle control, and differentiation, indicating its potential relevance to crop productivity. This study aims to identify and characterize KIX domains within the soybean (Glycine max L.) genome to predict their potential role in improving crop productivity. The conservation and evolutionary history of the KIX domains were explored in 59 plant species, confirming the presence of the KIX domains in diverse plants. Specifically, 13 KIX domains were identified within the soybean genome and classified into four main groups, namely GmKIX8/9, GmMED15, GmHAC, and GmRECQL, through sequence alignment, structural analysis, and phylogenetic tree construction. Association analysis was performed between KIX domain haplotypes and soybean seed-related agronomic traits using re-sequencing data from a core collection of 422 accessions. The results revealed correlations between SNP variations observed in GmKIX8-3 and GmMED15-4 and soybean seed phenotypic traits. Additionally, transcriptome analysis confirmed significant expression of the KIX domains during the early stages of soybean seed development. This study provides the first characterization of the structural, expression, genomic haplotype, and molecular features of the KIX domain in soybean, offering a foundation for functional analysis of the KIX domain in soybean and other plants.

Genetic Diversity of Korean Wild Soybean Core Collections and Genome-Wide Association Study for Days to Flowering.

The utilization of wild soybean germplasms in breeding programs increases genetic diversity, and they contain the rare alleles of traits of interest. Understanding the genetic diversity of wild germplasms is essential for determining effective strategies that can improve the economic traits of soybeans. Undesirable traits make it challenging to cultivate wild soybeans. This study aimed to construct a core subset of 1467 wild soybean accessions of the total population and analyze their genetic diversity to understand their genetic variations. Genome-wild association studies were conducted to detect the genetic loci underlying the time to flowering for a core subset collection, and they revealed the allelic variation in E genes for predicting maturity using the available resequencing data of wild soybean. Based on principal component and cluster analyses, 408 wild soybean accessions in the core collection covered the total population and were explained by 3 clusters representing the collection regions, namely, Korea, China, and Japan. Most of the wild soybean collections in this study had the E1e2E3 genotype according to association mapping and a resequencing analysis. Korean wild soybean core collections can provide helpful genetic resources to identify new flowering and maturity genes near the E gene loci and genetic materials for developing new cultivars, facilitating the introgression of genes of interest from wild soybean.

A Genome-Wide Association Study of Protein, Oil, and Amino Acid Content in Wild Soybean (Glycine soja).

Soybean (Glycine max L.) is a globally important source of plant proteins, oils, and amino acids for both humans and livestock. Wild soybean (Glycine soja Sieb. and Zucc.), the ancestor of cultivated soybean, could be a useful genetic source for increasing these components in soybean crops. In this study, 96,432 single-nucleotide polymorphisms (SNPs) across 203 wild soybean accessions from the 180K Axiom® Soya SNP array were investigated using an association analysis. Protein and oil content exhibited a highly significant negative correlation, while the 17 amino acids exhibited a highly significant positive correlation with each other. A genome-wide association study (GWAS) was conducted on the protein, oil, and amino acid content using the 203 wild soybean accessions. A total of 44 significant SNPs were associated with protein, oil, and amino acid content. Glyma.11g015500 and Glyma.20g050300, which contained SNPs detected from the GWAS, were selected as novel candidate genes for the protein and oil content, respectively. In addition, Glyma.01g053200 and Glyma.03g239700 were selected as novel candidate genes for nine of the amino acids (Ala, Asp, Glu, Gly, Leu, Lys, Pro, Ser, and Thr). The identification of the SNP markers related to protein, oil, and amino acid content reported in the present study is expected to help improve the quality of selective breeding programs for soybeans.

Fine genetic mapping of the genomic region controlling leaflet shape and number of seeds per pod in the soybean.

Narrow leaflet cultivars tend to have more seeds per pod than broad leaflet cultivars in soybean [Glycine max (L.) Merr.], which suggests that the leaflet-shape trait locus is tightly linked to or cosegregates with the trait locus controlling the number of seeds per pod (NSPP). Here, we attempted to further elucidate the relationship between leaflet shape and NSPP. A BC(3)F(2) population from a cross between the 'Sowon' (narrow leaflets and high NSPP) and 'V94-5152' (broad leaflets and low NSPP) variants was used. The results of the molecular genetic analyses indicated that, although the NSPP characteristic, in particular, the occurrence of 4-seeded pods, is governed by additional modifying genes that are likely present in Sowon, the two traits cosegregate in the BC(3)F(2) population. The mapping results generated using public markers demonstrated that the narrow leaflet-determining gene in Sowon is an allele of the previously highly studied ln gene on chromosome 20. A high-resolution map delimited the genomic region controlling both the leaflet shape and NSPP traits to a sequence length of 66 kb, corresponding to 0.7 cM. Among the three genes annotated in this 66 kb region, Glyma20g25000.1 appeared to be a good candidate for the Ln-encoding gene, owing to its 47.8% homology with the protein encoding for the JAGGED gene that regulates lateral organ development in Arabidopsis. Taken together, our results suggested that phenotypic variations for narrow leaflet and NSPP are predominantly from the pleiotropic effects of the ln gene. Thus, our results should provide a molecular framework for soybean breeding programs with the objective of improving soybean yield.