Interactions of gene expression, alternative splicing, and DNA methylation in determining nodule identity.

Soybean nodulation is a highly controlled process that involves complex gene regulation at both transcriptional and post-transcriptional levels. In the present study, we profiled gene expression changes, alternative splicing events, and DNA methylation patterns during nodule formation, development, and senescence. The transcriptome data uncovered key transcription patterns of nodule development that included 9669 core genes and 7302 stage-specific genes. Alternative splicing analysis uncovered a total of 2323 genes that undergo alternative splicing events in at least one nodule developmental stage, with activation of exon skipping and repression of intron retention being the most common splicing events in nodules compared to roots. Approximately 40% of the differentially spliced genes were also differentially expressed at the same nodule developmental stage, implying a substantial association between gene expression and alternative splicing. Genome-wide-DNA methylation analysis revealed dynamic changes in nodule methylomes that were specific to each nodule stage, occurred in a sequence-specific manner, and impacted the expression of 1864 genes. An attractive hypothesis raised by our data is that increased DNA methylation may contribute to the efficiency of alternative splicing. Together, our results provide intriguing insights into the associations between gene expression, alternative splicing, and DNA methylation that may shape transcriptome complexity and proteome specificity in developing soybean nodules.

Identification of introduced and stably inherited DNA methylation variants in soybean associated with soybean cyst nematode parasitism.

DNA methylation is a widespread epigenetic mark that contributes to transcriptome reprogramming during plant-pathogen interactions. However, the distinct role of DNA methylation in establishing resistant and susceptible responses remains largely unexplored. Here, we developed and used a pair of near-isogenic lines (NILs) to characterize DNA methylome landscapes of soybean roots during the susceptible and resistant interactions with soybean cyst nematode (SCN; Heterodera glycines). We also compared the methylomes of the NILs and their parents to identify introduced and stably inherited methylation variants. The genomes of the NILs were substantially differentially methylated under uninfected conditions. This difference was associated with differential gene expression that may prime the NIL responses to SCN infection. In response to SCN infection, the susceptible line exhibited reduced global methylation levels in both protein-coding genes and transposable elements, whereas the resistant line showed the opposite response, increased global methylation levels. Heritable and novel nonparental differentially methylated regions overlapping with genes associated with soybean response to SCN infection were identified and validated using transgenic hairy root system. Our analyses indicate that DNA methylation patterns associated with the susceptible and resistant interactions are highly specific and that novel and stably inherited methylation variants are of biological significance.

Pathotype Grouping of Cercospora sojina Isolates on Soybean and Sensitivity to QoI Fungicides.

Frogeye leaf spot (FLS), caused by Cercospora sojina, is a common disease of soybean in the southern and northern United States and causes significant yield loss. The use of the current race scheme for classification for C. sojina does not take into account the range of disease severity reactions within each differential. The objective of this research was to better understand the diversity among C. sojina isolates through the development and use of pathogenicity groups. In this study, 83 isolates acquired from 2006 to 2009 were screened using 12 soybean (Glycine max) differentials. Disease severity on the 12 differentials ranged from 0 to 9, where 0 is immune and 9 is very susceptible. The average severity for each isolate across differentials ranged from 1 to 7. The 83 isolates were grouped into five pathogenicity groups (PG): PG1, PG2, PG3, PG4, and PG5, reflecting the severity grouping. Using the 12 differentials, PG1 isolates were differentiated by the lack of infection on Davis, Peking, Kent, Palmetto, Hood, CNS, Tracy, and Richland. PG2 had a range of infections on a scale of 1 to 2 on all differentials except on Davis; PG3 isolates had severity ranging from 3 to 4 except on Davis. PG4 isolates caused no infection on Davis, a maximum disease severity of 5 on Peking, while the rest of differentials had severities from 5 to 6. PG5 isolates caused no infection on Davis, severity of 7 on CNS, and severity of 8 on Kent, Hood, and Palmetto. The remaining seven differentials had severities of 9. Across the geographical locations, the predominant pathotypes were PG3 and PG4 and represented 84% of the tested isolates. Azoxystrobin fungicide sensitivity tests showed that 88% of the isolates were sensitive and dominated the population, while only 6% had a high level of fungicide resistance, suggesting that FLS resistance to the QoI fungicide group was not yet completely developed and had not spread to other areas at the time when these isolates were acquired. The overall virulence profile of the isolates indicated that there was variation in disease severity, suggesting that selection of resistance for each PG may produce lines with more precisely defined interactions to specific pathotypes of C. sojina. This may improve the screening and selection of useful resistance genes that could be pyramided for resistance to each pathogenicity group.

A novel picornavirus-like genome from transcriptome sequencing of sugar beet cyst nematode represents a new putative genus.

Analysis of transcriptome sequence data from eggs and second-stage juveniles (J2s) of sugar beet cyst nematode (SBCN, Heterodera schachtii) identified the full-length genome of a positive-sense single-stranded RNA virus, provisionally named sugar beet cyst nematode virus 1 (SBCNV1). The SBCNV1 sequence was detected in both eggs and J2s, indicating its possible vertical transmission. The 9503-nucleotide genome sequence contains a single long open reading frame, which was predicted to encode a polyprotein with conserved domains for picornaviral structural proteins proximal to its amino terminus and RNA helicase, cysteine proteinase and RNA-dependent RNA polymerase (RdRp) conserved domains proximal to its carboxyl terminus, hallmarks of viruses belonging to the order Picornavirales. Phylogenetic analysis of the predicted SBCNV1 RdRp amino acid sequence indicated that the SBCNV1 sequence is most closely related to members of the family Secoviridae, which includes genera of nematode-transmitted plant-infecting viruses. SBCNV1 represents the first fully sequenced viral genome from SBCN.

Development and Phenotypic Screening of an Ethyl Methane Sulfonate Mutant Population in Soybean.

Soybean is an important oil-producing crop in the Fabaceae family and there are increasing demands for soybean oil and other soybean products. Genetic improvement of soybean is needed to increase its production. In order to provide genetic diversity and resources for identifying important genes, a new ethyl methane sulfonate (EMS) mutagenized soybean population was generated using the newly released germplasm, JTN-5203 (maturity group V). Treatment of soybean seeds with 60 mM EMS concentration was found to be suitable for inducing mutation. A total of 1,820 M1 individuals were produced from 15,000 treated seeds. The resulting M2 population was planted in the field for phenotyping. After harvest, seed traits including total oil, protein, starch, moisture content, fatty acid and amino acid compositions were measured by NIR. Phenotypic variations observed in this population include changes in leaf morphology, plant architecture, seed compositions, and yield. Of most interest, we identified plants with increased amounts of total protein (50% vs. 41% for control) and plants with higher amounts of total oil (25% vs. 21.2% control). Similarly, we identified plants with increases in oleic acid content and decreases in linoleic acid and linolenic acid. This EMS mutant population will be used for further studies including screening for various traits such as amino acid pathways, allergens, phytic acids, and other important soybean agronomic traits. In addition, these mutant individuals will be evaluated in the next generation to assess the heritability. Beneficial traits from these mutants can be exploited for future soybean breeding programs. This germplasm can also be used for discovering novel mutant alleles and for functional gene expression analysis using reverse genetics tools such as TILLING.

An (E,E)-α-farnesene synthase gene of soybean has a role in defence against nematodes and is involved in synthesizing insect-induced volatiles.

Plant terpene synthase genes (TPSs) have roles in diverse biological processes. Here, we report the functional characterization of one member of the soybean TPS gene family, which was designated GmAFS. Recombinant GmAFS produced in Escherichia coli catalysed the formation of a sesquiterpene (E,E)-α-farnesene. GmAFS is closely related to (E,E)-α-farnesene synthase gene from apple, both phylogenetically and structurally. GmAFS was further investigated for its biological role in defence against nematodes and insects. Soybean cyst nematode (SCN) is the most important pathogen of soybean. The expression of GmAFS in a SCN-resistant soybean was significantly induced by SCN infection compared with the control, whereas its expression in a SCN-susceptible soybean was not changed by SCN infection. Transgenic hairy roots overexpressing GmAFS under the control of the CaMV 35S promoter were generated in an SCN-susceptible soybean line. The transgenic lines showed significantly higher resistance to SCN, which indicates that GmAFS contributes to the resistance of soybean to SCN. In soybean leaves, the expression of GmAFS was found to be induced by Tetranychus urticae (two-spotted spider mites). Exogenous application of methyl jasmonate to soybean plants also induced the expression of GmAFS in leaves. Using headspace collection combined with gas chromatography-mass spectrometry analysis, soybean plants that were infested with T. urticae were shown to emit a mixture of volatiles with (E,E)-α-farnesene as one of the most abundant constituents. In summary, this study showed that GmAFS has defence roles in both below-ground and above-ground organs of soybean against nematodes and insects, respectively.

Transgenic soybean overexpressing GmSAMT1 exhibits resistance to multiple-HG types of soybean cyst nematode Heterodera glycines.

Soybean (Glycine max (L.) Merr.) salicylic acid methyl transferase (GmSAMT1) catalyses the conversion of salicylic acid to methyl salicylate. Prior results showed that when GmSAMT1 was overexpressed in transgenic soybean hairy roots, resistance is conferred against soybean cyst nematode (SCN), Heterodera glycines Ichinohe. In this study, we produced transgenic soybean overexpressing GmSAMT1 and characterized their response to various SCN races. Transgenic plants conferred a significant reduction in the development of SCN HG type 1.2.5.7 (race 2), HG type 0 (race 3) and HG type 2.5.7 (race 5). Among transgenic lines, GmSAMT1 expression in roots was positively associated with SCN resistance. In some transgenic lines, there was a significant decrease in salicylic acid titer relative to control plants. No significant seed yield differences were observed between transgenics and control soybean plants grown in one greenhouse with 22 °C day/night temperature, whereas transgenic soybean had higher yield than controls grown a warmer greenhouse (27 °C day/23 °C night) temperature. In a 1-year field experiment in Knoxville, TN, there was no significant difference in seed yield between the transgenic and nontransgenic soybean under conditions with negligible SCN infection. We hypothesize that GmSAMT1 expression affects salicylic acid biosynthesis, which, in turn, attenuates SCN development, without negative consequences to soybean yield or other morphological traits. Thus, we conclude that GmSAMT1 overexpression confers broad resistance to multiple SCN races, which would be potentially applicable to commercial production.

SNP identification and marker assay development for high-throughput selection of soybean cyst nematode resistance.

BACKGROUND: Soybean cyst nematode (SCN) is the most economically devastating pathogen of soybean. Two resistance loci, Rhg1 and Rhg4 primarily contribute resistance to SCN race 3 in soybean. Peking and PI 88788 are the two major sources of SCN resistance with Peking requiring both Rhg1 and Rhg4 alleles and PI 88788 only the Rhg1 allele. Although simple sequence repeat (SSR) markers have been reported for both loci, they are linked markers and limited to be applied in breeding programs due to accuracy, throughput and cost of detection methods. The objectives of this study were to develop robust functional marker assays for high-throughput selection of SCN resistance and to differentiate the sources of resistance. RESULTS: Based on the genomic DNA sequences of 27 soybean lines with known SCN phenotypes, we have developed Kompetitive Allele Specific PCR (KASP) assays for two Single nucleotide polymorphisms (SNPs) from Glyma08g11490 for the selection of the Rhg4 resistance allele. Moreover, the genomic DNA of Glyma18g02590 at the Rhg1 locus from 11 soybean lines and cDNA of Forrest, Essex, Williams 82 and PI 88788 were fully sequenced. Pairwise sequence alignment revealed seven SNPs/insertion/deletions (InDels), five in the 6th exon and two in the last exon. Using the same 27 soybean lines, we identified one SNP that can be used to select the Rhg1 resistance allele and another SNP that can be employed to differentiate Peking and PI 88788-type resistance. These SNP markers have been validated and a strong correlation was observed between the SNP genotypes and reactions to SCN race 3 using a panel of 153 soybean lines, as well as a bi-parental population, F5-derived recombinant inbred lines (RILs) from G00-3213xLG04-6000. CONCLUSIONS: Three functional SNP markers (two for Rhg1 locus and one for Rhg4 locus) were identified that could provide genotype information for the selection of SCN resistance and differentiate Peking from PI 88788 source for most germplasm lines. The robust KASP SNP marker assays were developed. In most contexts, use of one or two of these markers is sufficient for high-throughput marker-assisted selection of plants that will exhibit SCN resistance.

Computational discovery of soybean promoter cis-regulatory elements for the construction of soybean cyst nematode-inducible synthetic promoters.

Computational methods offer great hope but limited accuracy in the prediction of functional cis-regulatory elements; improvements are needed to enable synthetic promoter design. We applied an ensemble strategy for de novo soybean cyst nematode (SCN)-inducible motif discovery among promoters of 18 co-expressed soybean genes that were selected from six reported microarray studies involving a compatible soybean-SCN interaction. A total of 116 overlapping motif regions (OMRs) were discovered bioinformatically that were identified by at least four out of seven bioinformatic tools. Using synthetic promoters, the inducibility of each OMR or motif itself was evaluated by co-localization of gain of function of an orange fluorescent protein reporter and the presence of SCN in transgenic soybean hairy roots. Among 16 OMRs detected from two experimentally confirmed SCN-inducible promoters, 11 OMRs (i.e. 68.75%) were experimentally confirmed to be SCN-inducible, leading to the discovery of 23 core motifs of 5- to 7-bp length, of which 14 are novel in plants. We found that a combination of the three best tools (i.e. SCOPE, W-AlignACE and Weeder) could detect all 23 core motifs. Thus, this strategy is a high-throughput approach for de novo motif discovery in soybean and offers great potential for novel motif discovery and synthetic promoter engineering for any plant and trait in crop biotechnology.

Tillage, Fungicide, and Cultivar Effects on Frogeye Leaf Spot Severity and Yield in Soybean.

Frogeye leaf spot (FLS) of soybean, caused by Cercospora sojina, has been a problem in the southern United States for many years but has become an increasing problem in the northern United States more recently, causing significant yield losses. This increase in disease severity in the northern United States has been attributed to increased utilization of no-till planting and changes in climate. A field study was conducted at the University of Tennessee, Research and Education Center in Milan, TN from 2007 to 2010 to determine severity in tilled and no-till plots treated with or without fungicide at R3 and R5 growth stages. Three FLS-susceptible cultivars, one each in Maturity Groups III, IV, and V, were treated with pyraclostrobin (Headline) fungicide. Analysis of variance using the area under the disease progress curve (AUDPC) indicated no significant difference (P ≤ 0.05) in disease severity between tilled and no-till plots without fungicide. Fungicide did not significantly reduce disease under no-till, but did under tilled plots. This is the first study showing that no-till plots did not reduce or enhance the severity of FLS when no fungicide was applied. Fungicide application significantly reduced (P ≤ 0.05) disease severity and AUDPC and increased yield in tilled plots. The yield gains in tilled, fungicide-treated plots ranged from 1 to 17%. When fungicide was applied, disease severity was not reduced as significantly in no-till as in treated tilled plots, suggesting that fungicide programs under a no-till system may require further study to minimize the risk of FLS severity.

Overexpression of a soybean salicylic acid methyltransferase gene confers resistance to soybean cyst nematode.

Salicylic acid plays a critical role in activating plant defence responses after pathogen attack. Salicylic acid methyltransferase (SAMT) modulates the level of salicylic acid by converting salicylic acid to methyl salicylate. Here, we report that a SAMT gene from soybean (GmSAMT1) plays a role in soybean defence against soybean cyst nematode (Heterodera glycines Ichinohe, SCN). GmSAMT1 was identified as a candidate SCN defence-related gene in our previous analysis of soybean defence against SCN using GeneChip microarray experiments. The current study started with the isolation of the full-length cDNAs of GmSAMT1 from a SCN-resistant soybean line and from a SCN-susceptible soybean line. The two cDNAs encode proteins of identical sequences. The GmSAMT1 cDNA was expressed in Escherichia coli. Using in vitro enzyme assays, E. coli-expressed GmSAMT1 was confirmed to function as salicylic acid methyltransferase. The apparent Km value of GmSAMT1 for salicylic acid was approximately 46 µM. To determine the role of GmSAMT1 in soybean defence against SCN, transgenic hairy roots overexpressing GmSAMT1 were produced and tested for SCN resistance. Overexpression of GmSAMT1 in SCN-susceptible backgrounds significantly reduced the development of SCN, indicating that overexpression of GmSAMT1 in the transgenic hairy root system could confer resistance to SCN. Overexpression of GmSAMT1 in transgenic hairy roots was also found to affect the expression of selected genes involved in salicylic acid biosynthesis and salicylic acid signal transduction.

The receptor like kinase at Rhg1-a/Rfs2 caused pleiotropic resistance to sudden death syndrome and soybean cyst nematode as a transgene by altering signaling responses.

BACKGROUND: Soybean (Glycine max (L. Merr.)) resistance to any population of Heterodera glycines (I.), or Fusarium virguliforme (Akoi, O'Donnell, Homma & Lattanzi) required a functional allele at Rhg1/Rfs2. H. glycines, the soybean cyst nematode (SCN) was an ancient, endemic, pest of soybean whereas F. virguliforme causal agent of sudden death syndrome (SDS), was a recent, regional, pest. This study examined the role of a receptor like kinase (RLK) GmRLK18-1 (gene model Glyma_18_02680 at 1,071 kbp on chromosome 18 of the genome sequence) within the Rhg1/Rfs2 locus in causing resistance to SCN and SDS. RESULTS: A BAC (B73p06) encompassing the Rhg1/Rfs2 locus was sequenced from a resistant cultivar and compared to the sequences of two susceptible cultivars from which 800 SNPs were found. Sequence alignments inferred that the resistance allele was an introgressed region of about 59 kbp at the center of which the GmRLK18-1 was the most polymorphic gene and encoded protein. Analyses were made of plants that were either heterozygous at, or transgenic (and so hemizygous at a new location) with, the resistance allele of GmRLK18-1. Those plants infested with either H. glycines or F. virguliforme showed that the allele for resistance was dominant. In the absence of Rhg4 the GmRLK18-1 was sufficient to confer nearly complete resistance to both root and leaf symptoms of SDS caused by F. virguliforme and provided partial resistance to three different populations of nematodes (mature female cysts were reduced by 30-50%). In the presence of Rhg4 the plants with the transgene were nearly classed as fully resistant to SCN (females reduced to 11% of the susceptible control) as well as SDS. A reduction in the rate of early seedling root development was also shown to be caused by the resistance allele of the GmRLK18-1. Field trials of transgenic plants showed an increase in foliar susceptibility to insect herbivory. CONCLUSIONS: The inference that soybean has adapted part of an existing pathogen recognition and defense cascade (H.glycines; SCN and insect herbivory) to a new pathogen (F. virguliforme; SDS) has broad implications for crop improvement. Stable resistance to many pathogens might be achieved by manipulation the genes encoding a small number of pathogen recognition proteins.

Gene expression profiling of resistant and susceptible soybean lines infected with soybean cyst nematode.

Soybean cyst nematode (SCN) is the most devastating pathogen of soybean. Information about the molecular basis of soybean-SCN interactions is needed to assist future development of effective management tools against this pathogen. Toward this end, soybean transcript abundance was measured using the Affymetrix Soybean Genome Array in a susceptible and a resistant reaction of soybean to SCN infection. Two genetically related soybean sister lines TN02-226 and TN02-275, which are resistant and susceptible, respectively, to the SCN race 2 infection were utilized in these experiments. Pairwise comparisons followed by false discovery rate analysis indicated that the expression levels of 162 transcripts changed significantly in the resistant line, of which 84 increased while 78 decreased. However, in the susceptible line, 1,694 transcripts changed significantly, of which 674 increased while 1,020 decreased. Comparative analyses of these transcripts indicated that a total of 51 transcripts were in common between resistance and susceptible responses. In this set, 42 transcripts increased in the resistant line, but decreased in the susceptible line. Quantitative real-time reverse-transcription polymerase chain reaction confirmed the results of microarray analysis. Of the transcripts to which a function could be assigned, genes were associated with metabolism, cell wall modification, signal transduction, transcription, and defense. Microarray analyses examining two genetically related soybean lines against the same SCN population provided additional insights into the specific changes in gene expression of a susceptible and a resistant reaction beneficial for identification of genes involved in defense.

Iso-lines and inbred-lines confirmed loci that underlie resistance from cultivar 'Hartwig' to three soybean cyst nematode populations.

Soybean [Glycine max (L.) Merr.] cultivars varied in their resistance to different populations of the soybean cyst nematode (SCN), Heterodera glycines, called HG Types. The rhg1 locus on linkage group G was necessary for resistance to all HG types. However, the loci for resistance to H. glycines HG Type 1.3- (race 14) and HG Type 1.2.5- (race 2) of the soybean cyst nematode have varied in their reported locations. The aims were to compare the inheritance of resistance to three nematode HG Types in a population segregating for resistance to SCN and to identify the underlying quantitative trait loci (QTL). 'Hartwig', a soybean cultivar resistant to most SCN HG Types, was crossed with the susceptible cultivar 'Flyer'. A total of 92 F5-derived recombinant inbred lines (RILs; or inbred lines) and 144 molecular markers were used for map development. The rhg1 associated QTL found in earlier studies were confirmed and shown to underlie resistance to all three HG Types in RILs (Satt309; HG Type 0, P = 0.0001 R (2) = 22%; Satt275; HG Type 1.3, P = 0.001, R (2) = 14%) and near isogeneic lines (NILs; or iso-lines; Satt309; HG Type 1.2.5-, P = 0.001 R (2) = 24%). A new QTL underlying resistance to HG Type 1.2.5- was detected on LG D2 (Satt574; P = 0.001, R (2) = 11%) among 14 RILs resistant to the other HG types. The locus was confirmed in a small NIL population consisting of 60 plants of ten genotypes (P = 0.04). This QTL (cqSCN-005) is located in an interval previously associated with resistance to both SDS leaf scorch from 'Pyramid' and 'Ripley' (cqSDS-001) and SCN HG Type 1.3- from Hartwig and Pyramid. The QTL detected will allow marker assisted selection for multigenic resistance to complex nematode populations in combination with sudden death syndrome resistance (SDS) and other agronomic traits.

Molecular marker diversity of SCN-resistant sources in soybean.

Soybean cyst nematode (SCN) (Heterodera glycines Ichinohe; HG) is one of the most destructive pests of soybean (Glycine max (L.) Merr.) in the United States. Over 100 SCN-resistant accessions within the USDA Soybean Germplasm Collection have been identified, but little is known about the genetic diversity of this SCN-resistant germplasm. The objective of this research was to evaluate the genetic variation and determine the genetic relationships among SCN-resistant accessions. One hundred twenty-two genotypes were evaluated by 85 simple sequence repeat (SSR) markers from 20 linkage groups. Non-hierarchical (VARCLUS) and hierarchical (Ward's) clustering were combined with multidimensional scaling (MDS) to determine relationships among tested lines. The 85 SSR markers produced 566 allelic fragments with a mean polymorphic information content (PIC) value of 0.35. The 122 lines were grouped into 7 clusters by 2 different clustering methods and the MDS results consistently corresponded to the assigned clusters. Assigned clusters were dominated by genotypes that possess one or more unique SCN resistance genes and were associated with geographical origins. The results of analysis of molecular variance (AMOVA) showed that the variation differences among clusters and individual lines were significant, but the differences among individuals within clusters were not significant.