A retrotransposon insertion in the Mao1 promoter results in erect pubescence and higher yield in soybean.

Adaptive changes in crops contribute to the diversity of agronomic traits, which directly or indirectly affect yield. The change of pubescence form from appressed to erect is a notable feature during soybean domestication. However, the biological significance and regulatory mechanism underlying this transformation remain largely unknown. Here, we identified a major-effect locus, PUBESCENCE FORM 1 (PF1), the upstream region of Mao1, that regulates pubescence form in soybean. The insertion of a Ty3/Gypsy retrotransposon in PF1 can recruit the transcription factor GAGA-binding protein to a GA-rich region, which up-regulates Mao1 expression, underpinning soybean pubescence evolution. Interestingly, the proportion of improved cultivars with erect pubescence increases gradually with increasing latitude, and erect-pubescence cultivars have a higher yield possibly through a higher photosynthetic rate and photosynthetic stability. These findings open an avenue for molecular breeding through either natural introgression or genome editing toward yield improvement and productivity.

Tandemly duplicated MYB genes are functionally diverged in the regulation of anthocyanin biosynthesis in soybean.

Gene duplications have long been recognized as a driving force in the evolution of genes, giving rise to novel functions. The soybean (Glycine max) genome is characterized by a large number of duplicated genes. However, the extent and mechanisms of functional divergence among these duplicated genes in soybean remain poorly understood. In this study, we revealed that 4 MYB genes (GmMYBA5, GmMYBA2, GmMYBA1, and Glyma.09g235000)-presumably generated by tandem duplication specifically in the Phaseoleae lineage-exhibited a stronger purifying selection in soybean compared to common bean (Phaseolus vulgaris). To gain insights into the diverse functions of these tandemly duplicated MYB genes in anthocyanin biosynthesis, we examined the expression, transcriptional activity, induced metabolites, and evolutionary history of these 4 MYB genes. Our data revealed that Glyma.09g235000 is a pseudogene, while the remaining 3 MYB genes exhibit strong transcriptional activation activity, promoting anthocyanin biosynthesis in different soybean tissues. GmMYBA5, GmMYBA2, and GmMYBA1 induced anthocyanin accumulation by upregulating the expression of anthocyanin pathway-related genes. Notably, GmMYBA5 showed a lower capacity for gene induction compared to GmMYBA2 and GmMYBA1. Metabolomics analysis further demonstrated that GmMYBA5 induced distinct anthocyanin accumulation in Nicotiana benthamiana leaves and soybean hairy roots compared to GmMYBA2 and GmMYBA1, suggesting their functional divergence leading to the accumulation of different metabolites accumulation following gene duplication. Together, our data provide evidence of functional divergence within the MYB gene cluster following tandem duplication, which sheds light on the potential evolutionary directions of gene duplications during legume evolution.

Identification of an important QTL for seed oil content in soybean.

Seed oil content is one of the most important quantitative traits in soybean (Glycine max) breeding. Here, we constructed a high-density single nucleotide polymorphism linkage map using two genetically similar parents, Heinong 84 and Kenfeng 17, that differ dramatically in their seed oil contents, and performed quantitative trait loci (QTL) mapping of seed oil content in a recombinant inbred line (RIL) population derived from their cross. We detected five QTL related to seed oil content distributed on five chromosomes. The QTL for seed oil content explained over 10% of the phenotypic variation over two years. This QTL was mapped to an interval containing 20 candidate genes, including a previously reported gene, soybean RING Finger 1a (RNF1a) encoding an E3 ubiquitin ligase. Notably, two short sequences were inserted in the GmRNF1a coding region of KF 17 compared to that of HN 84, resulting in a longer protein variant in KF 17. Our results thus provide information for uncovering the genetic mechanisms determining seed oil content in soybean, as well as identifying an additional QTL and highlighting GmRNF1a as candidate gene for modulating seed oil content in soybean. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-023-01384-2.

Identification of ST1 reveals a selection involving hitchhiking of seed morphology and oil content during soybean domestication.

Seed morphology and quality of cultivated soybean (Glycine max) have changed dramatically during domestication from their wild relatives, but their relationship to selection is poorly understood. Here, we describe a semi-dominant locus, ST1 (Seed Thickness 1), affecting seed thickness and encoding a UDP-D-glucuronate 4-epimerase, which catalyses UDP-galacturonic acid production and promotes pectin biosynthesis. Interestingly, this morphological change concurrently boosted seed oil content, which, along with up-regulation of glycolysis biosynthesis modulated by ST1, enabled soybean to become a staple oil crop. Strikingly, ST1 and an inversion controlling seed coat colour formed part of a single selective sweep. Structural variation analysis of the region surrounding ST1 shows that the critical mutation in ST1 existed in earlier wild relatives of soybean and the region containing ST1 subsequently underwent an inversion, which was followed by successive selection for both traits through hitchhiking during selection for seed coat colour. Together, these results provide direct evidence that simultaneously variation for seed morphology and quality occurred earlier than variation for seed coat colour during soybean domestication. The identification of ST1 thus sheds light on a crucial phase of human empirical selection in soybeans and provides evidence that our ancestors improved soybean based on taste.

Pathogenicity and genome-wide sequence analysis reveals relationships between soybean mosaic virus strains.

Soybean mosaic virus (SMV) is the most prevalent viral pathogen in soybean. In China, the SMV strains SC and N are used simultaneously in SMV resistance assessments of soybean cultivars, but the pathogenic relationship between them is unclear. In this study, SMV strains N1 and N3 were found to be the most closely related to SC18. Moreover, N3 was found to be more virulent than N1. A global pathotype classification revealed the highest level of genetic diversity in China. The N3 type was the most frequent and widespread worldwide, implying that SMV possibly originated in China and spread across continents through the dissemination of infected soybean. It also suggests that the enhanced virulence of N3 facilitated its spread and adaptability in diverse geographical and ecological regions worldwide. Phylogenetic analysis revealed prominent geographical associations among SMV strains/isolates, and genomic nucleotide diversity analysis and neutrality tests demonstrated that the whole SMV genome is under negative selection, with the P1 gene being under the greatest selection pressure. The results of this study will facilitate the nationwide use of SMV-resistant soybean germplasm and could provide useful insights into the molecular variability, geographical distribution, phylogenetic relationships, and evolutionary history of SMV around the world.

MALE STERILITY 3 encodes a plant homeodomain-finger protein for male fertility in soybean.

Male-sterile plants are used in hybrid breeding to improve yield in soybean (Glycine max (L.) Merr.). Developing the capability to alter fertility under different environmental conditions could broaden germplasm resources and simplify hybrid production. However, molecular mechanisms potentially underlying such a system in soybean were unclear. Here, using positional cloning, we identified a gene, MALE STERILITY 3 (MS3), which encodes a nuclear-localized protein containing a plant homeodomain (PHD)-finger domain. A spontaneous mutation in ms3 causing premature termination of MS3 translation and partial loss of the PHD-finger. Transgenetic analysis indicated that MS3 knockout resulted in nonfunctional pollen and no self-pollinated pods, and RNA-seq analysis revealed that MS3 affects the expression of genes associated with carbohydrate metabolism. Strikingly, the fertility of mutant ms3 can restore under long-d conditions. The mutant could thus be used to create a new, more stable photoperiod-sensitive genic male sterility line for two-line hybrid seed production, with significant impact on hybrid breeding and production.

A transposon-mediated reciprocal translocation promotes environmental adaptation but compromises domesticability of wild soybeans.

Large structural variations frequently occur in higher plants; however, the impact of such variations on plant diversification, adaptation and domestication remains elusive. Here, we mapped and characterised a reciprocal chromosomal translocation in soybeans and assessed its effects on diversification and adaptation of wild (Glycine soja) and semiwild (Glycine gracilis) soybeans, and domestication of cultivated soybean (Glycine max), by tracing the distribution of the translocation in the USDA Soybean Germplasm Collection and population genetics analysis. We demonstrate that the translocation occurred through CACTA transposon-mediated chromosomal breakage in wild soybean c. 0.34 Ma and is responsible for semisterility in translocation heterozygotes and reduces their reproductive fitness. The translocation has differentiated Continental (i.e. China and Russia) populations from Maritime (i.e. Korea and Japan) populations of G. soja and predominately adapted to cold and dry climates. Further analysis revealed that the divergence of G. max from G. soja predates the translocation event and that G. gracilis is an evolutionary intermediate between G. soja and G. max. Our results highlight the effects of a chromosome rearrangement on the processes leading to plant divergence and adaptation, and provides evidence that suggests G. gracilis, rather than G. soja, as the ancestor of cultivated soybean.

Pedigree-based genetic dissection of quantitative loci for seed quality and yield characters in improved soybean.

As soybean plays an indispensable role in the supply of vegetable oil and protein, balancing the relationship between seed quality and yield traits according to human demand has become an important breeding goal for soybean improvement. Here, 256 intraspecific recombinant inbred lines (RILs), derived from a cross between Qi Huang No.34 (QH34) and Ji Dou No.17 (JD17), were used for quantitative trait loci (QTLs) mapping with remarkable four chemical and physical properties with a purpose for exploring the distribution of excellent alleles in germplasm resources in China. A total of 25 QTLs were detected, of which 10 QTLs inherited the alleles from the parent QH34. Pedigree research on favorable alleles on these QTLs showed the process of excellent alleles pyramided into QH34. Meta-analysis of the 25 QTLs by comparing with existed QTLs in previous study identified 17 novel QTLs. QTLs with pleiotropic effects have been detected. Furthermore, three representative elite recombinant inbred lines in different locations that have great potential in soybean breeding were selected, and finally, four seed weight-related candidate genes were identified. The discovery of these QTLs provides a new guidance for combining the diversity and rarity of germplasm resources, which can effectively increase population genetic diversity and broaden genetic basis of varieties. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-021-01211-6.

Genetic dissection of seed appearance quality using recombinant inbred lines in soybean.

Soybean seed appearance quality greatly affects the marketability. The objective of this study was to identify the quantitative trait loci (QTLs) that control the appearance quality of soybean seeds. A total of 256 recombinant inbred lines from Qi Huang No.34 × Ji Dou No.17 were utilized for QTL mapping. We innovatively applied a machine vision system to quantify the seed appearance of each line. As a result of QTL mapping, a total of 145 QTLs for the machine vision parameters were detected across three environments. We integrated QTLs mapped overlapped and obtained 16 QTL hotspots in total. Of these hotspots, hotspot-4-1 was suggested to be a major locus controlling seed size, and hotspot-15 was identified to affect the seed color and texture. The mapping for principal components of the seed appearance also supported it. This study comprehensively dug up the QTLs for seed appearance quality of soybean cultivars while providing an efficient method for phenotyping of seed appearance. These results would contribute to dissecting the genetic bases of seed appearance quality for the improvement of soybean. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-021-01262-9.

Genetic dissection of domestication-related traits in soybean through genotyping-by-sequencing of two interspecific mapping populations.

KEY MESSAGE: A total of 132 domestication-related QTLs, of which 41 were novel, were identified through genotyping-by-sequencing of two Glycine max × Glycine soja populations. Soybean [Glycine max (L.) Merr.] was domesticated in East Asia from the wild progenitor Glycine soja. The domestication process led to many distinct morphological changes that adapt it to cultivation. These include larger seeds, erect growth, larger stem diameter, reduced pod shattering, and altered growth habit. The objective of this study was to identify QTLs controlling key domestication-related traits (DRTs) using interspecific mapping populations. A total of 151 RILs from Williams 82 × PI 468916 and 510 RILs from Williams 82 × PI 479752 were utilized for QTL mapping. These lines were genotyped using a genotyping-by-sequencing protocol which resulted in approximately 5000 polymorphic SNP markers. The number of QTLs detected for each of the eleven DRTs ranged between 0-4 QTLs in the smaller Williams 82 × PI 468916 population and 3-16 QTLs in the larger Williams 82 × PI 479752 population. A total of 132 QTLs were detected, of which 51 are associated with selective sweeps previously related to soybean domestication. These QTLs were detected across all 20 chromosomes within 42 genomic regions. This study identifies 41 novel QTLs not detected in previous studies using smaller populations while also confirming the quantitative nature for several of the important DRTs in soybeans. These results would enable more effective use of the wild germplasm for soybean improvement.

Plasticity and innovation of regulatory mechanisms underlying seed oil content mediated by duplicated genes in the palaeopolyploid soybean.

Many plants have undergone whole genome duplication (WGD). However, how regulatory networks underlying a particular trait are reshaped in polyploids has not been experimentally investigated. Here we show that the regulatory pathways modulating seed oil content, which involve WRINKLED1 (WRI1), LEAFY COTYLEDON1 (LEC1), and LEC2 in Arabidopsis, have been modified in the palaeopolyploid soybean. Such modifications include functional reduction of GmWRI1b of the GmWRI1a/GmWRI1b homoeologous pair relevant to WRI1, complementary non-allelic dosage effects of the GmLEC1a/GmLEC1b homoeologous pair relevant to LEC1, pseudogenization of the singleton GmLEC2 relevant to LEC2, and the rise of the LEC2-like function of GmABI3b, contrasting to its homoeolog GmABI3a, which maintains the ABSCISIC ACID INSENSITIVE 3 (ABI3)-like function in modulating seed maturation and dormancy. The function of GmABI3b in modulating seed oil biosynthesis was fulfilled by direct binding to a RY (CATGCA) cis-regulatory element in the GmWRI1a promoter, which was absent in the GmWRI1b promoter, resulting in reduction of the GmWRI1b expression. Nevertheless, the three regulators each exhibited similar intensities of purifying selection to their respective duplicates since these pairs were formed by a WGD event that is proposed to have occurred approximately 13 million years ago (mya), suggesting that the differentiation in spatiotemporal expression between the duplicated genes is more likely to be the outcome of neutral variation in regulatory sequences. This study thus exemplifies the plasticity, dynamics, and novelty of regulatory networks mediated by WGD.

Dt2 is a gain-of-function MADS-domain factor gene that specifies semideterminacy in soybean.

Similar to Arabidopsis thaliana, the wild soybeans (Glycine soja) and many cultivars exhibit indeterminate stem growth specified by the shoot identity gene Dt1, the functional counterpart of Arabidopsis TERMINAL FLOWER1 (TFL1). Mutations in TFL1 and Dt1 both result in the shoot apical meristem (SAM) switching from vegetative to reproductive state to initiate terminal flowering and thus produce determinate stems. A second soybean gene (Dt2) regulating stem growth was identified, which, in the presence of Dt1, produces semideterminate plants with terminal racemes similar to those observed in determinate plants. Here, we report positional cloning and characterization of Dt2, a dominant MADS domain factor gene classified into the APETALA1/SQUAMOSA (AP1/SQUA) subfamily that includes floral meristem (FM) identity genes AP1, FUL, and CAL in Arabidopsis. Unlike AP1, whose expression is limited to FMs in which the expression of TFL1 is repressed, Dt2 appears to repress the expression of Dt1 in the SAMs to promote early conversion of the SAMs into reproductive inflorescences. Given that Dt2 is not the gene most closely related to AP1 and that semideterminacy is rarely seen in wild soybeans, Dt2 appears to be a recent gain-of-function mutation, which has modified the genetic pathways determining the stem growth habit in soybean.

Fine mapping and candidate gene analysis of two loci conferring resistance to Phytophthora sojae in soybean.

KEY MESSAGE: RpsUN1 and RpsUN2 were fine mapped to two genomic regions harboring disease resistance-like genes. The haplotypes and instability of the regions and candidate genes for the two resistance loci were characterized. Phytophthora root and stem rot caused by Phytophthora sojae, is one of the most destructive diseases of soybean. Deploying soybean cultivars carrying race-specific resistance conferred by Rps genes is the most practical approach to managing this disease. Previously, two Rps genes, RpsUN1 and RpsUN2 were identified in a landrace PI 567139B and mapped to a 6.5 cM region on chromosome 3 and a 3.0 cM region on chromosome 16, corresponding to 1387 and 423 kb of the soybean reference genome sequences. By analyzing recombinants defined by genotypic and phenotypic screening of the 826 F2:3 families derived from two reciprocal crosses between the two parental lines, RpsUN1 and RpsUN2, were further narrowed to a 151 kb region that harbors five genes including three disease resistance (R)-like genes, and a 36 kb region that contains four genes including five R-like genes, respectively, according to the reference genome. Expressional changes of these nine genes before and after inoculation with the pathogen, as revealed by RNA-seq, suggest that Glyma.03g034600 in the RpsUN1 region and Glyma.16g215200 and Glyma.16g214900 in the RpsUN2 region of PI 567139B may be associated with the resistance to P. sojae. It is also suggested that unequal recombination between/among R-like genes may have occurred, resulting in the formation of two recombinants with inconsistent genotypic and phenotypic observations. The haplotype variation of genomic regions where RpsUN1 and RpsUN2 reside in the entire soybean germplasm deposited in the US soybean germplasm collection suggests that RpsUN1 and RpsUN2 are most likely novel genes.

Molecular response to the pathogen Phytophthora sojae among ten soybean near isogenic lines revealed by comparative transcriptomics.

BACKGROUND: Phytophthora root and stem rot (PRR) of soybean, caused by Phytophthora sojae, is controlled by Rps genes. However, little is known regarding the Rps-induced molecular responses to P. sojae and how they actually overlap. We thus sequenced, analyzed, and compared the transcriptomes of 10 near isogenic lines (NILs), each with a unique Rps gene/allele, and the susceptible parent Williams, pre- and post-inoculation with the pathogen. RESULTS: A total of 4,330 differentially expressed genes (DEGs) were identified in Williams versus 2,014 to 5,499 DEGs in individual NILs upon inoculation with the pathogen. Comparisons of the DEGs between the NILs and Williams identified incompatible interaction genes (IIGs) and compatible interaction genes (CIGs). Hierarchical cluster and heatmap analyses consistently grouped the NILs into three clusters: Cluster I (Rps1-a), Cluster II (Rps1-b, 1-c and 1-k) and Cluster III (Rps3-a, 3-b, 3-c, 4, 5, and 6), suggesting an overlap in Rps-induced defense signaling among certain NILs. Gene ontology (GO) analysis revealed associations between members of the WRKY family and incompatible reactions and between a number of phytohormone signaling pathways and incompatible/compatible interactions. These associations appear to be distinguished according to the NIL clusters. CONCLUSIONS: This study characterized genes and multiple branches of putative regulatory networks associated with resistance to P. sojae in ten soybean NILs, and depicted functional "fingerprints" of individual Rps-mediated resistance responses through comparative transcriptomic analysis. Of particular interest are dramatic variations of detected DEGs, putatively involved in ethylene (ET)-, jasmonic acid (JA)-, (reactive oxygen species) ROS-, and (MAP-kinase) MAPK- signaling, among these soybean NILs, implicating their important roles of these signaling in differentiating molecular defense responses. We hypothesize that different timing and robustness in defense signaling to the same pathogen may be largely responsible for such variations.

GS6, a member of the GRAS gene family, negatively regulates grain size in rice.

Grain size is an important yield-related trait in rice. Intensive artificial selection for grain size during domestication is evidenced by the larger grains of most of today's cultivars compared with their wild relatives. However, the molecular genetic control of rice grain size is still not well characterized. Here, we report the identification and cloning of Grain Size 6 (GS6), which plays an important role in reducing grain size in rice. A premature stop at the +348 position in the coding sequence (CDS) of GS6 increased grain width and weight significantly. Alignment of the CDS regions of GS6 in 90 rice materials revealed three GS6 alleles. Most japonica varieties (95%) harbor the Type I haplotype, and 62.9% of indica varieties harbor the Type II haplotype. Association analysis revealed that the Type I haplotype tends to increase the width and weight of grains more than either of the Type II or Type III haplotypes. Further investigation of genetic diversity and the evolutionary mechanisms of GS6 showed that the GS6 gene was strongly selected in japonica cultivars. In addition, a "ggc" repeat region identified in the region that encodes the GRAS domain of GS6 played an important historic role in the domestication of grain size in rice. Knowledge of the function of GS6 might aid efforts to elucidate the molecular mechanisms that control grain development and evolution in rice plants, and could facilitate the genetic improvement of rice yield.

A reference-grade genome assembly for Astragalus mongholicus and insights into the biosynthesis and high accumulation of triterpenoids and flavonoids in its roots.

Astragalus membranaceus var. mongholicus (AMM), a member of the Leguminosae, is one of the most important medicinal plants worldwide. The dried roots of AMM have a wide range of pharmacological effects and are a traditional Chinese medicine. Here, we report the first chromosome-level reference genome of AMM, comprising nine pseudochromosomes with a total size of 1.47 Gb and 27 868 protein-encoding genes. Comparative genomic analysis reveals that AMM has not experienced an independent whole-genome duplication (WGD) event after the WGD event shared by the Papilionoideae species. Analysis of long terminal repeat retrotransposons suggests a recent burst of these elements at approximately 0.13 million years ago, which may explain the large size of the AMM genome. Multiple gene families involved in the biosynthesis of triterpenoids and flavonoids were expanded, and our data indicate that tandem duplication has been the main driver for expansion of these families. Among the expanded families, the phenylalanine ammonia-lyase gene family was primarily expressed in the roots of AMM, suggesting their roles in the biosynthesis of phenylpropanoid compounds. The functional versatility of 2,3-oxidosqualene cyclase genes in cluster III may play a critical role in the diversification of triterpenoids in AMM. Our findings provide novel insights into triterpenoid and flavonoid biosynthesis and can facilitate future research on the genetics and medical applications of AMM.

Enhanced Microcystis Aeruginosa removal and novel flocculation mechanisms using a novel continuous co-coagulation flotation (CCF).

Co-coagulation flotation (CCF) is a novel flotation technology that renders more efficient algal removal compared to traditional mechanical coagulation flotation (MCF) due to a short residence time (< 30 s) and fast rising behavior of algal flocs (> 250 m·h-1). This study compared the algal removal performance using continuous CCF and MCF using water samples taken from Lake Dianchi with severe Microcystis aeruginosa blooms. Removal efficiency, dosage of coagulant/flocculant, rising velocity and structural characteristics of the resulting flocs in the two processes were systematically compared. The results show that CCF could save >50 % polyaluminum chloride (PAC) and polyacrylamide (PAM) compared with MCF when the removal efficiency was both over 95 %. The average rising velocity of flocs in CCF could reach 254.3 m·h-1, much higher than that in MCF (154.5 m·h-1). In the respective optimal coagulation conditions, the flocs formed in CCF (G = 164.8 s-1) were larger (1843 ± 128 µm) and more spherical with a higher fractal dimension (Df = 1.85 ± 0.01) than those generated in MCF (G = 34.1 s-1). The Stokes's Law was found to correctly predict the rising velocity of spherical flocs with large fractal dimensions (Df > 1.7). In contrast, the Haarhoff and Edzwald's extended equation was more suitable for calculating the rising velocity of irregular flocs with small fractal dimension. This study provides new insights into the mechanisms of the enhanced algal removal by CCF and lays foundation for developing cost-efficient algal mitigation processes.

TH1, a DUF640 domain-like gene controls lemma and palea development in rice.

The developmental regulation of grasses lemma and palea and their relationship to the floral organs in dicots had been variously explicated and extensively debated. Here, we characterized a triangular hull mutant th1-1 from EMS-mutagenized Oryza sativa ssp. indica cv. 93-11. The th1-1 mutant exhibited obviously triangular hull with tortuous and slender lemma/palea. Using a map-based cloning strategy, the TH1 gene was narrowed down to a 60-kb region on the long arm of chromosome 2. Sequence verification revealed that the th1-1 mutant harbored 1-bp deletion in exon 2 of LOC_Os02g56610 which resulted in a frame-shift mutation. The RNA-interference transgenic plants of LOC_Os02g56610 displayed a similar phenotype to the th1 mutant. Consequently, LOC_Os02g56610 was identified as the TH1 gene which encoded 248 amino acids and contained a DUF640 domain. RT-PCR analysis and GUS staining showed that the transcripts of TH1 mainly accumulated in young inflorescence, lemma and palea of spikelet. These results suggested that TH1 was an important gene controlling the lemma and palea development in rice.

Identification of Potential Gene Regulatory Pathways Affecting the Ratio of Four-Seed Pod in Soybean.

The number of four-seed pods is one of the most important agronomic traits affected by gene and environment that can potentially improve soybean (Glycine max) yield. However, the gene regulatory network that affects the ratio of four-seed pod (the ratio of the number of four-seed pods to the total number of pods in each individual plant) is yet unclear. Here, we performed bulked segregant RNA sequencing (BSR-seq) on a series of recombinant inbred lines (RILs) derived from hybrid progenies between Heinong 48 (HN48), a cultivar with a high ratio of four-seed pod, and Henong 64 (HN64), a cultivar with a low ratio of four-seed pod. Two tissues, flower bud and young pod, at two different growth stages, R1 and R3, were analyzed under the ratios of four-seed pod at less than 10% and greater than 30%, respectively. To identify the potential gene regulation pathways associated with the ratio of soybean four-seed pod, we performed differentially expressed analysis on the four bulked groups. A differentially expressed gene (DEG) encoding a photosystem II 5-kDa protein had the function of participating in the energy conversion of photosynthesis. In addition, 79 common DEGs were identified at different developmental stages and under different ratios of four-seed pod. Among them, four genes encoding calcium-binding proteins and a WRKY transcription factor were enriched in the plant-pathogen interaction pathway, and they showed a high level of expression in roots. Moreover, 10 DEGs were identified in the reported quantitative trait locus (QTL) interval of four-seed pod, and two of them were significantly enriched in the pentose and glucuronate interconversion pathway. These findings provide basic insights into the understanding of the underlying gene regulatory network affected by specific environment and lay the foundation for identifying the targets that affect the ratio of four-seed pod in soybean.

Elevation of soybean seed oil content through selection for seed coat shininess.

Many leguminous species have adapted their seed coat with a layer of powdery bloom that contains hazardous allergens and makes the seeds less visible, offering duel protection against potential predators 1 . Nevertheless, a shiny seed surface without bloom is desirable for human consumption and health, and is targeted for selection under domestication. Here we show that seed coat bloom in wild soybeans is mainly controlled by Bloom1 (B1), which encodes a transmembrane transporter-like protein for biosynthesis of the bloom in pod endocarp. The transition from the 'bloom' to 'no-bloom' phenotypes is associated with artificial selection of a nucleotide mutation that naturally occurred in the coding region of B1 during soybean domestication. Interestingly, this mutation not only 'shined' the seed surface, but also elevated seed oil content in domesticated soybeans. Such an elevation of oil content in seeds appears to be achieved through b1-modulated upregulation of oil biosynthesis in pods. This study shows pleiotropy as a mechanism underlying the domestication syndrome 2 , and may pave new strategies for development of soybean varieties with increased seed oil content and reduced seed dust.