key: cord-0026286-y1iilu5s
authors: Xu, Yingwei; Yang, Lianzhi; Wang, Yaping; Zhu, Zhuoying; Yan, Jizhou; Qin, Si; Chen, Lanming
title: Prophage-encoded gene VpaChn25_0734 amplifies ecological persistence of Vibrio parahaemolyticus CHN25
date: 2022-01-22
journal: Curr Genet
DOI: 10.1007/s00294-022-01229-z
sha: 35ef95bbf11967dbc3c713893e9070e7411f86b0
doc_id: 26286
cord_uid: y1iilu5s

Vibrio parahaemolyticus is a waterborne pathogen that can cause acute gastroenteritis, wound infection, and septicemia in humans. The molecular basis of its pathogenicity is not yet fully understood. Phages are found most abundantly in aquatic environments and play a critical role in horizontal gene transfer. Nevertheless, current literature on biological roles of prophage-encoded genes remaining in V. parahaemolyticus is rare. In this study, we characterized one such gene VpaChn25_0734 (543-bp) in V. parahaemolyticus CHN25 genome. A deletion mutant ΔVpaChn25_0734 (543-bp) was obtained by homologous recombination, and a revertant ΔVpaChn25_0734-com (543-bp) was also constructed. The ΔVpaChn25_0734 (543-bp) mutant was defective in growth and swimming mobility particularly at lower temperatures and/or pH 7.0–8.5. Cell surface hydrophobicity and biofilm formation were significantly decreased in the ΔVpaChn25_0734 (543-bp) mutant (p < 0.05). Based on the in vitro Caco-2 cell model, the deletion of VpaChn25_0734 (543-bp) gene significantly reduced the cytotoxicity of V. parahaemolyticus CHN25 to human intestinal epithelial cells (p < 0.05). Comparative secretomic and transcriptomic analyses revealed a slightly increased extracellular proteins, and thirteen significantly changed metabolic pathways in the ΔVpaChn25_0734 (543-bp) mutant, showing down-regulated carbon source transport and utilization, biofilm formation, and type II secretion system (p < 0.05), consistent with the observed defective phenotypes. Taken, the prophage-encoded gene VpaChn25_0734 (543-bp) enhanced V. parahaemolyticus CHN25 fitness for survival in the environment and the host. The results in this study facilitate better understanding of pathogenesis and genome evolution of V. parahaemolyticus, the leading sea foodborne pathogen worldwide. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s00294-022-01229-z.

Vibrio parahaemolyticus is a Gram-negative bacterium that resides in aquatic environments worldwide (Letchumanan et al. 2014; Su and Chen 2020) . Ingestion of raw, undercooked, or mishandled seafood contaminated by pathogenic V. parahaemolyticus can cause acute gastroenteritis in humans and even death (Ghenem et al. 2017; Baker-Austin et al. 2018) . The bacterium was first identified from semidried juvenile sardines in 1950 in Japan, which caused 272 cases of acute diarrhea and 20 death (Fujino et al. 1953) . Communicated by Michael Polymenis. Afterward, outbreaks of human gastroenteritis caused by V. parahaemolyticus occurred in many Asian nations and consequently around the world (Elmahdi et al. 2016; Meparambu Prabhakaran et al. 2020) . It was estimated that more than 50% of foodborne gastrointestinal-Vibriosis cases were caused by V. parahaemolyticus in the United States (Karan et al. 2021) , which led to 45,000 illnesses each year (https:// www. cdc. gov/ Vibrio/, accessed on 26 April 2021). Sporadic outbreaks have also been reported in coastal European nations (Meparambu Prabhakaran et al. 2020) . China is one of the 15 countries with high incidence of diarrhea disease (Walker et al. 2013) . Recently, a national surveillance of 152,792 patients of all ages with acute diarrhea was administered in 217 hospitals and 93 reference laboratories in 31 provinces of China by the Chinese Center for Diseases Control and Prevention in 2009 (Wang et al. 2021a . Total 13 common bacterial pathogens were tested, among which V. parahaemolyticus ranked third and contributed to 10.83% of all positive detection. Pathogenic V. parahaemolyticus strains produce thermostable direct hemolysin (TDH) and/ or the TDH-related hemolysin (TRH) (Leoni et al. 2016; Cai and Zhang 2018) . However, some clinical V. parahaemolyticus strains do not produce the two major toxins but display virulence, indicating that putative toxic determinants exist. Their pathogenicity might be achieved with different strategies employed by different strains (Li et al. 2020) . Therefore, it is imperative to identify such risk factors in V. parahaemolyticus to assure food safety and human health.

Phages are found most abundantly in aquatic environments worldwide (Tan et al. 2021) . Temperate phages can integrate into bacterial chromosomes with carried genes and drive bacterial genome evolution (Garin-Fernandez and Wichels 2020; Nuidate et al. 2021; Wendling et al. 2021) . Vibrio spp. virulence-associated genes can be combined at high frequencies by horizontal gene transfer (HGT), leading to the emergence of pandemic or pathogenic clones (Castillo et al. 2018) . For example, in early 1996, an atypical outbreak of V. parahaemolyticus infection in India was associated with a new serotype O3:K6 that carried genetic markers tdh, toxRS new , and orf8. Dispersal of this new pathogenic clone has been reported in many countries around five continents, suggesting that V. parahaemolyticus from the serotype O3:K6 was pandemic (Santos et al. 2021 ). The orf8 gene was possibly acquired from a filamentous phage f237, which endorsed the bacterium more adhesive to host intestinal cells (Santos et al. 2021) . Recently, Grain-Fernandez et al. reported a filamentous phage vB_VpaI_VP-3218 in V. parahaemolyticus . It integrates into the host genome with potential zonula occludens toxin genes . Another newly isolated phage vB_ VpaP_MGD2 from clam Meretrix meretrix was reported to be a candidate biocontrol agent against early mortality syndrome/acute hepatopancreatic necrosis disease caused by multidrug-resistant V. parahaemolyticus in shrimp (Cao et al. 2021) .

In our previous research, V. parahaemolyticus CHN25 (serotype: O5: KUT) of aquatic animal origin was identified and characterized (Song et al. 2013; Sun et al. 2014; He et al. 2015; Zhu et al. 2017 Zhu et al. , 2020 Yang et al. 2020) . The bacterial genome (5,443,401 bp, 45 .2% G+C) contains five prophage gene clusters (6.5-36.6 kb) (Zhu et al. 2017) , the largest of which showed sequence similarity to a Vibrio phage Martha 12B12 (Zhu et al. 2017) . Within the cluster, a gene VpaChn25_0724 (294-bp) encoding a hypothetical protein has been identified recently ). In the present study, biological function of another prophageencoded gene VpaChn25_0734 (543-bp) within the same gene cluster was characterized. It showed sequence identity with components of a Phage_tail_S superfamily (Lu et al. 2020) . The objectives of this study were (1) to knock out the VpaChn25_0734 (543-bp) gene from V. parahaemolyticus CHN25 genome by homologous recombination to contract a deletion mutant ΔVpaChn25_0734 (543-bp). Meanwhile, a revertant ΔVpaChn25_0734-com (543-bp) was also constructed;

(2) to examine growth, mobility, biofilm formation, and cell toxicity of the ΔVpaChn25_0734 (543-bp) mutant compared with V. parahaemolyticus CHN25 wild type (WT) and ΔVpaChn25_0734-com (543-bp) strains; (3) to decipher the possible molecular mechanism underlying altered phenotypes of the ΔVpaChn25_0734 (543-bp) mutant by comparative secretomic and transcriptomic analysis. The results in this study facilitate better understanding of biological function of prophage-encoded genes in V. parahaemolyticus genomes.

Vibrio parahaemolyticus CHN25 (Song et al. 2013; Sun et al. 2014; He et al. 2015; Zhu et al. 2017 Zhu et al. , 2020 Yang et al. 2020 ) was stored at − 80 °C freezer in our research group. Escherichia coli DH5α λpir (BEINUO Biotech (Shanghai) Co. Ltd., China), E. coli β2155 λpir (kindly provided by Professor Weicheng Bei) and plasmid pDS132 (kindly provided by Professor Dominique Schneider) were used as a host strain for gene cloning, and a donor strain and a suicide vector for conjugation and gene knockout experiments, respectively (Zhu et al. 2017; Yang et al. 2020) . Plasmid pMMB207 (Biovector Science Lab, Inc., China) was used as a gene expression vector (Zhu et al. 2017; Yang et al. 2020) . V. parahaemolyticus and E. coli strains were incubated in the same media and conditions as described previously (Zhu et al. 2017; Yang et al. 2020 ).

Gene deletion mutant and revertant of V. parahaemolyticus CHN25 were constructed according to the methods described previously (Zhu et al. 2017; Yang et al. 2020) . Briefly, two primer pairs VpaChn25_0734 (543-bp)-up-F/R and VpaChn25_0734 (543-bp)-down-F/R (Table 1) were designed to target upstream (486 bp) and downstream (494 bp) sequences of the VpaChn25_0734 (543-bp) gene in V. parahaemolyticus CHN25 genome, respectively. Restriction endonucleases XbaI and SacI (TaKaRa, Japan) were introduced into polymerase chain reaction (PCR) products, which were inserted into corresponding cloning sites of the pDS132. The obtained recombinant plasmid pDS132 + VpaChn25_0734 (543-bp) was then transformed into diaminopimelic acid (DAP) auxotroph E. coli β2155 incubated in LB medium supplemented with 0.3 mM DAP (Sigma-Aldrich, USA) (Zhu et al. 2017; Yang et al. 2020) . Double-crossover deletions of the VpaChn25_0734 (543-bp) gene were screened using the colony PCR assay (Zhu et al. 2017; Yang et al. 2020) with the primer pair VpaChn25_0734 (543-bp)-up-ex-F/R (Table 1) . The obtained ΔVpaChn25_0734 (543-bp) mutant was confirmed as described previously (Zhu et al. 2017; Yang et al. 2020) .

Additionally, the VpaChn25_0734 (543-bp) gene was obtained by PCR using the primer pair VpaChn25_0734com (543-bp)-F/-R (Table 1 ). The PCR product was inserted into the pMMB207, and then transformed into E. coli DH5α. The obtained recombinant plasmid pMMB207 + VpaChn25_0734 (543-bp) was electro-transformed into the ΔVpaChn25_0734 (543-bp) mutant (Zhu et al. 2017; Yang et al. 2020) . Positive electro-transformants ΔVpaChn25_0734-com (543-bp) were screened, and confirmed (Zhu et al. 2017; Yang et al. 2020) . DNA sequencing was conducted by Shanghai Sangon Biological Engineering Technology and Services Co., Ltd., Shanghai, China.

Growth of V. parahaemolyticus strains in the LB or TSB media at different temperatures (37 °C, 25 °C, and 15 °C) for different time (24-60 h) was individually measured using Bioscreen Automated Growth Curve Analyzer (BioTek, . Growth curves of V. parahaemolyticus strains in the TSB medium in different pH conditions (pH 5.5-8.5) were also examined, respectively (Sun et al. 2014; Zhu et al. 2017; Yang et al. 2020 ).

Swimming motility of V. parahaemolyticus strains on semisolid TSB (pH 8.5, 3% NaCl) agar plates (0.25% agar) was individually measured as described previously . The TSB agar plates were incubated at 37 °C, 25 °C, and 15 °C for 12 h, 24 h, and 72 h, respectively. Diameters at bacterial colonies were measured and recorded. Biofilm formation of V. parahaemolyticus strains was individually examined using the crystalline violet staining method . The phosphate-buffered saline (PBS, pH 7.2-7.4), and crystal violet were purchased from Sangon (China).

Outer membrane permeability was determined according to the method by Wang et al. (Wang et al. 2021c ). Briefly, a 200 μL/well of bacterial suspension in TSB medium (pH 8.5, 3% NaCl) at mid-logarithmic growth phase (mid-LGP) at 37 °C was mixed with a 2 μL/well of 10 mm N-phenyl-1-naphthylamine (NPN) (Sangon, China). The change of fluorescence intensity per well was measured using BioTek Synergy 2 (BioTek, USA). The excitation and emission wavelengths were 350 nm and 420 nm, respectively (Wang et al. 2021c ). Inner membrane permeability was also determined (Huang et al. 2021 ). Briefly, a 200 μL/well of bacterial suspension was mixed with a 2.5 μL/well of 10 mM O-nitrophenyl-β-d-galactopyranoside (O-nitrophenyl)-β-dgalactopyranoside (ONPG) (Sangon, China). The mixture was incubated at 37 °C, and absorbance at OD 415 nm of each well was measured using BioTek synergy 2 (BioTek, USA) every 30 min (Huang et al. 2021) . Cell surface hydrophobicity and membrane fluidity assays were performed as described previously ).

Extraction of extracellular proteins of V. parahaemolyticus strains grown in TSB medium (pH 8.5, 3% NaCl) at 37 °C in static culture condition, and two-dimensional gel electrophoresis (2D-GE) were performed according to the methods described previously . Briefly, a 20 mg of each sample was applied to pH gradient gel strips (pH 4-7, 7 cm, Bio-Rad, Inc., USA) at 17 °C for passive rehydration for 16 h. Isoelectric focusing (IEF) for the first dimensional separation was performed using a six-step procedure, while sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE, 12.5% separation gels) was performed for the second dimensional separation . Amino acid sequences of protein spots were determined using liquid chromatography-tandem mass spectrometry (LC-MS/ MS) by Shanghai Houji Biological Co., Shanghai, China .

The in vitro cell modal assay was performed according to the method described previously 

Total RNA of each of V. parahaemolyticus strains grown in TSB medium (pH 8.5, 3% NaCl) to mid-LGP at 37 °C were individually extracted, and purified as described previously . Three independently prepared RNA samples were subjected for Illumina RNA-sequencing, which was carried out by Shanghai Majorbio Bio-pharm Technology Co. Ltd., China using Illumina HiSeq 2500 platform . Quality filtration of raw RNA-seq data, clean read aligning, and defining of differentially expressed genes (DEGs) were conducted as described previously . Gene set enrichment analysis (GSEA) of DEGs was also conducted . Representative DEGs were examined using quantitative reverse transcription-PCR (RT-PCR) assay ) with the primers listed in Table 1 . The data were analyzed using SPSS version 17.0 software (SPSS Inc., USA). All tests in this study were performed in at least triplicate. Phylogenetic tree was constructed using MEGA 5.0 software (Kumar et al. 2018 ).

Cell structure of V. parahaemolyticus strains grown in TSB medium (pH 8.5, 3% NaCl) to mid-LGP at 37 °C were observed and recorded using Schottky Field Emission Scanning Electron Microscope (SU5000, Japan, 5.0 kV × 35,000) at College of Food Science and Technology, Shanghai Ocean University (Shanghai, China).

The largest prophage gene cluster contains 24 predicted open read frames (ORFs) in V. parahaemolyticus CHN25 chromosome 1 (3,416,467 bp, from 816,554 to 846,961 bp) (Zhu et al. 2017 ). It has sequence similarity to the Vibrio phage Martha 12B12 (30,408 bp, GenBank accession no. NC_021070) . The VpaChn25_0734 (543bp) gene is the sixteenth ORF encodes a predicted phage virion morphogenetic protein with conserved structural domains belonging to the Phage_tail_S superfamily. In the cellular component catalogue, it was involved in constituting the intracellular membrane-bounded organelle (GO: 0043231) with a predicted GO-Score of 0.42.

To characterize biological roles of the prophageencoded gene VpaChn25_0734 (543-bp) remaining in V. parahaemolyticus CHN25 genome, an untagged in-frame gene deletion mutant ΔVpaChn25_0734 (543-bp) was obtained by homologous recombination (see the Materials and Methods). The VpaChn25_0734 (543-bp) gene was knocked out from the V. parahaemolyticus CHN25 genome, which was verified by PCR and DNA sequencing (data not shown), RT-qPCR, as well as transcriptome analysis (see below).

Additionally, a reverse mutant ΔVpaChn25_0734com (543-bp) was also successfully obtained. The VpaChn25_0734 (543-bp) gene was amplified, and cloned into the expression vector pMMB207, resulting in the recombinant vector pMMB207 + VpaChn25_0734 (543-bp). This recombinant vector was electro-transformed into the ΔVpaChn25_0734 (543-bp) mutant to obtain the revertant ΔVpaChn25_0734-com (543-bp), which was confirmed by the methods described above.

To investigate the influence of VpaChn25_0734 (543-bp) gene deletion on V. parahaemolyticus CHN25 survival, we determined growth curves of the WT, ΔVpaChn25_0734 (543-bp), and ΔVpaChn25_0734-com (543-bp) strains at 37 °C, 25 °C, and 15 °C, which V. parahaemolyticus experiences during life cycle (Maje et al. 2020; Yang et al. 2020; Ali et al. 2021) . At the optional growth temperature of 37 °C, the ΔVpaChn25_0734 (543-bp) mutant grew in TSB medium (pH 8.5, 3% NaCl) with a retardation phase (RP) of 2 h, when compared with the WT strain ( Fig. 1a) . Similarly, at 25 °C and 15 °C, the RPs of ΔVpaChn25_0734 (543-bp) were sixfold and 2.67-fold longer than those of WT, respectively (Fig. 1b, c) . Additionally, the revertant ΔVpaChn25_0734-com (543-bp) appeared partially complement the defective phenotype of the ΔVpaChn25_0734 (543-bp) mutant. These results indicated that the VpaChn25_0734 (543-bp) gene enhanced V. parahaemolyticus CHN25 fitness for thriving particularly at lower temperatures.

Acid tolerance is critical for foodborne pathogens to survive against host acid gastric fluid (pH 1-3, but above 6.0 after food consumption) (Sun et al. 2014 ). Therefore, we examined the growth of the three strains in TSB medium (3% NaCl) with pH values in the range from 5.5 to 8.0, and the results were illustrated in Fig. 2a -f. In acidic conditions (pH 5.5-6.5), the growth of WT, ΔVpaChn25_0734 (543-bp), and ΔVpaChn25_0734-com (543-bp) strains were greatly inhibited, and the maximum OD 600 nm values at stationary grow phage (SGP) were below 0.5 ( Fig. 2a-c) ; in the neutral condition (pH 7.0), the growth of ΔVpaChn25_0734 (543bp) and ΔVpaChn25_0734-com (543-bp) strains were also obviously inhibited compared to the WT strain (Fig. 2d ); in alkaline conditions (pH 7.5-8.0), ΔVpaChn25_0734 (543-bp) still grew poorly with a longer RP and lower biomass than WT (Fig. 2e , f). Additionally, the plasmid-born ΔVpaChn25_0734-com (543-bp) did not fully complement the defective phenotype of the ΔVpaChn25_0734 (543-bp) mutant. The results indicated that the VpaChn25_0734 (543-bp) gene amplified V. parahaemolyticus CHN25 environmental persistance at pH 7.0-8.0.

Vibrio parahaemolyticus is found to grow as swimming cells in liquid environments (Freitas et al. 2020) . In this study, swimming mobility of the WT, ΔVpaChn25_0734 (543-bp), and ΔVpaChn25_0734-com (543-bp) strains were examined at different temperatures. As shown in Fig. 3a-1 to c-1, no significant difference in swimming circles was observed among the three strains grown in semi-solid TSB agar medium (pH 8.5, 3% NaCl, 0.25% agar) at 37 °C (p > 0.05). However, at 25 °C, the ΔVpaChn25_0734 (543bp) mutant swam significantly slower (8.00 ± 0.71 mm) than the WT (12.75 ± 2.86 mm), and ΔVpaChn25_0734-com (543-bp) (9.00 ± 1.2 mm) strains, respectively (p < 0.05) ( Fig. 3a -2 to c-2). The similar case was observed at 15 °C ( Fig. 3a -3 to c-3). These results indicated that the deletion of VpaChn25_0734 (543-bp) gene significantly inhibited swimming mobility of V. parahaemolyticus CHN25 at the lower temperatures. Fig. 2 Survival of the WT, ΔVpaChn25_0734 (543-bp), and ΔVpaChn25_0734-com (543-bp) strains in different pH conditions. a-f pH 5.5 (a), pH 6.0 (b), pH 6.5 (c), pH 7.0 (d), pH 7.5 (e), and pH 8.0 (f), respectively

Given the growth inhibition at the lower temperatures, the dynamic process of biofilm formation of the three strains was examined at 37 °C for 60 h using the crystalline violet staining method. All three strains established biofilms at three different stages (development, maturation, and diffusion) when grown in TSB medium (pH 8.5, 3% NaCl) in static culture condition, which was consistent with the previous report . Nevertheless, the biofilm biomass formed by the ΔVpaChn25_0734 (543-bp) mutant was significantly less than the WT strain at development (24 h), maturation (48 h), and diffusion (60 h) stages (p < 0.05) (Fig. 4) . 

To investigate the impact of VpaChn25_0734 (543-bp) gene deletion on cell structure of V. parahaemolyticus CHN25, we further examined cell membrane permeability and fluidity, and surface hydrophobicity of the three strains. As shown in Fig. 5 , no significant difference in membrane permeability and fluidity was observed among the three strains grown in TSB medium (pH 8.5, 3% NaCl) at 37 °C (p > 0.05) (Fig. 5a-c) . However, the cell surface hydrophobicity of the ΔVpaChn25_0734 (543-bp) mutant was significantly decreased (2.52-fold) when compared to the WT (p < 0.01) (Fig. 5d ).

The ΔVpaChn25_0734 (543-bp) mutant also grew more slowly than the WT strain when cultured in the TSB medium (pH 8.5, 3% NaCl) in static culture condition at 37 °C (data not shown), consistent with the results of the biofilm formation experiments. As shown in Fig. 6 , slightly changed secretome profiles among the three strains were observed, based on the 2D-GE analysis. Three independent parallel 2D-GE gels of each sample produced identical results (data not shown). The different protein spots were excised from the 2D-GE gels and further identified by LC-MS/MS analysis. As shown in Fig. 6 , two differential extracellular proteins secreted by the ΔVpaChn25_0734 (543-bp) mutant were identified: spot 34-b-1 was identified as a 30S ribosomal protein S1 (RpsA); and Spot 34-b-2 was identified as a DNAdirected RNA polymerase subunit alpha (RpoA) ( Table 2) .

As illustrated in Fig. 7 , when compared to the infection with the WT strain, the viability of Caco-2 cells significantly increased 1.09-fold after infected with the ΔVpaChn25_0734 (543-bp) mutant at 37 °C for 4 h (p < 0.01) (Fig. 7a) . Meanwhile, Caco-2 cells were double-stained with membranelinked protein V-FITC and PI, and analyzed by flow cytometry. The results showed that the ΔVpaChn25_0734 (543-bp) mutant induced apoptosis of Caco-2 cells at a significantly reduced rate (0.87-fold) than the WT strain after the infection for 4 h (p < 0.01) (Fig. 7b) . These results indicated a decrease in cytotoxicity of V. parahaemolyticus CHN25 in the absence of the VpaChn25_0734 (543-bp) gene.

To investigate global-level gene expression change mediated by the VpaChn25_0734 (543-bp) gene deletion, transcriptomes of the WT, ΔVpaChn25_0734 (543-bp), and ΔVpaChn25_0734-com (543-bp) strains were determined, when incubated in TSB medium (pH 8.5, 3% NaCl) at 37 °C. Comparative transcriptome analysis showed that approximately 17.03% of the genes in the ΔVpaChn25_0734 (543-bp) mutant were differentially expressed, which were classified into various gene functional catalogues (Fig. 8) . A complete list of the DEGs in the three strains was deposited in the NCBI SRA database (http:// www. ncbi. nlm. nih. gov/ sra/) under the accession number PRJNA733855. Expression of representative DEGs in the ΔVpaChn25_0734 (543-bp) mutant were examined by the RT-PCR assay, and the resulting data were consistent with the transcriptomic analysis (Table S1 ).

Approximately thirteen significantly changed metabolic pathways were identified in the VpaChn25_0734 (543-bp) mutant, including the bacterial secretion system; terpenoid backbone biosynthesis; biofilm formation; pentose and glucuronate interconversions; fructose and mannose metabolism; purine metabolism; valine, leucine, and isoleucine degradation; pyruvate metabolism; peptidoglycan biosynthesis; valine, leucine and isoleucine biosynthesis, phosphotransferase system (PTS), two-component system (T-CS), and quorum sensing (QS) (Fig. 8) .

Remarkably, approximately 81 DEGs involved in 10 of the 13 metabolic pathways (except the terpenoid backbone biosynthesis, purine metabolism, and T-CS) were all significantly down-regulated in the ΔVpaChn25_0734 (543-bp) mutant (0.051-to 0.500-fold) (p < 0.05), when compared with the WT and ΔVpaChn25_0734-com (543-bp) strains (Table 3 ).

In the pentose and glucuronate interconversions, expression of two DEGs were significantly repressed at the transcriptional level (0.423-, and 0.462-fold) (p < 0.05), which encoded an l-arabinose isomerase (VpaChn25_RS23490), and an oligogalacturonate lyase family protein (VpaChn25_ RS16200), respectively. In the fructose and mannose metabolism, 3 DEGs encoding key enzymes were also significantly inhibited (0.132-to 0.431-fold) (p < 0.05), including a phosphoenolpyruvate-protein phosphotransferase (VpaChn25_ RS01915), a phospho-sugar mutase (VpaChn25_RS21580), and a mannitol-1-phosphate 5-dehydrogenase (VpaChn25_ RS01930). In the pyruvate metabolism, expression of 15 DGEs, e.g., a key pyruvate kinase (VpaChn25_RS10085) and a 2-hydroxyacid dehydrogenase (VpaChn25_RS16545), were significantly depressed (0.140-to 0.488-fold). Additionally, 21 DGEs involved in T-CS were significantly downregulated (0.131-to 0.498-fold). All these metabolic pathways were involved in carbohydrate metabolism, and the overall down-regulation trend of the DGEs sets suggested under-activity of carbon source transport and/or utilization, which may have consequently led to a certain impairment of energy supply.

In the peptidoglycan biosynthesis, expression of 6 DGEs were significantly down-regulated (0.374-to 0.498-fold), including a phospho-N-acetylmuramoylpentapeptide-transferase (VpaChn25_RS02225), a UDP-N-acetylmuramoyl-l-alanine-d-glutamate ligase (VpaChn25_RS02230), a UDP-N-acetylmuramoyl-tripeptide-d-alanyl-d-alanine ligase (VpaChn25_RS02220), a PBP1A family penicillin-binding protein (VpaChn25_ RS14065), a d -alanyl-d-alanine carboxypeptidase (VpaChn25_RS23100), and a UDP-N-acetylmuramatel-alanine ligase (VpaChn25_RS02245). Inhibition of the peptidoglycan biosynthetic may have compromised cell integrity and inhibited the normal growth and division of the ΔVpaChn25_0734 (543-bp) mutant.

In bacterial secretion systems, remarkably, all 13 DEGs involved in type II secretion system (T2SS) were significantly down-regulated (0.197-to 0.420-fold) (p < 0.05), including T2SS minor pseudopilins GspH, GspI, GspG, GspK, and GspJ (VpaChn25_RS00645, VpaChn25_ RS00650, VpaChn25_RS00640, VpaChn25_RS00660, and VpaChn25_RS00655); T2SS proteins M (VpaChn25_ RS00670) and GspL (VpaChn25_ RS00665); T2SS ATPase GspE (VpaChn25_RS00630); T2SS inner membrane protein GspF (VpaChn25_RS00635); T2SS secretin GspD (VpaChn25_RS00625); and T2SS proteins (Vpachn25_ RS13760, Vpachn25_RS13775, and Vpachn25_RS13765). Unexpectedly, expression of 3 DGEs involved in type VI secretion system (T6SS) were enhanced (2.166-to 3.688fold) (p < 0.05), including T6SS protein TssL (VpaChn25_ RS20595), ATPase TssH (VpaChn25_RS20535), and lipoprotein TssJ (VpaChn25_RS20605). T6SS secrets bacterial toxins and regulates biofilm formation or stress sensing (Aschtgen et al. 2010) . It was speculated that the depressed T2SS induced by the deletion of VpaChn25_0734 (543-bp) may have triggered T6SS for more support in material transport and nutrient uptake. In the biofilm formation, expression of 9 DGEs, e.g., a catalase (VpaChn25_RS22250), a porin (VpaChn25_ RS12510), and a carbon storage regulator CsrA (VpaChn25_ RS12935), was significantly depressed (0.051-to 0.499fold) (p < 0.05). Additionally, 6 DGEs involved in the QS (0.103-to 0.485-fold) also showed a decrease in transcriptional levels in the ΔVpaChn25_0734 (543-bp) mutant. The SEM observation of the three strains provided cell structure evidence for the above results, in which more cytocysts were obviously present in the ΔVpaChn25_0734 (543-bp) mutant than the WT and ΔVpaChn25_0734-com (543-bp) strains (Fig. 9) . Comparative transcriptomic analysis also revealed a few metabolic pathways that were significantly up-regulated in the ΔVpaChn25_0734 (543-bp) mutant (p < 0.05), such as purine metabolism, and terpenoid backbone biosynthesis. Moreover, there were several DEGs with greatly changed transcription levels (12.578-to 25.569-fold) mediated by the VpaChn25_0734 (543-bp) gene deletion. For example, expression of 22 DEGs in the purine metabolism was significantly up-regulated (2.133-to 9.555-fold) (p < 0.05). Additionally, a highly up-regulated gene VpaChn25_RS18565 (17.794-fold) coded for a thiolase that catalyzes the reverse Claisen condensation reaction, the first step of the biosynthesis of sterols and ketone bodies in lipid metabolism (Anbazhagan et al. 2014) . The gene (VpaChn25_RS16830, 25.569-fold) encoding a malate synthase, a key enzyme responsible for malic acid synthesis in the glyoxylate cycle (Pua et al. 2003) , was also greatly up-regulated. Taken, these DEGs elicited by the deletion of VpaChn25_0734 (543-bp) likely impacted multiple metabolic pathways and consequently resulted in the defective phenotypes of the ΔVpaChn25_0734 (543-bp) mutant.

A total of 119 V. parahaemolyticus isolates from aquatic products collected in Shanghai, China (Su and Chen 2020) were tested for the VpaChn25_0734 (543-bp) gene by the PCR assay. The results revealed that approximately 0.84% (n = 1) of the V. parahaemolyticus isolates carried the VpaChn25_0734 (543-bp) homologue. Moreover, sequence analysis against the GenBank database showed that VpaChn25_0734 (543-bp) homologues were present in one Vibrio phage (vB_VpaM_VP-3212) (marine metagenome genome assembly, Genbank accession no.: LR700235); two Vibrio species, including Vibrio campbellii (Genbank accession no.: CP020077) and V. parahaemolyticus (Genbank accession no.: CP046787); and non-Vibrio genus Shewanella oneidensis (Genbank accession no.: CP053946). Based on these identified homologues, a phylogenetic tree was constructed (Fig. 10) . This analysis revealed that the VpaChn25_0734 (543-bp) gene was phylogenetically close to a Vc3S01_A0696 gene of V. campbellii 2013 0629003S1. V. campbellii is an emerging aquaculture pathogen that causes Vibriosis in farmed shrimp (Nuidate et al. 2021) . These results indicated that the VpaChn25_0734 (543-bp) gene exists in V. parahaemolyticus population, within the Vibrio genus, and even across bacterial genera.

Vibrio parahaemolyticus is frequently isolated from seafoods at water temperatures above 15 °C (Su and Chen 2020) . Huge amounts of phages in aquatic ecosystems allow genetic material to jump between strains and species (Lerminiaux and Cameron 2019) . Currently, biological function of prophage-encoded genes in V. parahaemolyticus is not yet fully understood. In this study, one such gene VpaChn25_0734 (543-bp) encoding a predicted phage morphogenetic protein in V. parahaemolyticus CHN25 genome was for the first time subjected to a systematic study. We successfully constructed the ΔVpaChn25_0734 (543-bp) mutant and obtained its revertant ΔVpaChn25_0734-com (543-bp). Our data indicated that the deletion of VpaChn25_0734 (543-bp) not only resulted in the growth hindrance of V. parahaemolyticus CHN25, particularly at lower temperatures, but also reduced the bacterial tolerance at pH 7.0-8.5. Flagella play a crucial role in bacterial pathogenesis linked to colonization, chemotaxis, biofilm formation, and virulence (Chevance and Hughes 2008; Erhardt 2016) . Flagellar motility is one of the high-energy-utilizing cellular processes (Khan et al. 2020) . In this study, our results indicated that the swimming mobility of ΔVpaChn25_0734 (543-bp) was significantly inhibited at lower temperatures. Moreover, the ΔVpaChn25_0734 (543-bp) mutant had a defect to a certain extent in all three stages of biofilm formation. These results indicated that the prophage-encoded VpaChn25_0734 (543bp) gene amplified the environmental persistence of V. parahaemolyticus CHN25.

Bacteria have evolved several types of secretion systems to secrete relevant substrates and effectors for bacterial adaptation to the environment and interaction with hosts (Costa et al. 2015; Wang et al. 2021b) . In this study, comparative Fig. 10 Phylogenetic relationships between the VpaChn25_0734 (543-bp) gene and its homologues secretome analysis revealed a few differentially extracellular proteins secreted by the ΔVpaChn25_0734 (543-bp) mutant. For example, Spot 34-b-1 was identified as a 30S ribosomal protein S1, which is mainly involved in protein synthesis in bacteria. Despite being present intracellularly, however, it was also detected in secreted proteins of Streptococcus suis and Haemophilus parasuis (Wu et al. 2008; Wei et al. 2014) . Spot 34-b-2 was a DNA-directed RNA polymerase subunit alpha protein, and its dimerization is the first step in the sequential assembly of subunits to form the holoenzyme (Degen et al. 2014) . Additionally, another protein Spot 34-c-1 was identified as an aspartate kinase in the biosynthesis of the aspartate family of amino acids including methionine, threonine, and lysine (Zou et al. 2020) . V. parahaemolyticus infection resulted in autophagy, cell rounding, and cell lysis, which led to proinflammatory death of infected host cells (Shimohata et al. 2011) . In this study, the in vitro cell model analysis revealed that the viability of Caco-2 cells was significantly increased after the ΔVpaChn25_0734 (543bp) mutant infected at 37 °C for 4 h, when compared with the WT strain (p < 0.05). Moreover, ΔVpaChn25_0734 induced the apoptosis of Caco-2 cells at a lower rate than WT. Taken, the differentially secreted proteins by ΔVpaChn25_0734 and or the expression of VpaChn25_0734 (543-bp) gene benefited to V. parahaemolyticus CHN25 for its infection to the host cells.

Comparative transcriptome analysis revealed 13 significantly changed metabolic pathways in the VpaChn25_0734 (543-bp) mutant. Remarkably, all the DEGs associated with the T2SS were significantly inhibited (0.197-to 0.42-fold) (p < 0.05), including GspH, GspI, GspG, GspK, GspJ, GspL, GspE, GspF, GspD, Vpachn25_RS00670, Vpachn25_RS13760, Vpachn25_RS13775, Vpachn25_ RS13765. The main function of bacterial T2SS is to obtain nutrients via exoproteins such as hydrolases to degrade biopolymers. T2SS also promote the secretion of toxins, adhesins, mucin, or cytochromes involved in respiration, motility, or biofilm formation (Nivaskumar and Francetic 2014; Yu et al. 2015; Matsuda et al. 2019) . For example, in this study, expression of T2SS minor pseudopilin GspK (VpaChn25_RS00660) was significantly decreased (0.335fold) in the ΔVpaChn25_0734 mutant. In T2SS, GspK binds on GspE N-domain via GspL and promotes GspE binding to cardiolipin to stimulate the initial rounds of ATP hydrolysis, which is required for assembly of GspG (Camberg et al. 2007; Cisneros et al. 2012) . In this study, expression of GspH (Vpachn25_RS00645) was also highly repressed (0.197-fold). In T2SS, GspH binds the initiating tip complex in vitro via its globular domain, which could ensure a transition between initiation and ATPase-catalyzed elongation required for maximal efficiency of protein secretion under native conditions (Douzi et al. 2009; Korotkov and Sandkvist 2019) . Similarly, three DEGs involved in the secretory-signal recognition particle (Sec-SRP) system were significantly down-regulated (0.188-to 0.5-fold) (p < 0.05), including a protein translocase subunit SecF (VpaChn25_RS18155), a protein translocase subunit SecD (VpaChn25_RS18150), and a protein translocase subunit SecF (VpaChn25_RS02895). These results suggested inactive nutrient uptake and/or biofilm formation of the ΔVpaChn25_0734 (543-bp) mutant. Moreover, the DEG (VpaChn25_RS22250) encoding a catalase in the biofilm system was greatly down-regulated (0.051-fold). Catalase prevents the accumulation of hydrogen peroxide and protects cellular organelles from damage by peroxide (Goyal and Basak 2010) . Additionally, the DEG (VpaChn25_RS12510) in the porin synthesis, which is a part of the biofilm system, was also decreased (0.359-fold). Porins are non-specific protein channels in bacterial outer membrane that enable the influx of hydrophilic solutes (Ades 2004; Liu et al. 2012 ). It will be interesting to further investigate cell location of the VpaChn25_0734 (543-bp)-encoding protein in V. parahaemolyticus CHN25 in the future research.

The PTS functions in carbohydrate transport and regulates numerous cellular processes (Schauder et al. 1998; Galinier and Deutscher 2017) . In this study, comparative transcriptome data revealed that the VpaChn25_RS19520, VpaChn25_RS13650, and VpaChn25_RS22265 genes in PTS were significantly down-regulated (0.18, 0.257, 0.362fold) in the VpaChn25_0734 (543-bp) mutant. Meanwhile, 20 DEGs involved in carbohydrate metabolism also showed significantly down-regulated transcription (0.132-to 0.488fold). For example, the DEG encoding a pyruvate kinase (PK, VpaChn25_RS10085) in the pyruvate metabolism was highly down-regulated (0.140-fold). PK catalyzes the last step of glycolysis, in which phosphoenolpyruvate is converted to pyruvate with the production of adenosine triphosphate (ATP) (Israelsen and Vander Heiden 2015) . Additionally, 21 DGEs involved in T-CS were also significantly down-regulated (0.131-to 0.498-fold). ATP-dependent ABC transporters actively transport molecules across the lipid membrane (Moussatova et al. 2008 ). The overall down-regulation trend of the DGEs sets in carbohydrate metabolism suggested that the deletion of VpaChn25_0734 (543-bp) gene negatively impacted nutrient acquisition and utilization, and consequently affected energy production and consumption of V. parahaemolyticus CHN25, consistent with defective phenotypes observed in this study.

Additionally, some DGEs were highly up-regulated in the VpaChn25_0734 (543-bp) mutant. For example, expression of an inosine monophosphate (IMP) dehydrogenase (VpaChn25_RS03010) was increased by 4.645-fold. IMP dehydrogenase is the key enzyme in the de novo biosynthesis of guanosine monophosphate (GMP) (Salomon et al. 2004) . In this study, expression of an octaprenyl diphosphate synthase (ispB, VpaChn25_RS01730) in the terpenoid backbone biosynthesis was increased by 2.392-fold, which was essential for the growth of E. coli (Okada et al. 1997; Helfrich et al. 2019) . Additionally, the DEG encoding a 4-(cytidine 5′-diphospho)-2-C-methyl-d-erythritol kinase (VpaChn25_ RS03685, 2.642-fold) was also significantly up-regulated. This enzyme is essential for isoprenoid biosynthesis in some pathogenic microorganisms (Wada et al. 2003) . The interaction between VpaChn25_0734 (543-bp) and these DEGs should be further investigated in future research.

In the present study, we characterized for the first time the prophage-encoded gene VpaChn25_0734 (543-bp) in V. parahaemolyticus CHN25 genome. The ΔVpaChn25_0734 (543-bp) mutant was obtained and the revertant ΔVpaChn25_0734-com (543-bp) was also constructed. The ΔVpaChn25_0734 (543-bp) mutant was defective in growth and swimming mobility particularly at lower temperatures, as well as in tolerance at pH 7.0-8.5. In the absence of the VpaChn25_0734 (543-bp) gene, cell surface hydrophobicity and biofilm formation at all stages of V. parahaemolyticus CHN25 were significantly decreased (p < 0.05). These results indicated that the prophage-encoded VpaChn25_0734 (543-bp) gene amplified environmental persistence of V. parahaemolyticus CHN25.

Comparative secretomic analysis revealed a slightly increased extracellular proteins secreted by the ΔVpaChn25_0734 (543-bp) mutant, when compared with the WT and ΔVpaChn25_0734-com (543-bp) strains. Based on the Caco-2 cell model in vitro, the deletion of VpaChn25_0734 (543-bp) gene significantly reduced the cytotoxicity of V. parahaemolyticus CHN25 to human intestinal epithelial cells (p < 0.05), which indicated that the VpaChn25_0734 (543-bp) gene enhanced V. parahaemolyticus CHN25 fitness for surviving in the host.

Comparative transcriptome analysis revealed 13 significantly changed metabolic pathways in the ΔVpaChn25_0734 (543-bp) mutant, showing significantly down-regulated carbon source transport and/or utilization, biofilm formation, and T2SS (p < 0.05), consistent with its defective phenotypes. Overall, the prophage-encoded VpaChn25_0734 (543bp) gene amplified ecological persistence of V. parahaemolyticus CHN25. The results in this study facilitate better understanding of pathogenicity and evolution of V. parahaemolyticus, a leading sea foodborne pathogen worldwide.

The online version contains supplementary material available at https:// doi. org/ 10. 1007/ s00294-022-01229-z.

Author contributions Y.X., major experiments, data analysis, and writing-original draft; L.Y., construction of the ΔVpaChn25_0734-com (543-bp) revertant, phenotype and RT-PCR analysis; Y.W., construction of the ΔVpaChn25_0734 (543-bp) mutant; Z.Z., secretome analysis; J.Y. supervision of the in vitro cell mode assay; S.Q. supervision, and discussion; and L.C., funding acquisition, conceptualization, and writing-review and editing. All authors have read and agreed to the published version of the manuscript.

Funding This study was supported by National Natural Science Foundation of China, Grant number 31671946; and Science and Technology Commission of Shanghai Municipal, Grant number 17050502200.

Data availability A complete list of the DEGs is available in the NCBI SRA database (http:// www. ncbi. nlm. nih. gov/ sra/) under the accession number PRJNA733855.

The authors declare no conflict of interest.

Ethics approval Not applicable.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Control of the alternative sigma factor sigmaE in Escherichia coli

Diversity of Vibrio parahaemolyticus in marine fishes of Bangladesh

Phylogenetic relationships and classification of thiolases and thiolase-like proteins of Mycobacterium tuberculosis and Mycobacterium smegmatis

Anchoring the type VI secretion system to the peptidoglycan: TssL, TagL

Vibrio spp. infections. Nat Rev Dis Primers 4:8

Structure, function and regulation of the thermostable direct hemolysin (TDH) in pandemic Vibrio parahaemolyticus

Synergistic stimulation of EpsE ATP hydrolysis by EpsL and acidic phospholipids

Characterization of vB_VpaP_ MGD2, a newly isolated bacteriophage with biocontrol potential against multidrug-resistant Vibrio parahaemolyticus

Exploring the genomic traits of non-toxigenic Vibrio parahaemolyticus strains isolated in Southern Chile

Coordinating assembly of a bacterial macromolecular machine

Minor pseudopilin self-assembly primes type II secretion pseudopilus elongation

Secretion systems in Gram-negative bacteria: structural and mechanistic insights

Transcription inhibition by the depsipeptide antibiotic salinamide A. Elife 3:e02451

The XcpV/GspI pseudopilin has a central role in the assembly of a quaternary complex within the T2SS pseudopilus

Antibiotic resistance of Vibrio parahaemolyticus and Vibrio vulnificus in various countries: a review

Strategies to block bacterial pathogenesis by interference with motility and chemotaxis

The release of a distinct cell type from swarm colonies facilitates dissemination of Vibrio parahaemolyticus in the environment

On the bacteriological examination of shirasu-food poisoning

Sophisticated regulation of transcriptional factors by the bacterial phosphoenolpyruvate: sugar phosphotransferase system

Looking for the hidden: characterization of iysogenic phages in potential pathogenic Vibrio species from the north sea

Genomic characterization of filamentous phage vB_VpaI_VP-3218, an inducible prophage of Vibrio parahaemolyticus

Vibrio parahaemolyticus: a review on distribution, pathogenesis, virulence determinants and epidemiology

Human catalase: looking for complete identity

Comparative secretomics reveals novel virulence-associated factors of Vibrio parahaemolyticus

Bacterial terpene biosynthesis: challenges and opportunities for pathway engineering

Rapid bacteria capturing and killing by AgNPs/N-CD @ ZnO hybrids strengthened photo-responsive xerogel for rapid healing of bacteria-infected wounds

Pyruvate kinase: function, regulation and role in cancer

Immunogenic characterization and protective efficacy of recombinant CsgA, major subunit of curli fibers, against Vibrio parahaemolyticus

Motility of Vibrio spp.: regulation and controlling strategies

Architecture, function, and substrates of the Type II secretion system

MEGA X: molecular evolutionary genetics analysis across computing platforms

Trh (tdh -/trh + ) gene analysis of clinical, environmental and food isolates of Vibrio parahaemolyticus as a tool for investigating pathogenicity

Horizontal transfer of antibiotic resistance genes in clinical environments

Vibrio parahaemolyticus: a review on the pathogenesis, prevalence, and advance molecular identification techniques

Foodborne Vibrio parahaemolyticus in China: prevalence, antibiotic susceptibility, and genetic characterization

Complete genome sequence of Klebsiella pneumoniae subsp. pneumoniae HS11286, a multidrug-resistant strain isolated from human sputum

CDD/SPARCLE: the conserved domain database in 2020

Characterisation of Vibrio species from surface and drinking water sources and assessment of biocontrol potentials of their bacteriophages

Export of a Vibrio parahaemolyticus toxin by the Sec and type III secretion machineries in tandem

Genetic and virulence characterisation of Vibrio parahaemolyticus isolated from Indian coast

ATP-binding cassette transporters in Escherichia coli

Type II secretion system: a magic beanstalk or a protein escalator

Induction and genome analysis of HY01, a newly reported prophage from an emerging shrimp pathogen Vibrio campbellii

The ispB gene encoding octaprenyl diphosphate synthase is essential for growth of Escherichia coli

Malate synthase gene expression during fruit ripening of Cavendish banana (Musa acuminata cv. Williams)

4-genomics of Candida albicans

Genotypic diversity and pathogenic potential of clinical and environmental Vibrio parahaemolyticus isolates from Brazil

Crystal structure of the IIB subunit of a fructose permease (IIBLev) from Bacillus subtilis

Vibrio parahaemolyticus infection induces modulation of IL-8 secretion through dual pathway via VP1680 in Caco-2 cells

The mosaic accessory gene structures of the SXT/R391-like integrative and conjugative elements derived from Vibrio spp. isolated from aquatic products and environment in the Yangtze River Estuary. China

Virulence, resistance, and genetic diversity of Vibrio parahaemolyticus recovered from commonly consumed aquatic products in Shanghai. China

Insights into Vibrio parahaemolyticus CHN25 response to artificial gastric fluid stress by transcriptomic analysis

Isolation and characterization of six Vibrio parahaemolyticus lytic bacteriophages from seafood samples

Crystal structure of 4-(cytidine 5′-diphospho)-2-C-methyl-d-erythritol kinase, an enzyme in the non-mevalonate pathway of isoprenoid synthesis

Global burden of childhood pneumonia and diarrhoea

Etiological, epidemiological, and clinical features of acute diarrhea in China

Pseudomonas aeruginosa T6SS-mediated molybdate transport contributes to bacterial competition during anaerobiosis

Engineering the permeability of Halomonas bluephagenesis enhanced its chassis properties

Comparative proteome analysis of the extracellular proteins of two Haemophilus parasuis strains Nagasaki and SW114

Fitness benefits to bacteria of carrying prophages and prophage-encoded antibiotic-resistance genes peak in different environments

Comparative proteome analysis of secreted proteins of Streptococcus suis serotype 9 isolates from diseased and healthy pigs

Prophage-related gene VpaChn25_0724 contributes to cell membrane integrity and growth of Vibrio parahaemolyticus CHN25

VgrG2 of type VI secretion system 2 of Vibrio parahaemolyticus induces autophagy in macrophages

Genomic and transcriptomic analyses reveal distinct biological functions for cold shock proteins (VpaCspA and VpaCspD) in Vibrio parahaemolyticus CHN25 during low-temperature survival

Comparative proteomics and secretomics revealed virulence and antibiotic resistance-associated factors in Vibrio parahaemolyticus recovered from commonly consumed aquatic products

Feng Z (2020) Pathway construction and metabolic engineering for fermentative production of β-alanine in Escherichia coli

Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations