Next Article in Journal
BMP3 Affects Cortical and Trabecular Long Bone Development in Mice
Next Article in Special Issue
The Regulation of Nitrate Reductases in Response to Abiotic Stress in Arabidopsis
Previous Article in Journal
Induction of Isochromanones by Co-Cultivation of the Marine Fungus Cosmospora sp. and the Phytopathogen Magnaporthe oryzae
Previous Article in Special Issue
Bidirectional lncRNA Transfer between Cuscuta Parasites and Their Host Plant
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 2
10.3390/ijms23020784
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Co-Silencing of the Voltage-Gated Calcium Channel β Subunit and High-Voltage Activated α1 Subunit by dsRNA Soaking Resulted in Enhanced Defects in Locomotion, Stylet Thrusting, Chemotaxis, Protein Secretion, and Reproduction in Ditylenchus destructor

by
Mingwei An
†,
Xueling Chen
†,
Zhuhong Yang
,
Jianyu Zhou
,
Shan Ye
* and
Zhong Ding
*
College of Plant Protection, Hunan Agricultural University, Changsha 410128, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Int. J. Mol. Sci. 2022, 23(2), 784; https://doi.org/10.3390/ijms23020784
Submission received: 25 November 2021 / Revised: 7 January 2022 / Accepted: 8 January 2022 / Published: 11 January 2022
(This article belongs to the Collection Recent Advances in Plant Molecular Science in China 2021)
Browse Figures
Review Reports Versions Notes

Abstract

:
The voltage-gated calcium channel (VGCC) β subunit (Cavβ) protein is a kind of cytosolic auxiliary subunit that plays an important role in regulating the surface expression and gating characteristics of high-voltage-activated (HVA) calcium channels. Ditylenchus destructor is an important plant-parasitic nematode. In the present study, the putative Cavβ subunit gene of D. destructor, namely, DdCavβ, was subjected to molecular characterization. In situ hybridization assays showed that DdCavβ was expressed in all nematode tissues. Transcriptional analyses showed that DdCavβ was expressed during each developmental stage of D. destructor, and the highest expression level was recorded in the third-stage juveniles. The crucial role of DdCavβ was verified by dsRNA soaking-mediated RNA interference (RNAi). Silencing of DdCavβ or HVA Cavα1 alone and co-silencing of the DdCavβ and HVA Cavα1 genes resulted in defective locomotion, stylet thrusting, chemotaxis, protein secretion and reproduction in D. destructor. Co-silencing of the HVA Cavα1 and Cavβ subunits showed stronger interference effects than single-gene silencing. This study provides insights for further study of VGCCs in plant-parasitic nematodes.

1. Introduction

Calcium channels are glycoproteins or lateral glycosylated protein complexes that are embedded in membrane lipids. These channels are widely present on the surface of the cell membrane and form a functional unit or pore. They have a high selective permeability for Ca2+ and control ion transport and cell function. Voltage-gated calcium channels (VGCCs) are important calcium channels that participate in many excitatory processes of cells, such as muscle contraction regulation, neurotransmitter release, gene regulation, and neuronal migration [1].
L-, N-, P/Q-, and R-type VGCCs consist of a central pore-forming α1 subunit and α2δ, β, and γ related auxiliary subunits [2,3]. The α1 subunit (Cavα1) comprises of four homologous domains constituted by six putative transmembrane segments (S1′S6) [4,5]. The β subunit (Cavβ) is a cytosolic auxiliary subunit that binds to the proximal part of the intracellular loop between domains Ⅰ and Ⅱ of the α1 subunit [6,7]. Cavα1 and Cavβ associate through a high-affinity interaction between the α-interaction domain (AID) and AID-binding pocket (ABP) [8].
Cavβ consists of five typical regions. The amino and carboxyl ends vary in both length and amino acid sequence. The Src homology 3 (SH3) and guanylate kinase (GK) domains are highly conserved and are connected by a weakly conserved HOOK domain, forming an N-SH3-HOOK-GK-C structure. The SH3-HOOK-GK core structure recapitulates many critical functions of Cavβ subunits, resulting in their classification as members of the membrane-associated GK protein family. Moreover, there are residues in the GK domain that interact with the AID [3,9,10]. As an auxiliary subunit, Cavβ has many functions; especially in regulating the surface expression and gating characteristics of high-voltage-activated (HVA) calcium channels [3,11]. Although Cavα1 is the principal component of VGCCs, for the proper trafficking and functioning, the channels require Cavβ [12].
There are four subfamilies of Cavβs (β1–β4) that are encoded by four distinct genes, each with splice variants, that have been identified in vertebrates. In Caenorhabditis elegans, two genes encoding β subunits, namely CCB-1 and CCB-2, were identified by genetic analysis and a homology search [13]. CCB-1 is necessary for viability [14]. However, there are few detailed reports on the Cavβ gene in plant-parasitic nematodes. Therefore, the study of Cavβ in Ditylenchus destructor will help us understand the role of VGCCs in plant-parasitic nematodes.
D. destructor develops and reproduces at 5–34 °C and each female lays 100–200 eggs in a lifetime [15]. It is easy to be cultured on fungi or potatoes and its life cycle is relatively short. Under suitable temperature and humidity conditions, the nematodes can complete a generation in 20 to 30 days [16]. D. destructor is one of the important pathogens of sweet potato (Ipomoea batatas) and causes significant yield loss of 20% to 50%, and up to 100% under heavy infestations [17,18]. In previous studies, HVA L-type (DdCα1D) and non-L-type (DdCα1A) VGCC α1 subunits of D. destructor were cloned, and their expression and the effects of their knockdown on nematode activity, including locomotion, chemotaxis, and reproduction, were studied [19]. Although some progress has been made in research on the function of Cavα1 subunits in D. destructor, the function of HVA Cavα1 in nematode stylet thrusting, protein secretion, and interaction with Cavβ remains unclear.
In the current study, we cloned the full-length cDNA encoding Cavβ of D. destructor and investigated the expression level and localization of Cavβ via qPCR and in situ hybridization. Using RNA interference (RNAi), we also further explored the effects of HVA Cavα1 and Cavβ gene silencing on locomotion, chemotaxis, and reproduction of D. destructor and determined the effects of these genes on stylet thrusting and protein secretion. All the data that were generated provided additional clues to improve our understanding of VGCCs in D. destructor.

2. Results

2.1. Primary Structures of Cavβ of D. destructor

The full-length cDNA sequence of the D. destructor Cavβ, named DdCavβ, was obtained via PCR amplification (GenBank number: MN867027). The DdCavβ cDNA was 2036 bp long with a 171-bp 5′ untranslated region (UTR), a 1530-bp open reading frame (ORF), and a 335-bp 3′ UTR. The sequence was predicted to encode a protein of 509 amino acids, with a molecular mass of 56.88 kDa and a pI of 8.157. It contained an SL1 sequence at its 5’ end and a poly-A tail at its 3’ end (Figure 1A). The structure of DdCavβ protein was predicted by a voltage-dependent L-type calcium channel subunit beta-1 as a homology template, and the homology was 67.99% (Figure 1B).
DdCavβ was similar to the conventional β-subunit that is reported in vertebrates and invertebrates, including C. elegans, and contained five typical regions: amino (NH2) and carboxyl (COOH) termini with variable lengths and amino acid sequences, a highly conserved SH3 domain, a GK domain, and a variable and flexible HOOK region (Figure 1). As shown in Figure 1, DdCavβ had residues that were involved in interactions with the AID of Cavα1. The DdCavβ protein was more homologous to nematode (C. elegans) CCB-1 (68.5%) and filaria (Brugia malayi) BMA-CCB-1 (67.5%) than to Cavβ in vertebrates (43.9–48.6%) (Table 1). DdCavβ showed 13.3% sequence identity with C. elegans CCB-2.

2.2. Homology and Phylogenetic Analysis of DdCavβ

To investigate the evolutionary relationships of Cavβ, the neighbor-joining method in MEGA 7.0 software was used to construct phylogenetic trees comprising of the amino acid sequences of DdCavβ and other Cavβ protein sequences from invertebrates and vertebrates. As shown in the phylogenetic tree, conventional and variant β subunits were clustered into two large groups (Figure 2). The D. destructor Cavβ grouped with C. elegans CCB-1 and the Loa loa, B. malayi, Drosophila melanogaster, Heterololigo bleekeri, and Lymnaea stagnalis β subunits to form a large clade with conventional β-subunits that are found in many invertebrates. The Cavβ of vertebrates formed a large clade with the variant β subunits of Homo sapiens (β1–β4), Oryctolagus cuniculus, Rattus norvegicus, and Xenopus laevis. C. elegans CCB-2 formed a single, small clade.

2.3. Tissue Localization of DdCavβ

To examine the site of DdCavβ mRNA expression in the tissue of D. destructor, in situ hybridization was performed. As shown in Figure 3, hybridization signals were detected in all the nematode tissues (Figure 3E–H), and hybridization with the control sense probe yielded no hybridization signals (Figure 3A–D). The wide distribution of DdCavβ indicates the distribution pattern of HVA calcium channels in D. destructor, as the β subunit is always connected to the HVA Cavα1 subunit [3].

2.4. Expression of DdCavβ in Nematodes at Different Stages

The stage-specific expression of the DdCavβ transcript was analyzed via qPCR. DdCavβ was differentially expressed in D. destructor at different stages. With the expression level in second-stage juveniles (J2s) used as a standard, the expression of the β subunit gene was significantly upregulated from the egg stage to the third-stage juveniles (J3s) and then significantly downregulated in the fourth-stage juveniles (J4s). The expression in the J3 stage was the highest, and that in the J4 stage was the lowest (Figure 4).

2.5. Efficacy of Silencing HVA Cavα1 and Cavβ Genes by dsRNA Soaking

To further explore the role of the HVA Cavα1 and Cavβ, effective and specific silencing of these genes by RNAi was performed. We first evaluated the individual silencing of DdCα1D, DdCα1A, and DdCavβ and used a mixture of dsRNAs to target both DdCα1D or DdCα1A, and DdCavβ to evaluate the function of VGCCs.
Soaking with dsCavβ specifically silenced the expression of DdCavβ by 47.0% compared to the control, while the expression of DdCα1D and DdCα1A was not changed (Figure 5A). Similar results were observed in the dsCα1D and dsCα1A treatments. The expression of DdCα1D and DdCα1A was reduced by 53.1% and 47.6%, respectively, compared to that in the dsgfp treatment, while the expression levels of other subunit genes were not changed (Figure 5B,C). When the nematodes were treated with a mixture of dsCα1D and dsCavβ, the relative expression levels of DdCα1D and DdCavβ were significantly downregulated by 57.4% and 50.9%, respectively, compared to those in the dsgfp treatment, while the expression of DdCα1A showed no significant change (Figure 5D). Under co-treatment with dsCα1A and dsCavβ, the relative expression levels of DdCα1A and DdCavβ were significantly downregulated by 49.8% and 55.2%, respectively, compared to those in the dsgfp treatment, while the expression of DdCα1D was not significantly changed (Figure 5E).

2.6. Analysis of the Knockdown Phenotype

2.6.1. Observation of Locomotion Activity after RNAi

Nematode locomotion behaviors were tested via the migration assay in which the numbers of worms migrating through the columns and into the collection vials were counted at 6 h, 14 h, and 24 h. During the whole observation period, the rate of migration of nematodes that were treated with each target gene dsRNA was significantly lower than that of nematodes that were treated with dsgfp. The lowest nematode migration rates were found under the simultaneous treatment with dsCα1D and dsCavβ, with values of 12%, 33%, and 47% at 6 h, 14 h, and 24 h, respectively. In contrast, the nematode migration rates under dsgfp treatment at 6 h, 14 h, and 24 h were 32%, 65%, and 83%, respectively. The results showed that DdCα1D and DdCavβ could be important for the locomotion ability of D. destructor and that the simultaneous knockdown of both genes had a stronger effect on the locomotion ability (Figure 6).

2.6.2. Attraction Rate Assay of D. destructor after RNAi

To evaluate the effect of knockdown of the α1 and β subunits on the chemotaxis of D. destructor, an in vitro chemotactic bioassay was carried out with sweet potato blocks in water agar plates. As shown in Figure 7, the attraction rate of nematodes that were treated with dsCα1D and dsCα1D + dsCavβ decreased significantly to 7.13% and 3.88%, respectively, compared with 22.50% in the control. The silencing of the DdCα1D gene affected the chemotaxis of D. destructor toward the sweet potato, and this effect was more significant with co-silencing of the DdCavβ gene.

2.6.3. Stylet Thrusting of D. destructor after RNAi

The neurotransmitter serotonin is a conserved regulator of various behaviors in animals. The application of serotonin induces stylet thrusting in some plant-parasitic nematode species in the absence of a host [20,21,22]. To evaluate the effect of knockdown of the α1 and β subunits on the stylet thrusting of D. destructor, dsRNA-treated nematodes were soaked in serotonin (5 mM/L) for 20 min and the frequency of stylet thrusting per minute was observed under a microscope. As shown in Figure 8, the frequency of stylet thrusting per minute of the nematodes that were treated with dsCα1D and dsCα1D + dsCavβ was significantly lower than that of the nematodes that were treated with dsgfp. The dsCα1D + dsCavβ-treated nematodes had the lowest frequency of stylet thrusting per minute at 13.5, compared to 53.0 for dsgfp-treated nematodes. The results showed that silencing the DdCα1D gene had a strong impact on the stylet thrusting of D. destructor, and this effect was more obvious with co-silencing of the DdCavβ gene.

2.6.4. Detection of Secreted Proteins of D. destructor after RNAi

Based on pioneering work, resorcinol was used to induce esophageal gland secretion by D. destructor [23,24]. After soaking J3s in target dsRNAs and then incubating with 0.1% resorcinol for 16 h, the supernatant was used to detect the proteins that were secreted by the nematodes. For the worms that were treated with dsCα1A, the secreted protein content was significantly reduced to 19,740 μg/mL compared to the controls that were treated with dsgfp (20,246.67 μg/mL). The secreted protein content of dsCα1A + dsCavβ-treated nematodes was the lowest at 19,280 μg/mL (Figure 9). The results showed that silencing of the DdCα1A gene had an effect on the protein secretion by D. destructor, and co-silencing of the DdCavβ gene enhanced this effect.

2.6.5. Reproduction Rate of D. destructor after RNAi

As shown in Figure 10, after 25 days, the reproduction rate of nematodes in each treatment group was significantly lower than that of nematodes that were treated with dsgfp. Compared with the reproduction rate of nematodes that were treated with dsgfp (79.1%), the reproduction rate of nematodes that were treated with both dsCα1A and dsCavβ was the lowest at only 13.6%. The results showed that gene silencing affected the reproduction of D. destructor, and simultaneous silencing of the DdCα1A and DdCavβ genes enhanced this effect.

3. Discussion

VGCCs are macromolecular complexes that are composed of a pore-forming α1 subunit and related α2/δ, β, and γ auxiliary subunits and are embedded in the plasma membrane of most excitable cells [2]. Cavβ is a kind of cellular solute auxiliary subunit that consists of five distinct regions, namely, the NH2 terminus, an SH3 domain, a HOOK region, a GK domain, and the COOH terminus. The NH2, HOOK, and COOH regions are variable in length and amino acid sequence, and the SH3 and GK domains are highly conserved [9,12]. In this study, a type of Cavβ was identified for the first time in D. destructor. Through analysis of the amino acid sequence of DdCavβ, we found that it has a typical β-domain. In addition, sequence alignment showed that its SH3 domain and GK domain are highly conserved. The phylogenetic tree indicated that DdCavβ, CCB-1, and the Cavβs of L. loa and B. malayi were clustered together. According to these results, we preliminarily determined that the gene that was cloned in this study was the Cavβ of D. destructor.
C. elegans has two putative VGCC β subunits, namely, CCB-1 and CCB-2. We compared the amino acid sequences of CCB-1 and CCB-2 with the genome sequence of D. destructor in the WormBase database and identified a suspected β subtype [25]. After cloning, we found that the gene was highly homologous (68.48%) to CCB-1 and only 13.27% homologous to CCB-2. CCB-1, which is the main auxiliary subunit of the C. elegans VGCC, has contains two conserved subunits, namely, the SH3 and GK domains. CCB-1 can enhance the current of calcium channels, while CCB-2 does not affect this current and has no GK domain [12,14].
It has been reported that in the late growth stage of C. elegans and some plant-parasitic nematodes, the body wall muscle begins to undergo atrophy where the muscle mass and muscle function gradually decrease over time, and mobility gradually decreases [26,27]. In this study, DdCavβ was detected at different developmental stages and the expression level gradually increased from the egg stage to the J3 stage, followed by a significant decrease in the J4 stage. This may be related to atrophy of the body wall muscle in the later stage of stem nematode growth. In the early stage of nematode growth, the eggs and juveniles require more nutrients to complete the development and infection processes and body wall muscle mass and muscle function are continuously enhanced and improved. However, at the J4 stage, the body wall muscle of nematodes may begin to atrophy, and muscle mass and muscle function gradually decrease over time, leading to a decrease in DdCavβ gene expression.
Cavβ is expressed in muscle, neurons, and other tissues of vertebrates [12,28,29,30]. CCB-1 is widely expressed in most neuronal and muscle types in C. elegans [14]. Ye et al. showed that DdCα1D was expressed within body wall muscles and DdCα1A was expressed in the esophageal gland, vulva, and vas deferens of D. destructor [19]. In this study, we confirmed that DdCavβ was mainly expressed in the muscles. The results are similar to those regarding the expression of Cavβ in vertebrates and C. elegans, and as the main auxiliary subunit, it is also consistent with the results of Ye et al.
Cavα1 is the main functional subunit of VGCCs and plays an important role in the pharyngeal and body wall muscles of C. elegans [31,32,33]. Ye et al., by RNAi, confirmed that the DdCα1D gene plays a key role in the modulation of the cell wall muscle and normal locomotion in D. destructor and that the DdCα1A gene plays an important role in reproduction regulation in nematodes [19]. In this study, we silenced the HVA Cavα1 and β subunits, individually and in combination, to examine the effects on the locomotion behaviors and reproduction rate of the nematodes. The results showed that the locomotion ability of D. destructor was greatly affected by DdCα1D gene silencing, and that the reproduction rate of D. destructor was affected by silencing of the DdCα1A gene, especially when the DdCavβ gene and α1 subunit were co-silenced. It is suggested that the Cavα1 subunit plays an important role in the movement and reproduction of D. destructor. As the main auxiliary subunit, the Cavβ subunit plays an important auxiliary role with Cavα1. At the same time, these results further confirmed the functions of the L-type and non-L-type α1 subunits in the movement and reproduction of D. destructor.
Chemotaxis is movement in the direction of higher concentrations of semichemicals, such as plant chemical signals [34]. Chemical sensing by nematodes is an important part of their host-seeking behavior. In chemosensory neurons, odor concentration information is determined by the time integration of the increased intracellular calcium concentration of L-type VGCCs in a pair of olfactory neurons [35]. To date, only egl-19 (L-form) and unc-2 (non-L-form) have been shown to affect the chemical sensory signals of C. elegans [31,36,37]. Among the olfactory neurons of C. elegans, AWA neurons can guide nematodes to find potential food sources through unstable signals that are produced by bacteria and activate calcium-mediated action potentials that are initiated by egl-19 [38,39,40]. In our study, nematode chemotaxis was significantly inhibited by DdCα1D (L-type) gene silencing, especially with co-silencing of the DdCavβ gene. However, when the DdCα1A (non-L-type) gene was silenced, the chemotactic inhibition effect on the nematode was not as obvious as that which was observed upon silencing of the DdCα1D gene. The results showed that the L-type and non-L-type α1 subunits affected the chemotaxis of D. destructor toward sweet potato, and the L-type α1 had a strong effect on D. destructor movement. The results are similar to those for C. elegans. At the same time, the results show that the β subunit, as the main auxiliary subunit of voltage-gated ion channels, plays an important auxiliary role with the α1 subunit.
Serotonin is a neuroregulator that regulates feeding behavior in almost all phyla of the animal kingdom [41]. It also acts as a neuroregulator in C. elegans and is related to physiological functions such as pharyngeal pumping, egg laying, and locomotion [41,42,43]. It has been shown that, in C. elegans, serotonin activates pharyngeal pumping by the SER-7 serotonin receptor (a G protein-coupled receptor) in MC motor neurons, and SER-7 activates the downstream Gsα signaling pathway and then stimulates cholinergic transmission from MCs to the pharyngeal muscle [44]. To stimulate MC motor neurons, the CCA-1 (T-type) channel is involved in the initiation of action potentials, helping the membrane reach the threshold for activating EGL-19 (L-type) channels after the excitatory postsynaptic potential of the MC motoneurons and allowing reliable and rapid depolarization and contraction of the pharyngeal muscle [45]. In this study, soaking the treated nematodes with serotonin at 5 mM for 20 min significantly decreased the frequency of stylet thrusting of the nematodes after silencing of the DdCα1D gene, especially with co-silencing of the DdCavβ gene. The frequency of needle twitching of the nematodes did not change significantly after DdCα1A gene knockdown. In addition, D. destructor and C. elegans do not have voltage-gated Na+ channels, and the action potential depends on VGCCs [46,47]. In the pharyngeal muscles, the L-type VGCCs contribute to shaping the action polarization phase [48]. This suggests that the L-type α1 subunit plays a key role in the contraction of D. destructor stylet muscles.
In plant-parasitic nematodes, proteins that are synthesized in the esophageal gland (subventral and dorsal glands) and secreted through the stylet play important roles in the host-nematode relationship. By detection of the secreted nematode proteins after gene silencing, we found that the DdCα1A gene had a strong effect on protein secretion by D. destructor, especially when the DdCα1A and DdCavβ genes were co-silenced. We speculate that the non-L-type α1 subunit plays an important role in regulating protein secretion in the esophageal glands of D. destructor. At the same time, this result not only echoed the results of the expression localization study of the non-L-type α1 subunit in our previous work but also indicated the important auxiliary role of the β subunit with the α1 subunit.

4. Materials and Methods

4.1. Nematode Culture

D. destructor, isolated from diseased sweet potato and cultured on potato dextrose agar (PDA) plates that were inoculated with Fusarium semitectum, was employed in this study. After 30 days of incubation at 25 °C, D. destructor worms in mixed life stages were washed off the plates with distilled water and collected by the Baermann funnel technique [49]. The nematode eggs were collected via density gradient centrifugation, inoculated on PDA plates, and collected every 7 days to obtain nematodes at different stages [50].

4.2. RNA Extraction

Total RNA was extracted from D. destructor in mixed stages and different individual life stages by using TRIzol reagent (Invitrogen, Waltham, MA, USA) according to the manufacturer’s instructions. First-strand cDNA was synthesized from the total RNA by using the Super-Script Ⅲ First-Strand Synthesis Kit (Invitrogen) with OligodT primers according to the manufacturer’s protocol.

4.3. Cloning the Cavβ of D. destructor

On the basis of the protein sequence of Cavβ of C. elegans, the putative Cavβ gene sequence of D. destructor was identified in the WormBase database (https://parasite.wormbase.org/index.html, accessed on 24 November 2021), and primers (MF, MR) were designed for cloning the conserved sequence (Table 2). The PCR conditions used were as follows: 94 °C for 5 min; 30 cycles of 94 °C for 30 s, 48 °C for 30 s and 72 °C for 2 min; and a final extension at 72 °C for 7 min. The 3′-end cDNAs of the β subunit gene were obtained by rapid amplification of cDNA ends (RACE) PCR in conjunction with the use of a PrimeScript RTase Kit (TaKaRa, Japan). The 3′ RACE outer primer and 3′ RACE inner primer were provided with the kits. The 5′-end cDNAs of the β subunit gene were amplified with a primer upstream of the spliced leader (SL1) sequence, which is a sequence that is specific to nematode mRNA [51]. Gene-specific primers (3′ β-F1, 3′ β-F2, 5′ β-R1, and 5′ β-R2) (Table 2) were designed for 3′ and 5′ RACE amplification based on the conserved sequence of the β subunit gene; these sequences were obtained from previous sequencing results. After initial confirmation via agarose gel electrophoresis, all the resulting nested PCR products were cloned into the pMD™ 19-T vector (TaKaRa, Japan), and a single colony was sequenced (Sangon Biotech Co., Ltd., Shanghai, China).

4.4. Gene Expression Analysis by qPCR

qRT-PCR was used to assess the gene expression patterns of Cavβ. The total RNA was isolated from nematodes at different developmental stages (egg, J2, J3, and J4). The concentration of each sample was analyzed by a microspectrophotometer, and the samples were diluted to the same concentration with ddH2O. The RNA was then reverse-transcribed to cDNA and the relative gene expression of Cavβ was quantified via qRT-PCR. 18S rRNA was used as an endogenous control and the primers that were used are listed in Table 2. qRT-PCR was performed using SYBR Green mix according to the manufacturer’s instructions (Bio-Rad, Hercules, CA, USA). The qRT-PCR conditions were as follows: 95 °C for 2 min, followed by 39 cycles of 95 °C for 5 s and 59 °C for 40 s. qRT-PCR data from three biological and technical replicates were analyzed using Bio-Rad CFX Manager™ software. The relative transcript levels of each sample were calculated using the 2−∆∆Ct method [52]. The experiment included three biological repeats and three technical repetitions.

4.5. In Situ Hybridization of Cavβ in D. destructor

DIG-labeled sense and antisense probes were amplified via DIG RNA Labeling Mix (Roche Applied Science, Penzberg, Germany) together with BamHIF- and HindIIIR-specific primers (Table 2). In situ hybridization was performed as previously described [53], with some modifications. Briefly, D. destructor nematodes were fixed in 4% paraformaldehyde at 4 °C for 18 h, followed by an additional incubation at 22 °C for 4 h. Hybridization signals in the nematodes were detected with alkaline phosphatase-conjugated anti-DIG antibodies and substrate (Sangon Biotech Co., Ltd., Shanghai, China). Finally, the nematodes were placed on slides and then observed and imaged under a microscope (Carl Zeiss, Germany).

4.6. In Vitro RNA Interference Targeting VGCCs in D. destructor

Total D. destructor cDNA was used as a template for double-stranded RNA (dsRNA) synthesis using the MEGA Script RNAi Kit (Ambion, Austin, TX, USA) according to the manufacturer’s protocol. T7-labeled gene-specific primers (Table 2) were designed to amplify regions of the Cavβ and Cavα1 (L-type and non-L-type α1 subunit) transcripts. A nonendogenous control GFP dsRNA was synthesized using the GFP-F and GFP-R primers (Table 2). An RNAi assay for nematodes was conducted as described previously [54]. Approximately 10,000 J3-stage nematodes were soaked in 0.4 μg/μL target dsRNA in 200 μL of soaking buffer for 24 h in the dark on a slowly moving rotator at room temperature. In addition, J3s that were incubated in dsgfp and in soaking buffer (without dsRNA) served as the control. The nematodes that were treated with dsRNA were washed three times with ddH2O, RNA was extracted, and qPCR was performed using the methods that are described above to analyze the suppression of Cavβ and Cavα1 mRNA expression in D. destructor. FITC was used to trace the efficacy of the uptake of dsRNA by D. destructor. This experiment was independently repeated three times.

4.7. Analysis of Knockdown Phenotypes

After soaking for 24 h, the third-stage juveniles (J3s) were collected by centrifugation at 5000× g for 2 min and washed 3 times with ddH2O. The effects of Cavβ RNAi on the locomotion, chemotaxis, and reproduction rates of D. destructor were evaluated by the method of Shan Ye et al. [19]. In brief, worm motility was evaluated using a sand column functional migration assay. Approximately 100 dsRNA-treated J3s were added to the top of moistened sand columns and the columns were placed vertically in collection vials. The numbers of worms that were migrating through the columns and into the collection vials were counted at 6 h, 14 h, and 24 h. A total of 100 dsRNA-treated worms were inoculated on PDA with Fusarium semitectum. After 25 days, the nematodes were isolated and collected from the Petri dish and the total number of reproductions was counted to calculate the reproduction rate. The chemotaxis test was carried out on a 1% water agar plate in 90 mm glass dishes. Approximately 200 treated and control nematodes were added to the circular hole (1 cm diameter) in the center of the plate. Then, sweet potato slices with a diameter of 1 cm were placed 3.5 cm from the center, and the Petri dish was sealed with fresh-keeping film. After being placed in a dark incubator at 25 °C for 36 h, the nematodes within 2 cm of the sweet potato chips were isolated and counted and the attraction rate of sweet potatoes for nematodes was calculated. There were three biological repeats and three technical repeats in each detection experiment. The attraction rate was calculated as follows [55]:
Attraction rate (%) = Number of nematodes induced/Total number of nematodes input × 100.
The bioassay for nematode stylet thrusting was adapted from methodology that was described previously [56]. Approximately 50 J3s were incubated with 20 µL of serotonin solution (5 mM/L) in a 1.5 mL centrifuge tube for 20 min. The frequency of stylet thrusting per minute was observed under a microscope. A total of 10 nematodes were selected for each treatment. To detect the protein that was secreted by the esophageal gland of the nematode, approximately 2000 J3s were added to a 1.5 mL centrifuge tube containing 100 mL of 0.1% resorcinol neurotransmitter solution. After incubation at 25 °C for 16 h, the supernatant was used for determination of the secreted protein content with a Modified BCA Protein Assay Kit (Sangon Biotech Co., Ltd., Shanghai, China). Each process was performed for three biological repeats and three technical repeats.

4.8. Bioinformatic Analysis

To search for homologs, Cavβ sequences were compared using BLAST against the NCBI database (https://blast.ncbi.nlm.nih.gov/Blast.cgi, accessed on 12 October 2020). The Cavβ sequences were aligned using the DNAMAN software package (version 5.2.2, Lynnon Biosoft, San Ramon, CA, USA). Phylogenetic analyses of the Cavβ amino acid sequences and homologous sequences were obtained from the NCBI database (http://www.ncbi.nlm.nih.gov/blast/Blast.cgi, accessed on 12 October 2020) and prediction of its protein tertiary structure using the SWISS-MODEL online tool (https://www.swissmodel.expasy.org/, accessed on 24 November 2021).

5. Conclusions

Together, our results indicate that the Cavβ subunit plays an important auxiliary role associated with the Cavα1 subunit in D. destructor and further validates the role of the L-type and non-L-type α1 subunits in D. destructor. All the data provide additional clues for improving our understanding of the VGCCs of D. destructor.

Author Contributions

Conception of the work: Z.D. Collection of data: M.A., X.C., Z.Y. and J.Z. Analysis of data: Z.D., S.Y., M.A. and X.C. Writing of manuscript: M.A., S.Y. and Z.D. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Grant No. 31640063; No. 31872038; No. 32001879).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

We have no conflict of interest to declare.

References

  1. Arikkath, J.; Campbell, K.P. Auxiliary subunits: Essential components of the voltage-gated calcium channel complex. Curr. Opin. Neurobiol. 2003, 13, 298–307. [Google Scholar] [CrossRef]
  2. Karunasekara, Y.; Dulhunty, A.F.; Casarotto, M.G. The voltage-gated calcium-channel β subunit: More than just an accessory. Eur. Biophys. J. 2009, 39, 75–81. [Google Scholar] [CrossRef]
  3. Buraei, Z.; Yang, J. Structure and function of the β subunit of voltage-gated Ca2+ channels. Biochim. Biophys. Acta (BBA) Biomembr. 2012, 1828, 1530–1540. [Google Scholar] [CrossRef] [Green Version]
  4. Catterall, W.A. Structure and Regulation of Voltage-Gated Ca2+ Channels. Annu. Rev. Cell Dev. Biol. 2000, 16, 521–555. [Google Scholar] [CrossRef] [PubMed]
  5. Catterall, W.A. Voltage-Gated Calcium Channels. Cold Spring Harb. Perspect. Biol. 2011, 3, a003947. [Google Scholar] [CrossRef] [PubMed]
  6. Hidalgo, P.; Neely, A. Multiplicity of protein interactions and functions of the voltage-gated calcium channel β-subunit. Cell Calcium 2007, 42, 389–396. [Google Scholar] [CrossRef] [PubMed]
  7. Pragnell, M.; De Waard, M.; Mori, Y.; Tanabe, T.; Snutch, T.P.; Campbell, K. Calcium channel β-subunit binds to a conserved motif in the I–II cytoplasmic linker of the α1-subunit. Nature 1994, 368, 67–70. [Google Scholar] [CrossRef] [PubMed]
  8. Van Petegem, F.; Duderstadt, K.E.; Clark, K.A.; Wang, M.; Minor, D.L. Alanine-Scanning Mutagenesis Defines a Conserved Energetic Hotspot in the CaVα1 AID-CaVβ Interaction Site that Is Critical for Channel Modulation. Structure 2008, 16, 280–294. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. He, L.-L.; Zhang, Y.; Chen, Y.-H.; Yamada, Y.; Yang, J. Functional Modularity of the β-Subunit of Voltage-Gated Ca2+ Channels. Biophys. J. 2007, 93, 834–845. [Google Scholar] [CrossRef] [Green Version]
  10. Rima, M.; Daghsni, M.; Fajloun, Z.; M’Rad, R.; Brusés, J.L.; Ronjat, M.; De Waard, M. Protein partners of the calcium channel β subunit highlight new cellular functions. Biochem. J. 2016, 473, 1831–1844. [Google Scholar] [CrossRef]
  11. Dolphin, A.C. Calcium channel auxiliary α2δ and β subunits: Trafficking and one step beyond. Nat. Rev. Neurosci. 2012, 13, 542–555. [Google Scholar] [CrossRef] [PubMed]
  12. Buraei, Z.; Yang, J. The β Subunit of Voltage-Gated Ca2+ Channels. Physiol. Rev. 2010, 90, 1461–1506. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Bargmann, C.I. Neurobiology of the Caenorhabditis elegans Genome. Science 1998, 282, 2028–2033. [Google Scholar] [CrossRef]
  14. Lainé, V.; Frøkjær-Jensen, C.; Couchoux, H.; Jospin, M. The α1 Subunit EGL-19, the α2/δ Subunit UNC-36, and the β Subunit CCB-1 Underlie Voltage-dependent Calcium Currents in Caenorhabditis elegans Striated Muscle. J. Biol. Chem. 2011, 286, 36180–36187. [Google Scholar] [CrossRef] [Green Version]
  15. Sturhan, D.; Brzeski, M.W. Stem and bulb nematodes, Ditylenchus spp. In Manual of Agricultural Nematology, 1st ed.; Nickle, W.R., Ed.; CRC Press: Boca Raton, FL, USA, 1991; pp. 423–464. [Google Scholar]
  16. Faulkner, L.R.; Darling, H.M. Pathological histology, hosts, and culture of the potato rot nematode. Phytopathology 1961, 51, 778–786. [Google Scholar]
  17. Xu, Z.; Zhao, Y.-Q.; Yang, D.-J.; Sun, H.-J.; Zhang, C.-L.; Xie, Y.-P. Attractant and repellent effects of sweet potato root exudates on the potato rot nematode, Ditylenchus destructor. Nematology 2015, 17, 117–124. [Google Scholar] [CrossRef]
  18. Fan, W.; Wei, Z.; Zhang, M.; Ma, P.; Liu, G.; Zheng, J.; Guo, X.; Zhang, P. Resistance to Ditylenchus destructor Infection in Sweet Potato by the Expression of Small Interfering RNAs Targeting unc-15, a Movement-Related Gene. Phytopathology 2015, 105, 1458–1465. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Ye, S.; Zeng, R.; Zhou, J.; An, M.; Ding, Z. Molecular characterization of Ditylenchus destructor voltage-gated calcium channel α1 subunits and analysis of the effect of their knockdown on nematode activity. Biochimie 2020, 171–172, 91–102. [Google Scholar] [CrossRef]
  20. Masler, E. Responses of Heterodera glycines and Meloidogyne incognita to exogenously applied neuromodulators. J. Helminthol. 2007, 81, 421–427. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Han, Z.; Boas, S.; Schroeder, N.E. Serotonin Regulates the Feeding and Reproductive Behaviors of Pratylenchus penetrans. Phytopathology 2017, 107, 872–877. [Google Scholar] [CrossRef] [Green Version]
  22. Crisford, A.; Calahorro, F.; Ludlow, E.; Marvin, J.M.C.; Hibbard, J.K.; Lilley, C.J.; Kearn, J.; Keefe, F.; Johnson, P.; Harmer, R.; et al. Identification and characterisation of serotonin signalling in the potato cyst nematode Globodera pallida reveals new targets for crop protection. PLoS Pathog. 2020, 16, e1008884. [Google Scholar] [CrossRef]
  23. Jaubert, S.; Laffaire, J.B.; Piotte, C.; Abad, P.; Rosso, M.-N.; Ledger, T.N. Direct identification of stylet secreted proteins from root-knot nematodes by a proteomic approach. Mol. Biochem. Parasitol. 2002, 121, 205–211. [Google Scholar] [CrossRef]
  24. Bellafiore, S.; Shen, Z.; Rosso, M.-N.; Abad, P.; Shih, P.; Briggs, S.P. Direct Identification of the Meloidogyne incognita Secretome Reveals Proteins with Host Cell Reprogramming Potential. PLoS Pathog. 2008, 4, e1000192. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Zheng, J.; Peng, D.; Chen, L.; Liu, H.; Chen, F.; Xu, M.; Ju, S.; Ruan, L.; Sun, M. The Ditylenchus destructor genome provides new insights into the evolution of plant parasitic nematodes. Proc. R. Soc. B Boil. Sci. 2016, 283, 20160942. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Herndon, L.A.; Schmeissner, P.J.; Dudaronek, J.M.; Brown, P.A.; Listner, K.M.; Sakano, Y.; Paupard, M.C.; Hall, D.H.; Driscoll, M. Stochastic and genetic factors influence tissue-specific decline in ageing C. elegans. Nature 2002, 419, 808–814. [Google Scholar] [CrossRef]
  27. Han, Z.; Thapa, S.; Reuter-Carlson, U.; Reed, H.; Gates, M.; Lambert, K.N.; Schroeder, N.E. Immobility in the sedentary plant-parasitic nematode H. glycines is associated with remodeling of neuromuscular tissue. PLoS Pathog. 2018, 14, e1007198. [Google Scholar] [CrossRef]
  28. Ruth, P.; Röhrkasten, A.; Biel, M.; Bosse, E.; Regulla, S.; Meyer, H.E.; Flockerzi, V.; Hofmann, F. Primary Structure of the β Subunit of the DHP-Sensitive Calcium Channel from Skeletal Muscle. Science 1989, 245, 1115–1118. [Google Scholar] [CrossRef] [Green Version]
  29. Takahashi, M.; Seagar, M.J.; Jones, J.F.; Reber, B.F.; Catterall, W.A. Subunit structure of dihydropyridine-sensitive calcium channels from skeletal muscle. Proc. Natl. Acad. Sci. USA 1987, 84, 5478–5482. [Google Scholar] [CrossRef] [Green Version]
  30. Castellano, A.; Wei, X.; Birnbaumer, L.; Perez-Reyes, E. Cloning and expression of a neuronal calcium channel beta subunit. J. Biol. Chem. 1993, 268, 12359–12366. [Google Scholar] [CrossRef]
  31. Lee, R.Y.; Lobel, L.; Hengartner, M.; Horvitz, H.; Avery, L. Mutations in the alpha1 subunit of an L-type voltage-activated Ca2+ channel cause myotonia in Caenorhabditis elegans. EMBO J. 1997, 16, 6066–6076. [Google Scholar] [CrossRef] [Green Version]
  32. Jospin, M.; Jacquemond, V.; Mariol, M.-C.; Ségalat, L.; Allard, B. The L-type voltage-dependent Ca2+ channel EGL-19 controls body wall muscle function in Caenorhabditis elegans. J. Cell Biol. 2002, 159, 337–348. [Google Scholar] [CrossRef]
  33. Avery, L.; You, Y.-J. WormBook. WormBook 2012, 1–23. [Google Scholar] [CrossRef]
  34. Reynolds, A.M.; Dutta, T.K.; Curtis, R.H.C.; Powers, S.J.; Gaur, H.S.; Kerry, B.R. Chemotaxis can take plant-parasitic nematodes to the source of a chemo-attractant via the shortest possible routes. J. R. Soc. Interface 2010, 8, 568–577. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Tanimoto, Y.; Yamazoe-Umemoto, A.; Fujita, K.; Kawazoe, Y.; Miyanishi, Y.; Yamazaki, S.J.; Fei, X.; Busch, K.E.; Gengyo-Ando, K.; Nakai, J.; et al. Calcium dynamics regulating the timing of decision-making in C. elegans. eLife 2017, 6. [Google Scholar] [CrossRef] [PubMed]
  36. Schafer, W.; Kenyon, C.J. A calcium-channel homologue required for adaptation to dopamine and serotonin in Caenorhabditis elegans. Nature 1995, 375, 73–78. [Google Scholar] [CrossRef] [PubMed]
  37. Ferkey, D.M.; Sengupta, P.; L’Etoile, N.D. Chemosensory signal transduction in Caenorhabditis elegans. Genetics 2021, 217. [Google Scholar] [CrossRef]
  38. Choi, J.I.; Yoon, K.-H.; Kalichamy, S.S.; Yoon, S.-S.; Lee, J.I. A natural odor attraction between lactic acid bacteria and the nematode Caenorhabditis elegans. ISME J. 2015, 10, 558–567. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Worthy, S.E.; Haynes, L.; Chambers, M.; Bethune, D.; Kan, E.; Chung, K.; Ota, R.; Taylor, C.J.; Glater, E.E. Identification of attractive odorants released by preferred bacterial food found in the natural habitats of C. elegans. PLoS ONE 2018, 13, e0201158. [Google Scholar] [CrossRef]
  40. Liu, Q.; Kidd, P.B.; Dobosiewicz, M.; Bargmann, C.I. C. elegans AWA Olfactory Neurons Fire Calcium-Mediated All-or-None Action Potentials. Cell 2018, 175, 57–70.e17. [Google Scholar] [CrossRef] [Green Version]
  41. Horvitz, H.R.; Chalfie, M.; Trent, C.; Sulston, J.E.; Evans, P.D. Serotonin and Octopamine in the Nematode Caenorhabditis elegans. Science 1982, 216, 1012–1014. [Google Scholar] [CrossRef]
  42. Loer, C.M.; Kenyon, C.J. Serotonin-deficient mutants and male mating behavior in the nematode Caenorhabditis elegans. J. Neurosci. 1993, 13, 5407–5417. [Google Scholar] [CrossRef] [Green Version]
  43. Sawin, E.R.; Ranganathan, R.; Horvitz, H.C. elegans Locomotory Rate Is Modulated by the Environment through a Dopaminergic Pathway and by Experience through a Serotonergic Pathway. Neuron 2000, 26, 619–631. [Google Scholar] [CrossRef] [Green Version]
  44. Song, B.-M.; Avery, L. Serotonin Activates Overall Feeding by Activating Two Separate Neural Pathways in Caenorhabditis elegans. J. Neurosci. 2012, 32, 1920–1931. [Google Scholar] [CrossRef] [Green Version]
  45. Steger, K.A.; Shtonda, B.B.; Thacker, C.; Snutch, T.P.; Avery, L. The C. elegans T-type calcium channel CCA-1 boosts neuromuscular transmission. J. Exp. Biol. 2005, 208, 2191–2203. [Google Scholar] [CrossRef] [Green Version]
  46. Goodman, M.B.; Hall, D.H.; Avery, L.; Lockery, S.R. Active Currents Regulate Sensitivity and Dynamic Range in C. elegans Neurons. Neuron 1998, 20, 763–772. [Google Scholar] [CrossRef] [Green Version]
  47. Hobert, O. The neuronal genome of Caenorhabditis elegans. In WormBook: The Online Review of C. Elegans Biology; 2013; pp. 1–106. Available online: http://www.wormbook.org/index.html (accessed on 24 November 2021). [CrossRef]
  48. Shtonda, B.; Avery, L. CCA-1, EGL-19 and EXP-2 currents shape action potentials in theCaenorhabditis eleganspharynx. J. Exp. Biol. 2005, 208, 2177–2190. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Viglierchio, D.R.; Schmitt, R.V. On the methodology of nematode extraction from field samples: Baermann funnel modifica-tions. J. Nematol. 1983, 15, 438–444. [Google Scholar] [PubMed]
  50. Schaad, N.W.; Walker, J.T. The Use of Density-Gradient Centrifugation for the Purification of Eggs of Meloidogyne spp. J. Nematol. 1975, 7, 203–204. [Google Scholar]
  51. Blaxter, M.; Liu, L. Nematode spliced leaders–ubiquity, evolution and utility. Int. J. Parasitol. 1996, 26, 1025–1033. [Google Scholar] [PubMed]
  52. Livak, K.J.; Schmittgen, T.D. Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2−ΔΔCT Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  53. De Boer, J.M.; Yan, Y.; Smant, G.; Davis, E.L.; Baum, T.J. In-situ Hybridization to Messenger RNA in Heterodera glycines. J. Nematol. 1998, 30, 309–312. [Google Scholar]
  54. Urwin, P.E.; Lilley, C.J.; Atkinson, H.J. Ingestion of Double-Stranded RNA by Preparasitic Juvenile Cyst Nematodes Leads to RNA Interference. Mol. Plant-Microbe Interact. 2002, 15, 747–752. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. O’Halloran, D.M.; Burnell, A.M. An investigation of chemotaxis in the insect parasitic nematodeHeterorhabditis bacteriophora. Parasitology 2003, 127, 375–385. [Google Scholar] [CrossRef] [PubMed]
  56. Dutta, T.K.; Powers, S.J.; Gaur, H.S.; Birkett, M.; Curtis, R.H. Effect of small lipophilic molecules in tomato and rice root exudates on the behaviour of Meloidogyne incognita and M. graminicola. Nematology 2012, 14, 309–320. [Google Scholar] [CrossRef]
Ijms 23 00784 g001 550
Figure 1. DdCavβ sequence analysis. (A), The DdCavβ cDNA sequence and deduced amino acid sequence. The start codon (ATG) and the stop codon (TAA) are indicated in boxes. The fragment that is not translated before the start codon is the 3’ UTR, and the fragment that is not translated after the termination codon is the 5’ UTR. The underlined area represents SL1 at the front end of the 5’ UTR and the poly-A tail at the end of the 3’ UTR. The asterisk (*) indicates the stop codon (TAA). The light blue and brown colors indicate the NH2 terminus and the COOH terminus, respectively; yellow, the SH3 domain; purple, the HOOK region; and green, the GK domain. The residues that are involved in interactions with the AID are marked with the “•” symbol. (B), The predictive structure of DdCavβ.
Figure 1. DdCavβ sequence analysis. (A), The DdCavβ cDNA sequence and deduced amino acid sequence. The start codon (ATG) and the stop codon (TAA) are indicated in boxes. The fragment that is not translated before the start codon is the 3’ UTR, and the fragment that is not translated after the termination codon is the 5’ UTR. The underlined area represents SL1 at the front end of the 5’ UTR and the poly-A tail at the end of the 3’ UTR. The asterisk (*) indicates the stop codon (TAA). The light blue and brown colors indicate the NH2 terminus and the COOH terminus, respectively; yellow, the SH3 domain; purple, the HOOK region; and green, the GK domain. The residues that are involved in interactions with the AID are marked with the “•” symbol. (B), The predictive structure of DdCavβ.
Ijms 23 00784 g001
Ijms 23 00784 g002 550
Figure 2. Molecular phylogenetic analysis of DdCavβ via the neighbor-joining method. The species that were used in phylogenetic tree construction are summarized in Table 1. The phylogram was constructed based on the amino acid sequences of 16 Cavβ proteins using MEGA 7.0. The numbers below the branches indicate the bootstrap values.
Figure 2. Molecular phylogenetic analysis of DdCavβ via the neighbor-joining method. The species that were used in phylogenetic tree construction are summarized in Table 1. The phylogram was constructed based on the amino acid sequences of 16 Cavβ proteins using MEGA 7.0. The numbers below the branches indicate the bootstrap values.
Ijms 23 00784 g002
Ijms 23 00784 g003 550
Figure 3. Localization of the expressed DdCavβ gene in D. destructor via in situ hybridization. (A–D) In situ hybridization of sense probes for DdCavβ in D. destructor. (E–H) In situ hybridization of antisense probes for DdCavβ in D. destructor. st: stylet; mb: median bulb; e.g., esophageal gland; an: anus; vd: vas deferens; sp: spicules. Scale bar = 20 μm.
Figure 3. Localization of the expressed DdCavβ gene in D. destructor via in situ hybridization. (A–D) In situ hybridization of sense probes for DdCavβ in D. destructor. (E–H) In situ hybridization of antisense probes for DdCavβ in D. destructor. st: stylet; mb: median bulb; e.g., esophageal gland; an: anus; vd: vas deferens; sp: spicules. Scale bar = 20 μm.
Ijms 23 00784 g003
Ijms 23 00784 g004 550
Figure 4. The relative expression levels of DdCavβ at different developmental stages of D. destructor. The expression level in J2s was used as the standard value. The data are presented as the mean ± s.d of three biological replicates and three technical replicates (n = 9). Asterisks indicate a significant difference at the level of p < 0.01, as tested by Duncan’s new complex difference method.
Figure 4. The relative expression levels of DdCavβ at different developmental stages of D. destructor. The expression level in J2s was used as the standard value. The data are presented as the mean ± s.d of three biological replicates and three technical replicates (n = 9). Asterisks indicate a significant difference at the level of p < 0.01, as tested by Duncan’s new complex difference method.
Ijms 23 00784 g004
Ijms 23 00784 g005 550
Figure 5. The relative expression levels of DdCavβ, DdCα1D, and DdCα1A in D. destructor after soaking with specific dsRNA. (A) The relative expression levels of the DdCavβ, DdCα1D, and DdCα1A genes after dsCavβ treatment. (B) The relative expression levels of the DdCavβ, DdCα1D, and DdCα1A genes after dsCα1D treatment; (C) The relative expression levels of the DdCavβ, DdCα1D, and DdCα1A genes after dsCα1A treatment; (D) The relative expression levels of the DdCavβ, DdCα1D, and DdCα1A genes after dsCα1D + dsCavβ treatment; (E) The relative expression levels of the DdCavβ, DdCα1D, and DdCα1A genes after dsCα1A + dsCavβ treatment. Significant differences between the treatment and control are indicated with a line with asterisks (** p < 0.01; Student’s t test). “NS” indicates that there was no significant difference between the samples.
Figure 5. The relative expression levels of DdCavβ, DdCα1D, and DdCα1A in D. destructor after soaking with specific dsRNA. (A) The relative expression levels of the DdCavβ, DdCα1D, and DdCα1A genes after dsCavβ treatment. (B) The relative expression levels of the DdCavβ, DdCα1D, and DdCα1A genes after dsCα1D treatment; (C) The relative expression levels of the DdCavβ, DdCα1D, and DdCα1A genes after dsCα1A treatment; (D) The relative expression levels of the DdCavβ, DdCα1D, and DdCα1A genes after dsCα1D + dsCavβ treatment; (E) The relative expression levels of the DdCavβ, DdCα1D, and DdCα1A genes after dsCα1A + dsCavβ treatment. Significant differences between the treatment and control are indicated with a line with asterisks (** p < 0.01; Student’s t test). “NS” indicates that there was no significant difference between the samples.
Ijms 23 00784 g005
Ijms 23 00784 g006 550
Figure 6. Silencing of Cavβ or HVA Cavα1 and co-silencing of Cavβ and HVA Cavα1 affected the locomotory activity of D. destructor. Nematodes that were treated with gfp dsRNA were used as controls. Different lowercase letters indicate a significant difference at the level of p < 0.05, as tested by Duncan’s new complex difference method.
Figure 6. Silencing of Cavβ or HVA Cavα1 and co-silencing of Cavβ and HVA Cavα1 affected the locomotory activity of D. destructor. Nematodes that were treated with gfp dsRNA were used as controls. Different lowercase letters indicate a significant difference at the level of p < 0.05, as tested by Duncan’s new complex difference method.
Ijms 23 00784 g006
Ijms 23 00784 g007 550
Figure 7. Silencing of Cavβ or HVA Cavα1 and co-silencing of Cavβ and HVA Cavα1 caused defects in the chemotaxis of D. destructor toward sweet potato slices. Nematodes that were treated with gfp dsRNA were used as controls. Different lowercase letters indicate a significant difference at the level of p < 0.05, as tested by Duncan’s new complex difference method.
Figure 7. Silencing of Cavβ or HVA Cavα1 and co-silencing of Cavβ and HVA Cavα1 caused defects in the chemotaxis of D. destructor toward sweet potato slices. Nematodes that were treated with gfp dsRNA were used as controls. Different lowercase letters indicate a significant difference at the level of p < 0.05, as tested by Duncan’s new complex difference method.
Ijms 23 00784 g007
Ijms 23 00784 g008 550
Figure 8. Silencing of Cavβ or HVA Cavα1 and co-silencing of Cavβ and HVA Cavα1 affected the stylet thrusting of D. destructor. Nematodes that were treated with gfp dsRNA were used as controls. Different lowercase letters indicate a significant difference at the level of p < 0.05, as tested by Duncan’s new complex difference method.
Figure 8. Silencing of Cavβ or HVA Cavα1 and co-silencing of Cavβ and HVA Cavα1 affected the stylet thrusting of D. destructor. Nematodes that were treated with gfp dsRNA were used as controls. Different lowercase letters indicate a significant difference at the level of p < 0.05, as tested by Duncan’s new complex difference method.
Ijms 23 00784 g008
Ijms 23 00784 g009 550
Figure 9. Silencing of Cavβ or HVA Cavα1 and co-silencing of Cavβ and HVA Cavα1 affected the protein secretion by D. destructor. Nematodes that were treated with gfp dsRNA were used as controls. Different lowercase letters indicate a significant difference at the level of p < 0.05, as tested by Duncan’s new complex difference method.
Figure 9. Silencing of Cavβ or HVA Cavα1 and co-silencing of Cavβ and HVA Cavα1 affected the protein secretion by D. destructor. Nematodes that were treated with gfp dsRNA were used as controls. Different lowercase letters indicate a significant difference at the level of p < 0.05, as tested by Duncan’s new complex difference method.
Ijms 23 00784 g009
Ijms 23 00784 g010 550
Figure 10. Silencing of Cavβ or HVA Cavα1 and co-silencing of Cavβ and HVA Cavα1 affected the reproduction rate of D. destructor. Nematodes that were treated with gfp dsRNA were used as controls. Different lowercase letters indicate a significant difference at the level of p < 0.05, as tested by Duncan’s new complex difference method.
Figure 10. Silencing of Cavβ or HVA Cavα1 and co-silencing of Cavβ and HVA Cavα1 affected the reproduction rate of D. destructor. Nematodes that were treated with gfp dsRNA were used as controls. Different lowercase letters indicate a significant difference at the level of p < 0.05, as tested by Duncan’s new complex difference method.
Ijms 23 00784 g010
Table
Table 1. Cavβ information used for sequence alignment and phylogenic analysis.
Table 1. Cavβ information used for sequence alignment and phylogenic analysis.
SpeciesMolecule Name/Accession NumberIdentity (%)
Brugia malayiBMA-CCB-1/CRZ23597.167.5
Loa loabeta 2a/XP_020301925.167.2
Caenorhabditis elegansCCB-1/NP_491193.268.5
CCB-2/NP_001293380.113.3
Ditylenchus destructorMN867027
Drosophila melanogasterAAF21096.149.2
Heterololigo bleekeriBAB88219.150.2
Lymnaea stagnalisAAO83844.150.5
Homo sapiensbeta-1/NP_954856.144.8
beta-2/NP_000715.246.5
beta3/NP_000716.247.92
beta-4/NP_000717.248.7
Oryctolagus cuniculusbeta-subunit/AAA31180.144.9
CaB2b/CAA45576.143.9
CaB3/CAA45578.146.6
Rattus norvegicusbeta4/A4598248.6
Xenopus laevisAAA75519.146.4
Table
Table 2. Polymerase chain reaction (PCR) primers.
Table 2. Polymerase chain reaction (PCR) primers.
Name of PrimerSequence (5′-3′)Purpose
MFGATGGCAAAGTATAACCGAAGPrimers used for DdCavβ cloning
MRTTATCGGAAATCTTCTTCGTCGTAG
3′ RACE outer primerTACCGTCGTTCCACTAGTGATTT
3′ RACE inner primerCGCGGATCCTCCACTAGTGATTTCACTATAGG
3′ β-F1ACCTCACTTGCTCCTATTCATGTAT
3′ β-F2CCGTCCAGTCCACGATATGGCGAT
SL1GGTTTAATTACCCAAGTTTGAG
5′ β-R1TTTGGAGCTCTTAGAGGTTTG
5′ β-R2TTAGAGGTTTGCTGTCGAATG
18SFCTGATTAGCGATTCTTACGGAPrimers for real-time PCR analysis
18SRAGAAGCATGCCACCTTTGA
q-β-FAGCCGCTCAGCAATTGGACA
q-β-RTGAAAGACACTGCGCAGCCA
q-L-FGACCCGTTATTGTTGAGCCA
q-L-RACGTTCCTTCGAGATGAGA
q-NL-FTAGAAAACAGGCGAGACTTCC
q-NL-RCTCATCCGTTGTTCGATCCTC
BamHI FCGGGATCCGGAACGAGCAAACTCCAGGTCPrimers for ISH analysis
HindⅢ RCCCAAGCTTATGCTCACATGCGTCCTCAAG
dβ-FTAATACGACTCACTATAGGGAAGTTCCCCGAAAGTTCTCCAGPrimers used for synthesizing dsRNA
dβ-RTAATACGACTCACTATAGGGAGGTCCGCCCTTGTCATAATC
dL-FTAATACGACTCACTATAGGGAGGAAGATGACCTCTTGTTAG
dL-RTAATACGACTCACTATAGGGCCCAATATATGACCGTCTTTG
dNL-FTAATACGACTCACTATAGGGCGCAACACGTACCAAACTC
dNL-RTAATACGACTCACTATAGGGCTCATCTGAATCGCTAAGAGG
GFP-FTAATACGACTCACTATAGGGTACATCGCTCTTTCTTCACCG
GFP-RTAATACGACTCACTATAGGGACCAACAAGATGAAGAGCACC
The restriction enzyme sites and the T7 promoter sequences are underlined.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

An, M.; Chen, X.; Yang, Z.; Zhou, J.; Ye, S.; Ding, Z. Co-Silencing of the Voltage-Gated Calcium Channel β Subunit and High-Voltage Activated α1 Subunit by dsRNA Soaking Resulted in Enhanced Defects in Locomotion, Stylet Thrusting, Chemotaxis, Protein Secretion, and Reproduction in Ditylenchus destructor. Int. J. Mol. Sci. 2022, 23, 784. https://doi.org/10.3390/ijms23020784

AMA Style

An M, Chen X, Yang Z, Zhou J, Ye S, Ding Z. Co-Silencing of the Voltage-Gated Calcium Channel β Subunit and High-Voltage Activated α1 Subunit by dsRNA Soaking Resulted in Enhanced Defects in Locomotion, Stylet Thrusting, Chemotaxis, Protein Secretion, and Reproduction in Ditylenchus destructor. International Journal of Molecular Sciences. 2022; 23(2):784. https://doi.org/10.3390/ijms23020784

Chicago/Turabian Style

An, Mingwei, Xueling Chen, Zhuhong Yang, Jianyu Zhou, Shan Ye, and Zhong Ding. 2022. "Co-Silencing of the Voltage-Gated Calcium Channel β Subunit and High-Voltage Activated α1 Subunit by dsRNA Soaking Resulted in Enhanced Defects in Locomotion, Stylet Thrusting, Chemotaxis, Protein Secretion, and Reproduction in Ditylenchus destructor" International Journal of Molecular Sciences 23, no. 2: 784. https://doi.org/10.3390/ijms23020784

APA Style

An, M., Chen, X., Yang, Z., Zhou, J., Ye, S., & Ding, Z. (2022). Co-Silencing of the Voltage-Gated Calcium Channel β Subunit and High-Voltage Activated α1 Subunit by dsRNA Soaking Resulted in Enhanced Defects in Locomotion, Stylet Thrusting, Chemotaxis, Protein Secretion, and Reproduction in Ditylenchus destructor. International Journal of Molecular Sciences, 23(2), 784. https://doi.org/10.3390/ijms23020784

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

Article Metrics

Back to TopTop