Deciphering the genetic basis of resistance to soybean cyst nematode combining IBD and association mapping.

KEY MESSAGE: IBD analysis clarified the dynamics of chromosomal recombination during the ZP pedigree breeding process and identified ten genomic regions resistant to SCN race3 combining association mapping. Soybean cyst nematode (SCN, Heterodera glycines Ichinohe) is one of the most devastating pathogens for soybean production worldwide. The cultivar Zhongpin03-5373 (ZP), derived from SCN-resistant progenitor parents, Peking, PI 437654 and Huipizhi Heidou, is an elite line with high resistance to SCN race3. In the current study, a pedigree variation map was generated for ZP and its ten progenitors using 3,025,264 high-quality SNPs identified from an average of 16.2 × re-sequencing for each genome. Through identity by decent (IBD) tracking, we showed the dynamic change of genome and detected important IBD fragments, which revealed the comprehensively artificial selection of important traits during ZP breeding process. A total of 2,353 IBD fragments related to SCN resistance including SCN-resistant genes rhg1, rhg4 and NSFRAN07 were identified based on the resistant-related genetic paths. Moreover, 23 genomic regions underlying resistance to SCN race3 were identified by genome-wide association study (GWAS) in 481 re-sequenced cultivated soybeans. Ten common loci were found by both IBD tracking and GWAS analysis. Haplotype analysis of 16 potential candidate genes suggested a causative SNP (C/T, - 1065) located in the promoter of Glyma.08G096500 and encoding a predicted TIFY5b-related protein on chr8 was highly correlated with SCN race3 resistance. Our results more thoroughly elucidated the dynamics of genomic fragments during ZP pedigree breeding and the genetic basis of SCN resistance, which will provide useful information for gene cloning and the development of resistant soybean cultivars using a marker-assisted selection approach.

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.

Time series canopy phenotyping enables the identification of genetic variants controlling dynamic phenotypes in soybean.

Advances in plant phenotyping technologies are dramatically reducing the marginal costs of collecting multiple phenotypic measurements across several time points. Yet, most current approaches and best statistical practices implemented to link genetic and phenotypic variation in plants have been developed in an era of single-time-point data. Here, we used time-series phenotypic data collected with an unmanned aircraft system for a large panel of soybean (Glycine max (L.) Merr.) varieties to identify previously uncharacterized loci. Specifically, we focused on the dissection of canopy coverage (CC) variation from this rich data set. We also inferred the speed of canopy closure, an additional dimension of CC, from the time-series data, as it may represent an important trait for weed control. Genome-wide association studies (GWASs) identified 35 loci exhibiting dynamic associations with CC across developmental stages. The time-series data enabled the identification of 10 known flowering time and plant height quantitative trait loci (QTLs) detected in previous studies of adult plants and the identification of novel QTLs influencing CC. These novel QTLs were disproportionately likely to act earlier in development, which may explain why they were missed in previous single-time-point studies. Moreover, this time-series data set contributed to the high accuracy of the GWASs, which we evaluated by permutation tests, as evidenced by the repeated identification of loci across multiple time points. Two novel loci showed evidence of adaptive selection during domestication, with different genotypes/haplotypes favored in different geographic regions. In summary, the time-series data, with soybean CC as an example, improved the accuracy and statistical power to dissect the genetic basis of traits and offered a promising opportunity for crop breeding with quantitative growth curves.

Soybean genetic resources contributing to sustainable protein production.

KEY MESSAGE: Genetic resources contributes to the sustainable protein production in soybean. Soybean is an important crop for food, oil, and forage and is the main source of edible vegetable oil and vegetable protein. It plays an important role in maintaining balanced dietary nutrients for human health. The soybean protein content is a quantitative trait mainly controlled by gene additive effects and is usually negatively correlated with agronomic traits such as the oil content and yield. The selection of soybean varieties with high protein content and high yield to secure sustainable protein production is one of the difficulties in soybean breeding. The abundant genetic variation of soybean germplasm resources is the basis for overcoming the obstacles in breeding for soybean varieties with high yield and high protein content. Soybean has been cultivated for more than 5000 years and has spread from China to other parts of the world. The rich genetic resources play an important role in promoting the sustainable production of soybean protein worldwide. In this paper, the origin and spread of soybean and the current status of soybean production are reviewed; the genetic characteristics of soybean protein and the distribution of resources are expounded based on phenotypes; the discovery of soybean seed protein-related genes as well as transcriptomic, metabolomic, and proteomic studies in soybean are elaborated; the creation and utilization of high-protein germplasm resources are introduced; and the prospect of high-protein soybean breeding is described.

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.

QTL Mapping of Palmitic Acid Content Using Specific-Locus Amplified Fragment Sequencing (SLAF-Seq) Genotyping in Soybeans (Glycine max L.).

Soybeans are essential crops that supply protein and oil. The composition and contents of soybean fatty acids are relevant to human health and have a significant relationship with soybean oil processing and applications. Identifying quantitative trait locus (QTL) genes related to palmitic acid could facilitate the development of a range of nutritive soybean cultivars using molecular marker-assisted selection. In this study, we used a cultivar with higher palmitic acid content, 'Dongnong42', and a lower palmitic acid content cultivar, 'Hobbit', to establish F2:6 recombinant inbred lines. A high-density genetic map containing 9980 SLAF markers was constructed and distributed across 20 soybean chromosomes. The genetic map contained a total genetic distance of 2602.58 cM and an average genetic distance of 0.39 cM between adjacent markers. Two QTLs related to palmitic acid content were mapped using inclusive composite interval mapping, explaining 4.2-10.1% of the phenotypic variance in three different years and environments, including the QTL included in seed palmitic 7-3, which was validated by developing SSR markers. Based on the SNP/Indel and significant differential expression analyses of Dongnong42 and Hobbit, two genes, Glyma.15g119700 and Glyma.15g119800, were selected as candidate genes. The high-density genetic map, QTLs, and molecular markers will be helpful for the map-based cloning of palmitic acid content genes. These could be used to accelerate breeding for high nutritive value cultivars via molecular marker-assisted breeding.

Molecular characterization of common bean curly stunt virus: a novel recombinant geminivirus in China.

A new geminivirus was identified in common bean (Phaseolus vulgaris) showing severe stunt and leaf curling symptoms in Heilongjiang province of China, via sequencing and assembly of small RNAs. The genome of this geminivirus comprises 2,959 nucleotides (nt) and shares 21.77-54.97% nt sequence identity with other geminiviruses. The coat protein (CP) shares the highest amino acid (aa) sequence identity (23.5%) with that of sesame curly top virus (SeYMV; genus Turncurtovirus), whereas the C1 (Rep) shares the highest aa sequence identity (66.5%) with that of beet severe curly top virus (BSCTV; genus Curtovirus). This geminivirus neighbors the turncurtoviruses in phylogenetic trees based on the full genome sequence or the amino acid sequence of the Rep protein, but it forms a distinct clade in the phylogenetic tree based on the coat protein. Recombination analysis showed that parts of the C1 coding region of this geminivirus were recombined from a curtovirus or turncurtovirus. Based on these results, the name "common bean curly stunt virus" (CBCSV) is proposed for this virus.

Genome-wide signatures of the geographic expansion and breeding of soybean.

Soybean is a leguminous crop that provides oil and protein. Exploring the genomic signatures of soybean evolution is crucial for breeding varieties with improved adaptability to environmental extremes. We analyzed the genome sequences of 2,214 soybeans and proposed a soybean evolutionary route, i.e., the expansion of annual wild soybean (Glycine soja Sieb. & Zucc.) from southern China and its domestication in central China, followed by the expansion and local breeding selection of its landraces (G. max (L.) Merr.). We observed that the genetic introgression in soybean landraces was mostly derived from sympatric rather than allopatric wild populations during the geographic expansion. Soybean expansion and breeding were accompanied by the positive selection of flowering time genes, including GmSPA3c. Our study sheds light on the evolutionary history of soybean and provides valuable genetic resources for its future breeding.

Fine Mapping the Soybean Mosaic Virus Resistance Gene in Soybean Cultivar Heinong 84 and Development of CAPS Markers for Rapid Identification.

Heinong 84 is one of the major soybean varieties growing in Northeast China, and is resistant to the infection of all strains of soybean mosaic virus (SMV) in the region including the most prevalent strain, N3. However, the resistance gene(s) in Heinong 84 and the resistant mechanism are still elusive. In this study, genetic and next-generation sequencing (NGS)-based bulk segregation analysis (BSA) were performed to map the resistance gene using a segregation population from the cross of Heinong 84 and a susceptible cultivar to strain N3, Zhonghuang 13. Results show that the resistance of Heinong 84 is controlled by a dominant gene on chromosome 13. Further analyses suggest that the resistance gene in Heinong 84 is probably an allele of Rsv1. Finally, two pairs of single-nucleotide-polymorphism (SNP)-based primers that are tightly cosegregated with the resistance gene were designed for rapidly identifying resistant progenies in breeding via the cleaved amplified polymorphic sequence (CAPS) assay.

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.

Potential of marker selection to increase prediction accuracy of genomic selection in soybean (Glycine max L.).

Genomic selection is a promising molecular breeding strategy enhancing genetic gain per unit time. The objectives of our study were to (1) explore the prediction accuracy of genomic selection for plant height and yield per plant in soybean [Glycine max (L.) Merr.], (2) discuss the relationship between prediction accuracy and numbers of markers, and (3) evaluate the effect of marker preselection based on different methods on the prediction accuracy. Our study is based on a population of 235 soybean varieties which were evaluated for plant height and yield per plant at multiple locations and genotyped by 5361 single nucleotide polymorphism markers. We applied ridge regression best linear unbiased prediction coupled with fivefold cross-validations and evaluated three strategies of marker preselection. For plant height, marker density and marker preselection procedure impacted prediction accuracy only marginally. In contrast, for grain yield, prediction accuracy based on markers selected with a haplotype block analyses-based approach increased by approximately 4 % compared with random or equidistant marker sampling. Thus, applying marker preselection based on haplotype blocks is an interesting option for a cost-efficient implementation of genomic selection for grain yield in soybean breeding.

Heterosis and combining ability estimates in isoflavone content using different parental soybean accessions: wild soybean, a valuable germplasm for soybean breeding.

Isoflavone, a group of secondary metabolites in soybean, is beneficial to human health. Improving isoflavone content in soybean seeds has become one of the most important breeding objectives. However, the narrow genetic base of soybean cultivars hampered crop improvement. Wild soybean is an extraordinarily important gene pool for soybean breeding. In order to select an optimal germplasm for breeding programs to increase isoflavone concentration, 36 F1 soybean progenies from different parental accessions (cultivars, wild, Semi-wild and Interspecific) with various total isoflavone (TIF) concentration (High, Middle, Low) were analyzed for their isoflavone content. Results showed that male parents, except for Cultivars, showed positive GCA effects. In particular, wild soybean had higher positive GCA effects for TIF concentration. Both MP and BP heterosis value declined in the hybrid in which male parents were wild soybean, semi-wild soybean, interspecific offspring and cultivar in turn. In general, combining ability and heterosis in hybrids which had relative higher TIF concentration level parents showed better performance than those which had lower TIF concentration level parents. These results indicated characteristics of isoflavone content were mainly governed by additive type of gene action, and wild relatives could be utilized for breeding of soybean cultivars with this trait. A promising combination was found as the best potential hybrid for isoflavone content improvement.

Rapid quantitative detection of Vibrio parahaemolyticus in seafood by MPN-PCR.

This study aimed to adopt MPN-PCR (most probable number-polymerase chain reaction) for rapid detection of the quantity of Vibrio parahaemolyticus in seafood. V. parahaemolyticus in seafood could be quantitated by MPN statistics according to PCR products. The sensitivity of MPN-PCR was 100 times higher than that of direct PCR. Of 225 seafood samples from Qingdao, 165 were positive for the presence of V. parahaemolyticus, with an MPN value of >719 per gram, and about 41.5% of samples were positive for tdh gene-possessing cells. Eighty muscle tissues from the 225 seafood samples were investigated by direct PCR and MPN-PCR, but no V. parahaemolyticus was detected. The MPN-PCR test could be completed in less than 16 h from the time of sample preparation. It was rapid, sensitive, and reliable for comprehensive detection and quick quantitative determination of V. parahaemolyticus in seafood and it revealed the potential risk of illness associated with their consumption.

Expression and characterization of a metalloprotease from a Vibrio parahaemolyticus isolate.

The extracellular zinc metalloprotease from Vibrio parahaemolyticus (VPM) is a putative virulence factor for host infection. It is synthesized from the vpm gene of V. parahaemolyticus as a polypeptide of 814 amino acids with an estimated molecular mass of 89,833 Da, containing a zinc metalloprotease HEXXH consensus motif. To investigate the enzymatic properties of V. parahaemolyticus metalloprotease, the mature vpm gene was overexpressed in Escherichia coli, and the recombinant protein (rVPM) was purified by a His-binding metal affinity column (>95% purity). The activity of the recombinant protease produced in E. coli was examined by gelatin activity staining and proteolytic activity assays using gelatin and azocasein as substrates. rVPM showed maximum activity at about 37 degrees C and pH 8. The cytotoxicity against flounder gill cells and fish pathogenicity indicated a potential role in pathogenesis.