Whole-genome sequencing and intensive analysis of the undomesticated soybean (Glycine soja Sieb. and Zucc.) genome

Shared phylogeographic patterns and environmental responses of co-distributed soybean pests: Insights from comparative phylogeographic studies of Riptortus pedestris and Riptortus linearis in the subtropics of East Asia

Comparative phylogeographic studies of closely related species sharing co-distribution areas can elucidate the role of shared historical factors and environmental changes in shaping their phylogeographic pattern. The bean bugs, Riptortus pedestris and Riptortus linearis, which both inhabit subtropical regions in East Asia, are recognized as highly destructive soybean pests. Many previous studies have investigated the biological characteristics, pheromones, chemicals and control mechanisms of these two pests, but few studies have explored their phylogeographic patterns and underlying factors. In this study, we generated a double-digest restriction site-associated DNA sequencing (ddRAD-seq) dataset to investigate phylogeographic patterns and construct ecological niche models (ENM) for both Riptortus species. Our findings revealed similar niche occupancies and population genetic structures between the two species, with each comprising two phylogeographic lineages (i.e., the mainland China and the Indochina Peninsula clades) that diverged approximately 0.1 and 0.3 million years ago, respectively. This divergence likely resulted from the combined effects of temperatures variation and geographical barriers in the mountainous regions of Southwest China. Further demographic history and ENM analyses suggested that both pests underwent rapid expansion prior to the Last Glacial Maximum (LGM). Furthermore, ENM predicts a northward shift of both pests into new soybean-producing regions due to global warming. Our study indicated that co-distribution soybean pests with overlapping ecological niches and similar life histories in subtropical regions of East Asia exhibit congruent phylogeographic and demographic patterns in response to shared historical biogeographic drivers.

In Silico Comparison of WRKY Transcription Factors in Wild and Cultivated Soybean and Their Co-expression Network Arbitrating Disease Resistance

WRKY Transcription factors (TFs) play critical roles in plant defence mechanisms that are activated in response to biotic and abiotic stresses. However, information on the Glycine soja WRKYs (GsoWRKYs) is scarce. Owing to its importance in soybean breeding, here we identified putative WRKY TFs in wild soybean, and compared the results with Glycine max WRKYs (GmaWRKYs) by phylogenetic, conserved motif, and duplication analyses. Moreover, we explored the expression trends of WRKYs in G. max (oomycete, fungi, virus, bacteria, and soybean cyst nematode) and G. soja (soybean cyst nematode), and identified commonly expressed WRKYs and their co-expressed genes. We identified, 181 and 180 putative WRKYs in G. max and G. soja, respectively. Though the number of WRKYs in both studied species is almost the same, they differ in many ways, i.e., the number of WRKYs on corresponding chromosomes, conserved domain structures, WRKYGQK motif variants, and zinc-finger motifs. WRKYs in both species grouped in three major clads, i.e., I-III, where group-II had sub-clads IIa-IIe. We found that GsoWRKYs expanded mostly through segmental duplication. A large number of WRKYs were expressed in response to biotic stresses, i.e., Phakospora pachyrhizi, Phytoplasma, Heterodera glycines, Macrophomina phaseolina, and Soybean mosaic virus; 56 GmaWRKYs were commonly expressed in soybean plants infected with these diseases. Finally, 30 and 63 GmaWRKYs and GsoWRKYs co-expressed with 205 and 123 non-WRKY genes, respectively, indicating that WRKYs play essential roles in biotic stress tolerance in Glycine species.

Regulation of seed traits in soybean

Soybean (Glycine max) is an essential economic crop that provides vegetative oil and protein for humans, worldwide. Increasing soybean yield as well as improving seed quality is of great importance. Seed weight/size, oil and protein content are the three major traits determining seed quality, and seed weight also influences soybean yield. In recent years, the availability of soybean omics data and the development of related techniques have paved the way for better research on soybean functional genomics, providing a comprehensive understanding of gene functions. This review summarizes the regulatory genes that influence seed size/weight, oil content and protein content in soybean. We also provided a general overview of the pleiotropic effect for the genes in controlling seed traits and environmental stresses. Ultimately, it is expected that this review will be beneficial in breeding improved traits in soybean.

Global analysis of seed transcriptomes reveals a novel PLATZ regulator for seed size and weight control in soybean

Seed size and weight are important factors that influence soybean yield. Combining the weighted gene co-expression network analysis (WGCNA) of 45 soybean accessions and gene dynamic changes in seeds at seven developmental stages, we identified candidate genes that may control the seed size/weight. Among these, a PLATZ-type regulator overlapping with 10 seed weight QTLs was further investigated. This zinc-finger transcriptional regulator, named as GmPLATZ, is required for the promotion of seed size and weight in soybean. The GmPLATZ may exert its functions through direct binding to the promoters and activation of the expression of cyclin genes and GmGA20OX for cell proliferation. Overexpression of the GmGA20OX enhanced seed size/weight in soybean. We further found that the GmPLATZ binds to a 32-bp sequence containing a core palindromic element AATGCGCATT. Spacing of the flanking sequences beyond the core element facilitated GmPLATZ binding. An elite haplotype Hap3 was also identified to have higher promoter activity and correlated with higher gene expression and higher seed weight. Orthologues of the GmPLATZ from rice and Arabidopsis play similar roles in seeds. Our study reveals a novel module of GmPLATZ-GmGA20OX/cyclins in regulating seed size and weight and provides valuable targets for breeding of crops with desirable agronomic traits.

A loss-of-function mutant allele of a glycosyl hydrolase gene has been co-opted for seed weight control during soybean domestication

The resultant DNA from loss-of-function mutation can be recruited in biological evolution and development. Here, we present such a rare and potential case of "to gain by loss" as a neomorphic mutation during soybean domestication for increasing seed weight. Using a population derived from a chromosome segment substitution line of Glycine max (SN14) and Glycine soja (ZYD06), a quantitative trait locus (QTL) of 100-seed weight (qHSW) was mapped on chromosome 11, corresponding to a truncated β-1, 3-glucosidase (βGlu) gene. The novel gene hsw results from a 14-bp deletion, causing a frameshift mutation and a premature stop codon in the βGlu. In contrast to HSW, the hsw completely lost βGlu activity and function but acquired a novel function to promote cell expansion, thus increasing seed weight. Overexpressing hsw instead of HSW produced large soybean seeds, and surprisingly, truncating hsw via gene editing further increased the seed size. We further found that the core 21-aa peptide of hsw and its variants acted as a promoter of seed size. Transcriptomic variation in these transgenic soybean lines substantiated the integration hsw into cell and seed size control. Moreover, the hsw allele underwent selection and expansion during soybean domestication and improvement. Our work cloned a likely domesticated QTL controlling soybean seed weight, revealed a novel genetic variation and mechanism in soybean domestication, and provided new insight into crop domestication and breeding, and plant evolution.

Genomic Variations Explorer (GenVarX): a toolset for annotating promoter and CNV regions using genotypic and phenotypic differences

The rapid growth of sequencing technology and its increasing popularity in biology-related research over the years has made whole genome re-sequencing (WGRS) data become widely available. A large amount of WGRS data can unlock the knowledge gap between genomics and phenomics through gaining an understanding of the genomic variations that can lead to phenotype changes. These genomic variations are usually comprised of allele and structural changes in DNA, and these changes can affect the regulatory mechanisms causing changes in gene expression and altering the phenotypes of organisms. In this research work, we created the GenVarX toolset, that is backed by transcription factor binding sequence data in promoter regions, the copy number variations data, SNPs and Indels data, and phenotypes data which can potentially provide insights about phenotypic differences and solve compelling questions in plant research. Analytics-wise, we have developed strategies to better utilize the WGRS data and mine the data using efficient data processing scripts, libraries, tools, and frameworks to create the interactive and visualization-enhanced GenVarX toolset that encompasses both promoter regions and copy number variation analysis components. The main capabilities of the GenVarX toolset are to provide easy-to-use interfaces for users to perform queries, visualize data, and interact with the data. Based on different input windows on the user interface, users can provide inputs corresponding to each field and submit the information as a query. The data returned on the results page is usually displayed in a tabular fashion. In addition, interactive figures are also included in the toolset to facilitate the visualization of statistical results or tool outputs. Currently, the GenVarX toolset supports soybean, rice, and Arabidopsis. The researchers can access the soybean GenVarX toolset from SoyKB via https://soykb.org/SoybeanGenVarX/, rice GenVarX toolset, and Arabidopsis GenVarX toolset from KBCommons web portal with links https://kbcommons.org/system/tools/GenVarX/Osativa and https://kbcommons.org/system/tools/GenVarX/Athaliana, respectively.

Genetic analysis of protein content and oil content in soybean by genome-wide association study

Soybean seed protein content (PC) and oil content (OC) have important economic value. Detecting the loci/gene related to PC and OC is important for the marker-assisted selection (MAS) breeding of soybean. To detect the stable and new loci for PC and OC, a total of 320 soybean accessions collected from the major soybean-growing countries were used to conduct a genome-wide association study (GWAS) by resequencing. The PC ranged from 37.8% to 46.5% with an average of 41.1% and the OC ranged from 16.7% to 22.6% with an average of 21.0%. In total, 23 and 29 loci were identified, explaining 3.4%-15.4% and 5.1%-16.3% of the phenotypic variations for PC and OC, respectively. Of these, eight and five loci for PC and OC, respectively, overlapped previously reported loci and the other 15 and 24 loci were newly identified. In addition, nine candidate genes were identified, which are known to be involved in protein and oil biosynthesis/metabolism, including lipid transport and metabolism, signal transduction, and plant development pathway. These results uncover the genetic basis of soybean protein and oil biosynthesis and could be used to accelerate the progress in enhancing soybean PC and OC.

Genetic mapping and functional genomics of soybean seed protein

Soybean is an utterly important crop for high-quality meal protein and vegetative oil. Soybean seed protein content has become a key factor in nutrients for livestock feed as well as human dietary consumption. Genetic improvement of soybean seed protein is highly desired to meet the demands of rapidly growing world population. Molecular mapping and genomic analysis in soybean have identified many quantitative trait loci (QTL) underlying seed protein content control. Exploring the mechanisms of seed storage protein regulation will be helpful to achieve the improvement of protein content. However, the practice of breeding higher protein soybean is challenging because soybean seed protein is negatively correlated with seed oil content and yield. To overcome the limitation of such inverse relationship, deeper insights into the property and genetic control of seed protein are required. Recent advances of soybean genomics have strongly enhanced the understandings for molecular mechanisms of soybean with better seed quality. Here, we review the research progress in the genetic characteristics of soybean storage protein, and up-to-date advances of molecular mappings and genomics of soybean protein. The key factors underlying the mechanisms of the negative correlation between protein and oil in soybean seeds are elaborated. We also briefly discuss the future prospects of breaking the bottleneck of the negative correlation to develop high protein soybean without penalty of oil and yield.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01373-5.

The Allele Catalog Tool: a web-based interactive tool for allele discovery and analysis

BACKGROUND

The advancement of sequencing technologies today has made a plethora of whole-genome re-sequenced (WGRS) data publicly available. However, research utilizing the WGRS data without further configuration is nearly impossible. To solve this problem, our research group has developed an interactive Allele Catalog Tool to enable researchers to explore the coding region allelic variation present in over 1,000 re-sequenced accessions each for soybean, Arabidopsis, and maize.

RESULTS

The Allele Catalog Tool was designed originally with soybean genomic data and resources. The Allele Catalog datasets were generated using our variant calling pipeline (SnakyVC) and the Allele Catalog pipeline (AlleleCatalog). The variant calling pipeline is developed to parallelly process raw sequencing reads to generate the Variant Call Format (VCF) files, and the Allele Catalog pipeline takes VCF files to perform imputations, functional effect predictions, and assemble alleles for each gene to generate curated Allele Catalog datasets. Both pipelines were utilized to generate the data panels (VCF files and Allele Catalog files) in which the accessions of the WGRS datasets were collected from various sources, currently representing over 1,000 diverse accessions for soybean, Arabidopsis, and maize individually. The main features of the Allele Catalog Tool include data query, visualization of results, categorical filtering, and download functions. Queries are performed from user input, and results are a tabular format of summary results by categorical description and genotype results of the alleles for each gene. The categorical information is specific to each species; additionally, available detailed meta-information is provided in modal popups. The genotypic information contains the variant positions, reference or alternate genotypes, the functional effect classes, and the amino-acid changes of each accession. Besides that, the results can also be downloaded for other research purposes.

CONCLUSIONS

The Allele Catalog Tool is a web-based tool that currently supports three species: soybean, Arabidopsis, and maize. The Soybean Allele Catalog Tool is hosted on the SoyKB website ( https://soykb.org/SoybeanAlleleCatalogTool/ ), while the Allele Catalog Tool for Arabidopsis and maize is hosted on the KBCommons website ( https://kbcommons.org/system/tools/AlleleCatalogTool/Zmays and https://kbcommons.org/system/tools/AlleleCatalogTool/Athaliana ). Researchers can use this tool to connect variant alleles of genes with meta-information of species.

Genome-Wide Identification and Characterization of the PPO Gene Family in Cotton (Gossypium) and Their Expression Variations Responding to Verticillium Wilt Infection

Polyphenol oxidases (PPOs) are copper-binding metalloproteinases encoded by nuclear genes, ubiquitously existing in the plastids of microorganisms, plants, and animals. As one of the important defense enzymes, PPOs have been reported to participate in the resistant processes that respond to diseases and insect pests in multiple plant species. However, PPO gene identification and characterization in cotton and their expression patterns under Verticillium wilt (VW) treatment have not been clearly studied. In this study, 7, 8, 14, and 16 PPO genes were separately identified from Gossypium arboreum, G. raimondii, G. hirsutum, and G. barbadense, respectively, which were distributed within 23 chromosomes, though mainly gathered in chromosome 6. The phylogenetic tree manifested that all the PPOs from four cotton species and 14 other plants were divided into seven groups, and the analyses of the conserved motifs and nucleotide sequences showed highly similar characteristics of the gene structure and domains in the cotton PPO genes. The dramatically expressed differences were observed among the different organs at various stages of growth and development or under the diverse stresses referred to in the published RNA-seq data. Quantitative real-time PCR (qRT-PCR) experiments were also performed on the GhPPO genes in the roots, stems, and leaves of VW-resistant MBI8255 and VW-susceptible CCRI36 infected with Verticillium dahliae V991, proving the strong correlation between PPO activity and VW resistance. A comprehensive analysis conducted on cotton PPO genes contributes to the screening of the candidate genes for subsequent biological function studies, which is also of great significance for the in-depth understanding of the molecular genetic basis of cotton resistance to VW.

Understandings and future challenges in soybean functional genomics and molecular breeding

Soybean (Glycine max) is a major source of plant protein and oil. Soybean breeding has benefited from advances in functional genomics. In particular, the release of soybean reference genomes has advanced our understanding of soybean adaptation to soil nutrient deficiencies, the molecular mechanism of symbiotic nitrogen (N) fixation, biotic and abiotic stress tolerance, and the roles of flowering time in regional adaptation, plant architecture, and seed yield and quality. Nevertheless, many challenges remain for soybean functional genomics and molecular breeding, mainly related to improving grain yield through high-density planting, maize-soybean intercropping, taking advantage of wild resources, utilization of heterosis, genomic prediction and selection breeding, and precise breeding through genome editing. This review summarizes the current progress in soybean functional genomics and directs future challenges for molecular breeding of soybean.

The dicing activity of DCL3 and DCL4 is negatively affected by flavonoids

The dicing activities of DCL3 and DCL4 are inhibited by accumulated metabolites in soybean leaves. Epicatechin and 7,4'-dihydroxyflavone inhibited Arabidopsis DCL3 and DCL4 in vitro. Flavonoids are major secondary metabolites in plants, and soybean (Glycine max L.) is a representative plant that accumulates flavonoids, including isoflavonoids, to high levels. Naturally-occurring RNA interference (RNAi) against the chalcone synthase (CHS) gene represses flavonoid (anthocyanin) biosynthesis in an organ-specific manner, resulting in a colorless (yellow) seed coat in many soybean cultivars. To better understand seed coat-specific naturally-occurring RNAi in soybean, we characterized soybean Dicer-like (DCL) 3 and 4, which play critical roles in RNAi. Using a previously established dicing assay, two dicing activities producing 24- and 21-nt siRNAs, corresponding to DCL3 and DCL4, respectively, were detected in soybean. Dicing activity was detected in colorless seed coats where RNAi against CHS genes was found, but no dicing activity was detected in leaves where CHS expression was prevalent. Biochemical analysis revealed that soybean leaves contained two types of inhibitors effective for Arabidopsis Dicers (AtDCL3 and AtDCL4), one of which was a heat-labile high molecular weight compound of 50 to 100 kD while another was a low molecular weight substance. We found that some flavonoids, such as epicatechin and 7,4'-dihydroxyflavone, inhibited both AtDCL3 and AtDCL4, but AtDCL4 was more sensitive to these flavonoids than AtDCL3. These results suggest that flavonoids inhibit the dicing activity of DCL4 and thereby attenuate RNAi in soybean leaves.

Comprehensive Metabolite Profiling of Four Different Beans Fermented by Aspergillus oryzae

Fermented bean products are used worldwide; most of the products are made using only a few kinds of beans. However, the metabolite changes and contents in the beans generally used during fermentation are unrevealed. Therefore, we selected four different beans (soybean, Glycine max, GM; wild soybean, Glycine soja, GS; common bean, Phaseolus vulgaris, PV; and hyacinth bean, Lablab purpureus, LP) that are the most widely consumed and fermented with Aspergillus oryzae. Then, metabolome and multivariate statistical analysis were performed to figure out metabolite changes during fermentation. In the four beans, carbohydrates were decreased, but amino acids and fatty acids were increased in the four beans as they fermented. The relative amounts of amino acids were relatively abundant in fermented PV and LP as compared to other beans. In contrast, isoflavone aglycones (e.g., daidzein, glycitein, and genistein) and DDMP-conjugated soyasaponins (e.g., soyasaponins βa and γg) were increased in GM and GS during fermentation. Notably, these metabolite changes were more significant in GS than GM. In addition, the increase of antioxidant activity in fermented GS was significant compared to other beans. We expect our research provides a basis to extend choice for bean fermentation for consumers and food producers.

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.

Coupled and decoupled legumes and cereals in prehistoric northern and southern China

Legumes and cereals, which provide different nutrients, are cultivated as coupled crops in most centers of plant domestication worldwide. However, as the only legume domesticated in China, the spatio-temporal distribution of soybeans and its status in the millet- and rice-based agricultural system of the Neolithic and Bronze Ages remains elusive. Here, archaeobotanical evidence of soybeans (n=254), millet (n=462), rice (n=482), and zooarchaeological evidence of fish (n=138) were synthesized to elucidate the phenomenon of coupled or decoupled cereals and legumes in prehistoric China. During the Neolithic and Bronze Ages, soybeans was mostly confined to northern China and rarely found in southern China, serving as a companion to millet. In contrast, fish remains have been widely found in southern China, indicating a continuous reliance on fish as a staple food besides rice. Thus, an antipodal pattern of millet-soybeans and rice-fish agricultural systems may have been established in northern and southern China since the late Yangshao period (6000-5000 cal BP) respectively. These two agricultural systems were not only complementary in terms of diet, but they also exhibited positive interactions and feedback in the coculture system. Consequently, these two systems enabled the sustainable intensification of agriculture and served as the basis for the emergence of complex societies and early states in the Yellow and Yangtze Rivers.

Genome-wide identification and analysis of LOX genes in soybean cultivar “Zhonghuang 13”

Lipoxygenases (LOXs; EC1.13.11.12) are a family of iron- or manganese-containing dioxygenases that catalyze the oxygenation of polyunsaturated fatty acids (PUFAs) and play important roles in plant growth, development, and stress response. In this study, a total of 36 LOX gene family members were identified and annotated in Zhonghuang 13, a soybean cultivar bred by Chinese scientists in 2001. Sanger sequencing of the GmLOX1-coding sequence and colorimetric assays for the GmLOX1 protein showed that Zhonghuang 13 possessed the GmLOX1 gene. These LOX genes are divided into three subfamilies: 9-LOX, type Ⅰ 13-LOX and type II 13-LOX. In the 13-LOX group, the number of GmLOX members was the highest. These GmLOX genes are unevenly distributed on chromosomes 3, 7, 8, 10, 11, 12, 13, 15, 16, 19, and 20. Most of the 13-LOX genes exist in the form of gene clusters, indicating that these genes may originate from tandem duplications. The analysis of duplicated gene pairs showed that GmLOX genes underwent purifying selective pressure during evolution. The gene structures and conserved functional domains of these genes are quite similar. Compared to the orthologous gene pairs of LOX genes between wild soybean (Glycine soja W05) and Zhonghuang 13, the sequences of most gene pairs are relatively conserved. Many cis-elements are present in the promoter region and are involved in stress response, growth and development, hormone response and light response. The tissue-specific gene expression of GmLOX genes was evaluated. Represented by GmLOX1, GmLOX2, and GmLOX3, which were expressed at extremely high levels in seeds, they showed the characteristics of specific expression. This study provides detailed information on soybean lipoxygenase gene family members in Zhonghuang 13, which lays a foundation for further research.

Whole-Genome Identification of APX and CAT Gene Families in Cultivated and Wild Soybeans and Their Regulatory Function in Plant Development and Stress Response

Plants coevolved with their antioxidant defense systems, which detoxify and adjust levels of reactive oxygen species (ROS) under multiple plant stresses. We performed whole-genome identification of ascorbate peroxidase (APX) and catalase (CAT) families in cultivated and wild soybeans. In cultivated and wild soybean genomes, we identified 11 and 10 APX genes, respectively, whereas the numbers of identified CAT genes were four in each species. Comparative phylogenetic analysis revealed more homology among cultivated and wild soybeans relative to other legumes. Exon/intron structure, motif and synteny blocks are conserved in cultivated and wild species. According to the Ka/Ks value, purifying selection is a major force for evolution of these gene families in wild soybean; however, the APX gene family was evolved by both positive and purifying selection in cultivated soybean. Segmental duplication was a major factor involved in the expansion of APX and CAT genes. Expression patterns revealed that APX and CAT genes are differentially expressed across fourteen different soybean tissues under water deficit (WD), heat stress (HS) and combined drought plus heat stress (WD + HS). Altogether, the current study provides broad insights into these gene families in soybeans. Our results indicate that APX and CAT gene families modulate multiple stress response in soybeans.

Identification of genetic loci conferring seed coat color based on a high-density map in soybean

Seed coat color is a typical evolutionary trait. Identification of the genetic loci that control seed coat color during the domestication of wild soybean could clarify the genetic variations between cultivated and wild soybean. We used 276 F₁₀ recombinant inbred lines (RILs) from the cross between a cultivated soybean (JY47) and a wild soybean (ZYD00321) as the materials to identify the quantitative trait loci (QTLs) for seed coat color. We constructed a high-density genetic map using re-sequencing technology. The average distance between adjacent markers was 0.31 cM on this map, comprising 9,083 bin markers. We identified two stable QTLs (qSC08 and qSC11) for seed coat color using this map, which, respectively, explained 21.933 and 26.934% of the phenotypic variation. Two candidate genes (CHS3C and CHS4A) in qSC08 were identified according to the parental re-sequencing data and gene function annotations. Five genes (LOC100786658, LOC100801691, LOC100806824, LOC100795475, and LOC100787559) were predicted in the novel QTL qSC11, which, according to gene function annotations, might control seed coat color. This result could facilitate the identification of beneficial genes from wild soybean and provide useful information to clarify the genetic variations for seed coat color in cultivated and wild soybean.

Unraveling Origin, History, Genetics, and Strategies for Accelerated Domestication and Diversification of Food Legumes

Domestication is a dynamic and ongoing process of transforming wild species into cultivated species by selecting desirable agricultural plant features to meet human needs such as taste, yield, storage, and cultivation practices. Human plant domestication began in the Fertile Crescent around 12,000 years ago and spread throughout the world, including China, Mesoamerica, the Andes and Near Oceania, Sub-Saharan Africa, and eastern North America. Indus valley civilizations have played a great role in the domestication of grain legumes. Crops, such as pigeon pea, black gram, green gram, lablab bean, moth bean, and horse gram, originated in the Indian subcontinent, and Neolithic archaeological records indicate that these crops were first domesticated by early civilizations in the region. The domestication and evolution of wild ancestors into today’s elite cultivars are important contributors to global food supply and agricultural crop improvement. In addition, food legumes contribute to food security by protecting human health and minimize climate change impacts. During the domestication process, legume crop species have undergone a severe genetic diversity loss, and only a very narrow range of variability is retained in the cultivars. Further reduction in genetic diversity occurred during seed dispersal and movement across the continents. In general, only a few traits, such as shattering resistance, seed dormancy loss, stem growth behavior, flowering–maturity period, and yield traits, have prominence in the domestication process across the species. Thus, identification and knowledge of domestication responsive loci were often useful in accelerating new species’ domestication. The genes and metabolic pathways responsible for the significant alterations that occurred as an outcome of domestication might aid in the quick domestication of novel crops. Further, recent advances in “omics” sciences, gene-editing technologies, and functional analysis will accelerate the domestication and crop improvement of new crop species without losing much genetic diversity. In this review, we have discussed about the origin, center of diversity, and seed movement of major food legumes, which will be useful in the exploration and utilization of genetic diversity in crop improvement. Further, we have discussed about the major genes/QTLs associated with the domestication syndrome in pulse crops and the future strategies to improve the food legume crops.

Expanding the gene pool for soybean improvement with its wild relatives

Genetic diversity is a cornerstone of crop improvement, However, cultivated soybean (Glycine max) has undergone several genetic bottlenecks, including domestication in China, the introduction of landraces to other areas of the world and, latterly, selective breeding, leading to low genetic diversity the poses a major obstacle to soybean improvement. By contrast, there remains a relatively high level of genetic diversity in soybean's wild relatives, especially the perennial soybeans (Glycine subgenus Glycine), which could serve as potential gene pools for improving soybean cultivars. Wild soybeans are phylogenetically diversified and adapted to various habitats, harboring resistance to various biotic and abiotic stresses. Advances in genome and transcriptome sequencing enable alleles associated with desirable traits that were lost during domestication of soybean to be discovered in wild soybean. The collection and conservation of soybean wild relatives and the dissection of their genomic features will accelerate soybean breeding and facilitate sustainable agriculture and food production.

Multi-Omics Techniques for Soybean Molecular Breeding

Soybean is a major crop that provides essential protein and oil for food and feed. Since its origin in China over 5000 years ago, soybean has spread throughout the world, becoming the second most important vegetable oil crop and the primary source of plant protein for global consumption. From early domestication and artificial selection through hybridization and ultimately molecular breeding, the history of soybean breeding parallels major advances in plant science throughout the centuries. Now, rapid progress in plant omics is ushering in a new era of precision design breeding, exemplified by the engineering of elite soybean varieties with specific oil compositions to meet various end-use targets. The assembly of soybean reference genomes, made possible by the development of genome sequencing technology and bioinformatics over the past 20 years, was a great step forward in soybean research. It facilitated advances in soybean transcriptomics, proteomics, metabolomics, and phenomics, all of which paved the way for an integrated approach to molecular breeding in soybean. In this review, we summarize the latest progress in omics research, highlight novel findings made possible by omics techniques, note current drawbacks and areas for further research, and suggest that an efficient multi-omics approach may accelerate soybean breeding in the future. This review will be of interest not only to soybean breeders but also to researchers interested in the use of cutting-edge omics technologies for crop research and improvement.

Development of a versatile resource for post-genomic research through consolidating and characterizing 1500 diverse wild and cultivated soybean genomes

BACKGROUND

With advances in next-generation sequencing technologies, an unprecedented amount of soybean accessions has been sequenced by many individual studies and made available as raw sequencing reads for post-genomic research.

RESULTS

To develop a consolidated and user-friendly genomic resource for post-genomic research, we consolidated the raw resequencing data of 1465 soybean genomes available in the public and 91 highly diverse wild soybean genomes newly sequenced. These altogether provided a collection of 1556 sequenced genomes of 1501 diverse accessions (1.5 K). The collection comprises of wild, landraces and elite cultivars of soybean that were grown in East Asia or major soybean cultivating areas around the world. Our extensive sequence analysis discovered 32 million single nucleotide polymorphisms (32mSNPs) and revealed a SNP density of 30 SNPs/kb and 12 non-synonymous SNPs/gene reflecting a high structural and functional genomic diversity of the new collection. Each SNP was annotated with 30 categories of structural and/or functional information. We further identified paired accessions between the 1.5 K and 20,087 (20 K) accessions in US collection as genomic "equivalent" accessions sharing the highest genomic identity for minimizing the barriers in soybean germplasm exchange between countries. We also exemplified the utility of 32mSNPs in enhancing post-genomics research through in-silico genotyping, high-resolution GWAS, discovering and/or characterizing genes and alleles/mutations, identifying germplasms containing beneficial alleles that are potentially experiencing artificial selection.

CONCLUSION

The comprehensive analysis of publicly available large-scale genome sequencing data of diverse cultivated accessions and the newly in-house sequenced wild accessions greatly increased the soybean genome-wide variation resolution. This could facilitate a variety of genetic and molecular-level analyses in soybean. The 32mSNPs and 1.5 K accessions with their comprehensive annotation have been made available at the SoyBase and Ag Data Commons. The dataset could further serve as a versatile and expandable core resource for exploring the exponentially increasing genome sequencing data for a variety of post-genomic research.

The chromosome-level genome assembly of Astragalus sinicus and comparative genomic analyses provide new resources and insights for understanding legume-rhizobial interactions

The legume species Astragalus sinicus (Chinese milk vetch [CMV]) has been widely cultivated for centuries in southern China as one of the most important green manures/cover crops for improving rice productivity and preventing soil degeneration. In this study, we generated the first chromosome-scale reference genome of CMV by combining PacBio and Illumina sequencing with high-throughput chromatin conformation capture (Hi-C) technology. The CMV genome was 595.52 Mb in length, with a contig N50 size of 1.50 Mb. Long terminal repeats (LTRs) had been amplified and contributed to genome size expansion in CMV. CMV has undergone two whole-genome duplication (WGD) events, and the genes retained after the WGD shared by Papilionoideae species shaped the rhizobial symbiosis and the hormonal regulation of nodulation. The chalcone synthase (CHS) gene family was expanded and was expressed primarily in the roots of CMV. Intriguingly, we found that resistance genes were more highly expressed in roots than in nodules of legume species, suggesting that their expression may be increased to bolster plant immunity in roots to cope with pathogen infection in legumes. Our work sheds light on the genetic basis of nodulation and symbiosis in CMV and provides a benchmark for accelerating genetic research and molecular breeding in the future.

Sequencing the USDA core soybean collection reveals gene loss during domestication and breeding

The gene content of plants varies between individuals of the same species due to gene presence/absence variation, and selection can alter the frequency of specific genes in a population. Selection during domestication and breeding will modify the genomic landscape, though the nature of these modifications is only understood for specific genes or on a more general level (e.g., by a loss of genetic diversity). Here we have assembled and analyzed a soybean (Glycine spp.) pangenome representing more than 1,000 soybean accessions derived from the USDA Soybean Germplasm Collection, including both wild and cultivated lineages, to assess genomewide changes in gene and allele frequency during domestication and breeding. We identified 3,765 genes that are absent from the Lee reference genome assembly and assessed the presence/absence of all genes across this population. In addition to a loss of genetic diversity, we found a significant reduction in the average number of protein-coding genes per individual during domestication and subsequent breeding, though with some genes and allelic variants increasing in frequency associated with selection for agronomic traits. This analysis provides a genomic perspective of domestication and breeding in this important oilseed crop.

A de novo regulation design shows an effectiveness in altering plant secondary metabolism

•

Arabidopsis PAP1 and TT8 transcription factors and cis-regulatory elements of tobacco JAZs are for the first time selected as molecular tools for a De Novo regulation design.

•

This design creates a distant pathway-cross regulation (DPCR) that effectively downregulates the alkaloid biosynthesis in tobacco.

•

The DPCR significantly reduces carcinogenic compounds.

•

This design unearths two novel regulatory functions of PAP1 and TT8 and their complex and a new regulation mechanism of four NtJAZs.

Introduction

Transcription factors (TFs) and cis-regulatory elements (CREs) control gene transcripts involved in various biological processes. We hypothesize that TFs and CREs can be effective molecular tools for De Novo regulation designs to engineer plants.

Objectives

We selected two Arabidopsis TF types and two tobacco CRE types to design a De Novo regulation and evaluated its effectiveness in plant engineering.

Methods

G-box and MYB recognition elements (MREs) were identified in four Nicotiana tabacum JAZs (NtJAZs) promoters. MRE-like and G-box like elements were identified in one nicotine pathway gene promoter. TF screening led to select Arabidopsis Production of Anthocyanin Pigment 1 (PAP1/MYB) and Transparent Testa 8 (TT8/bHLH). Two NtJAZ and two nicotine pathway gene promoters were cloned from commercial Narrow Leaf Madole (NL) and KY171 (KY) tobacco cultivars. Electrophoretic mobility shift assay (EMSA), cross-linked chromatin immunoprecipitation (ChIP), and dual-luciferase assays were performed to test the promoter binding and activation by PAP1 (P), TT8 (T), PAP1/TT8 together, and the PAP1/TT8/Transparent Testa Glabra 1 (TTG1) complex. A DNA cassette was designed and then synthesized for stacking and expressing PAP1 and TT8 together. Three years of field trials were performed by following industrial and GMO protocols. Gene expression and metabolic profiling were completed to characterize plant secondary metabolism.

Results

PAP1, TT8, PAP1/TT8, and the PAP1/TT8/TTG1 complex bound to and activated NtJAZ promoters but did not bind to nicotine pathway gene promoters. The engineered red P + T plants significantly upregulated four NtJAZs but downregulated the tobacco alkaloid biosynthesis. Field trials showed significant reduction of five tobacco alkaloids and four carcinogenic tobacco specific nitrosamines in most or all cured leaves of engineered P + T and PAP1 genotypes.

Conclusion

G-boxes, MREs, and two TF types are appropriate molecular tools for a De Novo regulation design to create a novel distant-pathway cross regulation for altering plant secondary metabolism.

Progress in soybean functional genomics over the past decade

Soybean is one of the most important oilseed and fodder crops. Benefiting from the efforts of soybean breeders and the development of breeding technology, large number of germplasm has been generated over the last 100 years. Nevertheless, soybean breeding needs to be accelerated to meet the needs of a growing world population, to promote sustainable agriculture and to address future environmental changes. The acceleration is highly reliant on the discoveries in gene functional studies. The release of the reference soybean genome in 2010 has significantly facilitated the advance in soybean functional genomics. Here, we review the research progress in soybean omics (genomics, transcriptomics, epigenomics and proteomics), germplasm development (germplasm resources and databases), gene discovery (genes that are responsible for important soybean traits including yield, flowering and maturity, seed quality, stress resistance, nodulation and domestication) and transformation technology during the past decade. At the end, we also briefly discuss current challenges and future directions.

An SGSGeneloss-Based Method for Constructing a Gene Presence-Absence Table Using Mosdepth

Presence-absence variants (PAV) are genomic regions present in some individuals of a species, but not others. PAVs have been shown to contribute to genomic diversity, especially in bacteria and plants. These structural variations have been linked to traits and can be used to track a species' evolutionary history. PAVs are usually called by aligning short read sequence data from one or more individuals to a reference genome or pangenome assembly, and then comparing coverage. Regions where reads do not align define absence in that individual, and the regions are classified as PAVs. The method below details how to align sequence reads to a reference and how to use the sequencing-coverage calculator Mosdepth to identify PAVs and construct a PAV table for use in downstream comparative genome analysis.

A transposon-mediated reciprocal translocation promotes environmental adaptation but compromises domesticability of wild soybeans

Large structural variations frequently occur in higher plants; however, the impact of such variations on plant diversification, adaptation and domestication remains elusive. Here, we mapped and characterised a reciprocal chromosomal translocation in soybeans and assessed its effects on diversification and adaptation of wild (Glycine soja) and semiwild (Glycine gracilis) soybeans, and domestication of cultivated soybean (Glycine max), by tracing the distribution of the translocation in the USDA Soybean Germplasm Collection and population genetics analysis. We demonstrate that the translocation occurred through CACTA transposon-mediated chromosomal breakage in wild soybean c. 0.34 Ma and is responsible for semisterility in translocation heterozygotes and reduces their reproductive fitness. The translocation has differentiated Continental (i.e. China and Russia) populations from Maritime (i.e. Korea and Japan) populations of G. soja and predominately adapted to cold and dry climates. Further analysis revealed that the divergence of G. max from G. soja predates the translocation event and that G. gracilis is an evolutionary intermediate between G. soja and G. max. Our results highlight the effects of a chromosome rearrangement on the processes leading to plant divergence and adaptation, and provides evidence that suggests G. gracilis, rather than G. soja, as the ancestor of cultivated soybean.

Eight soybean reference genome resources from varying latitudes and agronomic traits

Comparative analysis of multiple reference genomes representing diverse genetic backgrounds is critical for understanding the role of key alleles important in domestication and genetic breeding of important crops such as soybean. To enrich the genetic resources for soybean, we describe the generation, technical assessment, and preliminary genomic variation analysis of eight de novo reference-grade soybean genome assemblies from wild and cultivated accessions. These resources represent soybeans cultured at different latitudes and exhibiting different agronomical traits. Of these eight soybeans, five are from new accessions that have not been sequenced before. We demonstrate the usage of these genomes to identify small and large genomic variations affecting known genes as well as screening for genic PAV regions for identifying candidates for further functional studies.

Analysis of genome variants in dwarf soybean lines obtained in F6 derived from cross of normal parents (cultivated and wild soybean)

Plant height is an important component of plant architecture and significantly affects crop breeding practices and yield. We studied DNA variations derived from F5 recombinant inbred lines (RILs) with 96.8% homozygous genotypes. Here, we report DNA variations between the normal and dwarf members of four lines harvested from a single seed parent in an F6 RIL population derived from a cross between Glycine max var. Peking and Glycine soja IT182936. Whole genome sequencing was carried out, and the DNA variations in the whole genome were compared between the normal and dwarf samples. We found a large number of DNA variations in both the dwarf and semi-dwarf lines, with one single nucleotide polymorphism (SNP) per at least 3.68 kb in the dwarf lines and 1 SNP per 11.13 kb of the whole genome. This value is 2.18 times higher than the expected DNA variation in the F6 population. A total of 186 SNPs and 241 SNPs were discovered in the coding regions of the dwarf lines 1282 and 1303, respectively, and we discovered 33 homogeneous nonsynonymous SNPs that occurred at the same loci in each set of dwarf and normal soybean. Of them, five SNPs were in the same positions between lines 1282 and 1303. Our results provide important information for improving our understanding of the genetics of soybean plant height and crop breeding. These polymorphisms could be useful genetic resources for plant breeders, geneticists, and biologists for future molecular biology and breeding projects.

Validation of a Triplex Quantitative Polymerase Chain Reaction Assay for Detection and Quantification of Traditional Protein Sources, Pisum sativum L. and Glycine max (L.) Merr., in Protein Powder Mixtures

Several botanicals have been traditionally used as protein sources, including the leguminous Pisum sativum L. and Glycine max (L.) Merr. While a rich history exists of cultivating these plants for their whole, protein-rich grain, modern use as powdered supplements present a new challenge in material authentication. The absence of clear morphological identifiers of an intact plant and the existence of long, complex supply chains behoove industry to create quick, reliable analytical tools to identify the botanical source of a protein product (many of which contain multiple sources). The utility of molecular tools for plant-based protein powder authentication is gaining traction, but few validated tools exist. Multiplex quantitative polymerase chain reaction (qPCR) can provide an economical means by which sources can be identified and relative proportions quantified. We followed established guidelines for the design, optimization, and validation of qPCR assay, and developed a triplex qPCR assay that can amplify and quantify pea and soy DNA targets, normalized by a calibrator. The assay was evaluated for analytical specificity, analytical sensitivity, efficiency, precision, dynamic range, repeatability, and reproducibility. We tested the quantitative ability of the assay using pea and soy DNA mixtures, finding exceptional quantitative linearity for both targets - 0.9983 (p < 0.0001) for soy and 0.9915 (p < 0.0001) for pea. Ratios based on mass of protein powder were also tested, resulting in non-linear patterns in data that suggested the requirement of further sample preparation optimization or algorithmic correction. Variation in fragment size within different lots of commercial protein powder samples was also analyzed, revealing low SD among lots. Ultimately, this study demonstrated the utility of qPCR in the context of protein powder mixtures and highlighted key considerations to take into account for commercial implementation.

GmCCD4 controls carotenoid content in soybeans

To better understand the mechanisms regulating plant carotenoid metabolism in staple crop, we report the map-based cloning and functional characterization of the Glycine max carotenoid cleavage dioxygenase 4 (GmCCD4) gene, which encodes a carotenoid cleavage dioxygenase enzyme involved in metabolizing carotenoids into volatile β-ionone. Loss of GmCCD4 protein function in four Glycine max increased carotenoid content (gmicc) mutants resulted in yellow flowers due to excessive accumulation of carotenoids in flower petals. The carotenoid contents also increase three times in gmicc1 seeds. A genome-wide association study indicated that the GmCCD4 locus was one major locus associated with carotenoid content in natural population. Further analysis indicated that the haplotype-1 of GmCCD4 gene was positively associated with higher carotenoid levels in soybean cultivars and accumulated more β-carotene in engineered E. coli with ectopic expression of different GmCCD4 haplotypes. These observations uncovered that GmCCD4 was a negative regulator of carotenoid content in soybean, and its various haplotypes provide useful resources for future soybean breeding practice.

Genomic resources in plant breeding for sustainable agriculture

Climate change during the last 40 years has had a serious impact on agriculture and threatens global food and nutritional security. From over half a million plant species, cereals and legumes are the most important for food and nutritional security. Although systematic plant breeding has a relatively short history, conventional breeding coupled with advances in technology and crop management strategies has increased crop yields by 56 % globally between 1965-85, referred to as the Green Revolution. Nevertheless, increased demand for food, feed, fiber, and fuel necessitates the need to break existing yield barriers in many crop plants. In the first decade of the 21st century we witnessed rapid discovery, transformative technological development and declining costs of genomics technologies. In the second decade, the field turned towards making sense of the vast amount of genomic information and subsequently moved towards accurately predicting gene-to-phenotype associations and tailoring plants for climate resilience and global food security. In this review we focus on genomic resources, genome and germplasm sequencing, sequencing-based trait mapping, and genomics-assisted breeding approaches aimed at developing biotic stress resistant, abiotic stress tolerant and high nutrition varieties in six major cereals (rice, maize, wheat, barley, sorghum and pearl millet), and six major legumes (soybean, groundnut, cowpea, common bean, chickpea and pigeonpea). We further provide a perspective and way forward to use genomic breeding approaches including marker-assisted selection, marker-assisted backcrossing, haplotype based breeding and genomic prediction approaches coupled with machine learning and artificial intelligence, to speed breeding approaches. The overall goal is to accelerate genetic gains and deliver climate resilient and high nutrition crop varieties for sustainable agriculture.

Soybean TILLING-by-Sequencing+ reveals the role of novel GmSACPD members in unsaturated fatty acid biosynthesis while maintaining healthy nodules

Developing soybean lines with high levels of stearic acid is a primary goal of the soybean industry. Most high-stearic-acid soybeans carry different GmSACPD-C mutated alleles. However, due to the dual role of GmSACPD-C in seeds and nodule development, all derived deleterious GmSACPD-C mutant alleles are of extremely poor agronomic value because of defective nodulation. The soybean stearoyl-acyl carrier protein desaturase (GmSACPD) gene family is composed of five members. Comparative genomics analysis indicated that SACPD genes were duplicated and derived from a common ancestor that is still present in chlorophytic algae. Synteny analysis showed the presence of segment duplications between GmSACPD-A/GmSACPD-B, and GmSACPD-C/GmSACPD-D. GmSACPD-E was not contained in any duplicated segment and may be the result of tandem duplication. We developed a TILLING by Target Capture Sequencing (Tilling-by-Sequencing+) technology, a versatile extension of the conventional TILLING by sequencing, and successfully identified 12, 14, and 18 ethyl methanesulfonate mutants at the GmSACPD-A, GmSACPD-B, and GmSACPD-D genes, respectively. Functional analysis of all identified mutants revealed an unprecedented role of GmSACPD-A, GmSACPD-B, and GmSACPD-D in unsaturated fatty acid biosynthesis without affecting nodule development and structure. This discovery will positively impact the development of high-stearic-acid lines to enhance soybean nutritional value without potential developmental tradeoffs.

The Effects of Gene Duplication Modes on the Evolution of Regulatory Divergence in Wild and Cultivated Soybean

Regulatory changes include divergence in both cis-elements and trans-factors, which play roles in organismal evolution. Whole genome duplications (WGD) followed by diploidization are a recurrent feature in the evolutionary history of angiosperms. Prior studies have shown that duplicated genes have different evolutionary fates due to variable selection constraints and results in genomic compositions with hallmarks of paleopolyploidy. The recent sequential WGDs and post-WGD evolution in the common ancestor of cultivated soybean (Glycine max) and wild soybean (Glycine soja), together with other models of gene duplication, have resulted in a highly duplicated genome. In this study, we investigated the transcriptional changes in G. soja and G. max. We identified a sizable proportion of interspecific differentially expressed genes (DEGs) and found parental expression level dominance of G. max in their F1 hybrids. By classifying genes into different regulatory divergence types, we found the trans-regulatory changes played a predominant role in transcriptional divergence between wild and cultivated soybean. The same gene ontology (GO) and protein family (Pfam) terms were found to be over-represented in DEGs and genes of cis-only between JY47 and GS, suggesting the substantial contribution of cis-regulatory divergences to the evolution of wild and cultivated soybeans. By further dissecting genes into five different duplication modes, we found genes in different duplication modes tend to accumulate different types of regulatory differences. A relatively higher proportion of cis-only regulatory divergences was detected in singleton, dispersed, proximal, and tandem duplicates than WGD duplicates and genome-wide level, which is in line with the prediction of gene balance hypothesis for the differential fates of duplicated genes post-WGD. The numbers of cis-only and trans-only regulated genes were similar for singletons, whereas there were more genes of trans-only than cis-only in the rest duplication types, especially in WGD in which there were two times more trans-only genes than that in cis-only type. Tandem duplicates showed the highest proportion of trans-only genes probably due to some special features of this class. In summary, our results demonstrate that genes in different duplication modes have different fates in transcriptional evolution underpinned by cis- or trans-regulatory divergences in soybean and likely in other paleopolyploid higher organisms.

Genome-wide identification and expression analyses of nitrate transporter family genes in wild soybean (Glycine soja)

Nitrate transporters (NRTs) are important channel proteins facilitating cross-membrane movement of small molecules like NO3- which is a critical nutrient for all life. However, the classification and evolution of nitrate transporters in the legume plants are still elusive. In this study, we surveyed the wild soybean (G. soja) genomic databases and identified 120 GsNRT1 and 5 GsNRT2 encoding genes. Phylogenetic analyses show that GsNRT1 subfamily is consisted of eight clades (NPF1 to NPF8), while GsNRT2 subfamily has only one clade. Gene chromosomal location and evolutionary historic analyses indicate that GsNRT genes are unevenly distributed on 19 out of 20 G. soja chromosomes and segmental duplications may take a major part in the expansion of GsNRT family. Investigations of gene structure and protein motif compositions suggest that GsNRT family members are highly conserved in structures of both gene and protein levels. In addition, we analyzed the spatial expression patterns of representative GsNRT genes and their responses to exogenous nitrogen and carbon supplies and different abiotic stresses. The qRT-PCR data indicated that 16 selected GsNRT genes showed various expression levels in the roots, stems, leaves, and pods of young G. soja plants, and these genes were regulated by not only nitrogen and carbohydrate nutrients but also NaCl, NaHCO3, abscisic acid (ABA), and salicylic acid (SA). These results suggest that GsNRT genes may be involved in the regulation of plant growth, development, and adaptation to environmental stresses, and the study will shed light on functional dissection of plant nitrate transporter proteins in the future.

Soybean BARCSoySNP6K: An assay for soybean genetics and breeding research

The limited number of recombinant events in recombinant inbred lines suggests that for a biparental population with a limited number of recombinant inbred lines, it is unnecessary to genotype the lines with many markers. For genomic prediction and selection, previous studies have demonstrated that only 1000-2000 genome-wide common markers across all lines/accessions are needed to reach maximum efficiency of genomic prediction in populations. Evaluation of too many markers will not only increase the cost but also generate redundant information. We developed a soybean (Glycine max) assay, BARCSoySNP6K, containing 6000 markers, which were carefully chosen from the SoySNP50K assay based on their position in the soybean genome and haplotype block, polymorphism among accessions and genotyping quality. The assay includes 5000 single nucleotide polymorphisms (SNPs) from euchromatic and 1000 from heterochromatic regions. The percentage of SNPs with minor allele frequency >0.10 was 95% and 91% in the euchromatic and heterochromatic regions, respectively. Analysis of progeny from two large families genotyped with SoySNP50K versus BARCSoySNP6K showed that the position of the common markers and number of unique bins along linkage maps were consistent based on the SNPs genotyped with the two assays; however, the rate of redundant markers was dramatically reduced with the BARCSoySNP6K. The BARCSoySNP6K assay is proven as an excellent tool for detecting quantitative trait loci, genomic selection and assessing genetic relationships. The assay is commercialized by Illumina Inc. and being used by soybean breeders and geneticists and the list of SNPs in the assay is an ideal resource for SNP genotyping by targeted amplicon sequencing.

Progresses, Challenges, and Prospects of Genome Editing in Soybean (Glycine max)

Soybean is grown worldwide for oil and protein source as food, feed and industrial raw material for biofuel. Steady increase in soybean production in the past century mainly attributes to genetic mediation including hybridization, mutagenesis and transgenesis. However, genetic resource limitation and intricate social issues in use of transgenic technology impede soybean improvement to meet rapid increases in global demand for soybean products. New approaches in genomics and development of site-specific nucleases (SSNs) based genome editing technologies have expanded soybean genetic variations in its germplasm and have potential to make precise modification of genes controlling the important agronomic traits in an elite background. ZFNs, TALENS and CRISPR/Cas9 have been adapted in soybean improvement for targeted deletions, additions, replacements and corrections in the genome. The availability of reference genome assembly and genomic resources increases feasibility in using current genome editing technologies and their new development. This review summarizes the status of genome editing in soybean improvement and future directions in this field.

Connecting signals and benefits through partner choice in plant-microbe interactions

Stabilizing mechanisms in plant-microbe symbioses are critical to maintaining beneficial functions, with two main classes: host sanctions and partner choice. Sanctions are currently presumed to be more effective and widespread, based on the idea that microbes rapidly evolve cheating while retaining signals matching cooperative strains. However, hosts that effectively discriminate among a pool of compatible symbionts would gain a significant fitness advantage. Using the well-characterized legume-rhizobium symbiosis as a model, we evaluate the evidence for partner choice in the context of the growing field of genomics. Empirical studies that rely upon bacteria varying only in nitrogen-fixation ability ignore host-symbiont signaling and frequently conclude that partner choice is not a robust stabilizing mechanism. Here, we argue that partner choice is an overlooked mechanism of mutualism stability and emphasize that plants need not use the microbial services provided a priori to discriminate among suitable partners. Additionally, we present a model that shows that partner choice signaling increases symbiont and host fitness in the absence of sanctions. Finally, we call for a renewed focus on elucidating the signaling mechanisms that are critical to partner choice while further aiming to understand their evolutionary dynamics in nature.

Genomic diversity analysis and identification of novel SSR markers in four tobacco varieties by high-throughput resequencing

Genome resequencing was carried out on two varieties of flue-cured tobacco (LY1306 and Qinyan 96), one variety of sun-cured tobacco (Wanmao 3), and one variety of air-cured Maryland tobacco (Wufeng 1), for a comparative analysis of genomic variation across the four varieties. Single nucleotide polymorphisms (SNPs), insertions and deletions (InDels), structural variations (SVs), and copy-number variations (CNVs) were then identified in each tobacco variety. Furthermore, a functional analysis of mutated genes was carried out. Through in-depth comparative analysis of genomes of different tobacco varieties, we identified genome variations in a number of SNPs, InDels, SVs, and CNVs, respectively. Computational analysis to predict the function of mutated genes containing these differential SNPs, InDels, SVs, and CNVs showed that they were mainly involved in different functions, such as carbohydrate metabolism and secondary metabolites biosynthesis. We mainly focused on genes that were involved in the biosynthesis of secondary metabolites and nicotine metabolism. In addition, we identified five simple sequence repeat (SSR)-based markers and verified them by PCR amplification in 10 tobacco varieties. Taken together, our study increases the understanding of genetic differences between tobacco types or varieties and identifies five SSR markers to classify tobacco varieties or types.

Impacts of genomic research on soybean improvement in East Asia

It has been commonly accepted that soybean domestication originated in East Asia. Although East Asia has the historical merit in soybean production, the USA has become the top soybean producer in the world since 1950s. Following that, Brazil and Argentina have been the major soybean producers since 1970s and 1990s, respectively. China has once been the exporter of soybean to Japan before 1990s, yet she became a net soybean importer as Japan and the Republic of Korea do. Furthermore, the soybean yield per unit area in East Asia has stagnated during the past decade. To improve soybean production and enhance food security in these East Asian countries, much investment has been made, especially in the breeding of better performing soybean germplasms. As a result, China, Japan, and the Republic of Korea have become three important centers for soybean genomic research. With new technologies, the rate and precision of the identification of important genomic loci associated with desired traits from germplasm collections or mutants have increased significantly. Genome editing on soybean is also becoming more established. The year 2019 marked a new era for crop genome editing in the commercialization of the first genome-edited plant product, which is a high-oleic-acid soybean oil. In this review, we have summarized the latest developments in soybean breeding technologies and the remarkable progress in soybean breeding-related research in China, Japan, and the Republic of Korea.

Genetic Architecture of Early Vigor Traits in Wild Soybean

A worldwide food shortage has been projected as a result of the current increase in global population and climate change. In order to provide sufficient food to feed more people, we must develop crops that can produce higher yields. Plant early vigor traits, early growth rate (EGR), early plant height (EPH), inter-node length, and node count are important traits that are related to crop yield. Glycine soja, the wild counterpart to cultivated soybean, Glycine max, harbors much higher genetic diversity and can grow in diverse environments. It can also cross easily with cultivated soybean. Thus, it holds a great potential in developing soybean cultivars with beneficial agronomic traits. In this study, we used 225 wild soybean accessions originally from diverse environments across its geographic distribution in East Asia. We quantified the natural variation of several early vigor traits, investigated the relationships among them, and dissected the genetic basis of these traits by applying a Genome-Wide Association Study (GWAS) with genome-wide single nucleotide polymorphism (SNP) data. Our results showed positive correlation between all early vigor traits studied. A total of 12 SNPs significantly associated with EPH were identified with 4 shared with EGR. We also identified two candidate genes, Glyma.07G055800.1 and Glyma.07G055900.1, playing important roles in influencing trait variation in both EGR and EPH in G. soja.

The Litsea genome and the evolution of the laurel family

The laurel family within the Magnoliids has attracted attentions owing to its scents, variable inflorescences, and controversial phylogenetic position. Here, we present a chromosome-level assembly of the Litsea cubeba genome, together with low-coverage genomic and transcriptomic data for many other Lauraceae. Phylogenomic analyses show phylogenetic discordance at the position of Magnoliids, suggesting incomplete lineage sorting during the divergence of monocots, eudicots, and Magnoliids. An ancient whole-genome duplication (WGD) event occurred just before the divergence of Laurales and Magnoliales; subsequently, independent WGDs occurred almost simultaneously in the three Lauralean lineages. The phylogenetic relationships within Lauraceae correspond to the divergence of inflorescences, as evidenced by the phylogeny of FUWA, a conserved gene involved in determining panicle architecture in Lauraceae. Monoterpene synthases responsible for production of specific volatile compounds in Lauraceae are functionally verified. Our work sheds light on the evolution of the Lauraceae, the genetic basis for floral evolution and specific scents.

Litsea cubeba belongs to the Lauraceae family within the Magnoliids clade. Here, the authors assemble its genome and reveal divergence of inflorescence and sexual differentiation, the phylogenetic relationships across the Lauraceae and related species, and biosynthetic genes related to essential oil synthesis.

Updates on Legume Genome Sequencing

Legumes (Fabaceae) are agronomically and economically one of the most important crops. Because legumes serve as a source of food, feed, and industrial materials, many studies in the field of legume genomics, including genome sequencing, have been conducted over the last decade. Here, we update the progress in genome sequencing of legume crops, including soybean (Glycine max [L.] Merr.), mung bean (V. radiata var. radiata), adzuki bean (V. angularis var. angularis), common bean (Phaseolus vulgaris L.), pigeon pea (Cajanus cajan), chickpea (Cicer arietinum), and peanut (Arachis hypogaea). Since the publication of the first reference genome sequence of each species, many accessions have been resequenced to study genetic diversity, speciation, and polyploidization in the legume lineage.

Identification of quantitative trait loci for increased α-tocopherol biosynthesis in wild soybean using a high-density genetic map

BACKGROUND

Soybean is one of the most important crop sources of tocopherols (Toc). However, the content of α-Toc, an isoform with the highest vitamin E activity in humans, is low in most cultivars. With the aim of broadening genetic variability, we performed quantitative trait locus (QTL) analysis for a high seed α-Toc trait detected in a wild soybean and characterized the sequence polymorphisms and expression profiles of γ-tocopherol methyltransferase (γ-TMT) genes as potential candidates.

RESULTS

A recombinant inbred line population was developed from a cross between the low α-Toc breeding line TK780 and the high α-Toc wild accession B04009. The α-Toc content in seeds correlated strongly with the ratio of α-Toc to γ-Toc contents. QTL analysis using a high-density map constructed with 7710 single nucleotide polymorphisms (SNPs) generated by restriction site-associated DNA sequencing detected six QTLs involved in α-Toc biosynthesis. Of these, three in chromosomes (Chr) 9, 11, and 12 produced consistent effects during a 2-year trial. B04009 allele at QTLs in Chr9 and Chr12 and TK780 allele at the QTL in Chr11 each promoted the conversion of γ-Toc to α-Toc, which elevated the seed α-Toc content. SNPs and indels were detected between the parents in three γ-TMT genes (γ-TMT1, γ-TMT2, and γ-TMT3) co-located in the QTLs in Chr9 and Chr12, of which some existed in the cis-regulatory elements associated with seed development and functions. In immature cotyledons, γ-TMT3 was expressed at higher levels in B04009 than TK780, irrespective of two thermal conditions tested, whereas the expression of γ-TMT2 was markedly upregulated under higher temperatures, particularly in B04009.

CONCLUSIONS

We identified QTLs consistently controlling α-Toc biosynthesis in wild soybean seeds in 2-year trials. The QTL on Chr9 had been previously identified in soybean, whereas the QTLs on Chr11 and Chr12 were novel. Further molecular dissections and characterization of the QTLs may facilitate the use of high α-Toc alleles from wild soybean in soybean breeding and an understanding of the molecular mechanisms underlying α-Toc biosynthesis in soybean seeds.

Analysis of Soybean Somatic Embryogenesis Using Chromosome Segment Substitution Lines and Transcriptome Sequencing

Soybean is an important cash crop that is widely used as a source of vegetable protein and edible oil. The regeneration ability of soybean directly affects the application of biotechnology. In this study, we used the exogenous hormone 2,4-D to treat immature embryos. Different levels of somatic incidence were selected from the chromosome segment substitution lines (CSSLs) constructed by SN14 and ZYD00006. Transcriptome sequencing of extreme materials was performed, and 2666 differentially expressed genes were obtained. At the same time, a difference table was generated by combining the data on CSSL rearrangement. In the extreme materials, a total of 93 differentially expressed genes were predicted and were then analyzed by cluster analysis and Gene Ontology (GO) annotation. After screening and annotating the target genes, three differentially expressed genes with hormone pathways were identified. The expression patterns of the target genes were verified by real-time quantitative PCR (qRT-PCR). Haplotype polymorphism detection and linkage disequilibrium analysis were performed on the candidate gene Glyma.09g248200. This study provided more information on the regulation network of soybean somatic embryogenesis and regeneration processes, and further identified important genes in the soybean regeneration process and provided a theoretical basis for accelerating the application of biotechnology to soybean for improving its breeding efficiency.

Symbiotic incompatibility between soybean and Bradyrhizobium arises from one amino acid determinant in soybean Rj2 protein

Cultivated soybean (Glycine max) carrying the Rj2 allele restricts nodulation with specific Bradyrhizobium strains via host immunity, mediated by rhizobial type III secretory protein NopP and the host resistance protein Rj2. Here we found that the single isoleucine residue I490 in Rj2 is required for induction of symbiotic incompatibility. Furthermore, we investigated the geographical distribution of the Rj2-genotype soybean in a large set of germplasm by single nucleotide polymorphism (SNP) genotyping using a SNP marker for I490. By allelic comparison of 79 accessions in the Japanese soybean mini-core collection, we suggest substitution of a single amino acid residue (R490 to I490) in Rj2 induces symbiotic incompatibility with Bradyrhizobium diazoefficiens USDA 122. The importance of I490 was verified by complementation of rj2-soybean by the dominant allele encoding the Rj2 protein containing I490 residue. The Rj2 allele was also found in Glycine soja, the wild progenitor of G. max, and their single amino acid polymorphisms were associated with the Rj2-nodulation phenotype. By SNP genotyping against 1583 soybean accessions, we detected the Rj2-genotype in 5.4% of G. max and 7.7% of G. soja accessions. Distribution of the Rj2-genotype soybean plants was relatively concentrated in the temperate Asian region. These results provide important information about the mechanism of host genotype-specific symbiotic incompatibility mediated by host immunity and suggest that the Rj2 gene has been maintained by environmental conditions during the process of soybean domestication.

Utilization of Interspecific High-Density Genetic Map of RIL Population for the QTL Detection and Candidate Gene Mining for 100-Seed Weight in Soybean

Seed-weight is one of the most important traits determining soybean yield. Hence, it is prerequisite to have detailed understanding of the genetic basis regulating seed-weight for the development of improved cultivars. In this regard, the present study used high-density interspecific linkage map of NJIR4P recombinant inbred population evaluated in four different environments to detect stable Quantitative trait loci (QTLs) as well as mine candidate genes for 100-seed weight. In total, 19 QTLs distributed on 12 chromosomes were identified in all individual environments plus combined environment, out of which seven were novel and eight are stable identified in more than one environment. However, all the novel QTLs were minor (R 2 < 10%). The remaining 12 QTLs detected in this study were co-localized with the earlier reported QTLs with narrow genomic regions, and out of these only 2 QTLs were major (R 2 > 10%) viz., qSW-17-1 and qSW-17-4. Beneficial alleles of all identified QTLs were derived from cultivated soybean parent (Nannong493-1). Based on Protein ANalysis THrough Evolutionary Relationships, gene annotation information, and literature search, 29 genes within 5 stable QTLs were predicted to be possible candidate genes that might regulate seed-weight/size in soybean. However, it needs further validation to confirm their role in seed development. In conclusion, the present study provides better understanding of trait genetics and candidate gene information through the use high-density inter-specific bin map, and also revealed considerable scope for genetic improvement of 100-seed weight in soybean using marker-assisted breeding.