Next Article in Journal
Monitoring of Khorasan (Triticum turgidum ssp. Turanicum) and Modern Kabot Spring Wheat (Triticum aestivum) Varieties by UAV and Sensor Technologies under Different Soil Tillage
Previous Article in Journal
Identification and Genomic Characterization of Pathogenic Bacillus altitudinis from Common Pear Trees in Morocco
Previous Article in Special Issue
Effects of Nitrogen and Water Deficiency on Agronomic Properties, Root Characteristics and Expression of Related Genes in Soybean
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Copy Number Quantification for the Soybean Cyst Nematode Resistance Locus rhg1 in the Soybean Varieties of Taiwan

1
Department of Plant Pathology and Microbiology, National Taiwan University, Taipei City 10617, Taiwan
2
Kaohsiung District Agricultural Research and Extension Station, Council of Agriculture, Kaohsiung 908126, Taiwan
3
Taoyuan District Agricultural Research and Extension Station, Council of Agriculture, Taoyuan 327005, Taiwan
*
Author to whom correspondence should be addressed.
Submission received: 24 May 2021 / Revised: 16 June 2021 / Accepted: 24 June 2021 / Published: 30 June 2021
(This article belongs to the Special Issue Soybean Breeding, Agronomic Practices, and Production Systems)

Abstract

:
Disease resistance is one of the most successful strategies in crop protection. For example, the implementation of PI 88788 type resistance, which contains high copy numbers of Resistance to Heterodera glycines 1 (rhg1) loci, into the commercial soybean varieties of the United States has significantly reduced the yield losses caused by soybean cyst nematode (SCN, H. glycines). Vegetable soybean, or edamame, has become a major exporting agricultural product in Taiwan with an annual revenue over $80 million USD since 2017. Several local varieties have been developed to fulfill the market needs such as the traits of flavor and sweetness. However, it remains unclear if the historical breeding programs ever incorporated rhg1 resistance into the varieties of Taiwan. This study applied the TaqMan qPCR method to measure the fluorescent signals specific to the rhg1 locus on the chromosome 18 of soybean, and the ratio of VIC and FAM signals were analyzed to predict the rhg1 copy number in the 21 soybean varieties of Taiwan. The results indicated the copy number and the single nucleotide polymorphisms of the 21 soybean varieties were identical to the susceptible soybean variety ‘Williams 82’. As importation of soybean will be continuously needed to fulfill the market and because SCN is absent in the soybean fields of Taiwan, lacking rhg1 resistance in the local soybean varieties may put the edamame industry at risk and early implementation of SCN resistance in the breeding program, alongside the application of quarantine regulations, will be the key to maintain the SCN-free status and to sustain the edamame industry in Taiwan.

1. Introduction

Soybean cyst nematode (SCN, Heterodera glycines) is the major yield-limiting biotic stress for soybean production worldwide [1]. The earliest report of cyst nematodes in soybean was documented in Japan [2] as Heterodera schachtii Schmidt [3], and was not re-classified as Heterodera glycines Ichinohe until 1952 [4]. Since then, the obligate endoparasitic plant nematode has been globally distributed in the major soybean-producing countries including the US, Argentina, Brazil, and China [5]. SCN infects soybean roots, leading to the stunting and yellowing of the host. Moreover, SCN infection often increases soybean susceptibility to fungal infection, which synergistically reduces the yield drastically. SCN caused annual yield losses of about 100 million bushels during 2006 to 2014 [6,7] and it is responsible for approximately over $1 billion loss per year in the United States.
The life cycle of SCN is composed of the egg stage, four juvenile stages, and the adult stage. The second stage juvenile (J2) SCN is the infection stage. The J2 juvenile uses its stylet and enzymes to penetrate soybean roots, and it migrates inside soybean roots until finding a proper feeding site close to the vascular cylinder. SCN secretes a variety of effectors to induce a specialized multinucleate feeding structure called syncytium [8]. Syncytia supports SCN growth including the sexual dimorphism at the J3 stage and the sexual maturity at the J4 stage. While the female SCN becomes enlarged and exposed on the root surface, male SCN remains reniform and migrates outside of the roots to fertilize the female SCN. After fertilization, the female SCN dies and becomes a cyst with hundreds of eggs inside. The mature and pigmented cyst provides the eggs a shelter from unfavorable environments, and a new life cycle of SCN begins upon encountering plant hosts and suitable environments [9]. The life cycle of SCN can be completed within thirty days if all conditions are optimized [10]. Therefore, SCN may be considered as a polycyclic disease in a single season of soybean production. Reducing the initial population in soils and preventing the spread of SCN to non-infected areas are important aspects of the SCN management. Tillage, removal of root debris in soil, and rotation with non-host crops such as corn (Zea mays L.) may effectively decrease nematode population [11,12]. Chemical management is commonly used for soil fumigation to reduce nematode population in soils [13]. More recently, seed treatments using the succinate dehydrogenase inhibitor (SDHI) fungicide fluopyram have been proposed for effective SCN-control without eliminating beneficial microorganisms in soils [14]. Nonetheless, the most widely applied and effective management strategy relies on the host resistance.
One of the well-known examples is the introduction of rhg1 (resistance to Heterodera glycines) locus from a resistant donor germplasm line, the plant introduction (PI) 88788, into soybean varieties such as ‘Fayette’ and about 90% commercial soybean varieties over decades [15]. The rhg1 locus has been mapped on the chromosome 18 [16], and the rhg1 haplotype of PI 88788 was defined as rhg1-b [15]. The resistance mechanism of rhg1-b has been shown to depend on the copy number variation, and the resistant PI 88788 contains nine copies of a 31.2 Kb cluster repeat with four candidate genes, including an amino acid transporter (Glyma18g02570, Glyma.18G022400), an α-SNAP (soluble NSF attachment protein) protein (Glyma18g02590, Glyma.18G022500), and a PLAC domain protein with unknown function (Glyma18g02600, Glyma.18G022600), and a wound-inducible protein 12 (Glyma18g02610, Glyma.18G022700). Among these genes, the Glyma.18G022400, Glyma.18G022500, and Glyma.18G022700 have been shown to be involved in SCN resistance by RNA silencing, and their expressions were up-regulated when PI 88788 is affected by SCN [15]. On the other hand, the susceptible variety ‘Williams 82’ contains only one copy of rhg1-c [15,17], and therefore the high copy number in PI 88788 may provide SCN resistance via enhanced expression of these three genes. Among these three genes, only the molecular function of α-SNAP protein in SCN resistance has been clearly studied. It has been shown that rhg1-b type of α-SNAP protein (α-SNAPRhg1HC) has amino acid polymorphism different from the ‘Williams 82’ rhg1-c type of α-SNAP protein (α-SNAPRhg1WT). The α-SNAPRhg1HC exhibited cytotoxicity and low affinity to the N-ethylmaleimide-sensitive factor (NSF), which lead to the disruption of SNARE-dependent vesicle trafficking [18]. The over-expression and hyper-accumulation of α-SNAPRhg1HC may disable the sustainability of feeding cells, which leads to syncytia death and SCN resistance. In addition, α-SNAP protein was shown to interact with two syntaxins, which were initially found in the QTL-mapping studies [19,20] and recently confirmed to be important for SCN resistance [21].
Different from the rhg1-b of PI 88788, the rhg1-a haplotype contains 2.3 to 3 copies in the resistant variety ‘Peking’ and the α-SNAP protein derived from ‘Peking’ represents the low copy version (α-SNAPRhg1LC) [15,17]. Nonetheless, both α-SNAPRhg1HC and α-SNAPRhg1LC were expressed in syncytia and both exhibited cytotoxicity in high dosage [18]. More intriguingly, both types displayed a higher affinity to the NSF on chromosome 7 (NSFRAN07, Glyma.07G195900), which contains five amino acid polymorphisms derived from a resistant rhg1-a type variety ‘Forrest’. Their interaction may also explain the segregation disorder observed for the rhg1 on the chromosome 18 and NSF on the chromosome 7 [22]. Other than copy number difference between rhg1-b and rhg1-a type resistant varieties, a retrotransposon (RAC, Rhg1 α-SNAP copia) was found to insert only in α-SNAPRhg1LC, but not α-SNAPRhg1HC or α-SNAPRhg1WT [23]. It is suggested that RAC has promoter activities, and it is possible that RAC provides additional regulation for rhg1-a.
Other than the rhg1 locus, Rhg4 is another reproducible locus from several mapping studies, and the Rhg4 resistance mostly co-occurs with the rhg1-a type. There are three allelic types of Rhg4, which includes the resistant ‘Peking’ type Rhg4-a, the susceptible ‘Williams 82’ type Rhg4-b, and the resistant PI 437654 type Rhg4-c. Copy number variation has also been found for the Rhg4 locus, and ‘Peking’ was found to have 2.3 copies of Rhg4-a while PI 437654 was found to have 4.3 copies of Rhg4-c [17]. Moreover, the sequence polymorphisms in the promoter regions of resistant types rhg1 and Rhg4 loci were found to regulate higher gene expression that together leads to SCN resistance [17]. Map-based cloning has been applied on the Rhg4 locus, which identified a serine hydroxymethyltransferase (SHMT) with amino acid polymorphisms between the resistant ‘Forrest’ and the susceptible ‘Essex’. The amino acid difference in SHMT may interfere its function in carbon folate metabolism, which leads to hypersensitive response-like cell death in syncytia and provides SCN resistance [24].
Studies on the nature of rhg1 and Rhg4 have provide foundation to develop molecular markers to assist breeding programs. While a competitive allele-specific PCR (KASP) assay has been used to detect amino acid polymorphisms in the α-SNAPRhg1HC, α-SNAPRhg1LC, and SHMT [25,26], a TaqMan-based quantitative polymerase chain reaction (qPCR) was developed to assess not only the allelic difference, but also the copy number variation [27]. Using the homeolog of Glyma.18G022600 on chromosome 11 (Glyma.11G234400) and the ‘Williams 82’ genome as a single copy reference, Lee et al. categorized the rhg1 locus into three subtypes, including the F subtype represented by the PI 88788 rhg1-b with nine copies, the P subtype represented by the ‘Peking’ rhg1-a with three copies, and the W subtype represented by the ‘Williams 82’ rhg1-c with one copy [27]. By estimating the qPCR fluorescence ratio, the qPCR assay has been applied to estimate many USDA soybean germplasms for their Rhg1 copy number, such as the PI 548402 (‘Peking’) has 2.8 ± 0.05, the PI 79609 has 7.4 ± 1.57, the PI 88788 was estimated at 9.4 ± 0.78, and the highest in record PI 209332 (‘Fayette’) has 9.8 ± 0.47 copies [28]. As the Taqman qPCR method has shown its merit in the USDA germplasm, the application to soybean from different regions or countries would extend the understanding of the rhg1 copy number variation in the global soybean varieties.
Soybean production in Taiwan locates mostly in the tropical regions and sub-tropical regions, and several endemic soybean varieties have been bred to improve agronomic traits suitable for the climate of Taiwan. The majority of soybean produced in Taiwan is consumed as vegetable soybean, or edamame, which is a classic appetizer in Japanese restaurants. Edamame produced in Taiwan is mostly exported to Japan, United States, and Canada with over $80 million dollars annually since 2017. However, a comprehensive examination on the SCN risk is lacking, and therefore, this study aims to apply the TaqMan qPCR to document the copy number variation and allelic types of rhg1 in the soybean varieties of Taiwan. The genetic insights will contribute to the SCN risk assessment of the potential disease outbreak.

2. Materials and Methods

2.1. Plant Materials and DNA Extraction

A total of 21 soybean (Glycine max (L.) Merr.) varieties were used in this study (Table 1). Soybean seeds were stored in 4 °C with desiccants to maintain the seed vigor. Soybean seeds were germinated on wet and sterilized tissue papers for five days at room temperature, and the taproots were spliced into a 2 mL screw-cap tube. Each biological replicate (screw-cap tube) contained taproot tissues collected from five soybean individuals. Each soybean variety was represented by four independent biological replicates; in other words, the genotype of each soybean variety was averaged from and represented by 20 soybean individuals. The screw-cap tubes were frozen in −20 °C freezer prior to DNA extraction.
The soybean taproots were mechanically homogenized by the CLUBIO Prep-CB24 homogenizer (Bioman, New Taipei City, Taiwan) at 6.0 speed for 15 s, followed by immediately adding 800 μL GPX1 buffer and 5 μL RNase A of the Genomic DNA Mini Kit (Plant) (Geneaid, New Taipei City, Taiwan). After 10 min of incubation at 60 °C, 200 μL GP2 buffer was added into the mixture and the tube was vortexed and incubated for 3 min on ice. A brief centrifugation of 14,000 rpm for 3 min was performed to precipitate plant tissues debris, and the supernatant was subjected to the subsequent DNA extraction procedure following the manufactural instruction. The final DNA was eluted by distilled deionized water, and the DNA concentration was quantified by the Nanodrop One instrument (Thermo Fisher Scientific, Waltham, MA, USA).

2.2. Quantitative Polymerase Chain Reaction (qPCR)

The quantification of rhg1 locus was modified from the methods developed by Lee et al. [27]. The original method applied 12 ng genomic DNA in 5 μL reaction, and the quantification was performed in 384-well plates using the ABI system (Foster City, CA, USA). The modified qPCR condition in this study included 20 ng genomic DNA in 15 μL reaction, and the quantification was performed in 96-well plates using the BioRad CFX Connect qPCR system (Hercules, CA, USA). The sequences of forward primer (5′-GCA GCT GTT GGA ATC ATT CTT TGT T), reverse primer (5′-AAA TTG AAC CTC TTT CTC ATT TGG ATC CT), and the Glyma.18G-specific probe (5′-VIC-BHQ2) as well as the Glyma.11F-specific probe (5′-FAM-BHQ2) followed the design of the original study [27]. The VIC-labeled probe and FAM-labeled probe were synthesized with HPLC cleaning (BIOTOOLS, New Taipei City, Taiwan), and the TOOLS Easy 2X Probe qPCR Mix (BIOTOOLS) was used in this study. Because of the difference in TaqMan qPCR conditions and instruments, the probe ratio of FAM to VIC was adjusted to 10:1 to enable the detection of FAM fluorescence signal in soybean varieties PI 88788 and PI 209332, both contain high VIC fluorescence signal that overwhelms the detection of FAM. In other words, the final concentration of both forward primer and reverse primer was 20 μM, and the concentration of FAM-labeled probe and the VIC-labeled probe was adjusted to 0.75 μM and 0.075 μM in each reaction, respectively. Each biological replicate was technically repeated for three times, and the 21 soybean varieties of Taiwan were measured in separate plates with ‘Williams 82’ (one copy), ‘Peking’ (three copies), and PI 88788 (nine copies) and/or PI 209332 (ten copies) included in each plate. Relative fluorescence units (RFU) estimated by the BioRad CFX Maestro version 2.0 software were collected in R environment 4.0.2 to assess the variations among plates and soybean varieties. Technical measurements of each four biological replicate were averaged, and the four biological replicates were used to calculate the mean and variance representing each soybean variety. The 99% confidence interval was estimated using the R package “ellipse” version 0.4.2, and the RFU scatterplot was graphed in R to visualize the similarity among the soybean varieties.

2.3. Copy Number Analysis

The qPCR signal division of VIC (Glyma.18G-specific probe) and FAM (Glyma.11F-specific probe) of each variety was used in a linear regression analysis and linear model prediction. Repeated measurements on ‘Williams 82’, ‘Peking’, PI 88788, and PI 209332 were performed to ensure the variation among plates were negligible, and a total of 24 measurements on ‘Williams 82’ from six plates, 24 measurements on ‘Peking’, 12 measurements on PI 88788 from three plates, and 20 measurements on PI 209332 from five plates were included to build the linear model. For Gaussian linear regression, the Shapiro–Wilk test was used to evaluate the assumption of normality. Box–Cox transformation was performed to fulfill the normality requirements in modeling linear regression. Considering the rhg1 copy numbers present as discrete values rather than continuous values, the generalized linear models using Poisson regression and negative binomial regression were applied. The Akaike information criterion (AIC) and Bayesian information criterion (BIC) were computed between Gaussian model, Poisson model, and negative binomial model. The model with the less AIC and BIC values was regarded as the best model and used for predicting the rhg1 copy number in the 21 soybean varieties of Taiwan. The prediction values of each soybean variety of Taiwan were pairwise compared to the prediction values of ‘Williams 82’, ‘Peking’, PI 88788, and PI 209332 using the students’ t test, and the significance was claimed at Bonferroni-adjusted p value of 0.01.

2.4. Sequence Polymorphisms of the rhg1 Locus in the Soybean Varieties of Taiwan

A forward primer (5′-CTA GTT AGA GCA TGA ACT GC) and a reverse primer (5′-GTA GTA ACA GGG CTA TCA C) were used to amplify the genomic region harboring two single nucleotide polymorphisms (SNPs) at the 10,978 and 10,995 position to distinguish the ‘Williams 82’ type rhg1-c, the ‘Peking’ type rhg1-a, or the PI 88788 type of rhg1-b locus [26]. The PCR was conducted using the high fidelity Phusion polymerase with HF buffer (New England Biolabs, Ipswich, MA, USA), and the PCR results were confirmed to be a single amplicon using 2% gel electrophoresis before direct submission to Sanger sequencing by the Tri-I Biotech Inc. (New Taipei City, Taiwan). The sequences of ‘Williams 82’, ‘Peking’, PI 88788, PI 209332, and the 21 soybean varieties of Taiwan were aligned using the Clustal Omega software.

3. Results

3.1. Relative Fluorescent Units (RFU) of rhg1 Copy Numbers in the Soybean Varieties of Taiwan

The modified qPCR conditions were used throughout the study to fit the instrument difference. While previous studies have confirmed the rhg1 copy number of PI 88788 and PI 209332 were nine and ten, respectively [15], this study resulted in a 99% confidence interval overlapping for these two varieties. Similar phenomenon was reported in the original publication and the fluorescence difference between PI 88788 and PI 209332 was not significantly different from each other [27]. Although these results indicated the resolution limits of TaqMan qPCR on soybean varieties with high copies, the method was applicable to distinguish the F type (PI 88788, nine copies of rhg1-b), P type (“Peking,” three copies of rhg1-a), and the W type (“Williams 82,” one copy of rhg1-c), for which the 99% confidence intervals separated clearly from each group (Figure 1). The modified TaqMan method was applied to the 21 soybean varieties of Taiwan, and the RFUs of the 21 soybean varieties scattered closely to the W type locus, and none of the 99% confidence interval overlapped with “Peking” or PI 88788 or PI 209332. Accordingly, based on the relative fluorescence of VIC and FAM, the results indicated the possible grouping of the 21 soybean varieties of Taiwan to the susceptible variety “Williams 82.”

3.2. Linear Modelling and Prediction of rhg1 Copy Numbers in the Soybean Varieties of Taiwan

In order to estimate the rhg1 copy number in the 21 soybean varieties of Taiwan, a Gaussian linear model using the ratio of VIC to FAM was built according to their copy number: one for ‘Williams 82’, three for ‘Peking’, nine for PI 88788, and ten for PI 209332. Due to the violation of normality (Shapiro–Wilk test, p = 2.535 × 10−7), the raw ratios of VIC to FAM were Box–Cox transformed and the transformed data were used to build the Gaussian linear model (Shapiro–Wilk test, p = 0.577). As the transformation greatly resolved the problems in normality (Figure 2A), the Box–Cox-transformed linear model was accepted for rhg1 copy number prediction. In addition, considering the values of copy number are discrete integers rather than continuous values with all possible decimals, Poisson and negative binomial regression were modelled to better meet the statistical assumptions. The Gaussian, Poisson, and negative binomial models were selected based on the AIC and BIC statistics, and the results indicated that Box–Cox-transformed Gaussian linear model was the best regression based on the lower AIC and BIC value (Table 2).
Using the Box–Cox-transformed Gaussian linear model to predict the rhg1 copy number in the 21 soybean varieties of Taiwan, the results indicated all of these 21 varieties have rhg1 locus approximately to one copy, and none of these varieties showed significant difference to the ‘Williams 82’ (0.675 ± 0.247) in the student’s t test. On the other hand, the prediction on the ‘Peking’ (4.161 ± 0.372), PI 88788 (9.341 ± 1.581), and PI 209332 (8.792 ± 1.287) was all significantly different from the ‘Williams 82’ as well as the 21 soybean varieties of Taiwan. In summary, the linear prediction of rhg1 copy number in the 21 soybean varieties of Taiwan to be one single copy, which displayed no significant difference from the susceptible variety ‘Williams 82’ (Figure 2B).

3.3. Sequence Polymorphisms of the rhg1 Locus in the Soybean Varieties of Taiwan

Although a single copy of rhg1 indicates the likelihood of SCN susceptibility, sequence polymorphism could still present in the single copy to differentiate the wild type rhg1-c (‘Williams 82’), rhg1-a (‘Peking’), and rhg1-b type (PI 88788). The possibility of a single rhg1-a or rhg1-b, rather than the wild type rhg1-c, presents in the soybean varieties of Taiwan should not be ruled out. In order to assess the sequence polymorphism, the genomic region flanking two single nucleotide polymorphisms (SNPs) of rhg1 locus, the SNP-p of ‘Peking’ representing the rhg1-a, and the SNP-f of PI 88788 representing the rhg1-b (Lee et al., 2015), were amplified and sequenced. Unfortunately, the single copy rhg1 locus in all 21 soybean varieties of Taiwan were found to be identical to the susceptible variety ‘Williams 82’ and therefore the genotypes were determined as rhg1-c (Figure S1).

4. Discussion

This study applied the Taqman qPCR method to quantify 21 soybean varieties of Taiwan, and the results revealed that these 21 varieties of Taiwan have only one copy of rhg1 locus. Although SCN resistance can also be mediated by low copy numbers of rhg1, the allelic types of rhg1 appears to be critical in incorporating with the Rhg4 resistance. For example, in the resistant soybean varieties with the rhg1-a locus such as ‘Peking’, PI 567516C, and PI 437564, the copy numbers were estimated to be 1.9 to 3.5 incorporating with 1 to 4.3 copies of Rhg4-a, Rhg4-b, or Rhg4-c [17]. On the other hand, for soybean varieties with the rhg1-b locus such as PI 88788 and PI 548316, it was suggested that a minimal of 5.6 copies of rhg1-b would be needed [17]. There has been no soybean variety containing a single copy of rhg1-c found to be SCN resistant; and therefore, the results of this study indicated that the 21 soybean varieties of Taiwan are likely susceptible to SCN.
Due to the absence of SCN in Taiwan, genotyping becomes the practical method to evaluate SCN resistance for the soybean varieties developed in Taiwan. While there were some publications claimed the presence of SCN in Taiwan [29,30,31,32], these literatures seemed to misunderstand an old literature [33] probably due to the language inconvenience. In the Hung’s publication (Supplementary File 1), a table on the page 95 and the corresponding notes on the page 96 summarized several nematodes found in different fields of Pingtung County in Taiwan, which were based on the morphological characteristics. These nematodes included Aphelenchus sp., Criconemoides rusticum, Helicotylenchus sp., Meloidogyne sp. Pratylenchus sp., Rotylenchus sp., Trichodorus sp., and Xiphinema sp. The key words of SCN and Heterodera glycines showed up in the page 27 for discussion, where it stated that “Heterodera glycines Ichinohe is absent in Pingtung, and the only sample collected in the Hsinchu Agricultural Research and Extension Station by Mr./Ms. Chang-Hsuan Hung of the Taiwan Agricultural Research Institute was a specimen” (translated from Mandarin). Accordingly, there appears to have no living record of SCN ever being confirmed in Taiwan. More recently, an article entitled “The occurrence and management of plant parasitic nematodes in recent years (translated from Mandarin)” stated the absence of SCN in soybean fields even after several intents of field scouting [34]. Consistently, a thorough nematode survey conducted during 2015–2019 in Taiwan did not discover SCN [35], and field scouting for soybean diseases in Hualien County, Taoyuan County, Kaohsiung City, and Pingtung County in 2019 and 2020 found no SCN.
Because nematode are microscopic worms, the species identification traditionally is based on morphology and body measurements only. For the identification of cyst nematode, the key characters are the vulval cone of the cyst, the length of stylet, tail and hyaline tail terminus of the second stage juvenile [36]. SCN belongs to the schachtii group of Heterodera, and is morphologically very similar to H. schachtii, thus accurate identification requires considerable experience. Therefore, some regions are left with unconfirmed status (CABI, May 24 2021; https://www.cabi.org/cpc) for SCN distribution. For example, SCN was reported in Egypt [37] but there is no further study to validate this early documentation. The record of Egypt in the Invasive Species Compendium of CABI is currently “Absent, Unconfirmed presence record(s), “which is the same situation for Taiwan. The absence of SCN in Taiwan together with the finding that there is no rhg1 resistance in the soybean varieties of Taiwan highlights the risk of SCN invasion and destruction that should not be negligible.

5. Conclusions

In the past five years, importation of soybean from SCN-presenting countries including the US, Brazil, and Canada dominated the market needs in Taiwan (Supplementary Table S1). With the continuing demand, nematode quarantine becomes an important task to minimize the chance of SCN invasion to Taiwan. Besides maintaining SCN on the quarantine pest list, further active and systemic approaches from pre-boarder to boarder, such as requiring appropriate sanitation plan on imported soybean seeds, developing sensitive and specific molecular detection assays for border inspection, as well as establishing effective quarantine measures for emergent nematode mitigation, will reduce the risk of SCN invasion and population establishment in Taiwan. As this study confirmed the absence of rhg1 resistance locus in the 21 soybean varieties of Taiwan, future breeding programs should focus on the introduction of high copy numbers of rhg1-b, rhg1-a with Rhg4. With the successful integration management in quarantine and resistance breeding, production of edamame/soybean in Taiwan could be sustainable and free from SCN devastation.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/agronomy11071346/s1, Table S1. Soybean seeds importation from countries with SCN infestation. Figure S1. Nucleotide polymorphism of the rhg1 locus among 21 soybean varieties of Taiwan. Supplementary File 1. The Publication of Hung, Y., 1958. The red block indicates discussion section related to soybean cyst nematode.

Author Contributions

H.-X.C. conceived the project and established the analysis system. C.-C.H. executed DNA extraction and qPCR measurement. C.-C.H. and H.-X.C. performed the data analyses. K.-L.C. and C.-H.L. helped with the field scouting, soybean materials, and pedigree information. J.-i.Y. and H.-X.C. wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This project is supported by the Ministry of Education, Taiwan (the Yushan Young Scholar Program NTU-109V0601) to Hao-Xun Chang.

Acknowledgments

We thank Chao-Hui Wu from the Tainan District Agricultural Research and Extension Station (Council of Agriculture, Tainan, Taiwan), Chang-Sheng Wang from the Department of Agronomy (National Chung Hsing University, Taichung City, Taiwan), and Da-Gin Lin from the Biotechnology Division and Plant Germplasm Division (Taiwan Agricultural Research Institute, Taichung City, Taiwan) for supporting with the soybean varieties and germplasms.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Niblack, T.L.; Lambert, K.N.; Tylka, G.L. A model plant pathogen from the kingdom Animalia: Heterodera glycines, the soybean cyst nematode. Annu. Rev. Phytopathol. 2006, 44, 283–303. [Google Scholar] [CrossRef]
  2. Hori, S. Phytopathological notes: Sick soil of soybean caused by a nematode. J. Plant Prot. 1915, 2, 927–930. (In Japanese) [Google Scholar]
  3. Ishikawa, T. Principal insect pests and diseases of plants occurring during the fourth year of Taisho. II. Occurrence of “Moon night” and wilt disease of soybean. J. Plant Prot. 1916, 3, 197. (In Japanese) [Google Scholar]
  4. Ichinohe, M. On the soybean nematode, Heterodera glycines n. sp., from Japan. Mag. Appl. Zool. 1952, 17, 1–4. [Google Scholar]
  5. Mitchum, M.G. Soybean resistance to soybean cyst nematode Heterodera glycines: An update. Phytopathology 2016, 106, 1444–1450. [Google Scholar] [CrossRef] [Green Version]
  6. Allen, T.W.; Bradley, C.A.; Sisson, A.J.; Byamukama, E.; Chilvers, M.I.; Coker, C.M.; Collins, A.A.; Damicone, J.P.; Dorrance, A.E.; Dufault, N.S.; et al. Soybean Yield Loss Estimates Due to Diseases in the United States and Ontario, Canada, from 2010 to 2014. Plant Health Prog. 2017, 18, 19–27. [Google Scholar] [CrossRef] [Green Version]
  7. Koenning, S.R.; Wrather, J.A. Suppression of soybean yield potential in the continental United States by plant diseases from 2006 to 2009. Plant Health Prog. 2010, 11, 5. [Google Scholar] [CrossRef] [Green Version]
  8. Mitchum, M.G.; Hussey, R.S.; Baum, T.J.; Wang, X.; Elling, A.A.; Wubben, M.; Davis, E.L. Nematode effector proteins: An emerging paradigm of parasitism. New Phytol. 2013, 199, 879–894. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Shaibu, A.S.; Li, B.; Zhang, S.; Sun, J. Soybean cyst nematode-resistance: Gene identification and breeding strategies. Crop J. 2020, 8, 892–904. [Google Scholar] [CrossRef]
  10. Anand, S.C.; Matson, K.W.; Sharma, S.B. Effect of Soil Temperature and pH on Resistance of Soybean to Heterodera glycines. J. Nematol. 1995, 27, 478–482. [Google Scholar]
  11. Chen, S.Y.; Porter, P.M.; Reese, C.D.; Stienstra, W.C. Crop sequence effects on soybean cyst nematode and soybean and corn yields. Crop Sci. 2001, 41, 1843–1849. [Google Scholar] [CrossRef]
  12. Miller, D.R.; Chen, S.Y.; Porter, P.M.; Johnson, G.A.; Wyse, D.L.; Stetina, S.R.; Klossner, L.D.; Nelson, G.A. Rotation crop evaluation for management of the soybean cyst nematode in Minnesota. Agron. J. 2006, 98, 569–578. [Google Scholar] [CrossRef] [Green Version]
  13. Hartwig, J.; Sommer, H.; Müller, F. Nematicides. Ullmann’s Encycl. Ind. Chem. 2008, 24, 1–5. [Google Scholar]
  14. Roth, M.G.; Jacobs, J.L.; Napieralski, S.; Byrne, A.M.; Stouffer-Hopkins, A.; Warner, F.; Chilvers, M.I. Fluopyram suppresses population densities of Heterodera glycines in field and greenhouse studies in Michigan. Plant Dis. 2020, 104, 1305–1311. [Google Scholar] [CrossRef]
  15. Cook, D.E.; Lee, T.G.; Guo, X.; Melito, S.; Wang, K.; Bayless, A.M.; Wang, J.; Hughes, T.J.; Willis, D.K.; Clemente, T.E.; et al. Copy number variation of multiple genes at Rhg1 mediates nematode resistance in soybean. Science 2012, 338, 1206–1209. [Google Scholar] [CrossRef] [Green Version]
  16. Concibido, V.C.; Diers, B.W.; Arelli, P.R. A decade of QTL mapping for cyst nematode resistance in soybean. Crop Sci. 2004, 44, 1121–1131. [Google Scholar] [CrossRef] [Green Version]
  17. Patil, G.B.; Lakhssassi, N.; Wan, J.; Song, L.; Zhou, Z.; Klepadlo, M.; Voung, T.D.; Stec, A.O.; Kahil, S.S.; Colantonio, V.; et al. Whole-genome re-sequencing reveals the impact of the interaction of copy number variants of the rhg1 and Rhg4 genes on broad-based resistance to soybean cyst nematode. Plant Biotechnol. J. 2019, 17, 1595–1611. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Bayless, A.M.; Smith, J.M.; Song, J.; McMinn, P.H.; Teillet, A.; August, B.K.; Bent, A.F. Disease resistance through impairment of α-SNAP–NSF interaction and vesicular trafficking by soybean Rhg1. Proc. Natl. Acad. Sci. USA 2016, 113, E7375–E7382. [Google Scholar] [CrossRef] [Green Version]
  19. Chang, W.; Dong, L.; Wang, Z.; Hu, H.; Han, Y.; Teng, W.; Zhang, H.; Guo, M.; Li, W. QTL underlying resistance to two HG types of Heterodera glycines found in soybean cultivar ‘L-10’. BMC Genom. 2011, 12, 1–10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Qiu, B.X.; Arelli, P.R.; Sleper, D.A. RFLP markers associated with soybean cyst nematode resistance and seed composition in a “Peking” × “Essex” population. Theor. Appl. Genet. 1999, 98, 356–364. [Google Scholar] [CrossRef]
  21. Dong, J.; Zielinski, R.E.; Hudson, M.E. t-SNAREs bind the Rhg1 α-SNAP and mediate soybean cyst nematode resistance. Plant J. 2020, 104, 318–331. [Google Scholar] [CrossRef]
  22. Bayless, A.M.; Zapotocny, R.W.; Grunwald, D.J.; Amundson, K.K.; Diers, B.W.; Bent, A.F. An atypical N-ethylmaleimide sensitive factor enables the viability of nematode-resistant Rhg1 soybeans. Proc. Natl. Acad. Sci. USA 2018, 115, E4512–E4521. [Google Scholar] [CrossRef] [Green Version]
  23. Bayless, A.M.; Zapotocny, R.W.; Han, S.; Grunwald, D.J.; Amundson, K.K.; Bent, A.F. The rhg1-a (Rhg1 low-copy) nematode resistance source harbors a copia-family retrotransposon within the Rhg1-encoded α-SNAP gene. Plant Direct 2019, 3, e00164. [Google Scholar] [CrossRef] [Green Version]
  24. Liu, S.; Kandoth, P.K.; Warren, S.D.; Yeckel, G.; Heinz, R.; Alden, J.; Yang, C.; Jamai, A.; El-Mellouki, T.; Juvale, P.S.; et al. A soybean cyst nematode resistance gene points to a new mechanism of plant resistance to pathogens. Nature 2012, 492, 256–260. [Google Scholar] [CrossRef] [PubMed]
  25. Liu, S.; Kandoth, P.K.; Lakhssassi, N.; Kang, J.; Colantonio, V.; Heinz, R.; Yeckel, G.; Zhou, Z.; Bekal, S.; Dapprich, J.; et al. The soybean GmSNAP18 gene underlies two types of resistance to soybean cyst nematode. Nat. Comm. 2017, 8, 14822. [Google Scholar] [CrossRef] [PubMed]
  26. Shi, Z.; Liu, S.; Noe, J.; Arelli, P.; Meksem, K.; Li, Z. SNP identification and marker assay development for high-throughput selection of soybean cyst nematode resistance. BMC Genom. 2015, 16, 314. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Lee, T.G.; Diers, B.W.; Hudson, M.E. An efficient method for measuring copy number variation applied to improvement of nematode resistance in soybean. Plant J. 2016, 88, 143–153. [Google Scholar] [CrossRef] [PubMed]
  28. Hudson, M.E.; Diers, B.W.; Lee, T.G. Copy Number Detection and Methods. U.S. Patent US20160265070A1, 2016. [Google Scholar]
  29. Bolla, R.I. The society of nematologists: Testerday, today, and tomorrow. Comp. Parasitol. 2002, 69, 151–156. [Google Scholar] [CrossRef]
  30. Epps, J.M.; Edwards, D.I.; Good, J.M.; Rebois, R.V. Annotated Bibliography of Nematodes of Soybean; Agricultural Research Service, United States Department of Agriculture: New York, NY, USA, 1973; Volume 29, pp. 1882–1968. [Google Scholar]
  31. Inagaki, H. Soybean nematodes. Tech. Bull. ASPAC Food Fertil. Technol. Cent. Taiwan 1977, 31, 17. [Google Scholar]
  32. Ko, H.-R.; Kang, H.; Park, E.-H.; Kim, E.-H.; Lee, J.-K. Identification of Heterodera glycines (Tylenchida; Heteroderidae) Using qPCR. Plant Pathol. J. 2019, 35, 654–661. [Google Scholar] [CrossRef]
  33. Hung, Y. A preliminary report on the plant-parasitic nematodes of soybean crop of the Pingtung district. Agric. Pest News 1958, 5, 1–5. [Google Scholar]
  34. Tsai, T.-T. Proceedings of symposium on the occurrence of important diseases in Taiwan in recent year and development of disease diagnosis, monitoring and control (In Mandarin). Spec. Publ. TARI 2012, 149, 63–81. [Google Scholar]
  35. Van den Hoogen, J.; Geisen, S.; Routh, D.; Ferris, H.; Traunspurger, W.; Wardle, D.A.; de Goede, R.G.M.; Adams, B.J.; Ahmad, W.; Andriuzzi, W.S.; et al. Soil nematode abundance and functional group composition at a global scale. Nature 2019, 572, 194–198. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. European and Mediterranean Plant Protection Organization. Heterodera glycines. Bull OEPP/EPPO Bull. 2009, 39, 59–61. [CrossRef]
  37. Diab, K.A. Occurrence of Heterodera glycines from the Golden Island, Giza, UAR. Nematologica 1968, 14, 148. [Google Scholar] [CrossRef]
Figure 1. Detection of relative fluorescence unit (RFU) among soybeans with different copy numbers of rhg1 locus. ‘Williams 82’ served as the susceptible soybean variety, which contains only one copy number of rhg1 locus. ‘Peking’, PI 88788, and PI 209332 served as the resistant soybean varieties, which contain three, nine, and ten copy numbers of rhg1 locus. All soybean varieties developed in Taiwan were measured to have RFUs undistinguishable from ‘Williams 82’. Each data point was averaged from three technical replicates in the qPCR measurement. Each variety was represented by four biological replicates, which included five plants in one DNA extraction; in other words, the four data points of each variety were represented by a total of 20 soybean plants. The ‘Williams 82’, ‘Peking’, and PI 209332 were measured in every plate as a control, and therefore, 16 data points were shown in the plot. The ellipse circles exhibit the 95% confidence interval for different soybean varieties.
Figure 1. Detection of relative fluorescence unit (RFU) among soybeans with different copy numbers of rhg1 locus. ‘Williams 82’ served as the susceptible soybean variety, which contains only one copy number of rhg1 locus. ‘Peking’, PI 88788, and PI 209332 served as the resistant soybean varieties, which contain three, nine, and ten copy numbers of rhg1 locus. All soybean varieties developed in Taiwan were measured to have RFUs undistinguishable from ‘Williams 82’. Each data point was averaged from three technical replicates in the qPCR measurement. Each variety was represented by four biological replicates, which included five plants in one DNA extraction; in other words, the four data points of each variety were represented by a total of 20 soybean plants. The ‘Williams 82’, ‘Peking’, and PI 209332 were measured in every plate as a control, and therefore, 16 data points were shown in the plot. The ellipse circles exhibit the 95% confidence interval for different soybean varieties.
Agronomy 11 01346 g001
Figure 2. Prediction of the rhg1 copy number in soybean varieties. (A) The quantile–quantile plot (Q-Q plot) showed the raw data of the relative fluorescent unit (RFU) ratio (VIC/FAM) violated the assumptions of normality. (B) The Q-Q plot showed that Box–Cox transformation evidently improves the normality violation, and the transformed data were subjected to the Gaussian linear modeling. (C) The rhg1 copy number prediction using the Gaussian linear model with Box–Cox transformed data successfully separated three resistant varieties including ‘Peking’, PI 88788, and PI 209332 from the susceptible variety ‘Williams 82’. None of the soybean varieties developed in Taiwan showed significant difference to ‘Williams 82’. The significance of student’s t test was declared based on the Bonferroni-adjusted p value at 0.01.
Figure 2. Prediction of the rhg1 copy number in soybean varieties. (A) The quantile–quantile plot (Q-Q plot) showed the raw data of the relative fluorescent unit (RFU) ratio (VIC/FAM) violated the assumptions of normality. (B) The Q-Q plot showed that Box–Cox transformation evidently improves the normality violation, and the transformed data were subjected to the Gaussian linear modeling. (C) The rhg1 copy number prediction using the Gaussian linear model with Box–Cox transformed data successfully separated three resistant varieties including ‘Peking’, PI 88788, and PI 209332 from the susceptible variety ‘Williams 82’. None of the soybean varieties developed in Taiwan showed significant difference to ‘Williams 82’. The significance of student’s t test was declared based on the Bonferroni-adjusted p value at 0.01.
Agronomy 11 01346 g002
Table 1. Twenty-one soybean varieties of Taiwan tested in this study.
Table 1. Twenty-one soybean varieties of Taiwan tested in this study.
VarietyDeveloper *Year *Pedigree *
Hualien No.1 (HL1)CW Tsai1979E32 × Wakajima
Hualien No.2 (HL2)FH Lin, PH Jan, and CS Chang1988Dalian-Tou × Wakajima
Tainan Select No.1 (TNS1)CH Wu, TJ Lien, and SS Wang2003Superior single plant selection from the Cha-Tou variety Edamame
Tainan No.3 (TN3)TJ Lien, TJ Yiu, CH Wu, CS Wu, and YC Wang1998HB (♂) ×
Local variety of Pingtung black soybean (♀)
Tainan No.4 (TN4)TJ Lien, CH Wu, SN Huang, TJ Yiu, and YC Wang1999Superior single plant selection from local variety Chin-Pi-Tou
Tainan No.5 (TN5)TJ Lien, TJ Yiu and CH Wu1999Local variety of Chaozhou black soybean (♂) ×
Local variety of Pingtung black soybean (♀)
Tainan No.10 (TN10)CH Wu and HL Huang2014Tainan No.2(♂) ×
95-(1)-4 (Ryokukou × Tainan No.8) (♀)
Tainan No.11 (TN11)CH Wu201792-(1)-6 (KLG10591 x Tanbaguro) (♂) ×
Tainan No. 3 (♀)
Kaohsiung No.5 (KS5)KF Chen and SC Cheng1995Superior single plant selection from Ryokukou
Kaohsiung No. 7 (KS7)KL Chou and SC Cheng2002Kaohsiung Select No.1 (♂) × Tanbaguro (♀)
Kaohsiung No. 9 (KS9)KL Chou2006Kaohsiung No.2 (♂) × Kaohsiung No.5 (♀)
Kaohsiung Select No.10 (KSS10)-1985SRF400 (♂) × SS (♀)
Kaohsiung No.11 (KS11)KL Chou2011KVS1329 (♂) × HC (♀)
Kaohsiung No.12 (KS12)KL Chou2012Kaohsiung No.6 (♂) × Sayamusume (♀)
Hengchun Black Soybean (HBS)--Germplasm
Hei-Wu-Yeh (HWY)--Germplasm
Shih Shih (SS)-1957Germplasm
Hsiang Chi (HC)--Germplasm
Japan Cha-Tou (JCT)--Germplasm
Chi-Ren Black Soybean (CRBS)--Germplasm
Chi-Ren-Hei-Long (CRHL)--Germplasm
* The “-“ in the table indicates untraceable developer, year, or pedigree.
Table 2. Model selection using the Akaike information criterion (AIC) and Bayesian information criterion (BIC).
Table 2. Model selection using the Akaike information criterion (AIC) and Bayesian information criterion (BIC).
Regression ModelAICBIC
Box–Cox transformed Gaussian linear model271.231278.377
Poisson model322.408327.172
Negative binomial model322.784329.930
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Huang, C.-C.; Yang, J.-i.; Chou, K.-L.; Lin, C.-H.; Chang, H.-X. Copy Number Quantification for the Soybean Cyst Nematode Resistance Locus rhg1 in the Soybean Varieties of Taiwan. Agronomy 2021, 11, 1346. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy11071346

AMA Style

Huang C-C, Yang J-i, Chou K-L, Lin C-H, Chang H-X. Copy Number Quantification for the Soybean Cyst Nematode Resistance Locus rhg1 in the Soybean Varieties of Taiwan. Agronomy. 2021; 11(7):1346. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy11071346

Chicago/Turabian Style

Huang, Cheng-Chun, Jiue-in Yang, Kuo-Lung Chou, Chen-Hsiang Lin, and Hao-Xun Chang. 2021. "Copy Number Quantification for the Soybean Cyst Nematode Resistance Locus rhg1 in the Soybean Varieties of Taiwan" Agronomy 11, no. 7: 1346. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy11071346

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop