ALBERT

Description: Plant height (PH) is an important agronomic trait that is closely related to soybean yield and quality. However, it is a complex quantitative trait governed by multiple genes and is influenced by environment. Unraveling the genetic mechanism involved in PH, and developing soybean cultivars with desirable PH is an imperative goal for soybean breeding. In this regard, the present study used high-density linkage maps of two related recombinant inbred line (RIL) populations viz., MT and ZM evaluated in three different environments to detect additive and epistatic effect quantitative trait loci (QTLs) as well as their interaction with environments for PH in Chinese summer planting soybean. A total of eight and 12 QTLs were detected by combining the composite interval mapping (CIM) and mixed-model based composite interval mapping (MCIM) methods in MT and ZM populations, respectively. Among these QTLs, nine QTLs viz., QPH-2, qPH-6-2MT, QPH-6, qPH-9-1ZM, qPH-10-1ZM, qPH-13-1ZM, qPH-16-1MT, QPH-17 and QPH-19 were consistently identified in multiple environments or populations, hence were regarded as stable QTLs. Furthermore, Out of these QTLs, three QTLs viz., qPH-4-2ZM, qPH-15-1MT and QPH-17 were novel. In particular, QPH-17 could detect in both populations, which was also considered as a stable and major QTL in Chinese summer planting soybean. Moreover, eleven QTLs revealed significant additive effects in both populations, and out of them only six showed additive by environment interaction effects, and the environment-independent QTLs showed higher additive effects. Finally, six digenic epistatic QTLs pairs were identified and only four additive effect QTLs viz., qPH-6-2MT, qPH-19-1MT/QPH-19, qPH-5-1ZM and qPH-17-1ZM showed epistatic effects. These results indicate that environment and epistatic interaction effects have significant influence in determining genetic basis of PH in soybean. These results would not only increase our understanding of the genetic control of plant height in summer planting soybean but also provide support for implementing marker assisted selection (MAS) in developing cultivars with ideal plant height as well as gene cloning to elucidate the mechanisms of plant height.

Description: Abnormal reactive oxygen species (ROS) may mediate cytoplasmic male sterility (CMS). To observe the effect of ROS on soybean CMS, metabolite content and antioxidant enzyme activity in the flower buds between soybean N8855-derived CMS line and its maintainer were compared. Of the 612 metabolites identified, a total of 74 metabolites were significantly differentiated in flower buds between CMS line and its maintainer. The differential metabolites involved 32 differential flavonoids, 13 differential phenolamides, and 1 differential oxidized glutathione (GSSG) belonging to a non-enzymatic ROS scavenging system. We observed lower levels of flavonoids and antioxidant enzyme activities in flower buds of the CMS line than in its maintainer. Our results suggest that deficiencies of enzymatic and non-enzymatic ROS scavenging systems in soybean CMS line cannot eliminate ROS in anthers effectively, excessive accumulation of ROS triggered programmed cell death and ultimately resulted in pollen abortion of soybean CMS line.

Description: Seed protein and oil content are the two important traits determining the quality and value of soybean. Development of improved cultivars requires detailed understanding of the genetic basis underlying the trait of interest. However, it is prerequisite to have a high-density linkage map for precisely mapping genomic regions, and therefore the present study used high-density genetic map containing 2267 recombination bin markers distributed on 20 chromosomes and spanned 2453.79 cM with an average distance of 1.08 cM between markers using restriction-site-associated DNA sequencing (RAD-seq) approach. A recombinant inbred line (RIL) population of 104 lines derived from a cross between Linhefenqingdou and Meng 8206 cultivars was evaluated in six different environments to identify main- and epistatic-effect quantitative trait loci (QTLs)as well as their interaction with environments. A total of 44 main-effect QTLs for protein and oil content were found to be distributed on 17 chromosomes, and 15 novel QTL were identified for the first time. Out of these QTLs, four were major and stable QTLs, viz., qPro-7-1, qOil-8-3, qOil-10-2 and qOil-10-4, detected in at least two environments plus combined environment with R2 values 〉10%. Within the physical intervals of these four QTLs, 111 candidate genes were screened for their direct or indirect involvement in seed protein and oil biosynthesis/metabolism processes based on gene ontology and annotation information. Based on RNA sequencing (RNA-seq) data analysis, 15 of the 111 genes were highly expressed during seed development stage and root nodules that might be considered as the potential candidate genes. Seven QTLs associated with protein and oil content exhibited significant additive and additive × environment interaction effects, and environment-independent QTLs revealed higher additive effects. Moreover, three digenic epistatic QTLs pairs were identified, and no main-effect QTLs showed epistasis. In conclusion, the use of a high-density map identified closely linked flanking markers, provided better understanding of genetic architecture and candidate gene information, and revealed the scope available for improvement of soybean quality through marker assisted selection (MAS).

Description: Soybean is one of the most important oil crops in the world. Revealing the molecular basis and exploring key candidate genes for seed oil synthesis has great significance for soybean improvement. In this study, we found that oil accumulation rates and gene expression levels changed dynamically during soybean seed development. The expression levels of genes in metabolic pathways such as carbon fixation, photosynthesis, glycolysis, and fatty acid biosynthesis were significantly up-regulated during the rapid accumulation of oil in developing soybean seeds. Through weighted correlation network analysis, we identified six co-expression modules associated with soybean seed oil content and the pink module was the most positively correlated (r = 0.83, p = 7 × 10−4) network. Through the integration of differential expression and co-expression analysis, we predicted 124 candidate genes potentially affecting soybean seed oil content, including seven genes in lipid metabolism pathway, two genes involved in glycolysis, one gene in sucrose metabolism, and 12 genes belonged to transcription factors as well as other categories. Among these, three genes (GmABI3b, GmNFYA and GmFAD2-1B) have been shown to control oil and fatty acid content in soybean seeds, and other newly identified candidate genes would broaden our knowledge to understand the molecular basis for oil accumulation in soybean seeds.

Description: Peroxidases play prominent roles in antioxidant responses and stress tolerance in plants; however, their functions in soybean tolerance to salt stress remain unclear. Here, we investigated the role of a peroxidase gene from the wild soybean (Glycine soja), GsPRX9, in soybean tolerance to salt stress. GsPRX9 gene expression was induced by salt treatment in the roots of both salt-tolerant and -sensitive soybean varieties, and its relative expression level in the roots of salt-tolerant soybean varieties showed a significantly higher increase than in salt-sensitive varieties after NaCl treatment, suggesting its possible role in soybean response to salt stress. GsPRX9-overexpressing yeast (strains of INVSc1 and G19) grew better than the control under salt and H2O2 stress, and GsPRX9-overexpressing soybean composite plants showed higher shoot fresh weight and leaf relative water content than control plants after NaCl treatment. Moreover, the GsPRX9-overexpressing soybean hairy roots had higher root fresh weight, primary root length, activities of peroxidase and superoxide dismutase, and glutathione level, but lower H2O2 content than those in control roots under salt stress. These findings suggest that the overexpression of the GsPRX9 gene enhanced the salt tolerance and antioxidant response in soybean. This study would provide new insights into the role of peroxidase in plant tolerance to salt stress.

Description: Drought is one of the most important factors affecting plant growth and productivity. The previous results on drought tolerance (DT) genetic system in soybean indicated a complex of genes not only few ones were involved in the trait. This study is featured with a relatively thorough identification of QTL-allele/candidate-gene system using an efficient restricted two-stage multi-locus multi-allele genome-wide association study, on two comprehensive DT indicators, membership index values of relative plant weight (MPW) and height (MPH), instead of a single biological characteristic, in a large sample (564 accessions) of the Chinese cultivated soybean population (CCSP). Based on 24,694 multi-allele markers, 75 and 64 QTL with 261 and 207 alleles (2–12/locus) were detected for MPW and MPH, explaining 54.7% and 47.1% of phenotypic variance, respectively. The detected QTL-alleles were organized into a QTL-allele matrix for each indicator, indicating DT is a super-trait conferred by two (even more) QTL-allele systems of sub-traits. Each CCSP matrix was separated into landrace (LR) and released cultivar (RC) sub-matrices, which showed significant differentiation in QTL-allele constitutions, with 58 LR alleles excluded and 16 new ones emerged in RC. Using the matrices, optimal crosses with great DT transgressive recombinants were predicted. From the detected QTL, 177 candidate genes were annotated and validated with quantitative Real-time PCR, and grouped into nine categories, with ABA and stress responders as the major parts. The key point of the above results is the establishment of relatively full QTL-allele matrices composed of numerous gene functions jointly conferring DT, therefore, demonstrates the complexity of DT genetic system and potential of CCSP in DT breeding.

Description: Aluminum (Al) toxicity is a major factor limiting crop productivity on acid soils. Soybean (Glycine max) is an important oil crop and there is great variation in Al tolerance in soybean germplasms. However, only a few Al-tolerance genes have been reported in soybean. Therefore, the purpose of this study was to identify candidate Al tolerance genes by comparative transcriptome analysis of two contrasting soybean varieties in response to Al stress. Two soybean varieties, M90-24 (M) and Pella (P), which showed significant difference in Al tolerance, were used for RNA-seq analysis. We identified a total of 354 Al-tolerance related genes, which showed up-regulated expression by Al in the Al-tolerant soybean variety M and higher transcript levels in M than P under Al stress. These genes were enriched in the Gene Ontology (GO) terms of cellular glucan metabolic process and regulation of transcription. Five out of 11 genes in the enriched GO term of cellular glucan metabolic process encode cellulose synthases, and one cellulose synthase gene (Glyma.02G205800) was identified as the key hub gene by co-expression network analysis. Furthermore, treatment of soybean roots with a cellulose biosynthesis inhibitor decreased the Al tolerance, indicating an important role of cellulose production in soybean tolerance to Al toxicity. This study provides a list of candidate genes for further investigation on Al tolerance mechanisms in soybean.

Description: Yield evaluation of breeding lines is the key to successful release of cultivars, which is becoming a serious issue due to soil heterogeneity in enlarged field tests. This study aimed at establishing plot-yield prediction models using unmanned aerial vehicle (UAV)-based hyperspectral remote sensing for yield-selection in large-scale soybean breeding programs. Three sets of soybean breeding lines (1103 in total) were tested in blocks-in-replication experiments for plot yield and canopy spectral reflectance on 454~950 nm bands at different growth stages using a UAV-based hyperspectral spectrometer (Cubert UHD185 Firefly). The four elements for plot-yield prediction model construction were studied respectively and concluded as: the suitable reflectance-sampling unit-size in a plot was its 20%–80% central part; normalized difference vegetation index (NDVI) and ration vegetation index (RVI) were the best combination of vegetation indices; the initial seed-filling stage (R5) was the best for a single stage prediction, while another was the best combination for a two growth-stage prediction; and multi-variate linear regression was suitable for plot-yield prediction. In model establishment for each material-set, a random half was used for modelling and another half for verification. Twenty-one two growth-stage two vegetation-index prediction models were established and compared for their modelling coefficient of determination (RM2) and root mean square error of the model (RMSEM), verification RV2 and RMSEV, and their sum RS2 and RMSES. Integrated with the coincidence rate between the model predicted and the practical yield-selection results, the models, MA1-2, MA4-2 and MA6-2, with coincidence rates of 56.8%, 58.5% and 52.4%, respectively, were chosen for yield-prediction in yield-test nurseries. The established model construction elements and methods can be used as local models for pre-harvest yield-selection and post-harvest integrated yield-selection in advanced breeding nurseries as well as yield potential prediction in plant-derived-line nurseries. Furthermore, multiple models can be used jointly for plot-yield prediction in soybean breeding programs.

Description: Aluminum (Al) toxicity is an important barrier to soybean (Glycine max (L.) Merr.) production in acid soils. However, little is known about the genes underlying Al tolerance in soybean. We aim to find the key candidate genes and investigate their roles in soybean tolerance to Al toxicity in this study. Comparative transcriptome analyses of the Al-tolerant (KF) and Al-sensitive (GF) soybean varieties under control and Al stress at 6, 12, and 24 h were investigated. A total of 1411 genes showed specific up-regulation in KF or more up-regulation in KF than in GF by Al stress, which were significantly enriched in the GO terms of peroxidase (POD) activity, transporter activity (including the known Al tolerance-related ABC transporter, ALMT, and MATE), and four families of transcription factors (AP2, C3H4, MYB, WRKY). The expression levels of seven POD genes were up-regulated by Al stress for at least one time point in KF. The H2O2 pretreatment significantly improved Al tolerance of KF, which is likely due to increased POD activity induced by H2O2. Our results suggest that PODs play important roles in soybean tolerance to Al toxicity. We also propose a list of candidate genes for Al tolerance in KF, which would provide valuable insights into the Al tolerance mechanisms in soybean.

Description: Annual wild soybean (G. soja) is the ancestor of the cultivated soybean (G. max). To reveal the genetic changes from soja to max, an improved wild soybean chromosome segment substitution line (CSSL) population, SojaCSSLP5, composed of 177 CSSLs with 182 SSR markers (SSR-map), was developed based on SojaCSSLP1 generated from NN1138-2(max)×N24852(soja). The SojaCSSLP5 was genotyped further through whole-genome resequencing, resulting in a physical map with 1366 SNPLDBs (SNP linkage-disequilibrium blocks), which are composed of more markers/segments, shorter marker length and more recombination breakpoints than the SSR-map and caused 721 new wild substituted segments. Using the SNPLDB-map, two loci co-segregating with seed-coat color (SCC) and six loci for days to flowering (DTF) with 88.02% phenotypic contribution were identified. Integrated with parental RNA-seq and DNA-resequencing, two SCC and six DTF candidate genes, including three previously cloned (G, E2 and GmPRR3B) and five newly detected ones, were predicted and verified at nucleotide mutant level, and then demonstrated with the consistency between gene-alleles and their phenotypes in SojaCSSLP5. In total, six of the eight genes were identified with the parental allele-pairs coincided to those in 303 germplasm accessions, then were further demonstrated by the consistency between gene-alleles and germplasm phenotypes. Accordingly, the CSSL population integrated with parental DNA and RNA sequencing data was demonstrated to be an efficient platform in identifying candidate wild vs. cultivated gene-alleles.