- 1State Key Laboratory of Ecological Pest Control for Fujian and Taiwan Crops, Fujian Agriculture and Forestry University, Fuzhou, China
- 2College of Plant Protection, Fujian Agriculture and Forestry University, Fuzhou, China
- 3Oil Crops Research Institute, College of Agriculture/Center of Legume Crop Genetics and Systems Biology, Fujian Agriculture and Forestry University, Fuzhou, China
- 4College of Life Science, Fujian Agriculture and Forestry University, Fuzhou, China
Due to its high genetic diversity and broad host range, Ralstonia solanacearum, the causative phytopathogen of the bacterial wilt (BW) disease, is considered a “species complex”. The R. solanacearum strain FJ1003 belonged to phylotype I, and was isolated from the Fuzhou City in Fujian Province of China. The pathogen show host specificity and infects tobacco, especially in the tropical and subtropical regions. To elucidate the pathogenic mechanisms of FJ1003 infecting tobacco, a complete genome sequencing of FJ1003 using single-molecule real-time (SMRT) sequencing technology was performed. The full genome size of FJ1003 was 5.90 Mb (GC%, 67%), containing the chromosome (3.7 Mb), megaplasmid (2.0 Mb), and small plasmid (0.2 Mb). A total of 5133 coding genes (3446 and 1687 genes for chromosome and megaplasmid, respectively) were predicted. A comparative genomic analysis with other strains having the same and different hosts showed that the FJ1003 strain had 90 specific genes, possibly related to the host range of R. solanacearum. Horizontal gene transfer (HGT) was widespread in the genome. A type Ⅲ effector protein (Rs_T3E_Hyp14) was present on both the prophage and genetic island (GI), suggesting that this gene might have been acquired from other bacteria via HGT. The Rs_T3E_Hyp14 was proved to be a virulence factor in the pathogenic process of R. solanacearum through gene knockout strategy, which affects the pathogenicity and colonization ability of R. solanacearum in the host. Therefore, this study will improve our understanding of the virulence of R. solanacearum and provide a theoretical basis for tobacco disease resistance breeding.
Introduction
Bacterial wilt disease (BWD) caused by Ralstonia solanacearum is a devastating soil-borne disease (Cheng et al., 2021). Owing to its large genetic diversity, it is also known as the R. solanacearum species complex (RSSC). RSSC damages numerous economically important crops, such as tobacco (Xiao et al., 2018), peanut (Chen et al., 2021), tomato (Im et al., 2020), potato (Tan et al., 2019), and pepper (Xianyang et al., 2018), thus resulting in substantial economic losses. Upon infection, RSSC can invade the host plant’s vascular tissues of roots through the intercellular spaces, where it multiplies and secretes exopolysaccharides that block the vascular tissue, thus causing the death of the host plant (Xiaoqiang, 2018). RSSC also lurks in some asymptomatic hosts to prepare for its subsequent transmission (Genin and Denny, 2012). According to classical taxonomic classification, RSSC can be classified into four phylotypes (Ι, II, III, and IV) (Paudel et al., 2020). However, based on genomics and proteomics, RSSC was reclassified into three phylotypes (I, II, and III) (Prior et al., 2016). Particularly, phylotypes I and III belong to the first type, while phylotypes II and IV belong to the second and third types, respectively. This classification has been accepted by several scholars and used for further experimentations.
Tobacco (Nicotiana tabacum L.) is an economically important global crop and a model plant for transgenic research. BWD is the third most damaging factor affecting tobacco productivity after black root rot and black shank disease (Liu et al., 2007). Ralstonia solanacearum invades tobacco through root wounds or natural cracks and then continues to infect xylem tissue, resulting in blockage of the host vascular bundle system and finally wilting and death (ZHOU et al., 2012). Tobacco bacterial wilt is widely occurring in the main tobacco-growing areas of China. Its incidence rate is about 80%, and in extreme cases, the crop is out of production, resulting in great economic losses to farmers (Wang, 2018).
The virulence factors of R. solanacearum responsible for promoting BWD mainly include exopolysaccharides, plant cell wall degrading enzymes, type III secretion system (T3SS), and type IV secretion system (T4SS) (Genin & Denny, 2012). Among them, the type III effector proteins (T3Es) secreted by T3SS are vital in the pathogenesis of R. solanacearum. Among bacteria, R. solanacearum has the most T3Es; with 60–75 of them existing in a single R. solanacearum genome (Genin and Denny, 2012; Peeters et al., 2013). Due to the large number of T3Es, R. solanacearum has a larger host range than other bacteria (Coll and Valls, 2013). Nonetheless, most T3Es are virulent, with only a few avirulent ones (Landry et al., 2020).
The whole-genome study of R. solanacearum gave novel insights into this complex species for further investigations. Previously, since the whole genome of the strain GMI1000 was sequenced first (Salanoubat et al., 2002), it is also known as the model strain. To date, 291 strains have complete genome assemblies and annotations in the NCBI database (https://www.ncbi.nlm.nih.gov/genome/browse/#!/prokaryotes/490/). Most of them consist of only one chromosome, with only a few having an additional small plasmid. Tobacco was the host of both Rs10 and CQPS-1 among the sequenced strains (https://www.ncbi.nlm.nih.gov/genome/browse/#!/prokaryotes/490/).
Therefore, this study reports the genome sequence of the FJ1003 strain of R. solanacearum isolated from tobacco in Fujian Province of China and analyzes its evolution, comparative genomics, and specific genes. The functions of Rs_T3E_Hyp14 were also investigated. Our findings provide a theoretical basis for the control and pathogenic mechanism of tobacco BWD and enrich the genomic information of R. solanacearum.
Materials and methods
Strain preparation and extraction of genomic DNA
The highly pathogenic R. solanacearum strain “FJ1003” (Zhang et al., 2017; Zhang et al., 2019) was provided by Professor Liu Bo from the Fujian Academy of Agricultural Sciences, Fuzhou, China (isolated from tobacco growing areas in Fujian Province). The strain was cultured on the bacteria peptone glucose (BG) solid medium (peptone 1.0%, yeast extract 0.1%, casamino acid 0.1%, glucose 1.0%, agar 1.4% and pH 7.4), and 50 ml BG liquid medium (peptone 1.0%, yeast extract 0.1%, casamino acid 0.1%, glucose 1.0%, pH 7.4) was used for single clone culturing at 28°C, 200 rpm overnight. The bacterial genome extraction kit (TIANGEN, Beijing, China) was used to extract the genomic DNA from the bacterial culture solution of R. solanacearum for subsequent sequencing experiments.
Genome sequencing and assembly
The constructed bacterial genome library was sequenced using the single-molecule real-time (SMRT) sequencing technology by Biomarker Biosciences Co., Ltd. (Beijing, China) (Eid et al., 2009; Faino et al., 2015). Low-quality filtering was performed on the original reads obtained by sequencing; reads with a read length less than 100 bp and average read quality less than 0.7 were removed. The longest reads were selected as seed reads, and the remaining reads were compared to seed reads. Seed reads are converted into highly accurate pre-assembled reads for genome assembly. The complete genome was obtained by assembling the filtered subreads using the MHAP software with the default parameter (Chin et al., 2013; Berlin et al., 2015). The splicing and assembly are carried out in an overlapping manner, and the result is in the form of a scaffold sequence. Mainly through the overlap-layout-consensus method, the rectified high-quality pre-assembled reads sequences were assembled de novo. The quality of the assembly was verified using a probabilistic algorithm to determine the final gene sequence, and the ends of the assembled sequence were trimmed to circularize the genome.
Genomic component analysis
We used the Repeat Masker (Tarailo Graovac and Chen, 2009) software to screen the repeated sequences in the bacterial genome. All coding genes and non-coding RNAs were predicted using the Prodigal software (Hyatt et al., 2010). GeneWise (Birney et al., 2004) was used to predict the presence of pseudogenes in the genome. The prophages and genomic islands were predicted using the PhiSpy (Akhter et al., 2012) and IslandPath-DIOMB (Bertelli and Brinkman, 2018) software, respectively. Furthermore, the PILLERCR software (Edgar, 2007) was used to predict the Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) in the genome.
Functional annotation
The predicted gene sequences were BLASTed using functional annotation databases like COG (Tatusov et al., 2000), KEGG (Kanehisa, 2004), GO (Ashburner et al., 2000), Swiss-Prot (Consortium, 2018), NR (Deng et al., 2006), and PHI (http://www.phi-base.org/index.jsp). The functions of the genes were analyzed using the COG, KEGG metabolic pathways, GO, and PHI enrichment analysis.
Gene family and unique gene analysis
OrthoMCL (Li et al., 2003) software was used to cluster the protein sequences of FJ1003, and the reference genome (GMI1000, CMR15, UW551, PSI07, CQPS_1, RS10) was used to discover the unique gene families of the strain. Based on sequence similarity, OrthoMCL can classify a set of proteins, such as genome-wide proteins, into ortholog groups. A total BLAST database was built using the protein sequences. All these sequences were run using BLASTP against the local database to get the results with an e-value < 1e−5. The percentage match length of the alignments was calculated using a 50% cut-off value. All ortholog, in-paralog, and co-ortholog pairs, and their normalized weight values, were entered into the MCL program for clustering.
Analysis of evolutionary relationship among species
The protein sequences of the endoglucanase genes from FJ1003 and the reference genome were obtained using the annotation information. All the protein sequences used to build the RS_T3E_Hyp14 phylogenetic tree were downloaded from the T3Es (https://iant.toulouse.inra.fr/bacteria/annotation/site/prj/T3Ev3/) and NCBI (https://www.ncbi.nlm.nih.gov/) database. The evolutionary history was inferred using the maximum likelihood method based on the JTT (Jones-Taylor-Thornton) matrix-based model (Jones et al., 1992). The tree with the highest log-likelihood was shown. The percentage of tree in which the related taxa were clustered together was shown next to the branches. Initial tree(s) for the heuristic search were automatically obtained by applying the neighbor-joining (NJ) and BioNJ algorithms to a matrix of pairwise distances estimated using the JTT model and then selecting the topology with a superior log-likelihood value. The tree was drawn to scale, with the branch lengths measured in the number of substitutions per site. All positions containing gaps and missing data were eliminated. Finally, evolutionary analysis was performed in MEGA7 (Kumar et al., 2016).
Collinearity analysis
Collinearity analysis was performed using the Multiple Collinearity Scan Toolkit (MCScanX) (Wang et al., 2012) to find the collinear genes between FJ1003 and the reference genome. Using the FJ1003 genomic protein sequences as a database and the reference genomic protein sequences as queries, the FJ1003 and other reference genomic protein sequences were compared by BLASTP. Finally, the file obtained after BLASTP was used as an input in MCScanX with the default parameters to visualize the results.
Construction of mutant and complementary strains of Ralstonia solanacearum
The construction of mutant and complementary strains was carried out following the method of GUI et al. (2021), with slight modifications. Using the genomic DNA of R. solanacearum as a template, the upstream 440 bp and downstream 469 bp fragments of Rs_T3E_Hyp14 were amplified by specific primers (F1:5′-TCAACCGCAGCTACATCAGCTGCATG-3′,R1:5′-GTTACCGGACAGCGAAGAGCCAGCATTTGTGCAACTCCGTTACTGAGATCG-3′,F2:5′-CGATCTCAGTAACGGAGTTGCACAAATGCTGGCTCTTCGCTGTCCGGTAAC-3′,R2:5′-AAGAGCCGCTCTTTGCACTCCGATCA-3′). The upstream and downstream fragments were linked by overlapping PCR, digested by restriction enzymes XbaI and SphI (Takara, Shiga, Japan), and then ligated into the PK18 vector (PK18-U-D) by T4 DNA ligase enzyme (Takara, Shiga, Japan) and transformed into R. Solanacearum FJ1003, by homologous recombination. The Rs_T3E_Hyp14 in R. solanacearum was removed, and the Rs_T3E_Hyp14 mutant (ΔRs_T3E_Hyp14) was obtained that was confirmed by specific primers. For the complementary strains of Rs_T3E_Hyp14, the fragment U-Rs_T3E_Hyp14-D was amplified from the genomic DNA of R. solanacearum by the F1 and R2 primers. The fragment “U-Rs_T3E_Hyp14-D” was ligated into the PK18 vector (PK18-U-Rs_T3E_Hyp14-D) by T4 DNA ligase enzyme and transformed into the mutant strain. The complementary strain (CRs_T3E_Hyp 14ΔRs_T3E_Hyp14) was obtained by homologous recombination.
Inoculation of Ralstonia Solanacearum and investigation of disease index
The R. solanacearum was cultured in BG medium until OD600 = 0.5 and then resuspended with sterile water. A susceptible tobacco variety, Honghua Dajinyuan (HD), was grown in plastic pots filled with nutrient soil to four-leaf stage. A sterilized scalpel was taken and inserted into the pot, causing mechanical damage to the roots. Then, 5 ml of bacterial solution was added to each pot and continued to culture in a greenhouse at 28°C. According to the classification and investigation methods of tobacco diseases and insect pests in China (GB/T23222—2008) (China, 2008), the disease index of infected tobacco plants was investigated at different time points. Three days after inoculation with R. solanacearum, 1 g of stem samples were ground in 1 ml sterile water, applied on BG medium after gradient dilution, and the number of colonies was calculated.
Results
Genome sequencing, assembly, and annotations
Using the SMRT sequencing technique, we sequenced the R. solanacearum FJ1003 strain isolated from tobacco at an average depth of 150X-250X, and obtained 1.38 Gb of high-quality sequencing data. We obtained three scaffolds, one bacterial chromosome, and two bacterial plasmids with a total length of 5.90 MB using the MHAP (Chin et al., 2013; Berlin et al., 2015) software (Table 1). Upon comparing the sequencing data to the assembled genome, the scaffolds were close to QV50 (Supplementary Figure S1), suggesting that we assembled a complete genome sequence.
The (G + C) content in prokaryotes is generally higher than in eukaryotes. The (G + C) content of the FJ1003 sequenced genome was ∼67%. The statistical results of the genomic characteristics are presented in Table 1. In short, the FJ1003 genome contained 5133 coding genes, of which 3446, 1564, and 123 genes were located on the chromosome, the megaplasmid, and the small plasmid, respectively. The FJ1003 genome also had 447 non-coding RNAs (4 miRNAs, 35 tRNAs, and 408 rRNAs) and 21 pseudogenes. CRISPR is an important component of the immune system in prokaryotes that helps resist foreign plasmids and phage sequences. It holds the potential to recognize and silence the invading functional elements (Zhong et al., 2021). The predicted CRISPR sequences in the sequenced species are described in Supplementary Table S1. At the same time, based on the sequence and assembly results of the genome, we developed a circular genome map (Figure 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/966092/fgene-13-966092-HTML/image_m/fgene-13-966092-g001.jpg)
FIGURE 1. Circular map showing the genome of the R. solanacearum strain FJ1003. The outermost circle indicates the marker of genome size, with each scale of 0.1 Mb; the second and third circles indicate the genes on the positive and negative strands of the genome, respectively, and different colors represent the different functional categories of COG; the fourth circle represents repetitive sequences, and the fifth circle is tRNAs. The innermost layer indicates the GC content.
Furthermore, we performed the functional annotation analysis of identified genes in the FJ1003 genome using the GO, NR, and COG databases. Among them, 3913 genes were annotated by the COG database, and 3855 genes by the GO database (Figure 2). In the NR database, 5064 genes were annotated, of which 4827 (94.81%) genes were highly conserved in R. solanacearum (Figure 3A). Similarly, the predicted results were used to BLAST against the PHI database to obtain the enrichment results (Figure 3B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/966092/fgene-13-966092-HTML/image_m/fgene-13-966092-g002.jpg)
FIGURE 2. COG and GO functional categories of protein-coding genes in the R. solanacearum FJ1003 genome. (A) Statistical results of COG classification; the abscissa indicates the content of COG classification; the ordinate is the number of genes. (B) Clustering results of GO annotation; the abscissa indicates the content of each GO classification; the left side of the ordinate indicates the percentage of the number of genes, and the right side indicates the number of genes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/966092/fgene-13-966092-HTML/image_m/fgene-13-966092-g003.jpg)
FIGURE 3. Enrichment analysis results based on the NR and PHI databases. (A) Statistical results of species distribution by the NR database. Different colors represent different homologous species. (B) The gene enrichment annotation to the PHI database. Different colors represent different types of disease-related genes.
Comparative analysis of the virulence factors in FJ1003
R. solanacearum has many virulence factors, among which the T3Es secreted by the T3SS, are considered the primary virulence factors (Landry et al., 2020). Based on the R. solanacearum T3Es (Peeters et al., 2013) and the NCBI database, we predicted 76 T3Es and compared them with T3Es from four other isolates (Table 2). Although RipBP and RS_T3E_Hyp6 were absent from the six other reference genomes, they were mainly found in the NCBI database searches of other R. solanacearum genomes. Notably, we found that five of the FJ1003 strain genes in the T3Es were absent from the model strain GMI1000. Although both CQPS-1 and FJ1003 infect tobacco, 19 T3Es were present in FJ1003 and were absent from CQPS-1 (Table 2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/966092/fgene-13-966092-HTML/image_m/fgene-13-966092-t002.jpg)
TABLE 2. Summary of type III effector proteins genes (Coverage%/Identity%) in FJ1003 and other strains.
The R. solanacearum T3SS family is considered to contain more than 20 structural genes. Therefore, we compared the FJ1003 T3SS with GM1000 T3SS and other virulence factors such as exopolysaccharides, cell wall degrading enzymes, and swimming ability-related genes (Supplementary Table S2). These virulence factors were conserved within the bacterial genome. The pehB, FliO, FliN, FliS, FliT, and FliK genes were 99% identical, while the remaining genes were 100% identical (Supplementary Table S2).
Evolutionary relationship among species
The detailed genome information of the genome sequences of the different phylotypes of R. solanacearum obtained from the NCBI database is presented in Supplementary Table S3. The phylogenetic tree (Figure 4) was drawn based on the sequence similarity among the endoglucanase genes from different genomes (Villa et al., 2005). Phylogenetic analysis revealed that the FJ1003 strain belonged to the same branch as the other phylotype I strain. Phylotypes I and III were closely related, while II and IV belonged to two independent branches, indicating that the “proteomics-genomics” method of classifying R. solanacearum into three classes is better suited for its classification. Despite belonging to the same phylotype, the genetic relationships between strains may vary depending on their host or geographical origin.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/966092/fgene-13-966092-HTML/image_m/fgene-13-966092-g004.jpg)
FIGURE 4. Molecular phylogenetic analysis between FJ1003 and other strains. This figure is mainly based on the sequence similarity of endoglucanase in different genomes. The percentage of tree in which the associated taxa were clustered together, is shown next to the branches. The tree is drawn to scale, with the corresponding branch lengths indicating the number of substitutions per site. The pink font on the far right is the name of the different species of R. solanacearum.
Some genes are specific to the FJ1003 genome
To find which genes are specific to the FJ1003 genome, we compared it with those of four other different strains from different phylotypes. The results showed that 73 genes were specific to the FJ1003 strain (Figure 5A; Supplementary Tables 4, S5-1). They included unknown functional proteins, transposable enzymes, twitching mobility, and plasmid replication initiation proteins. At the same time, we also analyzed the FJ1003 protein sequences with two other strains (CQPS-1 and RS10) with the same phylotype and host and found that 90 genes were unique to FJ1003 (Figure 5B; Supplementary Table S5-2). These unique genes included hypothetical proteins, hemagglutinin-like secreted proteins, membrane proteins, and integrases. Therefore, our results suggest that the genome of R. solanacearum is relatively complex and has high genetic diversity.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/966092/fgene-13-966092-HTML/image_m/fgene-13-966092-g005.jpg)
FIGURE 5. Venn diagrams showing the number of specific and shared genes among the R. solanacearum strain FJ1003 and other bacterial species. (A) Venn diagrams of the sequenced strains and reference strains from different phylotypes. (B) The hosts of the three strains are tobacco. Note: the overlapping regions in the figure indicate the number of gene families shared among the different species, while the regions that do not overlap with other species indicate the number of gene families specific to that species.
Genomic collinearity analysis
The analysis of collinear genes with other reference genomes (GMI1000, CMR15, PSI07, UW551, and CQPS_1, RS10) showed that there was similarity among the whole genomes of the different strains of R. solanacearum. However, there were some cases of local non-collinearity, indicating that there were a few inverted and ectopic genome rearrangements (Figure 6). We also found some blank regions in the genome, indicating that there was dissimilarity among the sequences of the various strains, possibly with the specific genes of a particular strain.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/966092/fgene-13-966092-HTML/image_m/fgene-13-966092-g006.jpg)
FIGURE 6. The collinear relationship between the FJ1003 strain and other strains. The genome location coordinates of the reference near-source species are measured on the left, while the assembled genome location coordinates are on the right. The lines in the figure represent the collinear regions between the two species. Lines of different colors represent different collinear areas between different chromosomes. The color of the lines is automatically assigned by the software.
Horizontal gene transfer is widespread in the genome
In terms of microbial evolution, horizontal gene transfer (HGT) is referred to as the transfer of genetic material between non-parental organisms, either interspecifically or intraspecifically, making it a fairly widespread phenomenon with a vital role in the evolutionary processes (Villa et al., 2005; Shutian, 2020). A genetic island (GI) acts as an integrated moving element on the host population gennome by facilitating gene exchanges within or between species. Using HGT between different hosts, the acceptor strains obtain large DNA fragments, increasing their versatility and adaptability. During the process of bacterial HGT, mild bacteriophages transfer the donor genes to the recipient, and this transferred genetic fragment which is integrated into the recipient, is called a prophage. The FJ1003 genome contains 23 GIs and five Prophages (Supplementary Table S6). RipP2, RipP3, RipBP, and RS_T3E_Hyp14 are present in the prophages, whereas Rs_T3E_Hyp14 and RS_T3E_Hyp6 are located on the GI (Supplementary Table S6).
RS-T3E-Hyp14 is involved in the pathogenicity of Ralstonia solanacearum
Notably, Rs_T3E_Hyp14 is present on prophage and GI both, proposing that this gene might have been acquired from other bacteria via HGT. We predict the tertiary structures of Rs_T3E_Hyp14 using SWISS-MODEL (https://swissmodel.expasy.org/interactive/PW9vQ4/models/) (Figure 7A) and found it to be a typical feature of T3E AvrRpt2 (a C70 family cysteine protease) from Erwinia amylovora. At the same time, through NCBI and T3E database comparison, we found that 14 strains of Ralstonia solanacearum had Rs_T3E_Hyp14 gene, and by blastp analysis, the Rs_T3E_Hyp14 gene of five strains was 100% similar to that of FJ1003 (Supplementary Table S7). We constructed an evolutionary tree of Rs_T3E_Hyp14 of these 14 strains using the ML method and found that the genetic distances of Rs_T3E_Hyp14 in most of the phylotype I strains were similar (Supplementary Table S7, Figure 7B). There is no report on the function of RS-T3E-Hyp14. The pathogenicity of R. solanacearum on tobacco is significantly reduced after we knockout RS-T3E-Hyp14, and the colonization of R. solanacearum on a host is also affected (Figures 7C–E). The phenotype of complementary strains is consistent with that of wild strains (Figures 7C–E). These results suggest that RS-T3E-Hyp14 is a virulence factor in the pathogenic process of R. Solanacearum FJ1003, which affects its pathogenicity.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/966092/fgene-13-966092-HTML/image_m/fgene-13-966092-g007.jpg)
FIGURE 7. Functional identification of RS-T3E-Hyp14 (A) Tertiary structures of Rs_T3E_Hyp14 in FJ1003. (B) Molecular phylogenetic analysis of predicted RS-T3E-Hyp14 effectors. The phylogenetic tree was constructed by the maximum likelihood method. (C) Symptoms of tobacco after 8 days of inoculation with Ralstonia solanacearum. (D) Disease index of Ralstonia solanacearum at different time after inoculation (E) Colonization of Ralstonia solanacearum in tobacco stems after inoculation for 3 days. Values are means ± SE for three replicates. (WT: wild type strains FJ1003, M: mutant strains ΔRs_T3E_Hyp14, C: complementary strains CRs_T3E_Hyp 14ΔRs_T3E_Hyp14).
Discussion
Bacterial wilt caused by R. solanacearum is a common bacterial disease affecting tobacco production. It is rapidly spreading and very destructive during all growing stages of tobacco (Lin, 2018). R. solanacearum demonstrates high genetic diversity and host specificity. FJ1003 is a strain isolated from tobacco and has not been reported to affect other host species. To date, 291 strains of R. solanacearum have been sequenced and assembled, among which both the RS10 and CQPS-1 strains were isolated from tobacco. Upon comparing the protein sequences of FJ1003 and four other reference strains, 73 genes were found to be FJ1003-specific, while the other two tobacco isolates (RS10 and CQPS-1) all had their specific genes. These results may further prove that R. solanacearum is genetically diverse. Moreover, subsequent whole-genome sequencing of the R. solanacearum isolated from tobacco from different geographical regions could shed light on this point of view for future investigations.
Using the syringe-like T3SS, R. solanacearum can inject multiple toxic effector proteins into the host cells, thereby disturbing the host immune system and finally causing host disease. Effector proteins have a wide variety of functions, including the effector proteins (RipM, RipS1, RipQ, RipG3, RipD, RipAD, and RipAU) that inhibit the plant immune response triggered by the bacterial flg22 (Landry et al., 2020). Upon injection of the bacterial effector protein RipTPS into the plant cell, trehalose-6-phosphate is synthesized in planta to provide nutrition for the reproduction of R. solanacearum in the host (Poueymiro et al., 2014). RipN can adversely alter the NADH/NAD + ratio and inhibit the immune response triggered by pathogen-associated molecular patterns (PAMPs) (Sun et al., 2019). RipI can modulate plant metabolism to enhance γ-aminobutyric acid (GABA) production (a signalling molecule) in plant cells, which R. solanacearum uses as a nutrition source during its reproductive process in plants, thus finally causing plant disease (Xian et al., 2020). Few effector proteins, like RipP1 and RipP2, function as avirulence proteins (Landry et al., 2020). Thus, these results indicate that FJ1003 contains 76 T3Es, with the function and pathogenesis of most effector proteins still being unknown, and further research is needed. The functional annotation by the PHI database indicates that some effector proteins related to the host-pathogen interaction may play a key role in the pathogenesis of R. solanacearum.
HGT, the exchange of genetic material between individuals, is an important factor behind microbial evolution and diversity (Syvanen, 1985). We detected five prophages (a total of 316 genes) and twenty-three GIs (a total of 252 genes) in the R. Solanacearum FJ1003 genome. Interestingly, RS-T3E-Hyp14 is present in both prophages and genomic islands, suggesting that this gene is associated with HGT. The function of RS-T3E-Hyp14 has not been reported. We identified the function of RS-T3E-Hyp14 and found that it could significantly reduce the pathogenicity of Ralstonia solanacearum and play an important role in the pathogenic process of R. solanacearum. Further studies on the pathogenesis of RS-T3E-Hyp14 and its interaction with tobacco are needed.
Most prokaryote genomes contain the CRISPR-Cas system, mainly containing CRISPR sequences (Louwen et al., 2014). This system is developed throughout the long-term struggle of prokaryotes against harmful foreign substances. It is a mechanism for protecting prokaryotes from harmful foreign agents. The current study predicted three copies of 100-bp CRISPR sequences in the FJ1003 genome. The widespread presence of HGT and CRISPR sequences in R. solanacearum may be one of the reasons for its wide host range, biological diversity, and strong pathogenicity.
The PHI database contains the pathogenic genes, virulence genes and effector protein genes of pathogens that have been verified by experiments or reported in the literature (http://www.phi-base.org/index.jsp). Inactivation or reduced expression of most of these genes can possibly result in the reduction or even the complete loss of the pathogenic ability of the pathogen in the corresponding host. The annotation of 869 genes in the FJ1003 genome by the PHI database falls under seven different modules (reduced virulence, unaffected pathogenicity, loss of pathogenicity, mixed outcome, lethal, increased virulence, effector, and chemistry target). These genes may play different roles in the host-pathogen interaction (ZANG et al., 2021). Therefore, further research is needed to provide an important theoretical basis for the screening and functional verification of pathogenic proteins.
Conclusion
In this study, we sequenced the genome of a plant pathogen, Ralstonia Solanacearum strain FJ1003, isolated from Fujian province, China. Analysis showed that the genome consists of a chromosome, a megaplasmid, and a small plasmid. The evolutionary analysis provided that FJ1003 belonged to phylotype I and this strain contained 76 T3Es. Horizontal gene transfer was widespread in this strain, and the RS-T3E-Hyp14 existed in both prophages and genomic islands. This study is the first report on the genetic functions of RS-T3E-Hyp14. The knockout mutants of RS-T3E-Hyp14 significantly reduced the pathogenicity of R. solanacearum and weakened its colonization in the host. Results provide valuable information for the control of bacterial wilt in tobacco plants.
Data availability statement
The genome sequences have been deposited in the NCBI SAR database under the accession number PRJNA780084, available at https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA780084.
Author contributions
The idea was convinced by WZ and KC was responsible for writing the manuscript and analyzing the data. YZ, LW, HL, HY, and MG took part in the data analysis and visualization. The experimental part of the manuscript was completed by ML, TL, JW, AR, CZ, and HD. The figures were drawn by YS, QY, HZ, and WZ lead the project by providing the technical guidlines for the experimental design, experimental procedure and data analysis, as well as the revision of the first draft of the manuscript. All authors have read and agreed to the published version of the manuscript.
Funding
This work was funded by the National Natural Science Foundation of China (grant numbers: U1705233, 32072103, and 31701463) and the Science and Technology Foundation of Fujian Province of China (2018N0004).
Acknowledgments
We are very grateful to Professor Liu Bo for providing us with Ralstonia solanacearum strain FJ1003, and also to the support provided by the National Natural Science Foundation of China (grant numbers: U1705233, 32072103, and 31701463). We would like to express our sincere gratitude to the Science and Technology Foundation of Fujian Province of China (2018N0004) for supporting this research. We would also like to thank TopEdit (www.topeditsci.com) for English language editing of this manuscript.
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.
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/fgene.2022.966092/full#supplementary-material
References
Akhter, S., Aziz, R. K., and Edwards, R. A. (2012). PhiSpy: a novel algorithm for finding prophages in bacterial genomes that combines similarity- and composition-based strategies. Nucleic Acids Res. 40, e126. doi:10.1093/nar/gks406
Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., and Butler, H. (2000). Gene ontology: tool for the unification of biology Nucleic Acids Res. 25, 25–29. doi:10.1038/75556
Berlin, K., Koren, S., Chin, C. S., Drake, J. P., Landolin, J. M., and Phillippy, A. M. (2015). Assembling large genomes with single-molecule sequencing and locality-sensitive hashing. Nat. Biotechnol. 33, 623–630. doi:10.1038/nbt.3238
Bertelli, C., and Brinkman, F. S. L. (2018). Improved genomic island predictions with IslandPath-DIMOB. BIOINFORMATICS 34, 2161–2167. doi:10.1093/bioinformatics/bty095
Birney, E., Clamp, M., and Durbin, R. (2004). GeneWise and genomewise. Genome Res. 14, 988–995. doi:10.1101/gr.1865504
Chen, K., Wang, L., Chen, H., Zhang, C., Wang, S., Chu, P., et al. (2021). Complete genome sequence analysis of the peanut pathogen Ralstonia solanacearum strain Rs-P.362200. BMC Microbiol. 21, 118. doi:10.1186/s12866-021-02157-7
Cheng, D., Zhou, D., Wang, Y., Wang, B., He, Q., Song, B., et al. (2021). Ralstonia solanacearum type III effector RipV2 encoding a novel E3 ubiquitin ligase (NEL) is required for full virulence by suppressing plant PAMP-triggered immunity. Biochem. Biophys. Res. Commun. 550, 120–126. doi:10.1016/j.bbrc.2021.02.082
Chin, C., Alexander, D. H., Marks, P., Klammer, A. A., Drake, J., Heiner, C., et al. (2013). Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data. Nat. Methods 10, 563–569. doi:10.1038/nmeth.2474
Coll, N. S., and Valls, M. (2013). Current knowledge on the Ralstonia solanacearum type III secretion system. Microb. Biotechnol. 6, 614–620. doi:10.1111/1751-7915.12056
Consortium, U. (2018). Corrigendum UniProt: the universal protein knowledgebase. Nucleic acids Res. 45, 158–169. doi:10.1093/nar/gky092/4841658
Deng, Y., Jianqi, L. I., Songfeng, W. U., Yunping, Z., Yaowen, C., Fuchu, H. E., et al. (2006). Integrated nr database in protein annotation system and its localization. Comput. Eng. 32, 71–73. doi:10.3969/j.issn.1000-3428.2006.05.026
Edgar, R. C. (2007). PILER-CR: fast and accurate identification of CRISPR repeats. BMC Bioinforma. 8, 18. doi:10.1186/1471-2105-8-18
Eid, J., Fehr, A., Gray, J., Rank, D., Baybayan, P., Bettman, B., et al. (2009). Real-time DNA sequencing from single polymerase molecules. SCIENCE 323, 133–138. doi:10.1126/science.1162986
Faino, L., Seidl, M. F., Datema, E., van den Berg, G. C., Janssen, A., Wittenberg, A. H., et al. (2015). Single-molecule real-time sequencing combined with optical mapping yields completely finished fungal genome. mBio 6, e00936-15. doi:10.1128/mBio.00936-15
Genin, S., and Denny, T. P. (2012). Pathogenomics of the Ralstonia solanacearum species complex. Annu. Rev. Phytopathol. 50, 67–89. doi:10.1146/annurev-phyto-081211-173000
Gui, C., Tahira, S., Xie, W., Zhuo, T., Hu, X., Zou, H., et al. (2021). The role of effector RipQ Q for virulence in Ralstonia solanacearum GMI1000. ACTA PHYTOPATHOL. SIN., 1–17. doi:10.13926/j.cnki.apps.000777
Hyatt, D., Chen, G., Locascio, P. F., Land, M. L., Larimer, F. W., Hauser, L. J., et al. (2010). Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinforma. 11, 119. doi:10.1186/1471-2105-11-119
Im, S. M., Yu, N. H., Joen, H. W., Kim, S. O., Park, H. W., Park, A. R., et al. (2020). Biological control of tomato bacterial wilt by oxydifficidin and difficidin-producing Bacillus methylotrophicus DR-08. Pestic. Biochem. Physiol. 163, 130–137. doi:10.1016/j.pestbp.2019.11.007
Jones, D. T., Taylor, W. R., and Thornton, J. M. (1992). The rapid generation of mutation data matrices from protein sequences. Comput. Appl. Biosci. 8, 275–282. doi:10.1093/bioinformatics/8.3.275
Kanehisa, M., Goto, S., Kawashima, S., Okuno, Y., and Hattori, M. (2004). The KEGG resource for deciphering the genome. Nucleic Acids Res. 32, 277D–280D. doi:10.1093/nar/gkh063
Kumar, S., Stecher, G., and Tamura, K. (2016). MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870–1874. doi:10.1093/molbev/msw054
Landry, D., González Fuente, M., Deslandes, L., and Peeters, N. (2020). The large, diverse, and robust arsenal of Ralstonia solanacearum type III effectors and their in planta functions. Mol. Plant Pathol. 21, 1377–1388. doi:10.1111/mpp.12977
Li, L., J, C., and Jr, S. (2003). OrthoMCL: Identification of ortholog groups for eukaryotic. Genome Res. 193, 2178–2189. doi:10.1101/gr.1224503
Lin, C. (2018). The antibacterial mechanism of MgONPs to Ralstonia solanacearum and its control effect on tobacco bacterial wilt. Master thesis. Chongqing: Southwest University.
Liu, Y., Qin, X., Wang, M., Liu, Y., Zhu, M., Ji, G., et al. (2007). The disease survey and pathogen isolation of tobacco bacterial wilt in yunnan province, China. Chin. Agric. Sci. Bull. 43, 311–314.
Louwen, R., Staals, R. H. J., Endtz, H. P., Baarlen, V. P., and Oost, V. D. J. (2014). The role of CRISPR-cas systems in virulence of pathogenic bacteria. Microbiol. Mol. Biol. Rev. 78, 74–88. doi:10.1128/MMBR.00039-13
Paudel, S., Dobhal, S., Alvarez, A. M., and Arif, M. (2020). Taxonomy and phylogenetic research on Ralstonia solanacearum species complex: A complex pathogen with extraordinary economic consequences. Pathogens 9, 886. doi:10.3390/pathogens9110886
Peeters, N., Carrère, S., Anisimova, M., Plener, L., Cazalé, A., and Genin, A. S. (2013). Repertoire, unified nomenclature and evolution of the Type III effector gene set in the Ralstonia solanacearum species complex. BMC Genomics 14, 859. doi:10.1186/1471-2164-14-859
Poueymiro, M., Cazalé, A. C., François, J. M., Parrou, J. L., Peeters, N., and Genin, S. (2014). A Ralstonia solanacearum type III effector directs the production of the plant signal metabolite trehalose-6-phosphate. mBio 5, e02065–14. doi:10.1128/mBio.02065-14
Prior, P., Ailloud, F., Dalsing, B. L., Remenant, B., Sanchez, B., and Allen, C. (2016). Genomic and proteomic evidence supporting the division of the plant pathogen Ralstonia solanacearum into three species. BMC GENOMICS 17, 90. doi:10.1186/s12864-016-2413-z
Salanoubat, M., Genin, S., Artiguenave, F., Gouzy, J., Mangenot, S., ArlatM., , et al. (2002). Genome sequence of the plant pathogen Ralstonia solanacearum. NATURE 415, 497–502. doi:10.1038/415497a
Shutian, L. (2020). Detecting the horizontal gene transfer for microbial genomes and clustering the biological sequences. Fuzhou, China: Fujian Normal University.
Sun, Y., Li, P., Shen, D., Wei, Q., He, J., and Lu, Y. (2019). The Ralstonia solanacearum effector RipN suppresses plant PAMP-triggered immunity, localizes to the endoplasmic reticulum and nucleus, and alters the NADH/NAD(+) ratio in Arabidopsis. Mol. Plant Pathol. 20, 533–546. doi:10.1111/mpp.12773
Syvanen, M. (1985). Cross-species gene transfer; implications for a new theory of evolution. J. Theor. Biol. 112, 333–343. doi:10.1016/s0022-5193(85)80291-5
Tan, X., Qiu, H., Li, F., Cheng, D., Zheng, X., Wang, B., et al. (2019). Complete genome sequence of sequevar 14M Ralstonia solanacearum strain HA4-1 reveals novel type III effectors acquired through horizontal gene transfer. Front. Microbiol. 10, 1893. doi:10.3389/fmicb.2019.01893
Tarailo Graovac, M., and Chen, N. (2009). Using repeat masker to identify repetitive elements in genomic sequences. Curr. Protoc. Bioinforma. 25, 410. doi:10.1002/0471250953.bi0410s25
Tatusov, R. L., Galperin, M. Y., Natale, D. A., and Koonin, E. V. (2000). The COG database: a tool for genome-scale analysis of protein functions and evolution. Nucleic Acids Res. 28, 33–36. doi:10.1093/nar/28.1.33
Villa, J. E., Tsuchiya, K., Horita, M., Natural, M., Opina, N., and Hyakumachi, M. (2005). Phylogenetic relationships of Ralstonia solanacearum species complex strains from Asia and other continents based on 16S rDNA, endoglucanase, and hrpB gene sequences. J. Gen. Plant Pathol. 71, 39–46. doi:10.1007/s10327-004-0156-1
Wang, S. (2018). Cloning and characterization of NtWRKY45 and NtRRS4 genes resistant to bacterial wilt in tobacco. Fuzhou, China: Fujian Agriculture and Forestry University.
Wang, Y., Tang, H., Debarry, J. D., Tan, X., Li, J., Wang, X., et al. (2012). MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 40, e49. doi:10.1093/nar/gkr1293
Xian, L., Yu, G., Wei, Y., Rufian, J. S., Li, Y., Zhuang, H., et al. (2020). A bacterial effector protein hijacks plant metabolism to support pathogen nutrition. Cell Host Microbe 28, 548–557. doi:10.1016/j.chom.2020.07.003
Xianyang, Z., Liuyi, P., Ming, C., Yongqi, F., Miaolian, X., and Xiaochun, Z. (2018). Inductive effect of methyl jasmonate to bacterial wilt and the effects on the activities of antioxidant enzymes in pepper seedlings. J. Plant Prot. 45, 1103–1111. doi:10.13802/j.cnki.zwbhxb.2018.2018128
Xiao, Y., Liu, X., Meng, D., Tao, J., Gu, Y., Yin, H., et al. (2018). The role of soil bacterial community during winter fallow period in the incidence of tobacco bacterial wilt disease. Appl. Microbiol. Biotechnol. 102, 2399–2412. doi:10.1007/s00253-018-8757-3
Xiaoqiang, C. (2018). Whole genome analysis of Ralstonia solanacearum FJAT-91 and function study on the extracellular polysaccharide. Fuzhou, China: Fujian Agriculture and Forestry University.
Zang, R., Song, L., Yin, X., Xu, C., Geng, Y., and Meng, Z. (2021). Genome-wide prediction and analysis of the secreted proteins of Botryosphaeria dothidea. ACTA PHYTOPATHOL. SIN. 51, 559–571. doi:10.13926/j.cnki.apps.000732
Zhang, C., Chen, H., Cai, T., Deng, Y., Zhuang, R., Zhang, N., et al. (2017). Overexpression of a novel peanut NBS-LRR geneAhRRS5 enhances disease resistance toRalstonia solanacearum in tobacco. Plant Biotechnol. J. 15, 39–55. doi:10.1111/pbi.12589
Zhang, C., Chen, H., Zhuang, R. R., Chen, Y. T., Deng, Y., Cai, T. C., et al. (2019). Overexpression of the peanut CLAVATA1-like leucine-rich repeat receptor-like kinase AhRLK1 confers increased resistance to bacterial wilt in tobacco. J. Exp. Bot. 70, 5407–5421. doi:10.1093/jxb/erz274
Zhong, Z., Rao, X., and Le, S. (2021). Molecular mechanisms of bacterialresistance to bacteriophage infection. Microbiology China 48, 3249–3260. doi:10.13344/j.microbiol.china.210483
Keywords: Ralstonia solanacearum, tobacco, bacterial wilt, disease resistance, effector proteins, genome sequencing, Rs_T3E_Hyp14
Citation: Chen K, Zhuang Y, Wang L, Li H, Lei T, Li M, Gao M, Wei J, Dang H, Raza A, Yang Q, Sharif Y, Yang H, Zhang C, Zou H and Zhuang W (2022) Comprehensive genome sequence analysis of the devastating tobacco bacterial phytopathogen Ralstonia solanacearum strain FJ1003. Front. Genet. 13:966092. doi: 10.3389/fgene.2022.966092
Received: 10 June 2022; Accepted: 18 July 2022;
Published: 22 August 2022.
Edited by:
Yuanhu Xuan, Shenyang Agricultural University, ChinaReviewed by:
Mei Qiong, Shenyang Agricultural University, ChinaElaine Chan, The University of Hong Kong, Hong Kong SAR, China
Copyright © 2022 Chen, Zhuang, Wang, Li, Lei, Li, Gao, Wei, Dang, Raza, Yang, Sharif, Yang, Zhang, Zou and Zhuang. 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: Huasong Zou, hszou@fafu.edu.cn; Weijian Zhuang, weijianz@fafu.edu.cn
†These authors have contributed equally to this work