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.

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.

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.

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.

Multi-omics-driven advances in the understanding of triacylglycerol biosynthesis in oil seeds.

Vegetable oils are rich sources of polyunsaturated fatty acids and energy as well as valuable sources of human food, animal feed, and bioenergy. Triacylglycerols, which are comprised of three fatty acids attached to a glycerol backbone, are the main component of vegetable oils. Here, we review the development and application of multiple-level omics in major oilseeds and emphasize the progress in the analysis of the biological roles of key genes underlying seed oil content and quality in major oilseeds. Finally, we discuss future research directions in functional genomics research based on current omics and oil metabolic engineering strategies that aim to enhance seed oil content and quality, and specific fatty acids components according to either human health needs or industrial requirements.

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.

MutL homolog 1 participates in interference-sensitive meiotic crossover formation in soybean.

MutL homolog 1 (MLH1), a member of the MutL-homolog family, is required for normal recombination in most organisms. However, its role in soybean (Glycine max) remains unclear to date. Here, we characterized the Glycine max female and male sterility 1 (Gmfms1) mutation that reduces pollen grain viability and increases embryo sac abortion in soybean. Map-based cloning revealed that the causal gene of Gmfms1 is Glycine max MutL homolog 1 (GmMLH1), and CRISPR/Cas9 knockout approach further validated that disruption of GmMLH1 confers the female-male sterility phenotype in soybean. Loss of GmMLH1 function disrupted bivalent formation, leading to univalent mis-segregation during meiosis and ultimately to female-male sterility. The Gmmlh1 mutant showed about a 78.16% decrease in meiotic crossover frequency compared to the wild type. The residual chiasmata followed a Poisson distribution, suggesting that interference-sensitive crossover formation was affected in the Gmmlh1 mutant. Furthermore, GmMLH1 could interact with GmMLH3A and GmMLH3B both in vivo and in vitro. Overall, our work demonstrates that GmMLH1 participates in interference-sensitive crossover formation in soybean, and provides additional information about the conserved functions of MLH1 across plant species.

Transcription factor GmAlfin09 regulates endoplasmic reticulum stress in soybean via peroxidase GmPRDX6.

Soybean (Glycine max [L.] Merr.) is a valuable oil crop but is also highly susceptible to environmental stress. Thus, developing approaches to enhance soybean stress resistance is vital to soybean yield improvement. In previous studies, transcription factor Alfin has been shown to serve as an epigenetic regulator of plant growth and development. However, no studies on Alfin have yet been reported in soybean. In this study, the endoplasmic reticulum (ER) stress- and reactive oxygen species (ROS)-related GmAlfin09 was identified. Screening of genes co-expressed with GmAlfin09 unexpectedly led to the identification of soybean peroxidase 6 (GmPRDX6). Further analyses revealed that both GmAlfin09 and GmPRDX6 were responsive to ER stress, with GmPRDX6 localizing to the ER under stress. Promoter binding experiments confirmed the ability of GmAlfin09 to bind to the GmPRDX6 promoter directly. When GmAlfin09 and GmPRDX6 were overexpressed in soybean, enhanced ER stress resistance and decreased ROS levels were observed. Together, these findings suggest that GmAlfin09 promotes the upregulation of GmPRDX6, and GmPRDX6 subsequently localizes to the ER, reduces ROS levels, promotes ER homeostasis, and ensures the normal growth of soybean even under ER stress. This study highlights a vital target gene for future molecular breeding of stress-resistant soybean lines.

Soybean ethylene response factors GmENS1 and GmENS2 promote nodule senescence.

The final phase in root nodule development is nodule senescence. The mechanism underlying the initiation of nodule senescence requires further elucidation. Here, we investigated the intrinsic signals governing soybean (Glycine max L. Merr.) nodule senescence, uncovering ethylene as a key signal in this intricate mechanism. Two AP2/ERF transcription factor genes, GmENS1 and GmENS2 (Ethylene-responsive transcription factors required for Nodule Senescence), exhibit heightened expression levels in both aged nodules and nodules treated with ethylene. Overexpression of either GmENS1 or GmENS2 accelerated senescence in soybean nodules, whereas the knockout or knockdown of both genes delayed senescence and enhanced nitrogenase activity. Furthermore, our findings indicated that GmENS1 and GmENS2 directly bind to the promoters of GmNAC039, GmNAC018, and GmNAC030, encoding three NAC transcription factors essential for activating soybean nodule senescence. Notably, the nodule senescence process mediated by GmENS1 or GmENS2 overexpression was suppressed in the soybean nac039/018/030 triple mutant compared with the wild-type control. These data indicate GmENS1 and GmENS2 as pivotal transcription factors mediating ethylene-induced nodule senescence through the direct activation of GmNAC039/GmNAC018/GmNAC030 expression in soybean.

Novel and multifaceted regulations of photoperiodic flowering by phytochrome A in soybean.

Photoperiod is an important environmental cue. Plants can distinguish the seasons and flower at the right time through sensing the photoperiod. Soybean is a sensitive short-day crop, and the timing of flowering varies greatly at different latitudes, thus affecting yields. Soybean cultivars in high latitudes adapt to the long day by the impairment of two phytochrome genes, PHYA3 and PHYA2, and the legume-specific flowering suppressor, E1. However, the regulating mechanism underlying phyA and E1 in soybean remains largely unknown. Here, we classified the regulation of the E1 family by phyA2 and phyA3 at the transcriptional and posttranscriptional levels, revealing that phyA2 and phyA3 regulate E1 by directly binding to LUX proteins, the critical component of the evening complex, to regulate the stability of LUX proteins. In addition, phyA2 and phyA3 can also directly associate with E1 and its homologs to stabilize the E1 proteins. Therefore, phyA homologs control the core flowering suppressor E1 at both the transcriptional and posttranscriptional levels, to double ensure the E1 activity. Thus, our results disclose a photoperiod flowering mechanism in plants by which the phytochrome A regulates LUX and E1 activity.

A high-performance computational workflow to accelerate GATK SNP detection across a 25-genome dataset.

BACKGROUND: Single-nucleotide polymorphisms (SNPs) are the most widely used form of molecular genetic variation studies. As reference genomes and resequencing data sets expand exponentially, tools must be in place to call SNPs at a similar pace. The genome analysis toolkit (GATK) is one of the most widely used SNP calling software tools publicly available, but unfortunately, high-performance computing versions of this tool have yet to become widely available and affordable. RESULTS: Here we report an open-source high-performance computing genome variant calling workflow (HPC-GVCW) for GATK that can run on multiple computing platforms from supercomputers to desktop machines. We benchmarked HPC-GVCW on multiple crop species for performance and accuracy with comparable results with previously published reports (using GATK alone). Finally, we used HPC-GVCW in production mode to call SNPs on a "subpopulation aware" 16-genome rice reference panel with ~ 3000 resequenced rice accessions. The entire process took ~ 16 weeks and resulted in the identification of an average of 27.3 M SNPs/genome and the discovery of ~ 2.3 million novel SNPs that were not present in the flagship reference genome for rice (i.e., IRGSP RefSeq). CONCLUSIONS: This study developed an open-source pipeline (HPC-GVCW) to run GATK on HPC platforms, which significantly improved the speed at which SNPs can be called. The workflow is widely applicable as demonstrated successfully for four major crop species with genomes ranging in size from 400 Mb to 2.4 Gb. Using HPC-GVCW in production mode to call SNPs on a 25 multi-crop-reference genome data set produced over 1.1 billion SNPs that were publicly released for functional and breeding studies. For rice, many novel SNPs were identified and were found to reside within genes and open chromatin regions that are predicted to have functional consequences. Combined, our results demonstrate the usefulness of combining a high-performance SNP calling architecture solution with a subpopulation-aware reference genome panel for rapid SNP discovery and public deployment.

Chromatograms and Mass Spectra of High-Mannose and Paucimannose N-Glycans for Rapid Isomeric Identifications.

N-Linked glycosylation is one of the most essential post-translational modifications of proteins. However, N-glycan structural determination remains challenging because of the small differences in structures between isomers. In this study, we constructed a database containing collision-induced dissociation MSn mass spectra and chromatograms of high-performance liquid chromatography for the rapid identification of high-mannose and paucimannose N-glycan isomers. These N-glycans include isomers by breaking of arbitrary numbers of glycosidic bonds at arbitrary positions of canonical Man9GlcNAc2 N-glycans. In addition, some GlcMannGlcNAc2 N-glycan isomers were included in the database. This database is particularly useful for the identification of the N-glycans not in conventional N-glycan standards. This study demonstrated the application of the database to structural assignment for high-mannose N-glycans extracted from bovine whey proteins, soybean proteins, human mammary epithelial cells, and human breast carcinoma cells. We found many N-glycans that are not expected to be generated by conventional biosynthetic pathways of multicellular eukaryotes.

An NBS-LRR protein in the Rpp1 locus negates the dominance of Rpp1-mediated resistance against Phakopsora pachyrhizi in soybean.

The soybean Rpp1 locus confers resistance to Phakopsora pachyrhizi, causal agent of rust, and resistance is usually dominant over susceptibility. However, dominance of Rpp1-mediated resistance is lost when a resistant genotype (Rpp1 or Rpp1b) is crossed with susceptible line TMG06_0011, and the mechanism of this dominant susceptibility (DS) is unknown. Sequencing the Rpp1 region reveals that the TMG06_0011 Rpp1 locus has a single nucleotide-binding site leucine-rich repeat (NBS-LRR) gene (DS-R), whereas resistant PI 594760B (Rpp1b) is similar to PI 200492 (Rpp1) and has three NBS-LRR resistance gene candidates. Evidence that DS-R is the cause of DS was reflected in virus-induced gene silencing of DS-R in Rpp1b/DS-R or Rpp1/DS-R heterozygous plants with resistance partially restored. In heterozygous Rpp1b/DS-R plants, expression of Rpp1b candidate genes was not significantly altered, indicating no effect of DS-R on transcription. Physical interaction of the DS-R protein with candidate Rpp1b resistance proteins was supported by yeast two-hybrid studies and in silico modeling. Thus, we conclude that suppression of resistance most likely does not occur at the transcript level, but instead probably at the protein level, possibly with Rpp1 function inhibited by binding to the DS-R protein. The DS-R gene was found in other soybean lines, with an estimated allele frequency of 6% in a diverse population, and also found in wild soybean (Glycine soja). The identification of a dominant susceptible NBS-LRR gene provides insight into the behavior of NBS-LRR proteins and serves as a reminder to breeders that the dominance of an R gene can be influenced by a susceptibility allele.

A natural non-synonymous single nucleotide polymorphism in GmbHLH113 negates its inhibitory effect on root hair elongation in soybean.

Root hair length (RHL) is an important character that affects nutrient acquisition in plants. The regulatory network in soybean controlling RHL is yet to be fully understood. In this study, we identified a quantitative trait locus (QTL) regulating RHL. One candidate causal gene in this QTL (GmbHLH113), preferentially expressed in root hairs, was annotated as encoding a basic helix-loop-helix transcription factor. In wild soybeans, the allelic type of GmbHLH113 with a glycine in the 13th residue, which was associated with a reduction in RHL, was shown to localize in the nucleus and activate gene transcription. Another allelic type with a single nucleotide polymorphism that resulted in a glutamate in the 13th residue is fixed in cultivated soybeans, and it lost the ability to localize to the nucleus or negatively regulate RHL. The ectopic expression of GmbHLH113 from W05 in Arabidopsis root hairs resulted in shorter RHL and reduced phosphorus (P) accumulation in shoots. Hence, a loss-of-function allele in cultivated soybeans might have been selected during domestication due to its association with a longer RHL and improved nutrient acquisition.

The Soybean Rpp3 Gene Encodes a TIR-NBS-LRR Protein that Confers Resistance to Phakopsora pachyrhizi.

Soybean rust is an economically significant disease caused by the fungus Phakopsora pachyrhizi that negatively impacts soybean (Glycine max [L.] Merr.) production throughout the world. Susceptible plants infected by P. pachyrhizi develop tan-colored lesions on the leaf surface that give rise to funnel-shaped uredinia as the disease progresses. While most soybean germplasm is susceptible, seven genetic loci (Rpp1 to Rpp7) that provide race-specific resistance to P. pachyrhizi (Rpp) have been identified. Rpp3 was first discovered and characterized in the soybean accession PI 462312 (Ankur), and it was also determined to be one of two Rpp genes present in PI 506764 (Hyuuga). Genetic crosses with PI 506764 were later used to fine-map the Rpp3 locus to a 371-kb region on chromosome 6. The corresponding region in the susceptible Williams 82 (Wm82) reference genome contains several homologous nucleotide binding site-leucine rich repeat (NBS-LRR) genes. To identify Rpp3, we designed oligonucleotide primers to amplify Rpp3 candidate (Rpp3C) NBS-LRR genes at this locus from PI 462312, PI 506764, and Wm82 using polymerase chain reaction (PCR). Five Rpp3C genes were identified in both Rpp3-resistant soybean lines, and co-silencing these genes compromised resistance to P. pachyrhizi. Gene expression analysis and sequence comparisons of the Rpp3C genes in PI 462312 and PI 506764 suggest that a single candidate gene, Rpp3C3, is responsible for Rpp3-mediated resistance. [Formula: see text] The author(s) have dedicated the work to the public domain under the Creative Commons CC0 "No Rights Reserved" license by waiving all of his or her rights to the work worldwide under copyright law, including all related and neighboring rights, to the extent allowed by law, 2024.

Overexpression of GmPAL Genes Enhances Soybean Resistance Against Heterodera glycines.

Soybean cyst nematode (Heterodera glycines, soybean cyst nematode [SCN]) disease adversely affects the yield of soybean and leads to billions of dollars in losses every year. To control the disease, it is necessary to study the resistance genes of the plant and their mechanisms. Isoflavonoids are secondary metabolites of the phenylalanine pathway, and they are synthesized in soybean. They are essential in plant response to biotic and abiotic stresses. In this study, we reported that phenylalanine ammonia-lyase (PAL) genes GmPALs involved in isoflavonoid biosynthesis, can positively regulate soybean resistance to SCN. Our previous study demonstrated that the expression of GmPAL genes in the resistant cultivar Huipizhi (HPZ) heidou are strongly induced by SCN. PAL is the rate-limiting enzyme that catalyzes the first step of phenylpropanoid metabolism, and it responds to biotic or abiotic stresses. Here, we demonstrate that the resistance of soybeans against SCN is suppressed by PAL inhibitor l-α-(aminooxy)-ß-phenylpropionic acid (L-AOPP) treatment. Overexpression of eight GmPAL genes caused diapause of nematodes in transgenic roots. In a petiole-feeding bioassay, we identified that two isoflavones, daidzein and genistein, could enhance resistance against SCN and suppress nematode development. This study thus reveals GmPAL-mediated resistance against SCN, information that has good application potential. The role of isoflavones in soybean resistance provides new information for the control of SCN. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Insights into Phakopsora pachyrhizi Effector-Effector Interactions.

The multifaceted role of pathogen-encoded effectors in plant-pathogen interactions is complex and not fully understood. Effectors operate within intricate host environments, interacting with host proteins and other effectors to modulate virulence. The complex interplay between effectors raises the concept of metaeffectors, wherein some effectors regulate the activity of others. While previous research has demonstrated the importance of effector repertoires in pathogen virulence, only a limited number of studies have investigated the interactions between these effectors. This study explores the interactions among Phakopsora pachyrhizi effector candidates (PpECs). P. pachyrhizi haustorial transcriptome analysis identified a collection of predicted PpECs. Among these, PpEC23 was found to interact with PpEC48, prompting further exploration into their potential interaction with other effectors. Here, we utilized a yeast two-hybrid screen to explore protein-protein interactions between PpECs. A split-luciferase complementation assay also demonstrated that these interactions could occur within soybean cells. Interestingly, PpEC48 displayed the ability to interact with several small cysteine-rich proteins (SCRPs), suggesting its affinity for this specific class of effectors. We show that these interactions involve a histidine-rich domain within PpEC48, emphasizing the significance of structural motifs in mediating effector interactions. The unique nature of PpEC48, showing no sequence matches in other organisms, suggests its relatively recent evolution and potential orphan gene status. Our work reveals insights into the intricate network of interactions among P. pachyrhizi effector-effector interactions. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

A comprehensive analysis of the WRKY family in soybean and functional analysis of GmWRKY164-GmGSL7c in resistance to soybean mosaic virus.

BACKGROUND: Soybean mosaic disease caused by soybean mosaic virus (SMV) is one of the most devastating and widespread diseases in soybean producing areas worldwide. The WRKY transcription factors (TFs) are widely involved in plant development and stress responses. However, the roles of the GmWRKY TFs in resistance to SMV are largely unclear. RESULTS: Here, 185 GmWRKYs were characterized in soybean (Glycine max), among which 60 GmWRKY genes were differentially expressed during SMV infection according to the transcriptome data. The transcriptome data and RT-qPCR results showed that the expression of GmWRKY164 decreased after imidazole treatment and had higher expression levels in the incompatible combination between soybean cultivar variety Jidou 7 and SMV strain N3. Remarkably, the silencing of GmWRKY164 reduced callose deposition and enhanced virus spread during SMV infection. In addition, the transcript levels of the GmGSL7c were dramatically lower upon the silencing of GmWRKY164. Furthermore, EMSA and ChIP-qPCR revealed that GmWRKY164 can directly bind to the promoter of GmGSL7c, which contains the W-box element. CONCLUSION: Our findings suggest that GmWRKY164 plays a positive role in resistance to SMV infection by regulating the expression of GmGSL7c, resulting in the deposition of callose and the inhibition of viral movement, which provides guidance for future studies in understanding virus-resistance mechanisms in soybean.

Elucidating the ecophysiology of soybean pod-sucking stinkbug Riptortus pedestris (Hemiptera: Alydidae) based on de novo genome assembly and transcriptome analysis.

Food security is important for the ever-growing global population. Soybean, Glycine max (L.) Merr., is cultivated worldwide providing a key source of food, protein and oil. Hence, it is imperative to maintain or to increase its yield under different conditions including challenges caused by abiotic and biotic stresses. In recent years, the soybean pod-sucking stinkbug Riptortus pedestris has emerged as an important agricultural insect pest in East, South and Southeast Asia. Here, we present a genomics resource for R. pedestris including its genome assembly, messenger RNA (mRNA) and microRNA (miRNA) transcriptomes at different developmental stages and from different organs. As insect hormone biosynthesis genes (genes involved in metamorphosis) and their regulators such as miRNAs are potential targets for pest control, we analyzed the sesquiterpenoid (juvenile) and ecdysteroid (molting) hormone biosynthesis pathway genes including their miRNAs and relevant neuropeptides. Temporal gene expression changes of these insect hormone biosynthesis pathways were observed at different developmental stages. Similarly, a diet-specific response in gene expression was also observed in both head and salivary glands. Furthermore, we observed that microRNAs (bantam, miR-14, miR-316, and miR-263) of R. pedestris fed with different types of soybeans were differentially expressed in the salivary glands indicating a diet-specific response. Interestingly, the opposite arms of miR-281 (-5p and -3p), a miRNA involved in regulating development, were predicted to target Hmgs genes of R. pedestris and soybean, respectively. These observations among others highlight stinkbug's responses as a function of its interaction with soybean. In brief, the results of this study not only present salient findings that could be of potential use in pest management and mitigation but also provide an invaluable resource for R. pedestris as an insect model to facilitate studies on plant-pest interactions.