Search

Soybean (Glycine max) produces seeds that are rich in unsaturated fatty acids and is an important oilseed crop worldwide. Seed oil content and composition largely determine the economic value of soybean. Due to natural genetic variation, seed oil content varies substantially across soybean cultivars. Although much progress has been made in elucidating the genetic trajectory underlying fatty acid metabolism and oil biosynthesis in plants, the causal genes for many quantitative trait loci (QTLs) regulating seed oil content in soybean remain to be revealed. In this study, we identified GmFATA1B as the gene underlying a QTL that regulates seed oil content and composition, as well as seed size in soybean. Nine extra amino acids in the conserved region of GmFATA1B impair its function as a fatty acyl-acyl carrier protein thioesterase, thereby affecting seed oil content and composition. Heterogeneously overexpressing the functional GmFATA1B allele in Arabidopsis thaliana increased both the total oil content and the oleic acid and linoleic acid contents of seeds. Our findings uncover a previously unknown locus underlying variation in seed oil content in soybean and lay the foundation for improving seed oil content and composition in soybean., (© 2023 Institute of Botany, Chinese Academy of Sciences.)

Plant resistance genes could affect rhizosphere microbiota, which in turn enhanced plant resistance to stresses. Our previous study found that overexpression of the GsMYB10 gene led to enhanced tolerance of soybean plants to aluminum (Al) toxicity. However, whether GsMYB10 gene could regulate rhizosphere microbiota to mitigate Al toxicity remains unclear. Here, we analyzed the rhizosphere microbiomes of HC6 soybean (WT) and transgenic soybean (trans-GsMYB10) at three Al concentrations, and constructed three different synthetic microbial communities (SynComs), including bacterial, fungal and cross-kingdom (bacteria and fungi) SynComs to verify their role in improving Al tolerance of soybean. Trans-GsMYB10 shaped the rhizosphere microbial communities and harbored some beneficial microbes, such as Bacillus, Aspergillus and Talaromyces under Al toxicity. Fungal and cross-kingdom SynComs showed a more effective role than the bacterial one in resistance to Al stress, and these SynComs helped soybean resist Al toxicity via affecting some functional genes that involved cell wall biosynthesis and organic acid transport etc. Overall, this study reveals the mechanism of soybean functional genes regulating the synergistic resistance of rhizosphere microbiota and plants to Al toxicity, and also highlights the possibility of focusing on the rhizobial microbial community as a potential molecular breeding target to produce crops., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier B.V. All rights reserved.)

Soybean is a major crop that produces valuable seed oil and protein for global consumption. Seed oil and protein are regulated by complex quantitative trait loci (QTLs) and have undergone intensive selections during the domestication of soybean. It is essential to identify the major genetic components and understand their mechanism behind seed oil and protein in soybean. We report that MOTHER-OF-FT-AND-TFL1 (GmMFT) is the gene of a classical QTL that has been reported to regulate seed oil and protein content in many studies. Mutation of MFT decreased seeds oil content and weight in both Arabidopsis and soybean, whereas increased expression of GmMFT enhanced seeds oil content and weight. Haplotype analysis showed that GmMFT has undergone selection, which resulted in the extended haplotype homozygosity in the cultivated soybean and the enriching of the oil-favorable allele in modern soybean cultivars. This work unraveled the GmMFT-mediated mechanism regulating seed oil and protein content and seed weight, and revealed a previously unknown function of MFT that provides new insights into targeted soybean improvement and breeding., (© 2023 The Authors. New Phytologist © 2023 New Phytologist Foundation.)

Salt stress has a detrimental impact on crop yield, quality, and profitability. The tau-like glutathione transferases (GSTs) represent a significant group of enzymes that play a crucial role in plant stress responses, including salt stress. In this study, we identified a tau-like glutathione transferase family gene from soybean named GmGSTU23 . Expression pattern analysis revealed that GmGSTU23 was predominantly expressed in the roots and flowers and exhibited a concentration-time-specific pattern in response to salt stress. Transgenic lines were generated and subjected to phenotypic characterization under salt stress. The transgenic lines exhibited increased salt tolerance, root length, and fresh weight compared to the wild type. Antioxidant enzyme activity and malondialdehyde content were subsequently measured, and the data revealed no significant differences between the transgenic and wild-type plants in the absence of salt stress. However, under salt stress, the wild-type plants exhibited significantly lower activities of SOD, POD, and CAT than the three transgenic lines, whereas the activity of APX and the content of MDA showed the opposite trend. We identified changes in glutathione pools and associated enzyme activity to gain insights into the underlying mechanisms of the observed phenotypic differences. Notably, under salt stress, the transgenic Arabidopsis's GST activity, GR activity, and GSH content were significantly higher than those of the wild type. In summary, our findings suggest that GmGSTU23 mediates the scavenging of reactive oxygen species and glutathione by enhancing the activity of glutathione transferase, thereby conferring enhanced tolerance to salt stress in plants.

Soybean is a major source of high-quality protein for humans and animals. The content of sulfur-containing amino acids (SAA) in soybean is insufficient, which has become the main factor limiting soybean nutrition. In this study, we used the high-density genetic maps derived from Guizao 1 and Brazil 13 to evaluate the quantitative trait loci of cysteine (Cys), methionine (Met), SAA, glycinin (7S), β-conglycinin (11S), ratio of glycinin to β-conglycinin (RGC), and protein content (PC). In genetic map linkage analysis, the major and stable 44 QTLs were detected, which shared nine bin intervals. Among them, the bin interval (bin157-bin160) on chromosome 5 was detected in multiple environments as a stable QTL, which was linked to 11S, 7S, RGC, and SSA. Based on the analysis of bioinformatics and RNA-sequencing data, 16 differentially expressed genes (DEGs) within these QTLs were selected as candidate genes. These results will help to elucidate the genetic mechanism of soybean SAA-related traits and provide the basis for the gene mining of sulfur-containing amino acids.

Cadmium (Cd) is a widely distributed heavy metal that is toxic to plants and humans. Although silicon (Si) has been reported to reduce Cd accumulation and toxicity in plants, evidence on the functions of Si and its mechanisms in the possible alleviation of soybean are limited. Therefore, a controlled experiment was conducted to investigate the impacts and mechanisms of Si on Cd retention in soybean. Here, we determined the growth index, Cd distribution, and antioxidant activity systems of Si, as well as expression levels of differentially expressed genes (DEGs) in Si under Cd stress, and conducted RNA-seq analysis. We not only found that Si can significantly promote soybean plant growth, increase plant antioxidant activities, and reduce the Cd translocation factor, but also revealed that a total of 636 DEGs were shared between CK and Cd, CK and Cd + Si, and Cd and Cd + Si. Moreover, several genes were significantly enriched in antioxidant systems and Cd distribution and transport systems. Therefore, the expression status of Si-mediated Cd stress response genes is likely involved in improving oxidative stress and changing Cd uptake and transport, as well as improving plant growth that contributes to Si alleviating Cd toxicity in plants. Moreover, numerous potential target genes were identified for the engineering of Cd-tolerant cultivars in soybean breeding programs., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

Aluminum (Al) toxicity is an essential factor that adversely limits soybean ( Glycine max (L.) Merr.) growth in acid soils. WRKY transcription factors play important roles in soybean responses to abiotic stresses. Here, GmWRKY81 was screened from genes that were differentially expressed under Al treatment in Al-tolerant soybean Baxi10 and Al-sensitive soybean Bendi2. We found that GmWRKY81 was significantly induced by 20 μM AlCl 3 and upregulated by AlCl 3 treatment for 2 h. In different tissues, the expression of GmWRKY81 was differentially induced. In 0-1 cm root tips, the expression of GmWRKY81 was induced to the highest level. The overexpression of GmWRKY81 in soybean resulted in higher relative root elongation, root weight, depth, root length, volume, number of root tips and peroxidase activity but lower root average diameter, malonaldehyde and H 2 O 2 contents, indicating enhanced Al tolerance. Moreover, RNA-seq identified 205 upregulated and 108 downregulated genes in GmWRKY81 transgenic lines. Fifteen of these genes that were differentially expressed in both AlCl 3 -treated and GmWRKY81 -overexpressing soybean had the W-box element, which can bind to the upstream-conserved WRKY domain. Overall, the combined functional analysis indicates that GmWRKY81 may improve soybean Al tolerance by regulating downstream genes participating in Al 3+ transport, organic acid secretion and antioxidant reactions.

Background: Deciphering the hereditary mechanism of seed iron (Fe) and zinc (Zn) content in soybean is important and sustainable to address the "hidden hunger" that presently affects approximately 2 billion people worldwide. Therefore, in order to detect genomic regions related to soybean seed Fe and Zn content, a recombinant inbred line (RIL) population with 248 lines was assessed in four environments to detect Quantitative Trait Loci (QTLs) related to soybean seed Fe and Zn content., Result: Wide variation was found in seed Fe and Zn content in four environments, and genotype, environment, and genotype × environment interactions had significant influences on both the seed Fe and Zn content. A positive correlation was observed between seed Fe content and seed Zn content, and broad-sense heritability (H 2 ) of seed Fe and Zn content were 0.73 and 0.75, respectively. In this study, five QTLs for seed Fe content were detected with 4.57 - 32.71% of phenotypic variation explained (PVE) and logarithm of odds (LOD) scores ranging from 3.60 to 33.79. Five QTLs controlling the seed Zn content were detected, and they individually explained 3.35 to 26.48% of the phenotypic variation, with LOD scores ranging from 3.64 to 20.4. Meanwhile, 409,541 high-quality single-nucleotide variants (SNVs) and 85,102 InDels (except intergenic regions) between two bi-parental lines were identified by whole genome resequencing. A total of 12 candidate genes were reported in one major QTL for seed Fe content and two major QTLs for seed Zn content, with the help of RNA-Seq analysis, gene ontology (GO) enrichment, gene annotation, and bi-parental whole genome sequencing (WGS) data., Conclusions: Limited studies were performed about microelement of soybean, so these results may play an important role in the biofortification of Fe and Zn and accelerate the development of marker-assisted selection (MAS) for breeding soybeans fortified with iron and zinc., (© 2022. The Author(s).)

The toxicity of aluminum (Al) in acidic soil limits global crop yield. The ATP-binding cassette (ABC) transporter-like gene superfamily has functions and structures related to transportation, so it responds to aluminum stress in plants. In this study, one half-size ABC transporter gene was isolated from wild soybeans ( Glycine soja ) and designated GsABCI1 . By real-time qPCR, GsABCI1 was identified as not specifically expressed in tissues. Phenotype identification of the overexpressed transgenic lines showed increased tolerance to aluminum. Furthermore, GsABCI1 transgenic plants exhibited some resistance to aluminum treatment by ion translocation or changing root components. This work on the GsABCI1 identified the molecular function, which provided useful information for understanding the gene function of the ABC family and the development of new aluminum-tolerant soybean germplasm.

Soybean, a typical short-day crop, is sensitive to photoperiod, which is a major limiting factor defining its north-to-south cultivation range. The long-juvenile (LJ) trait is controlled primarily by the J locus which has been used for decades by soybean breeders to delay flowering and improve grain yield in tropical regions. The J gene encodes an ortholog of the Arabidopsis Evening Complex (EC) component EARLY FLOWERING 3 (ELF3). To identify modifiers of J, we conducted a forward genetic screen and isolated a mutant (eoj57) that in combination with j has longer flowering delay compared with j single mutant plants. Map-based cloning and genome re-sequencing identified eoj57 (designated as GmLUX2) as an ortholog of the Arabidopsis EC component LUX ARRHYTHMO (LUX). To validate that GmLUX2 is a modifier of J, we used trans-complementation and identified a natural variant allele with a similar phenotype. We also show that GmLUX2 physically interacts with GmELF3a/b and binds DNA, whereas the mutant and natural variant are attenuated in both activities. Transcriptome analysis shows that the GmLUX2-GmELF3a complex co-regulates the expression of several circadian clock-associated genes and directly represses E1 expression. These results provide mechanistic insight into how GmLUX2-GmELF3 controls flowering time via synergistic regulation of gene expression. These novel insights expand our understanding of the regulation of the EC complex, and facilitate the development of soybean varieties adapted for growth at lower latitudes.

Acyl activating enzyme 3 (AAE3) was identified as being involved in the acetylation pathway of oxalate degradation, which regulates the responses to biotic and abiotic stresses in various higher plants. Here, we investigated the role of Glycine soja AAE3 ( GsAAE3 ) in Cadmium (Cd) and Aluminum (Al) tolerances. The recombinant GsAAE3 protein showed high activity toward oxalate, with a K m of 105.10 ± 12.30 μM and V max of 12.64 ± 0.34 μmol min -1 mg -1 protein, suggesting that it functions as an oxalyl-CoA synthetase. The expression of a GsAAE3-green fluorescent protein (GFP) fusion protein in tobacco leaves did not reveal a specific subcellular localization pattern of GsAAE3 . An analysis of the GsAAE3 expression pattern revealed an increase in GsAAE3 expression in response to Cd and Al stresses, and it is mainly expressed in root tips. Furthermore, oxalate accumulation induced by Cd and Al contributes to the inhibition of root growth in wild soybean. Importantly, GsAAE3 overexpression increases Cd and Al tolerances in A. thaliana and soybean hairy roots, which is associated with a decrease in oxalate accumulation. Taken together, our data provide evidence that the GsAAE3 -encoded protein plays an important role in coping with Cd and Al stresses.

This research aimed to identify stable quantitative trait loci (QTL) associated with oil and protein content in soybean. A population of 196 recombinant inbred lines (RILs) derived from Huachun 2 × Wayao was used to evaluate these target traits. A high-density genetic linkage map was constructed by using high-throughput genome-wide sequencing technology, which contained 3413 recombination bin markers and spanned 5400.4 cM with an average distance of 1.58 cM between markers. Eighteen stable QTLs controlling oil and protein content were detected. Among them, qOil-11-1 was identified for the first time as a novel QTL, while qOil-5-1 , qPro-10-1 , and qPro-14-1 were strong and stable QTLs with high log-likelihood (LOD) values. Sixteen differentially expressed genes (DEGs) within these four QTLs were shown to be highly expressed during seed development based on RNA sequencing (RNA-seq) data analysis. Our results may contribute toward gene mining and marker-assisted selection (MAS) for breeding a high-quality soybean in the future.

The MADS-box transcription factors (TFs) are essential in regulating plant growth and development, and conferring abiotic and metal stress resistance. This study aims to investigate GsMAS1 function in conferring tolerance to aluminum stress in Arabidopsis. The GsMAS1 from the wild soybean BW69 line encodes a MADS-box transcription factor in Glycine soja by bioinformatics analysis. The putative GsMAS1 protein was localized in the nucleus. The GsMAS1 gene was rich in soybean roots presenting a constitutive expression pattern and induced by aluminum stress with a concentration-time specific pattern. The analysis of phenotypic observation demonstrated that overexpression of GsMAS1 enhanced the tolerance of Arabidopsis plants to aluminum (Al) stress with larger values of relative root length and higher proline accumulation compared to those of wild type at the AlCl 3 treatments. The genes and/or pathways regulated by GsMAS1 were further investigated under Al stress by qRT-PCR. The results indicated that six genes resistant to Al stress were upregulated, whereas AtALMT1 and STOP2 were significantly activated by Al stress and GsMAS1 overexpression. After treatment of 50 μM AlCl 3 , the RNA abundance of AtALMT1 and STOP2 went up to 17-fold and 37-fold than those in wild type, respectively. Whereas the RNA transcripts of AtALMT1 and STOP2 were much higher than those in wild type with over 82% and 67% of relative expression in GsMAS1 transgenic plants, respectively. In short, the results suggest that GsMAS1 may increase resistance to Al toxicity through certain pathways related to Al stress in Arabidopsis.

The IREG (IRON REGULATED/ferroportin) family of genes plays vital roles in regulating the homeostasis of iron and conferring metal stress. This study aims to identify soybean IREG family genes and characterize the function of GmIREG3 in conferring tolerance to aluminum stress. Bioinformatics and expression analyses were conducted to identify six soybean IREG family genes. One GmIREG , whose expression was significantly enhanced by aluminum stress, GmIREG3 , was studied in more detail to determine its possible role in conferring tolerance to such stress. In total, six potential IREG-encoding genes with the domain of Ferroportin1 (PF06963) were characterized in the soybean genome. Analysis of the GmIREG3 root tissue expression patterns, subcellular localizations, and root relative elongation and aluminum content of transgenic Arabidopsis overexpressing GmIREG3 demonstrated that GmIREG3 is a tonoplast localization protein that increases transgenic Arabidopsis aluminum resistance but does not alter tolerance to Co and Ni. The systematic analysis of the GmIREG gene family reported herein provides valuable information for further studies on the biological roles of GmIREGs in conferring tolerance to metal stress. GmIREG3 contributes to aluminum resistance and plays a role similar to that of FeIREG3 . The functions of other GmIREG genes need to be further clarified in terms of whether they confer tolerance to metal stress or other biological functions.

Aluminum (Al3+) toxicity is a typical abiotic stress that severely limits crop production in acidic soils. In this study, an RIL (recombinant inbred line, F12) population derived from the cross of Zhonghuang 24 (ZH 24) and Huaxia 3 (HX 3) (160 lines) was tested using hydroponic cultivation. Relative root elongation (RRE) and apical Al3+ content (AAC) were evaluated for each line, and a significant negative correlation was detected between the two indicators. Based on a high-density genetic linkage map, the phenotypic data were used to identify quantitative trait loci (QTLs) associated with these traits. With composite interval mapping (CIM) of the linkage map, five QTLs that explained 39.65% of RRE and AAC variation were detected on chromosomes (Chrs) Gm04, Gm16, Gm17 and Gm19. Two new QTLs, qRRE_04 and qAAC_04, were located on the same region of bin93-bin94 on Chr Gm04, which explained 7.09% and 8.98% phenotypic variation, respectively. Furthermore, the results of the expression analysis of candidate genes in the five genetic regions of the QTLs showed that six genes (Glyma.04g218700, Glyma.04g212800, Glyma.04g213300, Glyma.04g217400, Glyma.04g216100 and Glyma.04g220600) exhibited significant differential expression between the Al3+ treatment and the control of two parents. The results of qRT-PCR analysis indicated that Glyma.04g218700 was upregulated by Al3+ treatment with the hundreds-fold increased expression level and may be a candidate gene with potential roles in the response to aluminum stress. Therefore, our efforts will enable future functional analysis of candidate genes and will contribute to the strategies for improvement of aluminum tolerance in soybean., Competing Interests: The authors have declared that no competing interests exist.

Key Message: The Mendelian locus conferring resistance to powdery mildew in soybean was precisely mapped using a combination of phenotypic screening, genetic analyses, and high-throughput genome-wide sequencing. Powdery mildew (PMD), caused by the fungus Microsphaera diffusa Cooke & Peck, leads to considerable yield losses in soybean [Glycine max (L.) Merr.] under favourable environmental conditions and can be controlled by identifying germplasm resources with resistance genes. In this study, resistance to M. diffusa among resistant varieties B3, Fudou234, and B13 is mapped as a single Mendelian locus using three mapping populations derived from crossing susceptible with resistant cultivars. The position of the PMD resistance locus in B3 is located between simple sequence repeat (SSR) markers GMES6959 and Satt_393 on chromosome 16, at genetic distances of 7.1 cM and 4.6 cM, respectively. To more finely map the PMD resistance gene, a high-density genetic map was constructed using 248 F 8 recombinant inbred lines derived from a cross of Guizao1 × B13. The final map includes 3748 bins and is 3031.9 cM in length, with an average distance of 0.81 cM between adjacent markers. This genotypic analysis resulted in the precise delineation of the B13 PMD resistance locus to a 188.06-kb genomic region on chromosome 16 that harbours 28 genes, including 17 disease resistance (R)-like genes in the reference Williams 82 genome. Quantitative real-time PCR assays of possible candidate genes revealed differences in the expression levels of 9 R-like genes between the resistant and susceptible parents. These results provide useful information for marker-assisted breeding and gene cloning for PMD resistance.

Background: The different leaf type associated traits of soybean (Glycine max L.) including leaf area, leaf length, leaf width, leaf shape and petiole length are considered to be associated with seed yield. In order to identify quantitative trait loci (QTLs) affecting leaf type traits, two advanced recombinant inbred line (RIL, ZH, Zhonghuang 24 × Huaxia 3; GB, Guizao 1 × Brazil 13) populations were introduced to score phenotypic values in plants across nine different environments (years, seasons, locations and soybean growth stages). Two restriction site-associated DNA sequencing (RAD-seq) based high-density genetic linkage maps with an average distance of 1.00 centimorgan (cM) between adjacent bin markers were utilized for QTL fine mapping., Results: Correlation analysis showed that most of the traits were correlated with each other and regulated both by hereditary and environmental factors. A total of 190 QTLs were identified for leaf type associated traits in the two populations, of which 14 loci were found to be environmentally stable. Moreover, these detected QTLs were categorized into 34 QTL hotspots, and four important QTL hotspots with phenotypic variance ranging from 3.89-23.13% were highlighted. Furthermore, Glyma04g05840, Glyma19g37820, Glyma14g07140 and Glyma19g39340 were predicted in the intervals of the stable loci and important QTL hotspots for leaf type traits by adopting Gene Ontology (GO) enrichment analysis., Conclusions: Our findings of the QTLs and the putative genes will be beneficial to gain new insights into the genetic basis for soybean leaf type traits and may further accelerate the breeding process for reasonable leaf type soybean.

Background: Multidrug and toxic compound extrusion (MATE) transporters, which exist widely in plants, function as crucial regulators in plant resistance to aluminum (Al) toxicity by inducing citrate efflux. However, the functions of most MATE family members in soybean (Glycine soja) remain to be elucidated., Results: Expression pattern analysis showed that GsMATE was constitutively expressed in different soybean organs, with the highest level in root compared with those in stem, leaf and cotyledon. In addition, Al stress induced expression of GsMATE in soybean. Temporal analysis indicated that GsMATE expression was greatly enhanced by increasing concentrations of aluminum [Al 3+ ] after short exposure, reaching the high levels detected in the BW69 (Al-resistant) and the JW81 (Al-sensitive) lines of Glycine soja of wild soybean at 6 h and 8 h, respectively. Furthermore, transient GsMATE expression in Arabidopsis protoplasts showed that GsMATE protein localized to the plasma membrane. Overexpression of GsMATE on an Arabidopsis columbia-0 (Col-0) background resulted in increased Al tolerance in transgenic plants. Analysis of hematoxylin staining showed that the roots of GsMATE transgenic lines were stained less intensely than those of the wild-type exposured to the same AlCl 3 concentrations. Therefore, GsMATE enhanced the resistance of transgenic plants to Al toxicity by reducing Al accumulation in Arabidopsis roots., Conclusions: In summary, our results indicate that GsMATE is responsive to aluminum stress and may participate in the regulation of sensitivity to Al toxicity in Arabidopsis. In addition, the GsMATE protein is an Al-induced citrate transporter of the MATE family and exerts an essential role in Al tolerance in Glycine soja.

Although sugarcane-soybean intercropping has been widely used to control disease and improve productivity in the field, the response of soil fungal communities to intercropping has not been fully understood. In this study, the rhizosphere fungal communities of sugarcane and soybean under monoculture and intercropping systems were investigated using Illumina MiSeq sequencing of ITS gene. Intercropping decreased the alpha-diversity and changed fungal community composition compared to monocultures. Taxonomic analyses showed that the dominant phyla were Ascomycota, Zygomycota and Basidiomycota. The abundance of Ascomycota decreased in intercropping sugarcane-grown soil compared to monoculture, while it increased in soybean-grown soil in the intercropping system. In addition, intercropping increased the abundance of important fungal genera, such as Trichoderma, Hypocreales and Fusarium but decreased the relative abundance of Gibberella and Chaetomium. The results of canonical correspondence analysis and automatic linear modelling indicated that fungal community compositions were closely associated with soil parameters such as total nitrogen (TN), soil organic matter (SOC), pH and NO 3 - , which suggests that the impacts of intercropping on the soil fungal community are linked to the alteration of soil chemical properties.

Key Message: Map-based cloning identified GmHAD1, a gene which encodes a HAD-like acid phosphatase, associated with soybean tolerance to low phosphorus stress. Phosphorus (P) deficiency in soils is a major limiting factor for crop growth worldwide. Plants may adapt to low phosphorus (LP) conditions via changes to root morphology, including the number, length, orientation, and branching of the principal root classes. To elucidate the genetic mechanisms for LP tolerance in soybean, quantitative trait loci (QTL) related to root morphology responses to LP were identified via hydroponic experiments. In total, we identified 14 major loci associated with these traits in a RIL population. The log-likelihood scores ranged from 2.81 to 7.43, explaining 4.23-13.98% of phenotypic variance. A major locus on chromosome 08, named qP8-2, was co-localized with an important P efficiency QTL (qPE8), containing phosphatase genes GmACP1 and GmACP2. Another major locus on chromosome 10 named qP10-2 explained 4.80-13.98% of the total phenotypic variance in root morphology. The qP10-2 contains GmHAD1, a gene which encodes an acid phosphatase. In the transgenic soybean hairy roots, GmHAD1 overexpression increased P efficiency by 8.4-16.5% relative to the control. Transgenic Arabidopsis plants had higher biomass than wild-type plants across both short- and long-term P reduction. These results suggest that GmHAD1, an acid phosphatase gene, improved the utilization of organic phosphate by soybean and Arabidopsis plants.

Key Message: Fifteen stable QTLs were identified using a high-density soybean genetic map across multiple environments. One major QTL, qIF5-1, contributing to total isoflavone content explained phenotypic variance 49.38, 43.27, 46.59, 45.15 and 52.50%, respectively. Soybeans (Glycine max L.) are a major source of dietary isoflavones. To identify novel quantitative trait loci (QTL) underlying isoflavone content, and to improve the accuracy of marker-assisted breeding in soybean, a valuable mapping population comprised of 196 F 7:8-10 recombinant inbred lines (RILs, Huachun 2 × Wayao) was utilized to evaluate individual and total isoflavone content in plants grown in four different environments in Guangdong. A high-density genetic linkage map containing 3469 recombination bin markers based on 0.2 × restriction site-associated DNA tag sequencing (RAD-seq) technology was used to finely map QTLs for both individual and total isoflavone contents. Correlation analyses showed that total isoflavone content, and that of five individual isoflavone, was significantly correlated across the four environments. Based on the high-density genetic linkage map, a total of 15 stable quantitative trait loci (QTLs) associated with isoflavone content across multiple environments were mapped onto chromosomes 02, 05, 07, 09, 10, 11, 13, 16, 17, and 19. Further, one of them, qIF5-1, localized to chromosomes 05 (38,434,171-39,045,620 bp) contributed to almost all isoflavone components across all environments, and explained 6.37-59.95% of the phenotypic variance, especially explained 49.38, 43.27, 46.59, 45.15 and 52.50% for total isoflavone. The results obtained in the present study will pave the way for a better understanding of the genetics of isoflavone accumulation and reveals the scope available for improvement of isoflavone content through marker-assisted selection.

Background: One of the overarching goals of soybean breeding is to develop lines that combine increased yield with improved quality characteristics. High-density-marker QTL mapping can serve as an effective strategy to identify novel genomic information to facilitate crop improvement. In this study, we genotyped a recombinant inbred line (RIL) population (Zhonghuang 24 × Huaxia 3) using a restriction-site associated DNA sequencing (RAD-seq) approach. A high-density soybean genetic map was constructed and used to identify several QTLs that were shown to influence six yield-related and two quality traits., Results: A total of 47,472 single-nucleotide polymorphisms (SNPs) were detected for the RILs that were integrated into 2639 recombination bin units, with an average distance of 1.00 cM between adjacent markers. Forty seven QTLs for yield-related traits and 13 QTLs for grain quality traits were found to be distributed on 16 chromosomes in the 2 year studies. Among them, 18 QTLs were stable, and were identified in both analyses. Twenty six QTLs were identified for the first time, with a single QTL (qNN19a) in a 56 kb region explaining 32.56% of phenotypic variation, and an additional 10 of these were novel, stable QTLs. Moreover, 8 QTL hotpots on four different chromosomes were identified for the correlated traits., Conclusions: With RAD-sequencing, some novel QTLs and important QTL clusters for both yield-related and quality traits were identified based on a new, high-density bin linkage map. Three predicted genes were selected as candidates that likely have a direct or indirect influence on both yield and quality in soybean. Our findings will be helpful for understanding common genetic control mechanisms of co-localized traits and to select cultivars for further analysis to predictably modulate soybean yield and quality simultaneously.

Key Message: Using a combination of phenotypic screening, genetic and statistical analyses, and high-throughput genome-wide sequencing, we have finely mapped a dominant Phytophthora resistance gene in soybean cultivar Wayao. Phytophthora root rot (PRR) caused by Phytophthora sojae is one of the most important soil-borne diseases in many soybean-production regions in the world. Identification of resistant gene(s) and incorporating them into elite varieties are an effective way for breeding to prevent soybean from being harmed by this disease. Two soybean populations of 191 F 2 individuals and 196 F 7:8 recombinant inbred lines (RILs) were developed to map Rps gene by crossing a susceptible cultivar Huachun 2 with the resistant cultivar Wayao. Genetic analysis of the F 2 population indicated that PRR resistance in Wayao was controlled by a single dominant gene, temporarily named RpsWY, which was mapped on chromosome 3. A high-density genetic linkage bin map was constructed using 3469 recombination bins of the RILs to explore the candidate genes by the high-throughput genome-wide sequencing. The results of genotypic analysis showed that the RpsWY gene was located in bin 401 between 4466230 and 4502773 bp on chromosome 3 through line 71 and 100 of the RILs. Four predicted genes (Glyma03g04350, Glyma03g04360, Glyma03g04370, and Glyma03g04380) were found at the narrowed region of 36.5 kb in bin 401. These results suggest that the high-throughput genome-wide resequencing is an effective method to fine map PRR candidate genes.

Background: The P1B-ATPase subfamily is an important group involved in transporting heavy metals and has been extensively studied in model plants, such as Arabidopsis thaliana and Oryza sativa. Emerging evidence indicates that one homolog in Glycine max is also involved in cadmium (Cd) stress, but the gene family has not been fully investigated in soybean., Results: Here, we identified 20 heavy metal ATPase (HMA) family members in the soybean genome, presented as 10 paralogous pairs, which is significantly greater than the number in Arabidopsis or rice, and was likely caused by the latest whole genome duplication event in soybean. A phylogenetic analysis divided the 20 members into six groups, each having conserved or divergent gene structures and protein motif patterns. The integration of RNA-sequencing and qRT-PCR data from multiple tissues provided an overall expression pattern for the HMA family in soybean. Further comparisons of expression patterns and the single nucleotide polymorphism distribution between paralogous pairs suggested functional conservation and the divergence of HMA genes during soybean evolution. Finally, analyses of the HMAs expressed in response to Cd stress provided evidence on how plants manage Cd tolerance, at least in the two contrasting soybean genotypes examined., Conclusions: The genome-wide identification, chromosomal distribution, gene structures, and evolutionary and expression analyses of the 20 HMA genes in soybean provide an overall insight into their potential involvement in Cd responses. These results will facilitate further research on the HMA gene family, and their conserved and divergent biological functions in soybean.

Background: The precise timing of flowering is fundamental to successful reproduction, and has dramatic significance for crop yields. Although prolonged low temperatures are not required for flowering induction in soybean, vernalization pathway genes have been retained during the evolution of this species. Little information is currently available in regarding these genes in soybean., Results: We were able to detect the expression of Glyma11g13220 in different organs at all monitored developmental stages in soybean. Glyma11g13220 expression was higher in leaves and pods than in shoot apexes and stems. In addition, Glyma11g13220 was responsive to photoperiod and low temperature in soybean. Furthermore, Glyma11g13220 was found to be a nuclear-localized protein. Over-expression of Glyma11g13220 in an Arabidopsis Columbia-0 (Col-0) background resulted in early flowering. Quantitative real-time PCR analysis revealed that transcript levels of flower repressor FLOWERING LOCUS C (FLC), and FD decreased significantly in transgenic Arabidopsis compared with wild-type Col-0, while the expression of VERNALIZATION INSENSITIVE 3 (VIN3) and FLOWERING LOCUS T (FT) noticeably increased., Conclusions: Our results suggest that Glyma11g13220, a homolog of Arabidopsis VRN1, is a functional protein. Glyma11g13220, which is responsive to photoperiod and low temperature in soybean, may participate in the vernalization pathway in Arabidopsis and help regulate flowering time. Arabidopsis VRN1 and Glyma11g13220 exhibit conserved as well as diverged functions.

Background: Soybean is one of the most important crops, providing large amounts of dietary proteins and edible oil, and is also an excellent model for studying evolution of duplicated genes. However, relative to the model plants Arabidopsis and rice, the present knowledge about soybean transcriptome is quite limited., Results: In this study, we employed RNA-seq to investigate transcriptomes of 11 soybean tissues, for genome-wide discovery of truly expressed genes, and novel and alternative transcripts, as well as analyses of conservation and divergence of duplicated genes and their functional implications. We detected a total of 54,132 high-confidence expressed genes, and identified 6,718 novel transcriptional regions with a mean length of 372 bp. We also provided strong evidence for alternative splicing (AS) events for ~15.9% of the genes with two or more exons. Among them, 1,834 genes exhibited stage-dependent AS, and 202 genes had tissue-biased exon-skipping events. We further defined the conservation and divergence in expression patterns between duplicated gene pairs from recent whole genome duplications (WGDs); differentially expressed genes, tissue preferentially expressed genes, transcription factors and specific gene family members were identified for shoot apical meristem and flower development., Conclusions: Our results significantly improved soybean gene annotation, and also provide valuable resources for functional genomics and studies of the evolution of duplicated genes from WGDs in soybean.

MicroRNAs (miRNAs) play crucial roles in regulating the expression of various stress responses genes in plants. To investigate soybean (Glycine max) miRNAs involved in the response to cadmium (Cd), microarrays containing 953 unique miRNA probes were employed to identify differences in the expression patterns of the miRNAs between different genotypes, Huaxia3 (HX3, Cd-tolerant) and Zhonghuang24 (ZH24, Cd-sensitive). Twenty six Cd-responsive miRNAs were identified in total. Among them, nine were detected in both cultivars, while five were expressed only in HX3 and 12 were only in ZH24. The expression of 16 miRNAs was tested by qRT-PCR and most of the identified miRNAs were found to have similar expression patterns with microarray. Three hundred and seventy six target genes were identified for 204 miRNAs from a mixture degradome library, which was constructed from the root of HX3 and ZH24 with or without Cd treatment. Fifty five genes were identified to be cleaved by 14 Cd-responsive miRNAs. Gene ontology (GO) annotations showed that these target transcripts are implicated in a broad range of biological processes. In addition, the expression patterns of ten target genes were validated by qRT-PCR. The characterization of the miRNAs and the associated target genes in response to Cd exposure provides a framework for understanding the molecular mechanism of heavy metal tolerance in plants.

Gibberellic acids (GAs) are plant hormones that play fundamental roles in plant growth and developmental processes. Previous studies have demonstrated that three key enzymes of GA20ox, GA3ox, and GA2ox are involved in GA biosynthesis. In this study, the Arabidopsis DREB1A gene driven by the CaMV 35S promoter was introduced into soybean plants by Agrobacterium- mediated transformation. The results showed that the transgenic soybean plants exhibited a typical phenotype of GA-deficient mutants, such as severe dwarfism, small and dark-green leaves, and late flowering compared to those of the non-transgenic plants. The dwarfism phenotype was rescued by the application of exogenous GA(3) once a week for three weeks with the concentrations of 144 µM or three times in one week with the concentrations of 60 µM. Quantitative RT-PCR analysis revealed that the transcription levels of the GA synthase genes were higher in the transgenic soybean plants than those in controls, whereas GA-deactivated genes except GmGA2ox4 showed lower levels of expression. The transcript level of GmGA2ox4 encoding the only deactivation enzyme using C(20)-GAs as the substrates in soybean was dramatically enhanced in transgenic plants compared to that of wide type. Furthermore, the contents of endogenous bioactive GAs were significantly decreased in transgenic plants than those of wide type. The results suggested that AtDREB1A could cause dwarfism mediated by GA biosynthesis pathway in soybean.

Soybean is a short-day plant and its flowering process can be reversed when switching from short-day to long-day conditions. Flowering reversion provides a useful system to study the flowering process in both forward and backward directions. In this study, we optimized a soybean flowering reversion system using a photoperiod-sensitive cultivar Zigongdongdou. Three types of terminal structures were found during flowering reversion: reversed terminal raceme (RTR), short terminal raceme (STR), and vegetative terminal (VT). The relative frequency of these terminal structures during flowering reversion under long day was dependent on the duration of the prior short day (SD) pretreatment. This process is phytochrome dependent and young plants were more susceptible to flowering reversion. Leaf removal increased the minimal SD period needed for the induction of STR. To demonstrate the application of this system, we studied the patterns of in situ expression of the GmNMH7 gene during flowering development and reversion. NMH7 family members encode MADS-box proteins and are unique in legume families since their expression can be detected in both developing flowers and nodules. In situ hybridization experiments using plants grown under different photoperiod cycles provided several lines of evidence supporting a close relationship between GmNMH7 gene expression and floral development in soybean. Furthermore, it seems that GmNMH7 may participate in flower development at different stages. Interestingly, the expression pattern of GmNMH7 in root nodules was also found to be regulated by photoperiod. These results support the notion that the photoperiod sensitive GmNMH7 gene may play multiple roles in growth and development in soybean.