Resequencing 302 wild and cultivated accessions identifies genes related to domestication and improvement in soybean

Genetic diversity and genomic strategies for improving drought and waterlogging tolerance in soybeans

Drought and its interaction with high temperature are the major abiotic stress factors affecting soybean yield and production stability. Ongoing climate changes are anticipated to intensify drought events, which will further impact crop production and food security. However, excessive water also limits soybean production. The success of soybean breeding programmes for crop improvement is dependent on the extent of genetic variation present in the germplasm base. Screening for natural genetic variation in drought- and flooding tolerance-related traits, including root system architecture, water and nitrogen-fixation efficiency, and yield performance indices, has helped to identify the best resources for genetic studies in soybean. Genomic resources, including whole-genome sequences of diverse germplasms, millions of single-nucleotide polymorphisms, and high-throughput marker genotyping platforms, have expedited gene and marker discovery for translational genomics in soybean. This review highlights the current knowledge of the genetic diversity and quantitative trait loci associated with root system architecture, canopy wilting, nitrogen-fixation ability, and flooding tolerance that contributes to the understanding of drought- and flooding-tolerance mechanisms in soybean. Next-generation mapping approaches and high-throughput phenotyping will facilitate a better understanding of phenotype-genotype associations and help to formulate genomic-assisted breeding strategies, including genomic selection, in soybean for tolerance to drought and flooding stress.

Candidate loci involved in domestication and improvement detected by a published 90K wheat SNP array

Selection is one of the most important forces in crop evolution. Common wheat is a major world food crop and a typical allopolyploid with a huge and complex genome. We applied four approaches to detect loci selected in wheat during domestication and improvement. A total of 7,984 candidate loci were detected, accounting for 23.3% of all 34,317 SNPs analysed, a much higher proportion than estimated in previous reports. We constructed a first generation wheat selection map which revealed the following new insights on genome-wide selection: (1) diversifying selection acted by increasing, decreasing or not affecting gene frequencies; (2) the number of loci under selection during domestication was much higher than that during improvement; (3) the contribution to wheat improvement by the D sub-genome was relatively small due to the bottleneck of hexaploidisation and diversity can be expanded by using synthetic wheat and introgression lines; and (4) clustered selection regions occur throughout the wheat genome, including the centromere regions. This study will not only help future wheat breeding and evolutionary studies, but will also accelerate study of other crops, especially polyploids.

Asymmetric subgenome selection and cis-regulatory divergence during cotton domestication

Comparative population genomics offers an excellent opportunity for unraveling the genetic history of crop domestication. Upland cotton (Gossypium hirsutum) has long been an important economic crop, but a genome-wide and evolutionary understanding of the effects of human selection is lacking. Here, we describe a variation map for 352 wild and domesticated cotton accessions. We scanned 93 domestication sweeps occupying 74 Mb of the A subgenome and 104 Mb of the D subgenome, and identified 19 candidate loci for fiber-quality-related traits through a genome-wide association study. We provide evidence showing asymmetric subgenome domestication for directional selection of long fibers. Global analyses of DNase I-hypersensitive sites and 3D genome architecture, linking functional variants to gene transcription, demonstrate the effects of domestication on cis-regulatory divergence. This study provides new insights into the evolution of gene organization, regulation and adaptation in a major crop, and should serve as a rich resource for genome-based cotton improvement.

Selection for Oil Content During Soybean Domestication Revealed by X-Ray Tomography of Ancient Beans

When and under what circumstances domestication related traits evolved in soybean (Glycine max) is not well understood. Seed size has been a focus of archaeological attention because increased soybean seed weight/size is a trait that distinguishes most modern soybeans from their ancestors; however, archaeological seed size analysis has had limited success. Modern domesticated soybean has a significantly higher oil content than its wild counterpart so oil content is potentially a source of new insight into soybean domestication. We investigated soybean oil content using X-ray computed tomography (CT; specifically, synchrotron radiation X-ray CT or SRX-CT) of charred, archaeological soybean seeds. CT identified holes in the specimens that are associated with oil content. A high oil content facilitates the development of small holes, whereas a high protein content results in larger holes. The volume of small holes increased slowly from 7,500 to 4,000 cal B.P. We infer that human selection for higher oil content began as early as 7,500 cal B.P. and that high oil content cultivars were well established by 4,000 cal B.P.

Potential Uses of Wild Germplasms of Grain Legumes for Crop Improvement

Challenged by population increase, climatic change, and soil deterioration, crop improvement is always a priority in securing food supplies. Although the production of grain legumes is in general lower than that of cereals, the nutritional value of grain legumes make them important components of food security. Nevertheless, limited by severe genetic bottlenecks during domestication and human selection, grain legumes, like other crops, have suffered from a loss of genetic diversity which is essential for providing genetic materials for crop improvement programs. Illustrated by whole-genome-sequencing, wild relatives of crops adapted to various environments were shown to maintain high genetic diversity. In this review, we focused on nine important grain legumes (soybean, peanut, pea, chickpea, common bean, lentil, cowpea, lupin, and pigeonpea) to discuss the potential uses of their wild relatives as genetic resources for crop breeding and improvement, and summarized the various genetic/genomic approaches adopted for these purposes.

Deploying Fourier Coefficients to Unravel Soybean Canopy Diversity

Soybean canopy outline is an important trait used to understand light interception ability, canopy closure rates, row spacing response, which in turn affects crop growth and yield, and directly impacts weed species germination and emergence. In this manuscript, we utilize a methodology that constructs geometric measures of the soybean canopy outline from digital images of canopies, allowing visualization of the genetic diversity as well as a rigorous quantification of shape parameters. Our choice of data analysis approach is partially dictated by the need to efficiently store and analyze large datasets, especially in the context of planned high-throughput phenotyping experiments to capture time evolution of canopy outline which will produce very large datasets. Using the Elliptical Fourier Transformation (EFT) and Fourier Descriptors (EFD), canopy outlines of 446 soybean plant introduction (PI) lines from 25 different countries exhibiting a wide variety of maturity, seed weight, and stem termination were investigated in a field experiment planted as a randomized complete block design with up to four replications. Canopy outlines were extracted from digital images, and subsequently chain coded, and expanded into a shape spectrum by obtaining the Fourier coefficients/descriptors. These coefficients successfully reconstruct the canopy outline, and were used to measure traditional morphometric traits. Highest phenotypic diversity was observed for roundness, while solidity showed the lowest diversity across all countries. Some PI lines had extraordinary shape diversity in solidity. For interpretation and visualization of the complexity in shape, Principal Component Analysis (PCA) was performed on the EFD. PI lines were grouped in terms of origins, maturity index, seed weight, and stem termination index. No significant pattern or similarity was observed among the groups; although interestingly when genetic marker data was used for the PCA, patterns similar to canopy outline traits was observed for origins, and maturity indexes. These results indicate the usefulness of EFT method for reconstruction and study of canopy morphometric traits, and provides opportunities for data reduction of large images for ease in future use.

A genome-wide association study of seed composition traits in wild soybean (Glycine soja)

BACKGROUND

Cultivated soybean (Glycine max) is a major agricultural crop that provides a crucial source of edible protein and oil. Decreased amounts of saturated palmitic acid and increased amounts of unsaturated oleic acid in soybean oil are considered optimal for human cardiovascular health and therefore there has considerable interest by breeders in discovering genes affecting the relative concentrations of these fatty acids. Using a genome-wide association (GWA) approach with nearly 30,000 single nucleotide polymorphisms (SNPs), we investigated the genetic basis of protein, oil and all five fatty acid levels in seeds from a sample of 570 wild soybeans (Glycine soja), the progenitor of domesticated soybean, to identify quantitative trait loci (QTLs) affecting these seed composition traits.

RESULTS

We discovered 29 SNPs located on ten different chromosomes that are significantly associated with the seven seed composition traits in our wild soybean sample. Eight SNPs co-localized with QTLs previously uncovered in linkage or association mapping studies conducted with cultivated soybean samples, while the remaining SNPs appeared to be in novel locations. Twenty-four of the SNPs significantly associated with fatty acid variation, with the majority located on chromosomes 14 (6 SNPs) and seven (8 SNPs). Two SNPs were common for two or more fatty acids, suggesting loci with pleiotropic effects. We also identified some candidate genes that are involved in fatty acid metabolism and regulation. For each of the seven traits, most of the SNPs produced differences between the average phenotypic values of the two homozygotes of about one-half standard deviation and contributed over 3% of their total variability.

CONCLUSIONS

This is the first GWA study conducted on seed composition traits solely in wild soybean populations, and a number of QTLs were found that have not been previously discovered. Some of these may be useful to breeders who select for increased protein/oil content or altered fatty acid ratios in the seeds. The results also provide additional insight into the genetic architecture of these traits in a large sample of wild soybean, and suggest some new candidate genes whose molecular effects on these traits need to be further studied.

Comparison of Genetic Diversity between Chinese and American Soybean (Glycine max (L.)) Accessions Revealed by High-Density SNPs

Soybean is one of the most important economic crops for both China and the United States (US). The exchange of germplasm between these two countries has long been active. In order to investigate genetic relationships between Chinese and US soybean germplasm, 277 Chinese soybean accessions and 300 US soybean accessions from geographically diverse regions were analyzed using 5,361 SNP markers. The genetic diversity and the polymorphism information content (PIC) of the Chinese accessions was higher than that of the US accessions. Population structure analysis, principal component analysis, and cluster analysis all showed that the genetic basis of Chinese soybeans is distinct from that of the USA. The groupings observed in clustering analysis reflected the geographical origins of the accessions; this conclusion was validated with both genetic distance analysis and relative kinship analysis. F_ST-based and EigenGWAS statistical analysis revealed high genetic variation between the two subpopulations. Analysis of the 10 loci with the strongest selection signals showed that many loci were located in chromosome regions that have previously been identified as quantitative trait loci (QTL) associated with environmental-adaptation-related and yield-related traits. The pattern of diversity among the American and Chinese accessions should help breeders to select appropriate parental accessions to enhance the performance of future soybean cultivars.

Genetic Diversity, Population Structure, and Linkage Disequilibrium of an Association-Mapping Panel Revealed by Genome-Wide SNP Markers in Sesame

The characterization of genetic diversity and population structure can be used in tandem to detect reliable phenotype-genotype associations. In the present study, we genotyped a set of 366 sesame germplasm accessions by using 89,924 single-nucleotide polymorphisms (SNPs). The number of SNPs on each chromosome was consistent with the physical length of the respective chromosome, and the average marker density was approximately 2.67 kb/SNP. The genetic diversity analysis showed that the average nucleotide diversity of the panel was 1.1 × 10^-3, with averages of 1.0 × 10^-4, 2.7 × 10^-4, and 3.6 × 10^-4 obtained, respectively for three identified subgroups of the panel: Pop 1, Pop 2, and the Mixed. The genetic structure analysis revealed that these sesame germplasm accessions were structured primarily along the basis of their geographic collection, and that an extensive admixture occurred in the panel. The genome-wide linkage disequilibrium (LD) analysis showed that an average LD extended up to ∼99 kb. The genetic diversity and population structure revealed in this study should provide guidance to the future design of association studies and the systematic utilization of the genetic variation characterizing the sesame panel.

Genome-Wide Association and Transcriptome Analyses Reveal Candidate Genes Underlying Yield-determining Traits in Brassica napus

Yield is one of the most important yet complex crop traits. To improve our understanding of the genetic basis of yield establishment, and to identify candidate genes responsible for yield improvement in Brassica napus, we performed genome-wide association studies (GWAS) for seven yield-determining traits [main inflorescence pod number (MIPN), branch pod number (BPN), pod number per plant (PNP), seed number per pod (SPP), thousand seed weight, main inflorescence yield (MIY), and branch yield], using data from 520 diverse B. napus accessions from two different yield environments. In total, we detected 128 significant single nucleotide polymorphisms (SNPs), 93 of which were revealed as novel by integrative analysis. A combination of GWAS and transcriptome sequencing on 21 haplotype blocks from samples pooled by four extremely high-yielding or low-yielding accessions revealed the differential expression of 14 crucial candiate genes (such as Bna.MYB83, Bna.SPL5, and Bna.ROP3) associated with multiple traits or containing multiple SNPs associated with the same trait. Functional annotation and expression pattern analyses further demonstrated that these 14 candiate genes might be important in developmental processes and biomass accumulation, thus affecting the yield establishment of B. napus. These results provide valuable information for understanding the genetic mechanisms underlying the establishment of high yield in B. napus, and lay the foundation for developing high-yielding B. napus varieties.

A Combined Comparative Transcriptomic, Metabolomic, and Anatomical Analyses of Two Key Domestication Traits: Pod Dehiscence and Seed Dormancy in Pea (Pisum sp.)

The origin of the agriculture was one of the turning points in human history, and a central part of this was the evolution of new plant forms, domesticated crops. Seed dispersal and germination are two key traits which have been selected to facilitate cultivation and harvesting of crops. The objective of this study was to analyze anatomical structure of seed coat and pod, identify metabolic compounds associated with water-impermeable seed coat and differentially expressed genes involved in pea seed dormancy and pod dehiscence. Comparative anatomical, metabolomics, and transcriptomic analyses were carried out on wild dormant, dehiscent Pisum elatius (JI64, VIR320) and cultivated, indehiscent Pisum sativum non-dormant (JI92, Cameor) and recombinant inbred lines (RILs). Considerable differences were found in texture of testa surface, length of macrosclereids, and seed coat thickness. Histochemical and biochemical analyses indicated genotype related variation in composition and heterogeneity of seed coat cell walls within macrosclereids. Liquid chromatography-electrospray ionization/mass spectrometry and Laser desorption/ionization-mass spectrometry of separated seed coats revealed significantly higher contents of proanthocyanidins (dimer and trimer of gallocatechin), quercetin, and myricetin rhamnosides and hydroxylated fatty acids in dormant compared to non-dormant genotypes. Bulk Segregant Analysis coupled to high throughput RNA sequencing resulted in identification of 770 and 148 differentially expressed genes between dormant and non-dormant seeds or dehiscent and indehiscent pods, respectively. The expression of 14 selected dormancy-related genes was studied by qRT-PCR. Of these, expression pattern of four genes: porin (MACE-S082), peroxisomal membrane PEX14-like protein (MACE-S108), 4-coumarate CoA ligase (MACE-S131), and UDP-glucosyl transferase (MACE-S139) was in agreement in all four genotypes with Massive analysis of cDNA Ends (MACE) data. In case of pod dehiscence, the analysis of two candidate genes (SHATTERING and SHATTERPROOF) and three out of 20 MACE identified genes (MACE-P004, MACE-P013, MACE-P015) showed down-expression in dorsal and ventral pod suture of indehiscent genotypes. Moreover, MACE-P015, the homolog of peptidoglycan-binding domain or proline-rich extensin-like protein mapped correctly to predicted Dpo1 locus on PsLGIII. This integrated analysis of the seed coat in wild and cultivated pea provides new insight as well as raises new questions associated with domestication and seed dormancy and pod dehiscence.

Characterization of Insect Resistance Loci in the USDA Soybean Germplasm Collection Using Genome-Wide Association Studies

Management of insects that cause economic damage to yields of soybean mainly rely on insecticide applications. Sources of resistance in soybean plant introductions (PIs) to different insect pests have been reported, and some of these sources, like for the soybean aphid (SBA), have been used to develop resistant soybean cultivars. With the availability of SoySNP50K and the statistical power of genome-wide association studies, we integrated phenotypic data for beet armyworm, Mexican bean beetle (MBB), potato leafhopper (PLH), SBA, soybean looper (SBL), velvetbean caterpillar (VBC), and chewing damage caused by unspecified insects for a comprehensive understanding of insect resistance in the United States Department of Agriculture Soybean Germplasm Collection. We identified significant single nucleotide (SNP) polymorphic markers for MBB, PLH, SBL, and VBC, and we highlighted several leucine-rich repeat-containing genes and myeloblastosis transcription factors within the high linkage disequilibrium region surrounding significant SNP markers. Specifically for soybean resistance to PLH, we found the PLH locus is close but distinct to a locus for soybean pubescence density on chromosome 12. The results provide genetic support that pubescence density may not directly link to PLH resistance. This study offers a novel insight of soybean resistance to four insect pests and reviews resistance mapping studies for major soybean insects.

Genome-Wide Association Study for Pre-harvest Sprouting Resistance in a Large Germplasm Collection of Chinese Wheat Landraces

Pre-harvest sprouting (PHS) is mainly caused by the breaking of seed dormancy in high rainfall regions, which leads to huge economic losses in wheat. In this study, we evaluated 717 Chinese wheat landraces for PHS resistance and carried out genome-wide association studies (GWAS) using to 9,740 DArT-seq and 178,803 SNP markers. Landraces were grown across six environments in China and germination testing of harvest-ripe grain was used to calculate the germination rate (GR) for each accession at each site. GR was highly correlated across all environments. A large number of landraces (194) displayed high levels of PHS resistance (i.e., mean GR < 0.20), which included nine white-grained accessions. Overall, white-grained accessions displayed a significantly higher mean GR (42.7-79.6%) compared to red-grained accessions (19.1-56.0%) across the six environments. Landraces from mesic growing zones in southern China showed higher levels of PHS resistance than those sourced from xeric areas in northern and north-western China. Three main quantitative trait loci (QTL) were detected by GWAS: one on 5D that appeared to be novel and two co-located with the grain color transcription factor Tamyb10 on 3A and 3D. An additional 32 grain color related QTL (GCR-QTL) were detected when the set of red-grained landraces were analyzed separately. GCR-QTL occurred at high frequencies in the red-grained accessions and a strong correlation was observed between the number of GCR-QTL and GR (R² = 0.62). These additional factors could be critical for maintaining high levels of PHS resistance and represent targets for introgression into white-grained wheat cultivars. Further, investigation of the origin of haplotypes associated with the three main QTL revealed that favorable haplotypes for PHS resistance were more common in accessions from higher rainfall zones in China. Thus, a combination of natural and artificial selection likely resulted in landraces incorporating PHS resistance in China.

Development and utilization of a new chemically-induced soybean library with a high mutation density

Mutagenized populations have provided important materials for introducing variation and identifying gene function in plants. In this study, an ethyl methanesulfonate (EMS)-induced soybean (Glycine max) population, consisting of 21,600 independent M2 lines, was developed. Over 1,000 M4 (5) families, with diverse abnormal phenotypes for seed composition, seed shape, plant morphology and maturity that are stably expressed across different environments and generations were identified. Phenotypic analysis of the population led to the identification of a yellow pigmentation mutant, gyl, that displayed significantly decreased chlorophyll (Chl) content and abnormal chloroplast development. Sequence analysis showed that gyl is allelic to MinnGold, where a different single nucleotide polymorphism variation in the Mg-chelatase subunit gene (ChlI1a) results in golden yellow leaves. A cleaved amplified polymorphic sequence marker was developed and may be applied to marker-assisted selection for the golden yellow phenotype in soybean breeding. We show that the newly developed soybean EMS mutant population has potential for functional genomics research and genetic improvement in soybean.

Identification of Substitutions and Small Insertion-Deletions Induced by Carbon-Ion Beam Irradiation in Arabidopsis thaliana

Heavy-ion beam irradiation is one of the principal methods used to create mutants in plants. Research on mutagenic effects and molecular mechanisms of radiation is an important subject that is multi-disciplinary. Here, we re-sequenced 11 mutagenesis progeny (M3) Arabidopsis thaliana lines derived from carbon-ion beam (CIB) irradiation, and subsequently focused on substitutions and small insertion-deletion (INDELs). We found that CIB induced more substitutions (320) than INDELs (124). Meanwhile, the single base INDELs were more prevalent than those in large size (≥2 bp). In details, the detected substitutions showed an obvious bias of C > T transitions, by activating the formation of covalent linkages between neighboring pyrimidine residues in the DNA sequence. An A and T bias was observed among the single base INDELs, in which most of these were induced by replication slippage at either the homopolymer or polynucleotide repeat regions. The mutation rate of 200-Gy CIB irradiation was estimated as 3.37 × 10^-7 per site. Different from previous researches which mainly focused on the phenotype, chromosome aberration, genetic polymorphism, or sequencing analysis of specific genes only, our study revealed genome-wide molecular profile and rate of mutations induced by CIB irradiation. We hope our data could provide valuable clues for explaining the potential mechanism of plant mutation breeding by CIB irradiation.

Back into the wild-Apply untapped genetic diversity of wild relatives for crop improvement

Deleterious effects of climate change and human activities, as well as diverse environmental stresses, present critical challenges to food production and the maintenance of natural diversity. These challenges may be met by the development of novel crop varieties with increased biotic or abiotic resistance that enables them to thrive in marginal lands. However, considering the diverse interactions between crops and environmental factors, it is surprising that evolutionary principles have been underexploited in addressing these food and environmental challenges. Compared with domesticated cultivars, crop wild relatives (CWRs) have been challenged in natural environments for thousands of years and maintain a much higher level of genetic diversity. In this review, we highlight the significance of CWRs for crop improvement by providing examples of CWRs that have been used to increase biotic and abiotic stress resistance/tolerance and overall yield in various crop species. We also discuss the surge of advanced biotechnologies, such as next-generation sequencing technologies and omics, with particular emphasis on how they have facilitated gene discovery in CWRs. We end the review by discussing the available resources and conservation of CWRs, including the urgent need for CWR prioritization and collection to ensure continuous crop improvement for food sustainability.

Recent breeding programs enhanced genetic diversity in both desi and kabuli varieties of chickpea (Cicer arietinum L.)

In order to understand the impact of breeding on genetic diversity and gain insights into temporal trends in diversity in chickpea, a set of 100 chickpea varieties released in 14 countries between 1948 and 2012 were re-sequenced. For analysis, the re-sequencing data for 29 varieties available from an earlier study was also included. Copy number variations and presence absence variations identified in the present study have potential to drive phenotypic variations for trait improvement. Re-sequencing of a large number of varieties has provided opportunities to inspect the genetic and genomic changes reflecting the history of breeding, which we consider as breeding signatures and the selected loci may provide targets for crop improvement. Our study also reports enhanced diversity in both desi and kabuli varieties as a result of recent chickpea breeding efforts. The current study will aid the explicit efforts to breed for local adaptation in the context of anticipated climate changes.

Identification and Analysis of NaHCO3 Stress Responsive Genes in Wild Soybean (Glycine soja) Roots by RNA-seq

Soil alkalinity is a major abiotic constraint to crop productivity and quality. Wild soybean (Glycine soja) is considered to be more stress-tolerant than cultivated soybean (G. max), and has considerable genetic variation for increasing alkalinity tolerance of soybean. In this study, we analyzed the transcriptome profile in the roots of an alkalinity tolerant wild soybean variety N24852 at 12 and 24 h after 90 mM NaHCO₃ stress by RNA-sequencing. Compared with the controls, a total of 449 differentially expressed genes (DEGs) were identified, including 95 and 140 up-regulated genes, and 108 and 135 down-regulated genes at 12 and 24 h after NaHCO₃ treatment, respectively. Quantitative RT-PCR analysis of 14 DEGs showed a high consistency with their expression profiles by RNA-sequencing. Gene Ontology (GO) terms related to transcription factors and transporters were significantly enriched in the up-regulated genes at 12 and 24 h after NaHCO₃ stress, respectively. Nuclear factor Y subunit A transcription factors were enriched at 12 h after NaHCO₃ stress, and high percentages of basic helix-loop-helix, ethylene-responsive factor, Trihelix, and zinc finger (C2H2, C3H) transcription factors were found at both 12 and 24 h after NaHCO₃ stress. Genes related to ion transporters such as ABC transporter, aluminum activated malate transporter, glutamate receptor, nitrate transporter/proton dependent oligopeptide family, and S-type anion channel were enriched in up-regulated DEGs at 24 h after NaHCO₃ treatment, implying their roles in maintaining ion homeostasis in soybean roots under alkalinity. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis showed "phenylpropanoid biosynthesis" and "phenylalanine metabolism" pathways might participate in soybean response to alkalinity. This study provides a foundation to further investigate the functions of NaHCO₃ stress-responsive genes and the molecular basis of soybean tolerance to alkalinity.

Two genomic regions associated with fiber quality traits in Chinese upland cotton under apparent breeding selection

Fiber quality is one of the most important agronomic traits of cotton, and understanding the genetic basis of its target traits will accelerate improvements to cotton fiber quality. In this study, a panel comprising 355 upland cotton accessions was used to perform genome-wide association studies (GWASs) of five fiber quality traits in four environments. A total of 16, 10 and 7 SNPs were associated with fiber length (FL), fiber strength (FS) and fiber uniformity (FU), respectively, based on the mixed linear model (MLM). Most importantly, two major genomic regions (MGR1 and MGR2) on chromosome D_t7 and four potential candidate genes for FL were identified. Analyzing the geographical distribution of favorable haplotypes (FHs) among these lines revealed that two favorable haplotype frequencies (FHFs) were higher in accessions from low-latitude regions than in accessions from high-latitude regions. However, the genetic diversity of lines from the low-latitude regions was lower than the diversity of lines from the high-latitude regions in China. Furthermore, the FHFs differed among cultivars developed during different breeding periods. These results indicate that FHs have undergone artificial selection during upland cotton breeding in recent decades in China and provide a foundation for the further improvement of fiber quality traits.

Advancements in breeding, genetics, and genomics for resistance to three nematode species in soybean

KEY MESSAGE

Integration of genetic analysis, molecular biology, and genomic approaches drastically enhanced our understanding of genetic control of nematode resistance and provided effective breeding strategies in soybeans. Three nematode species, including soybean cyst (SCN, Heterodera glycine), root-knot (RKN, Meloidogyne incognita), and reniform (RN, Rotylenchulus reniformis), are the most destructive pests and have spread to soybean growing areas worldwide. Host plant resistance has played an important role in their control. This review focuses on genetic, genomic studies, and breeding efforts over the past two decades to identify and improve host resistance to these three nematode species. Advancements in genetics, genomics, and bioinformatics have improved our understanding of the molecular and genetic mechanisms of nematode resistance and enabled researchers to generate large-scale genomic resources and marker-trait associations. Whole-genome resequencing, genotyping-by-sequencing, genome-wide association studies, and haplotype analyses have been employed to map and dissect genomic locations for nematode resistance. Recently, two major SCN-resistant loci, Rhg1 and Rhg4, were cloned and other novel resistance quantitative trait loci (QTL) have been discovered. Based on these discoveries, gene-specific DNA markers have been developed for both Rhg1 and Rhg4 loci, which were useful for marker-assisted selection. With RKN resistance QTL being mapped, candidate genes responsible for RKN resistance were identified, leading to the development of functional single nucleotide polymorphism markers. So far, three resistances QTL have been genetically mapped for RN resistance. With nematode species overcoming the host plant resistance, continuous efforts in the identification and deployment of new resistance genes are required to support the development of soybean cultivars with multiple and durable resistance to these pests.

Application of Genomic Technologies to the Breeding of Trees

The recent introduction of next generation sequencing (NGS) technologies represents a major revolution in providing new tools for identifying the genes and/or genomic intervals controlling important traits for selection in breeding programs. In perennial fruit trees with long generation times and large sizes of adult plants, the impact of these techniques is even more important. High-throughput DNA sequencing technologies have provided complete annotated sequences in many important tree species. Most of the high-throughput genotyping platforms described are being used for studies of genetic diversity and population structure. Dissection of complex traits became possible through the availability of genome sequences along with phenotypic variation data, which allow to elucidate the causative genetic differences that give rise to observed phenotypic variation. Association mapping facilitates the association between genetic markers and phenotype in unstructured and complex populations, identifying molecular markers for assisted selection and breeding. Also, genomic data provide in silico identification and characterization of genes and gene families related to important traits, enabling new tools for molecular marker assisted selection in tree breeding. Deep sequencing of transcriptomes is also a powerful tool for the analysis of precise expression levels of each gene in a sample. It consists in quantifying short cDNA reads, obtained by NGS technologies, in order to compare the entire transcriptomes between genotypes and environmental conditions. The miRNAs are non-coding short RNAs involved in the regulation of different physiological processes, which can be identified by high-throughput sequencing of RNA libraries obtained by reverse transcription of purified short RNAs, and by in silico comparison with known miRNAs from other species. All together, NGS techniques and their applications have increased the resources for plant breeding in tree species, closing the former gap of genetic tools between trees and annual species.

Genome-wide association study of 12 agronomic traits in peach

Peach (Prunus persica L.) is a highly valuable crop species and is recognized by molecular researchers as a model fruit for the Rosaceae family. Using whole-genome sequencing data generated from 129 peach accessions, here we perform a comprehensive genome-wide association study for 12 key agronomic traits. We show that among the 10 qualitative traits investigated, nine exhibit consistent and more precise association signals than previously identified by linkage analysis. For two of the qualitative traits, we describe candidate genes, one potentially involved in cell death and another predicted to encode an auxin-efflux carrier, that are highly associated with fruit shape and non-acidity, respectively. Furthermore, we find that several genomic regions harbouring association signals for fruit weight and soluble solid content overlapped with predicted selective sweeps that occurred during peach domestication and improvement. Our findings contribute to the large-scale characterization of genes controlling agronomic traits in peach.

Peach is both an economically important crop species and a model for Rosaceae fruit development research. Here, the authors perform genome-wide association analysis in peach and find candidate genes associated with variation in agronomically important fruit phenotypes.

Functional conservation and divergence of GmCHLI genes in polyploid soybean

Polyploidy is prevalent in nature. As the fate of duplicated genes becomes more complicated when the encoded proteins function as oligomers, functional investigations into duplicated oligomer-encoding genes in polyploid genomes will facilitate our understanding of how traits are expressed. In this study, we identified GmCHLI1, a gene encoding the I subunit of magnesium (Mg)-chelatase, which functions in hexamers as responsible for the semi-dominant etiolation phenotype in soybean. Four GmCHLI copies derived from two polyploidy events were identified in the soybean genome. Further investigation with regard to expression patterns indicated that these four copies have diverged into two pairs; mutation in the other copy of the pair that includes GmCHLI1 also resulted in a chlorophyll-deficient phenotype. Protein interaction assays showed that these four GmCHLIs can interact with each other, but stronger interactions were found with mutated subunits. The results indicate that, in polyploidy, deficiency in each copy of duplicated oligomer-encoding genes could result in a mutant phenotype due to hetero-oligomer formation, which is different from the model of allelic dosage or functional redundancy. In addition, we interestingly found an increase in isoflavonoids in the heterozygous etiolated plants, which might be useful for improving soybean seed quality.

Building strong relationships between conservation genetics and primary industry leads to mutually beneficial genomic advances

Several reviews in the past decade have heralded the benefits of embracing high-throughput sequencing technologies to inform conservation policy and the management of threatened species, but few have offered practical advice on how to expedite the transition from conservation genetics to conservation genomics. Here, we argue that an effective and efficient way to navigate this transition is to capitalize on emerging synergies between conservation genetics and primary industry (e.g., agriculture, fisheries, forestry and horticulture). Here, we demonstrate how building strong relationships between conservation geneticists and primary industry scientists is leading to mutually-beneficial outcomes for both disciplines. Based on our collective experience as collaborative New Zealand-based scientists, we also provide insight for forging these cross-sector relationships.

PGen: large-scale genomic variations analysis workflow and browser in SoyKB

BACKGROUND

With the advances in next-generation sequencing (NGS) technology and significant reductions in sequencing costs, it is now possible to sequence large collections of germplasm in crops for detecting genome-scale genetic variations and to apply the knowledge towards improvements in traits. To efficiently facilitate large-scale NGS resequencing data analysis of genomic variations, we have developed "PGen", an integrated and optimized workflow using the Extreme Science and Engineering Discovery Environment (XSEDE) high-performance computing (HPC) virtual system, iPlant cloud data storage resources and Pegasus workflow management system (Pegasus-WMS). The workflow allows users to identify single nucleotide polymorphisms (SNPs) and insertion-deletions (indels), perform SNP annotations and conduct copy number variation analyses on multiple resequencing datasets in a user-friendly and seamless way.

RESULTS

We have developed both a Linux version in GitHub ( https://github.com/pegasus-isi/PGen-GenomicVariations-Workflow ) and a web-based implementation of the PGen workflow integrated within the Soybean Knowledge Base (SoyKB), ( http://soykb.org/Pegasus/index.php ). Using PGen, we identified 10,218,140 single-nucleotide polymorphisms (SNPs) and 1,398,982 indels from analysis of 106 soybean lines sequenced at 15X coverage. 297,245 non-synonymous SNPs and 3330 copy number variation (CNV) regions were identified from this analysis. SNPs identified using PGen from additional soybean resequencing projects adding to 500+ soybean germplasm lines in total have been integrated. These SNPs are being utilized for trait improvement using genotype to phenotype prediction approaches developed in-house. In order to browse and access NGS data easily, we have also developed an NGS resequencing data browser ( http://soykb.org/NGS_Resequence/NGS_index.php ) within SoyKB to provide easy access to SNP and downstream analysis results for soybean researchers.

CONCLUSION

PGen workflow has been optimized for the most efficient analysis of soybean data using thorough testing and validation. This research serves as an example of best practices for development of genomics data analysis workflows by integrating remote HPC resources and efficient data management with ease of use for biological users. PGen workflow can also be easily customized for analysis of data in other species.

Soybean GmDREBL Increases Lipid Content in Seeds of Transgenic Arabidopsis

A DREB-type transcription factor gene GmDREBL has been characterized for its functions in oil accumulation in seeds. The gene is specifically expressed in soybean seeds. The GmDREBL is localized in nucleus and has transcriptional activation ability. Overexpression of GmDREBL increased the fatty acid content in the seeds of transgenic Arabidopsis plants. GmDREBL can bind to the promoter region of WRI1 to activate its expression. Several other genes in the fatty acid biosynthesis pathway were also enhanced in the GmDREBL-transgenic plants. The GmDREBL can be up-regulated by GmABI3 and GmABI5. Additionally, overexpression of GmDREBL significantly promoted seed size in transgenic plants compared to that of WT plants. Expression of the DREBL is at higher level on the average in cultivated soybeans than that in wild soybeans. The promoter of the DREBL may have been subjected to selection during soybean domestication. Our results demonstrate that GmDREBL participates in the regulation of fatty acid accumulation by controlling the expression of WRI1 and its downstream genes, and manipulation of the gene may increase the oil contents in soybean plants. Our study provides novel insights into the function of DREB-type transcription factors in oil accumulation in addition to their roles in stress response.

Molecular approaches for genetic improvement of seed quality and characterization of genetic diversity in soybean: a critical review

Soybean is an economically important leguminous crop. Genetic improvements of soybeans have focused on enhancement of seed and oil yield, development of varieties suited to different cropping systems, and breeding resistant/tolerant varieties for various biotic and abiotic stresses. Plant breeders have used conventional breeding techniques for the improvement of these traits in soybean. The conventional breeding process can be greatly accelerated through the application of molecular and genomic approaches. Molecular markers have proved to be a new tool in soybean breeding by enhancing selection efficiency in a rapid and time-bound manner. An overview of molecular approaches for the genetic improvement of soybean seed quality parameters, considering recent applications of marker-assisted selection and 'omics' research, is provided in this article.

A Global Analysis of the Polygalacturonase Gene Family in Soybean (Glycine max)

Polygalacturonase is one of the pectin hydrolytic enzymes involved in various developmental and physiological processes such as seed germination, organ abscission, pod and anther dehiscence, and xylem cell formation. To date, no systematic analysis of polygalacturonase incorporating genome organization, gene structure, and expression profiling has been conducted in soybean (Glycine max var. Williams 82). In this study, we identified 112 GmPG genes from the soybean Wm82.a2v1 genome. These genes were classified into three groups, group I (105 genes), group II (5 genes), and group III (2 genes). Fifty-four pairs of duplicate paralogous genes were preferentially identified from duplicated regions of the soybean genome, which implied that long segmental duplications significantly contributed to the expansion of the GmPG gene family. Moreover, GmPG transcripts were analyzed in various tissues using RNA-seq data. The results showed the differential expression of 64 GmPGs in the tissue and partially redundant expression of some duplicate genes, while others showed functional diversity. These findings suggested that the GmPGs were retained by substantial subfunctionalization during the soybean evolutionary processes. Finally, evolutionary analysis based on single nucleotide polymorphisms (SNPs) in wild and cultivated soybeans revealed that 107 GmPGs had selected site(s), which indicated that these genes may have undergone strong selection during soybean domestication. Among them, one non-synonymous SNP of GmPG031 affected floral development during selection, which was consistent with the results of RNA-seq and evolutionary analyses. Thus, our results contribute to the functional characterization of GmPG genes in soybean.

Genomic Signatures of North American Soybean Improvement Inform Diversity Enrichment Strategies and Clarify the Impact of Hybridization

Crop improvement represents a long-running experiment in artificial selection on a complex trait, namely yield. How such selection relates to natural populations is unclear, but the analysis of domesticated populations could offer insights into the relative role of selection, drift, and recombination in all species facing major shifts in selective regimes. Because of the extreme autogamy exhibited by soybean (Glycine max), many "immortalized" genotypes of elite varieties spanning the last century have been preserved and characterized using ∼50,000 single nucleotide polymorphic (SNP) markers. Also due to autogamy, the history of North American soybean breeding can be roughly divided into pre- and posthybridization eras, allowing for direct interrogation of the role of recombination in improvement and selection. Here, we report on genome-wide characterization of the structure and history of North American soybean populations and the signature of selection in these populations. Supporting previous work, we find that maturity defines population structure. Though the diversity of North American ancestors is comparable to available landraces, prehybridization line selections resulted in a clonal structure that dominated early breeding and explains many of the reductions in diversity found in the initial generations of soybean hybridization. The rate of allele frequency change does not deviate sharply from neutral expectation, yet some regions bare hallmarks of strong selection, suggesting a highly variable range of selection strengths biased toward weak effects. We also discuss the importance of haplotypes as units of analysis when complex traits fall under novel selection regimes.

A soybean quantitative trait locus that promotes flowering under long days is identified as FT5a, a FLOWERING LOCUS T ortholog

FLOWERING LOCUS T (FT) is an important floral integrator whose functions are conserved across plant species. In soybean, two orthologs, FT2a and FT5a, play a major role in initiating flowering. Their expression in response to different photoperiods is controlled by allelic combinations at the maturity loci E1 to E4, generating variation in flowering time among cultivars. We determined the molecular basis of a quantitative trait locus (QTL) for flowering time in linkage group J (Chromosome 16). Fine-mapping delimited the QTL to a genomic region of 107kb that harbors FT5a We detected 15 DNA polymorphisms between parents with the early-flowering (ef) and late-flowering (lf) alleles in the promoter region, an intron, and the 3' untranslated region of FT5a, although the FT5a coding regions were identical. Transcript abundance of FT5a was higher in near-isogenic lines for ef than in those for lf, suggesting that different transcriptional activities or mRNA stability caused the flowering time difference. Single-nucleotide polymorphism (SNP) calling from re-sequencing data for 439 cultivated and wild soybean accessions indicated that ef is a rare haplotype that is distinct from common haplotypes including lf The ef allele at FT5a may play an adaptive role at latitudes where early flowering is desirable.

Identification of candidate domestication regions in the radish genome based on high-depth resequencing analysis of 17 genotypes

This study provides high-quality variation data of diverse radish genotypes. Genome-wide SNP comparison along with RNA-seq analysis identified candidate genes related to domestication that have potential as trait-related markers for genetics and breeding of radish. Radish (Raphanus sativus L.) is an annual root vegetable crop that also encompasses diverse wild species. Radish has a long history of domestication, but the origins and selective sweep of cultivated radishes remain controversial. Here, we present comprehensive whole-genome resequencing analysis of radish to explore genomic variation between the radish genotypes and to identify genetic bottlenecks due to domestication in Asian cultivars. High-depth resequencing and multi-sample genotyping analysis of ten cultivated and seven wild accessions obtained 4.0 million high-quality homozygous single-nucleotide polymorphisms (SNPs)/insertions or deletions. Variation analysis revealed that Asian cultivated radish types are closely related to wild Asian accessions, but are distinct from European/American cultivated radishes, supporting the notion that Asian cultivars were domesticated from wild Asian genotypes. SNP comparison between Asian genotypes identified 153 candidate domestication regions (CDRs) containing 512 genes. Network analysis of the genes in CDRs functioning in plant signaling pathways and biochemical processes identified group of genes related to root architecture, cell wall, sugar metabolism, and glucosinolate biosynthesis. Expression profiling of the genes during root development suggested that domestication-related selective advantages included a main taproot with few branched lateral roots, reduced cell wall rigidity and favorable taste. Overall, this study provides evolutionary insights into domestication-related genetic selection in radish as well as identification of gene candidates with the potential to act as trait-related markers for background selection of elite lines in molecular breeding.

Environmental versus geographical effects on genomic variation in wild soybean (Glycine soja) across its native range in northeast Asia

A fundamental goal in evolutionary biology is to understand how various evolutionary factors interact to affect the population structure of diverse species, especially those of ecological and/or agricultural importance such as wild soybean (Glycine soja). G. soja, from which domesticated soybeans (Glycine max) were derived, is widely distributed throughout diverse habitats in East Asia (Russia, Japan, Korea, and China). Here, we utilize over 39,000 single nucleotide polymorphisms genotyped in 99 ecotypes of wild soybean sampled across their native geographic range in northeast Asia, to understand population structure and the relative contribution of environment versus geography to population differentiation in this species. A STRUCTURE analysis identified four genetic groups that largely corresponded to the geographic regions of central China, northern China, Korea, and Japan, with high levels of admixture between genetic groups. A canonical correlation and redundancy analysis showed that environmental factors contributed 23.6% to population differentiation, much more than that for geographic factors (6.6%). Precipitation variables largely explained divergence of the groups along longitudinal axes, whereas temperature variables contributed more to latitudinal divergence. This study provides a foundation for further understanding of the genetic basis of climatic adaptation in this ecologically and agriculturally important species.

Genome-Wide Association Study of Resistance to Soybean Cyst Nematode (Heterodera glycines) HG Type 2.5.7 in Wild Soybean (Glycine soja)

Soybean cyst nematode (SCN) is the most destructive soybean pest worldwide. Host plant resistance is the most environmentally friendly and cost-effective way of mitigating SCN damage to soybeans. However, overuse of the resistant soybean [Glycine max (L.) Merr.] cultivars from limited genetic resources has resulted in SCN race shifts in many soybean-growing areas. Thus, exploration of novel sources of SCN resistance and dissection of the genetic basis are urgently needed. In this study, we screened 235 wild soybean (Glycine soja Sieb. & Zucc.) accessions to identify genotypes resistant to SCN HG Type 2.5.7 (race 5), a less investigated type but is prevalent in the southeastern US. We also dissected the genetic basis of SCN resistance using a genome-wide association study with SNPs genotyped by SoySNP50k iSelect BeadChip. In total, 43 resistant accessions (female index < 30) were identified, with 10 SNPs being significantly associated with SCN HG 2.5.7 resistance in this wild species. Furthermore, four significant SNPs were localized to linked regions of the known quantitative trait locus (QTL) rhg1 on chromosome 18. The other four SNPs on chromosome 18 and two SNPs on chromosome 19 are novel. Genes encoding disease resistance-related proteins with a leucine-rich region, a mitogen-activated protein kinase (MAPK) on chromosome 18, and a MYB transcription factor on chromosome 19 were identified as promising candidate genes. The identified SNPs and candidate genes will benefit future marker-assisted breeding and dissection of the molecular mechanisms underlying the soybean-SCN interaction.

Environmental versus geographical effects on genomic variation in wild soybean (Glycine soja) across its native range in northeast Asia

A fundamental goal in evolutionary biology is to understand how various evolutionary factors interact to affect the population structure of diverse species, especially those of ecological and/or agricultural importance such as wild soybean (Glycine soja). G. soja, from which domesticated soybeans (Glycine max) were derived, is widely distributed throughout diverse habitats in East Asia (Russia, Japan, Korea, and China). Here, we utilize over 39,000 single nucleotide polymorphisms genotyped in 99 ecotypes of wild soybean sampled across their native geographic range in northeast Asia, to understand population structure and the relative contribution of environment versus geography to population differentiation in this species. A STRUCTURE analysis identified four genetic groups that largely corresponded to the geographic regions of central China, northern China, Korea, and Japan, with high levels of admixture between genetic groups. A canonical correlation and redundancy analysis showed that environmental factors contributed 23.6% to population differentiation, much more than that for geographic factors (6.6%). Precipitation variables largely explained divergence of the groups along longitudinal axes, whereas temperature variables contributed more to latitudinal divergence. This study provides a foundation for further understanding of the genetic basis of climatic adaptation in this ecologically and agriculturally important species.

Improvement in nitrogen fixation capacity could be part of the domestication process in soybean

Biological nitrogen fixation (BNF) in soybeans is a complex process involving the interplay between the plant host and the symbiotic rhizobia. As nitrogen supply has a crucial role in growth and development, higher nitrogen fixation capacity would be important to achieve bigger plants and larger seeds, which were important selection criteria during plant domestication by humans. To test this hypothesis, we monitored the nitrogen fixation-related performance in 31 cultivated and 17 wild soybeans after inoculation with the slow-growing Bradyrhizobium diazoefficiens sp. nov. USDA110 and the fast-growing Sinorhizobium (Ensifer) fredii CCBAU45436. Our results showed that, in general, cultivated soybeans gave better performance in BNF. Electron microscopic studies indicated that there was an exceptionally high accumulation of poly-β-hydroxybutyrate bodies in bacteroids in the nodules of all wild soybeans tested, suggesting that the C/N balance in wild soybeans may not be optimized for nitrogen fixation. Furthermore, we identified new quantitative trait loci (QTLs) for total ureides and total nodule fresh weight by employing a recombinant inbred population composed of descendants from a cross between a cultivated and a wild parent. Using nucleotide diversity (θπ), divergence index (Fst) and distribution of fixed single-nucleotide polymorphisms as parameters, we found that some regions in the total ureides QTL on chromosome 17 and the total nodule fresh weight QTL on chromosome 12 exhibited very low diversity among cultivated soybeans, suggesting that these were traits specially selected during the domestication and breeding process.

Development and Validation of a High-Density SNP Genotyping Array for African Oil Palm

High-density single nucleotide polymorphism (SNP) genotyping arrays are powerful tools that can measure the level of genetic polymorphism within a population. To develop a whole-genome SNP array for oil palms, SNP discovery was performed using deep resequencing of eight libraries derived from 132 Elaeis guineensis and Elaeis oleifera palms belonging to 59 origins, resulting in the discovery of >3 million putative SNPs. After SNP filtering, the Illumina OP200K custom array was built with 170 860 successful probes. Phenetic clustering analysis revealed that the array could distinguish between palms of different origins in a way consistent with pedigree records. Genome-wide linkage disequilibrium declined more slowly for the commercial populations (ranging from 120 kb at r(2) = 0.43 to 146 kb at r(2) = 0.50) when compared with the semi-wild populations (19.5 kb at r(2) = 0.22). Genetic fixation mapping comparing the semi-wild and commercial population identified 321 selective sweeps. A genome-wide association study (GWAS) detected a significant peak on chromosome 2 associated with the polygenic component of the shell thickness trait (based on the trait shell-to-fruit; S/F %) in tenera palms. Testing of a genomic selection model on the same trait resulted in good prediction accuracy (r = 0.65) with 42% of the S/F % variation explained. The first high-density SNP genotyping array for oil palm has been developed and shown to be robust for use in genetic studies and with potential for developing early trait prediction to shorten the oil palm breeding cycle.

First-generation HapMap in Cajanus spp. reveals untapped variations in parental lines of mapping populations

Whole genome re-sequencing (WGRS) was conducted on a panel of 20 Cajanus spp. accessions (crossing parentals of recombinant inbred lines, introgression lines, multiparent advanced generation intercross and nested association mapping population) comprising of two wild species and 18 cultivated species accessions. A total of 791.77 million paired-end reads were generated with an effective mapping depth of ~12X per accession. Analysis of WGRS data provided 5 465 676 genome-wide variations including 4 686 422 SNPs and 779 254 InDels across the accessions. Large structural variations in the form of copy number variations (2598) and presence and absence variations (970) were also identified. Additionally, 2 630 904 accession-specific variations comprising of 2 278 571 SNPs (86.6%), 166 243 deletions (6.3%) and 186 090 insertions (7.1%) were also reported. Identified polymorphic sites in this study provide the first-generation HapMap in Cajanus spp. which will be useful in mapping the genomic regions responsible for important traits.

1,135 Genomes Reveal the Global Pattern of Polymorphism in Arabidopsis thaliana

Arabidopsis thaliana serves as a model organism for the study of fundamental physiological, cellular, and molecular processes. It has also greatly advanced our understanding of intraspecific genome variation. We present a detailed map of variation in 1,135 high-quality re-sequenced natural inbred lines representing the native Eurasian and North African range and recently colonized North America. We identify relict populations that continue to inhabit ancestral habitats, primarily in the Iberian Peninsula. They have mixed with a lineage that has spread to northern latitudes from an unknown glacial refugium and is now found in a much broader spectrum of habitats. Insights into the history of the species and the fine-scale distribution of genetic diversity provide the basis for full exploitation of A. thaliana natural variation through integration of genomes and epigenomes with molecular and non-molecular phenotypes.

The Power of Inbreeding: NGS-Based GWAS of Rice Reveals Convergent Evolution during Rice Domestication

Low-coverage whole-genome sequencing is an effective strategy for genome-wide association studies in humans, due to the availability of large reference panels for genotype imputation. However, it is unclear whether this strategy can be utilized in other species without reference panels. Using simulations, we show that this approach is even more relevant in inbred species such as rice (Oryza sativa L.), which are effectively haploid, allowing easy haplotype construction and imputation-based genotype calling, even without the availability of large reference panels. We sequenced 203 rice varieties with well-characterized phenotypes from the United States Department of Agriculture Rice Mini-Core Collection at an average depth of 1.5× and used the data for mapping three traits. For the first two traits, amylose content and seed length, our approach leads to direct identification of the previously identified causal SNPs in the major-effect loci. For the third trait, pericarp color, an important trait underwent selection during domestication, we identified a new major-effect locus. Although known loci can explain color variation in the varieties of two main subspecies of Asian domesticated rice, japonica and indica, the new locus identified is unique to another domesticated rice subgroup, aus, and together with existing loci, can fully explain the major variation in pericarp color in aus. Our discovery of a unique genetic basis of white pericarp in aus provides an example of convergent evolution during rice domestication and suggests that aus may have a domestication history independent of japonica and indica.

Identification and Validation of Candidate Genes Associated with Domesticated and Improved Traits in Soybean

Soybean, an important source of vegetable oils and proteins for humans, has undergone significant phenotypic changes during domestication and improvement. However, there is limited knowledge about genes related to these domesticated and improved traits, such as flowering time, seed development, alkaline-salt tolerance, and seed oil content (SOC). In this study, more than 106,000 single nucleotide polymorphisms (SNPs) were identified by restriction site associated DNA sequencing of 14 wild, 153 landrace, and 119 bred soybean accessions, and 198 candidate domestication regions (CDRs) were identified via multiple genetic diversity analyses. Of the 1489 candidate domestication genes (CDGs) within these CDRs, a total of 330 CDGs were related to the above four traits in the domestication, gene ontology (GO) enrichment, gene expression, and pathway analyses. Eighteen, 60, 66, and 10 of the 330 CDGs were significantly associated with the above four traits, respectively. Of 134 trait-associated CDGs, 29 overlapped with previous CDGs, 11 were consistent with candidate genes in previous trait association studies, and 66 were covered by the domesticated and improved quantitative trait loci or their adjacent regions, having six common CDGs, such as one functionally characterized gene (). Of the 68 seed size (SS) and SOC CDGs, 37 were further confirmed by gene expression analysis. In addition, eight genes were found to be related to artificial selection during modern breeding. Therefore, this study provides an integrated method for efficiently identifying CDGs and valuable information for domestication and genetic research.

Genome engineering and plant breeding: impact on trait discovery and development

KEY MESSAGE

New tools for the precise modification of crops genes are now available for the engineering of new ideotypes. A future challenge in this emerging field of genome engineering is to develop efficient methods for allele mining. Genome engineering tools are now available in plants, including major crops, to modify in a predictable manner a given gene. These new techniques have a tremendous potential for a spectacular acceleration of the plant breeding process. Here, we discuss how genetic diversity has always been the raw material for breeders and how they have always taken advantage of the best available science to use, and when possible, increase, this genetic diversity. We will present why the advent of these new techniques gives to the breeders extremely powerful tools for crop breeding, but also why this will require the breeders and researchers to characterize the genes underlying this genetic diversity more precisely. Tackling these challenges should permit the engineering of optimized alleles assortments in an unprecedented and controlled way.

Genome-wide association study using whole-genome sequencing rapidly identifies new genes influencing agronomic traits in rice

A genome-wide association study (GWAS) can be a powerful tool for the identification of genes associated with agronomic traits in crop species, but it is often hindered by population structure and the large extent of linkage disequilibrium. In this study, we identified agronomically important genes in rice using GWAS based on whole-genome sequencing, followed by the screening of candidate genes based on the estimated effect of nucleotide polymorphisms. Using this approach, we identified four new genes associated with agronomic traits. Some genes were undetectable by standard SNP analysis, but we detected them using gene-based association analysis. This study provides fundamental insights relevant to the rapid identification of genes associated with agronomic traits using GWAS and will accelerate future efforts aimed at crop improvement.

Soybean (Glycine max) Mutant and Germplasm Resources: Current Status and Future Prospects

Genetic bottlenecks during domestication and modern breeding limited the genetic diversity of soybean (Glycine max (L.) Merr.). Therefore, expanding and diversifying soybean genetic resources is a major priority for the research community. These resources, consisting of natural and induced genetic variants, are valuable tools for improving soybean and furthering soybean biological knowledge. During the twentieth century, researchers gathered a wealth of genetic variation in the forms of landraces, Glycine soja accessions, Glycine tertiary germplasm, and the U.S. Department of Agriculture (USDA) Type and Isoline Collections. During the twenty-first century, soybean researchers have added several new genetic and genomic resources. These include the reference genome sequence, genotype data for the USDA soybean germplasm collection, next-generation mapping populations, new irradiation and transposon-based mutagenesis populations, and designer nuclease platforms for genome engineering. This paper briefly surveys the publicly accessible soybean genetic resources currently available or in development and provides recommendations for developing such genetic resources in the future. © 2016 by John Wiley & Sons, Inc.

The transcriptomic signature of developing soybean seeds reveals the genetic basis of seed trait adaptation during domestication

Cultivated soybean has undergone many transformations during domestication. In this paper we report a comprehensive assessment of the evolution of gene co-expression networks based on the analysis of 40 transcriptomes from developing soybean seeds in cultivated and wild soybean accessions. We identified 2680 genes that are differentially expressed during seed maturation and established two cultivar-specific gene co-expression networks. Through analysis of the two networks and integration with quantitative trait locus data we identified two potential key drivers for seed trait formation, GA20OX and NFYA. GA20OX encodes an enzyme in a rate-limiting step of gibberellin biosynthesis, and NFYA encodes a transcription factor. Overexpression of GA20OX and NFYA enhanced seed size/weight and oil content, respectively, in seeds of transgenic plants. The two genes showed significantly higher expression in cultivated than in wild soybean, and the increases in expression were associated with genetic variations in the promoter region of each gene. Moreover, the expression of GA20OX and NFYA in seeds of soybean accessions correlated with seed weight and oil content, respectively. Our study reveals transcriptional adaptation during soybean domestication and may identify a mechanism of selection by expression for seed trait formation, providing strategies for future breeding practice.

Genomic variation and DNA repair associated with soybean transgenesis: a comparison to cultivars and mutagenized plants

BACKGROUND

The safety of mutagenized and genetically transformed plants remains a subject of scrutiny. Data gathered and communicated on the phenotypic and molecular variation induced by gene transfer technologies will provide a scientific-based means to rationally address such concerns. In this study, genomic structural variation (e.g. large deletions and duplications) and single nucleotide polymorphism rates were assessed among a sample of soybean cultivars, fast neutron-derived mutants, and five genetically transformed plants developed through Agrobacterium based transformation methods.

RESULTS

On average, the number of genes affected by structural variations in transgenic plants was one order of magnitude less than that of fast neutron mutants and two orders of magnitude less than the rates observed between cultivars. Structural variants in transgenic plants, while rare, occurred adjacent to the transgenes, and at unlinked loci on different chromosomes. DNA repair junctions at both transgenic and unlinked sites were consistent with sequence microhomology across breakpoints. The single nucleotide substitution rates were modest in both fast neutron and transformed plants, exhibiting fewer than 100 substitutions genome-wide, while inter-cultivar comparisons identified over one-million single nucleotide polymorphisms.

CONCLUSIONS

Overall, these patterns provide a fresh perspective on the genomic variation associated with high-energy induced mutagenesis and genetically transformed plants. The genetic transformation process infrequently results in novel genetic variation and these rare events are analogous to genetic variants occurring spontaneously, already present in the existing germplasm, or induced through other types of mutagenesis. It remains unclear how broadly these results can be applied to other crops or transformation methods.

DNA Sequence Evolution and Rare Homoeologous Conversion in Tetraploid Cotton

Allotetraploid cotton species are a vital source of spinnable fiber for textiles. The polyploid nature of the cotton genome raises many evolutionary questions as to the relationships between duplicated genomes. We describe the evolution of the cotton genome (SNPs and structural variants) with the greatly improved resolution of 34 deeply re-sequenced genomes. We also explore the evolution of homoeologous regions in the AT- and DT-genomes and especially the phenomenon of conversion between genomes. We did not find any compelling evidence for homoeologous conversion between genomes. These findings are very different from other recent reports of frequent conversion events between genomes. We also identified several distinct regions of the genome that have been introgressed between G. hirsutum and G. barbadense, which presumably resulted from breeding efforts targeting associated beneficial alleles. Finally, the genotypic data resulting from this study provides access to a wealth of diversity sorely needed in the narrow germplasm of cotton cultivars.

Molecular and geographic evolutionary support for the essential role of GIGANTEAa in soybean domestication of flowering time

BACKGROUND

Flowering time is a domestication trait of Glycine max and varies in soybeans, yet, a gene for flowering time variation has not been associated with soybean domestication. GIGANTEA (GI) is a major gene involved in the control of flowering time in Arabidopsis, although three GI homologs complicate this model in the soybean genome.

RESULTS

In the present work, we revealed that the geographic evolution of the GIGANTEAa (GIa) haplotypes in G. max (GmGIa) and Glycine soja (GsGIa). Three GIa haplotypes (H1, H2, and H3) were found among cultivated soybeans and their wild relatives, yet an additional 44 diverse haplotypes were observed in wild soybeans. H1 had a premature stop codon in the 10(th) exon, whereas the other haplotypes encoded full-length GIa protein isoforms. In both wild-type and cultivated soybeans, H2 was present in the Southern region of China, and H3 was restricted to areas near the Northeast region of China. H1 was genetically derived from H2, and it was dominant and widely distributed among cultivated soybeans, whereas in wild populations, the ortholog of this domesticated haplotype H1 was only found in Yellow River basin with a low frequency. Moreover, this mutated GIa haplotype significantly correlated with early flowering. We further determined that the differences in gene expression of the three GmGIa haplotypes were not correlated to flowering time variations in cultivated soybeans. However, only the truncated GmGIa H1 could partially rescue gi-2 Arabidopsis from delayed flowering in transgenic plants, whereas both GmGIa H2 and H3 haplotypes could significantly repress flowering in transgenic Arabidopsis with a wild-type background.

CONCLUSIONS

Thus, GmGIa haplotype diversification may have contributed to flowering time adaptation that facilitated the radiation of domesticated soybeans. In light of the evolution of the GIa gene, soybean domestication history for an early flowering phenotype is discussed.

Environmental Association Analyses Identify Candidates for Abiotic Stress Tolerance in Glycine soja, the Wild Progenitor of Cultivated Soybeans

Natural populations across a species range demonstrate population structure owing to neutral processes such as localized origins of mutations and migration limitations. Selection also acts on a subset of loci, contributing to local adaptation. An understanding of the genetic basis of adaptation to local environmental conditions is a fundamental goal in basic biological research. When applied to crop wild relatives, this same research provides the opportunity to identify adaptive genetic variation that may be used to breed for crops better adapted to novel or changing environments. The present study explores an ex situ conservation collection, the USDA germplasm collection, genotyped at 32,416 SNPs to identify population structure and test for associations with bioclimatic and biophysical variables in Glycine soja, the wild progenitor of Glycine max (soybean). Candidate loci were detected that putatively contribute to adaptation to abiotic stresses. The identification of potentially adaptive variants in this ex situ collection may permit a more targeted use of germplasm collections.