Skip to main content
Have a personal or library account? Click to login
Evaluation of Glycine max genotypes (maturity groups IV and V) for resistance to Meloidogyne enterolobii via phenotypic and GWAS assays Cover

Evaluation of Glycine max genotypes (maturity groups IV and V) for resistance to Meloidogyne enterolobii via phenotypic and GWAS assays

Open Access
|Jul 2026

Full Article

Soybean (Glycine max) is a vital crop in North Carolina (NC) and national agriculture. Soybean is planted on over 1.6 million acres in NC, which generates up to 600 million USD annually (USDA-NASS, 2025). In the US, soybeans were planted over 87.2 million acres in 2021, accounting for 31% of soybean production globally (American Soybean Association, 2022). US soybean production is under threat from an economically devastating root-knot nematode, Meloidogyne enterolobii (guava root-knot nematode).

Meloidogyne enterolobii is a known pathogen of soybean and is an obligate parasite that enters the host plant roots to feed (Philbrick et al., 2020). The nematode uses a needle-like mouthpart named a stylet to probe and release cellulolytic and proteolytic enzymes that allow for the nematode to move toward the vascular cylinder of the host (Philbrick et al., 2020). Once it has reached the host’s vascular system, the nematode becomes sedentary and uses its stylet to create a feeding site composed of “giant cells.” These giant cells are 5–7 plant cells that surround the head of the nematode and are formed through the deregulation of the plant cell’s cycle (Favery et al., 2016). The giant cells are multinucleated, enlarged cells that have undergone nuclear division, but not cytokinesis. The juvenile nematodes undergo a series of molts while at the feeding site, reaching the mature adult stage after the final molt. Once the adult females are mature, they become swollen with eggs, which are then laid on the outside of the root in a gelatinous matrix. In a full lifecycle of 30–35 days, one M. enterolobii can typically produce 500–1,000 eggs via parthenogenesis (Philbrick et al., 2020). The giant cells stimulate the surrounding plant cells to divide rapidly creating large galls on the plant roots, which are detrimental to the overall root architecture and therefore the life of the plant (Favery et al., 2016). Symptoms produced by this nematode can include chlorosis, stunting, and most importantly for growers, a reduction in yield.

Meloidogyne enterolobii is currently found in 16 counties in NC (Bonyak et al., 2024). The first report of a field infestation of M. enterolobii in NC was on soybean and cotton plants in two counties in 2011, which shows the ability of M. enterolobii to parasitize and reproduce on soybean as well as its ability to rapidly spread (Ye et al., 2013). The geographic distribution of this nematode in the central Coastal Plains region of NC is a concern, as the population increases rapidly through quick reproduction and is easily spread through soil attached to machines, tools, even footwear (Philbrick et al., 2020). Meloidogyne enterolobii has a very broad range of host plants, including sweet potato, which is often rotated with soybean in this region (Gorny et al., 2023).

There are a number of management techniques that can mitigate the effects and population growth of Meloidogyne enterolobii in agricultural fields. Types of disease management for root-knot nematodes includes cultural, biocontrol, chemical, and plant resistance (Azlay et al., 2023). Plant host resistance is one of the most effective management strategies for root-knot nematodes while also being environmentally sustainable (Schwarz and Gorny, 2024). Resistance genes within the host plant can initiate a hypersensitive response (cell death) in the cells of feeding sites, decrease reproduction, or inhibit effectors released by the root-knot nematode (Saucet et al., 2016; Sato et al., 2019; Bali and Gleason, 2024). Studying host resistance is critically important for agricultural stability and proper management of M. enterolobii.

Once a number of plants from a resistance assay have been identified, genomic commonalities can be determined through analyses such as a genome-wide association study (GWAS). A GWAS examines plant cultivar genomes for genomic variants that could be potentially associated with phenotypic traits such as resistance to M. enterolobii. Identification of a genetic marker facilitates the transfer of resistance to agronomically valuable progeny (Balding et al., 2019). In 1994, a single gene Rmi1 was found to confer partial resistance in soybean to Meloidogyne incognita, but when tested against isolates of M. enterolobii, the gene was overcome and the resistance was not maintained (Luzzi et al., 1994; Schwarz and Gorny, 2024). Plant resistance is a key tool for regulating M. enterolobii and research that pursues new sources of resistance are vital for the continued management of this pathogen.

A study conducted by Schwarz and Gorny (2024) evaluated Glycine max and Glycine soja genotypes for novel sources of resistance to M. enterolobii. Glycine max is the preferred soybean species for commercial breeding programs, and has been cultivated to promote favorable phenotypic traits for growers, such as higher yield and vigor (Schwarz and Gorny, 2024). Wild soybean, Glycine soja, although not possessing the phenotypic traits required by growers, is an important pool of genetic diversity and may be a source of resistance genes to regulate nematode infections (Schwarz and Gorny, 2024). Of the 72 Glycine soja and 44 Glycine max lines tested by Schwarz and Gorny (2024), only five of the G. soja were found to be resistant and none of the G. max showed resistance to M. enterolobii. While this study was successful in analyzing the host status of these soybean lines, it was not able to find any resistance to M. enterolobii in the commercially viable G. max lines. By shifting the focus toward Glycine max genotypes and increasing the total number of soybean lines evaluated, there is potential to find resistance. The objective of this greenhouse screening trial was to assess diverse soybean genotypes (Glycine max) for resistance to M. enterolobii.

1
Materials and methods
1.1
Planting materials

The evaluation of 198 unique lines of Glycine max for resistance to M. enterolobii was conducted at NC State University Method Road Greenhouse facilities. Soybean seed was obtained from the USDA-ARS Germplasm Resources Information Network (GRIN). While the selection of lines for inclusion in the study was ultimately random, the criteria for selection of lines included those in the maturity groups IV or V (ideal for growing in the South and Mid-Atlantic regions) and coming from diverse geographical locations (to capture potential diverse sources of resistance).

1.2
Experimental design and inoculation

The lines were randomized and divided into four rounds (R1, R2, R3, and R4), where each line was tested twice in an initial screening experiment. Each round included five tomato plants (“Rutgers”), which are susceptible to M. enterolobii, as positive inoculation controls. Also included in each round of the initial screen were five Glycine max lines that were among the least susceptible genotypes identified in a previous study (Schwarz and Gorny 2024) (PI 646156, PI 614732, PI 562611, PI 527702, and PI 629013) to act as susceptible controls. Three seeds of each line were sown into a 1:1 mix of steam sterilized sand and soil in a 10.16 cm diameter plastic pot, then thinned to one plant per pot after emergence. Plants were fertilized once at planting using Osmocote Smart-Release Plant Food (The Scotts Company, Marysville, OH).

The plants were grown for 20 days in a greenhouse with temperatures maintained at 27–28℃. At 20 days after seeding, each plant was inoculated with 3,000 M. enterolobii eggs. The inoculum was produced from a NC isolate (Schwarz et al. 2020) maintained on tomato (“Rutgers”) and bell pepper (“California Wonder”) in the greenhouse. The egg inoculum was collected from the culture plants using the NaOCl extraction method of Hussey and Barker (1973). The roots of culture plants were removed from the pots, rinsed in water to remove excess soil, then massaged by hand in a 10% bleach solution (Clorox, The Clorox Company, Oakland, CA) to release the eggs from the roots. The egg and bleach solution was then poured over a set of stacked sieves (from bottom to top, 25, 75, and 250 μm) allowing for the collection of clean eggs on the bottom sieve and thorough washing of any remaining bleach solution. The eggs were decanted from the bottom sieve into a clean 50 mL Falcon tube along with 35 mL of water. A sucrose solution (70% w/v sugar) was added to the Falcon tube to bring the total volume up to 50 mL. The Falcon tube was centrifuged at 800  g for 5 min to allow the eggs to float in the sucrose solution while remaining soil sank to the bottom. The contents of the Falcon tube were poured over the 25 μm sieve, avoiding the debris pellet, and the eggs were rinsed to remove any remaining sucrose. The Falcon tube was rinsed and used to collect the cleaned eggs off the sieve. The eggs were tallied by taking three 100 μL aliquots from the total sample and counted using an inverted Nikon TMS microscope (Nikon Instruments, Melville, NY). The three counts were averaged and multiplied by 10 to determine the number of eggs per milliliter. The calculations were completed to determine the number of milliliters needed to inoculate the soybean lines with 3,000 eggs each.

Plants were inoculated by creating a small, shallow hole (1 cm × 2.5 cm) in the soil near the base of the plant, pipetting the inoculum solution into the hole, and then closing the hole. The plants were maintained through daily watering, increasing to twice daily watering in the summer, in the greenhouse for 60 days post inoculation. The initial screening experiments were carried out from May 2024 to November 2024.

1.3
Data collection

At 60 days post inoculation, the plants were cut at the soil line to remove foliage. The foliage for each plant was weighed. The root systems were removed from the pot and cleaned thoroughly to remove the soil. The gall severity ratings were determined by examining the clean root architecture and assigning a percentage using the modified rating scale of Bridge and Page (1980). The gall severity rating was used as one metric to quantify resistance; soybean lines with less than or equal to 10% galling were identified as possibly resistant and were included in the confirmation screening test. The roots of all plants were weighed, then underwent the NaOCl extraction method as described above to extract M. enterolobii eggs, resulting in a 50 mL egg solution from each plant. Three, 100 μL aliquots were taken from the egg solution and eggs were counted. The counts for the three aliquots were averaged, then multiplied by 10 to estimate the eggs present in 1 mL. This value was then multiplied by 50 to determine the eggs within the total 50 mL egg solution and estimate the total number of M. enterolobii eggs per root system. The eggs per gram of root for each plant was calculated by taking the total number of eggs and dividing by the root weight. Soybean lines with less than or equal to 500 eggs per gram of root tissue were identified as possibly resistant. The reproductive factor (RF) was used to quantify the ability of M. enterolobii to reproduce successfully on the host plant root system. To calculate this, the final population (total number of eggs extracted at the termination of the experiment) was divided by the initial population (number of eggs used for inoculation, herein 3,000). This ratio was utilized to define the host status of the soybean line, where an RF value <1 indicates host resistance and an RF value >1 indicates host susceptibility.

1.4
Confirmation screening experiment

Once all 198 lines were assessed twice for resistance status using the defined resistance metrics of eggs per gram of root (<500 eggs/g root), gall severity rating (<10%), and RF values (<1), the ten soybean lines exhibiting the greatest degree of resistance to M. enterolobii using these metrics were selected for a confirmation screening experiment. This experiment was performed as described above for the initial screening rounds, except soybean lines were tested in high replicate numbers (ten replicates per line). This round also included soybean lines PI 556852 and PI 548559; these lines had an average gall rating of 0% while still having other metrics within the susceptible range. The “Rutgers” tomato plants were included as positive controls in the confirmation screening. Two soybean lines identified as susceptible from a previous study (Schwarz and Gorny 2024), PI 562611 and PI 527702, and two soybean lines (PI 548330 and PI 548386) tested in the initial screening that had extremely high RF values were included as susceptible controls. The confirmation screening was conducted using the same protocols as above. An additional metric, the reproductive index (RI), was calculated to allow for a standardized value to be compared across trials and locations. The RI was calculated by taking the RF for a soybean line and dividing it by the average RF of the tomato controls (herein, RF = 11.45) and then multiplying by 100.

1.5
Analysis of greenhouse data

The statistical analysis was completed in RStudio (version 4.4.1; R Core Team, 2024). The data from the initial screening were not statistically analyzed and instead were observed for evidence of potentially phenotypic resistant lines for the confirmation screening. The data from the confirmation screening were assessed and found to have a non-normal distribution for the following variables: RF, root galling severity, and eggs per gram of root. Therefore, non-parametric methods were used to assess the effect of soybean genotype on these response variables. The effect of soybean line on RF, root galling severity, and eggs per gram of root were evaluated using the Kruskal–Wallis rank sum test. A significant result was found at the ɑ = 0.05 level for the Kruskal–Wallis test for each of these response variables. To separate the means, a Dunn’s Test was completed using the Holm correction. The Dunn’s Test computed the multiple comparisons between soybean and control lines and an adjusted p-value of 0.05 was used to determine significance between the pairwise comparisons. The RI for each soybean line was transformed (log10(x + 1)) to fit assumptions of normality and an ANOVA test indicated that the soybean line had a significant effect on RI at the ɑ = 0.05 level. To separate mean values, Fisher’s least significant difference (LSD) test was applied to identify which soybean lines were significantly different at the ɑ = 0.05 level.

1.6
GWAS

The USDA Soybean Germplasm collection has previously been genotyped with the SoySNP50K beadchip, and the marker data from this study are publicly available on soybase (Grant et al., 2009; Song et al., 2013). While these public data had high density SNP data for 152 of the soybean lines used in this study, they did not have marker data for 46 lines, for which we had the phenotypic data. To obtain genotypic data for the remaining 46 lines, 25 seeds of each line were germinated on germination paper for 5 days after which sprouted root tips were harvested and frozen in liquid nitrogen. Frozen root tissue was then lyophilized for 2 days in a VirTis Benchtop Pro freeze dryer, and DNA was then extracted from the lyophilized tissue using a Qiagen DNEasy Plant Mini DNA extraction kit (Qiagen Sciences Inc., Germantown, MD) following the manufacturers protocol. Extracted DNA was then genotyped with the BARCSoySNP6K genotyping assay at the Soybean Genomics and Improvement Lab in Beltsville, MD, USA (Song et al., 2014). Lower density marker data from these 46 lines were then imputed to the higher density marker set in the publicly available data using the AlphaPlantImpute2 software (Niehoff et al., 2022).

The rTASSEL software was used to perform GWAS analysis in R version 4.5 (Monier et al., 2022). Principal components were calculated with the pcaMethods package, and the first three principal components were used in the analysis to control for population structure (Stacklies et al., 2007). The presence of population structure was checked using quantile-quantile plots with rTASSEL. Total egg count, eggs per gram of root, gall rating, and reproduction factor phenotypes were analyzed in the GWAS analysis and all were log-transformed prior to analysis to normalize their distributions. Marker quality control was performed with rTASSEL. Markers with a minor allele frequency above 0.05, a maximum heterozygosity of 0.2, and a missing frequency below 0.15 were retained in the analysis. The centered IBS method was used to calculate a kinship matrix, and the mixed linear model method was used to perform the GWAS analysis. A Bonferroni correction was utilized to account for multiple testing.

2
Results
2.1
Initial screening host status

In the initial screening experiments (R1 through R4), eight soybean genotypes were found to be resistant to M. enterolobii (Tables 14). In In R1, 26 the soybean lines were foundto be resistant (RF < 1) including PI 548682, PI 548555, PI 639740, PI 548518, PI 548339, PI 548401, PI 553052, PI 556912, PI 548548, PI 548627, PI 595362, PI 619232, PI 559931, PI 548546, PI 506417, PI 548392, PI 548431, PI 598358, PI 548619, PI 612146, PI 548422, PI 556860, PI 559934, PI 548602, PI 595626, and PI 548466) and three of the susceptible control lines were found to be resistant (RF < 1, PI 646156, PI 562611, and PI 629013) (Table 1). In R2, 17 soybean lines were found to be resistant (RF < 1, PI 593256, PI 548667, PI 355070, PI 548598, PI 548679, PI 556906, PI 604100, PI 548559, PI 631122, PI 512039, PI 556852, PI 606748, PI 548309, PI 598222, PI 548622, PI 548359, and PI 543856) and one susceptible control line was found to be resistant (RF < 1, PI 527702) (Table 2). In R3, 11 soybean lines were determined to be resistant (RF < 1, PI 639740, PI 619232, PI 670461, PI 556546, PI 548392, PI 665036, PI 561218, PI 540555, PI 548532, PI 556871, and PI 548466) and one susceptible control line was found to be resistant (RF < 1, PI 629013) (Table 3). In R4, three soybean lines were found to be resistant (RF < 1, PI 548392, PI 548654, and PI 525454) (Table 4).

Table 1

Response of soybean lines to infection by Meloidogyne enterolobii in an initial greenhouse screening experiment, round 1. aAll soybean lines, including the susceptible controls, were planted once (n = 1) and the positive inoculation controls, tomato, had five replicates (n = 5). Plants were inoculated with 3,000 eggs of M. enterolobii and evaluated at 60 days post inoculation. The table includes the root galling severity as a percentage, eggs per gram of root tissue (Eggs g−1 of root), and RF for each line tested. Lines were scored as susceptible (S, RF > 1) or resistant (R, RF < 1)

Soybean lineRoot gall severity (0–100%)RFEggs g−1 of rootSusceptible or resistant
PI 54868200.2854R
PI 54855510.1116R
PI 63974010.1759R
PI 54851810.1741R
PI 54833910.0616R
PI 548401000R
PI 55305200.44145R
PI 55691200.1755R
PI 54854800.1745R
PI 54862720.22106R
PI 59536210.67131R
PI 61923210.56105R
PI 55993100.4459R
PI 553048103.33714S
PI 508269412.913,073S
PI 670461510.331,512S
PI 55654676.501,204S
PI 54854630.78182R
PI 59738437.94973S
PI 543793116.001,644S
PI 59738854.56607S
PI 632418107.001,250S
PI 55691328.781,017S
PI 54379473.33503S
PI 54841301.56188S
PI 54839023.22486S
PI 50641700.1631R
PI 54839210.1738R
PI 59093232.00216S
PI 66774056.39668S
PI 544354107.061,126S
PI 5568763018.502,382S
PI 5484755016.001,420S
PI 5484586010.501,221S
PI 5486964010.83778S
PI 5489779015.392,308S
PI 5336052516.831,048S
PI 54867833.721,074S
PI 355067354.33508S
PI 6650368014.111,647S
PI 5483422023.002,091S
PI 548633805.28701S
PI 54836424.83558S
PI 5483431029.504,116S
PI 614155805.78610S
PI 5565767511.17753S
PI 59365354.283,774S
PI 54843100.06208R
PI 548430554.22379S
PI 548415505.00593S
PI 66402633.72221S
PI 5106706013.721,486S
PI 540884608.11891S
PI 54867121.06434S
PI 561218759.781,281S
PI 5486455011.561,981S
PI 561219807.39972S
PI 5485636023.721,888S
PI 540555309.611,153S
PI 6073804015.442,032S
PI 59835800.1114R
PI 54861910.7884R
PI 5486137016.172,354S
PI 61214610.83102R
PI 548626205.22933S
PI 51866874.17510S
PI 556846104.61699S
PI 66402753.67595S
PI 548344501.39152S
PI 62088326.17826S
PI 57644013.83485S
PI 54842220.94140R
PI 55686010.44150R
PI 55993400.50167R
PI 54854711.67255S
PI 56119158.941,491S
PI 5483865040.6110,413S
PI 59738264.00510S
PI 54853271.72270S
PI 50826855.44770S
PI 5483501011.332,099S
PI 553051109.281,435S
PI 56020704.94693S
PI 54860210.06325R
PI 59562630.5685R
PI 55687122.28419S
PI 54846610.72141R
PI 59365415.83946S
PI 55937018.671,444S
PI 646156 (Control)00.1121R
PI 614732 (Control)53.94497S
PI 562611 (Control)00.1738R
PI 629013 (Control)00.1793R
Tomato Control120.37221S

aLines from round 1 that either did not germinate or died shortly after inoculation include: PI 556826, PI 561596, PI 518664 (seed no longer available at GRIN), PI 556820, PI 601983, PI 548464, PI 577798, PI 561599, PI 556930, PI 527702 (Control), Tomato 1.

Table 2

Response of soybean lines to infection by Meloidogyne enterolobii in an initial greenhouse screening experiment, round 2. bAll soybean lines, including the susceptible controls, were planted once (n = 1) and the positive inoculation controls, tomato, had five replicates (n = 5). Plants were inoculated with 3,000 eggs of M. enterolobii and evaluated at 60 days post inoculation. The table includes the root galling severity as a percentage, eggs per gram of root tissue (Eggs g−1 of root), and RF for each line tested. Lines were scored as susceptible (S, RF > 1) or resistant (R, RF < 1)

Soybean lineRoot gall severity (0–100%)RFEggs g−1 of rootSusceptible or resistant
PI 59325600.1544R
PI 548991204.11512S
PI 5483462011.222371S
PI 595363103.41451S
PI 6067492.53.86566S
PI 54866730.67107R
PI 35507000.0819R
PI 518673255.031376S
PI 548402202.72287S
PI 633970507.89568S
PI 54859800.0615R
PI 54867910.3982R
PI 556635301.3388S
PI 5486543.53.41629S
PI 5346465010.311844S
PI 548987101.44156S
PI 55690600.1723R
PI 60410010.56157R
PI 54855900.0613R
PI 54841052.11603S
PI 5484286014.671588S
PI 564525709.11882S
PI 54854971.56328S
PI 54843976.06804S
PI 54844116.001104S
PI 550734109.221390S
PI 59576555.72777S
PI 556841335.42714S
PI 63112250.5638R
PI 596540356.00684S
PI 561400253.78281S
PI 56020612.74596S
PI 51203900.1730R
PI 55685200.065R
PI 61111278.381497S
PI 60674800.2845R
PI 54830900.3661R
PI 550731204.00508S
PI 5530492512.281835S
PI 59822200.069R
PI 556711303.78290S
PI 615582209.601051S
PI 548517105.891077S
PI 5486220.330.57105R
PI 61480653.78640S
PI 55304352.9424S
PI 548586505.11518S
PI 548606106.441510S
PI 553038211.891518S
PI 5483305517.892736S
PI 612608212.333190S
PI 5722407.54.19713S
PI 548514153.78389S
PI 54868501.10198S
PI 572239507.89947S
PI 583366252.50199S
PI 5483593.50.64141R
PI 602496202.83659S
PI 54385600.3366R
PI 5489766012.561235S
PI 58698111.67177S
PI 5254548012.001040S
PI 646156 (Control)12.83500S
PI 614732 (Control)53.83383S
PI 562611 (Control)206.78992S
PI 527702 (Control)100.5683R
PI 629013 (Control)151.33165S
Tomato 1242.97688S

bLines from round 2 that either died or did not germinate include: PI 630984, PI 355068, PI 561597, PI 548674, PI 618809, PI 533654, PI 548460, PI 533654, PI 548603, PI 556773, PI 564849, PI 583288, PI 556834, PI 636463, PI 633567, PI 550732, PI 548541, PI 665035, PI 548429, PI 560307, PI 556506, PI 553049, PI 548669, PI 548301, PI 550733, PI 564082, PI 527701, PI 596414, PI 54853.

Table 3

Response of soybean lines to infection by Meloidogyne enterolobii in an initial greenhouse screening experiment, round 3. cAll soybean lines, including the susceptible controls, were planted once (n = 1) and the positive inoculation controls, tomato, had five replicates (n = 5). Plants were inoculated with 3,000 eggs of M. enterolobii and evaluated at 60 days post inoculation. The table includes the root galling severity as a percentage, eggs/g root tissue (Eggs g−1 of root), and RF for each line tested. Lines were scored as susceptible (S, RF > 1) or resistant (R, RF < 1)

Soybean lineRoot gall severity (0–100%)RFEggs g−1 of rootSusceptible or resistant
PI 54868212.50.81132S
PI 54855552.22179S
PI 639740100.5642R
PI 548401342.57387S
PI 548548604.61498S
PI 61923200.1115R
PI 559931151.44175S
PI 553048102.28594S
PI 50826958.941,214S
PI 670461101.0091R
PI 55654650.4433R
PI 548546507.721,114S
PI 543793151.78237S
PI 597388301.72146S
PI 632418200.7242S
PI 556913102.83478S
PI 543794202.94429S
PI 54841354.941,514S
PI 548390608.281,364S
PI 5064173017.171,565S
PI 54839200.89202R
PI 56159632.53.61808S
PI 667740101.11154S
PI 556876259.17986S
PI 548475251.62103S
PI 548458202.33149S
PI 548696103.89525S
PI 533605253.56544S
PI 5486786010.89981S
PI 355067103.17642S
PI 66503611.06110R
PI 548342101.83224S
PI 5486332511.221,433S
PI 548364402.06339S
PI 54834355.56947S
PI 614155152.33236S
PI 556576202.06147S
PI 593653103.28692S
PI 548431159.561,129S
PI 548430602.72391S
PI 548415255.83931S
PI 664026103.44517S
PI 51067012.33761S
PI 540884503.72395S
PI 548671605.28851S
PI 561218100.6150R
PI 54864513.33758S
PI 561219609.611,872S
PI 548563509.28910S
PI 54055500.50183R
PI 607380705.61399S
PI 598358609.612,827S
PI 548619102.83368S
PI 5486134011.831,320S
PI 61214617.54.14443S
PI 55682010.55.94828S
PI 548626503.33488S
PI 51866852.39326S
PI 556846202.44346S
PI 664027204.50431S
PI 6019838.31.20332S
PI 5483446012.44821S
PI 620883604.44430S
PI 576440101.75169S
PI 5484222010.612,166S
PI 556860403.39319S
PI 55993411.61250S
PI 54854752.56782S
PI 56119122.53.44474S
PI 54846415.281,667S
PI 548386457.081,079S
PI 597382205.89849S
PI 548532100.89139R
PI 50826852.17363S
PI 5483503.71.61143S
PI 548602307.78558S
PI 577798253.58378S
PI 5956263.71.02146S
PI 55687100.1769R
PI 561599452.80509S
PI 54846600.5518R
PI 593654152.94517S
PI 55693061.36153S
PI 5593706029.062,602S
PI 646156 (Control)02.33222S
PI 562611 (Control)12.53.03576S
PI 629013 (Control)00.2880R
Tomato903.391,640S

cLines from round 3 that either died or did not germinate include: PI 556826, PI 548518, PI 548339, PI 553052, PI 556912, PI 548627, PI 595362, PI 597384, PI 590932, PI 518664 (seed not available from GRIN), PI 544354, PI 548977, PI 553051, PI 560207, PI 614732 (Control), PI 527702 (Control).

Table 4

Response of soybean lines to infection by Meloidogyne enterolobii in an initial greenhouse screening experiment, round 4. dAll soybean lines, including the susceptible controls, were planted once (n = 1) and the positive inoculation controls, tomato, had 5 replicates (n = 5). Plants were inoculated with 3,000 eggs of M. enterolobii and evaluated at 60 days post inoculation. The table includes the root galling severity as a percentage, eggs/g root tissue (Eggs g−1 of root), and reproductive factor (RF) for each line tested. Lines were scored as susceptible (S, RF > 1) or resistant (R, RF < 1)

Soybean lineRoot gall severity (0–100%)RFEggs g−1 of rootSusceptible or resistant
PI 54851803.11915S
PI 548339204.44725S
PI 5484012043.727,287S
PI 556912105.56,1042S
PI 548627206.171,555S
PI 59536253.33667S
PI 5973843014.612,884S
PI 54839200.89202R
PI 590932158.441,189S
PI 544354153.17270S
PI 556930202.28240S
PI 59325619.171,536S
PI 63098413.315.331,622S
PI 548991154.671,157S
PI 5483462512.673,455S
PI 5953634021.172,318S
PI 35506812.53.11377S
PI 6067492010.001,685S
PI 548667157.78791S
PI 355070155.39624S
PI 518673606.941,532S
PI 56159711.73.39407S
PI 54840206.00687S
PI 633970205.78623S
PI 5486746013.616,282S
PI 61880954.721,889S
PI 54859853.67542S
PI 5486792510.33925S
PI 556635155.28519S
PI 54865410.6180R
PI 5346464012.172,897S
PI 548987209.002,213S
PI 55690609.392,965S
PI 60410013.721,283S
PI 54846037.502,073S
PI 54855902.67559S
PI 53365437.611,177S
PI 5486033015.671,787S
PI 5567735018.894,077S
PI 548410208.721,118S
PI 56484918.611,625S
PI 583288512.611,532S
PI 55683435.72698S
PI 6335674030.947,253S
PI 5484281010.39849S
PI 550732109.721,611S
PI 54854953.61686S
PI 55304435.111,357S
PI 54854152.50305S
PI 54843959.671,306S
PI 665035710.281,799S
PI 548441517.893,292S
PI 5484291012.502,072S
PI 5507341030.674,107S
PI 59576505.11431S
PI 560307209.831,085S
PI 631122103.11316S
PI 59654007.17881S
PI 56140004.44620S
PI 55650611.44175S
PI 56020657.50771S
PI 512039116.171,474S
PI 55685212.17323S
PI 561576102.39456S
PI 6111121010.391,430S
PI 60674812.89392S
PI 548309508.781,125S
PI 55073151.44476S
PI 553049257.063,159S
PI 556697159.671,330S
PI 59822219.061,269S
PI 55671111.06132S
PI 61558207.501,071S
PI 54851753.17519S
PI 614808514.173,512S
PI 54862256.611,349S
PI 61480658.501,244S
PI 548586154.06503S
PI 556851105.94529S
PI 5486692718.243,017S
PI 54830133.55269S
PI 548606103.89637S
PI 553038523.114,561S
PI 5483303022.782,680S
PI 55073356.561,009S
PI 612608153.50477S
PI 57224022.83379S
PI 548514309.11952S
PI 54868551.39105S
PI 56408252.67303S
PI 572239254.22370S
PI 5483591040.224,588S
PI 561578304.44460S
PI 602496106.44770S
PI 54385613.611,354S
PI 548976208.78823S
PI 596414605.89448S
PI 58698112.28304S
PI 52545400.78121R
PI 548538103.00209S
PI 562611 (Control)011.891,690S
PI 527702 (Control)108.33812S
PI 629013 (Control)05.33590S
Tomato304.241,487S

dLines from round 4 that either died or did not germinate include: PI 636463, PI 564525, PI 556841, PI 553043, PI 527701, and PI 583366.

2.2
Selection of soybean lines for confirmation screening

A total of ten suspected resistant lines were included in the confirmation screen. All but two lines were under the threshold of resistance for the following metrics: galling rating, RF value, and eggs per gram of root (PI 548682, PI 639740, PI 619232, PI 559931, PI 548392, PI 595626, PI 548466, and PI 561576). Two other soybean lines had RF values slightly above 1.0 but had indications of resistance to M. enterolobii for their gall ratings and eggs per gram of root values (PI 548555 and PI 559934).

2.3
Root galling severity

The Kruskal–Wallis test found that at least one soybean genotype had a significantly different root gall severity rating (P < 0.05). The Dunn Test indicated that soybean lines PI 548682, PI 639740, PI 559931, PI 548466, PI 548555, PI 556852, and PI 548386 were significantly different than the tomato control (P < 0.0001, Table 5) (Figs 1 and 2).

Table 5

Response of soybean lines to infection by Meloidogyne enterolobii in a confirmation screening greenhouse experiment. eAll soybean lines, including the susceptible controls, and the positive inoculation controls, tomato, had ten replicates (n = 10). Plants were inoculated with 3,000 eggs of M. enterolobii and evaluated at 60 days post inoculation. The table includes the root galling severity as a percentage, eggs/g root tissue (Eggs g−1 of root), and RF for each line tested. Values in the table are the average of ten replicates. Lines were scored as susceptible (S, RF > 1) or resistant (R, RF < 1)

Soybean lineRoot gall severity (0–100%)RFEggs g−1 of rootSusceptible or resistant
PI 54868211 abc0.85 ab108 abS
PI 6397408 abc0.79 ab115 abR
PI 61923236 abd1.81 abcd215 abcdeS
PI 5599317 ac1.18 abc216 abcdS
PI 54839245 bd8.14 cd1,339 cdeS
PI 55993423 abcd1.13 abc170 abdS
PI 59562635 abd3.39 bcd492 bcdeS
PI 5484666 ac1.84 abcd598 abcdeS
PI 56157627 abcd2.16 abcd321 abcdeS
PI 5485554 c0.62 a82 aR
PI 5568525 ac0.78 ab126 abR
PI 54855940 abd2.23 abcd497 abcdeS
PI 548330 (Control)46 bd7.85 cd1,682 ceS
PI 548386 (Control)10 abc4.55 bcd1,862 eS
PI 527702 (Control)55 abcd1.56 abcd186 abcdeS
PI 562611 (Control)31 abd1.86 abcd208 abcdS
Tomato71 d11.45 d3,177 eS

eLines from confirmation screening that either died or did not germinate include: PI 548466, PI 561576, PI 548555, PI 556852, PI 548559, PI 548330, PI 548386, and PI 527702 (Control).

Figure 1

(a) Resistant Glycine max, PI 639740, inoculated with Meloidogyne enterolobii. Small, infrequent galls and nitrogen-fixing nodules are present. (b) Susceptible tomato (“Rutgers”) root system inoculated with Meloidogyne enterolobii. Large galls and nitrogen-fixing nodules are present. (c) Susceptible soybean control, PI 527702, root system inoculated with Meloidogyne enterolobii. Large galls and nitrogen-fixing nodules are present. All images are root systems photographed at 60 days post inoculation.

Figure 2

(a) Resistant Glycine max, PI 548555, inoculated with Meloidogyne enterolobii. Very small galls and nitrogen-fixing nodules are present. (b) Susceptible tomato (“Rutgers”) root system inoculated with Meloidogyne enterolobii. Large galls are present. (c) Resistant Glycine max, PI 556852, inoculated with Meloidogyne enterolobii. Very small galls and nitrogen-fixing nodules are present. All images are root systems photographed at 60 days post inoculation.

2.4
RF

The Kruskal–Wallis test found that at least one soybean (Glycine max) genotype had a significantly different RF (P < 0.05). The Dunn Test indicated that soybean lines PI 548682, PI 639740, PI 559931, PI 559934, PI 548555, and PI 556852 were significantly different than the tomato control (P < 0.0001, Table 5).

2.5
Eggs per gram of root

The Kruskal–Wallis test found at least one soybean genotype to have a significantly different number of eggs per gram of root (P < 0.05). The Dunn test indicated that soybean lines PI 548682, PI 639740, PI 559931, PI 559934, PI 548555, PI 556852, and PI 562611 (control line) were significantly different than the tomato control (P < 0.0001, Table 5).

2.6
RI

An ANOVA test found that at least one soybean genotype had a significantly different RI than the others (P < 0.05). The LSD test found soybean lines PI 548386, PI 595626, PI 548559, PI 561576, PI 619232, PI 548466, PI 559934, PI 559931, PI 548682, PI 639740, PI 556852, PI 548555, PI 527702, and PI 562611 to be significantly different than the tomato control (P = 0.025).

2.7
GWAS

A total of 32,674 single nucleotide polymorphism (SNP) markers were used in the GWAS analysis following quality control filtering. A Bonferroni corrected threshold (α = 0.05) was used to determine if an association signal was significant. Following the application of these criteria, no significant associations were detected in the current data set for any of the phenotypes analyzed.

3
Discussion

In this study, the three resistant soybean lines that were under the resistant threshold values for all three-resistance metrics (RF, eggs per gram of root, and root gall severity) were PI 548555, PI 639740, and PI 556852 (Table 5). The soybean line PI 548555 is within the maturity group IV and was developed in Kansas, United States. The soybean line PI 639740 is also within the maturity group IV and was developed in Illinois, United States. The soybean line PI 556852 met these resistance requirements and was included in the confirmation screening due to low galling percentages, despite an elevated RF value in R4 of the initial screening. The soybean line PI 556852 is within the maturity group IV and was developed in the United States. The two soybean lines with high RF values in the initial screening, PI 548330 and PI 548386, that were included in the confirmation screening as susceptible controls were confirmed to be susceptible to M. enterolobii in the confirmation screening.

The five Glycine max lines that were among the least susceptible genotypes identified in a previous study (Schwarz and Gorny, 2024; PI 646156, PI 614732, PI 562611, PI 527702, and PI 629013) were included as susceptible controls in the initial screening trial. While the majority were found to be susceptible, a few were classified as resistant within the four rounds of the initial screen. The previous study by Schwarz and Gorny (2024) for which these results were compared to, was conducted differently. In this study, the soybean lines were planted into 10.16 cm diameter pots, whereas Schwarz and Gorny (2024) used cone containers that were only 3.8 cm in diameter. These cone containers may have restricted the growth of the soybean roots and allowed soil to dry out more rapidly when compared to the large diameter pots (Poorter et al., 2012). The restricted root growth and drought stress potentially put the plant into distress and may have allowed for M. enterolobii to reproduce more quickly than when compared to this study. Additionally, there was a noted decrease in non-germinating or dead plants as the initial screening progressed through the four rounds. The first three rounds were conducted over the warmest months for summer in NC (May through early September) and the soybean plants are particularly sensitive to heat, as they wilt and can enter distress very quickly. To prevent heat stress, watering was increased to twice a day, but the average air temperature during summer of 2024 was high, leading to multiple soybean plants entering distress. The final round of the initial screening was completed from the very end of August through November when the temperature was cooling off.

The confirmation screening used two of the five susceptible controls that were previously tested, PI 562611 and PI 527702. These controls had resistance metrics that identified them as susceptible to M. enterolobii in the confirmation screening, which is consistent with the results from the study by Schwarz and Gorny (2024). In the confirmation screening, these lines had a higher rate of galling (31 and 55%, respectively) when compared to the Schwarz and Gorny (2024) study (15.5 and 11.5%, respectively). Since these lines were grown in cone containers in the Schwarz and Gorny (2024) study, the restricted root architecture could explain the lack of galling (Poorter et al., 2012). This previous study also had higher RF values for these soybean controls (14.86 and 17.41, respectively) when compared to the present study’s confirmation screening (1.86 and 1.55, respectively). As addressed earlier, the increased stress on the plant could have allowed the nematodes to reproduce more successfully than in the present study (Poorter et al., 2012). Here, these soybean lines were tested in replicates of ten in the confirmation screening and in the Schwarz and Gorny (2024) study, only five replicates were used, which could also impact data analysis.

The GWAS indicated that there were no genomic variations that could be attributed to the phenotypic resistance to M. enterolobii seen for the three resistant lines PI 548555, PI 639740 and PI 556852, but this could be due to a number of reasons. The process of assessing plant cultivars for phenotypic resistance, that could indicate potential genotypic resistance, can be a lengthy experiment due to the large number of plants that must be tested. GWAS often include hundreds more cultivars than were used in this study, for logistical reasons involved with the processing of plants infected with M. enterolobii this would take several years to complete. When determining resistance to M. enterolobii, as completed in this study, it can be quite laborious particularly due to the pathosystem involving a nematode. Each plant must be manually processed in its entirety including the root system to properly quantify resistance, this can add countless hours of physical work. Due to the time constraints of this study, the total number of soybean cultivars was limited to 198 lines. While this is still a large number of lines to examine, it does not necessarily provide the power needed to distinguish between alleles that might be conferring resistance to M. enterolobii. Another potential reason for the lack of results from the GWAS could be explained by the low number of resistant lines found in this study. Out of the 198 lines, only 3 were found to be resistant, which is a small percentage of the genetic material analyzed for variations. Other candidate explanations could lie in the as yet unknown genetic control of this resistance as this specific analysis would be unlikely to uncover causal alleles if the resistance is controlled by many alleles with small effects. There is also no precedent for genetic material of soybean that is resistant to M. enterolobii that the genomes can be compared to or analyzed with as these are the first resistant Glycine max lines to M. enterolobii.

The aim of this screening study was to investigate Glycine max lines for resistance to M. enterolobii, and PI 548555, PI 639740, and PI 556852 were found to have resistance. While these lines did present phenotypic resistance, it is vital that these lines be further tested, especially tested out in fields, to ensure that resistance is consistent. In order to mitigate loss of soybean production due to M. enterolobii in NC and other soybean producing states, the development of resistant varieties is essential. While the GWAS for this specific study found no common genetic variation that could be associated with the resistance documented, there is still a future for these Glycine max lines. Resistant progeny are in development through crosses with these three resistant lines, PI 548555, PI 639740, and PI 556852. Once more genetic materials are available for Glycine max lines that are resistant to M. enterolobii, there is potential for future studies such as a GWAS, to investigate if a single gene or loci is regulating the documented resistance. The challenge of M. enterolobii is still present but the use of future resistant cultivars bred from those found in this screening study especially in tandem with other management methods such as the use of crop rotation with non-host crops and chemical management allow growers to get ahead of this nematode pathogen.

Acknowledgements

This work is part of a master’s thesis for the first author. The authors extend thanks to Jessica Dotray and P. Bennett Jeffreys (listed alphabetically by surname) for excellent technical assistance.

Funding information

Authors state no funding involved.

Author contribution statement

SC: Methodology, investigation, formal analysis, writing original draft preparation, review and editing; JG: Investigation, formal analysis, writing original draft preparation, review and editing; ET: Conceptualization, methodology, review and editing; LL: Conceptualization, review and editing; AMG: Conceptualization, supervision, project administration, resources, funding acquisition, review and editing.

Conflict of interest statement

Authors state no conflict of interest.

DOI: https://doi.org/10.2478/jofnem-2026-0018 | Journal eISSN: 2640-396X | Journal ISSN: 0022-300X
Language: English
Page range: 220 - 233
Submitted on: May 15, 2026
Accepted on: Jun 8, 2026
Published on: Jul 17, 2026
In partnership with: Paradigm Publishing Services
Publication frequency: 1 issue per year

© 2026 Sarah Cates, Jay Gillenwater, Earl Taliercio, LeAnn Lux, Adrienne M. Gorny, published by Society of Nematologists, Inc.
This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License.