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.

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.

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.

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.