<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Microbiol.</journal-id>
<journal-title>Frontiers in Microbiology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Microbiol.</abbrev-journal-title>
<issn pub-type="epub">1664-302X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fmicb.2021.632095</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Microbiology</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Integrated Proteomic and Transcriptomic Analyses Reveal the Roles of <italic>Brucella</italic> Homolog of BAX Inhibitor 1 in Cell Division and Membrane Homeostasis of <italic>Brucella suis</italic> S2</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name><surname>Zhang</surname> <given-names>Guangdong</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/998166/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Zhong</surname> <given-names>Fangli</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Chen</surname> <given-names>Lei</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Qin</surname> <given-names>Peipei</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Li</surname> <given-names>Junmei</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Zhi</surname> <given-names>Feijie</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Tian</surname> <given-names>Lulu</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Zhou</surname> <given-names>Dong</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/314290/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Lin</surname> <given-names>Pengfei</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/314338/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Chen</surname> <given-names>Huatao</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/86891/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Tang</surname> <given-names>Keqiong</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Liu</surname> <given-names>Wei</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1041714/overview"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Jin</surname> <given-names>Yaping</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/314417/overview"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Wang</surname> <given-names>Aihua</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<xref ref-type="corresp" rid="c002"><sup>&#x002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1149245/overview"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>College of Veterinary Medicine, Northwest A&#x0026;F University</institution>, <addr-line>Yangling</addr-line>, <country>China</country></aff>
<aff id="aff2"><sup>2</sup><institution>Key Laboratory of Animal Biotechnology, Ministry of Agriculture, Northwest A&#x0026;F University</institution>, <addr-line>Yangling</addr-line>, <country>China</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Marie-Joelle Virolle, Centre National de la Recherche Scientifique (CNRS), France</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Judith Smith, University of Wisconsin&#x2013;Madison, United States; William Doerrler, Louisiana State University, United States</p></fn>
<corresp id="c001">&#x002A;Correspondence: Aihua Wang, <email>aihuawang1966@126.com</email></corresp>
<corresp id="c002">Yaping Jin, <email>yapingjin@163.com</email></corresp>
<fn fn-type="other" id="fn004"><p>This article was submitted to Microbial Physiology and Metabolism, a section of the journal Frontiers in Microbiology</p></fn>
</author-notes>
<pub-date pub-type="epub">
<day>28</day>
<month>01</month>
<year>2021</year>
</pub-date>
<pub-date pub-type="collection">
<year>2021</year>
</pub-date>
<volume>12</volume>
<elocation-id>632095</elocation-id>
<history>
<date date-type="received">
<day>22</day>
<month>11</month>
<year>2020</year>
</date>
<date date-type="accepted">
<day>12</day>
<month>01</month>
<year>2021</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x00A9; 2021 Zhang, Zhong, Chen, Qin, Li, Zhi, Tian, Zhou, Lin, Chen, Tang, Liu, Jin and Wang.</copyright-statement>
<copyright-year>2021</copyright-year>
<copyright-holder>Zhang, Zhong, Chen, Qin, Li, Zhi, Tian, Zhou, Lin, Chen, Tang, Liu, Jin and Wang</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>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.</p></license>
</permissions>
<abstract>
<p>BAX inhibitor 1 (BI-1) is an evolutionarily conserved transmembrane protein first identified in a screening process for human proteins that suppress BAX-induced apoptosis in yeast cells. Eukaryotic BI-1 is a cytoprotective protein that suppresses cell death induced by multiple stimuli in eukaryotes. <italic>Brucella</italic>, the causative agent of brucellosis that threatens public health and animal husbandry, contains a conserved gene that encodes BI-1-like protein. To explore the role of the <italic>Brucella</italic> homolog of BI-1, BrBI, in <italic>Brucella suis</italic> S2, we constructed the <italic>brbI</italic> deletion mutant strain and its complemented strain. <italic>brbI</italic> deletion altered the membrane properties of <italic>Brucella suis</italic> S2 and decreased its resistance to acidic pH, H<sub>2</sub>O<sub>2</sub>, polymyxin B, and lincomycin. Additionally, deleting <italic>brbI</italic> led to defective growth, cell division, and viability in <italic>Brucella suis</italic> S2. We then revealed the effect of <italic>brbI</italic> deletion on the physiological characteristics of <italic>Brucella suis</italic> S2 via integrated transcriptomic and proteomic analyses. The integrated analysis showed that <italic>brbI</italic> deletion significantly affected the expression of multiple genes at the mRNA and/or protein levels. Specifically, the affected divisome proteins, FtsB, FtsI, FtsL, and FtsQ, may be the molecular basis of the impaired cell division of the <italic>brbI</italic> mutant strain, and the extensively affected membrane proteins and transporter-associated proteins were consistent with the phenotype of the membrane properties&#x2019; alterations of the <italic>brbI</italic> mutant strain. In conclusion, our results revealed that BrBI is a bacterial cytoprotective protein involved in membrane homeostasis, cell division, and stress resistance in <italic>Brucella suis</italic> S2.</p>
</abstract>
<kwd-group>
<kwd><italic>Brucella suis</italic> S2</kwd>
<kwd>BAX inhibitor 1 (BI-1)</kwd>
<kwd>cell viability</kwd>
<kwd>cell division</kwd>
<kwd>membrane homeostasis</kwd>
<kwd>stress resistance</kwd>
<kwd>proteomics</kwd>
<kwd>transcriptomics</kwd>
</kwd-group>
<contract-num rid="cn001">2018YFD0500900</contract-num>
<contract-sponsor id="cn001">National Key Research and Development Program of China<named-content content-type="fundref-id">10.13039/501100012166</named-content></contract-sponsor>
<counts>
<fig-count count="11"/>
<table-count count="1"/>
<equation-count count="1"/>
<ref-count count="46"/>
<page-count count="21"/>
<word-count count="0"/>
</counts>
</article-meta>
</front>
<body>
<sec id="S1">
<title>Introduction</title>
<p>BAX inhibitor 1 (BI-1) is an evolutionarily conserved transmembrane protein, first identified in a screening process for human protein candidates that functionally suppress BAX-induced apoptosis in yeast cells (<xref ref-type="bibr" rid="B21">Huckelhoven, 2004</xref>; <xref ref-type="bibr" rid="B19">Henke et al., 2011</xref>). Eukaryotic BI-1 is a cytoprotective protein located on the endoplasmic reticulum (ER) membrane with 6-7 transmembrane helices (<xref ref-type="bibr" rid="B10">Chae et al., 2004</xref>; <xref ref-type="bibr" rid="B6">Bultynck et al., 2012</xref>). In animals, BI-1 exerts its cytoprotective role in the intrinsic apoptotic pathway, which consists of the mitochondria-dependent and ER stress-dependent pathways (<xref ref-type="bibr" rid="B21">Huckelhoven, 2004</xref>). In mitochondria-dependent apoptosis, BI-1 interacts directly with the anti-apoptosis proteins, Bcl-2 and Bcl-XL, to inhibit BAX/BAK translocation to the mitochondria, thus protecting cells from mitochondrial dysfunction-induced apoptosis (<xref ref-type="bibr" rid="B21">Huckelhoven, 2004</xref>; <xref ref-type="bibr" rid="B30">Liu, 2017</xref>). In ER-stress-dependent apoptosis, BI-1 interacts directly with IRE1&#x03B1;, an ER-stress sensor, to reduce BAX and IRE1&#x03B1; binding, thus suppressing ER-stress-induced apoptosis (<xref ref-type="bibr" rid="B10">Chae et al., 2004</xref>; <xref ref-type="bibr" rid="B29">Lisbona et al., 2009</xref>; <xref ref-type="bibr" rid="B40">Shi et al., 2018</xref>). BI-1 also modulates Ca<sup>2+</sup> homeostasis of the ER, relieving the apoptosis triggered by Ca<sup>2+</sup> overload (<xref ref-type="bibr" rid="B6">Bultynck et al., 2012</xref>; <xref ref-type="bibr" rid="B30">Liu, 2017</xref>). Although the mechanism of apoptosis in plants differs from that in animals, plant BI-1 shares a highly conserved structure and function with animal BI-1. In yeast and animal cells, plant BI-1 inhibits mammalian BAX-induced cell death. In plant cells, plant BI-1 suppresses mammalian BAX-induced apoptosis and inhibits cell death triggered by biotic and abiotic stresses such as pathogens, acidic pH, oxidative stress, heat shock, and low nutrition (<xref ref-type="bibr" rid="B45">Watanabe and Lam, 2009</xref>).</p>
<p>Automated computational analysis using gene prediction method revealed that BI-1-like proteins also extensively exist in prokaryotes (<xref ref-type="bibr" rid="B16">Gamboa-Tuz et al., 2018</xref>). However, only two prokaryotic BI-1 homologs have been reported to date; both share amino acid sequence similarity and concordant hydrophobicity profiles with eukaryotic BI-1 (<xref ref-type="bibr" rid="B25">Kihara et al., 1998</xref>; <xref ref-type="bibr" rid="B11">Chang et al., 2014</xref>). <italic>Escherichia coli</italic> (<italic>E. coli</italic>) YccA is the first identified prokaryotic member of BI-1 family proteins, which is also a transmembrane protein located on the inner membrane. Although the precise physiological role of YccA in <italic>E. coli</italic> is poorly understood, <xref ref-type="bibr" rid="B25">Kihara et al. (1998</xref>, <xref ref-type="bibr" rid="B26">1999)</xref> reported that <italic>E. coli</italic> YccA is a proteolytic substrate of FtsH and shares the same recognition site with some other FtsH substrates. FtsH is a membrane-anchored ATP-dependent metalloprotease, conserved in both prokaryotes and eukaryotes (<xref ref-type="bibr" rid="B39">Schumann, 1999</xref>; <xref ref-type="bibr" rid="B27">Langklotz et al., 2012</xref>). In prokaryotes, its main role is to degrade redundant or abnormal membrane proteins as part of a quality control mechanism to maintain membrane homeostasis (<xref ref-type="bibr" rid="B39">Schumann, 1999</xref>; <xref ref-type="bibr" rid="B27">Langklotz et al., 2012</xref>; <xref ref-type="bibr" rid="B3">Bittner et al., 2017</xref>). <xref ref-type="bibr" rid="B44">van Stelten et al. (2009)</xref> revealed that YccA overexpression in <italic>E. coli</italic> can relieve the bacterial cell death triggered by jamming of the SecY translocation system, a translocator complex consisting of SecY, E, and G (<xref ref-type="bibr" rid="B5">Breukink, 2009</xref>). Protein translocation mediated by the SecY translocator is a fundamental process for bacteria (<xref ref-type="bibr" rid="B5">Breukink, 2009</xref>). Notably, SecY is an FtsH substrate that shares the same recognition site as YccA, suggesting a bacterial cytoprotective function for YccA (<xref ref-type="bibr" rid="B25">Kihara et al., 1998</xref>). <italic>Bacillus subtilis</italic> (<italic>B. subtilis</italic>) YetJ is another prokaryotic BI-1 (<xref ref-type="bibr" rid="B11">Chang et al., 2014</xref>). Chang et al. analyzed the crystal structure of <italic>B. subtilis</italic> YetJ and characterized it as a Ca<sup>2+</sup> homeostasis maintainer, implying that YetJ shares molecularly similar functions to those of eukaryotic BI-1 (<xref ref-type="bibr" rid="B11">Chang et al., 2014</xref>). To date, relevant studies on prokaryotic BI-1 remain rare, but suggest that prokaryotic BI-1 may be important for maintaining bacterial cell envelope properties as well as bacterial cell viability. As expected, <italic>Brucella</italic> spp. contain a eukaryotic BI-1 homolog. In the <italic>Brucella suis</italic> S2 (<italic>Brucella suis</italic> bv. 1 str. S2 [<italic>B. suis</italic> S2]) genome, <italic>BSS2_RS00410</italic> is annotated as a gene encoding a BI-1-like protein that remains undefined. Because this BI-1-like protein is 100% conserved in all <italic>Brucella</italic> species, we hereafter refer to this <italic><underline>Br</underline>ucella</italic> homologue of <underline>BI</underline>-1 as BrBI.</p>
<p><italic>Brucella</italic> spp. is Gram-negative, facultative intracellular bacteria belonging to the class Alphaproteobacteria, which can invade various host cells, including macrophages and trophoblasts (<xref ref-type="bibr" rid="B7">Byndloss and Tsolis, 2016</xref>). Once inside the host cell, <italic>Brucella</italic> initially suffers harsh intracellular environmental conditions such as low nutrition, acidic pH, and oxidative stress (<xref ref-type="bibr" rid="B17">Gomez et al., 2013</xref>; <xref ref-type="bibr" rid="B8">Celli, 2019</xref>). As a successful pathogen, <italic>Brucella</italic> has evolved numerous strategies, such as the type IV secretory system (T4SS) encoded by the <italic>virB</italic> operon, the quorum-sensing system, two-component systems, and modified lipopolysaccharide (LPS), to disturb the host immune system and adapt to intracellular environments (<xref ref-type="bibr" rid="B4">Boschiroli et al., 2002</xref>; <xref ref-type="bibr" rid="B9">Celli et al., 2003</xref>; <xref ref-type="bibr" rid="B43">Uzureau et al., 2007</xref>; <xref ref-type="bibr" rid="B34">Martinez-Nunez et al., 2010</xref>; <xref ref-type="bibr" rid="B23">Ke et al., 2015</xref>; <xref ref-type="bibr" rid="B7">Byndloss and Tsolis, 2016</xref>; <xref ref-type="bibr" rid="B1">Altamirano-Silva et al., 2018</xref>). Owing to these stealthy strategies, a fraction of internalized <italic>Brucella</italic> survives and multiplies in the host cell, eventually leading to long-term chronic infection. Brucellosis, caused by <italic>Brucella</italic> spp., is a severe globally distributed zoonosis characterized by abortion and infertility in animals and debilitating chronic infections in humans, which threatens public health and sustainable development of animal husbandry worldwide (<xref ref-type="bibr" rid="B7">Byndloss and Tsolis, 2016</xref>; <xref ref-type="bibr" rid="B8">Celli, 2019</xref>).</p>
<p>Given the membrane-associated role and the bacterial cytoprotective potential of known prokaryotic BI-1, we hypothesized that BrBI may have roles in modulating <italic>B. suis</italic> S2 viability and/or membrane properties. Consequently, we investigated the roles of BrBI in <italic>B. suis</italic> S2 stress tolerance, antibiotic resistance, cell viability, and cell division <italic>via</italic> integrated proteomic and transcriptomic analyses.</p>
</sec>
<sec id="S2" sec-type="materials|methods">
<title>Materials and Methods</title>
<sec id="S2.SS1">
<title>Biosafety Statement</title>
<p>All experiments were performed in accordance with the &#x201C;Regulations on biosafety of Pathogenic Microorganism Laboratory&#x201D; (2004) No. 424 set by the State Council of the People&#x2019;s Republic of China and approved by Biosafety Committee of Northwest A&#x0026;F University.</p>
</sec>
<sec id="S2.SS2">
<title>Bacterial Strains and Cultures</title>
<p>The bacterial strains used in this study were derived from the <italic>B. suis</italic> vaccine strain S2 (<xref ref-type="bibr" rid="B15">Di et al., 2016</xref>). The strains were cultured in tryptic soy broth (TSB) with shaking or on tryptic soy agar (TSA) at 37&#x00B0;C. To determine the concentrations of viable <italic>B. suis</italic> S2 strains, single bacterial colonies from freshly streaked TSA plates were inoculated into 50 ml TSB and grown at 37&#x00B0;C for 2 days with gentle shaking. Subsequently, the cultures were diluted 100-fold with TSB and incubated at 37&#x00B0;C to exponential phase. The cultures were harvested and resuspended in sterile phosphate-buffered saline (PBS), then serially diluted in sterile PBS and plated onto TSA plates to confirm the concentrations of viable bacteria. Antibiotics, when required, were added at the following concentrations: kanamycin, 50 &#x03BC;g/ml; ampicillin, 100 &#x03BC;g/ml. All work involving live <italic>B. suis</italic> S2 strains was performed in a biosafety level 3 (BSL-3) facility.</p>
</sec>
<sec id="S2.SS3">
<title>Generation of the <italic>brbI</italic> Deletion Mutant Strain and Its Complemented Strain</title>
<p>To obtain the <italic>brbI</italic> deletion mutant strain, a resistance gene replacement procedure was carried out. Briefly, the upstream homologous arm of <italic>brbI</italic> was amplified using the UP-F and UP-R primers, while the downstream homologous arm of <italic>brbI</italic> was amplified using the DW-F and DW-R primers. A 1375-bp kanamycin-resistance cassette (KanR) fragment was amplified from the pEGFP-C1 plasmid using the KanR-F and KanR-R primers. The KanR fragment and the two homologous arms were then inserted into the pUC18 vector to generate the recombinant pUC18-<italic>brbI</italic>-KanR plasmid which was transfected into electrocompetent <italic>B. suis</italic> S2 cells by electroporation (<xref ref-type="bibr" rid="B35">McQuiston et al., 1995</xref>). Transformants were screened in the presence of 50 &#x03BC;g/ml kanamycin. The mutant strain was further confirmed via polymerase chain reaction (PCR) using primers BrBI-TF and BrBI-TR and via quantitative real time polymerase chain reaction (qRT-PCR) using primers BrBI-qF and BrBI-qR.</p>
<p>To obtain the complemented strain, an expression plasmid containing a FLAG-tag (pBBARpc) was constructed (<xref ref-type="supplementary-material" rid="SM1">Supplementary Material</xref>). The intact <italic>brbI</italic> open reading frame (ORF) was amplified using the BrBI-F/BrBI-R primer pair, and the PCR product was cloned into pBBARpc to generate the recombinant pBBARpc-<italic>brbI</italic> plasmid. Electrocompetent cells of the mutant strain were prepared and electroporated according to standard procedures. Subsequently, transformants were screened in the presence of 100 &#x03BC;g/ml ampicillin. The complemented strain was further confirmed via western blot (WB) using an anti-FLAG antibody (abs830005, Absin, Shanghai, China) and via qRT-PCR using the primers, BrBI-qF and BrBI-qR. The sequences of all the primers used in this study are listed in <xref ref-type="table" rid="T1">Table 1</xref>.</p>
<table-wrap position="float" id="T1">
<label>TABLE 1</label>
<caption><p>Primers used in this study.</p></caption>
<table cellspacing="5" cellpadding="5" frame="hsides" rules="groups">
<thead>
<tr>
<td valign="top" align="left"><bold>Primer</bold></td>
<td valign="top" align="left"><bold>Sequence (5&#x2032;-3&#x2032;)</bold></td>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">UP-F</td>
<td valign="top" align="left"><underline>GCATGC</underline>CAGGTGACGGAGAACATT (<italic>Sph</italic> I, underlined)</td>
</tr>
<tr>
<td valign="top" align="left">UP-R</td>
<td valign="top" align="left"><underline>ACTAGT</underline>TGTGATTCATTATGGGGATTTC (<italic>Spe</italic> I, underlined)</td>
</tr>
<tr>
<td valign="top" align="left">DW-F</td>
<td valign="top" align="left"><underline>GTCGAC</underline>TGGGCAACCGTGAATAATC (<italic>Sal</italic> I, underlined)</td>
</tr>
<tr>
<td valign="top" align="left">DW-R</td>
<td valign="top" align="left"><underline>GAATTC</underline>AAGCGCAGCGAATGAAGAT (<italic>Eco</italic>R I, underlined)</td>
</tr>
<tr>
<td valign="top" align="left">KanR-F</td>
<td valign="top" align="left"><underline>GTCGAC</underline>TCAGGTGGCACTTTTCGGGGA (<italic>Sal</italic> I, underlined)</td>
</tr>
<tr>
<td valign="top" align="left">KanR-R</td>
<td valign="top" align="left"><underline>ACTAGT</underline>TTGGGCGTCGCTTGGTCGGT (<italic>Spe</italic> I, underlined)</td>
</tr>
<tr>
<td valign="top" align="left">BrBI-F</td>
<td valign="top" align="left"><underline>GTCGAC</underline>ATGGCTGACTTTCGTAAT (<italic>Sal</italic> I, underlined)</td>
</tr>
<tr>
<td valign="top" align="left">BrBI-R</td>
<td valign="top" align="left"><underline>GAATTC</underline>TTTCCGCTGCCCTTATCT (<italic>Eco</italic>R I, underlined)</td>
</tr>
<tr>
<td valign="top" align="left">BrBI-TF</td>
<td valign="top" align="left">CTTTCCAGCAGCTTCTTTTT</td>
</tr>
<tr>
<td valign="top" align="left">BrBI-TR</td>
<td valign="top" align="left">CAGCGTCGCCTTCATAGTA</td>
</tr>
<tr>
<td valign="top" align="left">BrBI-qF</td>
<td valign="top" align="left">CCTTATCCTTGCTTCGGTGG</td>
</tr>
<tr>
<td valign="top" align="left">BrBI-qR</td>
<td valign="top" align="left">TCCTTGATTTCCTGCGTATCG</td>
</tr>
<tr>
<td valign="top" align="left">orf2-qF</td>
<td valign="top" align="left">AACCATCACCGCCATCTATAC</td>
</tr>
<tr>
<td valign="top" align="left">orf2-qR</td>
<td valign="top" align="left">TTTGCCATTTCGTCCATTGTC</td>
</tr>
</tbody>
</table></table-wrap>
</sec>
<sec id="S2.SS4">
<title>Analysis of the Growth Phenotype of <italic>B. suis</italic> Strains</title>
<p>The three <italic>B. suis</italic> S2 strains were inoculated at equal densities (1 &#x00D7; 10<sup>7</sup> colony-forming units [CFUs]) into 10 ml of TSB and incubated at 37&#x00B0;C with shaking. The cultures were collected at specific times (0, 8, 16, 24, 32, 40, 48, 56, 64, 72, and 80 h), and the optical density (OD) at 600 nm was measured in a 96-well plate using a microplate reader (Bio-Rad, CA, United States). The three <italic>B. suis</italic> S2 strains were streaked on antibiotic-free TSA plates and incubated at 37&#x00B0;C. After 3 days, the colony sizes of the three <italic>B. suis</italic> S2 strains were measured.</p>
</sec>
<sec id="S2.SS5">
<title>Bacterial Aggregation Assay</title>
<p>The bacterial aggregation assay was performed as previously reported, with some modifications (<xref ref-type="bibr" rid="B31">Liu et al., 2017</xref>). Briefly, the three <italic>B. suis</italic> S2 strains were cultured to the exponential phase as described above, then the cultures were left standing on test-tube racks at room temperature for 24 h. To quantify the bacterial aggregation, the OD<sub>600 nm</sub> of 100 &#x03BC;l of the culture obtained before and after standing was measured in a 96-well microplate using a microplate reader (Bio-Rad). Percent bacterial aggregation was calculated using the following formula:</p>
<disp-formula id="S2.E1">
<label>(1)</label>
<mml:math id="M1">
<mml:mrow>
<mml:mo>%</mml:mo>
<mml:mi>b</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>o</mml:mi>
<mml:mpadded width="+3.3pt">
<mml:mi>n</mml:mi>
</mml:mpadded>
<mml:mo rspace="5.8pt">=</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mi>O</mml:mi>
<mml:mi>D</mml:mi>
<mml:mmultiscripts>
<mml:mo>-</mml:mo>
<mml:mprescripts/>
<mml:mrow>
<mml:mtext>total</mml:mtext>
</mml:mrow>
<mml:none/>
</mml:mmultiscripts>
<mml:mi>O</mml:mi>
<mml:mi>D</mml:mi>
<mml:mmultiscripts>
<mml:mo stretchy="false">]</mml:mo>
<mml:mprescripts/>
<mml:mrow>
<mml:mi>u</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:mi>p</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:mi>p</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:mi>e</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:mi>r</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:mi>p</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:mi>h</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:mi>a</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:mi>s</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:none/>
</mml:mmultiscripts>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:mi>O</mml:mi>
<mml:mi>D</mml:mi>
<mml:mmultiscripts>
<mml:mo rspace="5.8pt" stretchy="false">)</mml:mo>
<mml:mprescripts/>
<mml:mrow>
<mml:mtext>total</mml:mtext>
</mml:mrow>
<mml:none/>
</mml:mmultiscripts>
</mml:mrow>
<mml:mo rspace="5.8pt">&#x00D7;</mml:mo>
<mml:mn>100</mml:mn>
</mml:mrow>
</mml:math>
</disp-formula>
</sec>
<sec id="S2.SS6">
<title>Scanning Electron Microscopy (SEM) Assay</title>
<p>The three <italic>B. suis</italic> S2 strains were cultured to the exponential phase as described above, then inocula of the three strains were fixed in 2.5% glutaraldehyde solution for 3 h at room temperature. After washing three times with PBS, the bacterial samples were dehydrated in an ascending ethanol gradient (30, 50, 70, 80, 90, 95, and 100% for 15 min each). After critical point drying and gold coating, bacterial cells were observed and imaged using a scanning electron microscope (Regulus8100, Hitachi, Tokyo, Japan).</p>
</sec>
<sec id="S2.SS7">
<title>Analysis of <italic>B. suis</italic> S2 Strains Viability</title>
<p>The three <italic>B. suis</italic> S2 strains were prepared as described above and adjusted to the same OD<sub>600 nm</sub> with PBS. The concentrations of viable bacteria from the three <italic>B. suis</italic> S2 strains that had the same OD<sub>600 nm</sub> were measured by plating onto TSA plates as described above. Propidium iodide (PI) staining was also performed. 2 &#x00D7; 10<sup>7</sup> CFUs of each bacterial suspension was diluted into 100 &#x03BC;l with PBS and pipetted into separate wells of a 96-well, flat-bottom microplate, then incubated with isovolumetric PI stain solution (20 &#x03BC;M) at room temperature in the dark for 15 min. The fluorescence intensity was then measured (excitation wavelength: 535 nm; emission wavelength: 615 nm) using a fluorescence microplate reader (Spark 10M, TECAN, M&#x00E4;nnedorf, Switzerland).</p>
</sec>
<sec id="S2.SS8">
<title>Stress Assay</title>
<p>The three <italic>B. suis</italic> S2 strains were cultured and counted as described above. For the sodium dodecyl sulfate (SDS) sensitivity assay, 10 &#x03BC;l of 10-fold diluted samples from the three <italic>B. suis</italic> S2 strains (10<sup>4</sup>&#x2013;10<sup>6</sup> CFUs/ml) were added to a TSA plate containing 0.00325% SDS, then the plate was incubated for 3 days at 37&#x00B0;C.</p>
<p>For the antibiotic-resistance assay, bacterial samples were seeded into a 96-well microplate (2 &#x00D7; 10<sup>4</sup> CFUs/well) and treated with serially diluted antibiotics (polymyxin B: final concentrations of 1200, 600, 300, 150, 75, 37.5, 18.75, 9.325, and 0 &#x03BC;g/ml; lincomycin: final concentrations of 500, 375, 250, 125, 62.5, 31.25, 15.6, 7.8, 3.9, 1.95, and 0 &#x03BC;g/ml) at 37&#x00B0;C for 1h. The bacterial samples were then diluted 10-fold and plated onto TSA plates, then the plates were incubated for 3 days at 37 &#x00B0;C. The survival rate of each strain was calculated by dividing the antibiotic-treated bacterial CFUs by the untreated bacterial CFUs.</p>
<p>For the environmental stress tolerance assay, samples from the three <italic>B. suis</italic> S2 strains were inoculated into TSB containing H<sub>2</sub>O<sub>2</sub> (1.25 and 2.5 mM) for 1h and HCl (pH 3.5 and 4.5) for 30 min at 37&#x00B0;C. After treatment, a sample of each strain was diluted serially 10-fold and plated onto TSA plates. The plates were then incubated at 37&#x00B0;C for 3 days and the survival rate of each strain was calculated by dividing the treated bacterial CFUs by the untreated bacterial CFUs.</p>
</sec>
<sec id="S2.SS9">
<title>RNA Isolation and Quantitative Real-Time PCR</title>
<p>The <italic>B. suis</italic> S2 strains were grown in TSB to the exponential phase, then collected via centrifugation at 4&#x00B0;C. Total RNA was then isolated using TRIzol reagent (Invitrogen, Inc., Carlsbad, CA, United States). The resulting RNA samples were reverse-transcribed into cDNA using HiScript III RT SuperMix (Vazyme, Nanjing, China) per the manufacturer&#x2019;s recommended protocols. qRT-PCR was then performed using the ChamQ SYBR qPCR Master Mix (Vazyme, Nanjing, China) and Bio-Rad CFX96 Real-Time PCR Systems (Bio-Rad, United States). The 2<sup>&#x2013;&#x0394;&#x0394;<italic>Ct</italic></sup> method was used to analyze the relative transcription levels. The results for each target mRNA were normalized to 16S rRNA transcript levels and averaged.</p>
</sec>
<sec id="S2.SS10">
<title>Proteomic Profiling</title>
<sec id="S2.SS10.SSS1">
<title>Protein Preparation</title>
<p>The bacterial cultures (&#x223C;2 &#x00D7; 10<sup>7</sup> CFUs) were centrifuged at 12,000 &#x00D7; <italic>g</italic> for 5 min at 4&#x00B0;C. The bacterial pellets were then resuspended and homogenized in 300 &#x03BC;l of SDT buffer (4% SDS, 100 mM Dithiothreitol [DTT], 150 mM Tris-HCl pH 8.0). After centrifuging at 14,000 &#x00D7; <italic>g</italic> for 40 min, the supernatant was filtered with 0.22-&#x03BC;m filters and quantified with the BCA Protein Assay Kit (P0012, Beyotime).</p>
</sec>
<sec id="S2.SS10.SSS2">
<title>Protein Digestion and Tandem Mass Tag (TMT) Labeling</title>
<p>Proteins (200 &#x03BC;g per sample) were incorporated into 30 &#x03BC;l SDT buffer. The detergent, DTT and other low-molecular-weight components were removed using UA buffer (8 M urea, 150 mM Tris-HCl pH 8.5) by repeated ultrafiltration (Sartorius, 30 kDa). Next, 100 &#x03BC;l iodoacetamide (100 mM in UA buffer) was added to block reduced cysteine residues. The samples were incubated for 30 min in the dark. The filters were washed three times with 100 &#x03BC;l UA buffer, then twice with 100 &#x03BC;l 0.1 M triethylammonium bicarbonate (TEAB) buffer. Finally, the protein suspensions were digested with 4 &#x03BC;g trypsin (Promega) in 40 &#x03BC;l 0.1 M TEAB buffer overnight at 37&#x00B0;C, and the resulting peptides were collected as a filtrate. The peptide content was estimated via ultraviolet light spectral density at 280 nm using an extinctions coefficient of 1.1 of 0.1% (g/l) solution, which was calculated according to the frequencies of tryptophan and tyrosine in vertebrate proteins. Finally, 100 &#x03BC;g of peptide mixture per sample was labeled using TMT reagent according to the manufacturer&#x2019;s instructions (Thermo Fisher Scientific).</p>
</sec>
<sec id="S2.SS10.SSS3">
<title>Peptide Fractionation With Reversed Phase Chromatography</title>
<p>TMT-labeled peptides were fractionated by reversed phase chromatography using the Agilent 1260 Infinity II HPLC. The peptide mixture was diluted with buffer A (10 mM HCOONH4, 5% ACN, pH 10.0) and loaded onto a XBridge Peptide BEH C18 Column, 130 &#x00C5;, 5 &#x03BC;m, 4.6 mm &#x00D7; 100 mm column. The peptides were eluted at a flow rate of 1 ml/min with a gradient of 0%&#x2013;7% buffer B (10 mM HCOONH4, 85% ACN, pH 10.0) for 5 min, 7&#x2013;40% buffer B during 5&#x2013;40 min, 40&#x2013;100% buffer B during 45&#x2013;50 min, 100% buffer B during 50&#x2013;65 min. The elution was monitored at 214 nm based on the ultraviolet light trace, and fractions were collected every 1 min during 5&#x2013;50 min. The collected fractions were combined into 10 fractions and dried down via vacuum centrifugation at 45&#x00B0;C.</p>
</sec>
<sec id="S2.SS10.SSS4">
<title>Nano Liquid Chromatography-Tandem Mass Spectrometry (LC-MS/MS) Analysis</title>
<p>The peptide mixture was loaded onto the C18-reversed phase analytical column (Thermo Fisher Scientific, Acclaim PepMap RSLC 50 &#x03BC;m &#x00D7; 15 cm, nano viper, P/N164943) in buffer A (0.1% formic acid) and separated with a linear gradient of buffer B (80% acetonitrile and 0.1% formic acid) at flow 300 nl/min. The linear gradient was: 6% buffer B for 5 min, 6&#x2013;28% buffer B for 63 min, 28&#x2013;38% buffer B for 10 min, 38&#x2013;100% buffer B for 7 min, hold in 100% buffer B for 5 min.</p>
</sec>
<sec id="S2.SS10.SSS5">
<title>LC-MS/MS Analysis</title>
<p>LC-MS/MS analysis was performed on a Q Exactive plus mass spectrometer (Thermo Fisher Scientific) coupled to Easy nLC (Thermo Fisher Scientific) for 90 min. The mass spectrometer was operated in positive ion mode. MS data were acquired using a data-dependent top 10 method that dynamically chose the most abundant precursor ions from the survey scan (350&#x2013;1800 m/z) for HCD fragmentation. Survey scans were acquired at a resolution of 70000 at m/z 200 with an AGC target of 3e6 and a maxIT of 50 ms. MS2 scans were acquired at a resolution of 17500 for HCD spectra at m/z 200 with an AGC target of 2e5 and a maxIT of 45 ms. The isolation width was 2 m/z. Only ions with a charge state between 2 and 6 and a minimum intensity of 2e3 were selected for fragmentation. Dynamic exclusion for selected ions was 30 s. Normalized collision energy was 30 eV.</p>
</sec>
<sec id="S2.SS10.SSS6">
<title>Data Analysis</title>
<p>MS/MS raw files were processed using MASCOT engine (Matrix Science, London, United Kingdom; version 2.6) embedded into Proteome Discoverer 2.2, and searched against the NCBI database (Accession: <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="NZ_CP006961.1">NZ_CP006961.1</ext-link> and <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="NZ_CP006962.1">NZ_CP006962.1</ext-link>). The search parameters included trypsin as the enzyme used to generate peptides with a maximum of two missed cleavages permitted. A precursor mass tolerance of 10 ppm was specified and 0.05 Da tolerance for MS2 fragments. Except for TMT labels, carbamidomethyl (C) was set as a fixed modification. Variable modifications were Oxidation (M) and Acetyl (Protein N-term). A peptide and protein false discovery rate of 1% was enforced using a reverse database search strategy. Proteins with fold change (FC) &#x003E;1.5, <italic>p</italic>-value &#x003C; 0.05 (by Student&#x2019;s <italic>t</italic>-test), and the false discovery rate (FDR) &#x003C; 0.05 were considered differentially expressed.</p>
</sec>
</sec>
<sec id="S2.SS11">
<title>Transcriptomic Profiling</title>
<sec id="S2.SS11.SSS1">
<title>RNA Extraction</title>
<p>Total RNA was extracted from the tissue using TRIzol<sup>&#x00AE;</sup> Reagent according the manufacturer&#x2019;s instructions (Invitrogen) and genomic DNA was removed using DNase I (Takara). RNA quality was determined using the 2100 Bioanalyser (Agilent) and quantified using the ND-2000 (NanoDrop Technologies). A high-quality RNA sample (OD260 nm/280 nm = 1.8&#x223C;2.2, OD260 nm/230 nm &#x2265; 2.0, RIN &#x2265; 6.5, 28S:18S &#x2265; 1.0) was used to construct the sequencing library.</p>
</sec>
<sec id="S2.SS11.SSS2">
<title>Library Preparation, and Illumina HiSeq Sequencing</title>
<p>RNA-seq strand-specific libraries were prepared using the TruSeq RNA sample preparation kit from Illumina (San Diego, CA, United States), using 5 &#x03BC;g of total RNA. The rRNA by was removed using a RiboZero rRNA removal kit (Epicenter), and fragmented using fragmentation buffer. cDNA synthesis, end repair, A-base addition and ligation of the Illumina-index ed adaptors were performed according to Illumina&#x2019;s protocol. Libraries were then size selected for cDNA target fragments of 200&#x2013;300 bp on 2% Low Range Ultra Agarose followed by PCR amplification using Phusion DNA polymerase (NEB) for 15 PCR cycles. After quantified by TBS380, paired-end libraries were sequenced by Illumina NovaSeq 6000 sequencing (150 bp &#x00D7; 2).</p>
</sec>
<sec id="S2.SS11.SSS3">
<title>Read Quality Control and Mapping</title>
<p>The raw paired end reads were trimmed and quality controlled using Trimmomatic with the parameters, SLIDINGWINDOW: 4:15 MINLEN:75 (version 0.36<sup><xref ref-type="fn" rid="footnote1">1</xref></sup>). The clean reads were then separately aligned to the reference genome with orientation mode using Rockhopper software<sup><xref ref-type="fn" rid="footnote2">2</xref></sup> to perform a computational analysis of the bacterial RNA-seq data. Once input, Rockhopper uses RNA sequencing reads generated by high-throughput sequencing technology, which enables calculating gene expression levels with default parameters.</p>
</sec>
<sec id="S2.SS11.SSS4">
<title>Differential Expression Analysis and Functional Enrichment</title>
<p>To identify differentially expressed genes (DEGs) between the two samples, the expression level for each transcript was calculated using the fragments per kilobase of read per million mapped reads (RPKM) method. edgeR<sup><xref ref-type="fn" rid="footnote3">3</xref></sup> was used for differential expression analysis. The DEGs between the two samples were selected using the following criteria: (i) the logarithmic fold-change was &#x003E;2 and (ii) the FDR was &#x003C;0.05. To understand the functions of the DEGs, Gene Ontology (GO) functional enrichment analysis was performed using Goatools<sup><xref ref-type="fn" rid="footnote4">4</xref></sup>. DEGs were significantly enriched in GO terms when their Bonferroni-corrected <italic>p</italic>-value was &#x003C; 0.05.</p>
</sec>
</sec>
<sec id="S2.SS12">
<title>Bioinformatic Analysis</title>
<p>The amino acid sequences were aligned using Clustal Omega Webservers<sup><xref ref-type="fn" rid="footnote5">5</xref></sup> (<xref ref-type="bibr" rid="B33">Madeira et al., 2019</xref>). The protein transmembrane helices were predicted using TMHMM Server v. 2.0<sup><xref ref-type="fn" rid="footnote6">6</xref></sup>. The BrBI three-dimensional (3D) structure was predicted using the I-TASSER Webserver<sup><xref ref-type="fn" rid="footnote7">7</xref></sup> (<xref ref-type="bibr" rid="B37">Roy et al., 2010</xref>; <xref ref-type="bibr" rid="B46">Yang and Zhang, 2015</xref>).</p>
</sec>
<sec id="S2.SS13">
<title>Statistical Analysis</title>
<p>All experiments were independently repeated at least three times, and the data are expressed as the mean &#x00B1; standard deviation (SD) of representative experiments performed in triplicate. All data were analyzed with SPSS (SPSS, Inc., Cary, NC, United States). One-way analysis of variance (ANOVA) with Tukey&#x2019;s post-test was used to analyze differences between groups. A <italic>p</italic>-value &#x003C; 0.05 was considered statistically significant.</p>
</sec>
</sec>
<sec id="S3">
<title>Results</title>
<sec id="S3.SS1">
<title>BrBI Is Predicted to Be a Transmembrane Protein</title>
<p>As human BI-1 is the first and best characterized eukaryotic BI-1, and <italic>E. coli</italic> YccA and <italic>B. subtilis</italic> YetJ are relatively well-characterized prokaryotic homologs of eukaryotic BI-1, we compared BrBI with these three BI-1 family proteins using bioinformatic methods. BrBI shares 21.09% identity and 47.64% similarity to human BI-1, 21.12% identity and 59.36% similarity to <italic>E. coli</italic> YccA, and 21.43% identity and 47.01% similarity to <italic>B. subtilis</italic> YetJ in amino acid sequences (<xref ref-type="fig" rid="F1">Figure 1A</xref>). Additionally, BrBI contains seven predictive transmembrane helices with a 100% probability, which is concordant with human BI-1, <italic>E. coli</italic> YccA, and <italic>B. subtilis</italic> YetJ (<xref ref-type="fig" rid="F1">Figure 1B</xref>). As human BI-1, <italic>E. coli</italic> YccA, and <italic>B. subtilis</italic> YetJ have been identified as transmembrane proteins, we deduced that BrBI is also a transmembrane protein, and BI-1 family proteins may be structurally conserved. We further predicted the 3D structural model of BrBI. The predicted 3D structure of BrBI contained seven transmembrane segments (<xref ref-type="fig" rid="F1">Figure 1C</xref>) and is homologous to <italic>B. subtilis</italic> YetJ according to the Protein Data Bank library, further revealing that BrBI is a transmembrane protein homologous to other BI-1 family proteins.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption><p>Bioinformatic analysis of BrBI. <bold>(A)</bold> Multiple amino acid sequence alignment of BI-1 family proteins. Orange indicates residue G; yellow indicates residue P; green indicates residues T, S, N, and Q; blue indicates residues W, L, V, I, M, A, F, and C; cyan indicates residues H and Y; magenta indicates residues E and D; red indicates residues K and R; &#x201C;&#x002A;&#x201D; indicates fully conserved residues; &#x201C;:&#x201D; indicates strongly conserved residues; &#x201C;.&#x201D;indicates weakly conserved residues. <bold>(B)</bold> Prediction of BI-1 family protein transmembrane helices. The protein transmembrane helices were predicted using the TMHMM Server. <bold>(C)</bold> Predictive 3D structure of BrBI by I-TSAAER. <bold>(a)</bold> Red ribbons indicate transmembrane helices; green lines indicate loop structures; N indicates the N-terminus of BrBI; C indicates the C-terminus of BrBI. <bold>(b)</bold> Colored ribbons with numbers indicate 7 transmembrane helices of BrBI. <bold>(c,d)</bold> View of BrBI from the cytoplasm <bold>(c)</bold> and periplasmic space <bold>(d)</bold>. Colored ribbons with numbers indicate 7 transmembrane helices of BrBI. For <bold>(b&#x2013;d)</bold>, the ribbon was colored using the rainbow color scheme implemented in PyMOL.</p></caption>
<graphic xlink:href="fmicb-12-632095-g001.tif"/>
</fig>
</sec>
<sec id="S3.SS2">
<title>Construction of the <italic>brbI</italic> Mutant Strain and Its Complemented Strain</title>
<p>To explore the functions of BrBI, we generated a <italic>B. suis</italic> S2 <italic>brbI</italic> deletion mutant strain (&#x0394;<italic>brbI</italic>) via homologous recombination, replacing the intact <italic>brbI</italic> ORF with the kanamycin-resistance gene used as a selection marker (<xref ref-type="fig" rid="F2">Figure 2A</xref>). To further confirm successful generation of the &#x0394;<italic>brbI</italic> strain, we designed a primer pair BrBI-TF/TR residing within the <italic>brbI</italic> ORF (<xref ref-type="fig" rid="F2">Figure 2A</xref>). We randomly selected four kanamycin-resistant transformants for PCR verification and identified a transformant lacking <italic>brbI</italic> (<xref ref-type="fig" rid="F2">Figure 2B</xref>). This <italic>brbI</italic> mutant strain was further confirmed by qRT-PCR (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure S2A</xref>). Additionally, qRT-PCR indicated that the deletion of <italic>brbI</italic> did not affect the downstream gene expression (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure S2B</xref>). The complemented strain (&#x0394;<italic>brbI::brbI</italic>) was generated via exogenous expression of the intact <italic>brbI</italic> ORF in the &#x0394;<italic>brbI</italic> strain using an expression vector harboring a FLAG-tag. &#x0394;<italic>brbI::brbI</italic> strain generation was confirmed via WB using a mouse anti-FLAG antibody (<xref ref-type="fig" rid="F2">Figure 2C</xref>) and qRT-PCR (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure S2A</xref>).</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption><p>Verification of the &#x0394;<italic>brbI</italic> strain and its complemented strain. <bold>(A)</bold> Schematic of the &#x0394;<italic>brbI</italic> strain generation. <bold>(B)</bold> &#x0394;<italic>brbI</italic> strain verification via PCR using the BrBI-TF/TR primer pair. Lane 1 is the <italic>B. suis</italic> S2 strain; lanes 2&#x2013;5 are &#x0394;<italic>brbI</italic> strain candidates. Lane 5 is the confirmed mutant strain. <bold>(C)</bold> The &#x0394;<italic>brbI::brbI</italic> strain was verified via WB using a mouse anti-FLAG monoclonal antibody.</p></caption>
<graphic xlink:href="fmicb-12-632095-g002.tif"/>
</fig>
<sec id="S3.SS2.SSS1">
<title>Deletion of <italic>brbI</italic> Suppressed <italic>B. suis</italic> S2 Growth</title>
<p>When cultured in TSB, the &#x0394;<italic>brbI</italic> strain displayed a decreased growth rate compared with that of the <italic>B. suis</italic> S2 and &#x0394;<italic>brbI::brbI</italic> strains (<xref ref-type="fig" rid="F3">Figure 3A</xref>). When streaked on an antibiotic-free TSA plate, growth of the &#x0394;<italic>brbI</italic> strain remained significantly impaired, and the &#x0394;<italic>brbI</italic> strain exhibited smaller colonies than did the <italic>B. suis</italic> S2 and &#x0394;<italic>brbI::brbI</italic> strains (<xref ref-type="fig" rid="F3">Figure 3B</xref>). <italic>brbI</italic> deletion clearly suppressed the <italic>B. suis</italic> S2 growth capacity, indicating that deleting <italic>brbI</italic> may affect <italic>B. suis</italic> S2 division and/or survival.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption><p>Physiological phenotypes of the &#x0394;<italic>brbI</italic> strain. <bold>(A)</bold> Growth characteristics of the &#x0394;<italic>brbI</italic> strain in TSB. &#x201C;&#x002A;&#x201D; Indicates comparison of the <italic>B. suis</italic> S2 and &#x0394;<italic>brbI</italic> strains; &#x201C;#&#x201D; indicates comparison of the &#x0394;<italic>brbI::brbI</italic> and &#x0394;<italic>brbI</italic> strains. <bold>(B)</bold> Growth characteristics of the &#x0394;<italic>brbI</italic> strain on an antibiotic-free TSA plate. <bold>(C)</bold> Viable bacterial concentrations in the three <italic>B. suis</italic> S2 strains with the same OD<sub>600 nm</sub>. The <italic>Brucella</italic> samples were adjusted to the same OD600 nm, then the concentrations of viable bacteria from each sample were measured by plating them onto TSA plates. <bold>(D)</bold> Fluorescence intensity of the three <italic>B. suis</italic> S2 strain samples containing the same number of bacteria. Three <italic>B. suis</italic> S2 strain samples containing 2 &#x00D7; 10<sup>7</sup> CFUs each were taken and stained with PI, then the fluorescence intensity of each sample was measured. <bold>(E)</bold> Bacterial aggregation phenotype of the &#x0394;<italic>brbI</italic> strain. The three <italic>B. suis</italic> S2 strains were harvested at the exponential phase and left standing at room temperature for 24 h. The OD<sub>600 nm</sub> of 100 &#x03BC;l of the culture obtained before and after standing was measured to quantify the bacterial aggregation. <bold>(F)</bold> <italic>brbI</italic> deletion increases the sensitivity of <italic>B. suis</italic> S2 to SDS. Ten &#x03BC;l of <italic>Brucella</italic> samples were added to a TSA plate containing 0.00325% or no SDS and incubated for 24 h at 37&#x00B0;C. The data are expressed as the mean &#x00B1; SD of representative experiments performed in triplicate and analyzed using one-way ANOVA with Tukey&#x2019;s post-test. <sup>#</sup><italic>p</italic> &#x003C; 0.05,<sup>&#x2217;&#x2217;</sup>/##<italic>p</italic> &#x003C; 0.01, <sup>&#x002A;&#x002A;&#x002A;</sup>/###<italic>p</italic>&#x003C; 0.001.</p></caption>
<graphic xlink:href="fmicb-12-632095-g003.tif"/>
</fig>
</sec>
</sec>
<sec id="S3.SS3">
<title><italic>brbI</italic> Deletion Impaired <italic>B. suis</italic> S2 Viability</title>
<p>Because <italic>brbI</italic> deletion suppressed <italic>B. suis</italic> S2 growth, and because of the conserved cytoprotective capacity of eukaryotic BI-1 proteins and the bacterial cytoprotective potential of <italic>B. subtilis</italic> YetJ and <italic>E. coli</italic> YccA, we hypothesized that <italic>brbI</italic> deletion might lead to increased cell death in <italic>B. suis</italic> S2. To test this, we investigated how deleting <italic>brbI</italic> would affect the viability of <italic>B. suis</italic> S2. When adjusted to the same OD<sub>600 nm</sub> (representing the same bacterial cells), the &#x0394;<italic>brbI</italic> strain showed significantly reduced viable bacterial concentrations (represented as CFUs/ml) compared with those of the <italic>B. suis</italic> S2 and &#x0394;<italic>brbI::brbI</italic> strains (reduced &#x223C;4.8- and 3.3-fold, respectively, data not shown), indicating that <italic>brbI</italic> deletion suppressed the survival of <italic>B. suis</italic> S2 (<xref ref-type="fig" rid="F3">Figure 3C</xref>). We performed PI staining to further confirm this. PI is a red fluorescent nucleic acid stain that cannot penetrate intact biological membranes; therefore, the fluorescence intensity of bacteria stained with PI can indicate the cell membrane integrity status. If PI penetrates the membrane, the cell is usually assumed to be dead. When bacterial samples containing 2 &#x00D7; 10<sup>7</sup> CFUs were stained with PI, the &#x0394;<italic>brbI</italic> strain displayed considerably stronger fluorescence intensity than did either the <italic>B. suis</italic> S2 strain or the complemented strain (increased &#x223C;3.7- and 2.5-fold, respectively, data not shown), indicating a greater proportion of dead bacteria in the &#x0394;<italic>brbI</italic> strain samples (<xref ref-type="fig" rid="F3">Figure 3D</xref>). These results revealed that BrBI is crucial for <italic>B. suis</italic> S2 viability, suggesting that BrBI is a bacterially cytoprotective protein.</p>
</sec>
<sec id="S3.SS4">
<title><italic>brbI</italic> Deletion Affected the <italic>B. suis</italic> S2 Membrane Properties</title>
<p>In addition to growth suppression and viability defects, another noticeable phenotype was observed in the &#x0394;<italic>brbI</italic> strain. After standing at room temperature for 24 h, the <italic>B. suis</italic> S2 strain culture harvested at the exponential phase remained uniform, whereas conspicuous bacterial aggregation was observed in the &#x0394;<italic>brbI</italic> strain (<xref ref-type="fig" rid="F3">Figure 3E</xref>). Moreover, <italic>brbI</italic> complementation rescued the bacterial aggregation phenotype. We speculated that <italic>brbI</italic> deletion may have affected the membrane properties of <italic>B. suis</italic> S2. We then performed an SDS-sensitivity assay to determine the membrane properties of the &#x0394;<italic>brbI</italic> strain. As expected, <italic>brbI</italic> deletion significantly increased the sensitivity of <italic>B. suis</italic> S2 to SDS (<xref ref-type="fig" rid="F3">Figure 3F</xref>), suggesting that BrBI may be important for the membrane properties in <italic>B. suis</italic> S2.</p>
</sec>
<sec id="S3.SS5">
<title>SEM Analysis of the <italic>brbI</italic> Mutant Strain</title>
<p>Scanning electron microscopy analysis was conducted to assess the effect of <italic>brbI</italic> deletion on the <italic>B. suis</italic> S2 morphology (<xref ref-type="fig" rid="F4">Figure 4</xref>). The <italic>B. suis</italic> S2 strain displayed a typical coccobacillus morphology. Conversely, the cell morphology of the &#x0394;<italic>brbI</italic> strain was aberrant: some of the mutant strain showed Y-shaped or branched morphology, while some were intumesced mid-cell, indicating that <italic>brbI</italic> deletion impaired <italic>B. suis</italic> S2 cell division. Complementation of the <italic>brbI</italic> relieved the morphological irregularities. Together, BrBI may involve in the process of cell division in <italic>B. suis</italic> S2, which may contribute to the defective growth and viability in the &#x0394;<italic>brbI</italic> strain.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption><p>SEM analysis of the &#x0394;<italic>brbI</italic> strain. The <italic>B. suis</italic> S2 strain displayed a typical coccobacillus morphology (short-rod). The morphology of the &#x0394;<italic>brbI</italic> strain was aberrant, and manifested as Y-shaped or branched (arrows), intumesced or elongated bacterial cells (triangles), and fused bacterial cells (star). <italic>brbI</italic> complementation reversed the morphological irregularities. Bar = 4.00 &#x03BC;m.</p></caption>
<graphic xlink:href="fmicb-12-632095-g004.tif"/>
</fig>
</sec>
<sec id="S3.SS6">
<title><italic>brbI</italic> Deletion Reduced the Stress Resistance of <italic>B. suis</italic> S2</title>
<p>In bacteria, the cell envelope is the first layer to sense environmental stresses. Given the altered membrane properties of the &#x0394;<italic>brbI</italic> strain, we assessed whether <italic>brbI</italic> deletion would affect the sensitivity of <italic>B. suis</italic> S2 to environmental stresses. When exposed to an acidic environment, the survival rate was significantly lower in the &#x0394;<italic>brbI</italic> strain than in the <italic>B. suis</italic> S2 and &#x0394;<italic>brbI::brbI</italic> strains (<xref ref-type="fig" rid="F5">Figure 5A</xref>), indicating that BrBI is crucial for <italic>B. suis</italic> S2 to tolerate acidic stress. When stimulated with H<sub>2</sub>O<sub>2</sub>, cell viability was significantly impaired in the &#x0394;<italic>brbI</italic> strain (<xref ref-type="fig" rid="F5">Figure 5B</xref>), implying that BrBI also promotes <italic>B. suis</italic> S2 tolerance to oxidative stress.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption><p>Role of BrBI on stress resistance in <italic>B. suis</italic> S2. Bacterial samples were inoculated into TSB containing HCl <bold>(A)</bold> for 30 min and H<sub>2</sub>O<sub>2</sub> <bold>(B)</bold> for 1 h at 37&#x00B0;C. For antibiotic resistance, bacterial samples were seeded into a 96-well microplate and treated with serially diluted polymyxin B <bold>(C)</bold> and lincomycin <bold>(D)</bold> at 37&#x00B0;C for 1 h. After treatment, samples were diluted 10-fold and plated onto TSA plates to determine the survival rate of each strain by dividing the treated bacterial CFUs by the untreated bacterial CFUs. In panels C and D, &#x201C;&#x002A;&#x201D; indicates comparison of the <italic>B. suis</italic> S2 and &#x0394;<italic>brbI</italic> strains; &#x201C;#&#x201D; indicates comparison of the &#x0394;<italic>brbI::brbI</italic> and &#x0394;<italic>brbI</italic> strains. The data are expressed as the mean &#x00B1; SD of representative experiments performed in triplicate and analyzed using one-way ANOVA with Tukey&#x2019;s post-test. <sup>&#x002A;/#</sup><italic>p</italic> &#x003C; 0.05, <sup>&#x002A;&#x002A;/##</sup><italic>p</italic> &#x003C; 0.01, <sup>&#x002A;&#x002A;&#x002A;/###</sup><italic>p</italic> &#x003C; 0.001.</p></caption>
<graphic xlink:href="fmicb-12-632095-g005.tif"/>
</fig>
<p>Regarding antibiotic resistance, the bacterial cell envelope is both a barrier to antibiotic penetration and the target of several antibiotics. Hence, we evaluated the impact of <italic>brbI</italic> deletion on the resistance of <italic>B. suis</italic> S2 to polymyxin B and lincomycin. Polymyxin B is a peptide antibiotic that strongly interacts with bacterial membranes and is frequently used to assess the membrane characteristics of <italic>Brucella</italic>. Exposure to polymyxin B concentrations of up to 75 &#x03BC;g/ml significantly decreased the survival rate of the &#x0394;<italic>brbI</italic> strain compared with that of the complemented and <italic>B. suis</italic> S2 strains, indicating increased sensitivity of the &#x0394;<italic>brbI</italic> strain to polymyxin B (<xref ref-type="fig" rid="F5">Figure 5C</xref>). Lincomycin is a lincosamide antibiotic that inhibits bacterial protein synthesis by interacting with the ribosomal 50S subunits. When pre-treated with different lincomycin concentrations for 30 min at 37&#x00B0;C, the <italic>B. suis</italic> S2 and &#x0394;<italic>brbI::brbI</italic> strains displayed similar resistances to lincomycin; however, resistance of the &#x0394;<italic>brbI</italic> strain to lincomycin decreased significantly at 62.5 &#x03BC;g/ml of lincomycin (<xref ref-type="fig" rid="F5">Figure 5D</xref>). Further, <italic>brbI</italic> deletion significantly reduced the MICs of polymyxin B and lincomycin against <italic>B. suis</italic> S2 (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure S3</xref>). Thus, BrBI played a crucial role in modulating the sensitivity of <italic>B. suis</italic> S2 to antibiotics.</p>
</sec>
<sec id="S3.SS7">
<title>Proteomic Analysis</title>
<p>Given the prominent phenotypes induced by <italic>brbI</italic> deletion, we compared the proteomes of the three <italic>B. suis</italic> S2 strains. In total, 2270 proteins were identified in the three <italic>B. suis</italic> S2 strains, and 363 significantly affected proteins were differentially expressed in the &#x0394;<italic>brbI</italic> strain (<xref ref-type="fig" rid="F6">Figure 6A</xref>). Among these differentially expressed proteins (DEPs), 248 were down regulated, and 115 were up regulated. In the complemented strain, 138 proteins were differentially expressed, consisting of 94 down-regulated and 44 up-regulated proteins (<xref ref-type="fig" rid="F6">Figure 6B</xref>). In the &#x0394;<italic>brbI::brbI</italic> strain, complementation of <italic>brbI</italic> recovered 263 DEPs that were identified in the &#x0394;<italic>brbI</italic> strain but induced 38 &#x0394;<italic>brbI::brbI</italic>-specific DEPs, and 100 DEPs were overlapped with the &#x0394;<italic>brbI</italic> strain (<xref ref-type="fig" rid="F6">Figure 6C</xref>). This suggested that <italic>brbI</italic> may have extra functions when expressed beyond physiologic levels. <xref ref-type="fig" rid="F6">Figures 6C&#x2013;E</xref> show a slight discrepancy in the numbers of overlapped DEPs, attributed to two proteins, BrBI and a DUF1127 domain-containing protein (WP_002965878.1), which were up-regulated in the &#x0394;<italic>brbI::brbI</italic> strain but down-regulated in the &#x0394;<italic>brbI</italic> strain (<xref ref-type="supplementary-material" rid="SM1">Supplementary Table S1</xref>). Compared with the <italic>B. suis</italic> S2 strain, the remaining 98 overlapped DEPs were regulated in the same directions in the &#x0394;<italic>brbI</italic> and &#x0394;<italic>brbI::brbI</italic> strains. Although these 98 proteins were still differently expressed in the &#x0394;<italic>brbI::brbI</italic> strain compared with the <italic>B. suis</italic> S2 strain, 21 of them were significantly complemented (<xref ref-type="supplementary-material" rid="SM1">Supplementary Table S2</xref>), showing a 78.24% ([21 + 263]/363) complementation effect. Cluster analysis of 401 DEPs identified in the &#x0394;<italic>brbI</italic> and &#x0394;<italic>brbI::brbI</italic> strains showed a similar result (<xref ref-type="fig" rid="F6">Figure 6F</xref>); these DEPs were divided into &#x0394;<italic>brbI</italic>-specific, &#x0394;<italic>brbI::brbI</italic>-specific, significantly complemented, and uncomplemented DEPs, indicating an acceptable complementation of the complemented strain.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption><p>Comparison of the global proteomes of <italic>B. suis</italic> S2, &#x0394;<italic>brbI</italic>, and &#x0394;<italic>brbI::brbI</italic> by LC-MS/MS. <bold>(A)</bold> Volcano plot showing that 363 proteins were significantly differentially expressed in &#x0394;<italic>brbI</italic> compared with <italic>B. suis</italic> S2 with an absolute FC &#x003E; 1.5 and <italic>p</italic> &#x003C; 0.05. <bold>(B)</bold> Volcano plot showing that 115 proteins were significantly differentially expressed in &#x0394;<italic>brbI::brbI</italic> compared with <italic>B. suis</italic> S2 with an absolute FC &#x003E; 1.5 and <italic>p</italic> &#x003C; 0.05. <bold>(C)</bold> Venn diagram comparing the DEPs between &#x0394;<italic>brbI</italic> and &#x0394;<italic>brbI::brbI</italic>. The blue circle indicates DEPs between <italic>B. suis</italic> S2 and &#x0394;<italic>brbI</italic>; the pink circle indicates DEPs between <italic>B. suis</italic> S2 and &#x0394;<italic>brbI::brbI</italic>; the overlapped area indicates DEPs identified in both &#x0394;<italic>brbI</italic> and &#x0394;<italic>brbI::brbI</italic>. <bold>(D)</bold> Venn diagram comparing the up-regulated proteins between &#x0394;<italic>brbI</italic> and &#x0394;<italic>brbI::brbI</italic>. The blue circle indicates up-regulated proteins between <italic>B. suis</italic> S2 and &#x0394;<italic>brbI</italic>; the pink circle indicates up-regulated proteins between <italic>B. suis</italic> S2 and &#x0394;<italic>brbI::brbI</italic>; the overlapped area indicates up-regulated proteins in both &#x0394;<italic>brbI</italic> and &#x0394;<italic>brbI::brbI</italic>. <bold>(E)</bold> Venn diagram comparing the down-regulated proteins between &#x0394;<italic>brbI</italic> and &#x0394;<italic>brbI::brbI</italic>. The blue circle indicates down-regulated proteins between <italic>B. suis</italic> S2 and &#x0394;<italic>brbI</italic>; the pink circle indicates down-regulated proteins between <italic>B. suis</italic> S2 and &#x0394;<italic>brbI::brbI</italic>; the overlapped area indicates down-regulated proteins in both of &#x0394;<italic>brbI</italic> and &#x0394;<italic>brbI::brbI</italic>. <bold>(F)</bold> Heat map of the two-way hierarchical clustering (436 proteins from three strains satisfying FC &#x003E; 1.5 using <italic>Z</italic>-score for normalized values [log2 based]). <bold>(G)</bold> GO enrichment analysis of the DEPs between <italic>B. suis</italic> S2 and &#x0394;<italic>brbI</italic> according to the cellular component, molecular function, and biological process categories.</p></caption>
<graphic xlink:href="fmicb-12-632095-g006.tif"/>
</fig>
<p>Gene Ontology analysis was performed to further analyse the effects of <italic>brbI</italic> deletion in <italic>B. suis</italic> S2. The 363 DEPs between the &#x0394;<italic>brbI</italic> and <italic>B. suis</italic> S2 strains were enriched to GO terms at the second level in three GO categories: biological process, cellular component and molecular function (<xref ref-type="fig" rid="F6">Figure 6G</xref> and <xref ref-type="supplementary-material" rid="SM1">Supplementary Table S3</xref>). For biological process, <italic>brbI</italic> deletion mainly affected proteins associated with &#x201C;regulation of DNA-templated transcription,&#x201D; &#x201C;oxidation-reduction,&#x201D; &#x201C;transport,&#x201D; &#x201C;proteolysis,&#x201D; and &#x201C;protein folding.&#x201D; For cellular component, BrBI was primarily involved in the &#x201C;plasma membrane,&#x201D; &#x201C;outer membrane-bounded periplasmic space,&#x201D; &#x201C;cytoplasm,&#x201D; &#x201C;integral component of plasma membrane,&#x201D; and &#x201C;cytosol.&#x201D; For molecular function, BrBI mainly regulated proteins associated with &#x201C;protein binding,&#x201D; &#x201C;ATPase activity coupled to transmembrane movement of substances,&#x201D; &#x201C;transporter activity,&#x201D; &#x201C;catalytic activity,&#x201D; and &#x201C;DNA binding transcription factor activity.&#x201D; In summary, <italic>brbI</italic> deletion predominantly affected proteins associated with the membrane, transporter, and transcription.</p>
</sec>
<sec id="S3.SS8">
<title>Transcriptomic Analysis</title>
<p>Given the effect of <italic>brbI</italic> deletion on transcriptional activities and because 14 of 363 DEPs identified in the &#x0394;<italic>brbI</italic> strain were transcriptional factors (<xref ref-type="supplementary-material" rid="SM1">Supplementary Table S4</xref>), we compared the transcriptomes of the three <italic>B. suis</italic> S2 strains.</p>
<p>In total, we obtained 14,610,410, 11,764,696, and 12,634,380 clean reads from the <italic>B. suis</italic> S2 strain, 12,396,722, 9,776,270, and 9,215,932 clean reads from the &#x0394;<italic>brbI</italic> strain, and 10,609,342, 18,576,960, and 12,138,906 clean reads from the &#x0394;<italic>brbI::brbI</italic> strain. Further, 87.46&#x2013;93.46% clean reads were mapped to the <italic>B. suis</italic> S2 reference genome, and 2986 mRNAs were identified in the three <italic>B. suis</italic> S2 strains (<xref ref-type="supplementary-material" rid="SM1">Supplementary Table S5</xref>). Compared with the <italic>B. suis</italic> S2 strain, 677 genes were significantly differentially transcribed in the &#x0394;<italic>brbI</italic> strain, including 347 down-regulated and 330 up-regulated genes (<xref ref-type="fig" rid="F7">Figure 7A</xref> and <xref ref-type="supplementary-material" rid="SM1">Supplementary Table S6</xref>). Compared with the <italic>B. suis</italic> S2 strain, 144 up-regulated and 170 down-regulated genes were identified in the &#x0394;<italic>brbI::brbI</italic> strain (<xref ref-type="fig" rid="F7">Figure 7B</xref> and <xref ref-type="supplementary-material" rid="SM1">Supplementary Table S6</xref>). In the &#x0394;<italic>brbI::brbI</italic> strain, <italic>brbI</italic> complementation recovered 442 DEGs that were identified in the &#x0394;<italic>brbI</italic> strain and induced 79 &#x0394;<italic>brbI::brbI</italic>-specific DEGs, and 235 DEGs were overlapped between the &#x0394;<italic>brbI</italic> and &#x0394;<italic>brbI::brbI</italic> strains (<xref ref-type="fig" rid="F7">Figures 7C&#x2013;E</xref>). This also suggested that <italic>brbI</italic> may have extra functions when expressed beyond physiologic status. Compared with the <italic>B. suis</italic> S2 strain, the 235 overlapped DEGs were regulated in the same directions in the &#x0394;<italic>brbI</italic> and &#x0394;<italic>brbI::brbI</italic> strains, and 68 of them were significantly complemented in the &#x0394;<italic>brbI::brbI</italic> strain (<xref ref-type="supplementary-material" rid="SM1">Supplementary Table S7</xref>), displaying a 75.33% ([68+442]/677) complementation effect. Cluster analysis revealed that 756 DEGs identified in the &#x0394;<italic>brbI</italic> and &#x0394;<italic>brbI::brbI</italic> strains could be divided into &#x0394;<italic>brbI</italic>-specific, &#x0394;<italic>brbI::brbI</italic>- specific, significantly complemented, and uncomplemented DEGs (<xref ref-type="fig" rid="F7">Figure 7F</xref>). Considering the stronger sensitivity of transcriptomics to changes, the 75.33% complementation should be an acceptable effect in the complemented strain.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption><p>Comparison of the global transcriptomes of <italic>B. suis</italic> S2, &#x0394;<italic>brbI</italic>, and &#x0394;<italic>brbI::brbI</italic> by RNA-sequencing. <bold>(A)</bold> Volcano plot showing that 677 genes were significantly differentially expressed in &#x0394;<italic>brbI</italic> compared with <italic>B. suis</italic> S2 with an FC &#x003E; 2.0 and <italic>p</italic> &#x003C; 0.05. <bold>(B)</bold> Volcano plot showing that 314 genes were significantly differentially expressed in &#x0394;<italic>brbI::brbI</italic> compared with <italic>B. suis</italic> S2 with an FC &#x003E; 2.0 and <italic>p</italic> &#x003C; 0.05. <bold>(C)</bold> Venn diagram comparing the DEGs between &#x0394;<italic>brbI</italic> and &#x0394;<italic>brbI::brbI</italic>. The blue circle indicates DEGs between <italic>B. suis</italic> S2 and &#x0394;<italic>brbI</italic>; the pink circle indicates DEGs between <italic>B. suis</italic> S2 and &#x0394;<italic>brbI::brbI</italic>; the overlapped area indicates DEGs in both &#x0394;<italic>brbI</italic> and &#x0394;<italic>brbI::brbI</italic>. <bold>(D)</bold> Venn diagram comparing the up-regulated genes between &#x0394;<italic>brbI</italic> and &#x0394;<italic>brbI::brbI</italic>. The blue circle indicates up-regulated genes between <italic>B. suis</italic> S2 and &#x0394;<italic>brbI</italic>; the pink circle indicates up-regulated genes between <italic>B. suis</italic> S2 and &#x0394;<italic>brbI::brbI</italic>; the overlapped area indicates up-regulated genes in both of &#x0394;<italic>brbI</italic> and &#x0394;<italic>brbI::brbI</italic>. <bold>(E)</bold> Venn diagram comparing the down-regulated genes between &#x0394;<italic>brbI</italic> and &#x0394;<italic>brbI::brbI</italic>. The blue circle indicates down-regulated genes between <italic>B. suis</italic> S2 and &#x0394;<italic>brbI</italic>; the pink circle indicates down-regulated genes between <italic>B. suis</italic> S2 and &#x0394;<italic>brbI::brbI</italic>; the overlapped area indicates down-regulated genes in both of &#x0394;<italic>brbI</italic> and &#x0394;<italic>brbI::brbI</italic>. <bold>(F)</bold> Heat map of the two-way hierarchical clustering (754 genes from three strains satisfying an FC &#x003E; 2.0 using <italic>Z</italic>-score for normalized values [log2 based]). <bold>(G)</bold> GO enrichment analysis of the DEGs between the <italic>B. suis</italic> S2 and &#x0394;<italic>brbI</italic> strains according to the cellular component, molecular function, and biological process categories.</p></caption>
<graphic xlink:href="fmicb-12-632095-g007.tif"/>
</fig>
<p>Gene Ontology enrichment analysis showed that the highly enriched GO terms in the &#x0394;<italic>brbI</italic> strain were mainly associated with the biological processes of &#x201C;metabolic process,&#x201D; &#x201C;cellular process,&#x201D; &#x201C;organic substance metabolic process,&#x201D; &#x201C;cellular metabolic process,&#x201D; and &#x201C;primary metabolic process&#x201D;; the cellular component of &#x201C;cell part,&#x201D; &#x201C;cell,&#x201D; &#x201C;plasma membrane,&#x201D; &#x201C;cytoplasm,&#x201D; and &#x201C;membrane&#x201D;; the molecular function of &#x201C;catalytic activity,&#x201D; &#x201C;binding,&#x201D; and &#x201C;protein binding&#x201D; (<xref ref-type="fig" rid="F7">Figure 7G</xref> and <xref ref-type="supplementary-material" rid="SM1">Supplementary Table S8</xref>). These enriched terms were mostly consistent with those in the proteomic profiles with a slight difference in that the predominant genes enriched in the transcriptomic profiles were involved in multiple metabolic processes. This suggested that BrBI played extensive roles in <italic>B. suis</italic> S2 at both the protein and mRNA levels.</p>
</sec>
<sec id="S3.SS9">
<title>Integrated Analysis Revealed the Global Roles of BrBI in <italic>B. suis</italic> S2</title>
<p>According to the proteomic and transcriptomic profiles, <italic>brbI</italic> deletion extensively affected the physiology of <italic>B. suis</italic> S2 at both the protein and mRNA levels. Therefore, we integrated transcriptomic and proteomic analyses to explore the in-depth effects of <italic>brbI</italic> deletion on <italic>B. suis</italic> S2.</p>
<p>As mentioned, 2986 mRNAs and 2270 proteins were identified in the transcriptomic and proteomic profiles, respectively, and these identified proteins could all be assigned to the identified mRNAs (<xref ref-type="fig" rid="F8">Figure 8A</xref>). Thus, we focused on the 2270 genes that were identified in both the transcriptomic and proteomic profiles. Of the DEGs identified in the &#x0394;<italic>brbI</italic> strain, 422 were included in the 2270 genes (<xref ref-type="fig" rid="F8">Figure 8B</xref>). Therefore, 635 genes were significantly affected at the mRNA and/or protein levels in the &#x0394;<italic>brbI</italic> strain (<xref ref-type="fig" rid="F8">Figure 8C</xref>), and we thus targeted these genes going forward. Among the 635 genes, 213 were significantly affected at only the protein level (mRNA_no and protein_sig), 272 were significantly affected at only the mRNA level (mRNA_sig and protein_no), and 150 were significantly affected at both the mRNA and protein levels (mRNA_sig and protein_sig) (<xref ref-type="fig" rid="F8">Figure 8C</xref> and <xref ref-type="supplementary-material" rid="SM1">Supplementary Table S4</xref>).</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption><p>Integrated proteomic and transcriptomic analyses. <bold>(A)</bold> Venn diagram showing that the 2270 proteins identified in the proteomes can all be assigned to the mRNAs identified in the transcriptomes. The blue circle indicates all mRNAs identified in &#x0394;<italic>brbI</italic>; the pink circle indicates all proteins identified in &#x0394;<italic>brbI</italic>. <bold>(B)</bold> Venn diagram showing that 422 of 677 DEGs in the &#x0394;<italic>brbI</italic> strain can be assigned to the 2270 genes identified in both the proteomes and transcriptomes. The blue circle indicates DEGs in &#x0394;<italic>brbI</italic>, the pink circle indicates all proteins in &#x0394;<italic>brbI</italic>; the overlapped area indicates DEGs with corresponding proteins in the proteomic profiles. <bold>(C)</bold> Venn diagram showing the effects of <italic>brbI</italic> deletion on <italic>B. suis</italic> S2 at the mRNA and/or protein levels. The blue ellipse indicates genes significantly affected by <italic>brbI</italic> deletion in the transcriptomic profiles; the pink ellipse indicates genes significantly affected by <italic>brbI</italic> deletion in the proteomic profiles; the green ellipse indicates genes not affected by <italic>brbI</italic> deletion in the transcriptomic profiles; the purple ellipse indicates genes not affected by <italic>brbI</italic> deletion in the proteomic profiles.</p></caption>
<graphic xlink:href="fmicb-12-632095-g008.tif"/>
</fig>
<p>The existence of the mRNA_no and protein_sig genes in the &#x0394;<italic>brbI</italic> strain indicated that BrBI participated in the processes of protein translation, stability, and degradation. Because BrBI is a transmembrane protein, and its homolog, <italic>E. coli</italic> YccA, has been characterized as a protease activity-associated protein, we speculated that BrBI may initiate its extensive roles by regulating protein degradation. Consequently, the occurrence of 272 genes being affected at only the mRNA level suggested that some of the DEPs induced by <italic>brbI</italic> deletion should be transcription factors, which was revealed in the proteomic analysis. However, these mRNA_sig and protein_no genes were unaffected by <italic>brbI</italic> deletion at the protein level, indicating that these genes may have a powerful and possibly BrBI-independent protein regulation pathway to overcome the altered gene transcriptions and restore the corresponding proteins to a physiological level. Most of the mRNA_sig and protein_sig genes were regulated in the same directions at the mRNA and protein levels, indicating that these genes should be regulated by <italic>brbI</italic> deletion at only the mRNA level. However, 9 out of 150 mRNA_sig and protein_sig genes showed opposite regulation trends at the mRNA and protein levels. <italic>BSS2_RS10645</italic>, <italic>BSS2_RS10185</italic>, and <italic>BSS2_RS10195</italic> were down-regulated at the mRNA level but up-regulated at protein level, and <italic>BSS2_RS04340</italic>, <italic>BSS2_RS01465</italic>, <italic>BSS2_RS00810</italic>, <italic>BSS2_RS01210</italic>, <italic>BSS2_RS15335</italic>, and <italic>BSS2_RS06495</italic> were upregulated at the mRNA level but downregulated at the protein level (<xref ref-type="supplementary-material" rid="SM1">Supplementary Table S4</xref>), indicating that they should have a powerful, predominant, and BrBI-dependent protein regulation pathway that can reverse the transcription alterations induced by <italic>brbI</italic> deletion.</p>
<p>According to the known <italic>Brucella</italic> genome annotations, five groups of genes were summarized from the 635 genes discussed here: membrane protein-, cell-division protein-, transcriptional regulator-, peptidase/protease-, and transporter-associated protein-encoding genes (<xref ref-type="supplementary-material" rid="SM1">Supplementary Table S4</xref>). In total, 27 membrane protein-, 4 cell-division protein-, 31 transcriptional regulator-, 21 peptidase/protease-, and 81 transporter-associated protein-encoding genes were identified. The peptidase/protease- and transporter-associated protein-encoding genes were generally distributed into the mRNA_sig and protein_sig (4 genes and 26 genes, respectively), mRNA_no and protein_sig (14 genes and 25 genes, respectively), and mRNA_sig and protein_no (3 genes and 30 genes, respectively) groups. The transcriptional factor-encoding genes were distributed into the mRNA_sig and protein_no (14 genes) and mRNA_no and protein_sig (17 genes) groups. The membrane protein-encoding genes were mainly distributed into the mRNA_sig and protein_sig (seven genes) and mRNA_no &#x0026; protein_sig (20 genes) groups. The cell-division protein-encoding genes were specifically distributed into the mRNA_no and protein_sig group.</p>
<p>These results revealed that BrBI plays extensive roles in <italic>B. suis</italic> S2 via protein and/or transcriptional pathways. Importantly, the generally distributed transcriptional factor- and peptidase/protease-encoding genes could provide insights into the underlying mechanisms of BrBI. BrBI may initially modulate transcriptional factors or peptidase/proteases at the protein level, which may further trigger regulations on other regulatory or functional genes at the mRNA and/or protein levels. Notably, the transporter-associated protein-encoding genes were the most abundant and generally distributed, which may explain why <italic>brbI</italic> deletion led to the significantly increased sensitivity of <italic>B. suis</italic> S2 to environmental stresses and antibiotics. Most of the transporter-associated proteins were membrane-related, indicating that BrBI plays a crucial role in the <italic>B. suis</italic> S2 membrane. Further, 27 membrane protein-encoding genes were identified, and interestingly, most were mRNA_no and protein_sig genes. The identified cell-division protein-encoding genes were all mRNA_no and protein_sig genes, indicating a specific role of BrBI in <italic>B. suis</italic> S2 division. Therefore, the role of BrBI in the <italic>B. suis</italic> S2 division and membrane are discussed in detail in the following sections.</p>
</sec>
<sec id="S3.SS10">
<title>BrBI Is Involved in <italic>B. suis</italic> S2 Division</title>
<p>In Gram-negative bacteria, cell division is tightly mediated by a multi-protein complex called the divisome that is assembled to the mid-cell. FtsZ, FtsA, FtsB, FtsI, FtsL, FtsQ, and FtsW are the primary components of the divisome that execute bacterial cytokinesis (<xref ref-type="bibr" rid="B14">Dajkovic and Lutkenhaus, 2006</xref>; <xref ref-type="bibr" rid="B42">Szwedziak and Ghosal, 2017</xref>). The divisome is assembled in two stages. First, FtsZ polymers are attached to the mid-cell membrane by interacting with FtsA to form a Z-ring scaffold, which provides the constrictive force. Second, the late-division proteins, including FtsE, FtsK, FtsQ, FtsL and FtsB, FtsW, FtsI, and FtsN, are sequentially recruited to the Z-ring to form the complete septum to carry out cell division (<xref ref-type="bibr" rid="B14">Dajkovic and Lutkenhaus, 2006</xref>; <xref ref-type="bibr" rid="B42">Szwedziak and Ghosal, 2017</xref>). The early (Z-ring) and late divisome components are linked by the FtsQBL complex, which may have a structural role as a scaffold in assembling the divisome (<xref ref-type="bibr" rid="B12">Choi et al., 2018</xref>; <xref ref-type="bibr" rid="B13">Condon et al., 2018</xref>). In most bacteria, to guarantee that division occurs at the correct site, a Min system comprising MinC, MinD, and MinE prevents establishing the divisome at bacterial poles (<xref ref-type="bibr" rid="B42">Szwedziak and Ghosal, 2017</xref>; <xref ref-type="bibr" rid="B36">Ramm et al., 2019</xref>). As a spatial modulator of the divisome, defects in the Min system produce small anucleate minicells from the poles of rod-shaped mother cells, increasing the average cell length of the DNA-containing mother cells (<xref ref-type="bibr" rid="B32">Lutkenhaus, 2007</xref>).</p>
<p>As mentioned, <italic>brbI</italic> deletion significantly affected four cell-division protein-encoding genes, indicating a crucial role of BrBI in <italic>B. suis</italic> S2 division. Therefore, we integrated proteomic and transcriptomic analyses to thoroughly explore how <italic>brbI</italic> deletion would affect <italic>B. suis</italic> S2 division. First, <italic>brbI</italic> deletion did not affect FtsZ and FtsA expression at either the mRNA or protein levels, indicating a minor role of BrBI in Z-ring assembly (<xref ref-type="fig" rid="F9">Figure 9</xref> and <xref ref-type="supplementary-material" rid="SM1">Supplementary Table S9</xref>). Further, deleting <italic>brbI</italic> did not affect the late divisome proteins, FtsE, FtsK, and FtsW, at either the mRNA or protein levels; however, <italic>brbI</italic> deletion significantly upregulated FtsB, FtsI, FtsL, and FtsQ at only the protein level, possibly indicating disordered cytokinesis septum formation (<xref ref-type="fig" rid="F9">Figure 9</xref> and <xref ref-type="supplementary-material" rid="SM1">Supplementary Table S9</xref>). Moreover, <italic>brbI</italic> deletion did not affect MinC, MinD, and MinE expression at either the mRNA or protein levels, indicating a minor role of BrBI in the Min system (<xref ref-type="fig" rid="F9">Figure 9</xref> and <xref ref-type="supplementary-material" rid="SM1">Supplementary Table S9</xref>).</p>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption><p>Cluster analysis for cell division-associated genes. <bold>(A)</bold> Cluster analysis for cell division-associated genes at the protein level. <bold>(B)</bold> Cluster analysis for cell division-associated genes at the mRNA level.</p></caption>
<graphic xlink:href="fmicb-12-632095-g009.tif"/>
</fig>
<p>In summary, <italic>brbI</italic> deletion may affect proper formation of the division septum, which may further lead to aberrant <italic>B. suis</italic> S2 division. Specifically, the role of BrBI in <italic>B. suis</italic> S2 division may occur predominantly by regulating the translation/degradation/stability of FtsB, FtsI, FtsL, and FtsQ.</p>
</sec>
<sec id="S3.SS11">
<title>BrBI Is Crucial for <italic>B. suis</italic> S2 Membrane Homeostasis</title>
<p>The integrated proteomic and transcriptomic analyses showed that <italic>brbI</italic> deletion significantly affected the 27 identified membrane protein-encoding genes (excluding <italic>ftsB</italic>, <italic>ftsI</italic>, <italic>ftsL</italic>, and <italic>ftsQ</italic>); of these, seven were affected at both the mRNA and protein levels, and 20 were affected only at the protein level (<xref ref-type="supplementary-material" rid="SM1">Supplementary Table S4</xref>). Outer membrane proteins (OMPs) BamA and LptD were significantly down-regulated at only the protein level in the &#x0394;<italic>brbI</italic> strain (<xref ref-type="fig" rid="F10">Figure 10</xref> and <xref ref-type="supplementary-material" rid="SM1">Supplementary Table S9</xref>). In Gram-negative bacteria, the lipopolysaccharide transport (Lpt), &#x03B2;-barrel assembly machine (Bam), and the localization of lipoproteins (Lol) pathways are required to transport the outer membrane components (<xref ref-type="bibr" rid="B41">Sperandeo et al., 2019</xref>). Lpt, Bam, and Lol are necessary for respectively transporting LPS, &#x03B2;-barrel proteins, and lipoproteins across the periplasm to the outer membrane. Thus, down-regulation of BamA and LptD induced by <italic>brbI</italic> deletion indicated an alteration of the <italic>B. suis</italic> S2 outer membrane. Consistently, the OMP2a, OMP2b, OMP22, OMP25b, OMP31, and OMPW family protein (WP_002964666.1) and outer membrane &#x03B2;-barrel protein (WP_004689965.1) were downregulated at only the protein level with the exception of OMP25b, which was downregulated at both the mRNA and protein levels in the &#x0394;<italic>brbI</italic> strain (<xref ref-type="fig" rid="F10">Figure 10</xref> and <xref ref-type="supplementary-material" rid="SM1">Supplementary Table S9</xref>). In <italic>Brucella</italic>, OMP2a and OMP2b are group 2 OMPs, and OMP22, OMP25b, and OMP31 are group 3 OMPs, which are &#x03B2;-barrel proteins (<xref ref-type="bibr" rid="B18">Goolab et al., 2015</xref>). Therefore, down-regulation of OMP25b induced by <italic>brbI</italic> deletion may occur at the mRNA level; however, regulation of the remaining OMPs by <italic>brbI</italic> deletion may occur indirectly through BamA at the protein level.</p>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption><p>Cluster analysis for membrane protein encoding genes. <bold>(A)</bold> Cluster analysis for membrane protein encoding genes at the protein level. <bold>(B)</bold> Cluster analysis for membrane protein encoding genes at the mRNA level.</p></caption>
<graphic xlink:href="fmicb-12-632095-g010.tif"/>
</fig>
<p>The FtsH protease modulators, HflC and HflK, were also significantly up-regulated at only the protein level via <italic>brbI</italic> deletion (<xref ref-type="fig" rid="F10">Figure 10</xref> and <xref ref-type="supplementary-material" rid="SM1">Supplementary Table S9</xref>). FtsH is a membrane-anchored ATP-dependent metalloprotease that degrades redundant or abnormal membrane proteins to maintain membrane homeostasis in prokaryotes (<xref ref-type="bibr" rid="B27">Langklotz et al., 2012</xref>; <xref ref-type="bibr" rid="B3">Bittner et al., 2017</xref>). HflC and HflK form a membrane protein complex (HflKC) and further complex with FtsH. It has been reported that HflKC may act as a negative modulator for the proteolytic activity of FtsH (<xref ref-type="bibr" rid="B27">Langklotz et al., 2012</xref>). Consequently, we speculate that the upregulated HflC and HflK by <italic>brbI</italic> deletion may have suppressed the proteolytic activity of FtsH, subsequently leading to an accumulation of membrane proteins or cytosolic regulatory proteins that were meant to be degraded, which eventually perturbed <italic>B. suis</italic> S2 membrane homeostasis.</p>
<p>Moreover, <italic>brbI</italic> deletion significantly affected several inner membrane proteins and dozens of membrane-related transporter proteins, but the mechanism by which BrBI regulates these proteins remains unknown. Nevertheless, the integrated proteomic and transcriptomic analyses revealed the extensive roles of BrBI in <italic>B. suis</italic> S2 membrane homeostasis. Specifically, OMP was regulated by BrBI predominantly at the protein level.</p>
</sec>
</sec>
<sec id="S4">
<title>Discussion</title>
<p>In this study, we explored the role of BrBI in <italic>B. suis</italic> S2 via basic physiological tests and integrated proteomic and transcriptomic analyses. BrBI played extensive roles in the physiological processes of <italic>B. suis</italic> S2 at both the mRNA and protein levels, including membrane homeostasis, cell division, transport activity, transcription activity, and proteolysis.</p>
<p>Deleting <italic>brbI</italic> in <italic>B. suis</italic> S2 yielded the prominent phenotype of suppressed bacterial growth and viability. Eukaryotic BI-1 has been characterized as a highly conserved cytoprotective protein, although the mechanisms of apoptosis differ between animals and plants (<xref ref-type="bibr" rid="B21">Huckelhoven, 2004</xref>; <xref ref-type="bibr" rid="B45">Watanabe and Lam, 2009</xref>; <xref ref-type="bibr" rid="B19">Henke et al., 2011</xref>; <xref ref-type="bibr" rid="B6">Bultynck et al., 2012</xref>). Bacterial (prokaryotic) and eukaryotic BI-1 share a conserved structure (<xref ref-type="bibr" rid="B25">Kihara et al., 1998</xref>; <xref ref-type="bibr" rid="B16">Gamboa-Tuz et al., 2018</xref>), which is consistent with our finding that BrBI shared highly concordant transmembrane helices with <italic>E. coli</italic> YccA, <italic>B. subtilis</italic> YetJ, and human BI-1. The nature of the cytoprotective role of eukaryotic BI-1 is its function in suppressing cell death, which is not limited to apoptosis and induced by multiple stimuli (<xref ref-type="bibr" rid="B21">Huckelhoven, 2004</xref>). It has been hypothesized that BI-1 may be of ancient bacterial origin with a conserved cell viability-associated role (<xref ref-type="bibr" rid="B21">Huckelhoven, 2004</xref>; <xref ref-type="bibr" rid="B19">Henke et al., 2011</xref>). Our finding that BrBI impairs <italic>B. suis</italic> S2 viability supports this hypothesis. However, the mechanisms of cell death differ between bacteria and eukaryotes, and bacteria do not go through apoptosis. Thus, the mechanism underlying the bacterial cytoprotective role of BrBI remains unclear.</p>
<p>Scanning electron microscopy analysis showed that the &#x0394;<italic>brbI</italic> strain displayed defective cell division, which may have contributed to the <italic>brbI</italic> deletion-induced defects in <italic>B. suis</italic> S2 growth and survival. Our integrated proteomic and transcriptomic analyses revealed that <italic>brbI</italic> deletion did not affect the expressions of FtsA, FtsZ, FtsE, FtsK, FtsW, MinC, MinD, and MinE at either the mRNA or protein levels but specifically increased the protein levels of FtsB, FtsI, FtsL, and FtsQ. Because the FtsQBL complex acts as a scaffold to link the early and late divisome components, these up-regulated proteins may trigger redundant cytokinesis septum formation, subsequently leading to defective cell division (<xref ref-type="bibr" rid="B22">Jain et al., 2018</xref>). However, how BrBI regulates division-associated proteins remains unclear. Although defective cell division may be the molecular basis for the suppressed growth and viability of the &#x0394;<italic>brbI</italic> strain, whether BrBI has other mechanisms that modulate <italic>B. suis</italic> S2 viability and/or growth remains unclear.</p>
<p>The integrated proteomic and transcriptomic analyses also showed that <italic>brbI</italic> deletion significantly affected the expression of membrane-associated protein-encoding genes at the mRNA and/or protein levels, such as in the outer membrane protein-, inner membrane protein-, and transporter-associated protein-encoding genes, suggesting a crucial role of BrBI in the <italic>B. suis</italic> S2 membrane properties. The Gram-negative bacterial membrane is multi layered, comprising a symmetric phospholipid bilayer termed the inner membrane, an asymmetric outer membrane containing phospholipids and LPS, and an intermediate aqueous periplasm containing a thin peptidoglycan layer (<xref ref-type="bibr" rid="B20">Hinz et al., 2011</xref>; <xref ref-type="bibr" rid="B28">Lehman and Grabowicz, 2019</xref>). Regarding integral membrane proteins, inner membrane proteins span the membrane as &#x03B1;-helices, while most OMPs contain a &#x03B2;-barrel structure (<xref ref-type="bibr" rid="B41">Sperandeo et al., 2019</xref>). Our results showed that <italic>brbI</italic> deletion significantly affected two key proteins, BamA and LptD, which were down-regulated at only the protein level. In Gram-negative bacteria, Lpt, Bam, and Lol are the three pathways required to respectively transport LPS, &#x03B2;-barrel proteins, and lipoproteins across the periplasm to the outer membrane (<xref ref-type="bibr" rid="B41">Sperandeo et al., 2019</xref>). Thus, down-regulation of BamA and LptD by <italic>brbI</italic> deletion suggested a defective assembly of OMPs and LPS in the &#x0394;<italic>brbI</italic> strain. Consistently, seven OMPs were significantly down-regulated. OMP2a, OMP2b, OMP22, and OMP31 were &#x03B2;-barrel containing proteins, which were affected by <italic>brbI</italic> deletion at only the protein level, possibly indicating a link with BamA. Two other key proteins, HflK and HflC, were also significantly affected by <italic>brbI</italic> deletion at only the protein level. Because HflK/C is a negative regulator of FtsH, a critical membrane quality-control mechanism, up-regulation of HflK/C suggested that FtsH may be regulated to alter membrane proteins. However, this hypothesis requires further research. This raises other critical questions. How does BrBI regulate BamA, LptD, HflC, and HflC? Are there other mechanisms by which BrBI regulates membrane proteins? Is there an interaction between BamA and HflK/C? As FtsB, FtsI, FtsL, and FtsQ are also membrane proteins, are they regulated by BamA or HflK/C? Future studies are needed to answer these questions.</p>
<p><italic>brbI</italic> deletion decreased the resistance of <italic>B. suis</italic> S2 to acid, hydrogen peroxide, polymyxin B, and lincomycin. Although the mechanisms of these stressors in <italic>Brucella</italic> are unclear, the integrated analysis provides insight into how these stressors are affected by <italic>brbI</italic> deletion. First, because the cell envelope is the first layer by which bacteria sense and respond to extracellular stresses, <italic>brbI</italic> deletion-induced membrane disorder may contribute to the increased sensitivity of <italic>B. suis</italic> S2 to stresses. <italic>brbI</italic> deletion significantly downregulated universal stress protein (WP_006191191.1), which may also have led to the decreased resistance of <italic>B. suis</italic> S2 to multiple stresses. Specifically, <italic>brbI</italic> deletion downregulated a series of multidrug efflux RND transporter subunits (WP_011068960.1, WP_004690093.1, and WP_006189806.1), which may participate in reducing antibiotic resistances. Additionally, the acid-activated periplasmic chaperone, HdeA (WP_006191791.1), was significantly decreased in the &#x0394;<italic>brbI</italic> strain, which may contribute to reducing acid resistance.</p>
<p>As an inner membrane protein, how does BrBI play such an extensive regulatory role? In bacteria, FtsH is the sole ATP-dependent metalloprotease embedded in the plasma membrane (<xref ref-type="bibr" rid="B38">Saikawa et al., 2004</xref>) and plays a crucial role in quality control by degrading unneeded or damaged membrane proteins and cytosolic proteins (<xref ref-type="bibr" rid="B2">Bieniossek et al., 2006</xref>). Although the degradation mechanism of FtsH remains unclear, a molecular architecture-based model has been proposed. FtsH forms a crown-shaped hexameric molecule consisting of two rings; one is built by the C-terminal protease domains; the other is formed by the AAA (ATPases associated with diverse cellular activities) domains facing the cytosolic leaflet of the membrane (<xref ref-type="bibr" rid="B2">Bieniossek et al., 2006</xref>). After the AAA domains bind a recognition tag, ATP hydrolysis leads to translocation of the target polypeptide into the interior of the FtsH &#x201C;crown,&#x201D; followed by proteolysis (<xref ref-type="bibr" rid="B2">Bieniossek et al., 2006</xref>). HflK and HflC are cytoplasmic membrane proteins and form a complex (<xref ref-type="bibr" rid="B24">Kihara et al., 1996</xref>). In <italic>E. coli</italic>, the HflK/C complex directly interacts with FtsH to suppress degradation of several substrates, including SecY, which is a proteolytic substrate of FtsH that shares the same recognition site with YccA (<xref ref-type="bibr" rid="B24">Kihara et al., 1996</xref>, <xref ref-type="bibr" rid="B43">1998</xref>). It has been reported that <italic>E. coli</italic> YccA can cross-link with both FtsH and HflK/C, and these four proteins subsequently form an exceptionally large complex in the plasma membrane (<xref ref-type="bibr" rid="B34">Kihara et al., 1998</xref>, <xref ref-type="bibr" rid="B23">1999</xref>; <xref ref-type="bibr" rid="B7">Saikawa et al., 2004</xref>). YccA11 is a mutant lacking eight amino acids in the N-terminal hydrophilic region of YccA, which can still interact with HflK/C and FtsH but cannot be degraded by FtsH (<xref ref-type="bibr" rid="B1">Kihara et al., 1998</xref>, <xref ref-type="bibr" rid="B7">1999</xref>). YccA11 mutation can stabilize excess SecY as does an HflK/C mutation, but deleting the intact YccA cannot accelerate SecY degradation as HflK/C deletion does (<xref ref-type="bibr" rid="B8">Kihara et al., 1998</xref>). Moreover, overexpression of HflK/C significantly increased the cross-linking efficiency between YccA11 and FtsH. Kihara et al. concluded that HflK/C was required for YccA to interact with FtsH, and HflK/C acted as a negative regulator for FtsH to degrade specific substrates. Further, the authors speculated that negative regulation of HflK/C for FtsH may attribute to binding of HflK/C and FtsH substrates, which restrained the substrates accessing to the protease domains of FtsH (<xref ref-type="bibr" rid="B15">Kihara et al., 1998</xref>; <xref ref-type="bibr" rid="B35">Bieniossek et al., 2006</xref>). Combined with the molecular architecture-based model of FtsH proteolysis, the mechanism of HflK/C may be that the HflK/C complex combines with FtsH to increase the spatial proximity of HflK/C and FtsH substrates, then its acting-domain binds with FtsH substrates, restraining the substrates to be translocated into the interior of the FtsH &#x201C;crown.&#x201D; As BrBI is homologous to <italic>E. coli</italic> YccA, and FtsH, HflK, and HflC also have chromosomally encoded homologs in <italic>Brucella</italic> spp., we speculate that the possible mechanism of BrBI is that: (i) <italic>brbI</italic> deletion up-regulates HflK and HflC, and the increased HflK/C may suppress FtsH to degrade specific membrane and/or cytosolic proteins, possibly including proteases, transcriptional factors, and other regulatory or functional proteins (such as OMPs, transporter subunits, response regulators, and chaperone proteins); (ii) the accumulated regulatory proteins may further trigger downstream cascade regulation of multiple genes at the mRNA and/or protein levels, ultimately resulting in functional or structural protein alterations; and (iii) alterations of the functional or structural proteins are eventually reflected in their phenotypes (<xref ref-type="fig" rid="F11">Figure 11A</xref>).</p>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption><p>Schematic representation of the major findings. <bold>(A)</bold> The speculated mechanism of BrBI. <bold>(B)</bold> Effects of <italic>brbI</italic> deletion on outer membrane proteins, cell division proteins, and transporter proteins. &#x201C;&#x2013;&#x201D; Indicates negative transcriptional regulators, &#x201C;+&#x201D; indicates positive transcriptional regulators, the red arrow indicates the regulation direction of the protein, the double arrow indicates indirect regulation.</p></caption>
<graphic xlink:href="fmicb-12-632095-g011.tif"/>
</fig>
<p>In conclusion, mutagenesis and integrated proteomic and transcriptomic analyses revealed that BrBI plays crucial roles in <italic>B. suis</italic> S2 membrane homeostasis, cell viability, and cell division of. Besides, we speculated that as a membrane protein, BrBI may initially regulate specific membrane proteins (such as HflK, HflC, and BamA) at the protein level, then these proteins further regulate multiple gene expressions (may including regulatory, functional, or structural genes) at the mRNA and/or protein levels, and these cascade regulations are eventually reflected on the physiology of <italic>B. suis</italic> S2 (<xref ref-type="fig" rid="F11">Figure 11B</xref>). Nevertheless, because the prokaryotic members of BI-1 family proteins have rarely been studied, the detailed regulatory mechanism of BrBI requires more thorough investigations.</p>
</sec>
<sec id="S5">
<title>Data Availability Statement</title>
<p>RNA-seq raw data can be found in the NCBI database, the data was assigned an accession number of <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="GSE160911">GSE160911</ext-link>. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium (<ext-link ext-link-type="uri" xlink:href="http://proteomecentral.proteomexchange.org">http://proteomecentral.proteomexchange.org</ext-link>) via the iProX partner repository with the dataset identifier IPX0002588001/PXD022537.</p>
</sec>
<sec id="S6">
<title>Author Contributions</title>
<p>GZ, YJ, and AW conceived and designed the experiments. GZ, FlZ, LC, PQ, JL, FjZ, and LT performed the experiments and analyzed the data. GZ, YJ, and AW wrote and revised the manuscript. DZ, PL, HC, KT, and WL revised the manuscript. All the authors have read and approved the final manuscript.</p>
</sec>
<sec sec-type="COI-statement" id="conf1">
<title>Conflict of Interest</title>
<p>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.</p>
</sec>
</body>
<back>
<fn-group>
<fn fn-type="financial-disclosure">
<p><bold>Funding.</bold> This study was supported by the National Key R&#x0026;D Program of China (No. 2018YFD0500900) and the National Natural Science Foundation of China (No. 31672584).</p>
</fn>
</fn-group>
<ack>
<p>We are grateful to the LC-Bio Technologies (Hangzhou) Co., Ltd. for assisting in RNA-sequencing, LC-MS/MS and data analysis.</p>
</ack>
<sec id="S9" sec-type="supplementary material"><title>Supplementary Material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fmicb.2021.632095/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fmicb.2021.632095/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Data_Sheet_1.zip" id="SM1" mimetype="application/zip" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Altamirano-Silva</surname> <given-names>P.</given-names></name> <name><surname>Meza-Torres</surname> <given-names>J.</given-names></name> <name><surname>Castillo-Zeledon</surname> <given-names>A.</given-names></name> <name><surname>Ruiz-Villalobos</surname> <given-names>N.</given-names></name> <name><surname>Zuniga-Pereira</surname> <given-names>A. M.</given-names></name> <name><surname>Chacon-Diaz</surname> <given-names>C.</given-names></name><etal/></person-group> (<year>2018</year>). <article-title>Brucella abortus Senses the Intracellular Environment through the BvrR/BvrS Two-Component System, Which Allows B. abortus To Adapt to Its Replicative Niche.</article-title> <source><italic>Infect. Immun.</italic></source> <volume>86</volume> <fpage>713</fpage>&#x2013;<lpage>717 e</lpage>. <pub-id pub-id-type="doi">10.1128/IAI.00713-17</pub-id> <pub-id pub-id-type="pmid">29378792</pub-id></citation></ref>
<ref id="B2"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bieniossek</surname> <given-names>C.</given-names></name> <name><surname>Schalch</surname> <given-names>T.</given-names></name> <name><surname>Bumann</surname> <given-names>M.</given-names></name> <name><surname>Meister</surname> <given-names>M.</given-names></name> <name><surname>Meier</surname> <given-names>R.</given-names></name> <name><surname>Baumann</surname> <given-names>U.</given-names></name></person-group> (<year>2006</year>). <article-title>The molecular architecture of the metalloprotease FtsH.</article-title> <source><italic>Proc. Natl. Acad. Sci. U. S. A.</italic></source> <volume>103</volume> <fpage>3066</fpage>&#x2013;<lpage>3071</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.0600031103</pub-id> <pub-id pub-id-type="pmid">16484367</pub-id></citation></ref>
<ref id="B3"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bittner</surname> <given-names>L. M.</given-names></name> <name><surname>Arends</surname> <given-names>J.</given-names></name> <name><surname>Narberhaus</surname> <given-names>F.</given-names></name></person-group> (<year>2017</year>). <article-title>When, how and why? Regulated proteolysis by the essential FtsH protease in <italic>Escherichia coli</italic>.</article-title> <source><italic>Biol. Chem.</italic></source> <volume>398</volume> <fpage>625</fpage>&#x2013;<lpage>635</lpage>. <pub-id pub-id-type="doi">10.1515/hsz-2016-0302</pub-id> <pub-id pub-id-type="pmid">28085670</pub-id></citation></ref>
<ref id="B4"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Boschiroli</surname> <given-names>M. L.</given-names></name> <name><surname>Ouahrani-Bettache</surname> <given-names>S.</given-names></name> <name><surname>Foulongne</surname> <given-names>V.</given-names></name> <name><surname>Michaux-Charachon</surname> <given-names>S.</given-names></name> <name><surname>Bourg</surname> <given-names>G.</given-names></name> <name><surname>Allardet-Servent</surname> <given-names>A.</given-names></name><etal/></person-group> (<year>2002</year>). <article-title>The Brucella suis virB operon is induced intracellularly in macrophages.</article-title> <source><italic>Proc. Natl. Acad. Sci. U. S. A.</italic></source> <volume>99</volume> <fpage>1544</fpage>&#x2013;<lpage>1549</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.032514299</pub-id> <pub-id pub-id-type="pmid">11830669</pub-id></citation></ref>
<ref id="B5"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Breukink</surname> <given-names>E.</given-names></name></person-group> (<year>2009</year>). <article-title>Cell biology. Lethal traffic jam.</article-title> <source><italic>Science</italic></source> <volume>325</volume> <fpage>684</fpage>&#x2013;<lpage>685</lpage>. <pub-id pub-id-type="doi">10.1126/science.1178424</pub-id> <pub-id pub-id-type="pmid">19661407</pub-id></citation></ref>
<ref id="B6"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bultynck</surname> <given-names>G.</given-names></name> <name><surname>Kiviluoto</surname> <given-names>S.</given-names></name> <name><surname>Henke</surname> <given-names>N.</given-names></name> <name><surname>Ivanova</surname> <given-names>H.</given-names></name> <name><surname>Schneider</surname> <given-names>L.</given-names></name> <name><surname>Rybalchenko</surname> <given-names>V.</given-names></name><etal/></person-group> (<year>2012</year>). <article-title>The C terminus of Bax inhibitor-1 forms a Ca2+-permeable channel pore.</article-title> <source><italic>J. Biol. Chem.</italic></source> <volume>287</volume> <fpage>2544</fpage>&#x2013;<lpage>2557</lpage>. <pub-id pub-id-type="doi">10.1074/jbc.M111.275354</pub-id> <pub-id pub-id-type="pmid">22128171</pub-id></citation></ref>
<ref id="B7"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Byndloss</surname> <given-names>M. X.</given-names></name> <name><surname>Tsolis</surname> <given-names>R. M.</given-names></name></person-group> (<year>2016</year>). <article-title>Brucella spp. Virulence Factors and Immunity.</article-title> <source><italic>Annu. Rev. Anim. Biosci.</italic></source> <volume>4</volume> <fpage>111</fpage>&#x2013;<lpage>127</lpage>. <pub-id pub-id-type="doi">10.1146/annurev-animal-021815-111326</pub-id> <pub-id pub-id-type="pmid">26734887</pub-id></citation></ref>
<ref id="B8"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Celli</surname> <given-names>J.</given-names></name></person-group> (<year>2019</year>). <article-title>The Intracellular Life Cycle of Brucella spp.</article-title> <source><italic>Microbiol. Spectr.</italic></source> <volume>7</volume>:<issue>30848234</issue>. <pub-id pub-id-type="doi">10.1128/microbiolspec.BAI-0006-2019</pub-id> <pub-id pub-id-type="pmid">30848234</pub-id></citation></ref>
<ref id="B9"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Celli</surname> <given-names>J.</given-names></name> <name><surname>de Chastellier</surname> <given-names>C.</given-names></name> <name><surname>Franchini</surname> <given-names>D. M.</given-names></name> <name><surname>Pizarro-Cerda</surname> <given-names>J.</given-names></name> <name><surname>Moreno</surname> <given-names>E.</given-names></name> <name><surname>Gorvel</surname> <given-names>A. P.</given-names></name></person-group> (<year>2003</year>). <article-title>Brucella evades macrophage killing via VirB-dependent sustained interactions with the endoplasmic reticulum.</article-title> <source><italic>J. Exp. Med.</italic></source> <volume>198</volume> <fpage>545</fpage>&#x2013;<lpage>556</lpage>. <pub-id pub-id-type="doi">10.1084/jem.20030088</pub-id> <pub-id pub-id-type="pmid">12925673</pub-id></citation></ref>
<ref id="B10"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chae</surname> <given-names>H. J.</given-names></name> <name><surname>Kim</surname> <given-names>H. R.</given-names></name> <name><surname>Xu</surname> <given-names>C. Y.</given-names></name> <name><surname>Bailly-Maitre</surname> <given-names>B.</given-names></name> <name><surname>Krajewska</surname> <given-names>M.</given-names></name> <name><surname>Krajewski</surname> <given-names>S.</given-names></name><etal/></person-group> (<year>2004</year>). <article-title>BI-1 regulates an apoptosis pathway linked to endoplasmic reticulum stress.</article-title> <source><italic>Mol. Cell</italic></source> <volume>15</volume> <fpage>355</fpage>&#x2013;<lpage>366</lpage>. <pub-id pub-id-type="doi">10.1016/j.molcel.2004.06.038</pub-id> <pub-id pub-id-type="pmid">15304216</pub-id></citation></ref>
<ref id="B11"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chang</surname> <given-names>Y.</given-names></name> <name><surname>Bruni</surname> <given-names>R.</given-names></name> <name><surname>Kloss</surname> <given-names>B.</given-names></name> <name><surname>Assur</surname> <given-names>Z.</given-names></name> <name><surname>Kloppmann</surname> <given-names>E.</given-names></name> <name><surname>Rost</surname> <given-names>B.</given-names></name><etal/></person-group> (<year>2014</year>). <article-title>Structural basis for a pH-sensitive calcium leak across membranes.</article-title> <source><italic>Science</italic></source> <volume>344</volume> <fpage>1131</fpage>&#x2013;<lpage>1135</lpage>. <pub-id pub-id-type="doi">10.1126/science.1252043</pub-id> <pub-id pub-id-type="pmid">24904158</pub-id></citation></ref>
<ref id="B12"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Choi</surname> <given-names>Y.</given-names></name> <name><surname>Kim</surname> <given-names>J.</given-names></name> <name><surname>Yoon</surname> <given-names>H. J.</given-names></name> <name><surname>Jin</surname> <given-names>K. S.</given-names></name> <name><surname>Ryu</surname> <given-names>S.</given-names></name> <name><surname>Lee</surname> <given-names>H. H.</given-names></name></person-group> (<year>2018</year>). <article-title>Structural Insights into the FtsQ/FtsB/FtsL Complex, a Key Component of the Divisome.</article-title> <source><italic>Sci. Rep.</italic></source> <volume>8</volume>:<issue>18061</issue>. <pub-id pub-id-type="doi">10.1038/s41598-018-36001-2</pub-id> <pub-id pub-id-type="pmid">30584256</pub-id></citation></ref>
<ref id="B13"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Condon</surname> <given-names>S. G. F.</given-names></name> <name><surname>Mahbuba</surname> <given-names>D.-A.</given-names></name> <name><surname>Armstrong</surname> <given-names>C. R.</given-names></name> <name><surname>Diaz-Vazquez</surname> <given-names>G.</given-names></name> <name><surname>Craven</surname> <given-names>S. J.</given-names></name> <name><surname>LaPointe</surname> <given-names>L. M.</given-names></name><etal/></person-group> (<year>2018</year>). <article-title>The FtsLB subcomplex of the bacterial divisome is a tetramer with an uninterrupted FtsL helix linking the transmembrane and periplasmic regions.</article-title> <source><italic>J. Biol. Chem.</italic></source> <volume>293</volume> <fpage>1623</fpage>&#x2013;<lpage>1641</lpage>. <pub-id pub-id-type="doi">10.1074/jbc.RA117.000426</pub-id> <pub-id pub-id-type="pmid">29233891</pub-id></citation></ref>
<ref id="B14"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dajkovic</surname> <given-names>A.</given-names></name> <name><surname>Lutkenhaus</surname> <given-names>J.</given-names></name></person-group> (<year>2006</year>). <article-title>Z ring as executor of bacterial cell division.</article-title> <source><italic>J. Mol. Microbiol. Biotechnol.</italic></source> <volume>11</volume> <fpage>140</fpage>&#x2013;<lpage>151</lpage>. <pub-id pub-id-type="doi">10.1159/000094050</pub-id> <pub-id pub-id-type="pmid">16983191</pub-id></citation></ref>
<ref id="B15"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Di</surname> <given-names>D. D.</given-names></name> <name><surname>Jiang</surname> <given-names>H.</given-names></name> <name><surname>Tian</surname> <given-names>L. L.</given-names></name> <name><surname>Kang</surname> <given-names>J. L.</given-names></name> <name><surname>Zhang</surname> <given-names>W.</given-names></name> <name><surname>Yi</surname> <given-names>X. P.</given-names></name><etal/></person-group> (<year>2016</year>). <article-title>Comparative genomic analysis between newly sequenced Brucella suis Vaccine Strain S2 and the Virulent Brucella suis Strain 1330.</article-title> <source><italic>BMC Genom.</italic></source> <volume>17</volume>:<issue>741</issue>. <pub-id pub-id-type="doi">10.1186/s12864-016-3076-5</pub-id> <pub-id pub-id-type="pmid">27645563</pub-id></citation></ref>
<ref id="B16"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gamboa-Tuz</surname> <given-names>S. D.</given-names></name> <name><surname>Pereira-Santana</surname> <given-names>A.</given-names></name> <name><surname>Zhao</surname> <given-names>T.</given-names></name> <name><surname>Schranz</surname> <given-names>M. E.</given-names></name> <name><surname>Castano</surname> <given-names>E.</given-names></name> <name><surname>Rodriguez-Zapata</surname> <given-names>L. C.</given-names></name></person-group> (<year>2018</year>). <article-title>New insights into the phylogeny of the TMBIM superfamily across the tree of life: Comparative genomics and synteny networks reveal independent evolution of the BI and LFG families in plants.</article-title> <source><italic>Mol. Phylogenet. Evol.</italic></source> <volume>126</volume> <fpage>266</fpage>&#x2013;<lpage>278</lpage>. <pub-id pub-id-type="doi">10.1016/j.ympev.2018.04.032</pub-id> <pub-id pub-id-type="pmid">29702215</pub-id></citation></ref>
<ref id="B17"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gomez</surname> <given-names>G.</given-names></name> <name><surname>Adams</surname> <given-names>L. G.</given-names></name> <name><surname>Rice-Ficht</surname> <given-names>A.</given-names></name> <name><surname>Ficht</surname> <given-names>T. A.</given-names></name></person-group> (<year>2013</year>). <article-title>Host-Brucella interactions and the Brucella genome as tools for subunit antigen discovery and immunization against brucellosis.</article-title> <source><italic>Front. Cell. Infect. Microbiol</italic></source> <volume>3</volume>:<issue>17</issue>. <pub-id pub-id-type="doi">10.3389/fcimb.2013.00017</pub-id> <pub-id pub-id-type="pmid">23720712</pub-id></citation></ref>
<ref id="B18"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Goolab</surname> <given-names>S.</given-names></name> <name><surname>Roth</surname> <given-names>R. L.</given-names></name> <name><surname>van Heerden</surname> <given-names>H.</given-names></name> <name><surname>Crampton</surname> <given-names>M. C.</given-names></name></person-group> (<year>2015</year>). <article-title>Analyzing the molecular mechanism of lipoprotein localization in Brucella.</article-title> <source><italic>Front. Microbiol.</italic></source> <volume>6</volume>:<issue>1189</issue>. <pub-id pub-id-type="doi">10.3389/fmicb.2015.01189</pub-id> <pub-id pub-id-type="pmid">26579096</pub-id></citation></ref>
<ref id="B19"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Henke</surname> <given-names>N.</given-names></name> <name><surname>Lisak</surname> <given-names>D. A.</given-names></name> <name><surname>Schneider</surname> <given-names>L.</given-names></name> <name><surname>Habicht</surname> <given-names>J.</given-names></name> <name><surname>Pergande</surname> <given-names>M.</given-names></name> <name><surname>Methner</surname> <given-names>A.</given-names></name></person-group> (<year>2011</year>). <article-title>The ancient cell death suppressor BAX inhibitor-1.</article-title> <source><italic>Cell. Calcium</italic></source> <volume>50</volume> <fpage>251</fpage>&#x2013;<lpage>260</lpage>. <pub-id pub-id-type="doi">10.1016/j.ceca.2011.05.005</pub-id> <pub-id pub-id-type="pmid">21663964</pub-id></citation></ref>
<ref id="B20"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hinz</surname> <given-names>A.</given-names></name> <name><surname>Lee</surname> <given-names>S.</given-names></name> <name><surname>Jacoby</surname> <given-names>K.</given-names></name> <name><surname>Manoil</surname> <given-names>C.</given-names></name></person-group> (<year>2011</year>). <article-title>Membrane proteases and aminoglycoside antibiotic resistance.</article-title> <source><italic>J. Bacteriol.</italic></source> <volume>193</volume> <fpage>4790</fpage>&#x2013;<lpage>4797</lpage>. <pub-id pub-id-type="doi">10.1128/JB.05133-11</pub-id> <pub-id pub-id-type="pmid">21764915</pub-id></citation></ref>
<ref id="B21"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Huckelhoven</surname> <given-names>R.</given-names></name></person-group> (<year>2004</year>). <article-title>BAX Inhibitor-1, an ancient cell death suppressor in animals and plants with prokaryotic relatives.</article-title> <source><italic>Apoptosis</italic></source> <volume>9</volume> <fpage>299</fpage>&#x2013;<lpage>307</lpage>. <pub-id pub-id-type="doi">10.1023/b:appt.0000025806.71000.1c</pub-id></citation></ref>
<ref id="B22"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jain</surname> <given-names>P.</given-names></name> <name><surname>Malakar</surname> <given-names>B.</given-names></name> <name><surname>Khan</surname> <given-names>M. Z.</given-names></name> <name><surname>Lochab</surname> <given-names>S.</given-names></name> <name><surname>Singh</surname> <given-names>A.</given-names></name> <name><surname>Nandicoori</surname> <given-names>V. K.</given-names></name></person-group> (<year>2018</year>). <article-title>Delineating FtsQ-mediated regulation of cell division in Mycobacterium tuberculosis.</article-title> <source><italic>J. Biol. Chem.</italic></source> <volume>293</volume> <fpage>12331</fpage>&#x2013;<lpage>12349</lpage>. <pub-id pub-id-type="doi">10.1074/jbc.RA118.003628</pub-id> <pub-id pub-id-type="pmid">29903917</pub-id></citation></ref>
<ref id="B23"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ke</surname> <given-names>Y. H.</given-names></name> <name><surname>Wang</surname> <given-names>Y. F.</given-names></name> <name><surname>Li</surname> <given-names>W. F.</given-names></name> <name><surname>Chen</surname> <given-names>Z. L.</given-names></name></person-group> (<year>2015</year>). <article-title>Type IV secretion system of Brucella spp. and its effectors.</article-title> <source><italic>Front. Cell. Infect. Microbiol.</italic></source> <volume>5</volume>:<issue>72</issue>. <pub-id pub-id-type="doi">10.3389/famb.2015.00072</pub-id></citation></ref>
<ref id="B24"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kihara</surname> <given-names>A.</given-names></name> <name><surname>Akiyama</surname> <given-names>Y.</given-names></name> <name><surname>Ito</surname> <given-names>K.</given-names></name></person-group> (<year>1996</year>). <article-title>A protease complex in the <italic>Escherichia coli</italic> plasma membrane: HflKC (HflA) forms a complex with FtsH (HflB), regulating its proteolytic activity against SecY.</article-title> <source><italic>EMBO J.</italic></source> <volume>15</volume> <fpage>6122</fpage>&#x2013;<lpage>6131</lpage>. <pub-id pub-id-type="doi">10.1002/j.1460-2075.1996.tb01000.x</pub-id></citation></ref>
<ref id="B25"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kihara</surname> <given-names>A.</given-names></name> <name><surname>Akiyama</surname> <given-names>Y.</given-names></name> <name><surname>Ito</surname> <given-names>K.</given-names></name></person-group> (<year>1998</year>). <article-title>Different pathways for protein degradation by the FtsH/HflKC membrane-embedded protease complex: an implication from the interference by a mutant form of a new substrate protein.</article-title> <source><italic>YCCA. J. Mol. Biol.</italic></source> <volume>279</volume> <fpage>175</fpage>&#x2013;<lpage>188</lpage>. <pub-id pub-id-type="doi">10.1006/jmbi.1998.1781</pub-id> <pub-id pub-id-type="pmid">9636708</pub-id></citation></ref>
<ref id="B26"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kihara</surname> <given-names>A.</given-names></name> <name><surname>Akiyama</surname> <given-names>Y.</given-names></name> <name><surname>Ito</surname> <given-names>K.</given-names></name></person-group> (<year>1999</year>). <article-title>Dislocation of membrane proteins in FtsH-mediated proteolysis.</article-title> <source><italic>EMBO J.</italic></source> <volume>18</volume> <fpage>2970</fpage>&#x2013;<lpage>2981</lpage>. <pub-id pub-id-type="doi">10.1093/emboj/18.11.2970</pub-id> <pub-id pub-id-type="pmid">10357810</pub-id></citation></ref>
<ref id="B27"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Langklotz</surname> <given-names>S.</given-names></name> <name><surname>Baumann</surname> <given-names>U.</given-names></name> <name><surname>Narberhaus</surname> <given-names>F.</given-names></name></person-group> (<year>2012</year>). <article-title>Structure and function of the bacterial AAA protease FtsH.</article-title> <source><italic>Biochim. Biophys. Acta</italic></source> <volume>1823</volume> <fpage>40</fpage>&#x2013;<lpage>48</lpage>. <pub-id pub-id-type="doi">10.1016/j.bbamcr.2011.08.015</pub-id> <pub-id pub-id-type="pmid">21925212</pub-id></citation></ref>
<ref id="B28"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lehman</surname> <given-names>K. M.</given-names></name> <name><surname>Grabowicz</surname> <given-names>M.</given-names></name></person-group> (<year>2019</year>). <article-title>Countering Gram-Negative Antibiotic Resistance: Recent Progress in Disrupting the Outer Membrane with Novel Therapeutics.</article-title> <source><italic>Antibiot</italic></source> <volume>8</volume>:<issue>163</issue>. <pub-id pub-id-type="doi">10.3390/antibiotics8040163</pub-id> <pub-id pub-id-type="pmid">31554212</pub-id></citation></ref>
<ref id="B29"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lisbona</surname> <given-names>F.</given-names></name> <name><surname>Rojas-Rivera</surname> <given-names>D.</given-names></name> <name><surname>Thielen</surname> <given-names>P.</given-names></name> <name><surname>Zamorano</surname> <given-names>S.</given-names></name> <name><surname>Todd</surname> <given-names>D.</given-names></name> <name><surname>Martinon</surname> <given-names>F.</given-names></name><etal/></person-group> (<year>2009</year>). <article-title>BAX inhibitor-1 is a negative regulator of the ER stress sensor IRE1alpha.</article-title> <source><italic>Mol Cell</italic></source> <volume>33</volume> <fpage>679</fpage>&#x2013;<lpage>691</lpage>. <pub-id pub-id-type="doi">10.1016/j.molcel.2009.02.017</pub-id> <pub-id pub-id-type="pmid">19328063</pub-id></citation></ref>
<ref id="B30"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname> <given-names>Q.</given-names></name></person-group> (<year>2017</year>). <article-title>TMBIM-mediated Ca(2+) homeostasis and cell death.</article-title> <source><italic>Biochim. Biophys. Acta Mol. Cell Res.</italic></source> <volume>1864</volume> <fpage>850</fpage>&#x2013;<lpage>857</lpage>. <pub-id pub-id-type="doi">10.1016/j.bbamcr.2016.12.023</pub-id> <pub-id pub-id-type="pmid">28064000</pub-id></citation></ref>
<ref id="B31"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname> <given-names>Q.</given-names></name> <name><surname>Hu</surname> <given-names>M.</given-names></name> <name><surname>Yeo</surname> <given-names>W. S.</given-names></name> <name><surname>He</surname> <given-names>L.</given-names></name> <name><surname>Li</surname> <given-names>T.</given-names></name> <name><surname>Zhu</surname> <given-names>Y.</given-names></name><etal/></person-group> (<year>2017</year>). <article-title>Rewiring of the FtsH regulatory network by a single nucleotide change in saeS of Staphylococcus aureus.</article-title> <source><italic>Sci. Rep.</italic></source> <volume>7</volume>:<issue>8456</issue>. <pub-id pub-id-type="doi">10.1038/s41598-017-08774-5</pub-id> <pub-id pub-id-type="pmid">28814746</pub-id></citation></ref>
<ref id="B32"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lutkenhaus</surname> <given-names>J.</given-names></name></person-group> (<year>2007</year>). <article-title>Assembly dynamics of the bacterial MinCDE system and spatial regulation of the Z ring.</article-title> <source><italic>Annu. Rev. Biochem.</italic></source> <volume>76</volume> <fpage>539</fpage>&#x2013;<lpage>562</lpage>. <pub-id pub-id-type="doi">10.1146/annurev.biochem.75.103004.142652</pub-id> <pub-id pub-id-type="pmid">17328675</pub-id></citation></ref>
<ref id="B33"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Madeira</surname> <given-names>F.</given-names></name> <name><surname>Park</surname> <given-names>Y. M.</given-names></name> <name><surname>Lee</surname> <given-names>J.</given-names></name> <name><surname>Buso</surname> <given-names>N.</given-names></name> <name><surname>Gur</surname> <given-names>T.</given-names></name> <name><surname>Madhusoodanan</surname> <given-names>N.</given-names></name><etal/></person-group> (<year>2019</year>). <article-title>The EMBL-EBI search and sequence analysis tools APIs in 2019.</article-title> <source><italic>Nucleic Acids Res.</italic></source> <volume>47</volume> <fpage>W636</fpage>&#x2013;<lpage>W641</lpage>. <pub-id pub-id-type="doi">10.1093/nar/gkz268</pub-id> <pub-id pub-id-type="pmid">30976793</pub-id></citation></ref>
<ref id="B34"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Martinez-Nunez</surname> <given-names>C.</given-names></name> <name><surname>Altamirano-Silva</surname> <given-names>P.</given-names></name> <name><surname>Alvarado-Guillen</surname> <given-names>F.</given-names></name> <name><surname>Moreno</surname> <given-names>E.</given-names></name> <name><surname>Guzman-Verri</surname> <given-names>C.</given-names></name> <name><surname>Chaves-Olarte</surname> <given-names>E.</given-names></name></person-group> (<year>2010</year>). <article-title>The Two-Component System BvrR/BvrS Regulates the Expression of the Type IV Secretion System VirB in Brucella abortus.</article-title> <source><italic>J. Bacteriol.</italic></source> <volume>192</volume> <fpage>5603</fpage>&#x2013;<lpage>5608</lpage>. <pub-id pub-id-type="doi">10.1128/Jb.00567-10</pub-id> <pub-id pub-id-type="pmid">20833814</pub-id></citation></ref>
<ref id="B35"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>McQuiston</surname> <given-names>J. R.</given-names></name> <name><surname>Schurig</surname> <given-names>G. G.</given-names></name> <name><surname>Sriranganathan</surname> <given-names>N.</given-names></name> <name><surname>Boyle</surname> <given-names>S. M.</given-names></name></person-group> (<year>1995</year>). <article-title>Transformation of Brucella species with suicide and broad host-range plasmids.</article-title> <source><italic>Methods Mol. Biol.</italic></source> <volume>47</volume> <fpage>143</fpage>&#x2013;<lpage>148</lpage>. <pub-id pub-id-type="doi">10.1385/0-89603-310-4:143</pub-id></citation></ref>
<ref id="B36"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ramm</surname> <given-names>B.</given-names></name> <name><surname>Heermann</surname> <given-names>T.</given-names></name> <name><surname>Schwille</surname> <given-names>P.</given-names></name></person-group> (<year>2019</year>). <article-title>The E. coli MinCDE system in the regulation of protein patterns and gradients.</article-title> <source><italic>Cell. Mol. Life Sci.</italic></source> <volume>76</volume> <fpage>4245</fpage>&#x2013;<lpage>4273</lpage>. <pub-id pub-id-type="doi">10.1007/s00018-019-03218-x</pub-id> <pub-id pub-id-type="pmid">31317204</pub-id></citation></ref>
<ref id="B37"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Roy</surname> <given-names>A.</given-names></name> <name><surname>Kucukural</surname> <given-names>A.</given-names></name> <name><surname>Zhang</surname> <given-names>Y.</given-names></name></person-group> (<year>2010</year>). <article-title>I-TASSER: a unified platform for automated protein structure and function prediction.</article-title> <source><italic>Nat. Protoc.</italic></source> <volume>5</volume> <fpage>725</fpage>&#x2013;<lpage>738</lpage>. <pub-id pub-id-type="doi">10.1038/nprot.2010.5</pub-id> <pub-id pub-id-type="pmid">20360767</pub-id></citation></ref>
<ref id="B38"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Saikawa</surname> <given-names>N.</given-names></name> <name><surname>Akiyama</surname> <given-names>Y.</given-names></name> <name><surname>Ito</surname> <given-names>K.</given-names></name></person-group> (<year>2004</year>). <article-title>FtsH exists as an exceptionally large complex containing HflKC in the plasma membrane of <italic>Escherichia coli</italic>.</article-title> <source><italic>J. Struct. Biol.</italic></source> <volume>146</volume> <fpage>123</fpage>&#x2013;<lpage>129</lpage>. <pub-id pub-id-type="doi">10.1016/j.jsb.2003.09.020</pub-id> <pub-id pub-id-type="pmid">15037243</pub-id></citation></ref>
<ref id="B39"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Schumann</surname> <given-names>W.</given-names></name></person-group> (<year>1999</year>). <article-title>FtsH&#x2013;a single-chain charonin?</article-title> <source><italic>FEMS Microbiol. Rev.</italic></source> <volume>23</volume> <fpage>1</fpage>&#x2013;<lpage>11</lpage>. <pub-id pub-id-type="doi">10.1111/j.1574-6976.1999.tb00389.x</pub-id> <pub-id pub-id-type="pmid">10077851</pub-id></citation></ref>
<ref id="B40"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Shi</surname> <given-names>L.</given-names></name> <name><surname>Wang</surname> <given-names>Z.</given-names></name> <name><surname>Liu</surname> <given-names>X.</given-names></name> <name><surname>Li</surname> <given-names>M.</given-names></name> <name><surname>Zhang</surname> <given-names>S.</given-names></name> <name><surname>Song</surname> <given-names>X.</given-names></name></person-group> (<year>2018</year>). <article-title>Bax inhibitor-1 is required for resisting the Early Brain Injury induced by subarachnoid hemorrhage through regulating IRE1-JNK pathway.</article-title> <source><italic>Neurol. Res.</italic></source> <volume>40</volume> <fpage>189</fpage>&#x2013;<lpage>196</lpage>. <pub-id pub-id-type="doi">10.1080/01616412.2018.1424699</pub-id> <pub-id pub-id-type="pmid">29334839</pub-id></citation></ref>
<ref id="B41"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sperandeo</surname> <given-names>P.</given-names></name> <name><surname>Martorana</surname> <given-names>A. M.</given-names></name> <name><surname>Polissi</surname> <given-names>A.</given-names></name></person-group> (<year>2019</year>). <article-title>Lipopolysaccharide Biosynthesis and Transport to the Outer Membrane of Gram-Negative Bacteria.</article-title> <source><italic>Subcell. Biochem.</italic></source> <volume>92</volume> <fpage>9</fpage>&#x2013;<lpage>37</lpage>. <pub-id pub-id-type="doi">10.1007/978-3-030-18768-2_2</pub-id></citation></ref>
<ref id="B42"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Szwedziak</surname> <given-names>P.</given-names></name> <name><surname>Ghosal</surname> <given-names>D.</given-names></name></person-group> (<year>2017</year>). <article-title>FtsZ-ring Architecture and Its Control by MinCD.</article-title> <source><italic>Subcell. Biochem.</italic></source> <volume>84</volume> <fpage>213</fpage>&#x2013;<lpage>244</lpage>. <pub-id pub-id-type="doi">10.1007/978-3-319-53047-5_7</pub-id></citation></ref>
<ref id="B43"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Uzureau</surname> <given-names>S.</given-names></name> <name><surname>Godefroid</surname> <given-names>M.</given-names></name> <name><surname>Deschamps</surname> <given-names>C.</given-names></name> <name><surname>Lemaire</surname> <given-names>J.</given-names></name> <name><surname>De Bolle</surname> <given-names>X.</given-names></name> <name><surname>Letesson</surname> <given-names>J. J.</given-names></name></person-group> (<year>2007</year>). <article-title>Mutations of the quorum sensing-dependent regulator VjbR lead to drastic surface modifications in Brucella melitensis.</article-title> <source><italic>J. Bacteriol.</italic></source> <volume>189</volume> <fpage>6035</fpage>&#x2013;<lpage>6047</lpage>. <pub-id pub-id-type="doi">10.1128/JB.00265-07</pub-id> <pub-id pub-id-type="pmid">17557825</pub-id></citation></ref>
<ref id="B44"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>van Stelten</surname> <given-names>J.</given-names></name> <name><surname>Silva</surname> <given-names>F.</given-names></name> <name><surname>Belin</surname> <given-names>D.</given-names></name> <name><surname>Silhavy</surname> <given-names>T. J.</given-names></name></person-group> (<year>2009</year>). <article-title>Effects of antibiotics and a proto-oncogene homolog on destruction of protein translocator SecY.</article-title> <source><italic>Science</italic></source> <volume>325</volume> <fpage>753</fpage>&#x2013;<lpage>756</lpage>. <pub-id pub-id-type="doi">10.1126/science.1172221</pub-id> <pub-id pub-id-type="pmid">19661432</pub-id></citation></ref>
<ref id="B45"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Watanabe</surname> <given-names>N.</given-names></name> <name><surname>Lam</surname> <given-names>E.</given-names></name></person-group> (<year>2009</year>). <article-title>Bax Inhibitor-1, a Conserved Cell Death Suppressor, Is a Key Molecular Switch Downstream from a Variety of Biotic and Abiotic Stress Signals in Plants.</article-title> <source><italic>Int. J. Mol. Sci.</italic></source> <volume>10</volume> <fpage>3149</fpage>&#x2013;<lpage>3167</lpage>. <pub-id pub-id-type="doi">10.3390/ijms10073149</pub-id> <pub-id pub-id-type="pmid">19742129</pub-id></citation></ref>
<ref id="B46"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yang</surname> <given-names>J.</given-names></name> <name><surname>Zhang</surname> <given-names>Y.</given-names></name></person-group> (<year>2015</year>). <article-title>I-TASSER server: new development for protein structure and function predictions.</article-title> <source><italic>Nucleic Acids Res.</italic></source> <volume>43</volume> <fpage>W174</fpage>&#x2013;<lpage>W181</lpage>. <pub-id pub-id-type="doi">10.1093/nar/gkv342</pub-id> <pub-id pub-id-type="pmid">25883148</pub-id></citation></ref>
</ref-list><fn-group>
<fn id="footnote1">
<label>1</label>
<p><ext-link ext-link-type="uri" xlink:href="http://www.usadellab.org/cms/uploads/supplementary/Trimmomatic">http://www.usadellab.org/cms/uploads/supplementary/Trimmomatic</ext-link></p>
</fn>
<fn id="footnote2">
<label>2</label>
<p><ext-link ext-link-type="uri" xlink:href="http://cs.wellesley.edu/btjaden/Rockhopper/">http://cs.wellesley.edu/&#x223C;btjaden/Rockhopper/</ext-link></p></fn>
<fn id="footnote3">
<label>3</label>
<p><ext-link ext-link-type="uri" xlink:href="https://bioconductor.org/packages/release/bioc/html/edgeR.html">https://bioconductor.org/packages/release/bioc/html/edgeR.html</ext-link></p></fn>
<fn id="footnote4">
<label>4</label>
<p><ext-link ext-link-type="uri" xlink:href="https://github.com/tanghaibao/Goatools">https://github.com/tanghaibao/Goatools</ext-link></p></fn>
<fn id="footnote5">
<label>5</label>
<p><ext-link ext-link-type="uri" xlink:href="https://www.ebi.ac.uk/Tools/msa/clustalo/">https://www.ebi.ac.uk/Tools/msa/clustalo/</ext-link></p></fn>
<fn id="footnote6">
<label>6</label>
<p><ext-link ext-link-type="uri" xlink:href="https://services.healthtech.dtu.dk/service.php?TMHMM-2.0">https://services.healthtech.dtu.dk/service.php?TMHMM-2.0</ext-link></p></fn>
<fn id="footnote7">
<label>7</label>
<p><ext-link ext-link-type="uri" xlink:href="https://zhanglab.ccmb.med.umich.edu/I-TASSER/">https://zhanglab.ccmb.med.umich.edu/I-TASSER/</ext-link></p></fn>
</fn-group>
</back>
</article>