Molecular footprints of domestication and improvement in soybean revealed by whole genome re-sequencing

Phylogenetic and functional analysis of tiller angle control homeologs in allotetraploid cotton

Introduction

Plants can adapt their growth to optimize light capture in competitive environments, with branch angle being a crucial factor influencing plant phenotype and physiology. Decreased branch angles in cereal crops have been shown to enhance productivity in high-density plantings. The Tiller Angle Control (TAC1) gene, known for regulating tiller inclination in rice and corn, has been found to control branch angle in eudicots. Manipulating TAC1 in field crops like cotton offers the potential for improving crop productivity.

Methods

Using a homolog-based methodology, we examined the distribution of TAC1-related genes in cotton compared to other angiosperms. Furthermore, tissue-specific qPCR analysis unveiled distinct expression patterns of TAC1 genes in various cotton tissues. To silence highly expressed specific TAC1 homeologs in the stem, we applied CRISPR-Cas9 gene editing and Agrobacterium-mediated transformation, followed by genotyping and subsequent phenotypic validation of the mutants.

Results

Gene duplication events of TAC1 specific to the Gossypium lineage were identified, with 3 copies in diploid progenitors and 6 copies in allotetraploid cottons. Sequence analysis of the TAC1 homeologs in Gossypium hirsutum revealed divergence from other angiosperms with 1-2 copies, suggesting possible neo- or sub-functionalization for the duplicated copies. These TAC1 homeologs exhibited distinct gene expression patterns in various tissues over developmental time, with elevated expression of A11G109300 and D11G112200, specifically in flowers and stems, respectively. CRISPR-mediated loss of these TAC1 homeologous genes resulted in a reduction in branch angle and altered petiole angles, and a 5 to 10-fold reduction in TAC1 expression in the mutants, confirming their role in controlling branch and petiole angles. This research provides a promising strategy for genetically engineering branch and petiole angles in commercial cotton varieties, potentially leading to increased productivity.

Genome-wide association study on resistance of cultivated soybean to Fusarium oxysporum root rot in Northeast China

BACKGROUND

Fusarium oxysporum is a prevalent fungal pathogen that diminishes soybean yield through seedling disease and root rot. Preventing Fusarium oxysporum root rot (FORR) damage entails on the identification of resistance genes and developing resistant cultivars. Therefore, conducting fine mapping and marker development for FORR resistance genes is of great significance for fostering the cultivation of resistant varieties. In this study, 350 soybean germplasm accessions, mainly from Northeast China, underwent genotyping using the SoySNP50K Illumina BeadChip, which includes 52,041 single nucleotide polymorphisms (SNPs). Their resistance to FORR was assessed in a greenhouse. Genome-wide association studies utilizing the general linear model, mixed linear model, compressed mixed linear model, and settlement of MLM under progressively exclusive relationship models were conducted to identify marker-trait associations while effectively controlling for population structure.

RESULTS

The results demonstrated that these models effectively managed population structure. Eight SNP loci significantly associated with FORR resistance in soybean were detected, primarily located on Chromosome 6. Notably, there was a strong linkage disequilibrium between the large-effect SNPs ss715595462 and ss715595463, contributing substantially to phenotypic variation. Within the genetic interval encompassing these loci, 28 genes were present, with one gene Glyma.06G088400 encoding a protein kinase family protein containing a leucine-rich repeat domain identified as a potential candidate gene in the reference genome of Williams82. Additionally, quantitative real-time reverse transcription polymerase chain reaction analysis evaluated the gene expression levels between highly resistant and susceptible accessions, focusing on primary root tissues collected at different time points after F. oxysporum inoculation. Among the examined genes, only this gene emerged as the strongest candidate associated with FORR resistance.

CONCLUSIONS

The identification of this candidate gene Glyma.06G088400 improves our understanding of soybean resistance to FORR and the markers strongly linked to resistance can be beneficial for molecular marker-assisted selection in breeding resistant soybean accessions against F. oxysporum.

Evaluating the Response of Glycine soja Accessions to Fungal Pathogen Macrophomina phaseolina during Seedling Growth

Charcoal rot caused by the fungal pathogen Macrophomina phaseolina (Tassi) Goid is one of various devastating soybean (Glycine max (L.) Merr.) diseases, which can severely reduce crop yield. The investigation into the genetic potential for charcoal rot resistance of wild soybean (Glycine soja) accessions will enrich our understanding of the impact of soybean domestication on disease resistance; moreover, the identified charcoal rot-resistant lines can be used to improve soybean resistance to charcoal rot. The objective of this study was to evaluate the resistance of wild soybean accessions to M. phaseolina at the seedling stage and thereby select the disease-resistant lines. The results show that the fungal pathogen infection reduced the growth of the root and hypocotyl in most G. soja accessions. The accession PI 507794 displayed the highest level of resistance response to M. phaseolina infection among the tested wild soybean accessions, while PI 487431 and PI 483660B were susceptible to charcoal rot in terms of the reduction in root and hypocotyl growth. The mean values of the root and hypocotyl parameters in PI 507794 were significantly higher (p < 0.05) than those of PI 487431 and PI 483460B. A analysis of the resistance of wild soybean accessions to M. phaseolina using the root and hypocotyl as the assessment parameters at the early seedling stage provides an alternative way to rapidly identify potential resistant genotypes and facilitate breeding for soybean resistance to charcoal rot.

Pan-centromere reveals widespread centromere repositioning of soybean genomes

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

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

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

The Key to the Future Lies in the Past: Insights from Grain Legume Domestication and Improvement Should Inform Future Breeding Strategies

Crop domestication is a co-evolutionary process that has rendered plants and animals significantly dependent on human interventions for survival and propagation. Grain legumes have played an important role in the development of Neolithic agriculture some 12,000 years ago. Despite being early companions of cereals in the origin and evolution of agriculture, the understanding of grain legume domestication has lagged behind that of cereals. Adapting plants for human use has resulted in distinct morpho-physiological changes between the wild ancestors and domesticates, and this distinction has been the focus of several studies aimed at understanding the domestication process and the genetic diversity bottlenecks created. Growing evidence from research on archeological remains, combined with genetic analysis and the geographical distribution of wild forms, has improved the resolution of the process of domestication, diversification and crop improvement. In this review, we summarize the significance of legume wild relatives as reservoirs of novel genetic variation for crop breeding programs. We describe key legume features, which evolved in response to anthropogenic activities. Here, we highlight how whole genome sequencing and incorporation of omics-level data have expanded our capacity to monitor the genetic changes accompanying these processes. Finally, we present our perspective on alternative routes centered on de novo domestication and re-domestication to impart significant agronomic advances of novel crops over existing commodities. A finely resolved domestication history of grain legumes will uncover future breeding targets to develop modern cultivars enriched with alleles that improve yield, quality and stress tolerance.

An enhancing effect attributed to a nonsynonymous mutation in SOYBEAN SEED SIZE 1, a SPINDLY-like gene, is exploited in soybean domestication and improvement

Soybean (Glycine max) was domesticated from its wild relative Glycine soja. One-hundred-seed weight is one of the most important domesticated traits determining soybean yield; however, its underlying genetic basis remains elusive. We characterized a soybean seed size 1 (sss1) mutant featuring large seeds compared to its wild-type background. Positional cloning revealed that the candidate gene GmSSS1 encoded a SPINDLY homolog and was co-located in a well-identified quantitative trait locus (QTL)-rich region on chromosome 19. Knocking out GmSSS1 resulted in small seeds, while overexpressing GmSSS1/Gmsss1 induced large seeds. Modulating GmSSS1/Gmsss1 in transgenic plants can positively influence cell expansion and cell division. Relative to GmSSS1, one mutation leading to an E to Q substitution at the 182^nd residue in Gmsss1 conferred an enhancing effect on seed weight. GmSSS1 underwent diversification in wild-type and cultivated soybean, and the alleles encoding the Gmsss1-type substitution of 182^nd -Q, which originated along the central and downstream parts of the Yellow River, were selected and expanded during soybean domestication and improvement. We cloned the causative gene for the sss1 mutant, which is linked with a seed weight QTL, identified an elite allele of this gene for increasing seed weight, and provided new insights into soybean domestication and breeding.

The genomic signature of wild‐to‐crop introgression during the domestication of scarlet runner bean ( Phaseolus coccineus L.)

The scarlet runner bean (Phaseolus coccineus) is one of the five domesticated Phaseolus species. It is cultivated in small‐scale agriculture in the highlands of Mesoamerica for its dry seeds and immature pods, and unlike the other domesticated beans, P. coccineus is an open‐pollinated legume. Contrasting with its close relative, the common bean, few studies focusing on its domestication history have been conducted. Demographic bottlenecks associated with domestication might reduce genetic diversity and facilitate the accumulation of deleterious mutations. Conversely, introgression from wild relatives could be a source of variation. Using Genotyping by Sequencing data (79,286 single‐nucleotide variants) from 237 cultivated and wild samples, we evaluated the demographic history of traditional varieties from different regions of Mexico and looked for evidence of introgression between sympatric wild and cultivated populations. Traditional varieties have high levels of diversity, even though there is evidence of a severe initial genetic bottleneck followed by a population expansion. Introgression from wild to domesticated populations was detected, which might contribute to the recovery of the genetic variation. Introgression has occurred at different times: constantly in the center of Mexico; recently in the North West; and anciently in the South. Several factors are acting together to increase and maintain genetic diversity in P. coccineus cultivars, such as demographic expansion and introgression. Wild relatives represent a valuable genetic resource and have played a key role in scarlet runner bean evolution via introgression into traditional varieties.

Current overview on the genetic basis of key genes involved in soybean domestication

Modern crops were created through the domestication and genetic introgression of wild relatives and adaptive differentiation in new environments. Identifying the domestication-related genes and unveiling their molecular diversity provide clues for understanding how the domesticated variants were selected by ancient people, elucidating how and where these crops were domesticated. Molecular genetics and genomics have explored some domestication-related genes in soybean (Glycine max). Here, we summarize recent studies about the quantitative trait locus (QTL) and genes involved in the domestication traits, introduce the functions of these genes, clarify which alleles of domesticated genes were selected during domestication. A deeper understanding of soybean domestication could help to break the bottleneck of modern breeding by highlighting unused genetic diversity not selected in the original domestication process, as well as highlighting promising new avenues for the identification and research of important agronomic traits among different crop species.

High-throughput characterization, correlation, and mapping of leaf photosynthetic and functional traits in the soybean (Glycine max) nested association mapping population

Photosynthesis is a key target to improve crop production in many species including soybean [Glycine max (L.) Merr.]. A challenge is that phenotyping photosynthetic traits by traditional approaches is slow and destructive. There is proof-of-concept for leaf hyperspectral reflectance as a rapid method to model photosynthetic traits. However, the crucial step of demonstrating that hyperspectral approaches can be used to advance understanding of the genetic architecture of photosynthetic traits is untested. To address this challenge, we used full-range (500-2,400 nm) leaf reflectance spectroscopy to build partial least squares regression models to estimate leaf traits, including the rate-limiting processes of photosynthesis, maximum Rubisco carboxylation rate, and maximum electron transport. In total, 11 models were produced from a diverse population of soybean sampled over multiple field seasons to estimate photosynthetic parameters, chlorophyll content, leaf carbon and leaf nitrogen percentage, and specific leaf area (with R2 from 0.56 to 0.96 and root mean square error approximately <10% of the range of calibration data). We explore the utility of these models by applying them to the soybean nested association mapping population, which showed variability in photosynthetic and leaf traits. Genetic mapping provided insights into the underlying genetic architecture of photosynthetic traits and potential improvement in soybean. Notably, the maximum Rubisco carboxylation rate mapped to a region of chromosome 19 containing genes encoding multiple small subunits of Rubisco. We also mapped the maximum electron transport rate to a region of chromosome 10 containing a fructose 1,6-bisphosphatase gene, encoding an important enzyme in the regeneration of ribulose 1,5-bisphosphate and the sucrose biosynthetic pathway. The estimated rate-limiting steps of photosynthesis were low or negatively correlated with yield suggesting that these traits are not influenced by the same genetic mechanisms and are not limiting yield in the soybean NAM population. Leaf carbon percentage, leaf nitrogen percentage, and specific leaf area showed strong correlations with yield and may be of interest in breeding programs as a proxy for yield. This work is among the first to use hyperspectral reflectance to model and map the genetic architecture of the rate-limiting steps of photosynthesis.

SoySNP618K array: A high-resolution single nucleotide polymorphism platform as a valuable genomic resource for soybean genetics and breeding

Innovations in genomics have enabled the development of low-cost, high-resolution, single nucleotide polymorphism (SNP) genotyping arrays that accelerate breeding progress and support basic research in crop science. Here, we developed and validated the SoySNP618K array (618,888 SNPs) for the important crop soybean. The SNPs were selected from whole-genome resequencing data containing 2,214 diverse soybean accessions; 29.34% of the SNPs mapped to genic regions representing 86.85% of the 56,044 annotated high-confidence genes. Identity-by-state analyses of 318 soybeans revealed 17 redundant accessions, highlighting the potential of the SoySNP618K array in supporting gene bank management. The patterns of population stratification and genomic regions enriched through domestication were highly consistent with previous findings based on resequencing data, suggesting that the ascertainment bias in the SoySNP618K array was largely compensated for. Genome-wide association mapping in combination with reported quantitative trait loci enabled fine-mapping of genes known to influence flowering time, E2 and GmPRR3b, and of a new candidate gene, GmVIP5. Moreover, genomic prediction of flowering and maturity time in 502 recombinant inbred lines was highly accurate (>0.65). Thus, the SoySNP618K array is a valuable genomic tool that can be used to address many questions in applied breeding, germplasm management, and basic crop research.

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.

Assessment of the Breeding Potential of a Set of Genotypes Selected from a Natural Population of Akebia trifoliata (Three–Leaf Akebia)

Akebia trifoliata (three-leaf akebia) has long been used as a medicinal herb and has the potential to be used in diverse ways, especially as a fruit crop. However, efforts to domesticate and cultivate new varieties for commercial use are only in their infancy. Here, we evaluated the genetic diversity of 29 genotypes, which were previously selected from a natural population consisting of 1447 genotypes and exhibiting high resistance to fungal diseases and a smooth peel of A. trifoliata using 85 genome-specific single sequence repeat (SSR) markers. We also characterized variation in 19 phenotypic traits and nutritional components. Large variation in phenotypic traits and nutritional components was observed, especially in vitamin C, seed/pulp, and fruit color. Correlation analyses revealed that many phenotypic traits and nutritional components were significantly correlated. A principal component analysis identified five principal components, which explained 83.2% of the total variation in the data. The results of the SSR analysis revealed that 80 of the 85 SSR markers were polymorphic; the total number of alleles amplified was 532. The expected heterozygosity was 0.672, and Shannon’s information index was 1.328. A Ward dendrogram and unweighted pair group method with arithmetic mean dendrogram revealed high diversity among the 29 genotypes and suggested that the measured morphological and nutritional traits were genetically independent of disease resistance and texture traits, as well as SSR marker loci. Finally, our results suggest that additional rounds of selection from the selected population, despite its small size, could be effective for the development of new A. trifoliata fruit cultivars.

Domestication and improvement genes reveal the differences of seed size- and oil-related traits in soybean domestication and improvement

•

Due to reduced diversity, it is essential to map domesticated and improved genes.

•

13 known and 442 candidate genes were mined for seed size- and oil-related traits.

•

All the genes were used to explain trait changes in domestication and improvement.

•

56 domesticated and 15 improved genes may be valuable for future soybean breeding.

•

This study provides useful gene resources for future breeding and biology research.

To address domestication and improvement studies of soybean seed size- and oil-related traits, a series of domesticated and improved regions, loci, and candidate genes were identified in 286 soybean accessions using domestication and improvement analyses, genome-wide association studies, quantitative trait locus (QTL) mapping and bulked segregant analyses in this study. As a result, 534 candidate domestication regions (CDRs) and 458 candidate improvement regions (CIRs) were identified in this study and integrated with those in five and three previous studies, respectively, to obtain 952 CDRs and 538 CIRs; 1469 loci for soybean seed size- and oil-related traits were identified in this study and integrated with those in Soybase to obtain 433 QTL clusters. The two results were intersected to obtain 245 domestication and 221 improvement loci for the above traits. Around these trait-related domestication and improvement loci, 7 domestication and 7 improvement genes were found to be truly associated with these traits, and 372 candidate domestication and 87 candidate improvement genes were identified using gene expression, SNP variants in genome, miRNA binding, KEGG pathway, DNA methylation, and haplotype analysis. These genes were used to explain the trait changes in domestication and improvement. As a result, the trait changes can be explained by their frequencies of elite haplotypes, base mutations in coding region, and three factors affecting their expression levels. In addition, 56 domestication and 15 improvement genes may be valuable for future soybean breeding. This study can provide useful gene resources for future soybean breeding and molecular biology research.

Assessment of the Genetic Structure and Diversity of Soybean (Glycine&nbsp;max L.) Germplasm Using Diversity Array Technology and Single Nucleotide Polymorphism Markers

Knowledge of the genetic structure and diversity of germplasm collections is crucial for sustainable genetic improvement through hybridization programs and rapid adaptation to changing breeding objectives. The objective of this study was to determine the genetic diversity and population structure of 281 International Institute of Tropical Agriculture (IITA) soybean accessions using diversity array technology (DArT) and single nucleotide polymorphism (SNP) markers for the efficient utilization of these accessions. From the results, the SNP and DArT markers were well distributed across the 20 soybean chromosomes. The cluster and principal component analyses revealed the genetic diversity among the 281 accessions by grouping them into two stratifications, a grouping that was also evident from the population structure analysis, which divided the 281 accessions into two distinct groups. The analysis of molecular variance revealed that 97% and 98% of the genetic variances using SNP and DArT markers, respectively, were within the population. Genetic diversity indices such as Shannon's diversity index, diversity and unbiased diversity revealed the diversity among the different populations of the soybean accessions. The SNP and DArT markers used provided similar information on the structure, diversity and polymorphism of the accessions, which indicates the applicability of the DArT marker in genetic diversity studies. Our study provides information about the genetic structure and diversity of the IITA soybean accessions that will allow for the efficient utilization of these accessions in soybean improvement programs, especially in Africa.

The CCCH zinc finger family of soybean (Glycine max L.): genome-wide identification, expression, domestication, GWAS and haplotype analysis

BACKGROUND

The CCCH zinc finger (zf_CCCH) is a unique subfamily featured one or more zinc finger motif(s) comprising of three Cys and one His residues. The zf_CCCH family have been reported involving in various processes of plant development and adaptation.

RESULTS

In this study, the zf_CCCH genes were identified via a genome-wide search and were systematically analyzed. 116 Gmzf_CCCHs were obtained and classified into seventeen subfamilies. Gene duplication and expansion analysis showed that tandem and segmental duplications contributed to the expansion of the Gmzf_CCCH gene family, and that segmental duplication play the main role. The expression patterns of Gmzf_CCCH genes were tissue-specific. Eleven domesticated genes were detected involved in the regulation of seed oil and protein synthesis as well as growth and development of soybean through GWAS and haplotype analysis for Gmzf_CCCH genes among the 164 of 302 soybeans resequencing data. Among which, 8 genes play an important role in the synthesis of seed oil or fatty acid, and the frequency of their elite haplotypes changes significantly among wild, landrace and improved cultivars, indicating that they have been strongly selected in the process of soybean domestication.

CONCLUSIONS

This study provides a scientific foundation for the comprehensive understanding, future cloning and functional studies of Gmzf_CCCH genes in soybean, meanwhile, it was also helpful for the improvement of soybean with high oil content.

Identifying genomic regions targeted during eggplant domestication using transcriptome data

Identifying genes and traits that have diverged during domestication provides key information of importance for maintaining and even increasing yield and nutrients in existing crops. A “bottom-up” population genetics approach was used to identify signatures of selection across the eggplant genome, to better understand the process of domestication. RNA-seq data were obtained for 4 wild eggplants (Solanum insanum L.) and 16 domesticated eggplants (S. melongena L.) and mapped to the eggplant genome. Single-nucleotide polymorphism (SNPs) exhibiting signatures of selection in domesticates were identified as those exhibiting high F_ST between the 2 populations (evidence of significant divergence) and low π for the domesticated population (indicative of a selective sweep). Some of these regions appear to overlap with previously identified quantitative trait loci for domestication traits. Genes in regions of linkage disequilibrium surrounding these SNPs were searched against the Arabidopsis thaliana and tomato genomes to find orthologs. Subsequent gene ontology (GO) enrichment analysis identified over-representation of GO terms related to photosynthesis and response to the environment. This work reveals genomic changes involved in eggplant domestication and improvement, and how this compares to observed changes in the tomato genome, revealing shared chromosomal regions involved in the domestication of both species.

Multi-omic dissection of the drought resistance traits of soybean landrace LX

With diverse genetic backgrounds, soybean landraces are valuable resource for breeding programs. Herein, we apply multi-omic approaches to extensively characterize the molecular basis of drought tolerance in the soybean landrace LX. Initial screens established that LX performed better with PEG6000 treatment than control cultivars. LX germinated better than William 82 under drought conditions and accumulated more anthocyanin and flavonoids. Untargeted mass spectrometry in combination with transcriptomic analyses revealed the chemical diversity and genetic basis underlying the overall performance of LX landrace. Under control and drought conditions, significant differences in the expression of a suite of secondary metabolism genes, particularly those involved in the general phenylpropanoid pathway and flavonoid but not lignin biosynthesis, were seen in LX and William 82. The expression of these genes correlated with the corresponding metabolites in LX plants. Further correlation analysis between metabolites and transcripts identified pathway structural genes and transcription factors likely are responsible for the LX agronomic traits. The activities of some key biosynthetic genes or regulators were confirmed through heterologous expression in transgenic Arabidopsis and hairy root transformation in soybean. We propose a regulatory mechanism based on flavonoid secondary metabolism and adaptive traits of this landrace which could be of relevance to cultivated soybean.

Transcriptome Analyses Throughout Chili Pepper Fruit Development Reveal Novel Insights into the Domestication Process

Chili pepper (Capsicum spp.) is an important crop, as well as a model for fruit development studies and domestication. Here, we performed a time-course experiment to estimate standardized gene expression profiles with respect to fruit development for six domesticated and four wild chili pepper ancestors. We sampled the transcriptomes every 10 days from flowering to fruit maturity, and found that the mean standardized expression profiles for domesticated and wild accessions significantly differed. The mean standardized expression was higher and peaked earlier for domesticated vs. wild genotypes, particularly for genes involved in the cell cycle that ultimately control fruit size. We postulate that these gene expression changes are driven by selection pressures during domestication and show a robust network of cell cycle genes with a time shift in expression, which explains some of the differences between domesticated and wild phenotypes.

Whole-genome sequence diversity and association analysis of 198 soybean accessions in mini-core collections

We performed whole-genome Illumina resequencing of 198 accessions to examine the genetic diversity and facilitate the use of soybean genetic resources and identified 10 million single nucleotide polymorphisms and 2.8 million small indels. Furthermore, PacBio resequencing of 10 accessions was performed, and a total of 2,033 structure variants were identified. Genetic diversity and structure analysis congregated the 198 accessions into three subgroups (Primitive, World, and Japan) and showed the possibility of a long and relatively isolated history of cultivated soybean in Japan. Additionally, the skewed regional distribution of variants in the genome, such as higher structural variations on the R gene clusters in the Japan group, suggested the possibility of selective sweeps during domestication or breeding. A genome-wide association study identified both known and novel causal variants on the genes controlling the flowering period. Novel candidate causal variants were also found on genes related to the seed coat colour by aligning together with Illumina and PacBio reads. The genomic sequences and variants obtained in this study have immense potential to provide information for soybean breeding and genetic studies that may uncover novel alleles or genes involved in agronomically important traits.

Genomic analyses reveal the genetic basis of early maturity and identification of loci and candidate genes in upland cotton (Gossypium hirsutum L.)

Although upland cotton (Gossypium hirsutism L.) originated in the tropics, this early maturity cotton can be planted as far north as 46°N in China due to the accumulation of numerous phenotypic and physiological adaptations during domestication. However, how the genome of early maturity cotton has been altered by strong human selection remains largely unknown. Herein, we report a cotton genome variation map generated by the resequencing of 436 cotton accessions. Whole-genome scans for sweep regions identified 357 putative selection sweeps covering 4.94% (112 Mb) of the upland cotton genome, including 5184 genes. These genes were functionally related to flowering time control, hormone catabolism, ageing and defence response adaptations to environmental changes. A genome-wide association study (GWAS) for seven early maturity traits identified 307 significant loci, 22.48% (69) of which overlapped with putative selection sweeps that occurred during the artificial selection of early maturity cotton. Several previously undescribed candidate genes associated with early maturity were identified by GWAS. This study provides insights into the genetic basis of early maturity in upland cotton as well as breeding resources for cotton improvement.

Deep sampling and pooled amplicon sequencing reveals hidden genic variation in heterogeneous rye accessions

BACKGROUND

Loss of genetic variation negatively impacts breeding efforts and food security. Genebanks house over 7 million accessions representing vast allelic diversity that is a resource for sustainable breeding. Discovery of DNA variations is an important step in the efficient use of these resources. While technologies have improved and costs dropped, it remains impractical to consider resequencing millions of accessions. Candidate genes are known for most agronomic traits, providing a list of high priority targets. Heterogeneity in seed stocks means that multiple samples from an accession need to be evaluated to recover available alleles. To address this we developed a pooled amplicon sequencing approach and applied it to the out-crossing cereal rye (Secale cereale L.).

RESULTS

Using the amplicon sequencing approach 95 rye accessions of different improvement status and worldwide origin, each represented by a pooled sample comprising DNA of 96 individual plants, were evaluated for sequence variation in six candidate genes with significant functions on biotic and abiotic stress resistance, and seed quality. Seventy-four predicted deleterious variants were identified using multiple algorithms. Rare variants were recovered including those found only in a low percentage of seed.

CONCLUSIONS

We conclude that this approach provides a rapid and flexible method for evaluating stock heterogeneity, probing allele diversity, and recovering previously hidden variation. A large extent of within-population heterogeneity revealed in the study provides an important point for consideration during rye germplasm conservation and utilization efforts.

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.

Pan-Genome of Wild and Cultivated Soybeans

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

Polyphenol profile comparisons of seed coats of five pulse crops using a semi-quantitative liquid chromatography-mass spectrometric method

INTRODUCTION

Pulse crops are nutritious and therefore widely grown. Pulse seed coats are typically discarded, despite their high content of polyphenols that are known for their antioxidant properties and health benefits. A better understanding of polyphenol diversity and biochemical pathways will ultimately provide insight into how polyphenols are linked to health benefits, which will help to better utilise these seed coats.

OBJECTIVES

To explore polyphenol profiles among seed coats of diverse genotypes of five pulse crops using a targeted liquid chromatography mass spectrometry (LC-MS) method.

METHODS

Four genotypes of each of common bean, chickpea, pea, lentil and faba bean seed coats were selected for analysis. Following extraction, polyphenols were quantified using LC-MS.

RESULTS

An LC-MS method was developed to quantify 98 polyphenols from 13 different classes in 30 min. The low-tannin seed coats had the lowest concentrations of all polyphenols. Chickpea and pea seed coats had the most similar polyphenolic profiles. The black common bean showed the most diverse seed coat polyphenol profile, including several anthocyanins not detected in any of the other seed coats.

CONCLUSION

The LC-MS method reported herein was used to show polyphenol diversity within several polyphenol classes among the pulse crop seed coats. Detected in all seed coats, flavonols and hydroxybenzoic acids appear well-conserved in the edible Fabaceae. The presence of anthocyanins, flavan-3-ols and proanthocyanins in the coloured seed coats suggests that unique divergent branches were introduced in the flavonoid biosynthetic pathway, possibly in response to environmental stressors.

Toward a "Green Revolution" for Soybean

Soybean (Glycine max), as an economically important food and oilseedcrop, is a major source of plant proteins and oils. Although considerable progress has been made in increasing the yields of rice, wheat, and maize through the "Green Revolution", little improvements have been made for soybean. With the increasing demand of soybean production and the rapid development of crop breeding technologies, time has come for this important crop to undergo a Green Revolution. Here, we briefly summarize the history of crop breeding and Green Revolution in other crops. We then discuss the possible directions and potential approaches toward achieving a Green Revolution for soybean. We provide our views and perspectives on how to breed new soybean varieties with improved yield.

A Domestication-Associated Gene GmPRR3b Regulates the Circadian Clock and Flowering Time in Soybean

Improved soybean cultivars have been adapted to grow at a wide range of latitudes, enabling expansion of cultivation worldwide. However, the genetic basis of this broad adaptation is still not clear. Here, we report the identification of GmPRR3b as a major flowering time regulatory gene that has been selected during domestication and genetic improvement for geographic expansion. Through a genome-wide association study of a diverse soybean landrace panel consisting of 279 accessions, we identified 16 candidate quantitative loci associated with flowering time and maturity time. The strongest signal resides in the known flowering gene E2, verifying the effectiveness of our approach. We detected strong signals associated with both flowering and maturity time in a genomic region containing GmPRR3b. Haplotype analysis revealed that GmPRR3bH6 is the major form of GmPRR3b that has been utilized during recent breeding of modern cultivars. mRNA profiling analysis showed that GmPRR3bH6 displays rhythmic and photoperiod-dependent expression and is preferentially induced under long-day conditions. Overexpression of GmPRR3bH6 increased main stem node number and yield, while knockout of GmPRR3bH6 using CRISPR/Cas9 technology delayed growth and the floral transition. GmPRR3bH6 appears to act as a transcriptional repressor of multiple predicted circadian clock genes, including GmCCA1a, which directly upregulates J/GmELF3a to modulate flowering time. The causal SNP (Chr12:5520945) likely endows GmPRR3bH6 a moderate but appropriate level of activity, leading to early flowering and vigorous growth traits preferentially selected during broad adaptation of landraces and improvement of cultivars.

Exploring the genetic base of the soybean germplasm from Africa, America and Asia as well as mining of beneficial allele for flowering and seed weight

Genetic diversity is the foundation for any breeding program. The present study analyzed the genetic base of 163 soybean genotypes from three continents viz. Africa, America and Asia using 68 trait-linked simple sequence repeats (SSR) markers. The average number of alleles among the germplasm from the three continents followed the trend as Asia (9) > America (8) > Africa (7). Similar trends were observed for gene diversity (0.76 > 0.74 > 0.71) and polymorphism information content (PIC) (0.73 > 0.71 > 0.68). These findings revealed that soybean germplasm from Asia has wider genetic base followed by America, and least in Africa. The 163 genotypes were grouped into 4 clusters by phylogenetic analysis, whereas model-based population structure analysis also divided them into 4 subpopulations comprising 80.61% pure lines and 19.39% admixtures. The genotypes from Africa were easily distinguished from those of other two continents using phylogenetic analysis, indicating important role of geographyical differentiation for this genetic variability. Our results indicated that soybean germplasm has moved from Asia to America, and from America to Africa. Analysis of molecular variance (AMOVA) showed 8.41% variation among the four subpopulations, whereas 63.12% and 28.47% variation existed among and within individuals in the four subpopulations, respectively. Based on the association mapping, a total of 21 SSR markers showed significant association with days to flowering (DoF) and 100-seed weight (HSW). Two markers Satt365 and Satt581 on chromosome 6 and 10, respectively, showed pleiotropic effect or linkage on both traits. Genotype A50 (Gakuran Daizu/PI 506679) from Japan has 8 out of the 13 beneficial alleles for increased HSW. The diverse genotypes, polymorphic SSR markers and desirable alleles identified for DoF and HSW will be used in future breeding programs to improve reproductive, yield and quality traits.

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.

Identification of loci controlling adaptation in Chinese soya bean landraces via a combination of conventional and bioclimatic GWAS

Landraces often contain genetic diversity that has been lost in modern cultivars, including alleles that confer enhanced local adaptation. To comprehensively identify loci associated with adaptive traits in soya bean landraces, for example flowering time, a population of 1938 diverse landraces and 97 accessions of the wild progenitor of cultivated soya bean, Glycine soja was genotyped using tGBS® . Based on 99 085 high-quality SNPs, landraces were classified into three sub-populations which exhibit geographical genetic differentiation. Clustering was inferred from STRUCTURE, principal component analyses and neighbour-joining tree analyses. Using phenotypic data collected at two locations separated by 10 degrees of latitude, 17 trait-associated SNPs (TASs) for flowering time were identified, including a stable locus Chr12:5914898 and previously undetected candidate QTL/genes for flowering time in the vicinity of the previously cloned flowering genes, E1 and E2. Using passport data associated with the collection sites of the landraces, 27 SNPs associated with adaptation to three bioclimatic variables (temperature, daylength, and precipitation) were identified. A series of candidate flowering genes were detected within linkage disequilibrium (LD) blocks surrounding 12 bioclimatic TASs. Nine of these TASs exhibit significant differences in flowering time between alleles within one or more of the three individual sub-populations. Signals of selection during domestication and/or subsequent landrace diversification and adaptation were detected at 38 of the 44 flowering and bioclimatic TASs. Hence, this study lays the groundwork to begin breeding for novel environments predicted to arise following global climate change.

Dissection of soybean populations according to selection signatures based on whole-genome sequences

BACKGROUND

Domestication and improvement processes, accompanied by selections and adaptations, have generated genome-wide divergence and stratification in soybean populations. Simultaneously, soybean populations, which comprise diverse subpopulations, have developed their own adaptive characteristics enhancing fitness, resistance, agronomic traits, and morphological features. The genetic traits underlying these characteristics play a fundamental role in improving other soybean populations.

RESULTS

This study focused on identifying the selection signatures and adaptive characteristics in soybean populations. A core set of 245 accessions (112 wild-type, 79 landrace, and 54 improvement soybeans) selected from 4,234 soybean accessions was re-sequenced. Their genomic architectures were examined according to the domestication and improvement, and accessions were then classified into 3 wild-type, 2 landrace, and 2 improvement subgroups based on various population analyses. Selection and gene set enrichment analyses revealed that the landrace subgroups have selection signals for soybean-cyst nematode HG type 0 and seed development with germination, and that the improvement subgroups have selection signals for plant development with viability and seed development with embryo development, respectively. The adaptive characteristic for soybean-cyst nematode was partially underpinned by multiple resistance accessions, and the characteristics related to seed development were supported by our phenotypic findings for seed weights. Furthermore, their adaptive characteristics were also confirmed as genome-based evidence, and unique genomic regions that exhibit distinct selection and selective sweep patterns were revealed for 13 candidate genes.

CONCLUSIONS

Although our findings require further biological validation, they provide valuable information about soybean breeding strategies and present new options for breeders seeking donor lines to improve soybean populations.

Introgression of novel genetic diversity to improve soybean yield

KEY MESSAGE

Exotic soybean germplasm can be used to increase novel genetic diversity and yield potential of cultivars. Modern North American soybean (Glycine max [L.] Merr.) cultivars have been derived from only a few ancestors. The objectives of this research were to develop breeding lines with novel genetic diversity that were equivalent to the yield of a commercial cultivar parent and within those lines identify regions of novel genetic diversity that were not present in the Corteva Agriscience elite soybean germplasm pool. Nine lines created from diverse germplasm (USDA-ARS breeding program at the University of Illinois) were crossed to a RM34Elite parent to develop populations and sublines for yield testing. Across yield tests at 30 locations conducted between 2014 and 2016, eleven breeding lines were identified that were equivalent to or significantly higher in yield when compared to the RM34Elite parent. Among the eleven final lines, the introgressed novel haplotypes that were not present in current Corteva Agriscience soybean germplasm occupied an estimated 0.8-10.0% of the genome. JH-2665, the highest yielding line across 3 years of testing, yielded 280 kg/ha more than the RM34Elite parent and had an estimated 8.6% of the genome containing novel diversity haplotypes. JH-2665 had 96 regions of novel diversity introgression ranging from 1 to 12 cM in size, with six regions over 6 cM in length. The methods reported demonstrate how high-yielding lines with novel genetic diversity can be developed. This material will be useful for expanding the genetic diversity needed to improve genetic gain in future soybean cultivar development.

Genome-wide variation patterns between landraces and cultivars uncover divergent selection during modern wheat breeding

Genetic diversity, population structure, LD decay, and selective sweeps in 687 wheat accessions were analyzed, providing relevant guidelines to facilitate the use of the germplasm in wheat breeding. Common wheat (Triticum aestivum L.) is one of the most widely grown crops in the world. Landraces were subjected to strong human-mediated selection in developing high-yielding, good quality, and widely adapted cultivars. To investigate the genome-wide patterns of allelic variation, population structure and patterns of selective sweeps during modern wheat breeding, we tested 687 wheat accessions, including landraces (148) and cultivars (539) mainly from China and Pakistan in a wheat 90 K single nucleotide polymorphism array. Population structure analysis revealed that cultivars and landraces from China and Pakistan comprised three relatively independent genetic clusters. Cultivars displayed lower nucleotide diversity and a wider average LD decay across whole genome, indicating allelic erosion and a diversity bottleneck due to the modern breeding. Analysis of genetic differentiation between landraces and cultivars from China and Pakistan identified allelic variants subjected to selection during modern breeding. In total, 477 unique genome regions showed signatures of selection, where 109 were identified in both China and Pakistan germplasm. The majority of genomic regions were located in the B genome (225), followed by the A genome (175), and only 77 regions were located in the D genome. EigenGWAS was further used to identify key selection loci in modern wheat cultivars from China and Pakistan by comparing with global winter wheat and spring wheat diversity panels, respectively. A few known functional genes or loci found within these genome regions corresponded to known phenotypes for disease resistance, vernalization, quality, adaptability and yield-related traits. This study uncovered molecular footprints of modern wheat breeding and explained the genetic basis of polygenic adaptation in wheat. The results will be useful for understanding targets of modern wheat breeding, and in devising future breeding strategies to target beneficial alleles currently not pursued.

Genetic insights into the modification of the pre-fertilization mechanisms during plant domestication

Plant domestication is the process of adapting plants to human use by selecting specific traits. The selection process often involves the modification of some components of the plant reproductive mechanisms. Allelic variants of genes associated with flowering time, vernalization, and the circadian clock are responsible for the adaptation of crops, such as rice, maize, barley, wheat, and tomato, to non-native latitudes. Modifications in the plant architecture and branching have been selected for higher yields and easier harvests. These phenotypes are often produced by alterations in the regulation of the transition of shoot apical meristems to inflorescences, and then to floral meristems. Floral homeotic mutants are responsible for popular double-flower phenotypes in Japanese cherries, roses, camellias, and lilies. The rise of peloric flowers in ornamentals such as snapdragon and florists' gloxinia is associated with non-functional alleles that control the relative expansion of lateral and ventral petals. Mechanisms to force outcrossing such as self-incompatibility have been removed in some tree crops cultivars such as almonds and peaches. In this review, we revisit some of these important concepts from the plant domestication perspective, focusing on four topics related to the pre-fertilization mechanisms: flowering time, inflorescence architecture, flower development, and pre-fertilization self-incompatibility mechanisms.

Complementary gene interaction and xenia effect controls the seed coat colour in interspecific cross between Trifolium alexandrinum and T. apertum

Trifolium alexandrinum (Egyptian clover) is a widely cultivated winter annual fodder. Present work deals with inheritance of the seed coat colour in segregating progenies of the interspecific cross between T. alexandrinum and T. apertum. Although, both the parent species possessed yellow seed coat, the F₁ seeds were black coloured in the reciprocal cross (T. apertum × T. alexandrinum). Seeds borne on individual F₂ plants and the advancing generations segregated in yellow and black seed coat colour, which confirmed xenia effect. F₂ seeds collected from individual F₁ plants exhibited nine black and seven yellow segregation ratio. The segregation of the seed coat colour recorded from F₃ to F₅ generations revealed that yellow seed coat was true breeding (i.e. non-segregating) in this interspecific cross (including the reciprocal crosses). However, the black seeded progenies were either true breeding or segregated in nine black: seven yellow ratio or three black: one yellow ratio suggesting a complementary gene interaction or duplicate recessive epistasis. It indicated that the seed coat colour is controlled by complementary gene interaction along with xenia effect in interspecific crosses between T. alexandrinum and T. apertum. Occurrence of the complementary genes across the species could suggest T. apertum to be the progenitor of T. alexandrinum. Inheritance of seed coat colour in reference to its importance in Egyptian clover breeding is also discussed.

RNA-Seq analysis reveals transcript diversity and active genes after common cutworm (Spodoptera litura Fabricius) attack in resistant and susceptible wild soybean lines

BACKGROUND

Common cutworm (CCW) is highly responsible for destabilizing soybean productivity. Wild soybean is a resource used by breeders to discover elite defensive genes.

RESULTS

The transcriptomes of two wild accessions (W11 and W99) with different resistance to CCW were analyzed at early- and late-induction time points. After induction, the susceptible accession W11 differentially expressed 1268 and 508 genes at the early and late time points, respectively. Compared with W11, the resistant accession W99 differentially expressed 1270 genes at the early time point and many more genes (2308) at the late time point. In total, 3836 non-redundant genes were identified in both lines. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed that the differentially expressed genes (DEGs) in W99 at the late time point were mostly associated with specific processes and pathways. Among the non-redundant genes, 146 genes were commonly up-regulated in the treatment condition compared with the control condition at the early- and late-induction time points in both accessions used in this experiment. Approximately 40% of the common DEGs were related to secondary metabolism, disease resistance, and signal transduction based on their putative function. Excluding the common DEGs, W99 expressed more unique DEGs than W11. Further analysis of the 3836 DEGs revealed that the induction of CCW not only up-regulated defense-related genes, including 37 jasmonic acid (JA)-related genes, 171 plant-pathogen-related genes, and 17 genes encoding protease inhibitors, but also down-regulated growth-related genes, including 35 photosynthesis-related genes, 48 nutrition metabolism genes, and 28 auxin metabolism genes. Therefore, representative defense-related and growth-related genes were chosen for binding site prediction via co-expression of transcription factors (TFs) and spatial expression pattern analyses. In total, 53 binding sites of 28 TFs were identified based on 3 defense-related genes and 3 growth-related genes. Phosphate transporter PT1, which is a representative growth-related gene, was transformed into soybean, and the transgenic soybean plants were susceptible to CCW.

CONCLUSIONS

In summary, we described transcriptome reprograming after herbivore induction in wild soybean, identified the susceptibility of growth-related genes, and provided new resources for the breeding of herbivore-resistant cultivated soybeans.

Salinity stress response and 'omics' approaches for improving salinity stress tolerance in major grain legumes

Sustaining yield gains of grain legume crops under growing salt-stressed conditions demands a thorough understanding of plant salinity response and more efficient breeding techniques that effectively integrate modern omics knowledge. Grain legume crops are important to global food security being an affordable source of dietary protein and essential mineral nutrients to human population, especially in the developing countries. The global productivity of grain legume crops is severely challenged by the salinity stress particularly in the face of changing climates coupled with injudicious use of irrigation water and improper agricultural land management. Plants adapt to sustain under salinity-challenged conditions through evoking complex molecular mechanisms. Elucidating the underlying complex mechanisms remains pivotal to our knowledge about plant salinity response. Improving salinity tolerance of plants demand enriching cultivated gene pool of grain legume crops through capitalizing on 'adaptive traits' that contribute to salinity stress tolerance. Here, we review the current progress in understanding the genetic makeup of salinity tolerance and highlight the role of germplasm resources and omics advances in improving salt tolerance of grain legumes. In parallel, scope of next generation phenotyping platforms that efficiently bridge the phenotyping-genotyping gap and latest research advances including epigenetics is also discussed in context to salt stress tolerance. Breeding salt-tolerant cultivars of grain legumes will require an integrated "omics-assisted" approach enabling accelerated improvement of salt-tolerance traits in crop breeding programs.

Translational genomics and multi-omics integrated approaches as a useful strategy for crop breeding

Recent next generation sequencing-driven mass production of genomic data and multi-omics-integrated approaches have significantly contributed to broadening and deepening our knowledge on the molecular system of living organisms. Accordingly, translational genomics (TG) approach can play a pivotal role in creating an informational bridge between model systems and relatively less studied plants. This review focuses mainly on addressing recent advancement in omics-related technologies, a diverse array of bioinformatic resources and potential applications of TG for the crop breeding. To accomplish above objectives, information on omics data production, various DBs and high throughput technologies was collected, integrated, and used to analyze current status and future perspectives towards omics-assisted crop breeding. Various omics data and resources have been organized and integrated into the databases and/or bioinformatic infrastructures, and thereby serve as the ome's information center for cross-genome translation of biological data. Although the size of accumulated omics data and availability of reference genomes are different among plant families, translational approaches have been actively progressing to access particular biological characteristics. When multi-layered omics data are integrated in a synthetic manner, it will allow providing a stereoscopic view of dynamic molecular behavior and interacting networks of genes occurring in plants. Consequently, TG approach will lead us to broader and deeper insights into target traits for the plant breeding. Furthermore, such systems approach will renovate conventional breeding programs and accelerate precision crop breeding in the future.

Global status of 47 major wheat loci controlling yield, quality, adaptation and stress resistance selected over the last century

BACKGROUND

Wheat breeding over the last 100 years has increased productivity by adapting genotypes to local conditions, but the genomic changes and selection signals that caused phenotypic change during breeding are essentially unknown. Studying and understanding human selection of multiple important genes controlling key phenotypic traits will promote wheat molecular breeding.

RESULTS

A total of 1152 diverse global wheat materials were genotyped based on KASP markers from 47 genes controlling grain yield, grain quality, adaptation, and stress resistance. Significant phenotypic variations between landraces and modern cultivars were found in 11 adaptive and yield-related traits. Thirty-six improvement-selective favorable alleles, including 22 positive prolonged and 14 negative selection alleles, were identified through comparing frequency spectra. Sus1-7A-Hap-H, Sus1-7B-Hap-T, Sus2-2A-Hap-A, TGW6-A1a, Cwi-4A-Hap-C, vrn-A1, PHS1-PHS+ and Lr34+ were subjected to strong selection, and overwhelmingly strong selection had occurred before improvement selection at Psy-A1b, Psy-B1a or b, Psy-D1a and Cwi-5D-Hap-C. However, Rht-B1b, Rht-D1b and 1BL.1RS were rare or absent in Chinese landraces but present in modern Chinese cultivars and introduced accessions. Importantly, Lr68+, Fhb1+, Wx-B1b and Yr15+, currently existing at a low frequency, should be regarded as further major improvement targets in global wheat breeding. Gene flow analysis showed that introduced cultivars especially from the former USSR and Italy contributed to enriched genetic variation in modern Chinese cultivars.

CONCLUSIONS

This work objectively reports human selection on favorable alleles of multiple crucial genes in Asia, Europe, North America and CIMMYT, and traces the distribution of important genes in global wheat for molecular breeding.

Adapting legume crops to climate change using genomic approaches

Our agricultural system and hence food security is threatened by combination of events, such as increasing population, the impacts of climate change, and the need to a more sustainable development. Evolutionary adaptation may help some species to overcome environmental changes through new selection pressures driven by climate change. However, success of evolutionary adaptation is dependent on various factors, one of which is the extent of genetic variation available within species. Genomic approaches provide an exceptional opportunity to identify genetic variation that can be employed in crop improvement programs. In this review, we illustrate some of the routinely used genomics-based methods as well as recent breakthroughs, which facilitate assessment of genetic variation and discovery of adaptive genes in legumes. Although additional information is needed, the current utility of selection tools indicate a robust ability to utilize existing variation among legumes to address the challenges of climate uncertainty.

Identification of candidate domestication-related genes with a systematic survey of loss-of-function mutations

Domestication is an important key co-evolutionary process through which humans have extensively altered the genomic make-up and appearance of both plants and animals. The identification of domestication-related genes remains very arduous. In this study, we present a systematic analytical approach that harnesses two recent advances in genomics, whole-genome sequencing (WGS) and prediction of loss-of-function (LOF) mutations, to greatly facilitate the assembly of an enriched catalogue of domestication-related candidate genes. Using WGS data for 296 cultivated (Glycine max) and 64 wild soybean accessions, we identified 8699 LOF variants, and 116 genes that are uniquely fixed for one or more LOF allele(s) in domesticated soybeans. Existing soybean transcriptomic data led us to overcome analytical challenges associated with whole-genome duplications and to identify neo- or subfunctionalized genes. This systematic approach allowed us to identify 110 candidate domestication-related genes in an efficient and rapid way. This catalogue contains previously well characterized domestication genes in soybean, as well as some orthologs from other domesticated crop species. In addition, it comprises many promising candidate domestication genes. Overall, this collection of candidate domestication-related genes in soybean is almost twice as large as the sum of all previously reported candidate genes in all other crops. We believe this systematic approach could readily be used in wide range of species.

On the Road to Breeding 4.0: Unraveling the Good, the Bad, and the Boring of Crop Quantitative Genomics

Understanding the quantitative genetics of crops has been and will continue to be central to maintaining and improving global food security. We outline four stages that plant breeding either has already achieved or will probably soon achieve. Top-of-the-line breeding programs are currently in Breeding 3.0, where inexpensive, genome-wide data coupled with powerful algorithms allow us to start breeding on predicted instead of measured phenotypes. We focus on three major questions that must be answered to move from current Breeding 3.0 practices to Breeding 4.0: ( a) How do we adapt crops to better fit agricultural environments? ( b) What is the nature of the diversity upon which breeding can act? ( c) How do we deal with deleterious variants? Answering these questions and then translating them to actual gains for farmers will be a significant part of achieving global food security in the twenty-first century.

Transcriptomics profiling in response to cold stress in cultivated rice and weedy rice

Weedy rice is an important germplasm resource for rice improvement because it has useful genes for many abiotic stresses including cold tolerance. We identified the cold tolerance and cold sensitivity of two weedy rice lines (WR 03-35 and WR 03-26) and two cultivated rice lines (Kongyu 131 and 9311). During the seedling stage of these lines, we used RNA-seq to measure changes in weedy rice and cultivated rice whole-genome transcriptome before and after cold treatment. We identified 14,213 and 14,730 differentially expressed genes (DEGs) in cold-tolerant genotypes (WR 03-35, Kongyu 131), and 9219 and 720 DEGs were observed in two cold-sensitive genotypes (WR 03-26, 9311). Many common and special DEGs were analyzed in cold-tolerant and cold-sensitive genotypes, respectively. Some typical genes related to cold stress such as the basic helix-loop-helix (bHLH) gene and leucine-rich repeat (LRR) domain gene etc. The number of these DEGs in cold-tolerant genotypes is more than those found in cold-sensitive genotypes. The gene ontology (GO) enrichment analyses showed significantly enriched terms for biological processes, cellular components and molecular functions. In addition, some genes related to several plant hormones such as abscisic acid (ABA), gibberellic acid (GA), auxin and ethylene were identified. To confirm the RNA-seq data, semi-quantitative RT-PCR and qRT-PCR were performed on 12 randomly selected DEGs. The expression patterns of RNA-seq on these genes corresponded with the semi-quantitative RT-PCR and qRT-PCR method. This study suggests the gene resources related to cold stress from weedy rice could be valuable for understanding the mechanisms involved in cold stress and rice breeding for improving cold tolerance.

Genome-Wide Association Study Identifies Candidate Genes Related to Seed Oil Composition and Protein Content in Gossypium hirsutum L

Cotton (Gossypium spp.) is a leading natural fiber crop and an important source of vegetable protein and oil for humans and livestock. To investigate the genetic architecture of seed nutrients in upland cotton, a genome-wide association study (GWAS) was conducted in a panel of 196 germplasm resources under three environments using a CottonSNP80K chip of 77,774 loci. Relatively high genetic diversity (average gene diversity being 0.331) and phenotypic variation (coefficient of variation, CV, exceeding 3.9%) were detected in this panel. Correlation analysis revealed that the well-documented negative association between seed protein (PR) and oil may be to some extent attributable to the negative correlation between oleic acid (OA) and PR. Linkage disequilibrium (LD) was unevenly distributed among chromosomes and subgenomes. It ranged from 0.10-0.20 Mb (Chr19) to 5.65-5.75 Mb (Chr25) among the chromosomes and the range of Dt-subgenomes LD decay distances was smaller than At-subgenomes. This panel was divided into two subpopulations based on the information of 41,815 polymorphic single-nucleotide polymorphism (SNP) markers. The mixed linear model considering both Q-matrix and K-matrix [MLM(Q+K)] was employed to estimate the association between the SNP markers and the seed nutrients, considering the false positives caused by population structure and the kinship. A total of 47 SNP markers and 28 candidate quantitative trait loci (QTLs) regions were found to be significantly associated with seven cottonseed nutrients, including protein, total fatty acid, and five main fatty acid compositions. In addition, the candidate genes in these regions were analyzed, which included three genes, Gh_D12G1161, Gh_D12G1162, and Gh_D12G1165 that were most likely involved in the control of cottonseed protein concentration. These results improved our understanding of the genetic control of cottonseed nutrients and provided potential molecular tools to develop cultivars with high protein and improved fatty acid compositions in cotton breeding programs through marker-assisted selection.

DNA methylation footprints during soybean domestication and improvement

BACKGROUND

In addition to genetic variation, epigenetic variation plays an important role in determining various biological processes. The importance of natural genetic variation to crop domestication and improvement has been widely investigated. However, the contribution of epigenetic variation in crop domestication at population level has rarely been explored.

RESULTS

To understand the impact of epigenetics on crop domestication, we investigate the variation of DNA methylation during soybean domestication and improvement by whole-genome bisulfite sequencing of 45 soybean accessions, including wild soybeans, landraces, and cultivars. Through methylomic analysis, we identify 5412 differentially methylated regions (DMRs). These DMRs exhibit characters distinct from those of genetically selected regions. In particular, they have significantly higher genetic diversity. Association analyses suggest only 22.54% of DMRs can be explained by local genetic variations. Intriguingly, genes in the DMRs that are not associated with any genetic variation are enriched in carbohydrate metabolism pathways.

CONCLUSIONS

This study provides a valuable map of DNA methylation across diverse accessions and dissects the relationship between DNA methylation variation and genetic variation during soybean domestication, thus expanding our understanding of soybean domestication and improvement.

Genome-Wide Association Studies for Dynamic Plant Height and Number of Nodes on the Main Stem in Summer Sowing Soybeans

Plant height (PH) and the number of nodes on the main stem (NN) serve as major plant architecture traits affecting soybean seed yield. Although many quantitative trait loci for the two traits have been reported, their genetic controls at different developmental stages in soybeans remain unclear. Here, 368 soybean breeding lines were genotyped using 62,423 single nucleotide polymorphism (SNP) markers and phenotyped for the two traits at three different developmental stages over two locations in order to identify their quantitative trait nucleotides (QTNs) using compressed mixed linear model (CMLM) and multi-locus random-SNP-effect mixed linear model (mrMLM) approaches. As a result, 11 and 13 QTNs were found by CMLM to be associated with PH and NN, respectively. Among these QTNs, 8, 3, and 4 for PH and 6, 6, and 8 for NN were found at the three stages, and 3 and 6 were repeatedly detected for PH and NN. In addition, 34 and 30 QTNs were found by mrMLM to be associated with PH and NN, respectively. Among these QTNs, 11, 13, and 16 for PH and 11, 15, and 8 for NN were found at the three stages. A majority of these QTNs overlapped with the previously reported loci. Moreover, one QTN within the known E2 locus for flowering time was detected for the two traits at all three stages, and another that overlapped with the Dt1 locus for stem growth habit was also identified for the two traits at the mature stage. This may explain the highly significant correlation between the two traits. Our findings provide evidence for mixed major plus polygenes inheritance for dynamic traits and an extended understanding of their genetic architecture for molecular dissection and breeding utilization in soybeans.

De novo assembly of a Chinese soybean genome

Soybean was domesticated in China and has become one of the most important oilseed crops. Due to bottlenecks in their introduction and dissemination, soybeans from different geographic areas exhibit extensive genetic diversity. Asia is the largest soybean market; therefore, a high-quality soybean reference genome from this area is critical for soybean research and breeding. Here, we report the de novo assembly and sequence analysis of a Chinese soybean genome for "Zhonghuang 13" by a combination of SMRT, Hi-C and optical mapping data. The assembled genome size is 1.025 Gb with a contig N50 of 3.46 Mb and a scaffold N50 of 51.87 Mb. Comparisons between this genome and the previously reported reference genome (cv. Williams 82) uncovered more than 250,000 structure variations. A total of 52,051 protein coding genes and 36,429 transposable elements were annotated for this genome, and a gene co-expression network including 39,967 genes was also established. This high quality Chinese soybean genome and its sequence analysis will provide valuable information for soybean improvement in the future.

Gene-based SNP identification and validation in soybean using next-generation transcriptome sequencing

Gene-based molecular markers are increasingly used in crop breeding programs for marker-assisted selection. However, identification of genetic variants associated with important agronomic traits has remained a difficult task in soybean. RNA-Seq provides an efficient way, other than assessing global expression variations of coding genes, to discover gene-based SNPs at the whole genome level. In this study, RNA isolated from four soybean accessions each with three replications was subjected to high-throughput sequencing and a range of 44.2-65.9 million paired-end reads were generated for each library. A total of 75,209 SNPs were identified among different genotypes after combination of replications, 89.1% of which were located in expressed regions and 27.0% resulted in amino acid changes. GO enrichment analysis revealed that most significant enriched genes with nonsynonymous SNPs were involved in ribonucleotide binding or catalytic activity. Of 22 SNPs subjected to PCR amplification and Sanger sequencing, all of them were validated. To test the utility of identified SNPs, these validated SNPs were also assessed by genotyping a relative large population with 393 wild and cultivated soybean accessions. These SNPs identified by RNA-Seq provide a useful resource for genetic and genomic studies of soybean. Moreover, the collection of nonsynonymous SNPs annotated with their predicted functional effects also provides a valuable asset for further discovery of genes, identification of gene variants, and development of functional markers.