Free Standard AU & NZ Shipping For All Book Orders Over $80!
Register      Login
Crop and Pasture Science Crop and Pasture Science Society
Plant sciences, sustainable farming systems and food quality
RESEARCH ARTICLE

RNA-Seq analysis of resistant and susceptible soybean genotypes in response to stress from soybean cyst nematode (Heterodera glycines) HG type 1.2.3.5.7

Haipeng Jiang A , Fanshan Bu A , Shixin Yan A , Yi Li A , Tong Wu A , RuiYao Bai A , Aitong Xu A , Xue Zhao https://orcid.org/0000-0003-3362-1471 A and Yingpeng Han https://orcid.org/0000-0002-9829-6588 A *
+ Author Affiliations
- Author Affiliations

A Key Laboratory of Soybean Biology in Chinese Ministry of Education (Northeastern Key Laboratory of Soybean Biology and Genetics and Breeding in Chinese Ministry of Agriculture), Northeast Agricultural University, Harbin 150030, China.

* Correspondence to: hyp234286@aliyun.com

Handling Editor: Marta Santalla

Crop & Pasture Science 73(3) 238-248 https://doi.org/10.1071/CP21447
Submitted: 25 June 2021  Accepted: 7 September 2021   Published: 22 December 2021

© 2022 The Author(s) (or their employer(s)). Published by CSIRO Publishing

Abstract

The soybean cyst nematode (SCN, Heterodera glycine Ichinohe) is a major disease affecting soybean (Glycine max (L.) Merr.) production and yield. Breeding of new SCN-resistant cultivars and understanding their resistance mechanisms are valuable for improving SCN resistance of soybean. However, we still know little about resistance mechanisms to SCN. The purpose of our analysis was to understand the different resistance mechanisms of resistant and susceptible lines to SCN by dissecting their transcriptional changes during infection by SCN HG type 1.2.3.5.7. In this study, 119 recombinant inbred lines derived from a cross of cv. Dongnong L-204 (SCN resistant) and cv. Heinong 37 (SCN susceptible) were identified for resistance status to SCN HG type 1.2.3.5.7. Roots of the lines identified as extremely resistant or susceptible (L5, L89) were sequenced by transcriptome, and we obtained 66.61 Gb of data. Compared with the untreated control, there were 8394 and 6899 differentially expressed genes in SCN-infected roots of L5 and L89, respectively. Genes of metabolic pathways, biosynthesis of secondary metabolites, plant hormone signal transduction, and plant–pathogen interaction pathway were significantly expressed in both resistant and susceptible genotypes. The expression of genes of phenylpropanoid biosynthesis, flavonoid biosynthesis, thiamine metabolism, cutin, suberin and wax biosynthesis, and endocytosis pathway was significantly higher in the resistant line than the susceptible line. Transcription factor analysis showed that 88 transcription factors from 18 transcription factor families responded to SCN stress. Nine genes were identified by reverse transcriptase qPCR to be associated with SCN resistance. This study helps us to understand better the mechanism of soybean resistance to SCN.

Keywords: DEGs, HG type 1.2.3.5.7 (race 4), metabolic pathways, qPCR, resistance mechanisms, soybean cyst nematode, transcription factor, transcriptome.


References

Arelli PR, Young LD, Concibido VC (2009) Inheritance of resistance in soybean PI 567516C to LY1 nematode population infecting cv. Hartwig. Euphytica 165, 1–4.
Inheritance of resistance in soybean PI 567516C to LY1 nematode population infecting cv. Hartwig.Crossref | GoogleScholarGoogle Scholar |

Alkharouf NW, Khan R, Matthews B (2004) Analysis of expressed sequence tags from roots of resistant soybean infected by the soybean cyst nematode. Genome 47, 380–388.
Analysis of expressed sequence tags from roots of resistant soybean infected by the soybean cyst nematode.Crossref | GoogleScholarGoogle Scholar |

Alkharouf NW, Klink VP, Chouikha IB, Beard HS, MacDonald MH, Meyer S, Knap HT, Khan R, Matthews BF (2006) Timecourse microarray analyses reveal global changes in gene expression of susceptible Glycine max (soybean) roots during infection by Heterodera glycines (soybean cyst nematode). Planta 224, 838–852.
Timecourse microarray analyses reveal global changes in gene expression of susceptible Glycine max (soybean) roots during infection by Heterodera glycines (soybean cyst nematode).Crossref | GoogleScholarGoogle Scholar | 16575592PubMed |

Carter A, Tenuta A, Rajcan I, Welacky T, Woodrow L, Eskandari M (2018) Identification of quantitative trait loci for seed isoflavone concentration in soybean (Glycine max) against soybean cyst nematode stress. Plant Breeding 137, 721–729.
Identification of quantitative trait loci for seed isoflavone concentration in soybean (Glycine max) against soybean cyst nematode stress.Crossref | GoogleScholarGoogle Scholar |

Chen Y, Yang X, He K, Liu M, Li J, Gao Z, Lin Z, Zhang Y, Wang X, Qiu X (2006) The myb transcription factor superfamily of Arabidopsis: expression analysis and phylogenetic comparison with the rice myb family. Plant Molecular Biology 60, 107–124.
The myb transcription factor superfamily of Arabidopsis: expression analysis and phylogenetic comparison with the rice myb family.Crossref | GoogleScholarGoogle Scholar |

Cheng P, Gedling CR, Patil G, Vuong TD, Shannon JG, Dorrance AE, Nguyen HT (2017) Genetic mapping and haplotype analysis of a locus for quantitative resistance to Fusarium graminearum in soybean accession PI 567516C. Theoretical and Applied Genetics 130, 999–1010.
Genetic mapping and haplotype analysis of a locus for quantitative resistance to Fusarium graminearum in soybean accession PI 567516C.Crossref | GoogleScholarGoogle Scholar | 28275816PubMed |

Colgrove AL, Niblack TL (2008) Correlation of female indices from virulence assays on inbred lines and field populations of Heterodera glycines. Journal of Nematology 40, 39–45.

Cook DE, Lee TG, Guo X, Melito S, Wang K, Bayless AM, Wang JP, Hughes TJ, Willis DK, Clemente TE, Diers BW, Jiang J, Hudson ME, Bent AF (2012) Copy number variation of multiple genes at Rhg1 mediates nematode resistance in soybean. Science 338, 1206–1209.
Copy number variation of multiple genes at Rhg1 mediates nematode resistance in soybean.Crossref | GoogleScholarGoogle Scholar | 23065905PubMed |

Cook DE, Bayless AM, Wang K, Guo X, Song Q, Jiang J, Bent AF (2014) Distinct copy number, coding sequence, and locus methylation patterns underlie Rhg1-mediated soybean resistance to soybean cyst nematode. Plant Physiology 165, 630–647.
Distinct copy number, coding sequence, and locus methylation patterns underlie Rhg1-mediated soybean resistance to soybean cyst nematode.Crossref | GoogleScholarGoogle Scholar | 24733883PubMed |

Coquelle N, Glover JNM (2010) FHA domain pThr binding specificity: it’s all about me. Structure 18, 1549–1550.
FHA domain pThr binding specificity: it’s all about me.Crossref | GoogleScholarGoogle Scholar | 21134632PubMed |

Du L, Ali GS, Simons KA, Hou J, Yang T, Reddy ASN, Poovaiah BW (2009) Ca2+/calmodulin regulates salicylic-acid-mediated plant immunity. Nature 457, 1154–1158.
Ca2+/calmodulin regulates salicylic-acid-mediated plant immunity.Crossref | GoogleScholarGoogle Scholar | 19122675PubMed |

Foster PMD, Cattley RC, Mylchreest E (2000) Effects of di-n-butyl phthalate (DBP) on male reproductive development in the rat: implications for human risk assessment. Food and Chemical Toxicology 38, S97–S99.
Effects of di-n-butyl phthalate (DBP) on male reproductive development in the rat: implications for human risk assessment.Crossref | GoogleScholarGoogle Scholar |

Galon Y, Nave R, Boyce JM, Nachmias D, Knight MR, Fromm H (2008) Calmodulin-binding transcription activator (CAMTA) 3 mediates biotic defense responses in Arabidopsis. FEBS Letters 582, 943–948.
Calmodulin-binding transcription activator (CAMTA) 3 mediates biotic defense responses in Arabidopsis.Crossref | GoogleScholarGoogle Scholar | 18298954PubMed |

Gujjar RS, Akhtar M, Singh M (2014) Transcription factors in abiotic stress tolerance. Indian Journal of Plant Physiology 19, 306–316.
Transcription factors in abiotic stress tolerance.Crossref | GoogleScholarGoogle Scholar |

Han Y, Tan Y, Hu H, Chang W, Teng W (2017) Quantitative trait loci with additive and epistatic effects underlying resistance to two HG types of soybean cyst nematode. Plant Breeding 136, 720–727.
Quantitative trait loci with additive and epistatic effects underlying resistance to two HG types of soybean cyst nematode.Crossref | GoogleScholarGoogle Scholar |

He B, Shi P, Lv Y, Gao Z, Chen G (2020) Gene coexpression network analysis reveals the role of SRS genes in senescence leaf of maize (Zea mays L.). Journal of Genetics 99, 3
Gene coexpression network analysis reveals the role of SRS genes in senescence leaf of maize (Zea mays L.).Crossref | GoogleScholarGoogle Scholar | 32089522PubMed |

Jakoby M, Weisshaar B, Dröge-Laser W, Vicente-Carbajosa J, Tiedemann J, Kroj T, Parcy F (2002) bZIP transcription factors in Arabidopsis. Trends in Plant Science 7, 106–111.
bZIP transcription factors in Arabidopsis.Crossref | GoogleScholarGoogle Scholar | 11906833PubMed |

Jiang H, Bu F, Tian L, Sun Q, Bao D, Zhao X, Han Y (2020) RNA-Seq-based identification of potential resistance mechanism against the soybean cyst nematode (Heterodera glycines) HG Type 0 in soybean (Glycine max) cv. Dongnong L-204. Crop & Pasture Science 71, 539–551.
RNA-Seq-based identification of potential resistance mechanism against the soybean cyst nematode (Heterodera glycines) HG Type 0 in soybean (Glycine max) cv. Dongnong L-204.Crossref | GoogleScholarGoogle Scholar |

Jiang H, Tian L, Bu F, Sun Q, Zhao X, Han Y (2021) RNA-seq-based identification of potential resistance genes against the soybean cyst nematode (Heterodera glycines) HG Type 1.2.3.5.7 in ‘Dongnong L-10’. Physiological and Molecular Plant Pathology 114, 101627
RNA-seq-based identification of potential resistance genes against the soybean cyst nematode (Heterodera glycines) HG Type 1.2.3.5.7 in ‘Dongnong L-10’.Crossref | GoogleScholarGoogle Scholar |

Jiao Y, Vuong TD, Liu Y, Meinhardt C, Liu Y, Joshi T, Cregan PB, Xu D, Shannon JG, Nguyen HT (2015) Identification and evaluation of quantitative trait loci underlying resistance to multiple HG types of soybean cyst nematode in soybean PI 437655. Theoretical and Applied Genetics 128, 15–23.
Identification and evaluation of quantitative trait loci underlying resistance to multiple HG types of soybean cyst nematode in soybean PI 437655.Crossref | GoogleScholarGoogle Scholar | 25316311PubMed |

Klink VP, Overall CC, Alkharouf NW, MacDonald MH, Matthews BF (2007) Laser capture microdissection (LCM) and comparative microarray expression analysis of syncytial cells isolated from incompatible and compatible soybean (Glycine max) roots infected by the soybean cyst nematode (Heterodera glycines). Planta 226, 1389–1409.
Laser capture microdissection (LCM) and comparative microarray expression analysis of syncytial cells isolated from incompatible and compatible soybean (Glycine max) roots infected by the soybean cyst nematode (Heterodera glycines).Crossref | GoogleScholarGoogle Scholar | 17668236PubMed |

Liu S, Kandoth PK, Warren SD, Yeckel G, Heinz R, Alden J, Yang C, Jamai A, El-Mellouki T, Juvale PS, Hill J, Baum TJ, Cianzio S, Whitham SA, Korkin D, Mitchum MG, Meksem K (2012) A soybean cyst nematode resistance gene points to a new mechanism of plant resistance to pathogens. Nature 492, 256–260.
A soybean cyst nematode resistance gene points to a new mechanism of plant resistance to pathogens.Crossref | GoogleScholarGoogle Scholar | 23235880PubMed |

Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the CP21447_IE1.gif method. Methods 25, 402–408.
Analysis of relative gene expression data using real-time quantitative PCR and the CP21447_IE1.gif method.Crossref | method.&journal=Methods&volume=25&pages=402-408&publication_year=2001&author=KJ%20Livak&hl=en&doi=10.1006/meth.2001.1262" target="_blank" rel="nofollow noopener noreferrer" class="reftools">GoogleScholarGoogle Scholar | 11846609PubMed |

Martinelli F, Uratsu SL, Albrecht U, Reagan RL, Phu ML, Britton M, et al. (2012) Transcriptome profiling of citrus fruit response to Huanglongbing disease. PLoS ONE 7, e38039
Transcriptome profiling of citrus fruit response to Huanglongbing disease.Crossref | GoogleScholarGoogle Scholar | 22675433PubMed |

Matthews BF, Beard H, MacDonald MH, Kabir S, Youssef RM, Hosseini P, Brewer E (2013) Engineered resistance and hypersusceptibility through functional metabolic studies of 100 genes in soybean to its major pathogen, the soybean cyst nematode. Planta 237, 1337–1357.
Engineered resistance and hypersusceptibility through functional metabolic studies of 100 genes in soybean to its major pathogen, the soybean cyst nematode.Crossref | GoogleScholarGoogle Scholar | 23389673PubMed |

Matthews BF, Beard H, Brewer E, Kabir S, MacDonald MH, Youssef RM (2014) Arabidopsis genes, AtNPR1, AtTGA2 and AtPR-5, confer partial resistance to soybean cyst nematode (Heterodera glycines) when overexpressed in transgenic soybean roots. BMC Plant Biology 14, 96
Arabidopsis genes, AtNPR1, AtTGA2 and AtPR-5, confer partial resistance to soybean cyst nematode (Heterodera glycines) when overexpressed in transgenic soybean roots.Crossref | GoogleScholarGoogle Scholar | 24739302PubMed |

Miraeiz E, Chaiprom U, Afsharifar A, Karegar A, Drnevich JM, Hudson ME (2020) Early transcriptional responses to soybean cyst nematode HG Type 0 show genetic differences among resistant and susceptible soybeans. Theoretical and Applied Genetics 133, 87–102.
Early transcriptional responses to soybean cyst nematode HG Type 0 show genetic differences among resistant and susceptible soybeans.Crossref | GoogleScholarGoogle Scholar | 31570969PubMed |

Mitchum MG, Hussey RS, Davis EL, Baum TJ, Punja ZK, Boer SHD (2007) Application of biotechnology to understand pathogenesis in nematode plant pathogens. In ‘Biotechnology and plant disease management. Vol. 3’. (Eds ZK Punja, SH de Boer, H Sanfaçon) p. 58 (CABI: Wallingford, UK)

Pan G, Hanaoka T, Yoshimura M, Zhang S, Wang P, Tsukino H, Inoue K, Nakazawa H, Tsugane S, Takahashi K (2006) Decreased serum free testosterone in workers exposed to high levels of di-n-butyl phthalate (DBP) and di-2-ethylhexyl phthalate (DEHP): a cross-sectional study in China. Environmental Health Perspectives 114, 1643–1648.
Decreased serum free testosterone in workers exposed to high levels of di-n-butyl phthalate (DBP) and di-2-ethylhexyl phthalate (DEHP): a cross-sectional study in China.Crossref | GoogleScholarGoogle Scholar | 17107847PubMed |

Puthoff DP, Ehrenfried ML, Vinyard BT, Tucker ML (2007) GeneChip profiling of transcriptional responses to soybean cyst nematode, Heterodera glycines, colonization of soybean roots. Journal of Experimental Botany 58, 3407–3418.
GeneChip profiling of transcriptional responses to soybean cyst nematode, Heterodera glycines, colonization of soybean roots.Crossref | GoogleScholarGoogle Scholar | 17977850PubMed |

Suzuki M, Kao CY, McCarty DR (1997) The conserved B3 domain of VIVIPAROUS1 has a cooperative DNA binding activity. The Plant Cell 9, 799–807.
The conserved B3 domain of VIVIPAROUS1 has a cooperative DNA binding activity.Crossref | GoogleScholarGoogle Scholar | 9165754PubMed |

Udvardi MK, Kakar K, Wandrey M, Montanari O, Murray J, Andriankaja A, Zhang J-Y, Benedito V, Hofer JMI, Chueng F, Town CD (2007) Legume transcription factors: global regulators of plant development and response to the environment. Plant Physiology 144, 538–549.
Legume transcription factors: global regulators of plant development and response to the environment.Crossref | GoogleScholarGoogle Scholar | 17556517PubMed |

Usovsky M, Ye H, Vuong TD, Patil GB, Wan J, Zhou L, Nguyen HT (2021) Fine-mapping and characterization of qSCN18, a novel QTL controlling soybean cyst nematode resistance in PI 567516C. Theoretical and Applied Genetics 134, 621–631.
Fine-mapping and characterization of qSCN18, a novel QTL controlling soybean cyst nematode resistance in PI 567516C.Crossref | GoogleScholarGoogle Scholar | 33185711PubMed |

Vuong TD, Sleper DA, Shannon JG, Nguyen HT (2010) Novel quantitative trait loci for broad-based resistance to soybean cyst nematode (Heterodera glycines Ichinohe) in soybean PI 567516C. Theoretical and Applied Genetics 121, 1253–1266.
Novel quantitative trait loci for broad-based resistance to soybean cyst nematode (Heterodera glycines Ichinohe) in soybean PI 567516C.Crossref | GoogleScholarGoogle Scholar | 20559815PubMed |

Wan J, Vuong T, Jiao Y, Joshi T, Zhang H, Xu D, Nguyen HT (2015) Whole-genome gene expression profiling revealed genes and pathways potentially involved in regulating interactions of soybean with cyst nematode (Heterodera glycines Ichinohe). BMC Genomics 16, 148
Whole-genome gene expression profiling revealed genes and pathways potentially involved in regulating interactions of soybean with cyst nematode (Heterodera glycines Ichinohe).Crossref | GoogleScholarGoogle Scholar | 25880563PubMed |

Wu X, Blake S, Sleper DA, Shannon JG, Cregan P, Nguyen HT (2009) QTL, additive and epistatic effects for SCN resistance in PI 437654. Theoretical and Applied Genetics 118, 1093–1105.
QTL, additive and epistatic effects for SCN resistance in PI 437654.Crossref | GoogleScholarGoogle Scholar | 19184662PubMed |

Zhen-Hua W, Li-Bo S, Hai-Yan WU, Jing L, Xiu-Xia LI (2009) Distribution and developmental process of Heterodera glycines in soybean root. Scientia Agricultura Sinica 42, 3147–3153.