A reference-grade wild soybean genome.

Efficient crop improvement depends on the application of accurate genetic information contained in diverse germplasm resources. Here we report a reference-grade genome of wild soybean accession W05, with a final assembled genome size of 1013.2 Mb and a contig N50 of 3.3 Mb. The analytical power of the W05 genome is demonstrated by several examples. First, we identify an inversion at the locus determining seed coat color during domestication. Second, a translocation event between chromosomes 11 and 13 of some genotypes is shown to interfere with the assignment of QTLs. Third, we find a region containing copy number variations of the Kunitz trypsin inhibitor (KTI) genes. Such findings illustrate the power of this assembly in the analysis of large structural variations in soybean germplasm collections. The wild soybean genome assembly has wide applications in comparative genomic and evolutionary studies, as well as in crop breeding and improvement programs.

Whole-genome re-sequencing reveals the impact of the interaction of copy number variants of the rhg1 and Rhg4 genes on broad-based resistance to soybean cyst nematode.

Soybean cyst nematode (SCN) is the most devastating plant-parasitic nematode. Most commercial soybean varieties with SCN resistance are derived from PI88788. Resistance derived from PI88788 is breaking down due to narrow genetic background and SCN population shift. PI88788 requires mainly the rhg1-b locus, while 'Peking' requires rhg1-a and Rhg4 for SCN resistance. In the present study, whole genome re-sequencing of 106 soybean lines was used to define the Rhg haplotypes and investigate their responses to the SCN HG-Types. The analysis showed a comprehensive profile of SNPs and copy number variations (CNV) at these loci. CNV of rhg1 (GmSNAP18) only contributed towards resistance in lines derived from PI88788 and 'Cloud'. At least 5.6 copies of the PI88788-type rhg1 were required to confer SCN resistance, regardless of the Rhg4 (GmSHMT08) haplotype. However, when the GmSNAP18 copies dropped below 5.6, a 'Peking'-type GmSHMT08 haplotype was required to ensure SCN resistance. This points to a novel mechanism of epistasis between GmSNAP18 and GmSHMT08 involving minimum requirements for copy number. The presence of more Rhg4 copies confers resistance to multiple SCN races. Moreover, transcript abundance of the GmSHMT08 in root tissue correlates with more copies of the Rhg4 locus, reinforcing SCN resistance. Finally, haplotype analysis of the GmSHMT08 and GmSNAP18 promoters inferred additional levels of the resistance mechanism. This is the first report revealing the genetic basis of broad-based resistance to SCN and providing new insight into epistasis, haplotype-compatibility, CNV, promoter variation and its impact on broad-based disease resistance in plants.

Dissecting genomic hotspots underlying seed protein, oil, and sucrose content in an interspecific mapping population of soybean using high-density linkage mapping.

The cultivated [Glycine max (L) Merr.] and wild [Glycine soja Siebold & Zucc.] soybean species comprise wide variation in seed composition traits. Compared to wild soybean, cultivated soybean contains low protein, high oil, and high sucrose. In this study, an interspecific population was derived from a cross between G. max (Williams 82) and G. soja (PI 483460B). This recombinant inbred line (RIL) population of 188 lines was sequenced at 0.3× depth. Based on 91 342 single nucleotide polymorphisms (SNPs), recombination events in RILs were defined, and a high-resolution bin map was developed (4070 bins). In addition to bin mapping, quantitative trait loci (QTL) analysis for protein, oil, and sucrose was performed using 3343 polymorphic SNPs (3K-SNP), derived from Illumina Infinium BeadChip sequencing platform. The QTL regions from both platforms were compared, and a significant concordance was observed between bin and 3K-SNP markers. Importantly, the bin map derived from next-generation sequencing technology enhanced mapping resolution (from 1325 to 50 Kb). A total of five, nine, and four QTLs were identified for protein, oil, and sucrose content, respectively, and some of the QTLs coincided with soybean domestication-related genomic loci. The major QTL for protein and oil were mapped on Chr. 20 (qPro_20) and suggested negative correlation between oil and protein. In terms of sucrose content, a novel and major QTL were identified on Chr. 8 (qSuc_08) and harbours putative genes involved in sugar transport. In addition, genome-wide association using 91 342 SNPs confirmed the genomic loci derived from QTL mapping. A QTL-based haplotype using whole-genome resequencing of 106 diverse soybean lines identified unique allelic variation in wild soybean that could be utilized to widen the genetic base in cultivated soybean.

Genomic-assisted phylogenetic analysis and marker development for next generation soybean cyst nematode resistance breeding.

Soybean cyst nematode (SCN, Heterodera glycines Ichinohe) is a serious soybean pest. The use of resistant cultivars is an effective approach for preventing yield loss. In this study, 19,652 publicly available soybean accessions that were previously genotyped with the SoySNP50K iSelect BeadChip were used to evaluate the phylogenetic diversity of SCN resistance genes Rhg1 and Rhg4 in an attempt to identify novel sources of resistance. The sequence information of soybean lines was utilized to develop KASPar (KBioscience Competitive Allele-Specific PCR) assays from single nucleotide polymorphisms (SNPs) of Rhg1, Rhg4, and other novel quantitative trait loci (QTL). These markers were used to genotype a diverse set of 95 soybean germplasm lines and three recombinant inbred line (RIL) populations. SNP markers from the Rhg1 gene were able to differentiate copy number variation (CNV), such as resistant-high copy (PI 88788-type), low copy (Peking-type), and susceptible-single copy (Williams 82) numbers. Similarly, markers for the Rhg4 gene were able to detect Peking-type (resistance) genotypes. The phylogenetic information of SCN resistance loci from a large set of soybean accessions and the gene/QTL specific markers that were developed in this study will accelerate SCN resistance breeding programs.

A Sweetpotato Geranylgeranyl Pyrophosphate Synthase Gene, IbGGPS, Increases Carotenoid Content and Enhances Osmotic Stress Tolerance in Arabidopsis thaliana.

Sweetpotato highly produces carotenoids in storage roots. In this study, a cDNA encoding geranylgeranyl phyrophosphate synthase (GGPS), named IbGGPS, was isolated from sweetpotato storage roots. Green fluorescent protein (GFP) was fused to the C-terminus of IbGGPS to obtain an IbGGPS-GFP fusion protein that was transiently expressed in both epidermal cells of onion and leaves of tobacco. Confocal microscopic analysis determined that the IbGGPS-GFP protein was localized to specific areas of the plasma membrane of onion and chloroplasts in tobacco leaves. The coding region of IbGGPS was cloned into a binary vector under the control of 35S promoter and then transformed into Arabidopsis thaliana to obtain transgenic plants. High performance liquid chromatography (HPLC) analysis showed a significant increase of total carotenoids in transgenic plants. The seeds of transgenic and wild-type plants were germinated on an agar medium supplemented with polyethylene glycol (PEG). Transgenic seedlings grew significantly longer roots than wild-type ones did. Further enzymatic analysis showed an increased activity of superoxide dismutase (SOD) in transgenic seedlings. In addition, the level of malondialdehyde (MDA) was reduced in transgenics. qRT-PCR analysis showed altered expressions of several genes involved in the carotenoid biosynthesis in transgenic plants. These data results indicate that IbGGPS is involved in the biosynthesis of carotenoids in sweetpotato storage roots and likely associated with tolerance to osmotic stress.

Expression of root-related transcription factors associated with flooding tolerance of soybean (Glycine max).

Much research has been conducted on the changes in gene expression of the model plant Arabidopsis to low-oxygen stress. Flooding results in a low oxygen environment in the root zone. However, there is ample evidence that tolerance to soil flooding is more than tolerance to low oxygen alone. In this study, we investigated the physiological response and differential expression of root-related transcription factors (TFs) associated with the tolerance of soybean plants to soil flooding. Differential responses of PI408105A and S99-2281 plants to ten days of soil flooding were evaluated at physiological, morphological and anatomical levels. Gene expression underlying the tolerance response was investigated using qRT-PCR of root-related TFs, known anaerobic genes, and housekeeping genes. Biomass of flood-sensitive S99-2281 roots remained unchanged during the entire 10 days of flooding. Flood-tolerant PI408105A plants exhibited recovery of root growth after 3 days of flooding. Flooding induced the development of aerenchyma and adventitious roots more rapidly in the flood-tolerant than the flood-sensitive genotype. Roots of tolerant plants also contained more ATP than roots of sensitive plants at the 7th and 10th days of flooding. Quantitative transcript analysis identified 132 genes differentially expressed between the two genotypes at one or more time points of flooding. Expression of genes related to the ethylene biosynthesis pathway and formation of adventitious roots was induced earlier and to higher levels in roots of the flood-tolerant genotype. Three potential flood-tolerance TFs which were differentially expressed between the two genotypes during the entire 10-day flooding duration were identified. This study confirmed the expression of anaerobic genes in response to soil flooding. Additionally, the differential expression of TFs associated with soil flooding tolerance was not qualitative but quantitative and temporal. Functional analyses of these genes will be necessary to reveal their potential to enhance flooding tolerance of soybean cultivars.

Genome sequence of the palaeopolyploid soybean.

Soybean (Glycine max) is one of the most important crop plants for seed protein and oil content, and for its capacity to fix atmospheric nitrogen through symbioses with soil-borne microorganisms. We sequenced the 1.1-gigabase genome by a whole-genome shotgun approach and integrated it with physical and high-density genetic maps to create a chromosome-scale draft sequence assembly. We predict 46,430 protein-coding genes, 70% more than Arabidopsis and similar to the poplar genome which, like soybean, is an ancient polyploid (palaeopolyploid). About 78% of the predicted genes occur in chromosome ends, which comprise less than one-half of the genome but account for nearly all of the genetic recombination. Genome duplications occurred at approximately 59 and 13 million years ago, resulting in a highly duplicated genome with nearly 75% of the genes present in multiple copies. The two duplication events were followed by gene diversification and loss, and numerous chromosome rearrangements. An accurate soybean genome sequence will facilitate the identification of the genetic basis of many soybean traits, and accelerate the creation of improved soybean varieties.

Regulation of growth response to water stress in the soybean primary root. I. Proteomic analysis reveals region-specific regulation of phenylpropanoid metabolism and control of free iron in the elongation zone.

In water-stressed soybean primary roots, elongation was maintained at well-watered rates in the apical 4 mm (region 1), but was progressively inhibited in the 4-8 mm region (region 2), which exhibits maximum elongation in well-watered roots. These responses are similar to previous results for the maize primary root. To understand these responses in soybean, spatial profiles of soluble protein composition were analysed. Among the changes, the results indicate that region-specific regulation of phenylpropanoid metabolism may contribute to the distinct growth responses in the different regions. Several enzymes related to isoflavonoid biosynthesis increased in abundance in region 1, correlating with a substantial increase of isoflavonoid content in this region which could contribute to growth maintenance via various potential mechanisms. In contrast, caffeoyl-CoA O-methyltransferase, which is involved in lignin synthesis, was highly up-regulated in region 2. This response was associated with enhanced accumulation of lignin, which may be related to the inhibition of growth in this region. Several proteins that increased in abundance in both regions of water-stressed roots were related to protection from oxidative damage. In particular, an increase in the abundance of ferritin proteins effectively sequestered more iron and prevented excess free iron in the elongation zone under water stress.