Abstract
Introduction:
Pine wilt disease (PWD), caused by the pine wood nematode (PWN) Bursaphelenchus xylophilus, poses a significant threat to global pine forests. The sulfate permease (SULP) family is essential for sulfate transport, sulfur assimilation and cellular homeostasis, yet it remains uncharacterized in B. xylophilus. This study aimed to comprehensively identify all members of the SULP family in B. xylophilus and to elucidate their roles in nematode development and stress adaptation.
Methods:
Through genomic data analysis, we identified 10 members of the SULP family in B. xylophilus and conducted a comprehensive characterization of their physicochemical properties, conserved motifs, protein structures, and gene expression profiles across different developmental stages.
Results:
The results revealed Bx-sulps were located on 5 chromosomes of B. xylophilus. Phylogenetic analysis unveiled both conserved and divergent evolutionary patterns of these proteins compared to counterparts in other nematodes. Expression analysis demonstrated upregulation of Bx-sulps during the dauer third-instar larva (D3) stage, suggesting their involvement in stress response and diapause. Moreover, certain Bx-sulps exhibited high expression levels in adult stages, indicating a potential role in reproductive processes.
Discussion:
The study presents the first comprehensive examination of BxSULP family, shed light on its significance in nematode development and stress adaptation. These findings provide the groundwork for further functional investigations and may aid in the development of targeted strategies for managing PWD.
1 Introduction
Pine wilt disease (PWD), caused by the pine wood nematode (PWN) Bursaphelenchus xylophilus, is one of the greatest threats to pine (). First identified in Japan in the 20th century, PWD has become a major concern for global forestry (). In China, PWD has impacted 620 counties in 18 provinces, spreading rapidly to westward and northward (; ). At present, the management strategies for PWD focus on integration monitoring, early detection, and rapid response methods. Integrated pest management (IPM) practices, including cultural techniques (such as proper watering, mulching, and pruning to enhance tree health), biological control, and chemical interventions (insecticides and nematicides), are commonly employed to reduce PWD incidence (). B. xylophilus is primarily spread by vector insects, notably the genus Monochamus. These insects introduce B. xylophilus into healthy pines while laying eggs or feeding on infected ().
Under natural conditions, B. xylophilus alternates between two life-history modules: a propagative phase within living pines and a dispersal phase between hosts (). During propagation, eggs hatched into second instar larva (L2), which molt through third instar larva (L3) and fourth instar larva (L4) before maturing into adults (). Entirely confined to the resin canal system, this phase supported exponential population growth under favorable temperature and moisture. Declining temperature and host humidity triggered L2 to divert into the dispersal pathway, molting into dauer third-instar larva (D3) (). These accumulate lipid reserves, thicken their cuticles, and aggregated around pupal chambers of the M. alternatus (). After quiescence, D3 molted into dauer fourth-instar larva (D4), which were transported by emerging beetles and inoculated into new host trees during maturation feeding, initiating further infections ().
B. xylophilus is capable of proliferating and multiplying within the resin canals of its pine host, a trait closely linked to the adaptive evolution of its gene family (). Investigating the expression patterns and roles of individual gene family members of B. xylophilus can elucidate their contributions to the nematode’s growth and development, offering valuable insights for identifying molecular targets to manage B. xylophilus. Research of gene families in B. xylophilus has advanced considerably in recent decades. Genome analysis has revealed several gene families crucial for the adaptive evolution and pathogenicity of B. xylophilus. The cytochrome P450 (CYP450) family in B. xylophilus has undergone significant amplification, facilitating the detoxification and metabolism of pine defense chemicals by the nematode (). The flavin-containing monooxygenase (FMO) family contributes to detoxification processes, enabling the nematode to thrive in pine resin channels (). The glycoside hydrolase (GH) family aids in the degradation of pine cell walls, thereby enhancing the nematode’s pathogenicity (). Nevertheless, investigations into the Sulfate permease family in B. xylophilus remain relatively limited. Sulfur is essential for synthesizing sulfur-containing amino acids like cysteine and methionine, crucial for protein synthesis in nematodes (). This suggested that Bx-sulps might similarly determine PWN fitness within the host, making them high-value targets for RNAi- or HIGS-based control strategies.
In B. xylophilus, the SULP family may be linked to the evolutionary adaptation of nematodes to pine environments. Through enhanced sulfate uptake and utilization, B. xylophilus may exhibit improved adaptation to pine defense mechanisms. To date, the functional characterization of the SULP family in B. xylophilus remains unexplored. This study represents the first systematic investigation into the SULP family within B. xylophilus. The research aimed to identify all members of the SULP family in B. xylophilus and to elucidate their roles in nematode development and stress adaptation.
2 Materials and methods
2.1 Physicochemical properties analysis of BxSULPs
Firstly, we obtained the whole genome of B. xylophilus (SAMEA7282713) (https://www.ncbi.nlm.nih.gov/datasets/genome/GCA_904067135.1/) data from Wormbase (https://parasite.wormbase.org). To obtain a representative set of SULP orthologues, we first screened all publicly available nematode genomes deposited in NCBI database (https://www.ncbi.nlm.nih.gov/gene) () for the presence of complete, well-annotated SULP families. After applying strict criteria for sequence completeness and annotation quality, only three genomes (Brugia malayi, Caenorhabditis elegans, and Caenorhabditis briggsae) were retained. We used BLASTX to compare the queried protein sequence with the protein database of B. xylophilus to obtain the candidate BxSULPs sequence. BLAST searches were performed with an E-value cut-off of 1 × 10-5, ≥ 30% query coverage and ≥ 25% identity. After that, we queried the Pfam database (https://pfam.xfam.org/) to search for the conserved structural domains possessed by the SULP protein and downloaded its Hidden Markov Model (). The protein sequences with the Hidden Markov Model were then searched from the protein database of B. xylophilus by HMMER3.0 software. HMMER 3.0 was run using the trusted cut-off (TC) defined in Pfam for the Sulfate_transp and STAS domains, with an inclusion threshold of E ≤ 1 × 10–5 and a minimum alignment coverage of 70% of the Pfam model. Finally, the candidate sequences were characterized using NCBI-conserved domains (https://www.ncbi.nlm.nih.gov/Structure/cdd) () and InterPro (https://www.ebi.ac.uk/Interpro) () to eliminate duplicates and acquire the protein sequences of BxSULPs. The physicochemical properties of BxSULPs were analyzed using Expasy 3.0 (https://web.expasy.org/protparam) ().
2.2 Phylogenetic analysis of BxSULPs and analysis of conserved motifs
In order to evaluate the extent of phylogenetic similarity among the 10 protein sequences, we conducted phylogenetic and conserved motif analyses of BxSULPs. Evolutionary relationships of BxSULP were analyzed using MEGA 11.0 software through multiple sequence alignments of 10 BxSULPs and the optimal model for the maximum likelihood tree was determined. Based on the prediction results of the maximum likelihood tree in the best model, the LG model was selected to construct the phylogenetic evolutionary tree of BxSULP. This model had a BIC (Bayesian Information Criterion) value of 11645.256 and an AICc (Corrected Akaike Information Criterion) value of 11362.515. The model was chosen because it provided the lowest BIC value (11645.256), indicating that it is the best fit for the data among the tested models (). A maximum likelihood evolutionary tree was constructed with the Bootstrap test performed 1000 times (). Subsequently, the constructed phylogenetic evolutionary tree was landscaped using iTOL (https://itol.embl.de/tree/) (). Gene information for Bx-sulps was obtained from the gene structure annotation file, and gene structure analysis was conducted using the GSDS 2.0 (http://gsds.gao-lab.org/) (). Additionally, the distribution of conserved motifs in BxSULPs was examined using MEME Suite 5.5.6 (https://meme-suite.org/meme/) () with default parameters and a motif value of 8. The results from MEME were combined with the maximum likelihood tree to create gene structure maps with annotated conserved motifs.
2.3 Chromosomal distribution and protein structure analysis of BxSULPs
The gene density information of B. xylophilus from the annotation file of gene structure was extracted by Gene Density Profile function of TBtools. Then, we integrated the annotation file and the density file to map the distribution of BxSULP on chromosomes by Gene Location Visualize from GTF/GFF function (, ). Subsequent analysis included predicting the secondary structure of BxSULPs using online tools like Prabi (https://doua.prabi.fr/software/cap3) (), NetPhos 3.1 (https://services.healthtech.dtu.dk/services/NetPhos-3.1/) (), and Protter 1.0 (http://wlab.ethz.ch/protter/start/) (). Swiss-Model (https://swissmodel.expasy.org/interactive) () was employed to forecast the protein tertiary structure, and SAVES v6.1 (https://saves.mbi.ucla.edu) was used to create Ramachandran diagrams.
2.4 Isolation and identification of B. xylophilus
B. xylophilus was isolated from symptomatic Pinus thunbergii logs collected in Zhejiang Province, China. Xylem samples were cut into 2–3 cm thick disks and placed on a modified Baermann funnel (25 °C, 24 h) (). Nematodes that migrated into the collection tube were recovered and examined under a light microscope. Individual adult females were carefully picked out and identified to species according to the morphological diagnostic criteria (). Only those specimens that exhibit the typical morphological characteristics of B. xylophilus are retained for subsequent cultures and experiments.
2.5 Culture of B. xylophilus in different developmental stages
Synchronized eggs were obtained from embryos by incubating mixed-stage B. xylophilus in petri dishes at 25 °C in the absence of light for 10 minutes. This process facilitated the adhesion of eggs to the substrate through surface glycoproteins. The top layer containing water and nematodes was then removed to harvest the synchronized eggs. Subsequently, these eggs were allowed to hatch at 25 °C in the dark to yield second instar larvae (L2) in a nutrient-free setting. These synchronized L2 larvae were then transferred onto a Botrytis cinerea lawn on a PDA plate. Nematodes were harvested at 24, 48, and 72-hour intervals to obtain third instar larvae (L3), fourth instar larvae (L4), and adults, respectively. Dauer third-instar larvae (D3) and dauer fourth-instar larvae (D4) were isolated from infected pinewood by Baermann funnel. Male and female nematodes were manually sorted under microscope using worm pickers ().
2.6 Expression profiles of BxSULP family in different developmental stages
Gene expression in distinct developmental stages was analyzed by reverse-transcription quantitative polymerase chain reaction (RT-qPCR) by synchronizing B. xylophilus in different developmental stages. To explore the molecular response mechanisms of Bx-sulps to various treatments, expression levels of 10 Bx-sulps were determined across various developmental stages of B. xylophilus, including eggs. The statistical significance of the data was determined using a one-way analysis of variance (ANOVA) using GraphPad Prism 10.6.0.890 with a student’s t-test with p < 0.05 as the significance threshold.
2.7 Validation of gene expression by RT-qPCR
Differentially expressed genes were chosen for validation via RT-qPCR (). Primer sequences were generated with Primer 6.0 software and subsequently synthesized by Sangon Biotech (Shanghai, China) (Table 1). Amplification was performed on a QuantStudio 1 Plus instrument (Thermo Fisher Scientific, USA) following the manufacturer’s instructions for the One Step PrimeScript™ RT-PCR Kit (RR064A, Beijing), and relative gene expression was calculated using the -ΔΔCt method. Gene-expression data were normalized to the Bx28S rRNA reference gene.
Table 1
| Protein ID of B. xylophilus | Protein name | Primer sequences |
|---|---|---|
| BXYJ5.010268100 | BxSULP1 | F: CCCACTTTCAGGCTCTCCAG R: CCGAAACGACAGATTGGTGC |
| BXYJ5.020105600 | BxSULP2 | F: CTGGTTTCATTAACGGCGGC R: ACACCGTTAAGAGCAGGGTG |
| BXYJ5.020147800 | BxSULP3 | F: AACTTTTCGCTGTCTCCGCT R: AGACGGCGATGTGTTCTGAT |
| BXYJ5.030165800 | BxSULP4 | F: CACGCTTTCGAGCTATGTGC R: GAGACCGCCAATACTCTGGG |
| BXYJ5.030215700 | BxSULP5 | F: GGCCAAAATTCCAGCACCTG R: CCAGCAGATCCGAGCCATAG |
| BXYJ5.050106800 | BxSULP6 | F: CTTGGCTTCATGGGGACACT R: GACTCCACACTCAACGCTCA |
| BXYJ5.050140100 | BxSULP7 | F: AACTGCCCTGATGGTTGGAG R: CGGGGGCTACATCCATTCTC |
| BXYJ5.050219300 | BxSULP8 | F: TTTTGTCGCACCTGACCGTA R: GTCAATTCGTATACACTTTTTCCGT |
| BXYJ5.060009300 | BxSULP9 | F: TGATAGGGAAAGAGGGGCGA R: GCCAGAACGGAGTAGGCAAT |
| BXYJ5.060037700 | BxSULP10 | F: GAAGCCACTGCCAAACCAAG R: TCATCGGGTGGGAGTCTGAT |
| 28S | / | F: AACCGAACACGCGACAATAG R: GTGCGTATTCAGCCTTCTGG |
RT-qPCR primer sequences.
3 Results
3.1 Identification of BxSULPs
Following a BLASTX search in the Wormbase database and subsequent validation of conserved domains using Conserved Domains and InterPro, we identified 10 distinct SULP proteins in B. xylophilus, designated as BxSULP1 to BxSULP10 based according to their ascending order on chromosomes (Table 2). To further elucidate the evolutionary relationships, we conducted a comparative analysis of the protein sequences of SULP found in BxSULP1 to BxSULP10 with those of analogous proteins from other nematodes.
Table 2
| Protein ID of B. xylophilus | Protein name | Sequence ID in NCBI | Protein name in other nematodes | Other nematodes name | Scores of blastP | Similarity e-value | % identity |
|---|---|---|---|---|---|---|---|
| BXYJ5.010268100 | BxSULP1 | CBG_23242 | Protein CBR-SULP-4 | Caenorhabditis briggsae | 335 | 1.8×10-121 | 36.2 |
| BXYJ5.020105600 | BxSULP2 | CELE_F41D9.5 | Protein CE-SULP-3 | Caenorhabditis elegans | 302 | 1.5×10-140 | 58.4 |
| BXYJ5.020147800 | BxSULP3 | CELE_ZK287.2 | Protein CE-SULP-8 | Caenorhabditis elegans | 199 | 8.7×10-62 | 35.0 |
| BXYJ5.030165800 | BxSULP4 | CBG_23242 | Protein CBR-SULP-4 | Caenorhabditis briggsae | 224 | 2.6×10-37 | 28.4 |
| BXYJ5.030215700 | BxSULP5 | CELE_K12G11.2 | Protein CE-SULP-5 | Caenorhabditis elegans | 147 | 2.3×10-37 | 28.3 |
| BXYJ5.050106800 | BxSULP6 | CBG_23241 | Protein CBR-SULP-5 | Caenorhabditis briggsae | 234 | 2.2×10-101 | 47.0 |
| BXYJ5.050140100 | BxSULP7 | CELE_W01B11.2 | Protein CE-SULP-6 | Caenorhabditis elegans | 249 | 2.5×10-142 | 46.7 |
| BXYJ5.050219300 | BxSULP8 | CELE_ZK287.2 | Protein CE-SULP-8 | Caenorhabditis briggsae | 462 | 8.4×10-148 | 56.5 |
| BXYJ5.060009300 | BxSULP9 | CBG_23241 | Protein CBR-SULP-5 | Caenorhabditis elegans | 343 | 3.1×10-133 | 51.5 |
| BXYJ5.060037700 | BxSULP10 | CBG_23241 | Protein CBR-SULP-5 | Caenorhabditis briggsae | 303 | 1.9×10-129 | 48.9 |
Comparison of BxSULPs with other nematodes.
Comparative analysis of BxSULPs with homologs from nematodes like C. briggsae and C. elegans revealed both conserved and divergent evolutionary patterns. High BLASTP scores (above 100) and low E-values (below 1×10-10) indicated significant sequence similarity, implying shared conserved structural and functional domains, notably in the Sulfate Transporter and Anti-sigma factor antagonist (STAS) domain crucial for sulfate binding and transport. Notably, BxSULP8 exhibited a remarkable BLASTP score of 462 and an exceptionally low E-value of 8.4×10–148 compared to a STAS domain-containing protein in C. elegans. This high degree of sequence similarity suggested that BxSULP8 was highly conserved and likely played a crucial role in sulfate transport, similar to its orthologs in other nematodes. Given the critical function of the STAS domain in sulfate binding and transport, BxSULP8 might be particularly important for maintaining sulfate homeostasis in B. xylophilus. This warrants further functional investigations to elucidate its specific role in nematode biology. These differences underscored the importance of understanding the evolutionary and functional contexts of SULP across diverse nematode species, providing insights for targeted potential control strategies and targeted interventions.
3.2 Physicochemical analysis of BxSULPs
The average number of amino acids in BxSULPs was 202.2, with an average molecular mass of 22.71 kDa. The size of BxSULP varied, ranging from 17.63 kDa (BxSULP1) to 29.47 kDa (BxSULP2). Notably, BxSULP10, the largest protein, may exhibit more intricate functionalities or interactions. The isoelectric points (pI) of BxSULPs spanned from 6.51 (BxSULP6) to 10.19 (BxSULP9). Proteins with pI values near neutrality are anticipated to demonstrate enhanced stability under physiological conditions, whereas those with extreme pI values might serve specialized roles. BxSULP9 could participate in processes occurring in high pH environments. BxSULP6, characterized by the highest aliphatic index (119.51), may exhibit greater resistance to thermal denaturation compared to other BxSULPs, a critical feature for its potential function in extreme nematode environments. BxSULP2, BxSULP5, BxSULP9, and BxSULP10 were hydrophilic, while the remaining proteins were hydrophobic. BxSULP4, BxSULP7, and BxSULP8 demonstrated stability, contrasting with the instability observed in the other proteins. This suggested that protein may possess more robust structures, potentially ensuring sustained functionality over time (Table 3).
Table 3
| Protein name | Number of amino acids | Molecular weight/kDa | Isoelectric point1 | Aliphatic index2 | Hydropathicity3 | Instability index4 | Stability |
|---|---|---|---|---|---|---|---|
| BxSULP1 | 158 | 17.63 | 9.96 | 85.19 | 0.068 | 51.61 | Unstable |
| BxSULP2 | 257 | 29.47 | 9.60 | 94.79 | -0.125 | 46.28 | Unstable |
| BxSULP3 | 200 | 22.32 | 9.30 | 99.95 | 0.402 | 46.92 | Unstable |
| BxSULP4 | 167 | 18.75 | 8.58 | 105.03 | 0.100 | 34.51 | Stable |
| BxSULP5 | 173 | 20.23 | 10.09 | 82.89 | -0.124 | 53.37 | Unstable |
| BxSULP6 | 225 | 25.52 | 6.51 | 119.51 | 0.510 | 43.38 | Unstable |
| BxSULP7 | 181 | 19.76 | 8.22 | 114.64 | 0.893 | 37.56 | Stable |
| BxSULP8 | 227 | 24.71 | 9.90 | 109.91 | 0.763 | 31.44 | Stable |
| BxSULP9 | 172 | 19.50 | 10.19 | 72.67 | -0.247 | 59.18 | Unstable |
| BxSULP10 | 262 | 29.18 | 7.74 | 75.92 | -0.049 | 49.91 | Unstable |
Physicochemical properties analysis of BxSULPs.
1 Isoelectric point (pI): The isoelectric point (pI) is the pH at which a protein carries no net electrical charge.
2 Aliphatic indexes: The Aliphatic Index is a measure that describes the relative volume occupied by aliphatic side chains in a protein.
3 Hydropathicity: Hydropathicity is a relative value used to measure the hydrophobicity or hydrophilicity of a molecule. Positive values indicate hydrophobicity, while negative values indicate hydrophilicity.
4 Instability index: The Instability Index is a metric used to evaluate the stability of proteins. If the Instability Index is less than 40, the protein is considered stable. If the Instability Index is greater than 40, the protein is unstable.
3.3 Phylogenetic and conserved motif analysis of BxSULPs
We performed phylogenetic analysis on the full-length protein sequences. The conserved STAS and Sulfate_transp domains were subsequently mapped onto the alignment to verify motif conservation. The phylogenetic analysis revealed distinct clustering patterns among SULPs, indicating conserved and divergent evolutionary trajectories. The grouping of BxSULPs with homologs from closely related nematodes like C. briggsae and C. elegans suggested significant sequence conservation and shared ancestry. For example, BxSULP7 clustered closely with CbrSULP6 and CeSULP6, implying a strong evolutionary relationship, while functional equivalence remains to be experimentally determined. However, notable divergence was observed among SULPs from different nematode species. BxSULP2 exhibited a unique branching pattern, suggesting potential functional or adaptive divergence specific to B. xylophilus, possibly linked to its distinct ecological niche. The phylogenetic tree categorized BxSULP into distinct clades, denoted by different colors. Purple represents the SULP I family, red the SULP II family, and green the SULP III family, highlighting separate evolutionary lineages within the SULP family. SULP4, SULP5, SULP6, SULP8 of B. xylophilus, SULP3 of C. briggsae and SULP8 of B. malayi form a well-supported clade. Similarly, SULP1, SULP9 of B. xylophilus, SULP7 of C. elegans and SULP7 of C. briggsae clustered together in another distinct clade. SULP3, SULP7, SULP10 of B. xylophilus, SULP1, SULP2, SULP4, SULP5, SULP6, SULP8 of C. elegans, SULP2, SULP4, SULP5, SULP6, SULP8 of C. briggsae and SULP6 of B. malayi clustered in a separate branch, supporting the concept of gene duplication and subsequent functional diversification. These clades were also strongly supported by high bootstrap values, affirming their robustness (Figure 1).
Figure 1
In the context of sulfate permeases, the motifs are likely involved in critical functions such as sulfate binding, transmembrane transport, and protein stability (Supplementary Figure S1). Analysis determined the compositions of introns and exons in Bxsulps, revealing that the number of exons ranged from 7 to 16 across the 10 genes. Specifically, BxSULP3, BxSULP8 and BxSULP10 had only 7 exons, while BxSULP5 exhibited the highest number with 14 exons. Examination of structural motifs in the BxSULPs unveiled shared motifs among these proteins, as corroborated by the analysis of conserved motif structures (Supplementary Figure S2). The figure illustrated the conserved motifs identified within BxSULPs, emphasizing regions of sequence similarity across the proteins. These motifs are short, conserved sequences that are likely involved in critical functions such as sulfate binding, transmembrane transport, and protein stability. The clustering of these motifs indicates regions of high conservation, suggesting potential functional roles.
3.4 Chromosomal distribution of Bx-sulps
The chromosomal localization of the 10 Bx-sulps in B. xylophilus was determined using TBtools v2.056, precisely mapping the locations of these channel-encoding genes across various chromosomes (Figure 2). Analysis of chromosomal localization revealed that the genes encoding the 10 SULP proteins were distributed among five chromosomes of B. xylophilus. Specifically, on chromosome 1, Bx-sulp1 was situated in the 3′ region. Bx-sulp2 and Bx-sulp3 were positioned in the middle of chromosome 2. Furthermore, Bx-sulp4 and Bx-sulp5 were identified in the 3′ region of chromosome 3. Within chromosome 5, Bx-sulp6 and Bx-sulp7 were located in the central region, and Bx-sulp8 was found in the 3′ region. Additionally, Bx-sulp9 and Bx-sulp10 were situated in the 5′ region of chromosome 6, with no genes present on chromosome 4. The dispersion of Bx-sulps across multiple chromosomes indicates their non-clustered distribution throughout the genome. This dispersion had implications for gene regulation, as genes located on distinct chromosomes could be governed by diverse regulatory mechanisms and environmental factors.
Figure 2
3.5 Protein structure analysis of BxSULPs
Analysis of the three-dimensional structures of 10 BxSULPs revealed similar proportions of structural elements, including α-helices (32.56%–60.22%), β-folds (6.11%–23.12%), β-turns (0.60%–11.05%), and random coils (28.73%–53.82%). The transmembrane segment of BxSULPs predominantly comprised α-helices, facilitating their integration into the cell membrane to form channel structures for sulfate transport. Despite variations in the number and positioning of transmembrane segments across different proteins, a general conservation of these regions was observed. While BxSULPs exhibited sequence diversity, their three-dimensional structures displayed substantial similarity, indicating the preservation of core functional features essential for sulfate transport (Table 4; Figure 3). The Ramachandran plot was used to assess the three-dimensional structures of 10 proteins. The dihedral angles of BxSULPs residues were found to be located in the yellow favored regions, and the spatial structure was found to be over 90% stable, indicating a high degree of confidence in the tertiary structure. The Ramachandran plot analysis confirms the high stability of the three-dimensional structures of BxSULPs, with the majority of residues falling within the core regions, indicating a high degree of confidence in the predicted protein conformations (Figure 4).
Table 4
| Protein ID | Number of N-glycosylation sites | Number of amino acids | Percentage of number (%) | |||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| T | S | Y1 | α-helix2 | β-fold3 | β-turn4 | Random coil5 | α-helix | β-fold | β-turn | Random coil | ||
| BxSULP1 | 3 | 5 | 14 | 0 | 78 | 22 | 11 | 47 | 49.37 | 13.92 | 6.96 | 29.75 |
| BxSULP2 | 1 | 9 | 17 | 2 | 107 | 46 | 9 | 95 | 41.63 | 17.90 | 3.50 | 36.96 |
| BxSULP3 | 3 | 4 | 13 | 1 | 95 | 15 | 7 | 83 | 47.50 | 7.50 | 3.50 | 41.50 |
| BxSULP4 | 0 | 4 | 13 | 3 | 90 | 26 | 1 | 50 | 53.89 | 15.57 | 0.60 | 29.94 |
| BxSULP5 | 1 | 2 | 15 | 3 | 67 | 40 | 9 | 57 | 38.73 | 23.12 | 5.20 | 32.95 |
| BxSULP6 | 2 | 5 | 6 | 6 | 109 | 25 | 5 | 86 | 48.44 | 11.11 | 2.22 | 38.22 |
| BxSULP7 | 1 | 4 | 18 | 0 | 109 | 15 | 5 | 52 | 60.22 | 8.29 | 2.76 | 28.73 |
| BxSULP8 | 1 | 7 | 7 | 2 | 112 | 32 | 5 | 78 | 49.34 | 14.10 | 2.20 | 34.36 |
| BxSULP9 | 1 | 10 | 11 | 2 | 56 | 37 | 19 | 60 | 32.56 | 21.51 | 11.05 | 34.88 |
| BxSULP10 | 2 | 9 | 12 | 5 | 98 | 16 | 7 | 141 | 37.40 | 6.11 | 2.67 | 53.82 |
Secondary structure motifs of BxSULPs.
1 T, S, Y: T is Threonine, S is Serine, and Y is Tyrosine, quantity share is the number of amino acids that make up each secondary structure as a percentage of the total number of amino acids.
2 α-helix: The α-helix is a common secondary structure in proteins, formed by the right-handed coiling of a polypeptide chain around a central axis.
3 β-fold: The β-fold is a secondary structure formed by the parallel or antiparallel alignment of polypeptide chains.
4 β-turn: The β-turn is a common type of turn structure, typically composed of four amino acid residues.
5 Random coils: Random coil refers to the parts of a protein that cannot be classified as α-helices, β-fold, or other regular secondary structures.
Figure 3
Figure 4
The topological structure analysis indicated that BxSULP2, BxSULP3, BxSULP4, BxSULP6, BxSULP7, BxSULP8, and BxSULP10 were transmembrane proteins characterized by the presence of transmembrane domains enabling their integration into the cell membrane, primarily facilitating sulfate transmembrane transport. Conversely, BxSULP1, BxSULP5, and BxSULP9 were intra-membrane proteins. Although these proteins lacked transmembrane domains, they can be secreted outside the cell via signaling peptides. These signal peptides were short, hydrophobic regions at the N-terminus of the protein that guided the protein through the endoplasmic reticulum and into the extracellular space. These intra-membrane proteins operate within the extracellular matrix or engage in intercellular signaling upon secretion. Once secreted, these proteins can interact with other molecules in the extracellular matrix or participate in intercellular signaling pathways. Their functions might include binding to extracellular ligands, modulating cell-cell interactions, or influencing the local environment around the cell. Notably, BxSULP1 and BxSULP3 exhibited the highest number of N-glycosylation sites, with 3 sites each, enhancing their stability and functionality in membrane-associated or secreted contexts. This suggested a pivotal role in sulfate transport and interaction with the extracellular environment. The elongated signal peptide of BxSULP7 suggested more precise localization within the cell membrane, thereby enhancing its efficacy in sulfate transport. Despite variations in transmembrane structure, glycosylation patterns, and signal peptides among the BxSULPs, they exhibited a degree of structural similarity. These structural variations suggested functional partitioning. The multi-pass trans-membrane proteins were predicted to form the core sulfate-transport channels. In contrast, the secreted is formed that contain signal peptides might function as extracellular sulfate-binding or sulfate-donating molecules, modulate local sulfur availability or serving as signaling ligands during host invasion and diapause (Figure 5).
Figure 5
3.6 Expression profiles of Bx-sulps in different stages
To analyze gene expression in different developmental stages of B. xylophilus, we employed reverse-transcription quantitative polymerase chain reaction (RT-qPCR). By analyzing gene expression data from B. xylophilus under various nematode states and stress conditions, the expression data of 10 Bx-sulps were obtained. Significant variations in Bx-sulps expression across different nematode states (Figure 6). Specifically, Bx-sulps exhibited low expression levels during egg incubation and fourth instar larva (L4). Only two genes, Bx-sulp3 and Bx-sulp5, were highly expressed in second instar larva (L2). Bx-sulp1, Bx-sulp2 and Bx-sulp3 were highly expressed in third instar larva (L3). Notably, Bx-sulp4, Bx-sulp8, Bx-sulp9, and Bx-sulp10 showed heightened expression in dauer third-instar larva (D3), crucial for maintaining essential physiological functions during diapause and enhancing resistance. The overexpression of Bx-sulps potentially bolstered resistance during diapause, with sulfate likely playing a role in cell membrane stabilization and enhancing cell tolerance to extreme conditions. Moreover, Bx-sulp1, Bx-sulp6 and Bx-sulp7 were prominently expressed in female adult nematodes and male adult nematodes, suggesting involvement in sex differentiation. While L2, L3, and L4 stages exhibited similar sulp genes expression patterns, the transition to adulthood, either male or female, was marked by significant changes in gene expression profiles.
Figure 6
4 Discussion
In recent years, an increasing number of researchers had been investigating the molecular pathways of B. xylophilus by RNA-sequence analysis (). RNA-sequence analysis had primarily been employed to elucidate the pathogenicity, stress resistance, growth, and detoxification mechanisms of B. xylophilus (; ). For instance, Li conducted a comparative transcriptomic analysis to investigate the molecular interactions between the B. xylophilus and β-pinene, a secondary metabolite produced by pine as a defense against nematode infection (). The study revealed that a significant proportion of the differentially expressed genes were associated with cytochrome P450 (CYP), UDP-glucuronide transferase (UDT), and short-chain dehydrogenase (SDR) genes. Similarly, Lu analyzed transcriptome data of B. xylophilus following exposure to emamectin benzoate (). The analysis indicated that genes encoding pectate lyases and β-1,4-endoglucanases were down-regulated, while those related to glutamate-gated chloride channels (GluCls), γ-aminobutyric acid type β receptors (GABAB), and ATP-binding cassette transporters (ABC) were upregulated. These genes played crucial roles in the embryonic and larval development, reproduction, as well as the nervous and motor systems of B. xylophilus, making them potential targets for B. xylophilus control strategies. Further exploration of specific gene families within B. xylophilus will offer a solid theoretical foundation for the development of targeted control measures. This study provides theoretical basis for developing control strategies based on B. xylophilus by studying BxSUIPs.
The sulfate permease (SULP) family, commonly associated with sulfate transporters (SULTR) which is prevalent in plants and certain microorganisms. SULTR played a crucial role in sulfate absorption and transportation, serving as a carrier protein essential for active sulfate transport in plants (). The SULTR protein comprises 12 membrane-spanning domains (TMDs) at the N-terminal, a sulfur_transp domain at the N-terminal, and a STAS domain at the C-terminal region. These domains are directly involved in sulfate binding, transmembrane transport, and the maintenance of SULTR protein activity and stability (). In Heterodera glycines, an SULTR-type sulphate transporter (HgSULTR1;2) was transcriptionally up-regulated 5- to 12-fold within 24 h of root invasion and RNAi-mediated knock-down reduced both cyst formation and egg hatching by >60%, demonstrating that efficient sulphate uptake is essential for parasitic success (). Sulphate limitation triggers a global stress response in nematodes. In Caenorhabditis elegans, sulphate starvation induced dauer formation, implying that precise control of sulphate availability via SULP transporters can modulate developmental decisions and reproduction in B. xylophilus (). In other plant pathogens and symbionts, disruption of sulphate permeases reduced virulence. For example, deletion of the sulphate permease gene cysP in the bacterium Xanthomonas campestris attenuated growth in plant (). These findings suggested that sulp genes were evolutionarily conserved and functionally diversified across taxa.
Cysteine and methionine synthesis depend on a continuous supply of sulfate; any modulation of SULP-mediated transport could directly impact nematode protein metabolism. As a sulfur-containing compound, ergothioneine can significantly prolong the lifespan of C. elegans, improve its mobility and resistance to adverse environments (). The diapause stage is a survival strategy of B. xylophilus to endure unfavorable conditions like low temperatures and drought, allowing it to remain active across seasons and environments, thereby expanding its transmission range (). The high expression of Bx-sulps in the dauer third-instar larva (D3) of B. xylophilus indicated its importance in maintaining essential physiological functions and enhancing resistance during diapause, aligning with the notion that sulfur compounds can enhance nematode resilience to adverse environments. In the D3 stage, sulfate likely aids in stabilizing cell membranes and increasing cell tolerance to extreme conditions (). Moreover, the significant expression of Bx-sulp6 and Bx-sulp7 in female adult nematodes (FM) and male adult nematodes (M) suggested their involvement in sex differentiation, underscoring the potential role of sulfate transport in reproductive processes with implications for nematode population dynamics and pathogenicity. The varying expression patterns in different developmental stages emphasized the significance of elucidating the regulatory mechanisms that control Bx-sulps gene expression, this could facilitate the design of precise interventions to disrupt the nematode’s sulfate transport mechanisms. Future functional studies, including RNA interference (RNAi) and CRISPR-based gene editing, were warranted to dissect the individual contributions of each Bx-sulp genes to nematode development and pathogenicity.
The phylogenetic analysis displayed distinct clustering patterns, indicating that BxSULPs form well-supported clades with homologs from other nematodes, suggesting significant sequence conservation and shared ancestry. This observation implied that these proteins likely preserve essential functional features crucial for sulfate transport (). Conserved motifs were identified across the 10 BxSULPs, associated with key functions such as sulfate binding, transmembrane transport, and protein stability, further supporting their conserved functionality. The chromosomal localization of the 10 Bx-sulps revealed their dispersed distribution across five chromosomes. This dispersion had implications for gene regulation, genes situated on different chromosomes might be subject to diverse regulatory mechanisms and environmental influences (). The non-clustered distribution of Bx-sulps suggested a complex regulatory landscape, possibly involving diverse transcription factors and signaling pathways (). This intricate regulatory network might enable the nematode to finely regulate its sulfate transport mechanisms in response to varying environmental conditions and developmental stages (). The topological analysis indicated that BxSULP2, BxSULP3, BxSULP4, BxSULP6, BxSULP7, BxSULP8, and BxSULP10 were transmembrane proteins primarily facilitating sulfate transmembrane transport. In contrast, BxSULP1, BxSULP5, and BxSULP9 were intra-membrane proteins that might be secreted outside the cell through signaling peptides, localized within the extracellular matrix, or participate in intercellular signaling upon secretion. The hydropathicity analysis indicated that BxSULP5 and BxSULP9 exhibited hydrophilic properties, which can involve in extractable interactions or signaling processes, consistent with the results of topology analysis. The three-dimensional structural analysis of BxSULPs revealed that their transmembrane segments were primarily composed of α-helices. These α-helical structures facilitated protein integration into the cell membrane, crucial for forming functional channel configurations that facilitate efficient sulfate ion transport (). These findings provided a structural foundation for comprehending the function and regulation of BxSULPs.
The identification and comparative analysis of BxSULPs represent a significant advancement in comprehending the adaptive strategies of nematodes and their pathogenic interactions with pine hosts. This research elucidates the roles of the BxSULPs in the development and survival, enhancing the understanding of the molecular mechanisms governing development and stress responses of B. xylophilus. Disrupting the sulfate transport mechanisms of B. xylophilus can potentially impede its ability to thrive within the pine host environment, thereby mitigating its pathogenicity and spread. Furthermore, this study serves as a foundation for future investigations into the gene family of B. xylophilus and gene knockout experiments involving the SULP family.
5 Conclusions
In conclusion, this study represents the first comprehensive identification and characterization of the complete SULP gene family in B. xylophilus. The BxSULPs exhibit highly conserved structural features, indicative of evolutionary stability within nematodes, while minimal structural divergence among them suggests functional consistency and potential synergistic interactions. Expression profiling revealed that Bx-sulp4, Bx-sulp8, Bx-sulp9, and Bx-sulp10 were specifically up-regulated during the dauer (D3) stage, implicating these transporters in diapause maintenance and stress adaptation. Conversely, Bx-sulp1, Bx-sulp6 and Bx-sulp7 showed elevated expression in adult males or females, indicating a putative role in reproductive processes and sexual differentiation. The non-clustered genomic distribution and diverse membrane topologies of BxSULPs further provide a structural foundation for their stage-specific and environmentally modulated regulation. Collectively, these findings position the BxSULPs as promising molecular targets for RNAi-based or small-molecule intervention strategies aimed at disrupting critical biological processes such as diapause and fecundity, thereby offering novel prospects for the control of PWD. However, further investigation is required to elucidate the involvement of BxSULPs in the immune mechanisms of B. xylophilus. This study establishes a foundational framework for future functional analyses of BxSULPs and identifies potential targets for the development of targeted management strategies against B. xylophilus infestations.
Statements
Data availability statement
The raw sequencing data BXSULP1-10 generated from genes analyses in this study have been deposited in the NCBI database (PX242223-32).
Author contributions
HL: Data curation, Validation, Writing – original draft, Software, Formal Analysis. RW: Writing – original draft, Validation, Visualization, Methodology. NP: Writing – original draft, Validation, Visualization. SY: Project administration, Investigation, Funding acquisition, Resources, Writing – review & editing. JC: Visualization, Validation, Supervision, Writing – review & editing, Conceptualization. XH: Conceptualization, Writing – original draft.
Funding
The author(s) declare financial support was received for the research and/or publication of this article. The study supported by Southwest Forestry University Forestry major in Yunnan Province First-Class Construction Discipline (LXXK-2023M06); Yunnan Fundamental Research Projects (202401BD070001-115); Yunnan Provincial Plateau-Characteristic Science and Technology Project in Agriculture (202302AE090017).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2025.1630288/full#supplementary-material
References
1
AbelleiraA.PicoagaA.MansillaJ. P.AguinO. (2011). Detection of bursaphelenchus xylophilus, causal agent of pine wilt disease on pinus pinaster in Northwestern Spain. Plant Dis.95, 776. doi: 10.1094/PDIS-12-10-0902
2
AhmedW.YeW.PanJ.LiuS.JiW.ZhouS.et al. (2025). Evaluation the role of Luteibacter pinisoli DP2–30 in mitigating pine wilt disease caused by Bursaphelenchus xylophilus through modulation of host microbiome. Front. Plant Sci.16, 1515506. doi: 10.3389/fpls.2025.1515506
3
AliV.NozakiT. (2007). Current therapeutics, their problems, and sulfur-containing-amino-acid metabolism as a novel target against infections by “amitochondriate” protozoan parasites. Clin. Microbiol. Rev.20, 164–187. doi: 10.1128/CMR.00019-06
4
AnY.LiY.MaL.LiD.ZhangW.FengY.et al. (2023). Transcriptomic response of Pinus massoniana to infection stress from the pine wood nematode Bursaphelenchus xylophilus. Stress Biol.3, 50. doi: 10.1007/s44154-023-00131-z
5
BaileyT. L.JohnsonJ.GrantC. E.NobleW. S. (2015). The MEME suite. Nucleic Acids Res.43, W39–W49. doi: 10.1093/nar/gkv416
6
BailleulG.YangG.NicollC. R.AndreaM.FraaijeM. W.MascottiM. L. (2023). Evolution of enzyme functionality in the flavin-containing monooxygenases. Nat. Commun.14, 1042. doi: 10.1038/s41467-023-36756-x
7
BlomN.GammeltoftS.BrunakS. (1999). Sequence and structure-based prediction of eukaryotic protein phosphorylation sites 1–1 Edited by F. E. Cohen. J. Mol. Biol.294, 1351–1362. doi: 10.1006/jmbi.1999.3310
8
CantelliG.BatemanA.BrooksbankC.PetrovA. I.Malik-SheriffR. S.Ide-SmithM.et al. (2021). The european bioinformatics institute (EMBL-EBI) in 2021. Nucleic Acids Res.50, D11–D19. doi: 10.1093/nar/gkab1127
9
ChangY. N.JaumannE. A.ReichelK.HartmannJ.OliverD.HummerG.et al. (2019). Structural basis for functional interactions in dimers of SLC26 transporters. Nat. Commun.10, 2032. doi: 10.1038/s41467-019-10001-w
10
ChenC.ChenH.ZhangY.ThomasH. R.FrankM. H.HeY.et al. (2020). TBtools: an integrative toolkit developed for interactive analyses of big biological data. Molecluar Plant13, 1194–1202. doi: 10.1016/j.molp.2020.06.009
11
ChenJ.HaoX.TanR.LiY.WangB.PanJ.et al. (2022). Functional study on cytochrome P450 in response to L(-)-carvone stress in bursaphelenchus xylophilus. Genes (Basel)13, 1956. doi: 10.3390/genes13111956
12
ChenC.WuY.LiJ.WangX.ZengZ.XuJ.et al. (2023). TBtools-II: A “one for all, all for one” bioinformatics platform for biological big-data mining. Molecluar Plant16, 1733–1742. doi: 10.1016/j.molp.2023.09.010
13
ChenY.XiaoX.YangR.SunZ.YangS.ZhangH.et al. (2024). Genome-wide identification and expression-pattern analysis of sulfate transporter (SULTR) gene family in cotton under multiple abiotic stresses and fiber development. Funct. Integr. Genomics24, 108. doi: 10.1007/s10142-024-01387-y
14
DingY.ZhouX.ZuoL.WangH.YuD. (2016). Identification and functional characterization of the sulfate transporter gene GmSULTR1;2b in soybean. BMC Genomics17, 373. doi: 10.1186/s12864-016-2705-3
15
DuvaudS.GabellaC.LisacekF.StockingerH.IoannidisV.DurinxC. (2021). Expasy, the Swiss Bioinformatics Resource Portal, as designed by its users. Nucleic Acids Res.49, W216–w227. doi: 10.1093/nar/gkab225
16
El-GebaliS.MistryJ.BatemanA.EddyS. R.LucianiA.PotterS. C.et al. (2019). The Pfam protein families database in 2019. Nucleic Acids Res.47, D427–d432. doi: 10.1093/nar/gky995
17
GoekeS. P.BinderE. B. (2024). Linking environmental factors and gene regulation. eLife13, e96710. doi: 10.7554/eLife.96710
18
GoronzyI. N. (2025). Unravelling the complexity of gene regulation through multiplexed protein mapping. Nat. Rev. Mol. Cell Biol.26, 1. doi: 10.1038/s41580-025-00830-7
19
HachamY.KaplanA.CohenE.GalM.AmirR. (2025). Sulfur metabolism under stress: Oxidized glutathione inhibits methionine biosynthesis by destabilizing the enzyme cystathionine γ-synthase. J. Integr. Plant Biol.67, 87–100. doi: 10.1111/jipb.13799
20
HooverK.UzunovicA.GethingB.DaleA.LeungK.OstiguyN.et al. (2010). Lethal temperature for pinewood nematode, Bursaphelenchus xylophilus, in infested wood using microwave energy. J. Nematol.42, 101–110.
21
HoweK. L.BoltB. J.ShafieM.KerseyP.BerrimanM. (2016). WormBase ParaSite – a comprehensive resource for helminth genomics. Mol. Biochem. Parasitol.215, 2–10. doi: 10.1016/j.molbiopara.2016.11.005
22
HuB.JinJ.GuoA. Y.ZhangH.LuoJ.GaoG. (2015). GSDS 2.0: an upgraded gene feature visualization server. Bioinformatics31, 1296–1297. doi: 10.1093/bioinformatics/btu817
23
KangJ. S.LeeH.MoonI. S.LeeY.KohY. H.JeY. H.et al. (2009). Construction and characterization of subtractive stage-specific expressed sequence tag (EST) libraries of the pinewood nematode Bursaphelenchus xylophilus. Genomics94, 70–77. doi: 10.1016/j.ygeno.2009.03.001
24
KoichiroT.GlenS.SudhirK. (2021). MEGA11: molecular evolutionary genetics analysis version 11. Mol. Biol. Evol.38, 3022–3027. doi: 10.1093/molbev/msab120
25
LetunicI.BorkP. (2024). Interactive Tree of Life (iTOL) v6: recent updates to the phylogenetic tree display and annotation tool. Nucleic Acids Res.52, W78–W82. doi: 10.1093/nar/gkae268
26
LiY.FengY.WangX.CuiJ.DengX.ZhangX. (2020). Adaptation of pine wood nematode Bursaphelenchus xylophilus to β-pinene stress. BMC Genomics21, 478. doi: 10.1186/s12864-020-06876-5
27
LiH.WangR.PanJ.ChenJ.HaoX. (2025). Bioinformatics analysis of the glutamate-gated chloride channel family in bursaphelenchus xylophilus. Int. J. Mol. Sci.26, 3477. doi: 10.3390/ijms26083477
28
LiuX.LiangX.HanJ.CuiY.LeiM.WangB.et al. (2025). Genome-wide identification and transcriptome analysis of the cytochrome P450 genes revealed its potential role in the growth of Flammulina filiformis. BMC Genomics26, 346. doi: 10.1186/s12864-025-11555-4
29
LuF.GuoK.ChenA.ChenS.LinH.ZhouX. (2020a). Transcriptomic profiling of effects of emamectin benzoate on the pine wood nematode Bursaphelenchus xylophilus. Pest Manage. Sci.76, 747–757. doi: 10.1002/ps.5575
30
LuS.WangJ.ChitsazF.DerbyshireM. K.GeerR. C.GonzalesN. R.et al. (2020b). CDD/SPARCLE: the conserved domain database in 2020. Nucleic Acids Res.48, D265–d268. doi: 10.1093/nar/gkz991
31
LuoA.QiaoH.ZhangY.ShiW.HoS. Y.XuW.et al. (2010). Performance of criteria for selecting evolutionary models in phylogenetics: a comprehensive study based on simulated datasets. BMC Evolutionary Biol.10, 242. doi: 10.1186/1471-2148-10-242
32
MamiyaY. (1983). Pathology of the pine wilt disease caused by bursaphelenchus xylophilus. Annu. Rev. Phytopathol.21, 201–220. doi: 10.1146/annurev.py.21.090183.001221
33
MarcoB.StefanB.AndrewW.KonstantinA.GabrielS.TobiasS.et al. (2014). SWISS-MODEL: modelling protein tertiary and quaternary structure using evolutionary information. Nucleic Acids Res.42, W252–W258. doi: 10.1093/nar/gku340
34
MarietouA.RøyH.JørgensenB. B.KjeldsenK. U. (2018). Sulfate transporters in dissimilatory sulfate reducing microorganisms: A comparative genomics analysis. Front. Microbiol.9, 309. doi: 10.3389/fmicb.2018.00309
35
PetrovicD.SladeL.PaikopoulosY.D’AndreaD.SavicN.StancicA.et al. (2025). Ergothioneine improves healthspan of aged animals by enhancing cGPDH activity through CSE-dependent persulfidation. Cell Metab.37, 542–556. doi: 10.1016/j.cmet.2024.12.008
36
RinuK.KoushikR. G.NejraS.MarieD. L. (2020). ATP-binding cassette (ABC) transporter genes in plant-parasitic nematodes: an opinion for development of novel control strategy. Front. Plant Sci.11, 582424. doi: 10.3389/fpls.2020.582424
37
ShinyaR.KirinoH.MorisakaH.Takeuchi-KanekoY.FutaiK.UedaM. (2021). Comparative secretome and functional analyses reveal glycoside hydrolase family 30 and cysteine peptidase as virulence determinants in the pinewood nematode bursaphelenchus xylophilus. Front. Plant Sci.12, 640459. doi: 10.3389/fpls.2021.640459
38
ShinyaR.MorisakaH.KikuchiT.TakeuchiY.UedaM.FutaiK. (2013). Secretome analysis of the pine wood nematode bursaphelenchus xylophilus reveals the tangled roots of parasitism and its potential for molecular mimicry. PloS One8, e67377. doi: 10.1371/journal.pone.0067377
39
TanakaS. E.DayiM.MaedaY.TsaiI. J.TanakaR.BlighM.et al. (2019). Stage-specific transcriptome of Bursaphelenchus xylophilus reveals temporal regulation of effector genes and roles of the dauer-like stages in the lifecycle. Sci. Rep.9, 6080. doi: 10.1038/s41598-019-42570-7
40
TangJ.MaR.ZhuN.GuoK.GuoY.BaiL.et al. (2020). Bxy-fuca encoding α-L-fucosidase plays crucial roles in development and reproduction of the pathogenic pinewood nematode, Bursaphelenchus xylophilus. Pest Manage. Sci.76, 205–214. doi: 10.1002/ps.5497
41
ThibautG.StéphaneD.VincentN. (2015). VirHostNet 2.0: surfing on the web of virus/host molecular interactions data. Nucleic Acids Res.43, D583–D587. doi: 10.1093/nar/gku1121
42
TianH.KoskiT. M.ZhaoL.LiuZ.SunJ. (2022). Invasion History of the Pinewood Nematode Bursaphelenchus xylophilus Influences the Abundance of Serratia sp. in Pupal Chambers and Tracheae of Insect-Vector Monochamus alternatus. Front. Plant Sci.13, 856841. doi: 10.3389/fpls.2022.856841
43
TogashiK.AkbulutS.MatsunagaK.SugimotoH.YanagisawaK. (2019). Transmission of Bursaphelenchus xylophilus between Monochamus alternatus and Monochamus saltuarius through interspecific mating behaviour. J. Appl. Entomology143, 483–486. doi: 10.1111/jen.12604
44
UlrichO.AhrensC. H.SebastianM.BerndW. (2014). Protter: interactive protein feature visualization and integration with experimental proteomic data. Bioinf. (Oxford England)30, 884–886. doi: 10.1093/bioinformatics/btt607
45
VlaarL. E.BertranA.RahimiM.DongL.KammengaJ. E.HelderJ.et al. (2021). On the role of dauer in the adaptation of nematodes to a parasitic lifestyle. Parasit Vectors14, 554. doi: 10.1186/s13071-021-04953-6
46
WangL.ChenK.ZhouM. (2021). Structure and function of an Arabidopsis thaliana sulfate transporter. Nat. Commun.12, 4455. doi: 10.1038/s41467-021-24778-2
47
WangB.WenR.MaoX.ChenJ.HaoX. (2024). Unveiling the co-expression network and molecular targets behind rotenone resistance in the Bursaphelenchus xylophilus. Ecotoxicology Environ. Saf.287, 117293. doi: 10.1016/j.ecoenv.2024.117293
48
YuL.WangY.WangX.HanS.WangL.WangX. (2025). Transcriptomic, metabonomic and proteomic analyses reveal that terpenoids and flavonoids are required for Pinus koraiensis early defence against Bursaphelenchus xylophilus infection. BMC Plant Biol.25, 185. doi: 10.1186/s12870-025-06192-8
49
ZhangJ.LiY.YuanH.SunB.LiH. (2016). Biological control of the cereal cyst nematode (Heterodera filipjevi) by Achromobacter xylosoxidans isolate 09X01 and Bacillus cereus isolate 09B18. Biol. Control92, 1–6. doi: 10.1016/j.biocontrol.2015.08.004
50
ZhangW.YuH.LvY.BushleyK. E.WickhamJ. D.GaoS.et al. (2020). Gene family expansion of pinewood nematode to detoxify its host defence chemicals. Mol. Ecol.29, 940–955. doi: 10.1111/mec.15378
51
ZhangY.ZhaoR.JingT.LinS.DingX. (2024). Identification and transcriptome analysis of bursaphelenchus xylophilus with excellent low temperature resistance. Int. J. Mol. Sci.25, 13732. doi: 10.3390/ijms252413732
52
ZhouJ.DongY.LiuY.HuangY.JiangW.ZhengX.et al. (2024). Identification and expression analysis of sulfate transporter genes family and function analysis of gmSULTR3;1a from soybean. Int. J. Mol. Sci.25, 9080. doi: 10.3390/ijms25169080
Summary
Keywords
Bursaphelenchus xylophilus, sulfate permease family, gene and protein structure, bioinformatic analysis, gene expression
Citation
Li H, Wang R, Pu N, Yang S, Chen J and Hao X (2025) Analysis of the sulfate permease family in Bursaphelenchus xylophilus in the nematode development and stress adaptation. Front. Plant Sci. 16:1630288. doi: 10.3389/fpls.2025.1630288
Received
17 May 2025
Accepted
23 September 2025
Published
13 October 2025
Volume
16 - 2025
Edited by
Claudia Sl Vicente, University of Évora, Portugal
Reviewed by
Laith Khalil Tawfeeq Al-Ani, Universiti Sains Malaysia, Malaysia
Margarida Espada, University of Évora, Portugal
Santisree Parankusam, Sedia Biosceinces, United States
Updates
Copyright
© 2025 Li, Wang, Pu, Yang, Chen and Hao.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xin Hao, xinhao@swfu.edu.cn; Jie Chen, jiechen@swfu.edu.cn
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.