Identification of quantitative trait loci controlling soybean seed protein and oil content.

Soybean is a major source of seed protein and oil globally with an average composition of 40% protein and 20% oil in the seed. The goal of this study was to identify quantitative trait loci (QTL) conferring seed protein and oil content utilizing a population constructed by crossing an above average protein content line, PI 399084 to another line that had a low protein content value, PI 507429, both from the USDA soybean germplasm collection. The recombinant inbred line (RIL) population, PI 507429 x PI 399084, was evaluated in two replications over four years (2018-2021); the seeds were analyzed for seed protein and oil content using near-infrared reflectance spectroscopy. The recombinant inbred lines and the two parents were re-sequenced using genotyping by sequencing. A total of 12,761 molecular markers, which came from genotyping by sequencing, the SoySNP6k BeadChip and selected simple sequence repeat (SSR) markers from known protein QTL chromosomal regions were used for mapping. One QTL was identified on chromosome 2 explaining up to 56.8% of the variation for seed protein content and up to 43% for seed oil content. Another QTL identified on chromosome 15 explained up to 27.2% of the variation for seed protein and up to 41% of the variation for seed oil content. The protein and oil QTLs of this study and their associated molecular markers will be useful in breeding to improve nutritional quality in soybean.

Network Inference of Transcriptional Regulation in Germinating Low Phytic Acid Soybean Seeds.

The low phytic acid (lpa) trait in soybeans can be conferred by loss-of-function mutations in genes encoding myo-inositol phosphate synthase and two epistatically interacting genes encoding multidrug-resistance protein ATP-binding cassette (ABC) transporters. However, perturbations in phytic acid biosynthesis are associated with poor seed vigor. Since the benefits of the lpa trait, in terms of end-use quality and sustainability, far outweigh the negatives associated with poor seed performance, a fuller understanding of the molecular basis behind the negatives will assist crop breeders and engineers in producing variates with lpa and better germination rate. The gene regulatory network (GRN) for developing low and normal phytic acid soybean seeds was previously constructed, with genes modulating a variety of processes pertinent to phytic acid metabolism and seed viability being identified. In this study, a comparative time series analysis of low and normal phytic acid soybeans was carried out to investigate the transcriptional regulatory elements governing the transitional dynamics from dry seed to germinated seed. GRNs were reverse engineered from time series transcriptomic data of three distinct genotypic subsets composed of lpa soybean lines and their normal phytic acid sibling lines. Using a robust unsupervised network inference scheme, putative regulatory interactions were inferred for each subset of genotypes. These interactions were further validated by published regulatory interactions found in Arabidopsis thaliana and motif sequence analysis. Results indicate that lpa seeds have increased sensitivity to stress, which could be due to changes in phytic acid levels, disrupted inositol phosphate signaling, disrupted phosphate ion (Pi) homeostasis, and altered myo-inositol metabolism. Putative regulatory interactions were identified for the latter two processes. Changes in abscisic acid (ABA) signaling candidate transcription factors (TFs) putatively regulating genes in this process were identified as well. Analysis of the GRNs reveal altered regulation in processes that may be affecting the germination of lpa soybean seeds. Therefore, this work contributes to the ongoing effort to elucidate molecular mechanisms underlying altered seed viability, germination and field emergence of lpa crops, understanding of which is necessary in order to mitigate these problems.

Identification of Quantitative Disease Resistance Loci Toward Four Pythium Species in Soybean.

In this study, four recombinant inbred line (RIL) soybean populations were screened for their response to infection by Pythium sylvaticum, Pythium irregulare, Pythium oopapillum, and Pythium torulosum. The parents, PI 424237A, PI 424237B, PI 408097, and PI 408029, had higher levels of resistance to these species in a preliminary screening and were crossed with "Williams," a susceptible cultivar. A modified seed rot assay was used to evaluate RIL populations for their response to specific Pythium species selected for a particular population based on preliminary screenings. Over 2500 single-nucleotide polymorphism (SNP) markers were used to construct chromosomal maps to identify regions associated with resistance to Pythium species. Several minor and large effect quantitative disease resistance loci (QDRL) were identified including one large effect QDRL on chromosome 8 in the population of PI 408097 × Williams. It was identified by two different disease reaction traits in P. sylvaticum, P. irregulare, and P. torulosum. Another large effect QDRL was identified on chromosome 6 in the population of PI 408029 × Williams, and conferred resistance to P. sylvaticum and P. irregulare. These large effect QDRL will contribute toward the development of improved soybean cultivars with higher levels of resistance to these common soil-borne pathogens.

Mining germplasm panels and phenotypic datasets to identify loci for resistance to Phytophthora sojae in soybean.

Phytophthora sojae causes Phytophthora root and stem rot of soybean and has been primarily managed through deployment of qualitative Resistance to P. sojae genes (Rps genes). The effectiveness of each individual or combination of Rps gene(s) depends on the diversity and pathotypes of the P. sojae populations present. Due to the complex nature of P. sojae populations, identification of more novel Rps genes is needed. In this study, phenotypic data from previous studies of 16 panels of plant introductions (PIs) were analyzed. Panels 1 and 2 consisted of 448 Glycine max and 520 G. soja, which had been evaluated for Rps gene response with a combination of P. sojae isolates. Panels 3 and 4 consisted of 429 and 460 G. max PIs, respectively, which had been evaluated using individual P. sojae isolates with complex virulence pathotypes. Finally, Panels 5-16 (376 G. max PIs) consisted of data deposited in the USDA Soybean Germplasm Collection from evaluations with 12 races of P. sojae. Using these panels, genome-wide association (GWA) analyses were carried out by combining phenotypic and SoySNP50K genotypic data. GWA models identified two, two, six, and seven novel Rps loci with Panels 1, 2, 3, and 4, respectively. A total of 58 novel Rps loci were identified using Panels 5-16. Genetic and phenotypic dissection of these loci may lead to the characterization of novel Rps genes that can be effectively deployed in new soybean cultivars against diverse P. sojae populations.

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 transcriptional regulatory network of Rsv3-mediated extreme resistance against Soybean mosaic virus.

Resistance genes are an effective means for disease control in plants. They predominantly function by inducing a hypersensitive reaction, which results in localized cell death restricting pathogen spread. Some resistance genes elicit an atypical response, termed extreme resistance, where resistance is not associated with a hypersensitive reaction and its standard defense responses. Unlike hypersensitive reaction, the molecular regulatory mechanism(s) underlying extreme resistance is largely unexplored. One of the few known, naturally occurring, instances of extreme resistance is resistance derived from the soybean Rsv3 gene, which confers resistance against the most virulent Soybean mosaic virus strains. To discern the regulatory mechanism underlying Rsv3-mediated extreme resistance, we generated a gene regulatory network using transcriptomic data from time course comparisons of Soybean mosaic virus-G7-inoculated resistant (L29, Rsv3-genotype) and susceptible (Williams82, rsv3-genotype) soybean cultivars. Our results show Rsv3 begins mounting a defense by 6 hpi via a complex phytohormone network, where abscisic acid, cytokinin, jasmonic acid, and salicylic acid pathways are suppressed. We identified putative regulatory interactions between transcription factors and genes in phytohormone regulatory pathways, which is consistent with the demonstrated involvement of these pathways in Rsv3-mediated resistance. One such transcription factor identified as a putative transcriptional regulator was MYC2 encoded by Glyma.07G051500. Known as a master regulator of abscisic acid and jasmonic acid signaling, MYC2 specifically recognizes the G-box motif ("CACGTG"), which was significantly enriched in our data among differentially expressed genes implicated in abscisic acid- and jasmonic acid-related activities. This suggests an important role for Glyma.07G051500 in abscisic acid- and jasmonic acid-derived defense signaling in Rsv3. Resultantly, the findings from our network offer insights into genes and biological pathways underlying the molecular defense mechanism of Rsv3-mediated extreme resistance against Soybean mosaic virus. The computational pipeline used to reconstruct the gene regulatory network in this study is freely available at https://github.com/LiLabAtVT/rsv3-network.

Soybean mosaic virus: a successful potyvirus with a wide distribution but restricted natural host range.

TAXONOMY: Soybean mosaic virus (SMV) is a species within the genus Potyvirus, family Potyviridae, which includes almost one-quarter of all known plant RNA viruses affecting agriculturally important plants. The Potyvirus genus is the largest of all genera of plant RNA viruses with 160 species. PARTICLE: The filamentous particles of SMV, typical of potyviruses, are about 7500 Å long and 120 Å in diameter with a central hole of about 15 Å in diameter. Coat protein residues are arranged in helices of about 34 Å pitch having slightly less than nine subunits per turn. GENOME: The SMV genome consists of a single-stranded, positive-sense, polyadenylated RNA of approximately 9.6 kb with a virus-encoded protein (VPg) linked at the 5' terminus. The genomic RNA contains a single large open reading frame (ORF). The polypeptide produced from the large ORF is processed proteolytically by three viral-encoded proteinases to yield about 10 functional proteins. A small ORF, partially overlapping the P3 cistron, pipo, is encoded as a fusion protein in the N-terminus of P3 (P3N + PIPO). BIOLOGICAL PROPERTIES: SMV's host range is restricted mostly to two plant species of a single genus: Glycine max (cultivated soybean) and G. soja (wild soybean). SMV is transmitted by aphids non-persistently and by seeds. The variability of SMV is recognized by reactions on cultivars with dominant resistance (R) genes. Recessive resistance genes are not known. GEOGRAPHICAL DISTRIBUTION AND ECONOMIC IMPORTANCE: As a consequence of its seed transmissibility, SMV is present in all soybean-growing areas of the world. SMV infections can reduce significantly seed quantity and quality (e.g. mottled seed coats, reduced seed size and viability, and altered chemical composition). CONTROL: The most effective means of managing losses from SMV are the planting of virus-free seeds and cultivars containing single or multiple R genes. KEY ATTRACTIONS: The interactions of SMV with soybean genotypes containing different dominant R genes and an understanding of the functional role(s) of SMV-encoded proteins in virulence, transmission and pathogenicity have been investigated intensively. The SMV-soybean pathosystem has become an excellent model for the examination of the genetics and genomics of a uniquely complex gene-for-gene resistance model in a crop of worldwide importance.

Inference of Transcription Regulatory Network in Low Phytic Acid Soybean Seeds.

A dominant loss of function mutation in myo-inositol phosphate synthase (MIPS) gene and recessive loss of function mutations in two multidrug resistant protein type-ABC transporter genes not only reduce the seed phytic acid levels in soybean, but also affect the pathways associated with seed development, ultimately resulting in low emergence. To understand the regulatory mechanisms and identify key genes that intervene in the seed development process in low phytic acid crops, we performed computational inference of gene regulatory networks in low and normal phytic acid soybeans using a time course transcriptomic data and multiple network inference algorithms. We identified a set of putative candidate transcription factors and their regulatory interactions with genes that have functions in myo-inositol biosynthesis, auxin-ABA signaling, and seed dormancy. We evaluated the performance of our unsupervised network inference method by comparing the predicted regulatory network with published regulatory interactions in Arabidopsis. Some contrasting regulatory interactions were observed in low phytic acid mutants compared to non-mutant lines. These findings provide important hypotheses on expression regulation of myo-inositol metabolism and phytohormone signaling in developing low phytic acid soybeans. The computational pipeline used for unsupervised network learning in this study is provided as open source software and is freely available at https://lilabatvt.github.io/LPANetwork/.

Candidate Gene Sequence Analyses toward Identifying -Type Resistance to.

is one of three genetic loci conferring strain-specific resistance to (SMV). The locus has been mapped to a 154-kb region on chromosome 14, containing a cluster of five nucleotide-binding leucine-rich repeat (NB-LRR) resistance genes. High sequence similarity between the candidate genes challenges fine mapping of the locus. Among the five, Glyma14g38533 showed the highest transcript abundance in 1 to 3 h of SMV-G7 inoculation. Comparative sequence analyses were conducted with the five candidate NB-LRR genes from susceptible (-type) soybean [ (L.) Merr.] cultivar Williams 82, resistant (-type) cultivar Hwangkeum, and resistant lines L29 and RRR. Sequence comparisons revealed that Glyma14g38533 had far more polymorphisms than the other candidate genes. Interestingly, Glyma14g38533 gene from -type lines exhibited 150 single-nucleotide polymorphism (SNP and six insertion-deletion (InDel) markers relative to -type line, Furthermore, the polymorphisms identified in three -type lines were highly conserved. Several polymorphisms were validated in 18 -type resistant and six -type susceptible lines and were found associated with their disease response. The majority of the polymorphisms were located in LRR domain encoding region, which is involved in pathogen recognition via protein-protein interactions. These findings associating Glyma14g38533 with -type resistance to SMV suggest it is the most likely candidate gene for .

Identification of haplotypes at the Rsv4 genomic region in soybean associated with durable resistance to soybean mosaic virus.

KEY MESSAGE: Discovery of new germplasm sources and identification of haplotypes for the durable Soybean mosaic virus resistance gene, Rsv 4, provide novel resources for map-based cloning and genetic improvement efforts in soybean. The Soybean mosaic virus (SMV) resistance locus Rsv4 is of interest because it provides a durable type of resistance in soybean [Glycine max (L.) Merr.]. To better understand its molecular basis, we used a population of 309 BC3F2 individuals to fine-map Rsv4 to a ~120 kb interval and leveraged this genetic information in a second study to identify accessions 'Haman' and 'Ilpumgeomjeong' as new sources of Rsv4. These two accessions along with three other Rsv4 and 14 rsv4 accessions were used to examine the patterns of nucleotide diversity at the Rsv4 region based on high-depth resequencing data. Through a targeted association analysis of these 19 accessions within the ~120 kb interval, a cluster of four intergenic single-nucleotide polymorphisms (SNPs) was found to perfectly associate with SMV resistance. Interestingly, this ~120 kb interval did not contain any genes similar to previously characterized dominant disease resistance genes. Therefore, a haplotype analysis was used to further resolve the association signal to a ~94 kb region, which also resulted in the identification of at least two Rsv4 haplotypes. A haplotype phylogenetic analysis of this region suggests that the Rsv4 locus in G. max is recently introgressed from G. soja. This integrated study provides a strong foundation for efforts focused on the cloning of this durable virus resistance gene and marker-assisted selection of Rsv4-mediated SMV resistance in soybean breeding programs.

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.

The HC-Pro and P3 cistrons of an avirulent Soybean mosaic virus are recognized by different resistance genes at the complex Rsv1 locus.

The complex Rsv1 locus in soybean plant introduction (PI) 'PI96983' confers extreme resistance (ER) against Soybean mosaic virus (SMV) strain N but not SMV-G7 and SMV-G7d. Both the SMV helper-component proteinase (HC-Pro) and P3 cistrons can serve as avirulence factors recognized by Rsv1. To understand the genetics underlying recognition of the two cistrons, we have utilized two soybean lines (L800 and L943) derived from crosses between PI96983 (Rsv1) and Lee68 (rsv1) with distinct recombination events within the Rsv1 locus. L800 contains a single PI96983-derived member (3gG2) of an Rsv1-associated subfamily of nucleotide-binding leucine-rich repeat (NB-LRR) genes. In contrast, although L943 lacks 3gG2, it contains a suite of five other NB-LRR genes belonging to the same family. L800 confers ER against SMV-N whereas L943 allows limited replication at the inoculation site. SMV-N-derived chimeras containing HC-Pro from SMV-G7 or SMV-G7d gained virulence on L943 but not on L800 whereas those with P3 replacement gained virulence on L800 but not on L943. In reciprocal experiments, SMV-G7- and SMV-G7d-derived chimeras with HC-Pro replacement from SMV-N lost virulence on L943 but retained virulence on L800 whereas those with P3 replacement lost virulence on L800 while remaining virulent on L943. These data demonstrate that distinct resistance genes at the Rsv1 locus, likely belonging to the NB-LRR class, mediate recognition of HC-Pro and P3.

Evolution of a complex disease resistance gene cluster in diploid Phaseolus and tetraploid Glycine.

We used a comparative genomics approach to investigate the evolution of a complex nucleotide-binding (NB)-leucine-rich repeat (LRR) gene cluster found in soybean (Glycine max) and common bean (Phaseolus vulgaris) that is associated with several disease resistance (R) genes of known function, including Rpg1b (for Resistance to Pseudomonas glycinea1b), an R gene effective against specific races of bacterial blight. Analysis of domains revealed that the amino-terminal coiled-coil (CC) domain, central nucleotide-binding domain (NB-ARC [for APAF1, Resistance genes, and CED4]), and carboxyl-terminal LRR domain have undergone distinct evolutionary paths. Sequence exchanges within the NB-ARC domain were rare. In contrast, interparalogue exchanges involving the CC and LRR domains were common, consistent with both of these regions coevolving with pathogens. Residues under positive selection were overrepresented within the predicted solvent-exposed face of the LRR domain, although several also were detected within the CC and NB-ARC domains. Superimposition of these latter residues onto predicted tertiary structures revealed that the majority are located on the surface, suggestive of a role in interactions with other domains or proteins. Following polyploidy in the Glycine lineage, NB-LRR genes have been preferentially lost from one of the duplicated chromosomes (homeologues found in soybean), and there has been partitioning of NB-LRR clades between the two homeologues. The single orthologous region in common bean contains approximately the same number of paralogues as found in the two soybean homeologues combined. We conclude that while polyploidization in Glycine has not driven a stable increase in family size for NB-LRR genes, it has generated two recombinationally isolated clusters, one of which appears to be in the process of decay.

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.

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.