Pan-Genome of Wild and Cultivated Soybeans.

Soybean is one of the most important vegetable oil and protein feed crops. To capture the entire genomic diversity, it is needed to construct a complete high-quality pan-genome from diverse soybean accessions. In this study, we performed individual de novo genome assemblies for 26 representative soybeans that were selected from 2,898 deeply sequenced accessions. Using these assembled genomes together with three previously reported genomes, we constructed a graph-based genome and performed pan-genome analysis, which identified numerous genetic variations that cannot be detected by direct mapping of short sequence reads onto a single reference genome. The structural variations from the 2,898 accessions that were genotyped based on the graph-based genome and the RNA sequencing (RNA-seq) data from the representative 26 accessions helped to link genetic variations to candidate genes that are responsible for important traits. This pan-genome resource will promote evolutionary and functional genomics studies in soybean.

Integrative omics analysis elucidates the genetic basis underlying seed weight and oil content in soybean.

Synergistic optimization of key agronomic traits by traditional breeding has dramatically enhanced crop productivity in the past decades. However, the genetic basis underlying coordinated regulation of yield- and quality-related traits remains poorly understood. Here, we dissected the genetic architectures of seed weight and oil content by combining genome-wide association studies (GWAS) and transcriptome-wide association studies (TWAS) using 421 soybean (Glycine max) accessions. We identified 26 and 33 genetic loci significantly associated with seed weight and oil content by GWAS, respectively, and detected 5,276 expression quantitative trait loci (eQTLs) regulating expression of 3,347 genes based on population transcriptomes. Interestingly, a gene module (IC79), regulated by two eQTL hotspots, exhibited significant correlation with both seed weigh and oil content. Twenty-two candidate causal genes for seed traits were further prioritized by TWAS, including Regulator of Weight and Oil of Seed 1 (GmRWOS1), which encodes a sodium pump protein. GmRWOS1 was verified to pleiotropically regulate seed weight and oil content by gene knockout and overexpression. Notably, allelic variations of GmRWOS1 were strongly selected during domestication of soybean. This study uncovers the genetic basis and network underlying regulation of seed weight and oil content in soybean and provides a valuable resource for improving soybean yield and quality by molecular breeding.

The Myb73-GDPD2-GA2ox1 transcriptional regulatory module confers phosphate deficiency tolerance in soybean.

Phosphorus is indispensable in agricultural production. An increasing food supply requires more efficient use of phosphate due to limited phosphate resources. However, how crops regulate phosphate efficiency remains largely unknown. Here, we identified a major quantitative trait locus, qPE19, that controls 7 low-phosphate (LP)-related traits in soybean (Glycine max) through linkage mapping and genome-wide association studies. We identified the gene responsible for qPE19 as GLYCEROPHOSPHORYL DIESTER PHOSPHODIESTERASE2 (GmGDPD2), and haplotype 5 represents the optimal allele favoring LP tolerance. Overexpression of GmGDPD2 significantly affects hormone signaling and improves root architecture, phosphate efficiency and yield-related traits; conversely, CRISPR/Cas9-edited plants show decreases in these traits. GmMyb73 negatively regulates GmGDPD2 by directly binding to its promoter; thus, GmMyb73 negatively regulates LP tolerance. GmGDPD2 physically interacts with GA 2-oxidase 1 (GmGA2ox1) in the plasma membrane, and overexpressing GmGA2ox1 enhances LP-associated traits, similar to GmGDPD2 overexpression. Analysis of double mutants for GmGDPD2 and GmGA2ox1 demonstrated that GmGDPD2 regulates LP tolerance likely by influencing auxin and gibberellin dose-associated cell division in the root. These results reveal a regulatory module that plays a major role in regulating LP tolerance in soybeans and is expected to be utilized to develop phosphate-efficient varieties to enhance soybean production, particularly in phosphate-deficient soils.

Mitogen-activated protein kinases MPK3 and MPK6 phosphorylate receptor-like cytoplasmic kinase CDL1 to regulate soybean basal immunity.

Soybean cyst nematode (SCN; Heterodera glycines Ichinohe), one of the most devastating soybean (Glycine max) pathogens, causes significant yield loss in soybean production. Nematode infection triggers plant defense responses; however, the components involved in the upstream signaling cascade remain largely unknown. In this study, we established that a mitogen-activated protein kinase (MAPK) signaling module, activated by nematode infection or wounding, is crucial for soybeans to establish SCN resistance. GmMPK3 and GmMPK6 directly interact with CDG1-LIKE1 (GmCDL1), a member of the receptor-like cytoplasmic kinase (RLCK) subfamily VII. These kinases phosphorylate GmCDL1 at Thr-372 to prevent its proteasome-mediated degradation. Functional analysis demonstrated that GmCDL1 positively regulates immune responses and promotes SCN resistance in soybeans. GmMPK3-mediated and GmMPK6-mediated phosphorylation of GmCDL1 enhances GmMPK3 and GmMPK6 activation and soybean disease resistance, representing a positive feedback mechanism. Additionally, 2 L-type lectin receptor kinases, GmLecRK02g and GmLecRK08g, associate with GmCDL1 to initiate downstream immune signaling. Notably, our study also unveils the potential involvement of GmLecRKs and GmCDL1 in countering other soybean pathogens beyond nematodes. Taken together, our findings reveal the pivotal role of the GmLecRKs-GmCDL1-MAPK regulatory module in triggering soybean basal immune responses.

Enhanced weathering in the US Corn Belt delivers carbon removal with agronomic benefits.

Terrestrial enhanced weathering (EW) of silicate rocks, such as crushed basalt, on farmlands is a promising scalable atmospheric carbon dioxide removal (CDR) strategy that urgently requires performance assessment with commercial farming practices. We report findings from a large-scale replicated EW field trial across a typical maize-soybean rotation on an experimental farm in the heart of the United Sates Corn Belt over 4 y (2016 to 2020). We show an average combined loss of major cations (Ca2+ and Mg2+) from crushed basalt applied each fall over 4 y (50 t ha-1 y-1) gave a conservative time-integrated cumulative CDR potential of 10.5 ± 3.8 t CO2 ha-1. Maize and soybean yields increased significantly (P < 0.05) by 12 to 16% with EW following improved soil fertility, decreased soil acidification, and upregulation of root nutrient transport genes. Yield enhancements with EW were achieved with significantly (P < 0.05) increased key micro- and macronutrient concentrations (including potassium, magnesium, manganese, phosphorus, and zinc), thus improving or maintaining crop nutritional status. We observed no significant increase in the content of trace metals in grains of maize or soybean or soil exchangeable pools relative to controls. Our findings suggest that widespread adoption of EW across farming sectors has the potential to contribute significantly to net-zero greenhouse gas emissions goals while simultaneously improving food and soil security.

H2O2-dependent oxidation of the transcription factor GmNTL1 promotes salt tolerance in soybean.

Reactive oxygen species (ROS) play an essential role in plant growth and responses to environmental stresses. Plant cells sense and transduce ROS signaling directly via hydrogen peroxide (H2O2)-mediated posttranslational modifications (PTMs) on protein cysteine residues. Here, we show that the H2O2-mediated cysteine oxidation of NAC WITH TRANS-MEMBRANE MOTIF1-LIKE 1 (GmNTL1) in soybean (Glycine max) during salt stress promotes its release from the endoplasmic reticulum (ER) membrane and translocation to the nucleus. We further show that an oxidative posttranslational modification on GmNTL1 residue Cys-247 steers downstream amplification of ROS production by binding to and activating the promoters of RESPIRATORY BURST OXIDASE HOMOLOG B (GmRbohB) genes, thereby creating a feed-forward loop to fine-tune GmNTL1 activity. In addition, oxidation of GmNTL1 Cys-247 directly promotes the expression of CATION H+ EXCHANGER 1 (GmCHX1)/SALT TOLERANCE-ASSOCIATED GENE ON CHROMOSOME 3 (GmSALT3) and Na+/H+ Antiporter 1 (GmNHX1). Accordingly, transgenic overexpression of GmNTL1 in soybean increases the H2O2 levels and K+/Na+ ratio in the cell, promotes salt tolerance, and increases yield under salt stress, while an RNA interference-mediated knockdown of GmNTL1 elicits the opposite effects. Our results reveal that the salt-induced oxidation of GmNTL1 promotes its relocation and transcriptional activity through an H2O2-mediated posttranslational modification on cysteine that improves resilience of soybean against salt stress.

GmNAC039 and GmNAC018 activate the expression of cysteine protease genes to promote soybean nodule senescence.

Root nodules are major sources of nitrogen for soybean (Glycine max (L.) Merr.) growth, development, production, and seed quality. Symbiotic nitrogen fixation is time-limited, as the root nodule senesces during the reproductive stage of plant development, specifically during seed development. Nodule senescence is characterized by the induction of senescence-related genes, such as papain-like cysteine proteases (CYPs), which ultimately leads to the degradation of both bacteroids and plant cells. However, how nodule senescence-related genes are activated in soybean is unknown. Here, we identified 2 paralogous NAC transcription factors, GmNAC039 and GmNAC018, as master regulators of nodule senescence. Overexpression of either gene induced soybean nodule senescence with increased cell death as detected using a TUNEL assay, whereas their knockout delayed senescence and increased nitrogenase activity. Transcriptome analysis and nCUT&Tag-qPCR assays revealed that GmNAC039 directly binds to the core motif CAC(A)A and activates the expression of 4 GmCYP genes (GmCYP35, GmCYP37, GmCYP39, and GmCYP45). Similar to GmNAC039 and GmNAC018, overexpression or knockout of GmCYP genes in nodules resulted in precocious or delayed senescence, respectively. These data provide essential insights into the regulatory mechanisms of nodule senescence, in which GmNAC039 and GmNAC018 directly activate the expression of GmCYP genes to promote nodule senescence.

SoyMD: a platform combining multi-omics data with various tools for soybean research and breeding.

Advanced multi-omics technologies offer much information that can uncover the regulatory mechanisms from genotype to phenotype. In soybean, numerous multi-omics databases have been published. Although they cover multiple omics, there are still limitations when it comes to the types and scales of omics datasets and analysis methods utilized. This study aims to address these limitations by collecting and integrating a comprehensive set of multi-omics datasets. This includes 38 genomes, transcriptomes from 435 tissue samples, 125 phenotypes from 6686 accessions, epigenome data involving histone modification, transcription factor binding, chromosomal accessibility and chromosomal interaction, as well as genetic variation data from 24 501 soybean accessions. Then, common analysis pipelines and statistical methods were applied to mine information from these multi-omics datasets, resulting in the successful establishment of a user-friendly multi-omics database called SoyMD (https://yanglab.hzau.edu.cn/SoyMD/#/). SoyMD provides researchers with efficient query options and analysis tools, allowing them to swiftly access relevant omics information and conduct comprehensive multi-omics data analyses. Another notable feature of SoyMD is its capability to facilitate the analysis of candidate genes, as demonstrated in the case study on seed oil content. This highlights the immense potential of SoyMD in soybean genetic breeding and functional genomics research.

Dynamics of cis-regulatory sequences and transcriptional divergence of duplicated genes in soybean.

Transcriptional divergence of duplicated genes after whole genome duplication (WGD) has been described in many plant lineages and is often associated with subgenome dominance, a genome-wide mechanism. However, it is unknown what underlies the transcriptional divergence of duplicated genes in polyploid species that lack subgenome dominance. Soybean is a paleotetraploid with a WGD that occurred 5 to 13 Mya. Approximately 50% of the duplicated genes retained from this WGD exhibit transcriptional divergence. We developed accessible chromatin region (ACR) datasets from leaf, flower, and seed tissues using MNase-hypersensitivity sequencing. We validated enhancer function of several ACRs associated with known genes using CRISPR/Cas9-mediated genome editing. The ACR datasets were used to examine and correlate the transcriptional patterns of 17,111 pairs of duplicated genes in different tissues. We demonstrate that ACR dynamics are correlated with divergence of both expression level and tissue specificity of individual gene pairs. Gain or loss of flanking ACRs and mutation of cis-regulatory elements (CREs) within the ACRs can change the balance of the expression level and/or tissue specificity of the duplicated genes. Analysis of DNA sequences associated with ACRs revealed that the extensive sequence rearrangement after the WGD reshaped the CRE landscape, which appears to play a key role in the transcriptional divergence of duplicated genes in soybean. This may represent a general mechanism for transcriptional divergence of duplicated genes in polyploids that lack subgenome dominance.

Temporal and spatial resolution of distal protein motions that activate hydrogen tunneling in soybean lipoxygenase.

The enzyme soybean lipoxygenase (SLO) provides a prototype for deep tunneling mechanisms in hydrogen transfer catalysis. This work combines room temperature X-ray studies with extended hydrogen-deuterium exchange experiments to define a catalytically-linked, radiating cone of aliphatic side chains that connects an active site iron center of SLO to the protein-solvent interface. Employing eight variants of SLO that have been appended with a fluorescent probe at the identified surface loop, nanosecond fluorescence Stokes shifts have been measured. We report a remarkable identity of the energies of activation (Ea) for the Stokes shifts decay rates and the millisecond C-H bond cleavage step that is restricted to side chain mutants within an identified thermal network. These findings implicate a direct coupling of distal protein motions surrounding the exposed fluorescent probe to active site motions controlling catalysis. While the role of dynamics in enzyme function has been predominantly attributed to a distributed protein conformational landscape, the presented data implicate a thermally initiated, cooperative protein reorganization that occurs on a timescale faster than nanosecond and represents the enthalpic barrier to the reaction of SLO.

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.

Pan-centromere reveals widespread centromere repositioning of soybean genomes.

Centromere repositioning refers to a de novo centromere formation at another chromosomal position without sequence rearrangement. This phenomenon was frequently encountered in both mammalian and plant species and has been implicated in genome evolution and speciation. To understand the dynamic of centromeres on soybean genome, we performed the pan-centromere analysis using CENH3-ChIP-seq data from 27 soybean accessions, including 3 wild soybeans, 9 landraces, and 15 cultivars. Building upon the previous discovery of three centromere satellites in soybean, we have identified two additional centromere satellites that specifically associate with chromosome 1. These satellites reveal significant rearrangements in the centromere structures of chromosome 1 across different accessions, consequently impacting the localization of CENH3. By comparative analysis, we reported a high frequency of centromere repositioning on 14 out of 20 chromosomes. Most newly emerging centromeres formed in close proximity to the native centromeres and some newly emerging centromeres were apparently shared in distantly related accessions, suggesting their emergence is independent. Furthermore, we crossed two accessions with mismatched centromeres to investigate how centromere positions would be influenced in hybrid genetic backgrounds. We found that a significant proportion of centromeres in the S9 generation undergo changes in size and position compared to their parental counterparts. Centromeres preferred to locate at satellites to maintain a stable state, highlighting a significant role of centromere satellites in centromere organization. Taken together, these results revealed extensive centromere repositioning in soybean genome and highlighted how important centromere satellites are in constraining centromere positions and supporting centromere function.

GmEID1 modulates light signaling through the Evening Complex to control flowering time and yield in soybean.

Soybean (Glycine max) morphogenesis and flowering time are accurately regulated by photoperiod, which determine the yield potential and limit soybean cultivars to a narrow latitudinal range. The E3 and E4 genes, which encode phytochrome A photoreceptors in soybean, promote the expression of the legume-specific flowering repressor E1 to delay floral transition under long-day (LD) conditions. However, the underlying molecular mechanism remains unclear. Here, we show that the diurnal expression pattern of GmEID1 is opposite to that of E1 and targeted mutations in the GmEID1 gene delay soybean flowering regardless of daylength. GmEID1 interacts with J, a key component of circadian Evening Complex (EC), to inhibit E1 transcription. Photoactivated E3/E4 interacts with GmEID1 to inhibit GmEID1-J interaction, promoting J degradation resulting in a negative correlation between daylength and the level of J protein. Notably, targeted mutations in GmEID1 improved soybean adaptability by enhancing yield per plant up to 55.3% compared to WT in field trials performed in a broad latitudinal span of more than 24°. Together, this study reveals a unique mechanism in which E3/E4-GmEID1-EC module controls flowering time and provides an effective strategy to improve soybean adaptability and production for molecular breeding.

Effects of demographic history on recombination hotspots in soybean.

Recombination is the primary mechanism underlying genetic improvement in populations and allows plant breeders to create new allelic combinations for agronomic improvement. Soybean [Glycine max (L.) Merr.] has gone through multiple genetic bottlenecks that have significantly affected its genetic diversity, linkage disequilibrium, and altered allele frequencies. To investigate the impact of genetic bottlenecks on recombination hotspots in soybeans, historical recombination was studied in three soybean populations. The populations were wild soybean [Glycine soja (Sieb. and Zucc.)], landraces, and North American elite soybean cultivars that have been genotyped with the SoySNP50K BeadChip. While each population after a genetic bottleneck had an increased average haplotype block size, they did not have a significant difference in the number of hotspots between each population. Instead, the increase in observed haplotype block size is likely due to an elimination of individuals that contained historical recombination at hotspots which decreased the observed rate of recombination for the hotspot after each genetic bottleneck. Conversely, heterochromatic DNA which has an increased haplotype block size compared to euchromatic DNA had a significantly different number of hotspots but not a significant difference in the average hotspot recombination rate. Previously identified genomic motifs associated with hotspots were also associated with hotspots found in the historical populations suggesting a common mechanism. This characterization of historical recombination hotspots in soybeans provides further insights into the effect genetic bottlenecks and selection have on recombination hotspots.

GmNLP1 and GmNLP4 activate nitrate-induced CLE peptides NIC1a/b to mediate nitrate-regulated root nodulation.

Symbiotic nitrogen fixation is an energy-intensive process, to maintain the balance between growth and nitrogen fixation, high concentrations of nitrate inhibit root nodulation. However, the precise mechanism underlying the nitrate inhibition of nodulation in soybean remains elusive. In this study, CRISPR-Cas9-mediated knockout of GmNLP1 and GmNLP4 unveiled a notable nitrate-tolerant nodulation phenotype. GmNLP1b and GmNLP4a play a significant role in the nitrate-triggered inhibition of nodulation, as the expression of nitrate-responsive genes was largely suppressed in Gmnlp1b and Gmnlp4a mutants. Furthermore, we demonstrated that GmNLP1b and GmNLP4a can bind to the promoters of GmNIC1a and GmNIC1b and activate their expression. Manipulations targeting GmNIC1a and GmNIC1b through knockdown or overexpression strategies resulted in either increased or decreased nodule number in response to nitrate. Additionally, transgenic roots that constitutively express GmNIC1a or GmNIC1b rely on both NARK and hydroxyproline O-arabinosyltransferase RDN1 to prevent the inhibitory effects imposed by nitrate on nodulation. In conclusion, this study highlights the crucial role of the GmNLP1/4-GmNIC1a/b module in mediating high nitrate-induced inhibition of nodulation.

Light recovery after maize harvesting promotes soybean flowering in a maize-soybean relay strip intercropping system.

Moving from sole cropping to intercropping is a transformative change in agriculture, contributing to yield. Soybeans adapt to light conditions in intercropping by adjusting the onset of reproduction and the inflorescence architecture to optimize reproductive success. Maize-soybean strip intercropping (MS), maize-soybean relay strip intercropping (IS), and sole soybean (SS) systems are typical soybean planting systems with significant differences in light environments during growth periods. To elucidate the effect of changes in the light environment on soybean flowering processes and provide a theoretical basis for selecting suitable varieties in various planting systems to improve yields, field experiments combining planting systems (IS, MS, and SS) and soybean varieties (GQ8, GX7, ND25, and NN996) were conducted in 2021 and 2022. Results showed that growth recovery in the IS resulted in a balance in the expression of TERMINAL FLOWER 1 (TFL1) and FLOWERING LOCUS T (FT) in the meristematic tissues of soybeans, which promoted the formation of new branches or flowers. IS prolonged the flowering time (2-7 days) and increased the number of forming flowers compared with SS (93.0 and 169%) and MS (67.3 and 103.3%) at the later soybean flowering stage. The higher carbon and nitrogen content in the middle and bottom canopies of soybean contributed to decreased flower abscission by 26.7 and 30.2%, respectively, compared with SS. Canopy light environment recovery promoted branch and flower formation and transformation of flowers into pods with lower flower-pod abscission, which contributed to elevating soybean yields in late-maturing and multibranching varieties (ND25) in IS.

The acyl-acyl carrier protein thioesterases GmFATA1 and GmFATA2 are essential for fatty acid accumulation and growth in soybean.

Acyl-acyl carrier protein (ACP) thioesterases (FAT) hydrolyze acyl-ACP complexes to release FA in plastids, which ultimately affects FA biosynthesis and profiles. Soybean GmFATA1 and GmFATA2 are homoeologous genes encoding oleoyl-ACP thioesterases whose role in seed oil accumulation and plant growth has not been defined. Using CRISPR/Cas9 gene editing mutation of Gmfata1 or 2 led to reduced leaf FA content and growth defect at the early seedling stage. In contrast, no homozygous double mutants were obtained. Combined this indicates that GmFATA1 and GmFATA2 display overlapping, but not complete functional redundancy. Combined transcriptomic and lipidomic analysis revealed a large number of genes involved in FA synthesis and FA chain elongation are expressed at reduced level in the Gmfata1 mutant, accompanied by a lower triacylglycerol abundance at the early seedling stage. Further analysis showed that the Gmfata1 or 2 mutants had increased composition of the beneficial FA, oleic acid. The growth defect of Gmfata1 could be at least partially attributed to reduced acetyl-CoA carboxylase activity, reduced abundance of five unsaturated monogalactosyldiacylglycerol lipids, and altered chloroplast morphology. On the other hand, overexpression of GmFATA in soybean led to significant increases in leaf FA content by 5.7%, vegetative growth, and seed yield by 26.9%, and seed FA content by 23.2%. Thus, overexpression of GmFATA is an effective strategy to enhance soybean oil content and yield.

Mutations in the genes responsible for the synthesis of furan fatty acids resolve the light-induced off-odor in soybean oil.

Soybean oil is the second most produced edible vegetable oil and is used for many edible and industrial materials. Unfortunately, it has the disadvantage of 'reversion flavor' under photooxidative conditions, which produces an off-odor and decreases the quality of edible oil. Reversion flavor and off-odor are caused by minor fatty acids in the triacylglycerol of soybean oil known as furan fatty acids, which produce 3-methyl-2,4-nonanedione (3-MND) upon photooxidation. As a solution to this problem, a reduction in furan fatty acids leads to a decrease in 3-MND, resulting in a reduction in the off-odor induced by light exposure. However, there are no reports on the genes related to the biosynthesis of furan fatty acids in soybean oil. In this study, four mutant lines showing low or no furan fatty acid levels in soybean seeds were isolated from a soybean mutant library. Positional cloning experiments and homology search analysis identified two genes responsible for furan fatty acid biosynthesis in soybean: Glyma.20G201400 and Glyma.04G054100. Ectopic expression of both genes produced furan fatty acids in transgenic soybean hairy roots. The structure of these genes is different from that of the furan fatty acid biosynthetic genes in photosynthetic bacteria. Homologs of these two group of genes are widely conserved in the plant kingdom. The purified oil from the furan fatty acid mutant lines had lower amounts of 3-MND and reduced off-odor after light exposure, compared with oil from the wild-type.

Degenerate oligonucleotide primer MIG-seq: an effective PCR-based method for high-throughput genotyping.

Next-generation sequencing (NGS) library construction often involves using restriction enzymes to decrease genome complexity, enabling versatile polymorphism detection in plants. However, plant leaves frequently contain impurities, such as polyphenols, necessitating DNA purification before enzymatic reactions. To overcome this problem, we developed a PCR-based method for expeditious NGS library preparation, offering flexibility in number of detected polymorphisms. By substituting a segment of the simple sequence repeat sequence in the MIG-seq primer set (MIG-seq being a PCR method enabling library construction with low-quality DNA) with degenerate oligonucleotides, we introduced variability in detectable polymorphisms across various crops. This innovation, named degenerate oligonucleotide primer MIG-seq (dpMIG-seq), enabled a streamlined protocol for constructing dpMIG-seq libraries from unpurified DNA, which was implemented stably in several crop species, including fruit trees. Furthermore, dpMIG-seq facilitated efficient lineage selection in wheat and enabled linkage map construction and quantitative trait loci analysis in tomato, rice, and soybean without necessitating DNA concentration adjustments. These findings underscore the potential of the dpMIG-seq protocol for advancing genetic analyses across diverse plant species.

SoyDNGP: a web-accessible deep learning framework for genomic prediction in soybean breeding.

Soybean is a globally significant crop, playing a vital role in human nutrition and agriculture. Its complex genetic structure and wide trait variation, however, pose challenges for breeders and researchers aiming to optimize its yield and quality. Addressing this biological complexity requires innovative and accurate tools for trait prediction. In response to this challenge, we have developed SoyDNGP, a deep learning-based model that offers significant advancements in the field of soybean trait prediction. Compared to existing methods, such as DeepGS and DNNGP, SoyDNGP boasts a distinct advantage due to its minimal increase in parameter volume and superior predictive accuracy. Through rigorous performance comparison, including prediction accuracy and model complexity, SoyDNGP represents improved performance to its counterparts. Furthermore, it effectively predicted complex traits with remarkable precision, demonstrating robust performance across different sample sizes and trait complexities. We also tested the versatility of SoyDNGP across multiple crop species, including cotton, maize, rice and tomato. Our results showed its consistent and comparable performance, emphasizing SoyDNGP's potential as a versatile tool for genomic prediction across a broad range of crops. To enhance its accessibility to users without extensive programming experience, we designed a user-friendly web server, available at http://xtlab.hzau.edu.cn/SoyDNGP. The server provides two features: 'Trait Lookup', offering users the ability to access pre-existing trait predictions for over 500 soybean accessions, and 'Trait Prediction', allowing for the upload of VCF files for trait estimation. By providing a high-performing, accessible tool for trait prediction, SoyDNGP opens up new possibilities in the quest for optimized soybean breeding.