Transcriptome and metabolome profiling to elucidate the mechanism underlying the poor growth of Streptococcus suis serotype 2 after orphan response regulator CovR deletion

The deletion of orphan response regulator CovR reduces the growth rate of Streptococcus suis serotype 2 (S. suis 2). In this study, metabolome and transcriptome profiling were performed to study the mechanisms underlying the poor growth of S. suis 2 caused by the deletion of orphan response regulator CovR. By comparing S. suis 2 (ΔcovR) and S. suis 2 (SC19), 146 differentially accumulated metabolites (upregulated: 83 and downregulated: 63) and 141 differentially expressed genes (upregulated: 86 and downregulated: 55) were identified. Metabolome and functional annotation analysis revealed that the growth of ΔcovR was inhibited by the imbalance aminoacyl tRNA biosynthesis (the low contents of L-lysine, L-aspartic acid, L-glutamine, and L-glutamic acid, and the high content of L-methionine). These results provide a new insight into the underlying poor growth of S. suis 2 caused by the deletion of orphan response regulator CovR. Metabolites and candidate genes regulated by the orphan response regulator CovR and involved in the growth of S. suis 2 were reported in this study.


Introduction
Streptococcus suis (S. suis) is one of the main pathogenic bacteria on swine farms around the world (1).S. suis is also an emerging zoonotic pathogen that can cause a wide range of diseases (e.g., septicemia, endocarditis, arthritis, meningitis, and pneumonia) (2,3).S. suis mainly causes death in 5-to 10-week-old post-weaned piglets, resulting in significant economic losses to the pig farming industry (4).According to the differences in capsular polysaccharide (CPS) antigenicity, 35 serotypes of S. suis strains were found and defined (5).Among these 35 serotypes, the S. suis 2 strain not only has the highest isolation rate but also poses a serious threat to the development of pig farming and human health (6).
Two-component systems (TCSs) play a crucial role in bacterial infection and response to changes in the external living environment, which have been identified in the vast majority of bacteria (7).The vast majority of TCSs in bacteria consist of a histidine kinase (HK) that receives stimulus signals and a cognate response regulator (RR) that transfers stimulus signals.Once HK receives external stimulus signals, a histidine residue located on the HK will be phosphorylated, and then the phosphoryl group is transferred to the aspartic acid residue on the RR.Subsequently, the phosphorylated RR will interact with the promoter areas of the target genes associated with the response to the environmental stimuli (8).TCSs can determine the survival and pathogenicity of bacteria by regulating various cellular processes, including competition, metabolism, antibiotic resistance, biofilm formation, pathogenicity, and stress response.Therefore, TCSs are viewed as attractive candidates for the development of new strategies for preventing bacterial infections (7,9).
The CovR/S system consists of a sensor kinase gene (covS) and a RR gene (covR), which directly or indirectly regulate 15% of genes connected with bacterial pathogenesis in Streptococcus spp.(10)(11)(12).Due to the impact of CovR/S on a variety of critical virulence factors, it is considered the most essential TCS system for group A Streptococcus (GAS, Streptococcus pyogenes) (13).The CovR that can regulate the virulence of bacteria has been reported in many studies.In Streptococcus agalactia, for example, the covR strain has stronger hemolytic activity, adherence ability to epithelial cells, and higher lethality in the neonate rat sepsis model (14).Previous studies reported that CovR can bind to the promoter areas of the majority of virulence factor genes in GAS (13).These reports suggested that CovR can negatively regulate bacterial virulence.In bacteria, protein tyrosine kinases are linked to pathogenicity and exopolysaccharide synthesis (15).In Edwardsiella piscicida, adenylosuccinate synthase has been reported to have an inhibition effect on the host's nuclear factor-kappa B (NF-κB) signaling pathway, which is a common and important immune signaling pathway (16).Hemolysin is the key virulence component of S. suis 2, which creates holes in the target cell membrane and activates the inflammasome NLRP3 in order to release cytoplasmic K + efflux and cause Streptococcal toxic shock-like syndrome (17).Initiation factors IF1, IF2, and IF3 allow Shine-Dalgarno (SD) sequence-lead mRNAs initiate transcription (18).Ser/Thr kinase (STK) controls a number of cellular activities, including stress response, biofilm formation, membrane production, sporulation, metabolism, pathogenicity, and growth processes in bacteria.Recently, STK has been reported to be related to the virulence of S. suis (19).However, it is still unclear how the CovR/S system regulates genes and metabolites to promote the survivability of S. suis 2.
In recent years, transcriptomics and metabolomics have been extensively utilized to explain the survival and pathogenic mechanism of pathogenic bacteria.To reveal the mechanism by which the orphan RR factor CovR regulates the survivability of S. suis 2, liquid chromatography-tandem mass spectrometry (LC-MS) was employed to dig for the differentially expressed metabolites between S. suis 2 (wild type) and covR (mutant) strains.Furthermore, the assessment of the differentially expressed genes (DEGs) was carried out using RNA-Seq and qRT-PCR.This study will not only reveal the genes and metabolites that are regulated by the orphan RR factor CovR involved in the survivability of S. suis 2 but also provide some potential targets for preventing and controlling S. suis 2 infections.

Bacterial strains and growth condition
The S. suis 2 SC19 strain was donated by Prof. Chen Tan of Huazhong Agricultural University.SC19 is plated on Tryptose Soya agar (TSA) or cultured in Tryptose Soya broth (TSB) (Difco Laboratories, Detroit, MI, USA) with 10% (v/v) newborn bovine serum (Sijiqing Biological Engineering Materials Co., Ltd., Hangzhou, China) at 37 • C. covR mutant and complemented C covR strains are stored in our laboratory and cultured as described above (20).

Growth characteristics
The growth characteristics of wild-type SC19, covR mutant, and complemented C covR strains were determined by transferring overnight bacterial cultures into fresh preheated TSB supplemented with 10% newborn bovine serum at a ratio of 1:1000.The bacteria were grown at 37 • C, accompanied by shaking at 180 rpm for 12 h, and OD 600 was determined by an Automatic Microbial Growth Curve Analyzer (MGC-200, Ningbo Scientz Biotechnology Co., Ltd., Ningbo, China) for every hour.One-way ANOVA test was used to verify the significance of growth characteristics of SC19, covR, and C covR for every hour.

RNA extraction, library construction, and sequencing
The total RNA was extracted using the TRIzol reagent (Takara, Dalian, China) from SC19 (n = 3) and covR (n = 3) as described by Rio et al. (21), and the genomic DNA in the total RNA was removed by using DNase I (Takara).A total of 2 µg of RNA was used to build the RNA-seq transcriptome library by using the TruSeq TM RNA sample preparation Kit (Illumina, San Diego, CA).Briefly, the Ribo-Zero Magnetic Kit (epicenter) was used to exhaust ribosomal RNA (rRNA), and then the fragmented mRNA was broken into short fragments (200 bp) by using a fragmentation buffer.The mRNA fragment was performed by RT-PCR using the SuperScript Double-Stranded cDNA Synthesis Kit (Invitrogen, CA) with random hexamer primers (Illumina).The dTTP was replaced with dUTP when the second strand was synthesized.Then, after end-repair and phosphorylation, the synthesized cDNA was treated with "A" base addition according to the protocol of Illumina's introduction.The UNG enzyme was used to recognize and degrade those secondstrand cDNA with dUTP.The library was constructed using 200-bp cDNA target fragments of 2% low-range ultra-agarose.Subsequently, these fragments were amplified by performing 15 PCR cycles with Phusion DNA polymerase and quantified using TBS380.Finally, the paired-end RNA-seq sequencing library was constructed using Illumina Novaseq (2 × 150 bp read length).

Transcriptome data analysis
The raw data collected from Illumina Hiseq sequencing were converted into raw data using base calling, and the results were saved in FASTQ file format.For mapping assembly sequences and distinct genes to reference sequences (GenBank: CP020863.1),Bowtie with default parameters (http://bowtiebio.sourceforge.net/index.shtml) was employed.FPKM (fragments per kilobase million) was used to quantify the quantity of genes and transcripts (http://deWeylab.biostat.wisc.edu/rsem/,https:// ccb.jhu.edu/software/stringtie/). The transcriptome data were analyzed by the cloud platform of Shanghai Majorbio Bio-Pharm Technology Co., Ltd.(Shanghai, China).Additionally, the DEGs between covR and S. suis 2 strains were detected using the DESeq2 software.To control the probability of error with multiple test corrections (Benjamini/Hochberg, BH), the p-values were adjusted.DEGs were identified as their fold change ≥1.2 or ≤0.8 and adjusted p-values <0.05.The GO analysis for DEGs was performed using the Goatools software, while the KEGG analysis was performed using R script.

Extraction of metabolites from bacteria and quality control measurement
Overnight, SC19 (n = 4) and covR (n = 4) cultures were diluted at a 1:1000 ratio into TSB supplemented with 10% newborn bovine serum, and the mixture was then incubated to the logarithm phase.A total of 100 ml of the bacterium was transferred to a clean tube and centrifuged at 12,000g at 4 • C for 10 min, and then the supernatant was removed.The metabolites were extracted from the bacterium (SC19 and covR) by a 400-µl isovolumetric methanol and acetonitrile solution and then sonicated by an ultrasonic cleaner (SBL-10TD, Ningbo Xinzhi Biotechnology Co., Ltd.) at 40 kHz for 30 min at 5 • C.After being placed at −20 • C for 30 min, the mixture was centrifuged at 13,000g at 4 • C for 15 min.The supernatant was carefully transferred to a new microtube and evaporated to dry under a gentle stream of nitrogen.Then, the samples were resuspended with 100 µl of loading isovolumetric acetonitrile and water solution by short sonication using an ultrasonic cleaner (Ningbo Scientz Biotechnology Co., Ltd.) in a 5 • C water bath.The extracted metabolites were centrifuged at 13,000g for 15 min at 4 • C and then UHPLC-MS/MS analysis was performed.The same volumes of all samples were mixed to be used as a pooled quality control.

UHPLC-MS/MS analysis
The LC-MS System (UHPLC-Q Exactive HF-X System, Thermo Fisher Scientific, Shanghai, China) was used to analyze the extracted metabolites.First, the sample was injected into a Waters HSS T3 column (100 mm × 2.1 mm i.d., 1.8 µm), which was followed by mass spectrometry detection.The mobile phases consisted of two solvents: solvent A, which consisted of 95% 0.1% formic acid and 5% acetonitrile, and solvent B, which consisted of 47.5% 0.1% formic acid in acetonitrile with 47.5% isopropanol and 5% water.The solvent gradient employed the following conditions: First, the gradient changed from 0 B to 24.5% B at a flow rate of 0.4 ml/min for 3.5 min; then, the gradient changed from 24.5 B to 65% B for 1.5 min; the gradient changed from 65 B to 100% B for 0.5 min; the gradient remained at 100% B for 1.9 min with the flow rate being increased from 0.4 to 0.6 ml/min; the gradient changed from 100 B to 51.5% B for 0.2 min; the gradient changed from 51.5 B to 0% B for 0.2 min at a flow rate of 0.6 to 0.5 ml/min; the gradient remained at 0% B for 0.2 min at a flow rate of 0.5 to 0.4 ml/min; the system equilibration was performed as the gradient remained at 0% B for 1 min at a flow rate of 0.4 ml/min.The sample injection volume was 2 µl at a flow rate of 0.4 ml/min.The column temperature was maintained at 40 • C throughout the analysis, and the samples were stored at 4 • C. The mass spectrometric data were collected using continuous scanning in either positive or negative ion modes.Data acquisition was performed in the data-dependent acquisition mode.The full MS resolution was set at 60,000, and the MS/MS resolution was set at 7,500.The heater temperature and capillary temperature were set at 425 • C and 325 • C, respectively.The sheath gas flow rate was set at 50 arb and the aux gas was set at 13 arb.The ion-spray voltage floating was −3,500 V in negative mode and 3,500 V in positive mode.Normalized collision energy for MS/MS was set at 20-40-60 V rolling.The detection range for mass spectrometry was from 70 to 1,050 m/z.

Metabolite data and di erential metabolites analysis
A three-dimensional data matrix that included sample information, metabolite name, and mass spectral response intensity in CSV format was obtained by analyzing metabolite data with Progenesis QI (Waters Corporation, Milford, USA).Internal standard peaks, as well as any known false positive peaks that include noise, column bleed, and derivatized reagent peaks, were removed from the data matrix, which was then de-redundant and peak pooled.The metabolic mass spectra features were recognized by accurate mass, MS/MS fragment spectra, and isotope ratio differences were searched in reliable biochemical databases such as the Human Metabolome Database (HMDB) (http://www.hmdb.ca/) and the METLIN database (https://metlin.scripps.edu/).
After searching in the database, these data were uploaded to the Majorbio cloud platform (https://cloud.majorbio.com)for analysis.In each set of samples, at least 80% of the metabolic features detected were retained.The metabolic levels in specific samples that fell below the lower limit of quantitation were imputed with minimum metabolite values, and each metabolic feature was normalized by sum.The response intensity of the sample mass spectrum peaks was normalized using the sum normalization method to obtain the normalized data matrix.Those QC samples with variables that had a relative standard deviation >30% were removed, and a logarithmic transformation (log10) was applied to the final data matrix for subsequent analysis.
Next, principal component analysis (PCA) and orthogonal least squares partial discriminant analysis (OPLS-DA) were conducted on the preprocessed data array using the R package (Version 1.6.2).The stability of the model was evaluated through seven cycles of cross-validation.To identify significantly different metabolites, VIP values obtained from the OPLS-DA model and Student's t-test were used, with VIP >1 and p < 0.05 indicating significance.Metabolic pathway annotation was carried out using the KEGG data (https://www.kegg.jp/kegg/pathway.html) to determine the pathways with which the differential metabolites were associated.Pathway enrichment analysis was performed using the Python package scipy.stats,and Fisher's exact test was used to identify the most relevant biological pathways to the experimental treatments.The raw data were finally uploaded to MetaboLights, where they were given the unique identifier MTBLS8495 (www.ebi.ac.uk/metabolights/MTBLS8495).

Quantitative real-time polymerase chain reaction
In this study, the accuracy of RNA-seq data was validated by performing qRT-PCR on the selected DEGs.Among them, B9H01_RS09035 is glutamine ABC transporter substrate-binding protein; B9H01_RS06585 is glutamine ABC transporter permease; B9H01_RS06580 is glutamine ABC transporter permease; B9H01_RS05040 is peptide ABC transporter ATP-binding protein; B9H01_RS05035 is amino acid ABC transporter permease; B9H01_RS09790 is ABC transporter ATP-binding protein; and B9H01_RS01490, B9H01_RS06570, and B9H01_RS09040 are ABC transporters.To begin, RNA extraction was followed as described in the preceding part of the text, and 1 µg of RNA was reserved for transcription using a PrimeScript R RT reagent Kit and gDNA Eraser (Takara, Dalian, China).Next, a 7500 Fast Real-Time PCR System (Applied Biosystems, Foster City, CA, USA) was used to perform the PCR reaction using Taq Pro Universal SYBR qPCR Master Mix (Q712, Vazyme, Nanjing, China).The 16S RNA gene served as the control gene, and all mRNA transcription levels were reported as 2 − Ct .The primers in Table 1 were designed according to the genomic sequence of SC19 and used for qRT-PCR.

Statistical analyses
The data in this study were analyzed using SPSS (Version 17.0, SPSS Inc., Chicago, IL, USA) and GraphPad Prism 8. Unpaired t-tests between two groups or one-way ANOVA among multiple groups were used to obtain p-values ( * p < 0.05, * * p < 0.01, * * * p < 0.001; NS, not significant).

Results
The e ects of covR deletion on the growth characteristics of SC , covR, and c covR strains As shown in Figure 1, the growth of covR was significantly slower compared to that of SC19 and the complemented C covR strains after 3 h.

Metabolome profiling
For the purpose of determining the poor growth of mutant covR, metabolome profiling via untargeted LC-MS was performed on SC19 and mutant covR.The metabolites identified from LC-MS were presented in Supplementary Table S1.There was a total of 514 metabolites identified, which were classified into 11 categories.Supplementary Figure S1 provides the PCA of SC19, covR, and QC samples, which present significant separation.The coloration of each sample was analyzed using metabolite concentration data and shown in Supplementary Figure S1, which suggests that different samples were distinguished significantly.This means that the metabolome data were highly reliable.These results also suggested that SC19 and covR samples had substantially different metabolite profiles.

DAMs between SC and mutant covR strains
The uniquely mapped read ratio of all samples was higher than 96%, which means the reference sequence was reliable.To acquire Differentially accumulated metabolites (DAMs) between SC19 and mutant covR strains, parameters such as p-values <0.05, |Log2 Fold Change| >1, and variable importance in the projection (VIP) >1 were used to analyze all metabolite data.In total, 146 DAMs (83 were upregulated and 63 were downregulated) were acquired (Figures 2A, B, Supplementary Table S2).The results of the Kyoto Encyclopedia of Genes and Genomes (KEGG, Figure 2C) analysis showed that arginine biosynthesis (map00220), pyrimidine metabolism (map00240), aminoacyl-tRNA biosynthesis (map00970), alanine, aspartate and glutamate metabolism (map00250), ABC transporters (map02010), and lysine biosynthesis (map00300) were significantly enriched.Further analysis results showed that DMAs in these enrichment pathways were likely involved in aminoacyl-tRNA biosynthesis, energy metabolism, and amino acid biosynthesis and metabolism (Table 2).In addition, the 146 DMAs could be categorized into >10 classes (Figure 2D).These analysis results suggest that DAMs involved in aminoacyl-tRNA biosynthesis, energy metabolism, and amino acid biosynthesis and metabolism may play a key role in slowing down the growth rate of the mutant covR strain.

DEGs between SC and mutant covR strains
To explore genes responsible for poor growth of covR strains, DEGs (|Log2 Fold Change| >1.2 or <0.8) were identified between SC19 and mutant covR strains.A total of 141 DEGs (86 were upregulated and 55 were downregulated) were identified.The GO annotation analysis result showed that the 141 DEGs were classified into three categories (Figure 3B), among which the integral component of membrane (36, 25.5%) was the largest group; the second largest group was ATP binding (26, 18.4%); only a few were related to metal ion binding (3, 2.1%) and iron-sulfur cluster assembly (3, 2.1%).As shown in Figure 3C, the enrichment pathways in KEGG terms were ABC transporters, glycerolipid metabolism, propanoate metabolism, carbon fixation pathways in prokaryotes, valine, leucine, and isoleucine degradation, and glycosaminoglycan degradation (Supplementary Table S3).

The correlation analysis of DAMs and DEGs
To validate the accuracy of the identified DEGs by RNA-seq, nine DEGs were validated by qRT-PCR (Figure 4).Consistent with the RNA-Seq results (Supplementary Table S3), qRT-PCR results showed that the expression levels of the identified genes in the mutant covR strain were significantly increased compared with those in the SC19 strain.Further research showed that all the gathered data from metabolome and transcriptome profiles had high reliability according to the correlation analysis (Figure 5).We could conclude that variations in metabolite accumulation between SC19 and mutant covR strains were closely regulated by different gene expressions.

Discussion
Prior studies have reported the importance of the TCS CovR in S. suis 2 virulence and growth characteristics (14,(22)(23)(24)(25).However, the mechanism by which CovR deletion affects the growth of S. suis 2 is not yet clear.CovR is essential for the growth of S. suis First, we constructed the covR deletion mutant strain.We that the mutant covR grew significantly more slowly than the parental strain SC19 on the growth curve (Figure 1A) and it also growth slow on the TSA plate (Figure 1B).This phenomenon was also presented by Pan et al. (22).To further understand the mechanism of poor growth of the covR mutant, transcriptome and metabolome profiling were performed.the transcriptome profile showed that there were no significant changes in gene expression levels of acetyl-CoA and oxalosuccinic acid between mutant covR and SC19 strains.The reason for the up-accumulation of acetyl-CoA and oxalosuccinic acid may be due to less consumption rather than excessive production.The D-4'-Phosphopantothenate was down-accumulated in metabolism and down-expression in the transcriptome, which means there were fewer precursors of coenzyme A, but the level of acetyl-CoA was still up-accumulated.Less consumption is the main reason for the up-accumulation of acetyl-CoA and oxaloacetylsuccinic acid in the mutant covR strain.Therefore, the mutant consumed less acetyl-CoA compared to SC19 and caused less energy supply to the biosynthesis of complex macromolecules, which resulted in poor growth.
The poor growth of mutant covR was related to amino acid biosynthesis and metabolism Amino acid synthesis and metabolism play a key role in bacterial growth.There are over 300 naturally occurring amino acids, and only 22 constitute the monomer units of proteins, which can be classified as acidic or basic and hydrophilic or hydrophobic according to their respective side chain R (29).1-aminocyclopropanecarboxylic acid (ACC) is a cyclopropanecontaining non-proteinogenic amino acid isolated from various fruits and plant tissues and is considered the immediate precursor of the plant hormone ethylene (30).Recently, ACC has been found to perform similar synthesis machinery to plants in microorganisms (31,32).Citrulline and ornithine, as we all know, are the intermediate substances of the urea cycle and have recently been reported to be involved in the metabolism of arginine in different bacteria (33).N-acetylglutamate in prokaryotes is also the first intermediate in arginine biosynthesis (34).As we all know, the level of amino acids can influence the amount of aminoacyl-tRNA.Several reports have shown that, in E. coli and S. typhimurium, the levels of aminoacyl-tRNA synthetases change with the rate of growth in different media (35).Our study found that seven metabolites related to amino acid biosynthesis and metabolism were differently accumulated, including two up-accumulated (1aminocyclopropanecarboxylic acid and allysine) and five downaccumulated (citrulline, L-ornithine, N-acetylaspartate, N-acetyl-L-aspartic acid, and N-acetylglutamic acid) in mutant covR compared to S. suis 2 SC19.The down-accumulation of citrulline, ornithine, and N-acetylglutamate were all involved in arginine metabolism, which means that arginine may play a key role in the growth of SC19.
The argF, as a gene code ornithine carbamoyltransferase, was downregulated in mutant covR compared to SC19 in the transcriptome.Ornithine carbamoyltransferase can catalyze the condensation of ornithine with carbamoyl phosphate to form the amino acid citrulline, which is subsequently delivered to the cytosol (36).In bacteria, urea and ornithine produced by the decomposition of arginine can be further utilized as nitrogen and carbon sources (37).In other words, the mutant covR produced less arginine than S. suis 2 SC19, which decreased the availability of nitrogen for the biosynthesis of amino acids.Thus, the CovR deletion influenced the rate of amino acid biosynthesis, resulting in poor growth of S. suis 2.
Poor growth of mutant was related to aminoacyl-tRNA biosynthesis Aminoacyl-tRNA synthetases play a key role in the process of aminoacyl-tRNA biosynthesis, which produces aminoacyl-tRNA, AMP, and PPi by esterifying an amino acid to the 3 ′ end of a tRNA and hydrolyzing one molecule of ATP (38).Previous studies showed that translation elongation, which is a process of the aminoacyl-tRNA translocation of the ribosome, is an essential step in protein synthesis (39), and translation is the key to defining the growth rate (40).Therefore, the quantity of aminoacyl-tRNA will have an impact on the pace of bacterial growth.In our study, the aminoacyl-tRNA biosynthesis related to metabolism (L-methionine, L-lysine, L-aspartic L-glutamate, and L-glutamine) is down-accumulated except L-methionine in the mutant covR strain.Meanwhile, the amino acid ABC transporter substrate-binding protein genes (B9H01_RS09940) and the amino acid ABC transporter permease gene (B9H01_RS09935) were downregulated, and the amino acid ABC transporter ATPbinding protein genes (B9H01_RS06570, B9H01_RS09040, and B9H01_RS05040) were upregulated in mutant covR compared to SC19.The amino acid ABC transporter is a kind of active transporter for amino acid uptake (41).Therefore, we assume that mutant covR has less amino acid ABC transporters that involve in amino acid uptake, thus expressing more ATP-binding proteins to transport more amino acids to maintain homeostasis, but still making the imbalance of aminoacyl tRNA biosynthesis, ultimately leading to slower growth.

Conclusion
The mechanism of slow growth of covR was proposed in Figure 6 can be explained as follows: first, the mutant covR exhibits reduced consumption of acetyl-CoA, resulting in a decrease in ATP production for biosynthesis.This decrease in ATP production, in turn, leads to a shortage of nitrogen resources due to reduced arginine biosynthesis, which consequently slows down the rate of amino acid biosynthesis.Ultimately, the combination of decreased ATP and amino acid production results in a reduced availability of aminoacyl tRNA, leading to poor growth in mutant covR.
Of Hubei Key Laboratory of Animal Nutrition and Feed Science, Wuhan Polytechnic University (202314).

FIGURE
FIGURE Di erentially accumulated metabolites between SC and mutant covR strains.(A) Volcano plots displayed DAMs that showed significant change (upregulated and downregulated) or no significant change between SC and mutant covR strains.(B) The heat map of mutant covR and two SC strains.(C) KEGG enrichment analysis of DAMs.(D) Pie chart depicting the categories of DAMs.

FIGURE
FIGURE DEGs between mutant covR and SC strains.(A) The heat map of mutant covR and SC strains.(B) Go annotation analysis.(C) The KEGG enrichment pathways of DEGs.

FIGURE
FIGURECorrelation analysis of DEGs and DAMs in the mutant covR and SC strains.The levels of correlation of each sample were indicated by the color, from high (red) to low (blue), and the correlation coe cient between samples was represented the Z-score.*p < ., **p < ., ***p < ..

FIGURE
FIGURE Mechanisms of the poor growth of Streptococcus suis serotype after orphan response regulator CovR deletion.(A) Synthesis of aminoacyl tRNAs in SC in covR strain.(B) Synthesis of aminoacyl tRNAs in covR strain.
TABLE The list of DAMs identified between S. suis SC and mutant covR strains.