Genetic interactions regulating seed phytate and oligosaccharides in soybean (Glycine max L.).

Two low-phytate soybean (Glycine max (L.) Merr.) mutant lines- V99-5089 (mips mutation on chromosome 11) and CX-1834 (mrp-l and mrp-n mutations on chromosomes 19 and 3, respectively) have proven to be valuable resources for breeding of low-phytate, high-sucrose, and low-raffinosaccharide soybeans, traits that are highly desirable from a nutritional and environmental standpoint. A recombinant inbred population derived from the cross CX1834 x V99-5089 provides an opportunity to study the effect of different combinations of these three mutations on soybean phytate and oligosaccharides levels. Of the 173 recombinant inbred lines tested, 163 lines were homozygous for various combinations of MIPS and two MRP loci alleles. These individuals were grouped into eight genotypic classes based on the combination of SNP alleles at the three mutant loci. The two genotypic classes that were homozygous mrp-l/mrp-n and either homozygous wild-type or mutant at the mips locus (MIPS/mrp-l/mrp-n or mips/mrp-l/mrp-n) displayed relatively similar ~55% reductions in seed phytate, 6.94 mg g -1 and 6.70 mg g-1 respectively, as compared with 15.2 mg g-1 in the wild-type MIPS/MRP-L/MRP-N seed. Therefore, in the presence of the double mutant mrp-l/mrp-n, the mips mutation did not cause a substantially greater decrease in seed phytate level. However, the nutritionally-desirable high-sucrose/low-stachyose/low-raffinose seed phenotype originally observed in soybeans homozygous for the mips allele was reversed in the presence of mrp-l/mrp-n mutations: homozygous mips/mrp-l/mrp-n seed displayed low-sucrose (7.70%), high-stachyose (4.18%), and the highest observed raffinose (0.94%) contents per gram of dry seed. Perhaps the block in phytic acid transport from its cytoplasmic synthesis site to its storage site, conditioned by mrp-l/mrp-n, alters myo-inositol flux in mips seeds in a way that restores to wild-type levels the mips conditioned reductions in raffinosaccharides. Overall this study determined the combinatorial effects of three low phytic acid causing mutations on regulation of seed phytate and oligosaccharides in soybean.

A Method for Combining Isolates of Phytophthora sojae to Screen for Novel Sources of Resistance to Phytophthora Stem and Root Rot in Soybean.

Soybean cultivars with specific single resistance genes (Rps) are grown to reduce yield loss due to Phytophthora stem and root rot caused by the oomycete pathogen Phytophthora sojae. To identify novel Rps loci, soybean lines are often screened several times, each time with an isolate of P. sojae that differs in virulence on various Rps genes. The goal of this study was to determine whether several isolates of P. sojae that differ in virulence on Rps genes could be combined into a single source of inoculum and used to screen soybean lines for novel Rps genes. A set of 14 soybean differential lines, each carrying a specific Rps gene, was inoculated with three isolates of P. sojae, which differed in virulence on 6 to 10 Rps genes, individually or in a 1:1:1 mixture. Inoculum containing the 1:1:1 mixture of isolates was virulent on 13 Rps genes. The mixed-inoculum method was used to screen 1,019 soybean accessions in a blind assay for novel sources of resistance. In all, 17% of Glycine max accessions and 11% of G. soja accessions were resistant (≤30% dead plants), suggesting that these accessions may carry a novel Rps gene or genes. Advantages of combining isolates into a single source of inoculum include reduced cost, ability to screen soybean germplasm with inoculum virulent on all known Rps genes, and ease of identifying novel sources of resistance. This study is a precursor to identifying novel sources of resistance to P. sojae in soybean using RXLR effectors.

Genome-wide transcriptome analyses of developing seeds from low and normal phytic acid soybean lines.

BACKGROUND: Low phytic acid (lpa) crops are potentially eco-friendly alternative to conventional normal phytic acid (PA) crops, improving mineral bioavailability in monogastric animals as well as decreasing phosphate pollution. The lpa crops developed to date carry mutations that are directly or indirectly associated with PA biosynthesis and accumulation during seed development. These lpa crops typically exhibit altered carbohydrate profiles, increased free phosphate, and lower seedling emergence, the latter of which reduces overall crop yield, hence limiting their large-scale cultivation. Improving lpa crop yield requires an understanding of the downstream effects of the lpa genotype on seed development. Towards that end, we present a comprehensive comparison of gene-expression profiles between lpa and normal PA soybean lines (Glycine max) at five stages of seed development using RNA-Seq approaches. The lpa line used in this study carries single point mutations in a myo-inositol phosphate synthase gene along with two multidrug-resistance protein ABC transporter genes. RESULTS: RNA sequencing data of lpa and normal PA soybean lines from five seed-developmental stages (total of 30 libraries) were used for differential expression and functional enrichment analyses. A total of 4235 differentially expressed genes, including 512-transcription factor genes were identified. Eighteen biological processes such as apoptosis, glucan metabolism, cellular transport, photosynthesis and 9 transcription factor families including WRKY, CAMTA3 and SNF2 were enriched during seed development. Genes associated with apoptosis, glucan metabolism, and cellular transport showed enhanced expression in early stages of lpa seed development, while those associated with photosynthesis showed decreased expression in late developmental stages. The results suggest that lpa-causing mutations play a role in inducing and suppressing plant defense responses during early and late stages of seed development, respectively. CONCLUSIONS: This study provides a global perspective of transcriptomal changes during soybean seed development in an lpa mutant. The mutants are characterized by earlier expression of genes associated with cell wall biosynthesis and a decrease in photosynthetic genes in late stages. The biological processes and transcription factors identified in this study are signatures of lpa-causing mutations.

Genetic variants in root architecture-related genes in a Glycine soja accession, a potential resource to improve cultivated soybean.

BACKGROUND: Root system architecture is important for water acquisition and nutrient acquisition for all crops. In soybean breeding programs, wild soybean alleles have been used successfully to enhance yield and seed composition traits, but have never been investigated to improve root system architecture. Therefore, in this study, high-density single-feature polymorphic markers and simple sequence repeats were used to map quantitative trait loci (QTLs) governing root system architecture in an inter-specific soybean mapping population developed from a cross between Glycine max and Glycine soja. RESULTS: Wild and cultivated soybean both contributed alleles towards significant additive large effect QTLs on chromosome 6 and 7 for a longer total root length and root distribution, respectively. Epistatic effect QTLs were also identified for taproot length, average diameter, and root distribution. These root traits will influence the water and nutrient uptake in soybean. Two cell division-related genes (D type cyclin and auxin efflux carrier protein) with insertion/deletion variations might contribute to the shorter root phenotypes observed in G. soja compared with cultivated soybean. Based on the location of the QTLs and sequence information from a second G. soja accession, three genes (slow anion channel associated 1 like, Auxin responsive NEDD8-activating complex and peroxidase), each with a non-synonymous single nucleotide polymorphism mutation were identified, which may also contribute to changes in root architecture in the cultivated soybean. In addition, Apoptosis inhibitor 5-like on chromosome 7 and slow anion channel associated 1-like on chromosome 15 had epistatic interactions for taproot length QTLs in soybean. CONCLUSION: Rare alleles from a G. soja accession are expected to enhance our understanding of the genetic components involved in root architecture traits, and could be combined to improve root system and drought adaptation in soybean.

Physiological and transcriptomic characterization of submergence and reoxygenation responses in soybean seedlings.

Complete inundation at the early seedling stage is a common environmental constraint for soybean production throughout the world. As floodwaters subside, submerged seedlings are subsequently exposed to reoxygenation stress in the natural progression of a flood event. Here, we characterized the fundamental acclimation responses to submergence and reoxygenation in soybean at the seedling establishment stage. Approximately 90% of seedlings succumbed during 3 d of inundation under constant darkness, whereas 10 d of submergence were lethal to over 90% of seedlings under 12 h light/12 h dark cycles, indicating the significance of underwater photosynthesis in seedling survival. Submergence rapidly decreased the abundance of carbohydrate reserves and ATP in aerial tissue of seedlings although chlorophyll breakdown was not observed. The carbohydrate and ATP contents were recovered upon de-submergence, but sudden exposure to oxygen also induced lipid peroxidation, confirming that reoxygenation induced oxidative stress. Whole transcriptome analysis recognized genome-scale reconfiguration of gene expression that regulates various signalling and metabolic pathways under submergence and reoxygenation. Comparative analysis of differentially regulated genes in shoots and roots of soybean and other plants defines conserved, organ-specific and species-specific adjustments which enhance adaptability to submergence and reoxygenation through different metabolic pathways.

Amino acid changes in P3, and not the overlapping pipo-encoded protein, determine virulence of soybean mosaic virus on functionally immune Rsv1-genotype soybean.

A small open reading frame, termed 'pipo', is embedded in the P3 cistron of potyviruses. Currently, knowledge on pipo and its role(s) in the life cycle of potyviruses is limited. The P3 and helper-component proteinase (HC-Pro) cistrons of Soybean mosaic virus (SMV) harbour determinants affecting virulence on functionally immune Rsv1-genotype soybeans. Interestingly, a key virulence determinant of SMV on Rsv1-genotype soybeans (i.e. soybeans containing the Rsv1 resistance gene) that resides at polyprotein codon 947 overlaps both P3 and a pipo-encoded codon. This raises the question of whether PIPO or P3 is the virulence factor. To answer this question, the corresponding pipo of an avirulent and two virulent strains of SMV were studied by comparative genomics, followed by syntheses and analyses of site-directed mutants. Our data demonstrate that the virulence of SMV on Rsv1-genotype soybeans is affected by P3 and not the overlapping pipo-encoded protein.

Experimental adaptation of an RNA virus mimics natural evolution.

Identification of virulence determinants of viruses is of critical importance in virology. In search of such determinants, virologists traditionally utilize comparative genomics between a virulent and an avirulent virus strain and construct chimeras to map their locations. Subsequent comparison reveals sequence differences, and through analyses of site-directed mutants, key residues are identified. In the absence of a naturally occurring virulent strain, an avirulent strain can be functionally converted to a virulent variant via an experimental evolutionary approach. However, the concern remains whether experimentally evolved virulence determinants mimic those that have evolved naturally. To provide a direct comparison, we exploited a plant RNA virus, soybean mosaic virus (SMV), and its natural host, soybean. Through a serial in vivo passage experiment, the molecularly cloned genome of an avirulent SMV strain was converted to virulent variants on functionally immune soybean genotypes harboring resistance factor(s) from the complex Rsv1 locus. Several of the experimentally evolved virulence determinants were identical to those discovered through a comparative genomic approach with a naturally evolved virulent strain. Thus, our observations validate an experimental evolutionary approach to identify relevant virulence determinants of an RNA virus.

Infection and genotype remodel the entire soybean transcriptome.

BACKGROUND: High throughput methods, such as high density oligonucleotide microarray measurements of mRNA levels, are popular and critical to genome scale analysis and systems biology. However understanding the results of these analyses and in particular understanding the very wide range of levels of transcriptional changes observed is still a significant challenge. Many researchers still use an arbitrary cut off such as two-fold in order to identify changes that may be biologically significant. We have used a very large-scale microarray experiment involving 72 biological replicates to analyze the response of soybean plants to infection by the pathogen Phytophthora sojae and to analyze transcriptional modulation as a result of genotypic variation. RESULTS: With the unprecedented level of statistical sensitivity provided by the high degree of replication, we show unambiguously that almost the entire plant genome (97 to 99% of all detectable genes) undergoes transcriptional modulation in response to infection and genetic variation. The majority of the transcriptional differences are less than two-fold in magnitude. We show that low amplitude modulation of gene expression (less than two-fold changes) is highly statistically significant and consistent across biological replicates, even for modulations of less than 20%. Our results are consistent through two different normalization methods and two different statistical analysis procedures. CONCLUSION: Our findings demonstrate that the entire plant genome undergoes transcriptional modulation in response to infection and genetic variation. The pervasive low-magnitude remodeling of the transcriptome may be an integral component of physiological adaptation in soybean, and in all eukaryotes.

Differential accumulation of retroelements and diversification of NB-LRR disease resistance genes in duplicated regions following polyploidy in the ancestor of soybean.

The genomes of most, if not all, flowering plants have undergone whole genome duplication events during their evolution. The impact of such polyploidy events is poorly understood, as is the fate of most duplicated genes. We sequenced an approximately 1 million-bp region in soybean (Glycine max) centered on the Rpg1-b disease resistance gene and compared this region with a region duplicated 10 to 14 million years ago. These two regions were also compared with homologous regions in several related legume species (a second soybean genotype, Glycine tomentella, Phaseolus vulgaris, and Medicago truncatula), which enabled us to determine how each of the duplicated regions (homoeologues) in soybean has changed following polyploidy. The biggest change was in retroelement content, with homoeologue 2 having expanded to 3-fold the size of homoeologue 1. Despite this accumulation of retroelements, over 77% of the duplicated low-copy genes have been retained in the same order and appear to be functional. This finding contrasts with recent analyses of the maize (Zea mays) genome, in which only about one-third of duplicated genes appear to have been retained over a similar time period. Fluorescent in situ hybridization revealed that the homoeologue 2 region is located very near a centromere. Thus, pericentromeric localization, per se, does not result in a high rate of gene inactivation, despite greatly accelerated retrotransposon accumulation. In contrast to low-copy genes, nucleotide-binding-leucine-rich repeat disease resistance gene clusters have undergone dramatic species/homoeologue-specific duplications and losses, with some evidence for partitioning of subfamilies between homoeologues.

Replication of nonautonomous retroelements in soybean appears to be both recent and common.

Retrotransposons and their remnants often constitute more than 50% of higher plant genomes. Although extensively studied in monocot crops such as maize (Zea mays) and rice (Oryza sativa), the impact of retrotransposons on dicot crop genomes is not well documented. Here, we present an analysis of retrotransposons in soybean (Glycine max). Analysis of approximately 3.7 megabases (Mb) of genomic sequence, including 0.87 Mb of pericentromeric sequence, uncovered 45 intact long terminal repeat (LTR)-retrotransposons. The ratio of intact elements to solo LTRs was 8:1, one of the highest reported to date in plants, suggesting that removal of retrotransposons by homologous recombination between LTRs is occurring more slowly in soybean than in previously characterized plant species. Analysis of paired LTR sequences uncovered a low frequency of deletions relative to base substitutions, indicating that removal of retrotransposon sequences by illegitimate recombination is also operating more slowly. Significantly, we identified three subfamilies of nonautonomous elements that have replicated in the recent past, suggesting that retrotransposition can be catalyzed in trans by autonomous elements elsewhere in the genome. Analysis of 1.6 Mb of sequence from Glycine tomentella, a wild perennial relative of soybean, uncovered 23 intact retroelements, two of which had accumulated no mutations in their LTRs, indicating very recent insertion. A similar pattern was found in 0.94 Mb of sequence from Phaseolus vulgaris (common bean). Thus, autonomous and nonautonomous retrotransposons appear to be both abundant and active in Glycine and Phaseolus. The impact of nonautonomous retrotransposon replication on genome size appears to be much greater than previously appreciated.

Application of comparative genomics in developing molecular markers tightly linked to the virus resistance gene Rsv4 in soybean.

The Rsv4 gene confers resistance to all the known strain groups of soybean mosaic virus in soybean (Glycine max (L.) Merr.). To construct a fine genetic map near Rsv4 in soybean, we employed a comparative genomics approach that used genome sequence information of the model legume Lotus japonicus. Sequences of the soybean expressed sequence tags (ESTs) AI856415 and BF070293 mapping to one side of the Rsv4 gene showed high similarity with gene sequences of the transformation-competent artificial chromosome (TAC) clone LjT32P24 of Lotus. LjT32P24 is tightly linked to another sequenced TAC clone, LjT26I01, in Lotus. A new marker, AW307114A, developed from soybean EST AW307114, which is homologous to a Lotus gene within LjT26I01, was mapped to the other side of the Rsv4 gene. The identification of the microsyntenic relationship facilitated the development of additional 2 EST markers between BF070293-S and AW307114A bracketing the Rsv4 gene. Several other markers developed in this study were mapped to putative homoeologous or duplicated chromosomal regions in soybean. Alignment between the soybean maps indicated that Rsv4 is located near a local chromosomal rearrangement. This targeted comparative mapping serves to provide a foundation for marker-assisted selection and cloning of the Rsv4 gene.

High-resolution mapping of the barley leaf rust resistance gene Rph5 using barley expressed sequence tags (ESTs) and synteny with rice.

The rapidly growing expressed sequence tag (EST) resources of species representing the Poacea family and availability of comprehensive sequence information for the rice (Oryza sativa) genome create an excellent opportunity for comparative genome analysis. Extensive synteny between rice chromosome 1 and barley (Hordeum vulgare L.) chromosome 3 has proven extremely useful for saturation mapping of chromosomal regions containing target genes of large-genome barley with conserved orthologous genes from the syntenic regions of the rice genome. Rph5 is a gene conferring resistance to the barley leaf rust pathogen Puccinia hordei. It was mapped to chromosome 3HS, which is syntenic with rice chromosome 1S. The objective of this study was to increase marker density within the sub-centimorgan region around Rph5, using sequence-tagged site (STS) markers that were developed based on barley ESTs syntenic to the phage (P1)-derived artificial chromosome (PAC) clones comprising the distal region of rice chromosome 1S. Five rice PAC clones were used as queries in a blastn search to screen 375,187 barley ESTs. Ninety-four non-redundant EST sequences were identified from the EST database and used as templates to design 174 pairs of primer combinations. As a result, 9 barley EST-based STS markers were incorporated into the 'Bowman' x 'Magnif 102' high-resolution map of the Rph5 region. More importantly, six markers, including five EST-derived STS sequences, were found to co-segregate with Rph5. The results of this study demonstrate the usefulness of rice genomic resources for efficient deployment of barley ESTs for marker saturation of targeted barley genomic regions.

Recombination within a nucleotide-binding-site/leucine-rich-repeat gene cluster produces new variants conditioning resistance to soybean mosaic virus in soybeans.

The soybean Rsv1 gene for resistance to soybean mosaic virus (SMV; Potyvirus) has previously been described as a single-locus multi-allelic gene mapping to molecular linkage group (MLG) F. Various Rsv1 alleles condition different responses to the seven (G1-G7) described strains of SMV, including extreme resistance, localized and systemic necrosis, and mosaic symptoms. We describe the cloning of a cluster of NBS-LRR resistance gene candidates from MLG F of the virus-resistant soybean line PI96983 and demonstrate that multiple genes within this cluster interact to condition unique responses to SMV strains. In addition to cloning 3gG2, a strong candidate for the major Rsv1 resistance gene from PI96983, we describe various unique resistant and necrotic reactions coincident with the presence or absence of other members of this gene cluster. Responses of recombinant lines from a high-resolution mapping population of PI96983 (resistant) x Lee 68 (susceptible) demonstrate that more than one gene in this region of the PI96983 chromosome conditions resistance and/or necrosis to SMV. In addition, the soybean cultivars Marshall and Ogden, which carry other previously described Rsv1 alleles, are shown to possess the 3gG2 gene in a NBS-LRR gene cluster background distinct from PI96983. These observations suggest that two or more related non-TIR-NBS-LRR gene products are likely involved in the allelic response of several Rsv1-containing lines to SMV.

Genetic and Sequence Analysis of Markers Tightly Linked to the Soybean mosaic virus Resistance Gene, Rsv3.

Soybean mosaic virus (SMV) is a major viral pathogen, affecting soybean [Glycine max (L.) Merr.] production worldwide. The Rsv3 gene of soybean confers resistance to three of the most virulent strains (G5-G7) of SMV. The objectives of this study were to map Rsv3 and develop polymerase chain reaction (PCR) based markers for marker-assisted selection (MAS) purposes. Disease-response data were collected from two F(2) mapping populations, L29 (Rsv3) x Lee68 (rsv3) and Tousan 140 (Rsv3) x Lee68 (rsv3). Bulk segregant analysis based on amplified fragment length polymorphism (AFLP) markers demonstrated that the Rsv3 locus maps to the soybean molecular linkage group (MLG) B2 between restriction fragment length polymorphism (RFLP) markers A519 and Mng247. These two tightly linked RFLP markers were converted to PCR-based markers to expedite MAS. Sequence analysis of the Mng247 genomic region revealed similarity to the consensus sequence of a leucine-rich repeat (LRR) characteristic of the extracellular LRR class of disease resistance genes. Results from this study will be useful in pyramiding viral resistance genes and in cloning the Rsv3 gene.