Genetic Control and Geo-Climate Adaptation of Pod Dehiscence Provide Novel Insights into Soybean Domestication

Identification of quantitative trait loci of pod dehiscence in a collection of soybean grown in the southeast of Kazakhstan

Soybean [Glycine max (L.) Merr.] is one of the important crops that are constantly increasing their cultivation area in Kazakhstan. It is particularly significant in the southeastern regions of the country, which are currently predominant areas for cultivating this crop. One negative trait reducing yield in these dry areas is pod dehiscence (PD). Therefore, it is essential to understand the genetic control of PD to breed new cultivars with high yield potential. In this study, we evaluated 273 soybean accessions from different regions of the world for PD resistance in the conditions of southeastern regions of Kazakhstan in 2019 and 2021. The field data for PD suggested that 12 accessions were susceptible to PD in both studied years, and 32 accessions, in one of the two studied years. The genotyping of the collection using a DNA marker for the Pdh1 gene, a major gene for PD, revealed that 244 accessions had the homozygous R (resistant) allele, 14 had the homozygous S (susceptible) allele, and 15 accessions showed heterozygosity. To identify additional quantitative trait loci (QTLs), we applied an association mapping study using a 6K SNP Illumina iSelect array. The results suggested that in addition to major QTL on chromosome 16, linked to the physical location of Pdh1, two minor QTLs were identified on chromosomes 10 and 13. Both minor QTLs for PD were associated with calmodulin-binding protein, which presumably plays an important role in regulating PD in dry areas. Thus, the current study provided additional insight into PD regulation in soybean. The identified QTLs for PD can be efficiently employed in breeding for high-yield soybean cultivars.

Artificial selection of mutations in two nearby genes gave rise to shattering resistance in soybean

Resistance to pod shattering is a key domestication-related trait selected for seed production in many crops. Here, we show that the transition from shattering in wild soybeans to shattering resistance in cultivated soybeans resulted from selection of mutations within the coding sequences of two nearby genes - Sh1 and Pdh1. Sh1 encodes a C2H2-like zinc finger transcription factor that promotes shattering by repressing SHAT1-5 expression, thereby reducing the secondary wall thickness of fiber cap cells in the abscission layers of pod sutures, while Pdh1 encodes a dirigent protein that orchestrates asymmetric lignin distribution in inner sclerenchyma, creating torsion in pod walls that facilitates shattering. Integration analyses of quantitative trait locus mapping, genome-wide association studies, and allele distribution in representative soybean germplasm suggest that these two genes are primary modulators underlying this domestication trait. Our study thus provides comprehensive understanding regarding the genetic, molecular, and cellular bases of shattering resistance in soybeans.

QTL mapping for pod quality and yield traits in snap bean (Phaseolus vulgaris L.)

Pod quality and yield traits in snap bean (Phaseolus vulgaris L.) influence consumer preferences, crop adoption by farmers, and the ability of the product to be commercially competitive locally and globally. The objective of the study was to identify the quantitative trait loci (QTL) for pod quality and yield traits in a snap × dry bean recombinant inbred line (RIL) population. A total of 184 F6 RILs derived from a cross between Vanilla (snap bean) and MCM5001 (dry bean) were grown in three field sites in Kenya and one greenhouse environment in Davis, CA, USA. They were genotyped at 5,951 single nucleotide polymorphisms (SNPs), and composite interval mapping was conducted to identify QTL for 16 pod quality and yield traits, including pod wall fiber, pod string, pod size, and harvest metrics. A combined total of 44 QTL were identified in field and greenhouse trials. The QTL for pod quality were identified on chromosomes Pv01, Pv02, Pv03, Pv04, Pv06, and Pv07, and for pod yield were identified on Pv08. Co-localization of QTL was observed for pod quality and yield traits. Some identified QTL overlapped with previously mapped QTL for pod quality and yield traits, with several others identified as novel. The identified QTL can be used in future marker-assisted selection in snap bean.

Adaptation and the Geographic Spread of Crop Species

Crops are plant species that were domesticated starting about 11,000 years ago from several centers of origin, most prominently the Fertile Crescent, East Asia, and Mesoamerica. From their domestication centers, these crops spread across the globe and had to adapt to differing environments as a result of this dispersal. We discuss broad patterns of crop spread, including the early diffusion of crops associated with the rise and spread of agriculture, the later movement via ancient trading networks, and the exchange between the Old and New Worlds over the last ∼550 years after the European colonization of the Americas. We also examine the various genetic mechanisms associated with the evolutionary adaptation of crops to their new environments after dispersal, most prominently seasonal adaptation associated with movement across latitudes, as well as altitudinal, temperature, and other environmental factors.

Genetic study for seed germination and shattering in Euphorbia lagascae in response to different seed treatments

Euphorbia lagascae Spreng is a promising emerging oilseed crop, with its seed oil accounting for approximately 50% of the seed weight. Euphorbia oil contains a significant amount of vernolic acid, comprising two-thirds of its composition, which boasts various industrial applications, including acting as a stabilizer-plasticizer and natural dye. However, this species was known to have a high degree of seed-shattering and a low germination rate, which act as two important barriers to large-scale production and exploitation. Therefore, the objective of this study is to determine the genetic control of seed germination and seed-shattering traits in order to develop a reliable pipeline that would be applicable for industries and breeders to select superior E. lagascae lines and design a robust breeding scheme in a short time at reduced labor costs. For this objective, five different wild-type genotypes of E. lagascae that demonstrated high germination potential were crossed with an ethyl methanesulfonate (EMS) mutant genotype that produces non-shattering capsules. The F2 populations from two successful crosses (A and B) were separated into three different treated groups for seed germination evaluation and to study the segregation of 200 individuals per F2 population. The three treatments were: light, gibberellic acid (GA3), and control treatment. Consequently, plants treated with approximately 250 μmol/m2/s of light showed significant improvement in germination up to 75% in cross A and 82.4 % in cross B compared with the control plants and the group treated with 0.05% GA3. According to the chi-square test results, the inheritance pattern of seed germination in response to light treatment follows a 3:1 segregation ratio between germinated and non-germinated seeds, indicating a dominant gene action in the F2 generation. The same conclusion was followed for the shattering trait in the group treated with light, which was also simply inherited as a 3:1 ratio for shattering vs. non-shattering capsules. Our results emphasize the importance and significance of light treatment in producing uniform populations through acceptable germination and shattering resistance of the mutant genotypes of E. lagascae. This is the first report of light treatment that significantly improved seed germination of E. lagascae, which may enhance efforts in the development of this new industrial crop as a feedstock for vernolic acid production.

Parallel selection of loss-of-function alleles of Pdh1 orthologous genes in warm-season legumes for pod indehiscence and plasticity is related to precipitation

Pod dehiscence facilitates seed dispersal in wild legumes but results in yield loss in cultivated legumes. The evolutionary genetics of the legume pod dehiscence trait remain largely elusive. We characterized the pod dehiscence of chromosome segment substitution lines of Glycine max crossed with Glycine soja and found that the gene underlying the predominant quantitative trait locus (QTL) of soybean pod-shattering trait was Pod dehiscence 1 (Pdh1). A few rare loss-of-function (LoF) Pdh1 alleles were identified in G. soja, while only an allele featuring a premature stop codon was selected for pod indehiscence in cultivated soybean and spread to low-precipitation regions with increased frequency. Moreover, correlated interactions among Pdh1's haplotype, gene expression, and environmental changes for the developmental plasticity of the pod dehiscence trait were revealed in G. max. We found that orthologous Pdh1 genes specifically originated in warm-season legumes and their LoF alleles were then parallel-selected during the domestication of legume crops. Our results provide insights into the convergent evolution of pod dehiscence in warm-season legumes, facilitate an understanding of the intricate interactions between genetic robustness and environmental adaptation for developmental plasticity, and guide the breeding of new legume varieties with pod indehiscence.

Domestication of azuki bean and soybean in Japan: From the insight of archeological and molecular evidence

Domestication of azuki bean and soybean has enabled them to acquire non-dormant seeds, non-shattering pods, and larger seed size. Seed remains of the Jomon period recently discovered at archeological sites in the Central Highlands of Japan (6,000-4,000 BP) suggest that the use of azuki bean and soybean and their increase in seed size began earlier in Japan than in China and Korea; molecular phylogenetic studies indicate that azuki bean and soybean originated in Japan. Recent identification of domestication genes indicate that the domestication traits of azuki bean and soybean were established by different mechanisms. Analyses of domestication related genes using DNA extracted from the seed remains would reveal further details about their domestication processes.

Legume-wide comparative analysis of pod shatter locus PDH1 reveals phaseoloid specificity, high cowpea expression, and stress responsive genomic context

Pod dehiscence is a major source of yield loss in legumes, which is exacerbated by aridity. Disruptive mutations in "Pod indehiscent 1" (PDH1), a pod sclerenchyma-specific lignin biosynthesis gene, has been linked to significant reductions in dehiscence in several legume species. We compared syntenic PDH1 regions across 12 legumes and two outgroups to uncover key historical evolutionary trends at this important locus. Our results clarified the extent to which PDH1 orthologs are present in legumes, showing the typical genomic context surrounding PDH1 has only arisen relatively recently in certain phaseoloid species (Vigna, Phaseolus, Glycine). The notable absence of PDH1 in Cajanus cajan may be a major contributor to its indehiscent phenotype compared with other phaseoloids. In addition, we identified a novel PDH1 ortholog in Vigna angularis and detected remarkable increases in PDH1 transcript abundance during Vigna unguiculata pod development. Investigation of the shared genomic context of PDH1 revealed it lies in a hotspot of transcription factors and signaling gene families that respond to abscisic acid and drought stress, which we hypothesize may be an additional factor influencing expression of PDH1 under specific environmental conditions. Our findings provide key insights into the evolutionary history of PDH1 and lay the foundation for optimizing the pod dehiscence role of PDH1 in major and understudied legume species.

Genome-wide association analysis identified consistent QTL for seed yield in a soybean diversity panel tested across multiple environments

Improving seed yield is one of the main targets of soybean [Glycine max (L.) Merr.] breeding. Identification of loci that influence productivity and understanding their genetic mechanism will help marker-assisted trait introgression. The present study evaluated a diverse panel of 541 soybean genotypes consisting of three maturity groups (MGs III-V) in four environments in Kansas, U.S. Data on seed yield, seed weight, shattering resistance, days to maturity, and plant height showed significant genotype, environmental, and genotype × environment interaction variations. Seed yield and shattering had moderate broad-sense heritability (<85%), while the rest of the traits showed high broad-sense heritability (>90%). The SoySNP50K iSelect BeadChip dataset was used to identify significantly associated loci via genome-wide association studies (GWAS). A total of 19 single-nucleotide polymorphisms (SNPs) were significantly associated with seed yield. Particularly, two stable seed yield quantitative trait loci (QTL) on chromosomes 9 and 17 were consistently detected in at least three out of four environments. Candidate gene analysis surrounding seed yield QTL on chromosome 9 showed that Glyma.09G048900, an oxygen binding protein, was the closest to the QTL peak. Similarly, Glyma.17G090200 and Glyma.17G090400 were within 20-kb region of the seed yield QTL on chromosome 17. The candidate genes warrant further analysis to determine their functional mechanisms and develop markers for seed yield improvement.

Advancing Grain Legumes Domestication and Evolution Studies with Genomics

Grain legumes were domesticated in parallel with cereals in several regions of the world and formed the economic basis of early farming cultures. Since then, legumes have played a vital role in human and animal diets and in fostering agrobiodiversity. Increasing grain legume cultivation will be crucial to safeguard nutritional security and the resilience of agricultural ecosystems across the globe. A better understanding of the molecular underpinnings of domestication and crop evolution of grain legumes may be translated into practical approaches in modern breeding programs to stabilize yield, which is threatened by evolving pathogens and changing climates. During recent decades, domestication research in all crops has greatly benefited from the fast progress in genomic technologies. Yet still, many questions surrounding the domestication and diversification of legumes remain unanswered. In this review, we assess the potential of genomic approaches in grain legume research. We describe the centers of origin and the crucial domestication traits of grain legumes. In addition, we survey the effect of domestication on both above-ground and below-ground traits that have economic importance. Finally, we discuss open questions in grain legume domestication and diversification and outline how to bridge the gap between the preservation of historic crop diversity and their utilization in modern plant breeding.

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.

Loss of pod strings in common bean is associated with gene duplication, retrotransposon insertion and overexpression of PvIND

Fruit development has been central in the evolution and domestication of flowering plants. In common bean (Phaseolus vulgaris), the principal global grain legume staple, two main production categories are distinguished by fibre deposition in pods: dry beans, with fibrous, stringy pods; and stringless snap/green beans, with reduced fibre deposition, which frequently revert to the ancestral stringy state. Here, we identify genetic and developmental patterns associated with pod fibre deposition. Transcriptional, anatomical, epigenetic and genetic regulation of pod strings were explored through RNA-seq, RT-qPCR, fluorescence microscopy, bisulfite sequencing and whole-genome sequencing. Overexpression of the INDEHISCENT ('PvIND') orthologue was observed in stringless types compared with isogenic stringy lines, associated with overspecification of weak dehiscence-zone cells throughout the pod vascular sheath. No differences in DNA methylation were correlated with this phenotype. Nonstringy varieties showed a tandemly direct duplicated PvIND and a Ty1-copia retrotransposon inserted between the two repeats. These sequence features are lost during pod reversion and are predictive of pod phenotype in diverse materials, supporting their role in PvIND overexpression and reversible string phenotype. Our results give insight into reversible gain-of-function mutations and possible genetic solutions to the reversion problem, of considerable economic value for green bean production.

Elimination of an unfavorable allele conferring pod shattering in an elite soybean cultivar by CRISPR/Cas9

Pod shattering can lead to devastating yield loss of soybean and has been a negatively selected trait in soybean domestication and breeding. Nevertheless, a significant portion of soybean cultivars are still pod shattering-susceptible, limiting their regional and climatic adaptabilities. Here we performed genetic diagnosis on the shattering-susceptible trait of a national registered cultivar, Huachun6 (HC6), and found that HC6 carries the susceptible genotype of a candidate Pod dehiscence 1 (PDH1) gene, which exists in a significant portion of soybean cultivars. We next performed genome editing on PDH1 gene by clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated protein 9 (Cas9). In T₂ progenies, several transgene-free lines with pdh1 mutations were characterized without affecting major agronomic traits. The pdh1 mutation significantly improved the pod shattering resistance which is associated with aberrant lignin distribution in inner sclerenchyma. Our work demonstrated that precision breeding by genome editing on PDH1 holds great potential for precisely improving pod shattering resistance and adaptability of soybean cultivars.

A novel Synthetic phenotype association study approach reveals the landscape of association for genomic variants and phenotypes

INTRODUCTION

Genome-Wide Association Studies (GWAS) identify tagging variants in the genome that are statistically associated with the phenotype because of their linkage disequilibrium (LD) relationship with the causative mutation (CM). When both low-density genotyped accession panels with phenotypes and resequenced data accession panels are available, tagging variants can assist with post-GWAS challenges in CM discovery.

OBJECTIVES

Our objective was to identify additional GWAS evaluation criteria to assess correspondence between genomic variants and phenotypes, as well as enable deeper analysis of the localized landscape of association.

METHODS

We used genomic variant positions as Synthetic phenotypes in GWAS that we named "Synthetic phenotype association study" (SPAS). The extreme case of SPAS is what we call an "Inverse GWAS" where we used CM positions of cloned soybean genes. We developed and validated the Accuracy concept as a measure of the correspondence between variant positions and phenotypes.

RESULTS

The SPAS approach demonstrated that the genotype status of an associated variant used as a Synthetic phenotype enabled us to explore the relationships between tagging variants and CMs, and further, that utilizing CMs as Synthetic phenotypes in Inverse GWAS illuminated the landscape of association. We implemented the Accuracy calculation for a curated accession panel to an online Accuracy calculation tool (AccuTool) as a resource for gene identification in soybean. We demonstrated our concepts on three examples of soybean cloned genes. As a result of our findings, we devised an enhanced "GWAS to Genes" analysis (Synthetic phenotype to CM strategy, SP2CM). Using SP2CM, we identified a CM for a novel gene.

CONCLUSION

The SP2CM strategy utilizing Synthetic phenotypes and the Accuracy calculation of correspondence provides crucial information to assist researchers in CM discovery. The impact of this work is a more effective evaluation of landscapes of GWAS associations.

Genetic, anatomical, and environmental patterns related to pod shattering resistance in domesticated cowpea [Vigna unguiculata (L.) Walp]

Pod shattering, which causes the explosive release of seeds from the pod, is one of the main sources of yield losses in cowpea in arid and semi-arid areas. Reduction of shattering has therefore been a primary target for selection during domestication and improvement of cowpea, among other species. Using a mini-core diversity panel of 368 cowpea accessions, four regions with a statistically significant association with pod shattering were identified. Two genes (Vigun03g321100 and Vigun11g100600), involved in cell wall biosynthesis, were identified as strong candidates for pod shattering. Microscopic analysis was conducted on a subset of accessions representing the full spectrum of shattering phenotypes. This analysis indicated that the extent of wall fiber deposition was highly correlated with shattering. The results from this study also demonstrate that pod shattering in cowpea is exacerbated by arid environmental conditions. Finally, using a subset of West African landraces, patterns of historical selection for shattering resistance related to precipitation in the environment of origin were identified. Together, these results shed light on sources of resistance to pod shattering, which will, in turn, improve climate resilience of a major global nutritional staple.

Pod shattering in grain legumes: emerging genetic and environment-related patterns

A reduction in pod shattering is one of the main components of grain legume domestication. Despite this, many domesticated legumes suffer serious yield losses due to shattering, particularly under arid conditions. Mutations related to pod shattering modify the twisting force of pod walls or the structural strength of the dehiscence zone in pod sutures. At a molecular level, a growing body of evidence indicates that these changes are controlled by a relatively small number of key genes that have been selected in parallel across grain legume species, supporting partial molecular convergence. Legume homologs of Arabidopsis thaliana silique shattering genes play only minor roles in legume pod shattering. Most domesticated grain legume species contain multiple shattering-resistance genes, with mutants of each gene typically showing only partial shattering resistance. Hence, crosses between varieties with different genes lead to transgressive segregation of shattering alleles, producing plants with either enhanced shattering resistance or atavistic susceptibility to the trait. The frequency of these resistance pod-shattering alleles is often positively correlated with environmental aridity. The continued development of pod-shattering-related functional information will be vital for breeding crops that are suited to the increasingly arid conditions expected in the coming decades.

Comprehensive genomic resources related to domestication and crop improvement traits in Lima bean

Lima bean (Phaseolus lunatus L.), one of the five domesticated Phaseolus bean crops, shows a wide range of ecological adaptations along its distribution range from Mexico to Argentina. These adaptations make it a promising crop for improving food security under predicted scenarios of climate change in Latin America and elsewhere. In this work, we combine long and short read sequencing technologies with a dense genetic map from a biparental population to obtain the chromosome-level genome assembly for Lima bean. Annotation of 28,326 gene models show high diversity among 1917 genes with conserved domains related to disease resistance. Structural comparison across 22,180 orthologs with common bean reveals high genome synteny and five large intrachromosomal rearrangements. Population genomic analyses show that wild Lima bean is organized into six clusters with mostly non-overlapping distributions and that Mesomerican landraces can be further subdivided into three subclusters. RNA-seq data reveal 4275 differentially expressed genes, which can be related to pod dehiscence and seed development. We expect the resources presented here to serve as a solid basis to achieve a comprehensive view of the degree of convergent evolution of Phaseolus species under domestication and provide tools and information for breeding for climate change resiliency.

Functional characterization of soybean (Glycine max) DIRIGENT genes reveals an important role of GmDIR27 in the regulation of pod dehiscence

DIRIGENT (DIR) genes play important roles in regulating plant growth and development and have been studied in many plant species. However, information on DIR genes in soybean is limited. Here, we identified and characterized 54 GmDIRs and studied the characteristics of GmDIRs. Most of the GmDIRs contained a classical gene structure, one exon; 26 conserved motifs were found among these GmDIRs. The GmDIRs were grouped into four subfamilies, DIR-a, DIR-b, DIR-e and DIR-f, based on a phylogenetic analysis, and 24 duplicated gene pairs were identified. Differences in the cis-acting elements in the GmDIR promoter regions might result in distinct expression patterns of GmDIRs in different tissues. In addition, GmDIR27 had a close relationship with the pod dehiscence gene GmPdh1, and overexpression of GmDIR27 increased pod dehiscence by affecting several pod dehiscence-related gene expressions. Generally, our results provide essential information that aids future efforts to functionally characterize soybean GmDIR genes.