Large-scale data mining pipeline for identifying novel soybean genes involved in resistance against the soybean cyst nematode

The soybean cyst nematode (SCN) [Heterodera glycines Ichinohe] is a devastating pathogen of soybean [Glycine max (L.) Merr.] that is rapidly becoming a global economic issue. Two loci conferring SCN resistance have been identified in soybean, Rhg1 and Rhg4; however, they offer declining protection. Therefore, it is imperative that we identify additional mechanisms for SCN resistance. In this paper, we develop a bioinformatics pipeline to identify protein-protein interactions related to SCN resistance by data mining massive-scale datasets. The pipeline combines two leading sequence-based protein-protein interaction predictors, the Protein-protein Interaction Prediction Engine (PIPE), PIPE4, and Scoring PRotein INTeractions (SPRINT) to predict high-confidence interactomes. First, we predicted the top soy interacting protein partners of the Rhg1 and Rhg4 proteins. Both PIPE4 and SPRINT overlap in their predictions with 58 soybean interacting partners, 19 of which had GO terms related to defense. Beginning with the top predicted interactors of Rhg1 and Rhg4, we implement a "guilt by association" in silico proteome-wide approach to identify novel soybean genes that may be involved in SCN resistance. This pipeline identified 1,082 candidate genes whose local interactomes overlap significantly with the Rhg1 and Rhg4 interactomes. Using GO enrichment tools, we highlighted many important genes including five genes with GO terms related to response to the nematode (GO:0009624), namely, Glyma.18G029000, Glyma.11G228300, Glyma.08G120500, Glyma.17G152300, and Glyma.08G265700. This study is the first of its kind to predict interacting partners of known resistance proteins Rhg1 and Rhg4, forming an analysis pipeline that enables researchers to focus their search on high-confidence targets to identify novel SCN resistance genes in soybean.

The SoyaGen Project: Putting Genomics to Work for Soybean Breeders

The SoyaGen project was a collaborative endeavor involving Canadian soybean researchers and breeders from academia and the private sector as well as international collaborators. Its aims were to develop genomics-derived solutions to real-world challenges faced by breeders. Based on the needs expressed by the stakeholders, the research efforts were focused on maximizing realized yield through optimization of maturity and improved disease resistance. The main deliverables related to molecular breeding in soybean will be reviewed here. These include: (1) SNP datasets capturing the genetic diversity within cultivated soybean (both within a worldwide collection of > 1,000 soybean accessions and a subset of 102 short-season accessions (MG0 and earlier) directly relevant to this group); (2) SNP markers for selecting favorable alleles at key maturity genes as well as loci associated with increased resistance to key pathogens and pests (Phytophthora sojae, Heterodera glycines, Sclerotinia sclerotiorum); (3) diagnostic tools to facilitate the identification and mapping of specific pathotypes of P. sojae; and (4) a genomic prediction approach to identify the most promising combinations of parents. As a result of this fruitful collaboration, breeders have gained new tools and approaches to implement molecular, genomics-informed breeding strategies. We believe these tools and approaches are broadly applicable to soybean breeding efforts around the world.

Soybean (Glycine max) Haplotype Map (GmHapMap): a universal resource for soybean translational and functional genomics

Here, we describe a worldwide haplotype map for soybean (GmHapMap) constructed using whole-genome sequence data for 1007 Glycine max accessions and yielding 14.9 million variants as well as 4.3 M tag single-nucleotide polymorphisms (SNPs). When sampling random subsets of these accessions, the number of variants and tag SNPs plateaued beyond approximately 800 and 600 accessions, respectively. This suggests extensive coverage of diversity within the cultivated soybean. GmHapMap variants were imputed onto 21 618 previously genotyped accessions with up to 96% success for common alleles. A local association analysis was performed with the imputed data using markers located in a 1-Mb region known to contribute to seed oil content and enabled us to identify a candidate causal SNP residing in the NPC1 gene. We determined gene-centric haplotypes (407 867 GCHs) for the 55 589 genes and showed that such haplotypes can help to identify alleles that differ in the resulting phenotype. Finally, we predicted 18 031 putative loss-of-function (LOF) mutations in 10 662 genes and illustrated how such a resource can be used to explore gene function. The GmHapMap provides a unique worldwide resource for applied soybean genomics and breeding.

Soybean (Glycine max) Haplotype Map (GmHapMap): a universal resource for soybean translational and functional genomics

Here, we describe a worldwide haplotype map for soybean (GmHapMap) constructed using whole-genome sequence data for 1007 Glycine max accessions and yielding 14.9 million variants as well as 4.3 M tag single-nucleotide polymorphisms (SNPs). When sampling random subsets of these accessions, the number of variants and tag SNPs plateaued beyond approximately 800 and 600 accessions, respectively. This suggests extensive coverage of diversity within the cultivated soybean. GmHapMap variants were imputed onto 21 618 previously genotyped accessions with up to 96% success for common alleles. A local association analysis was performed with the imputed data using markers located in a 1-Mb region known to contribute to seed oil content and enabled us to identify a candidate causal SNP residing in the NPC1 gene. We determined gene-centric haplotypes (407 867 GCHs) for the 55 589 genes and showed that such haplotypes can help to identify alleles that differ in the resulting phenotype. Finally, we predicted 18 031 putative loss-of-function (LOF) mutations in 10 662 genes and illustrated how such a resource can be used to explore gene function. The GmHapMap provides a unique worldwide resource for applied soybean genomics and breeding.

The effects and interaction of soybean maturity gene alleles controlling flowering time, maturity, and adaptation in tropical environments

BACKGROUND

Soybean is native to the temperate zones of East Asia. Poor yields of soybean in West African countries may be partially attributed to inadequate adaptation of soybean to tropical environments. Adaptation will require knowledge of the effects of allelic combinations of major maturity genes (E1, E2, and E3) and stem architecture. The long juvenile trait (J) influences soybean flowering time in short, ~ 12 h days, which characterize tropical latitudes. Soybean plant architecture includes determinate or indeterminate stem phenotypes controlled by the Dt1 gene. Understanding the influence of these genetic components on plant development and adaptation is key to optimize phenology and improve soybean yield potential in tropical environments.

RESULTS

Soybean lines from five recombinant inbred populations were developed that varied in their combinations of targeted genes. The soybean lines were field tested in multiple environments and characterized for days to flowering (DTF), days to maturity (DTM), and plant height in locations throughout northern Ghana, and allelic combinations were determined for each line for associating genotype with phenotype. The results revealed significant differences based on genotype for DTF and DTM and allowed the comparison of different variant alleles of those genes. The mutant alleles of J and E1 had significant impact on DTF and DTM, and alleles of those genes interacted with each other for DTF but not DTM. The Dt1 gene significantly influenced plant height but not DTF or DTM.

CONCLUSIONS

This research identified major and minor effect alleles of soybean genes that can be combined to control DTF, DTM, and plant height in short day tropical environments in Ghana. These phenotypes contribute to adaptation to a low latitude environment that can be optimized in a soybean breeding program with targeted selection of desired allele combinations. The knowledge of the genetic control of these traits will enhance molecular breeding to produce optimally adapted soybean varieties targeted to tropical environments.

Effect of change from a determinate to a semi-determinate growth habit on the yield and lodging resistance of soybeans in the northeast region of Japan

Although an indeterminate growth habit is attractive to develop high-yield soybean varieties with higher number of pods (Glycine max (L). Merr.), lodging in indeterminate varieties remains a problem in Japan. As the semi-determinate varieties have shorter main stem length than the indeterminate varieties, this trait can be useful to improve varieties with high yield and low lodging risk. We introduced the genes Dt1 and Dt2, which regulate stem growth habit, into three determinate varieties by backcrossing and evaluated the resulting effects on yield and lodging tendency under four different growing environments. The yield and lodging degree of the semi-determinate and indeterminate lines were higher and more severe than those of the determinate lines. Despite the lower overall lodging score, the semi-determinate lines had marginally lower overall yield than that of the indeterminate lines. However, the effect of introduction of semi-determinate traits on yield and lodging degree was different in the three backgrounds, with the yield of semi-determinate lines being the highest and the difference in lodging degree between the semi-determinate and determinate lines being under 1.0 in one background. Therefore, semi-determinate growth habit has potential to develop high yielding varieties with low lodging risk.

Fatal Fusarium infection manifesting as osteomyelitis following previous treatment with amphotericin B in a multi-visceral transplant: Case report and review of Fusarium infections in solid organ transplantation

Invasive fusariosis in solid organ transplant is uncommon and usually presents as localized infection with favorable outcomes compared to hematologic malignancies or bone marrow transplants. We report the first case of Fusarium osteomyelitis in a patient following multi-visceral transplant and review Fusarium in organ transplant recipients and Fusarium bone and joint infections. Our case underscores the importance of early recognition and multidisciplinary approach to treatment and highlights potential failure to eradicate with amphotericin B monotherapy.

Comprehensive description of genomewide nucleotide and structural variation in short-season soya bean

Next-generation sequencing (NGS) and bioinformatics tools have greatly facilitated the characterization of nucleotide variation; nonetheless, an exhaustive description of both SNP haplotype diversity and of structural variation remains elusive in most species. In this study, we sequenced a representative set of 102 short-season soya beans and achieved an extensive coverage of both nucleotide diversity and structural variation (SV). We called close to 5M sequence variants (SNPs, MNPs and indels) and noticed that the number of unique haplotypes had plateaued within this set of germplasm (1.7M tag SNPs). This data set proved highly accurate (98.6%) based on a comparison of called genotypes at loci shared with a SNP array. We used this catalogue of SNPs as a reference panel to impute missing genotypes at untyped loci in data sets derived from lower density genotyping tools (150 K GBS-derived SNPs/530 samples). After imputation, 96.4% of the missing genotypes imputed in this fashion proved to be accurate. Using a combination of three bioinformatics pipelines, we uncovered ~92 K SVs (deletions, insertions, inversions, duplications, CNVs and translocations) and estimated that over 90% of these were accurate. Finally, we noticed that the duplication of certain genomic regions explained much of the residual heterozygosity at SNP loci in otherwise highly inbred soya bean accessions. This is the first time that a comprehensive description of both SNP haplotype diversity and SV has been achieved within a regionally relevant subset of a major crop.

Data on the epitope mapping of soybean A2 and A3 glycinin

The data information provided in this article relate to our research article “Using patient serum to epitope map soybean glycinins reveals common epitopes shared with many legumes and tree nuts” (Saeed et al., 2016) [1]. Here we provide western blot detection of glycinin subunits by soy-sensitive human sera, ELISA screens with overlapping synthetic peptides (epitope mapping), and various database/server epitope searches.

Using patient serum to epitope map soybean glycinins reveals common epitopes shared with many legumes and tree nuts

Soybean consumption is increasing in many Western diets; however, recent reviews suggest that the prevalence of soy allergy can be as high as 0.5% for the general population and up to 13% for children. The immunoglobulin-E (IgE) binding of sera from six soy-sensitive adult human subjects to soybean proteins separated by 2D gel electrophoresis was studied. Synthetic peptide sets spanning the mature glycinin subunit A2 and A3 primary sequences were used to map the IgE-binding regions. Putative epitopes identified in this study were also localized on glycinin hexamer models using bioinformatics software. We identified linear IgE-binding epitopes of the major storage protein Gly m 6 by screening individual soy-sensitive patient sera. These epitopes were then further analysed by 3D in silico model localization and compared to other plant storage protein epitopes. Web-based software applications were also used to study the ability to accurately predict epitopes with mixed results. A total of nine putative IgE-binding epitopes were identified in the glycinin A3 (A3.1-A3.3) and A2 (A2.1-A2.6) subunits. Most patients' sera IgE bound to only one or two epitopes, except for one patient's serum which bound to four different A2 epitopes. Two epitopes (A3.2 and A2.4) overlapped with a previously identified epitope hot spot of 11S globulins from other plant species. Most epitopes were predicted to be exposed on the surface of the 3D model of the glycinin hexamer. Amino acid sequence alignments of soybean acidic glycinins and other plant globulins revealed one dominant epitope hot spot among the four reported hot spots. This study may be helpful for future development of soy allergy immunotherapy and diagnosis.

Identification of loci governing eight agronomic traits using a GBS-GWAS approach and validation by QTL mapping in soya bean

Soya bean is a major source of edible oil and protein for human consumption as well as animal feed. Understanding the genetic basis of different traits in soya bean will provide important insights for improving breeding strategies for this crop. A genome-wide association study (GWAS) was conducted to accelerate molecular breeding for the improvement of agronomic traits in soya bean. A genotyping-by-sequencing (GBS) approach was used to provide dense genome-wide marker coverage (>47,000 SNPs) for a panel of 304 short-season soya bean lines. A subset of 139 lines, representative of the diversity among these, was characterized phenotypically for eight traits under six environments (3 sites × 2 years). Marker coverage proved sufficient to ensure highly significant associations between the genes known to control simple traits (flower, hilum and pubescence colour) and flanking SNPs. Between one and eight genomic loci associated with more complex traits (maturity, plant height, seed weight, seed oil and protein) were also identified. Importantly, most of these GWAS loci were located within genomic regions identified by previously reported quantitative trait locus (QTL) for these traits. In some cases, the reported QTLs were also successfully validated by additional QTL mapping in a biparental population. This study demonstrates that integrating GBS and GWAS can be used as a powerful complementary approach to classical biparental mapping for dissecting complex traits in soya bean.

First report of successful treatment of multidrug-resistant cytomegalovirus disease with the novel anti-CMV compound AIC246

We report the first case of cytomegalovirus (CMV) disease treated with AIC246, a novel anti-CMV compound which targets the viral terminase complex and remains active against virus resistant to DNA polymerase inhibitors. A lung transplant recipient developed refractory multidrug-resistant CMV disease involving the lungs, gastrointestinal tract and retina. His disease progressed despite treatment with all DNA polymerase inhibitors; multiple agents reported to have activity against CMV in case series, and reduction in his immunosuppressive medications. AIC246 which is in clinical development was obtained for emergency use, and combined with additional reduction in immunosuppression resulted in rapid clinical, virological and radiological resolution of disease. The patient has remained free of CMV disease or viremia off treatment for greater than 3 months. In summary AIC246, while still in development, may be a promising alternative to current therapies.

Mapping and validation of simple sequence repeat markers linked to a major gene controlling seed cadmium accumulation in soybean [Glycine max (L.) Merr]

Daily consumption of cadmium (Cd) contaminated foods poses a risk to human health. Cultivar selection is an important method to limit Cd uptake and accumulation; however, analyzing grain Cd concentration is costly and time-consuming. Developing markers for low Cd accumulation will facilitate marker assisted selection (MAS). Inheritance studies using a threshold value of 0.2 mg kg(-1) for low and high and an F(2:3) population showed that low Cd accumulation in soybean seed is under the control of a major gene (Cda1, proposed name) with the allele for low accumulation being dominant. A recombinant inbred line (RIL) population (F(6:8)) derived from the cross AC Hime (high Cd accumulation) and Westag-97 (low Cd accumulation) was used to identify the DNA markers linked to Cda gene(s) or quantitative trait loci (QTLs) controlling low Cd accumulation. We screened 171 simple sequence repeat (SSR) primers that showed polymorphism between parents on the 166 RILs. Of these, 40 primers were newly developed from the soybean genomic DNA sequence. Seven SSR markers, SatK138, SatK139, SatK140 (0.5 cM), SatK147, SacK149, SaatK150 and SattK152 (0.3 cM), were linked to Cda1 in soybean seed. All the linked markers were mapped to the same linkage group (LG) K. The closest flanking SSR markers linked to Cda1 were validated using a parallel population (RILs) involving Leo x Westag-97. Linked markers were also validated with diverse soybean genotypes differing in their seed Cd concentration and showed that SSR markers SatK147, SacK149, and SattK152 clearly differentiated the high and low Cd accumulating genotypes tested. To treat Cd uptake as a quantitative trait, QTL analysis using a linkage map constructed with 161 markers identified a major QTL associated with low Cd concentration in the seeds. The QTL was also mapped to the same location as Cda1 on LG-K. This QTL accounted for 57.3% of the phenotypic variation. Potential candidate genes (genes with known or predicted function that could influence the seed Cd concentration) like protein kinase, putative Adagio-like protein, and plasma membrane H(+)-ATPase were found to be located in the locus of interest. Of the four SSR markers located in the region, SattK152 was localized in the plasma membrane H(+)-ATPase gene. SSR markers closely linked to Cda1 in seeds of soybean were identified and have potential to be used for MAS to develop low Cd accumulating cultivars in a breeding program.