# Large-Scale Analysis of the Mycoplasma bovis Genome Identified Non-essential, Adhesion- and Virulence-Related Genes

Christoph Josi, Sibylle Bürki, Sara Vidal, Emilie Dordet-Frisoni, Christine Citti, Laurent Falquet, Paola Pilo

► To cite this version:

Christoph Josi, Sibylle Bürki, Sara Vidal, Emilie Dordet-Frisoni, Christine Citti, et al.. Large-Scale Analysis of the Mycoplasma bovis Genome Identified Non-essential, Adhesion- and Virulence-Related Genes. Frontiers in Microbiology, 2019, 10, pp.2085. 10.3389/fmicb.2019.02085. hal-02933215

HAL Id: hal-02933215

https://hal.inrae.fr/hal-02933215

Submitted on 20 Apr 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Large-Scale Analysis of the *Mycoplasma bovis* Genome Identified Non-essential, Adhesion- and Virulence-Related Genes

Christoph Josi<fnref n="1" />,², Sibylle Bürki<fnref n="1" />, Sara Vidal<fnref n="1" />, Emilie Dordet-Frisoni³, Christine Citti³, Laurent Falquet⁴† and Paola Pilo<fnref n="1" />*†

## OPEN ACCESS

**Edited by:**

Feng Gao,
Tianjin University, China

**Reviewed by:**

Yong-Qiang He,
Guangxi University, China
Angelika Lehner,
University of Zurich, Switzerland

*Correspondence:

Paola Pilo
pilop508@gmail.com

† These authors have contributed
equally to this work

**Specialty section:**

This article was submitted to
Evolutionary and Genomic
Microbiology,
a section of the journal
Frontiers in Microbiology

Received: 04 July 2019

Accepted: 23 August 2019

Published: 13 September 2019

**Citation:**

Josi C, Bürki S, Vidal S, Dordet-Frisoni E, Citti C, Falquet L and Pilo P (2019) Large-Scale Analysis of the Mycoplasma bovis Genome Identified Non-essential, Adhesion- and Virulence-Related Genes. Front. Microbiol. 10:2085. doi: 10.3389/fmicb.2019.02085

Mycoplasma bovis is an important pathogen of cattle causing bovine mycoplasmosis. Clinical manifestations are numerous, but pneumonia, mastitis, and arthritis cases are mainly reported. Currently, no efficient vaccine is available and antibiotic treatments are not always satisfactory. The design of new, efficient prophylactic and therapeutic approaches requires a better understanding of the molecular mechanisms responsible for M. bovis pathogenicity. Random transposon mutagenesis has been widely used in Mycoplasma species to identify potential gene functions. Such an approach can also be used to screen genomes and search for essential and non-essential genes for growth. Here, we generated a random transposon mutant library of M. bovis strain JF4278 containing approximately 4000 independent insertion sites. We then coupled high-throughput screening of this mutant library to transposon sequencing and bioinformatic analysis to identify M. bovis non-essential, adhesion- and virulence-related genes. Three hundred and fifty-two genes of M. bovis were assigned as essential for growth in rich medium. Among the remaining non-essential genes, putative virulence-related factors were subsequently identified. The complete mutant library was screened for adhesion using primary bovine mammary gland epithelial cells. Data from this assay resulted in a list of conditional-essential genes with putative adhesion-related functions by identifying non-essential genes for growth that are essential for host cell-adhesion. By individually assessing the adhesion capacity of six selected mutants, two previously unknown factors and the adhesin TrmFO were associated with a reduced adhesion phenotype. Overall, our study (i) uncovers new, putative virulence-related genes; (ii) offers a list of putative adhesion-related factors; and (iii) provides valuable information for vaccine design and for exploring M. bovis biology, pathogenesis, and host-interaction.

**Keywords:** Mycoplasma bovis, random transposon mutagenesis, non-essential genes, adhesion, virulence

<footnote n="1">¹ Department of Infectious Diseases and Pathobiology, Vetsuisse Faculty, Institute of Veterinary Bacteriology, University of Bern, Bern, Switzerland, ² Graduate School for Cellular and Biomedical Sciences, University of Bern, Bern, Switzerland, ³ UMR 1225, IHAP, Université de Toulouse, INRA, ENVT, Toulouse, France, ⁴ Department of Biology, Faculty of Science and Medicine, Swiss Institute of Bioinformatics, University of Fribourg, Fribourg, Switzerland</footnote>
# INTRODUCTION

*Mycoplasma bovis* is one of the major causative agents of bovine mycoplasmosis. This bacterium is emerging in industrialized countries, where it leads to high economic losses in the dairy and beef industries (Nicholas, 2011; Bürki et al., 2015). Currently, no effective vaccine is available and treatment with antibiotics is hardly efficient (Nicholas, 2011; Gautier-Bouchardon et al., 2014; Perez-Casal et al., 2017). Thus, the development of therapeutic and prophylactic measures to combat bovine mycoplasmosis is an ongoing challenge and requires in-depth understanding of the molecular mechanisms of pathogenicity of this bacterium (Calcutt et al., 2018).

Over the past years, some virulence factors of *M. bovis* have been described (Bürki et al., 2015). The most investigated ones are the variable surface lipoproteins (Vsps) that were shown to be involved in adhesion, antigenic variation, and immune evasion (Sachse et al., 1993, 1996, 2000; Rosengarten et al., 1994; Thomas et al., 2003; Buchenau et al., 2010). Four additional adhesion factors were characterized: the α-enolase (Song et al., 2012), the VpmaX protein (Zou et al., 2013), a NADH oxidase (NOX), which is also involved in H₂O₂ production (Zhao et al., 2017) and the TrmFO protein (Guo et al., 2017). Furthermore, two nucleases, the homologs MBOVPG45_0215 and MnuA, and the secretory nuclease MBOV_RS02825, were shown to be associated with cytotoxicity or with the degradation of neutrophil extracellular traps (Sharma et al., 2015; Zhang et al., 2016; Mitiku et al., 2018).

For a long time, the functional characterization of genes in mycoplasmas was hampered by the lack of genetic tools. Since mycoplasmas lack substantial parts of DNA recombination and repair mechanisms (Rocha et al., 2005), the frequency of successful targeted mutagenesis through homologous recombination is very low (Renaudin et al., 2014). For this reason, the fast generation of mutants by the use of transposons is widely used in the field of mycoplasmas (Hedreyda et al., 1993; Hutchison et al., 1999; Dybvig et al., 2000). Random transposon mutagenesis is now available in *M. bovis* (Chopra-Dewasthaly et al., 2005; Sharma et al., 2014). This approach has already led to the identification of the function of some genes in *M. bovis* (Sharma et al., 2015; Rasheed et al., 2017; Zhao et al., 2017). In the closely related species *M. agalactiae*, this technique was used to investigate growth phenotypes of individual or pools of mutants (Baranowski et al., 2010; Hegde et al., 2016).

Transposon mutagenesis has also been useful in identifying a number of mycoplasma genes essential for growth in culture media. In previous studies using other *Mycoplasma* spp. transposon mutant libraries with ≥1000 insertions/genome were built to identify essential genes for growth in rich medium (Glass et al., 2006; French et al., 2008; Dybvig et al., 2010; Lluch-Senar et al., 2015). Moreover, dispensable genes of *M. bovis* were previously identified by direct sequencing of several mutants (Sharma et al., 2014). Yet, a universal and consensual definition of essentiality versus non-essentiality is lacking (Hutchison et al., 2016). Approaches using random transposon mutagenesis take into account individual genes (Hutchison et al., 2016). Non-essentiality and essentiality assignment of genes with transposon

mutagenesis is based on the observation of whether an individual gene tolerates disruptions by a transposon (Glass et al., 2006; French et al., 2008; Dybvig et al., 2010; Sharma et al., 2014; Lluch-Senar et al., 2015). Thus, mutagenesis of single genes does not allow us to predict genes that could be simultaneously dispensable or mutants that are mutually exclusive (Hutchison et al., 2016). Moreover, the definition of essentiality versus non-essentiality is dependent on several factors including the selective environment, bacterial growth conditions, and makes use of an artificial environment potentially neglecting differences present in *in vivo* situations (Acevedo-Rocha et al., 2013; Tseng et al., 2013; Hegde et al., 2016; Hutchison et al., 2016).

In this study, we performed a large-scale analysis of the genome of *M. bovis* strain JF4278 to identify non-essential, adhesion- and virulence-related genes in rich culture medium and cell culture conditions. To achieve this aim, we coupled a high-throughput adhesion screening of a random transposon mutant library of *M. bovis* to bioinformatic analyses. Selected individual mutants were then tested for adhesion to confirm previous screening results.

# MATERIALS AND METHODS

## Bacterial Strain and Growth Conditions

*Mycoplasma bovis* strain JF4278 was isolated from the milk of a cow with pneumonia and mastitis in Switzerland in 2008 (Aebi et al., 2012). JF4278 was grown at 37°C in SP4 medium (Freundt, 1983) supplemented with 50 µg/mL cefoxitin sodium salt (Sigma-Aldrich, Buchs, Switzerland) for 24 h in broth medium or for 3-7 days on agar plates. For the growth of the transposon mutants, 15 µg/mL tetracycline hydrochloride (Sigma-Aldrich, Buchs, Switzerland) was added to the SP4 medium supplemented with cefoxitin. Mycoplasma cell concentrations were defined by CFU counting on agar plates using a stereomicroscope.

## Generation of a Random Transposon Mutant Library

A mutant library of JF4278 was generated by random transposon mutagenesis using a modified version of transposon Tn4001 (mini-Tn) inserted in the plasmid pMT85-Tet (Chopra-Dewasthaly et al., 2005; Zimmerman and Herrmann, 2005; Dordet Frisoni et al., 2013). *M. bovis* strain JF4278 was transformed using high concentrations of polyethylene glycol 8000 (PEG 8000; Sigma-Aldrich, Buchs, Switzerland) based on a protocol published for the transformation of *Ureaplasma parvum* (Aboklaish et al., 2014). Transformed *M. bovis* were plated on SP4 agar plates supplemented with cefoxitin and tetracycline. After 5-6 days single colonies were picked and suspended in 150 µL SP4 broth with cefoxitin and tetracycline in 96-well plates (TPP®, Trasadingen, Switzerland). The cultures were incubated for 24-48 h at 37°C. A total of approximately 3500-4000 individual mutants were collected and individually stored in 96-well plates. Pools of mutants were made from each plate by mixing 10 µL of individually grown mutants. These pools were combined in a single tube to collect all mutants, which were subsequently split into 1 mL aliquots and stored at -80°C.
## Primary Bovine Mammary Gland Epithelial Cells

Primary bovine mammary gland epithelial cells (bMec) were isolated from mammary gland tissues of cows directly after slaughter and cultured as previously described (Wellnitz and Kerr, 2004; Zbinden et al., 2015; Josi et al., 2018). No ethics approval was needed because primary cells were collected from organs of bovine carcasses at the slaughterhouse in accordance with the Swiss Federal Animal Protection Law, RS455. For experiments, cells were routinely seeded in T75 cell culture flasks or 24-well plates (TPP®, Trasadingen, Switzerland) and grown until confluency. Twenty hours before infection with *M. bovis*, the medium was changed to minimal essential medium (MEM)-Earle medium without addition of antibiotics but supplemented with 2.2 g/L NaHCO₃ (Biochrom, Berlin, Germany), 7% fetal bovine serum, and 1% L-glutamine (Biochrom, Berlin, Germany). At the time of infection, cell density reached approximately 1 × 10⁷ and 2.25 × 10⁵ cells per T75 flask or 24-well, respectively. For infection experiments, cell passages 3–5 were used.

## Screening of the Mutant Library Using an Adhesion Assay and DNA Extraction

The screening methodology of this study is summarized in **Figure 1**. A pool of the mutants, containing all collected mutants, was either directly lysed (input sample) or used for an adhesion assay (**Figure 1**). The full mutant library is referred as the input sample and the mutants recovered after the adhesion experiment are referred as the output sample (**Figure 1**). The input sample and output sample from one experiment were collected in parallel. The screening of the full random transposon mutant library was performed in three replicates leading to the collection of three input samples and three output samples. The DNA of input samples and output samples was sequenced by transposon sequencing (**Figure 1**). The sequencing data were processed to identify non-essential, adhesion- and virulence-related genes (**Figure 1**).

The adhesion assay was adapted in our laboratory from a previous protocol (Sachse et al., 1996; Sachse, 1998; Josi et al., 2018). For one replicate, one aliquot of the pooled complete mutant library was thawed, supplemented with 1 mL fresh SP4 medium, and incubated for 2 h at 37°C. Mycoplasma cultures were washed once and suspended in buffer A (0.05 M Tris-HCl, pH 7.2, 0.1 M NaCl, and 1 mM CaCl₂). Before infection, bMec were washed with buffer A and 6 mL buffer A was added in the T75 flask. Approximately 5 × 10⁷ mycoplasmas were added to infect bMec. *M. bovis* was allowed to adhere to bovine cells for 30 min on a rocker with 22 strokes/min and an amplitude of 3.5 cm at 37°C. bMec were washed three times with 15 mL PBS (Thermo Fisher Scientific) and manually rocked 30 times to remove unattached mycoplasmas. bMec with attached mycoplasmas were trypsinized and collected in 10 mL MEM-Earle medium. Twenty microliters of the output sample were used to determine the CFU. The remaining sample was centrifuged at 10,000 × g. Total DNA of the output sample was isolated using the DNeasy® Blood & Tissue Kit (QIAGEN GmbH, Hilden, Germany) for cultured cells according

to the manufacturer's protocol. To prevent saturation, three columns were used to extract DNA from each experiment. An initial RNase A treatment step was included and DNA from each column was eluted in 150 µL H₂O and pooled. In parallel, mycoplasmal concentration of the input sample was also determined by CFU counting. The input sample was then centrifuged and lysed with lysis solution [buffer A supplemented with 1% Tween®20 (Sigma-Aldrich, Buchs, Switzerland) and 0.24 mg/mL Proteinase K (Roche Diagnostics, Rotkreuz, Switzerland)]. The input sample was then heated for 1 h at 65°C and 15 min at 96°C. The percentage of *M. bovis* adherent to bMec was calculated by dividing the CFU calculated in the input sample by the CFU counted in the output sample. The assay was performed in triplicates leading to the collection of three input DNA samples and three output DNA samples. The concentration and purity of the three DNA samples from inputs and outputs were checked with a NanoDrop 1000 Spectrophotometer (Thermo Fisher Scientific) and stored at -20°C for further use.

## DNA Library Preparation for Transposon Sequencing on the Illumina NextSeq Platform

Input and output DNA samples were prepared for transposon sequencing on the Illumina NextSeq Platform. Two major steps were performed: (1) Enrichment of amplicons from the transposon-junctions and depletion of eukaryotic DNA and (2) Amplification of the library to add Illumina adapters. The method is briefly described here but the full and detailed protocol is in **Supplementary Table 1**. A method for amplicon library preparation was set up involving biotinylated primers specific for the inserted transposon present in *M. bovis* mutants. Purifications were done using streptavidin-coated magnetic beads to limit contamination with eukaryotic DNA in output samples. Amplicons of the transposon-junction sequences were prepared using a non-restrictive linear amplification-mediated PCR (nrLAM-PCR) approach (Gabriel et al., 2009). A previously published protocol was followed (Paruzynski et al., 2010) with the following modifications. Instead of a linker-ligation as described in the study of Paruzynski et al. (2010), a homopolymer-C tail of controlled length was added to linear PCR products needed for subsequent exponential PCRs (Lazinski and Camilli, 2013). The combination of nrLAM-PCR with homopolymer tail-mediated ligation PCR (HTML-PCR) was previously published (Dawoud et al., 2014). Amplicon libraries were sent to Microsynth (Microsynth AG, Balgach, Switzerland) for barcoding and adding Illumina adapters. The six libraries were gel purified, quantified, pooled, and sequenced in the same run with Illumina NextSeq 1 × 75 bp high-output.

## Bioinformatics Analysis
### Illumina Sequence Data Analysis

Transposon sequencing reads of input and output samples were obtained in FASTQ format. Adapter, primer, and transposon sequences were removed and reads were trimmed to 25 bp using Trimmomatic (Bolger et al., 2014). Raw reads and trimmed
<figure><img src="image_3.png" /></figure>

**FIGURE 1** | Schematic representation of the adhesion screening and bioinformatics analyses. The adhesion screening was performed with the full mutant library containing approximately 4000 mutants (input sample). After adhesion, adherent mutants (output sample) were recovered. The screening of the same full random transposon mutant library was performed in three replicates, leading to the collection of three input samples and three output samples. Transposon sequencing was carried out by parallel sequencing of the regions flanking each transposon to identify disrupted genes. Transposons sequencing read data from input samples from triplicates were used to categorize genes as essential and non-essential for growth on rich medium using the ESSENTIALS software (Zomer et al., 2012). Transposons sequencing read data of input and output samples from triplicates were then compared to identify conditional-essential genes (essential for adhesion) using the ESSENTIALS software (Zomer et al., 2012).

reads from input and output samples were deposited in the European Nucleotide Archive (ENA) (Bioproject: PRJEB31021). Quality control of trimmed reads was performed with FastQC and reads were mapped against the *M. bovis* strain JF4278 reference genome (GenBank accession no. NZ_LT578453.1) with Bowtie 2 (Langmead and Salzberg, 2012). Alignment parameters are shown in **Supplementary Table 1**. BAM files were generated after mapping and filtered for mapping quality ≥12 using BamTools (Barnett et al., 2011). These bioinformatics procedures were performed using the web-based scientific analysis platform Galaxy<fnref n="1" /> (Afgan et al., 2018). Filtered BAM files were used for visualization of the alignments in Geneious version R11.1.2 (Kearse et al., 2012).

## Essentiality Analysis

Essentiality analysis was performed with ESSENTIALS<fnref n="2" />, which is an open source, web-based software tool for insertion sequencing analysis (Zomer et al., 2012). Trimmed reads from input samples were uploaded as individual control samples from the same library. Reads were aligned directly to the *M. bovis* strain JF4278 reference genome with ESSENTIALS. Count data for read counts per insertion or gene were generated and

compared to the expected number of reads per gene calculated by ESSENTIALS. Count data were combined and normalized for differential expression analysis with EdgeR by ESSENTIALS (Robinson et al., 2010; Zomer et al., 2012). Parameters used for the essentiality analysis by ESSENTIALS are shown in **Supplementary Table 1**.

## Conditional-Essentiality Analysis

Conditional-essentiality analysis which corresponds to the detection of genes identified as non-essential for growth but essential for adhesion was performed using ESSENTIALS (Zomer et al., 2012). Compared to the essentiality analysis, output samples were included and uploaded as individual target samples from the same library. Input samples were uploaded as individual control samples from the same mutant library. Parameters used for the ESSENTIALS run were kept the same as explained above. In addition to the prediction of adhesion-related factors by ESSENTIALS, detection of putative adhesins with the SPAAN software (Sachdeva et al., 2005) was performed on the predicted *M. bovis* proteome.

## KEGG Pathway Analysis

Functional KEGG pathway annotation of JF4278-predicted proteins was achieved using the web-based version of KOBAS 3.0 against the annotated genome of *M. bovis* type strain

<footnote n="1">¹https://usegalaxy.org</footnote>

<footnote n="2">²http://bamics2.cmbi.ru.nl/websoftware/essentials/essentials_start.php</footnote>
PG45³ (GenBank accession no. NC_014760.1) (Wise et al., 2011; Xie et al., 2011). Since KEGG pathway annotation with KOBAS 3.0 was performed for protein sequences only, rRNAs and tRNAs were annotated manually to KEGG pathways via the *M. bovis* PG45 reference genome available on the KEGG database<fnref n="4" />. KEGG pathway enrichment of essential genes compared to the complete gene set of JF4278 was assessed. First, KEGG pathway fold enrichment (fold enrichment of proportion of JF4278 essential genes compared to the proportion of all JF4278 genes belonging to the respective KEGG pathway) was calculated. Secondly, the significance of overrepresentation of the respective KEGG pathway in the essential geneset compared to the complete geneset of JF4278 was calculated using Fisher's exact test.

## Comparison Between Essential Genes of JF4278 and Non-essential Genes in *M. bovis* PG45

*Mycoplasma bovis* PG45 non-essential protein-coding genes were blasted in NCBI protein-protein BLAST<fnref n="5" /> to find the respective JF4278 homologs (Sharma et al., 2014). Non-essential genes of PG45 that were found essential in JF4278 were compared to other *Mycoplasma* spp. (French et al., 2008; Dybvig et al., 2010; Lin and Zhang, 2011; Liu et al., 2012). The persistence of these genes in *Mycoplasma* spp. was assessed as described in the study of Liu et al. (2012). Since the latter study did not include *M. bovis* strains *M. agalactiae* strain PG2 (GenBank accession no. NC_009497.1) orthologs to JF4278 were first identified by protein-protein BLAST (Sirand-Pugnet et al., 2007b). Using PG2 orthologs, orthologs to other *Mycoplasma* spp. were then identified based on the study by Liu et al. (2012) and compared to the essential gene set of *M. pulmonis* strain UAB CTIP (French et al., 2008; Dybvig et al., 2010). *M. pulmonis* was selected for comparisons due to its thorough assignment of essential genes and because it belongs to the Hominis group like *M. bovis* (May et al., 2014). Finally, *M. pulmonis* orthologs were analyzed to define if they belong to the putative core set of essential genes in *Mycoplasma* genomes<fnref n="6" /> (Lin and Zhang, 2011).

## Additional Softwares

To show density distributions of essentiality values obtained with ESSENTIALS, Kernel density plots using an Epanechikov model were generated by a web-based tool<fnref n="7" />. The software TMPred<fnref n="8" /> was used to predict putative transmembrane helices (TM-Hs). Statistical analysis, dot plots, and box plots were made in R version 3.4.3.

## Isolation of Candidate Mutants Involved in Adhesion

Selected mutants identified as conditional-essential in the adhesion assays were recovered from the mutant library by PCR screening using a step-wise approach. PCR primers were

designed using the alignments obtained from the transposon sequencing data (Supplementary Table 1). PCRs were performed on the mutant pool of all collected 96-well plates (36 pools), to identify the 96-well plate containing each mutant. Thereafter, PCRs were carried out on pools of each row of the respective plate and of each well within the positive row. The culture obtained from the well assigned to the selected mutant was filtered with a 0.2-µm syringe filter and plated on SP4 agar plates containing tetracycline. A single colony was picked and grown in SP4 liquid medium containing tetracycline. A final PCR and sequencing of the obtained PCR product were performed to confirm the isolation of the selected mutants.

## Adhesion Assay With Selected Mutants

The adhesion assay for the relative quantification of adherent mycoplasmas was adapted from a previous protocol established in our laboratory (Josi et al., 2018). The quantification of adherent *M. bovis* was carried out using real-time qPCR (Rossetti et al., 2010). The adhesion assay was performed for each selected mutant individually as described in the section "Screening of the Mutant Library Using an Adhesion Assay and DNA Extraction" but in a 24-well plate setup. As described previously (Josi et al., 2018), lysates of the *M. bovis*-mutant inoculum and infected cell lysates were collected to quantify the relative amount of adherent mycoplasmas. The percentage of adherent *M. bovis* relative to the amount of *M. bovis* added for the infection was calculated as follows (adapted from Nicholson et al., 2012):

$$ \% \text{ of adherent } M. bovis = 2^{(a-b)} * 100 $$

where a = C<sub>t</sub> of mycoplasmas added for infection; and b = C<sub>t</sub> of adherent mycoplasmas. The assays were performed in duplicates in four independent experiments. Statistical comparison of each mutant to the wild type (wt) was performed by the two-sample Wilcoxon test.

# RESULTS AND DISCUSSION

## Mutant Library Transposon Sequencing Results and Alignment to *M. bovis* Strain JF4278 Reference Genome

The transposon sequencing results from the three replicates of input and output samples are shown in Table 1. Transposon sequencing runs yielded between 39 and 56 million reads for the different input and output samples (Table 1). The transposon sequencing reads were aligned to the *M. bovis* strain JF4278 genome retrieved from NCBI (GenBank accession no. NZ_LT578453.1). The alignment rate was between 43 and 67%, corresponding to between 22 and 28 million aligned reads for the different input and output samples (Table 1). The GC content of unaligned reads was 67.5%. In these reads, the G content increased steadily to nearly 100% toward the end of the read, reflecting sequencing of the added poly-C tail during amplicon library preparation. The remaining specific sequences at the start of unaligned reads were too short to be accurately aligned. On average 26 million transposon sequencing reads

<footnote n="3">³http://kobas.cbi.pku.edu.cn/index.php</footnote>

<footnote n="4">⁴https://www.genome.jp/kegg/</footnote>

<footnote n="5">⁵https://blast.ncbi.nlm.nih.gov/Blast.cgi</footnote>

<footnote n="6">⁶https://tubic.org/pdeg/core.php?page=1</footnote>

<footnote n="7">⁷https://www.wessa.net/rwasp_density.wasp</footnote>

<footnote n="8">⁸https://embnet.vital-it.ch/software/TMPRED_form.html</footnote>
<table><caption>**TABLE 1** | Summary of the transposon sequencing results from input and output samples of the adhesion assay.</caption><thead><tr><th>Sample name</th><th>Adhesion assay replicate<fnref n="1" /></th><th>Total number of transposon sequencing reads (million)</th><th>Number of reads aligned to JF4278 (million)</th><th>Alignment rate to JF4278 (%)</th></tr></thead><tbody><tr><td>Input 1</td><td>1</td><td>39</td><td>26</td><td>67</td></tr><tr><td>Output 1</td><td>1</td><td>46</td><td>26</td><td>57</td></tr><tr><td>Input 2</td><td>2</td><td>44</td><td>28</td><td>64</td></tr><tr><td>Output 2</td><td>2</td><td>56</td><td>28</td><td>50</td></tr><tr><td>Input 3</td><td>3</td><td>47</td><td>28</td><td>60</td></tr><tr><td>Output 3</td><td>3</td><td>51</td><td>22</td><td>43</td></tr><tr><td>Average</td><td>N/A</td><td>47</td><td>26</td><td>55</td></tr></tbody></table>

per sample could be aligned to the JF4278 reference genome. Therefore, with approximately 4000 independent transposon insertions sites within the JF4278 genome, the average sequencing depth (reads/transposon insertion) of each transposon-junction is approximately 6500 per sample. The size of the *M. bovis* strain JF4278 genome is 1,038,531 bp. With approximately 4000 independent transposon insertions sites within the JF4278 genome, the average density of our mutant library is one insertion event every 260 bp.

A schematic illustration of the alignment of the transposon sequencing reads from the three replicates of the full mutant library (input samples) is shown in **Figure 2**. Distribution of transposon insertion sites as well as read counts per transposon insertion (height of bars) among the three replicates of the full mutant library (input samples) is highly similar (Figure 2). The sequencing depth between different transposon insertion sites and the insertion density between different genes varies and this data were used to identify essential and non-essential genes of *M. bovis* strain JF4278 with the software ESSENTIALS (Figure 2). No obvious hot and cold spots for transposon *Tn4001* integration were observed (Figure 2). Compared to integration with *Tn916*, also commonly used in Mollicutes, *Tn4001* is not known to integrate preferably in hot spots in mycoplasmas (Renaudin et al., 2014). However, a preference of secondary insertions of *Tn4001* in regions rich in palindromic and cruciform elements was shown in *M. genitalium* (Glass et al., 2006).

## Essential Genes of *M. bovis* Strain JF4278 Can Be Identified With the Software ESSENTIALS

Genes essential for growth in SP4 medium were analyzed using the software ESSENTIALS (Zomer et al., 2012). Non-essentiality and essentiality assignment of genes is not based solely on the absence or presence of gene disruption but rather on the amount of disruption tolerated by a gene. Read data from input samples were used for the essentiality analysis. An expected number of reads per gene is automatically calculated by ESSENTIALS. This is based on the size of the used mutant library, the number of insertion sites per gene, and the sequencing depth of the transposon-junctions (Zomer et al., 2012). The obtained

read counts per gene of input samples were compared to the expected number of reads per gene calculated by ESSENTIALS (Figure 3). Thereby, genes are assigned a fold change (FC) value of read counts in log2 scale. Log2(FC) values for all JF4278 genes were obtained with ESSENTIALS (Figure 3 and Supplementary Table 2). The log2(FC) values of JF4278 genes show a bimodal distribution in the density plot (Figure 3). The log2(FC) value of non-essential genes is around zero, as the obtained read counts correspond to the expected number of reads per gene (Figure 3). On the other hand, genes with underrepresented read counts [strongly negative log2(FC) values] are more likely to be essential (Figure 3). The log2(FC) value of a gene can be interpreted similarly to a fitness score where significantly underrepresented (adj. P-value ≤ 0.05) genes below an automatically calculated log2(FC) cut-off value are considered essential. The log2(FC) cut-off suggested by ESSENTIALS for our dataset was <−4.64 (Figure 3).

According to these criteria [log2(FC) < −4.64 and adj. P-value ≤ 0.05], 352 out of 900 genes were found to be essential for growth in SP4 medium (Supplementary Table 2). For 36 out of the remaining 548 non-essential genes, log2(FC) and adjusted P-values could not be obtained. These genes mainly encode for proteins with DDE transposase and/or integrase core domain and were present in multiple copies (Supplementary Table 2). Reads from highly repetitive regions or from genes present in multiple copies cannot be assigned to a specific position in the genome and are excluded by the software ESSENTIALS (Zomer et al., 2012). However, if an insertion occurs in a non-repetitive segment of gene containing highly repetitive sequences, this read will be taken into account by the software ESSENTIALS. For this reason, not all DDE transposases are among the 36 genes without log2(FC) and adjusted P-values (Supplementary Table 2).

Distribution of transposon insertion sites, as well as the read counts in log10 scale (height of bars) from the three sequencing replicates of the full mutant library (input samples), is in line with the predicted essentiality of genes (Figure 2 and Supplementary Table 2). Genes with no or few transposon insertions and with limited numbers of reads were identified as essential (Figure 2, red arrows), while non-essential genes (Figure 2, yellow arrows) contained more insertions and higher read counts/gene. The high density of essential genes found in the locus between positions nt 679,886 and 729,966 corresponded to genes encoding small ribosomal proteins, ribosomal RNA, and other factors involved in translation (Figure 2 and Supplementary Table 2). Other prominent essential operons contained genes encoding proteins involved in oligopeptide permease (Opp) systems (nt 91,927–100,788) and subunits of the $F_1F_0$ ATPase (nt 484,670–477,651) (Beven et al., 2012).

## Essential Genes Are Involved in Genetic Information Processing and Metabolic Pathways

Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analysis of essential genes of strain JF4278 was performed. Using PG45 as reference genome, 277 genes of strain JF4278 were assigned to one or several of the 52 KEGG pathways

<footnote n="1">¹Input and output samples from one replicate were collected in parallel.</footnote>
<figure><img src="image_4.png" /></figure>

**FIGURE 2** | Location of transposon insertion sites in the genome of *M. bovis* strain JF4278. Genes and their orientation are indicated by arrows. Genes found essential with ESSENTIALS are filled in red, while non-essential genes are filled in yellow (Zomer et al., 2012). Mapping of the transposon integrations sites from the transposon sequencing of the full mutant library (input samples) is represented by black, green, and red bars (each corresponding to one replicate experiment). The height of the insertion bars represents the read count per insertion. Approximately 4000 independent transposon insertion sites were detected with an average sequencing depth (reads/transposon insertion) of 6800 per transposon-junction per replicate.

(Supplementary Table 3). Two hundred and nine genes out of the 277 were found to be essential with the ESSENTIALS software. First, KEGG pathway fold enrichment (fold enrichment of proportion of JF4278 essential genes compared to the proportion of all JF4278 genes belonging to the respective KEGG pathway) was calculated (Supplementary Table 3). Secondly, the significance of the overrepresentation of the respective KEGG pathway in the set of essential JF4278 genes compared to the complete set of JF4278 genes was calculated using Fisher's exact test (Supplementary Table 3). Eighteen out of the 52 KEGG pathways were significantly overrepresented (P-value ≤ 0.05) in the set of essential genes of strain JF4278 (Supplementary

**Table 3**). KEGG pathways with a highly significant enrichment (P-value ≤ 0.01) were examined in more detail (Figure 4). Genetic information processing pathways were among the top hits. Nearly all genes assigned to pathways “Ribosome,” “Aminoacyl-tRNA biosynthesis,” and “DNA replication” were found to be essential in JF4278. Only three genes were assigned to transcription in the pathway “RNA polymerase.” Even though these three genes were found to be essential, no significant enrichment was obtained (Supplementary Table 3). Genes involved in replication and translation belong to the basic cellular machinery and were also found to be essential in many bacterial species (Grazziotin et al., 2015). Furthermore, many
genes assigned to metabolic pathways were found to be essential in JF4278 (Figure 2). Among the top hits were genes involved in nucleotide metabolism (purine and pyrimidine metabolism), biosynthesis of secondary metabolites and amino acids. Genes assigned to oxidative phosphorylation encode for subunits of the F₁F₀ ATPase found in different *Mycoplasma* spp. However, in *Mycoplasma* spp. the F₁F₀ ATPase is involved in the maintenance of the electrochemical gradient by ATP hydrolysis rather than in the generation of ATP (Beven et al., 2012). While carbon metabolism genes were significantly overrepresented in the set of essential JF4278 genes, individual carbohydrate metabolism pathways were not (Figure 4 and Supplementary Table 3). In *Mycoplasma* spp., glycerol derived from host-phospholipids is a major source of carbon and energy (Blötz and Stülke, 2017). In the *M. bovis* strain JF4278, most of the genes assigned to “glycerophospholipid metabolism” were found to be essential (Supplementary Table 3). In *Mycoplasma* spp. many metabolic pathways are incomplete and mycoplasmas strongly depend on host carbon sources and nutrients (Razin et al., 1998; Masukagami et al., 2017). It can be expected that the remaining genes involved in metabolic pathways of mycoplasmas are essential to further metabolize nutrients. In *M. bovis* strain JF4278 the remaining parts of metabolic pathways were indeed found to be essential for growth on rich medium. Furthermore, the high proportion of essential ABC transporters found in the *M. bovis* strain JF4278 genome might reflect the strong dependency of mycoplasmas on nutrients uptake (Figure 4 and Supplementary Table 3).

## Nineteen Non-essential Genes of PG45 Are Essential in Strain JF4278

Protein-coding genes identified as non-essential in PG45 by Sharma et al. (2014) were compared to the results obtained for strain JF4278 with ESSENTIALS. In PG45, 319 independent transposon insertions were sequenced directly from genomic DNA. Based on the presence of at least one transposon insertion located after the first three codons and within the first 85% of the protein-coding sequence, 128 genes were assigned as non-essential in PG45 (Sharma et al., 2014). Using NCBI protein-protein BLAST, 118 out of the 128 non-essential PG45 were found to have homologs in JF4278 (Supplementary Table 4). Concerning the remaining 10 genes, no homolog was found in JF4278 (Supplementary Table 4).

Among these 118 genes, 99 were also found to be non-essential in JF4278 (Supplementary Table 4). Thereby, non-essentiality was confirmed for the vast majority of genes (99/118). However, 19 homologs were differentially assigned as essential in JF4278 compared to non-essential in PG45 (Supplementary Table 4). To assess this discrepancy, these 19 genes were compared to results previously obtained with *M. pulmonis* (Table 2). Ten out of these 19 genes had essential orthologs in *M. pulmonis* strain UAB CTIP (MYPU) (French et al., 2008; Dybvig et al., 2010; Table 2). Six out of the 19 differentially assigned genes were shown to be essential in other *Mycoplasma* spp. and even to belong to the core set of essential genes in *Mycoplasma* genomes (Table 2; Glass et al., 2006; French et al., 2008; Dybvig et al., 2010; Lin and Zhang, 2011; Sharma et al., 2014; Lluch-Senar et al., 2015).

<figure><img src="image_5.png" /><figcaption>**FIGURE 3** | Density plot of essentiality values obtained with the software ESSENTIALS. Read counts per gene of the input samples were compared to the expected number of reads per gene calculated by ESSENTIALS (Zomer et al., 2012). *M. bovis* strain JF4278 genes are assigned a fold change (FC) value of read counts in log2 scale compared to the expected number of reads per gene. Log2(FC) values (essentiality values) for all genes are shown on the x-axis of the density plot. Essentiality values show a bimodal distribution. The log2(FC) value of non-essential JF4278 genes is around zero. Genes more likely to be essential have negative log2(FC) values. The cut-off value for essential genes suggested by the software ESSENTIALS for this data set is log2(FC) < -4.64.</figcaption></figure>

Notably, among the non-essential genes identified in the study by Sharma et al. (2014), five genes were found to be essential in all previous gene essentiality studies in mycoplasmas. These five genes (*trmE*, *potB*, *metK*, *dnaJ*, and *smpB*) even belonged to the core set of essential genes in *Mycoplasma* genomes (Lin and Zhang, 2011). However, all of these five genes were found to be essential in *M. bovis* strain JF4278 (Table 2). In fact, *trmE* (MBOVJF4278_00029), *potB* (MBOVJF4278_00110), *metK* (MBOVJF4278_00655), *dnaJ* (MBOVJF4278_00841), and *smpB* (MBOVJF4278_00857) even belonged to the top candidates of JF4278 essential genes identified by ESSENTIALS [log2(FC) value ≤ -8.4, Supplementary Table 2]. Discrepancies in the essentiality assignment of *M. bovis* genes might be due to strain specificity as well as to culture conditions or to mutant library saturation or to the used essentiality analysis method. In consequence, the discrepancies in the essentiality assignment of genes belonging to the core set of essential genes in other *Mycoplasma* spp. should be considered with caution (Table 2). It should be also mentioned that some biological factors can influence the physical accessibility of transposases to some segments of the chromosome. For example, proteins binding DNA and the structure of the chromosomal DNA might modify the frequency of insertion and lead to the identification of false essential genes (Kimura et al., 2016).

## Non-essential Genes of *M. bovis* Strain JF4278

Genes identified as non-essential in JF4278 represent a list of individual genes tolerating disruption (Supplementary Table 2). The functional assignment of many *M. bovis* genes is often based on homology predictions and many genes are
<figure><img src="image_6.png" /><figcaption>**FIGURE 4** | KEGG pathway enrichment analysis of essential genes in strain JF4278. Different KEGG pathways from *M. bovis*, with the corresponding mbov ID number, are displayed. Significant (*P*-value ≤ 0.01) overrepresented KEGG pathways in the essential gene set compared to the complete gene set of strain JF4278 are shown. Gray columns correspond to the *P*-value for significance of enrichment in –log10 transformation. Spotted columns correspond to the fold enrichment of individual KEGG pathways between complete and essential gene set. Significance of enrichment was calculated with the Fisher's exact test.</figcaption></figure>

<table><caption>**TABLE 2** | Protein-coding genes found as non-essential in PG45 but essential in JF4278.</caption><thead><tr><th>Essential CDS in<br>JF4278<br>(MBOVJF4278_)</th><th>Product</th><th>Non-essential<br>PG45 homolog<fnref n="1" /><br>(MBOVPG45_)</th><th>PG2<br>ortholog²<br>(MAG_)</th><th>Gene<br>persistence³</th><th>Essential UAB<br>CTIP ortholog<fnref n="4" /><br>(MYPU_)</th><th>Core<br>mycoplasma<br>genes⁵ (CEMyc_)</th></tr></thead><tbody><tr><td>00029</td><td>tRNA modification GTPase MnmE</td><td>0060</td><td>0530</td><td>20</td><td>0130</td><td>0050</td></tr><tr><td>00034</td><td>50S ribosomal protein L34</td><td>0065</td><td>7500</td><td>20</td><td>1540</td><td>N/A</td></tr><tr><td>00088</td><td>Oligopeptide transport ATP-binding protein OppF</td><td>0116</td><td>1040</td><td>15</td><td>N/A</td><td>N/A</td></tr><tr><td>00110</td><td>Transport system permease protein PotB</td><td>0135</td><td>1260</td><td>20</td><td>4240</td><td>0750</td></tr><tr><td>00180</td><td>Hypothetical protein</td><td>0662</td><td>N/A</td><td>N/A</td><td>N/A</td><td>N/A</td></tr><tr><td>00300</td><td>SUA5-like translation suppressor</td><td>0534</td><td>4520</td><td>17</td><td>6130</td><td>N/A</td></tr><tr><td>00382</td><td>Potassium transporter trkA</td><td>0464</td><td>3750</td><td>17</td><td>1370</td><td>N/A</td></tr><tr><td>00448</td><td>Endonuclease MutS2</td><td>0404</td><td>3140</td><td>1</td><td>N/A</td><td>N/A</td></tr><tr><td>00463</td><td>Hypothetical protein</td><td>0390</td><td>3010</td><td>7</td><td>7700</td><td>N/A</td></tr><tr><td>00529</td><td>Hypothetical protein</td><td>0349</td><td>4660</td><td>11</td><td>N/A</td><td>N/A</td></tr><tr><td>00561</td><td>Putative endonuclease 4</td><td>0317</td><td>4990</td><td>20</td><td>6210*</td><td>01290</td></tr><tr><td>00655</td><td>S-adenosylmethionine synthase</td><td>0227</td><td>5800</td><td>20</td><td>7020</td><td>01380</td></tr><tr><td>00755</td><td>Modification methylase DpnIIIB</td><td>0766</td><td>6680</td><td>4</td><td>N/A</td><td>N/A</td></tr><tr><td>00792</td><td>Ribosome maturation factor RimP</td><td>0800</td><td>7000</td><td>7</td><td>N/A</td><td>N/A</td></tr><tr><td>00833</td><td>Putative RNA pseudouridine synthase</td><td>0831</td><td>7190</td><td>11</td><td>N/A</td><td>N/A</td></tr><tr><td>00841</td><td>Chaperone protein DnaJ</td><td>0839</td><td>7280</td><td>12</td><td>7330</td><td>01460</td></tr><tr><td>00846</td><td>Phenylalanine-tRNA ligase beta subunit</td><td>0845</td><td>7310</td><td>14</td><td>4860</td><td>N/A</td></tr><tr><td>00857</td><td>SsrA-binding protein</td><td>0855</td><td>7390</td><td>20</td><td>3520</td><td>0650</td></tr><tr><td>00865</td><td>Catabolite control protein A</td><td>0858</td><td>7400</td><td>5</td><td>N/A</td><td>N/A</td></tr></tbody></table>

annotated as hypothetical proteins with unknown function (Supplementary Table 2). Therefore, the list of non-essential JF4278 genes might be useful for future studies aiming to

generate *M. bovis* mutants of specific genes for subsequent functional characterization and broaden the knowledge of *M. bovis* biology.

<footnote n="1">¹Non-essential homolog in *M. bovis* PG45 (Sharma et al., 2014). ²*M. agalactiae* strain PG2 ortholog. ³Gene persistence out of 20 *Mycoplasma* spp. (Liu et al., 2012).</footnote>

<footnote n="4">⁴Essential gene in *M. pulmonis* strain UAB CTIP (French et al., 2008). ⁵Core set of essential genes in *Mycoplasma* genomes (Lin and Zhang, 2011). *Non-essential gene in *M. pulmonis* strain UAB CTIP (Dybvig et al., 2010).</footnote>
# Several Virulence-Related Factors Identified With the Virulence Factor Database Are Non-essential in *M. bovis* Strain JF4278

To investigate whether putative virulence-related factors are dispensable for *M. bovis*, all the putative proteins of *M. bovis* were analyzed using the Virulence Factor Database (VFDB) (Chen et al., 2016). The hypothetical proteins of *M. bovis* strain JF4278 were aligned against the full dataset of VFDB, containing characterized and predicted virulence-related factors. Seventy-eight genes in the genome of JF4278 were predicted to encode for virulence-related factors based on a BLAST score ≥80 (Supplementary Table 5). In Supplementary Table 5, the gene product names of the VFDB hit are given and not the product names of *M. bovis* strain JF4278. Additionally, previously characterized *M. bovis* virulence-related factors involved in adhesion, cytotoxicity, and H₂O₂ production were identified in JF4278 based on a BLAST score ≥80 (Wise et al., 2011; Song et al., 2012; Zou et al., 2013; Sharma et al., 2015; Zhang et al., 2016; Guo et al., 2017; Zhao et al., 2017; Mitiku et al., 2018; Supplementary Table 5). With the VFDB candidates and the known *M. bovis* virulence factors, 91 genes in the JF4278 genome were predicted to encode for virulence-related factors (Supplementary Table 5). As described above, a stricter non-essentiality assignment was applied to the dataset. Forty-six out of the 91 genes were found to have more than one insertion within the core (5–80%) of the gene (Supplementary Table 5). These genes represent putative virulence-related factors that are dispensable for *M. bovis* strain JF4278 growth in rich medium. However, many genes identified via BLAST against the VFDB data base show a very low BLAST score and *E*-value to the *M. bovis* gene and might not be homologs.

Top candidates of non-essential putative virulence-related factors are shown in Table 3 (BLAST score ≥ 200). Additionally, several previously characterized *M. bovis* virulence factors were found to be non-essential: the NADH oxidase (NOX; MBOVJF4278_00281) (Zhao et al., 2017); the TrmFO protein (MBOVJF4278_00132) (Guo et al., 2017), the VpmaX protein (MBOVJF4278_00429) (Zou et al., 2013), and the membrane nuclease MnuA protein (MBOVJF4278_00670) (Sharma et al., 2015; Mitiku et al., 2018; Table 3). Other known *M. bovis* virulence factors like the secretory nuclease (MBOVJF4278_00567; Zhang et al., 2016) or an α-enolase (MBOVJF4278_00443) (Song et al., 2012) were found to be essential or with only few insertions (Supplementary Table 5). The *vsp*-locus of *M. bovis* strain JF4278 comprises several genes (MBOVJF4278_00797–MBOVJF4278_00821) and was identified based on a protein BLAST against PG45 Vsps (Wise et al., 2011). All Vsps identified were found to be non-essential with ESSENTIALS (Supplementary Table 5). Since reads from highly repetitive regions of *vsp*-genes cannot be assigned to a specific position in the genome, they are not used by ESSENTIALS, resulting in *vsp*-genes showing few insertions (Supplementary Table 5). Whether just few insertions are tolerated or the complete *vsp*-locus of strain JF4278 is dispensable remains unknown. Moreover, due to the high frequency of phase

and size variation of Vsps, the expressed Vsp forms in the different mutants cannot be predicted (Behrens et al., 1994; Lysnyansky et al., 1996, 1999, 2001).

The exact function of putative virulence-related factors identified via VFDB in *M. bovis* remains speculative (Table 3). However, some notable factors were characterized in previous studies. The P48 lipoprotein (MBOVJF4278_00889 and _00150) stimulates the host humoral response in *M. agalactiae* (Table 3; Rosati et al., 2000). The endopeptidase Clp ATP-binding chain C (ClpC; MBOVJF4278_00700) is required for *Listeria monocytogenes* adhesion and invasion (Nair et al., 2000).

<table><caption>TABLE 3 | Non-essential virulence-related factors in *M. bovis* strain JF4278.</caption><thead><tr><th>Non-essential CDS in JF4278 (MBOVJF4278_)</th><th>VFDB hit and/or <i>M. bovis</i> virulence attributes</th><th>Blast score</th><th>Blast <i>E</i>-value</th></tr></thead><tbody><tr><td>00028</td><td>IS1634AV transposase [MmmSC]</td><td>1060</td><td>0</td></tr><tr><td>00094</td><td>IS1634AV transposase [MmmSC]</td><td>1059</td><td>0</td></tr><tr><td>00280</td><td>IS1634AV transposase [MmmSC]</td><td>1025</td><td>0</td></tr><tr><td>00861</td><td>IS1634AV transposase [MmmSC]</td><td>1018</td><td>0</td></tr><tr><td>00281</td><td>NADH oxidase [M. bovis]<sup>1</sup></td><td>926</td><td>0</td></tr><tr><td>00808</td><td>Variable surface lipoprotein N [M. bovis]<sup>2</sup>; (VpmaVprecursor [M. agalactiae])</td><td>870; (313)</td><td>0; (1E-85)</td></tr><tr><td>00132</td><td>TrmFO [M. bovis]<sup>3</sup></td><td>859</td><td>0</td></tr><tr><td>00670</td><td>Membrane nuclease MnuA [M. bovis]<sup>4</sup>; (membrane nuclease [M. pulmonis])</td><td>751; (173)</td><td>0; (2E-43)</td></tr><tr><td>00889</td><td>P48, predicted lipoprotein [M. agalactiae]</td><td>743</td><td>0</td></tr><tr><td>00150</td><td>P48, predicted lipoprotein [M. agalactiae]</td><td>627</td><td>1E-180</td></tr><tr><td>00700</td><td>Endopeptidase Clp ATP-binding chain C [Listeria monocytogenes]</td><td>546</td><td>1E-155</td></tr><tr><td>00006</td><td>Oligopeptide ABC transporter [MmmSC]</td><td>414</td><td>1E-116</td></tr><tr><td>00429</td><td>VpmaX [M. bovis]<sup>5</sup></td><td>353</td><td>2E-126</td></tr><tr><td>00819</td><td>Variable surface lipoprotein K [M. bovis]</td><td>289</td><td>2E-103</td></tr><tr><td>00390</td><td>PgPepO oligopeptidase [Mycobacterium sp. JLS]</td><td>280</td><td>3E-75</td></tr><tr><td>00450</td><td>Peptide methionine sulfoxide reductase [Neisseria meningitidis]</td><td>280</td><td>8E-76</td></tr><tr><td>00684</td><td>Lipid A ABC exporter [Haemophilus somnus]</td><td>273</td><td>2E-73</td></tr><tr><td>00453</td><td>Lipid A export ATP-binding protein MsbA [H. influenzae]</td><td>234</td><td>1E-61</td></tr><tr><td>00884</td><td>Putative lipoate protein ligase A [L. ivanovii subsp. ivanovii]</td><td>229</td><td>2E-60</td></tr><tr><td>00797</td><td>Variable surface lipoprotein H [M. bovis]<sup>2</sup></td><td>215</td><td>7E-74</td></tr><tr><td>00803</td><td>Variable surface lipoprotein G [M. bovis]<sup>2</sup></td><td>201</td><td>4E-66</td></tr></tbody></table>

Non-essential virulence-related proteins identified by VFDB for *M. bovis* JF4278 using BLAST. Known *M. bovis* virulence-related factors were identified by BLAST. ¹NADH oxidase (Zhao et al., 2017). ²Variable surface lipoproteins (Wise et al., 2011). ³TrmFO (Guo et al., 2017). ⁴Membrane nuclease MnuA (Sharma et al., 2015; Mitiku et al., 2018). ⁵VpmaX (Zou et al., 2013). Scores and *E*-values were generated by the BLAST algorithm. VFDB hits were filtered by BLAST scores >200.
Additionally, it promotes early escape from the phagosome of macrophages and intracellular survival by modulation of other factors at the transcription level (Rouquette et al., 1998; Nair et al., 2000; **Table 3**). Similarly, the lipoate protein ligase A (LplA1; MBOVJF4278_00884) was found to be necessary for efficient intracellular proliferation of *L. monocytogenes* (O'Riordan et al., 2003; **Table 3**). The PgPepO oligopeptidase (Zmp1; MBOVJF4278_00390) allows *Mycobacterium tuberculosis* survival in macrophages by interfering with phagosome maturation and inhibiting the formation of the inflammasome (Paolino et al., 2018; **Table 3**). The peptide methionine sulfoxide reductase (MsrA; MBOVJF4278_00450) might be involved in the protection from oxidative damage through exogenous as well as endogenous oxidants in *Neisseria meningitidis* and *M. genitalium* and adhesion in *M. genitalium* (Dhandayuthapani et al., 2001; Olry et al., 2002; **Table 3**). All these genes could have a potential role in the virulence of *M. bovis*. However, functional studies coupled with *in vitro* and *in vivo* infections are necessary to confirm or infirm these hypotheses. The identification of non-essential virulence-related factors could also lead to the design and generation of attenuated *M. bovis* strains for subsequent vaccine development. Many factors included in the VFDB are clearly linked to virulence and are well characterized as virulence attributes in the respective species. However, a clear and direct link to virulence is not obvious for all factors present in this database. For example, transposases and other mobile genetic elements are widely distributed in *Mycoplasma* spp. and are involved in genetic plasticity of genomes and horizontal gene transfer (Sirand-Pugnet et al., 2007a,b; Lysnyansky et al., 2009; **Table 3**). Indeed, insertion sequences and transposases can transfer antibiotic-resistant genes and genes encoding for virulence factors and indirectly impact the virulence of strains (Vandecraen et al., 2017; Partridge et al., 2018).

## A Set of *M. bovis* Genes Are Underrepresented After the Adhesion Assay, These Conditional-Essential Genes Might Be Adhesion-Related Factors

The full *M. bovis* strain JF4278 mutant library was screened for adhesion (**Figure 1**). The percentage of adherent *M. bovis* (output samples) relative to the total amount of mycoplasma cells used for infection (input samples) was calculated based on CFUs counted from input and output samples. The median value for mycoplasma adhesion to bMec cells from the three independent replicate experiments was 1.7% (**Figure 5**). Approximately, $5 \times 10^7$ mycoplasmas (CFU) including all the 4000 mutants were used in each replicate of the adhesion assay. With a median adhesion value of 1.7%, approximately $2 \times 10^2$ cells per mutant were recovered from the assay in each of the replicates. As mutants deficient in adhesion are expected to be more efficiently washed away during the adhesion assay, they are underrepresented in the output samples compared to other mutants. Based on this assumption, output samples were compared to input samples to identify

conditional-essential genes (genes identified as non-essential for growth but essential for adhesion).

The transposon-junctions of mutants recovered from the adhesion assay (output samples) were sequenced. Transposon sequencing data from input and output samples were used for the conditional-essentiality analysis using ESSENTIALS to identify genes putatively involved in adhesion. The read counts per gene of output samples were compared to data from input samples by ESSENTIALS. Thereby, genes from output samples are assigned a FC value of read counts in log2 scale (**Figure 6** and **Supplementary Table 6**). The log2(FC) values of JF4278 genes after the adhesion assay show a bimodal distribution in the density plot (**Figure 6**). The log2(FC) value of non-essential genes for growth and for adhesion is around zero (**Figure 6**). On the other hand, genes with underrepresented read counts after the adhesion assay [strongly negative log2(FC) values] are more likely to be conditional-essential (**Figure 6**). Significantly underrepresented (adj. *P*-value ≤ 0.05) genes below an automatically calculated log2(FC) cut-off value were considered as conditional-essential genes (Zomer et al., 2012). The log2(FC) cut-off for conditional-essential genes suggested by ESSENTIALS was <-3.38 (**Figure 6**).

According to these criteria [log2(FC) < -3.38 and adj. *P*-value ≤ 0.05], 196 out of 548 non-essential genes were found to be conditional-essential (**Supplementary Table 6**). On the other hand, transposon sequencing reads from 66 genes were significantly overrepresented after the adhesion assay (**Supplementary Table 6**). Distribution of transposon insertion sites as well as the read counts in log10 scale (height of bars) from the three sequencing replicates of the full mutant library

<figure><img src="image_7.png" /></figure>

**FIGURE 5** | Adhesion assay with the full *M. bovis* strain JF4278 mutant library. *M. bovis* adhesion to bMec cells after 30 min. The y-axis represents the percentage of adherent *M. bovis* relative to the added *M. bovis*, calculated by using the number of CFUs counted from input samples and output samples. Adhesion assay with the full *M. bovis* mutant library was performed in three independent experiments indicated by black dots. Median values are indicated as horizontal lines.
ESSENTIALS (Conditional)

Epanechnikov Kernel

<figure><img src="image_8.png" /></figure>

**FIGURE 6** | Density plot of conditional-essentiality values obtained with the software ESSENTIALS. Read counts per gene of output samples were compared to data from input samples by ESSENTIALS (Zomer et al., 2012). *M. bovis* strain JF4278 genes from output samples are assigned a fold change (FC) value of read counts in log2 scale compared to the read counts of input samples. Log2(FC) values for all genes are shown on the x-axis of the density plot. Log2(FC) values show a bimodal distribution. The log2(FC) value of non-essential JF4278 genes is around zero. Genes more likely to be conditional-essential (essential for adhesion) have negative log2(FC) values. The cut-off value for conditional-essential genes suggested by the software ESSENTIALS for this data set is log2(FC) < -3.38.

(input samples) and of the mutants recovered after the adhesion assay (output samples) are in line with the predicted conditional-essentiality of genes (Figure 7 and Supplementary Table 6). Genes with fewer transposon insertions and smaller numbers of reads in the output samples (Figure 7, Outputs) compared to input samples (Figure 7, Inputs) were identified as conditional-essential (Figure 7, purple arrows). In addition to the prediction of adhesion-related factors using ESSENTIALS, SPAAN analysis of the predicted proteome of *M. bovis* was performed (Sachdeva et al., 2005). Putative proteins were assigned a *P*<sub>ad</sub> value reflecting the probability of a protein to be an adhesin, with higher values reflecting higher probabilities (Supplementary Table 6).

The set of genes identified to be conditional-essential
represents a list of potential adhesion-related factors in *M. bovis*.
For example, the NADH oxidase (MBOVJF4278_00281) and
the TrmFO protein (MBOVJF4278_00132), which are known
*M. bovis* adhesins (Guo et al., 2017; Zhao et al., 2017),
were predicted to be conditional-essential and involved in
adhesion (Supplementary Table 6). Endopeptidase Clp ATP-
binding chain C (MBOVJF4278_00700) identified with the VFDB
and known to be involved in adhesion in *L. monocytogenes*
(Nair et al., 2000; Table 3) was assigned as a putative
adhesion-related factor in *M. bovis* (Supplementary Table 6).
However, ESSENTIALS failed to predict most of the Vsps
and VpmaX (MBOVJF4278_00429) to be conditional-essential
in the adhesion assay (Supplementary Table 6). This is
probably due to the repetitive regions of *vsp*-genes. Indeed,
sequence reads from repetitive regions cannot be assigned to
a specific position in the genome and are not considered by
ESSENTIALS. Moreover, the *M. bovis* adhesion factor α-enolase
(MBOVJF4278_00443) previously shown to be essential could

not be assigned to a conditional-essential phenotype (Song et al.,
2012; Supplementary Tables 2, 6).

The combined information from adhesion-related factors identified by ESSENTIALS and SPAAN is useful to identify candidate genes for future adhesion studies (Supplementary Table 6). Nevertheless, false positives cannot be ruled out and therefore candidate mutants for adhesion should be individually tested. As mentioned above, approximately 2 × 10² cells per mutant were recovered from the adhesion assay in each of the replicates. This rather small number of recovered cells per mutant represents a potential bottleneck during the adhesion assay. Especially mutants that grew slowly during the collection of individual mutants and are underrepresented in the full mutant library, even before the infection, could lead to false positive adhesion candidates. To prevent a potential bottleneck during the adhesion assay, the number of different mutants used could be reduced to recover more than 2 × 10² cells per mutant after adhesion. However, including fewer mutants in this study would have hampered the accurate assignment of essential and non-essential genes from input samples. For these reasons, conditional-essential genes with the lowest log2(FC) values (Supplementary Table 6) and with the highest adhesion probability values by SPAAN (Supplementary Table 6) are top candidate genes for adhesion to be tested in future studies.

Three Individual Mutants Showed a
Reduced Adhesion Phenotype

Six candidate mutants from the complete spectrum of putative
adhesion-related factors identified by ESSENTIALS and
<figure><img src="image_9.png" /><figcaption>**FIGURE 7** | Overview of non-essential, essential, and conditional-essential genes in the genome of *M. bovis* strain JF4278. Genes and their orientation are indicated by arrows. Genes found essential with ESSENTIALS are filled in red, while non-essential genes are filled in yellow (Zomer et al., 2012). Genes found to be conditional-essential (essential for adhesion) are filled in purple. Essentiality was assessed with transposon sequencing data of the full mutant library (input samples) with ESSENTIALS. Conditional-essentiality was assessed by comparing input samples to output samples with ESSENTIALS. Mapping of the transposon integrations sites from the transposon sequencing of the input samples (Inputs) and output samples (Outputs) are represented by black, green, and red bars (each corresponding to one replicate). The height of the insertion bars represents the read count per insertion.</figcaption></figure>

SPAAN were selected to be individually evaluated for adhesion: MBOVJF4278_00132, MBOVJF4278_00255, MBOVJF4278_00264, MBOVJF4278_00598, MBOVJF4278_00667, and

MBOVJF4278_00812 (Table 4). A mutant of the gene coding for the TrmFO protein (MBOVJF4278_00132) was selected because this protein was previously characterized as an adhesin
in *M. bovis* (Guo et al., 2017), despite the fact that it was not among the top candidates identified with ESSENTIALS and SPAAN analysis (Table 4 and Supplementary Table 6). According to TMPred TrmFO displays one putative TM-H with a score of 424 which is just below the suggested cut-off value of 500. The hypothetical protein MBOVJF4278_00255 is a homolog of the putative lipoprotein MBOVPG45_0586 of PG45 that belongs to the LppB family. MBOVJF4278_00255 was among the top candidates for adhesion identified with ESSENTIALS and SPAAN and was predicted to contain two TM-Hs (Table 4 and Supplementary Table 6). The hypothetical protein MBOVJF4278_00264 is a homolog of the putative lipoprotein MBOVPG45_0574 in PG45 and belongs to the PARCEL family. MBOVJF4278_00264 was identified with ESSENTIALS and SPAAN as an intermediate candidate for adhesion and was predicted to contain one TM-H (Table 4 and Supplementary Table 6). The putative phosphatase

MBOVJF4278_00598 is a homolog of MBOVPG45_0281 in PG45 and belongs to the HAD-superfamily hydrolase, subfamily IIB. Based on SPAAN and ESSENTIALS, this protein displayed a high probability to be involved in adhesion but no TM-H was predicted (Table 4 and Supplementary Table 6). The hypothetical protein MBOVJF4278_00667 was also among the top candidates for adhesion identified with ESSENTIALS. However, the SPAAN value was notably low (Table 4 and Supplementary Table 6). The TMPred software predicted one TM-H for this protein (Table 4). MBOVJF4278_00812, encoding for a hypothetical protein, was the only gene significantly underrepresented from the *vsp*-locus (Supplementary Table 6). The sequence obtained from the genome of JF4278 indicates that MBOVJF4278_00812 is a truncated form of the VspN in PG45 (MBOVPG45_0811). Even though adhesion-prediction by ESSENTIALS and SPAAN was low (Table 4 and Supplementary Table 6),

<figure><img src="image_10.png" /><figcaption>**FIGURE 8** | Adhesion assay of individual candidate mutants of *M. bovis* strain JF4278. *M. bovis* adhesion to bMec cells after 30 min. The y-axis represents the percentage of adherent *M. bovis* relative to the added *M. bovis*. Adhesion with the wild-type strain or single mutants was performed in four independent experiments as indicated by black dots. Median values are indicated as horizontal lines in the box plot. Statistical comparison of each mutant to the wild-type strain was performed by the two-sample Wilcoxon test. \* *P* < 0.05.</figcaption></figure>

<table><caption>**TABLE 4** | Features of the *M. bovis* potential adhesion-deficient candidate mutants.</caption><thead><tr><th>Candidate mutant (MBOVJF4278_)</th><th>Product</th><th>Tn insertion site in gene (%)</th><th>ESSENTIALS log2(FC)<sup>1</sup></th><th>Adhesion probability SPAAN<sup>2</sup></th><th>TM-helices<sup>3</sup></th></tr></thead><tbody><tr><td>00132</td><td>TrmFO</td><td>50.8</td><td>-4.0</td><td>0.23</td><td>0</td></tr><tr><td>00255</td><td>Hypothetical protein</td><td>19.8</td><td>-8.7</td><td>0.71</td><td>2</td></tr><tr><td>00264</td><td>Hypothetical protein</td><td>65.2</td><td>-7.9</td><td>0.48</td><td>1</td></tr><tr><td>00598</td><td>Putative phosphatase MPN_383</td><td>61.8</td><td>-8.1</td><td>0.70</td><td>0</td></tr><tr><td>00667</td><td>Hypothetical protein</td><td>25.6</td><td>-9.2</td><td>0.26</td><td>1</td></tr><tr><td>00812</td><td>Hypothetical protein</td><td>66.4</td><td>-5.0</td><td>0.54</td><td>1</td></tr></tbody></table>

<footnote n="1">¹Conditional-essentiality log2(FC) value by ESSENTIALS (Zomer et al., 2012). ²Adhesion probability value by SPAAN (Sachdeva et al., 2005). ³Number of transmembrane helices predicted by TMPred scores >500.</footnote>
MBOVJF4278_00812 was included as a candidate for adhesion in agreement with the Vsps being previously shown to be involved in adhesion (Sachse et al., 1993, 1996, 2000; Thomas et al., 2003).

Adhesion of the individual *M. bovis* mutants to bMec cells was compared to wt strain JF4278. The median adhesion was 1.8% for the wt strain JF4278 (Figure 8). Three out of six mutants showed significantly lower adhesion capacity to bMec cells compared to the wt: a mutant of the gene coding for the TrmFO protein (MBOVJF4278_00132) and mutants of two genes coding for hypothetical proteins (MBOVJF4278_00255 and MBOVJF4278_00667) (Figure 8). The MBOVJF4278_00132 mutant showed a median adhesion of 1.2% (Figure 8). Both of the two newly identified adhesion-related factors had a median adhesion of 0.8%, which is over 50% lower than the wt strain JF4278 (Figure 8). The median adhesion value of the other three selected mutants, MBOVJF4278_00264, MBOVJF4278_00598, and MBOVJF4278_0081, was not significantly different from the wt strain (Figure 8).

MBOVJF4278_00255 and MBOVJF4278_00667 were among the top candidates identified with ESSENTIALS (Supplementary Table 6). Moreover, the presence of at least one TM-H was predicted for both genes (Table 4). MBOVJF4278_00255 was also among the top candidates for adhesion by the SPAAN analysis (Supplementary Table 6). Therefore, a focus on top candidates and a combination of the results obtained with ESSENTIALS, SPAAN, and TMpred could improve the prediction of putative adhesins. Nevertheless, it is crucial to test single mutants to confirm results obtained via high-throughput approaches. The final assignment of an adhesin should then include experiments to complement mutants. Moreover, the surface exposure of these putative adhesins should be characterized *in vitro*. An additional point to consider is the expression of Vsps in the wt strain and in mutants. As the *vsp*-locus undergoes recombination events, the expression of Vsps might differ between a mutant and the wt strain, potentially affecting adhesion.

## CONCLUSION

Addressing gene essentiality is always dependent on the experimental settings, in particular the conditions used for preparing and propagating mutant libraries and the analysis method chosen. In recent years, methods to assign essentiality from transposon mutant libraries have developed, from looking at whether a gene was disrupted by a transposon, to applying rules for disruption, to high-throughput analysis approaches involving statistics for essentiality probabilities (Christen et al., 2011; Zomer et al., 2012; Lluch-Senar et al., 2015; Solaimanpour et al., 2015). Similar to ESSENTIALS, these methods also require high-density mutant libraries where few transposon insertions can be detected in essential genes, but this is counterbalanced by taking into account the amount of associated reads. Additionally, one recurrent question is whether the simultaneous disruption

of non-essential genes identified by transposon mutagenesis results in viable bacteria and whether some mutants are mutually exclusive (Hutchison et al., 2016). Data collected in this study resulted in proposing a list of individual *M. bovis* genes tolerating disruption among which several non-essential putative virulence-related factors were identified. Moreover, three clones showed reduced bovine cell adhesion. Overall, datasets of non-essential and conditional-essential genes are both a source of valuable information for further exploring *M. bovis* biology, pathogenesis, and interaction with the host as well as for vaccine design.

## DATA AVAILABILITY

The transposon sequencing datasets analyzed for this study can be found in the ENA (http://www.ebi.ac.uk/ena/data/view/PRJEB31021).

## AUTHOR CONTRIBUTIONS

CJ, SB, ED-F, CC, LF, and PP designed the experiments. CJ, SB, and SV performed the experiments. CJ, LF, and PP analyzed the data. CJ drafted the manuscript. All authors helped in writing the manuscript and critically revised it. All authors read and approved the final manuscript.

## FUNDING

This study was supported by the Swiss National Science Foundation (Reference No. 31003A_160159 to PP) and by the research fund of the Institute of Veterinary Bacteriology, University of Bern, Bern, Switzerland. The open access publication fees were supported by the Institute of Veterinary Bacteriology, University of Bern, Bern, Switzerland.

## ACKNOWLEDGMENTS

We are grateful to Aldert Zomer (Utrecht University, Netherlands), who was very helpful for the analysis with the software ESSENTIALS. We thank Olga Wellnitz and Denisa Dan (University of Bern, Switzerland) for the help in isolating primary bovine cells. We are also grateful to Lisandra Aguilar-Bultet (University of Basel, Switzerland) for the introduction into UNIX.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.02085/full#supplementary-material
REFERENCES

Aboklaish, A. F., Dordet-Frisoni, E., Citti, C., Toleman, M. A., Glass, J. I., and Spiller, O. B. (2014). Random insertion and gene disruption via transposon mutagenesis of *Ureaplasma parvum* using a mini-transposon plasmid. *Int. J. Med. Microbiol.* 304, 1218–1225. doi: 10.1016/j.ijmm.2014.09.003

Acevedo-Rocha, C. G., Fang, G., Schmidt, M., Ussery, D. W., and Danchin, A. (2013). From essential to persistent genes: a functional approach to constructing synthetic life. *Trends Genet.* 29, 273–279. doi: 10.1016/j.tig.2012.11.001

Aebi, M., Bodmer, M., Frey, J., and Pilo, P. (2012). Herd-specific strains of *Mycoplasma bovis* in outbreaks of mycoplasmal mastitis and pneumonia. *Vet. Microbiol.* 157, 363–368. doi: 10.1016/j.vetmic.2012.01.006

Afgan, E., Baker, D., Batut, B., Van Den Beek, M., Bouvier, D., Cech, M., et al. (2018). The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2018 update. *Nucleic Acids Res.* 46, W537–W544. doi: 10.1093/nar/gky379

Baranowski, E., Guiral, S., Sagne, E., Skapski, A., and Citti, C. (2010). Critical role of dispensable genes in *Mycoplasma agalactiae* interaction with mammalian cells. *Infect. Immun.* 78, 1542–1551. doi: 10.1128/IAI.01195-09

Barnett, D. W., Garrison, E. K., Quinlan, A. R., Stromberg, M. P., and Marth, G. T. (2011). BamTools: a C++ API and toolkit for analyzing and managing BAM files. *Bioinformatics* 27, 1691–1692. doi: 10.1093/bioinformatics/btr174

Behrens, A., Heller, M., Kirchhoff, H., Yogev, D., and Rosengarten, R. (1994). A family of phase- and size-variant membrane surface lipoprotein antigens (Vsps) of *Mycoplasma bovis*. *Infect. Immun.* 62, 5075–5084.

Beven, L., Charenton, C., Dautant, A., Bouyssou, G., Labroussaa, F., Skollermo, A., et al. (2012). Specific evolution of F1-Like ATPases in mycoplasmas. *PLoS One* 7:e38793. doi: 10.1371/journal.pone.0038793

Blötz, C., and Stülke, J. (2017). Glycerol metabolism and its implication in virulence in *Mycoplasma*. *FEMS Microbiol. Rev.* 41, 640–652. doi: 10.1093/femsre/fux033

Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. *Bioinformatics* 30, 2114–2120. doi: 10.1093 bioinformatics/btu170

Buchenau, I., Poumarat, F., Le Grand, D., Linkner, H., Rosengarten, R., and Hewicker-Trautwein, M. (2010). Expression of *Mycoplasma bovis* variable surface membrane proteins in the respiratory tract of calves after experimental infection with a clonal variant of *Mycoplasma bovis* type strain PG45. *Res. Vet. Sci.* 89, 223–229. doi: 10.1016/j.rvsc.2010.03.014

Bürki, S., Frey, J., and Pilo, P. (2015). Virulence, persistence and dissemination of *Mycoplasma bovis*. *Vet. Microbiol.* 179, 15–22. doi: 10.1016/j.vetmic.2015.02.024

Calcutt, M. J., Lysnyansky, I., Sachse, K., Fox, L. K., Nicholas, R. A. J., and Ayling, R. D. (2018). Gap analysis of *Mycoplasma bovis* disease, diagnosis and control: an aid to identify future development requirements. *Transbound. Emerg. Dis.* 65, 91–109. doi: 10.1111/tbed.12860

Chen, L., Zheng, D., Liu, B., Yang, J., and Jin, Q. (2016). VFDB 2016: hierarchical and refined dataset for big data analysis–10 years on. *Nucleic Acids Res.* 44, D694–D697. doi: 10.1093/nar/gkv1239

Chopra-Dewasthaly, R., Zimmermann, M., Rosengarten, R., and Citti, C. (2005). First steps towards the genetic manipulation of *Mycoplasma agalactiae* and *Mycoplasma bovis* using the transposon Tn4001mod. *Int. J. Med. Microbiol.* 294, 447–453. doi: 10.1016/j.ijmm.2004.09.010

Christen, B., Abeliuk, E., Collier, J. M., Kalogeraki, V. S., Passarelli, B., Coller, J. A., et al. (2011). The essential genome of a bacterium. *Mol. Syst. Biol.* 7:528. doi: 10.1038/msb.2011.58

Dawoud, T. M., Jiang, T., Mandal, R. K., Ricke, S. C., and Kwon, Y. M. (2014). Improving the efficiency of transposon mutagenesis in *Salmonella enteritidis* by overcoming host-restriction barriers. *Mol. Biotechnol.* 56, 1004–1010. doi: 10.1007/s12033-014-9779-4

Dhandayuthapani, S., Blaylock, M. W., Bebear, C. M., Rasmussen, W. G., and Baseman, J. B. (2001). Peptide methionine sulfoxide reductase (MsrA) is a virulence determinant in *Mycoplasma genitalium*. *J. Bacteriol.* 183, 5645–5650. doi: 10.1128/JB.183.19.5645-5650.2001

Dordet Frisoni, E., Marenda, M. S., Sagne, E., Nouvel, L. X., Guerillot, R., Glaser, P., et al. (2013). ICEA of *Mycoplasma agalactiae*: a new family of self-transmissible integrative elements that confers conjugative properties to the recipient strain. *Mol. Microbiol.* 89, 1226–1239. doi: 10.1111/mmi.12341

Dybvig, K., French, C. T., and Voelker, L. L. (2000). Construction and use of derivatives of transposon Tn4001 that function in *Mycoplasma pulmonis* and *Mycoplasma arthritidis*. *J. Bacteriol.* 182, 4343–4347. doi: 10.1128/Jb.182.15.4343-4347.2000

Dybvig, K., Lao, P., Jordan, D. S., and Simmons, W. L. (2010). Fewer essential genes in mycoplasmas than previous studies suggest. *FEMS Microbiol. Lett.* 311, 51–55. doi: 10.1111/j.1574-6968.2010.02078.x

French, C. T., Lao, P., Loraine, A. E., Matthews, B. T., Yu, H., and Dybvig, K. (2008). Large-scale transposon mutagenesis of *Mycoplasma pulmonis*. *Mol. Microbiol.* 69, 67–76. doi: 10.1111/j.1365-2958.2008.06262.x

Freundt, E. A. (1983). "Culture media for classic mycoplasmas," in *Methods in Mycoplasmosis*, eds S. Razin, and J. G. Tully, (New York, NY: Academic Press), 127–135. doi: 10.1016/b978-0-12-583801-6.50029-9

Gabriel, R., Eckenberg, R., Paruzynski, A., Bartholomae, C. C., Nowrouzi, A., Arens, A., et al. (2009). Comprehensive genomic access to vector integration in clinical gene therapy. *Nat. Med.* 15, 1431–1436. doi: 10.1038/nm.2057

Gautier-Bouchardon, A. V., Ferre, S., Le Grand, D., Paoli, A., Gay, E., and Poumarat, F. (2014). Overall decrease in the susceptibility of *Mycoplasma bovis* to antimicrobials over the past 30 years in France. *PLoS One* 9:e87672. doi: 10.1371/journal.pone.0087672

Glass, J. I., Assad-Garcia, N., Alperovich, N., Yooseph, S., Lewis, M. R., Maruf, M., et al. (2006). Essential genes of a minimal bacterium. *Proc. Natl. Acad. Sci. U.S.A.* 103, 425–430. doi: 10.1073/pnas.0510013103

Grazziotin, A. L., Vidal, N. M., and Venancio, T. M. (2015). Uncovering major genomic features of essential genes in Bacteria and a methanogenic Archaea. *FEBS J.* 282, 3395–3411. doi: 10.1111/febs.13350

Guo, Y., Zhu, H., Wang, J., Huang, J., Khan, F. A., Zhang, J., et al. (2017). TrmFO, a fibronectin-binding adhesin of *Mycoplasma bovis*. *Int. J. Mol. Sci.* 18:1732. doi: 10.3390/ijms18081732

Hedreyda, C. T., Lee, K. K., and Krause, D. C. (1993). Transformation of *Mycoplasma pneumoniae* with Tn4001 by electroporation. *Plasmid* 30, 170–175. doi: 10.1006/plas.1993.1047

Hegde, S., Zimmermann, M., Flöck, M., Brunthaler, R., Spergser, J., Rosengarten, R., et al. (2016). Genetic loci of *Mycoplasma agalactiae* involved in systemic spreading during experimental intramammary infection of sheep. *Vet. Res.* 47:106. doi: 10.1186/s13567-016-0387-0

Hutchison, C. A. III, Chuang, R. Y., Noskov, V. N., Assad-Garcia, N., Deerinck, T. J., Ellisman, M. H., et al. (2016). Design and synthesis of a minimal bacterial genome. *Science* 351:ad6253. doi: 10.1126/science.ad6253

Hutchison, C. A., Peterson, S. N., Gill, S. R., Cline, R. T., White, O., Fraser, C. M., et al. (1999). Global transposon mutagenesis and a minimal *Mycoplasma* genome. *Science* 286, 2165–2169. doi: 10.1126/science.286.5447.2165

Josi, C., Bürki, S., Stojiljkovic, A., Wellnitz, O., Stoffel, M. H., and Pilo, P. (2018). Bovine epithelial *in vitro* infection models for *Mycoplasma bovis*. *Front. Cell. Infect. Microbiol.* 8:329. doi: 10.3389/fcimb.2018.00329

Kearse, M., Moir, R., Wilson, A., Stones-Havas, S., Cheung, M., Sturrock, S., et al. (2012). Geneious basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. *Bioinformatics* 28, 1647–1649. doi: 10.1093/bioinformatics/bts199

Kimura, S., Hubbard, T. P., Davis, B. M., and Waldor, M. K. (2016). The nucleoid binding protein H-NS biases genome-wide transposon insertion landscapes. *Mbio* 7:e01351-16. doi: 10.1128/mBio.01351-16

Langmead, B., and Salzberg, S. L. (2012). Fast gapped-read alignment with Bowtie 2. *Nat. Methods* 9, 357–359. doi: 10.1038/nmeth.1923

Lazinski, D. W., and Camilli, A. (2013). Homopolymer tail-mediated ligation PCR: a streamlined and highly efficient method for DNA cloning and library construction. *Biotechniques* 54, 25–34. doi: 10.2144/000113981

Lin, Y., and Zhang, R. R. (2011). Putative essential and core-essential genes in *Mycoplasma* genomes. *Sci. Rep.* 1:53. doi: 10.1038/srep00053

Liu, W., Fang, L., Li, M., Li, S., Guo, S., Luo, R., et al. (2012). Comparative genomics of mycoplasma: analysis of conserved essential genes and diversity of the pan-genome. *PLoS One* 7:e35698. doi: 10.1371/journal.pone.0035698

Lluch-Senar, M., Delgado, J., Chen, W. H., Llorens-Rico, V., O'Reilly, F. J., Wodke, J. A., et al. (2015). Defining a minimal cell: essentiality of small ORFs and ncRNAs in a genome-reduced bacterium. *Mol. Syst. Biol.* 11:780. doi: 10.15252/msb.20145558

Lysnyansky, I., Calcutt, M. J., Ben-Barak, I., Ron, Y., Levisohn, S., Methe, B. A., et al. (2009). Molecular characterization of newly identified IS3, IS4 and IS30
insertion sequence-like elements in *Mycoplasma bovis* and their possible roles
in genome plasticity. *FEMS Microbiol. Lett.* 294, 172–182. doi: 10.1111/j.1574-6968.2009.01562.x

Lysnyansky, I., Ron, Y., Sachse, K., and Yogev, D. (2001). Intrachromosomal recombination within the vsp locus of *Mycoplasma bovis* generates a chimeric variable surface lipoprotein antigen. *Infect. Immun.* 69, 3703–3712. doi: 10.1128/IAI.69.6.3703-3712.2001

Lysnyansky, I., Rosengarten, R., and Yogev, D. (1996). Phenotypic switching of variable surface lipoproteins in *Mycoplasma bovis* involves high-frequency chromosomal rearrangements. *J. Bacteriol.* 178, 5395–5401. doi: 10.1128/jb.178.18.5395-5401.1996

Lysnyansky, I., Sachse, K., Rosenbusch, R., Levisohn, S., and Yogev, D. (1999). The vsp locus of *Mycoplasma bovis*: gene organization and structural features. *J. Bacteriol.* 181, 5734–5741.

Masukagami, Y., De Souza, D. P., Dayalan, S., Bowen, C., O'callaghan, S., Kouremenos, K., et al. (2017). Comparative metabolomics of *Mycoplasma bovis* and *Mycoplasma gallisepticum* reveals fundamental differences in active metabolic pathways and suggests novel gene annotations. *mSystems* 2:e00055-17. doi: 10.1128/mSystems.00055-17

May, M., Balish, M. F., and Blanchard, A. (2014). “The order *Mycoplasma matalense*,” in *The Prokaryotes: Firmicutes and Tenericutes*, eds E. Rosenberg, E. F. Delong, S. Lory, E. Stackebrandt, and F. Thompson, (Berlin: Springer), 515–550. doi: 10.1007/978-3-642-30120-9_289

Mitiku, F., Hartley, C. A., Sansom, F. M., Coombe, J. E., Mansell, P. D., Beggs, D. S., et al. (2018). The major membrane nuclease MnuA degrades neutrophil extracellular traps induced by *Mycoplasma bovis*. *Vet. Microbiol.* 218, 13–19. doi: 10.1016/j.vetmic.2018.03.002

Nair, S., Milohanic, E., and Berche, P. (2000). ClpC ATPase is required for cell adhesion and invasion of *Listeria monocytogenes*. *Infect. Immun.* 68, 7061–7068. doi: 10.1128/lai.68.12.7061-7068.2000

Nicholas, R. A. (2011). Bovine mycoplasmosis: silent and deadly. *Vet. Rec.* 168, 459–462. doi: 10.1136/vr.d2468

Nicholson, P., Joncourt, R., and Muhlemann, O. (2012). Analysis of nonsense-mediated mRNA decay in mammalian cells. *Curr. Protoc. Cell Biol.* 27, Unit 27.4. doi: 10.1002/0471143030.cb2704s55

Olyr, A., Boschi-Muller, S., Marraud, M., Sanglier-Cianferani, S., Van Dorselear, A., and Branlant, G. (2002). Characterization of the methionine sulfoxide reductase activities of PILB, a probable virulence factor from *Neisseria meningitidis*. *J. Biol. Chem.* 277, 12016–12022. doi: 10.1074/jbc.M112350200

O'Riordan, M., Moors, M. A., and Portnoy, D. A. (2003). *Listeria* intracellular growth and virulence require host-derived lipoic acid. *Science* 302, 462–464. doi: 10.1126/science.1088170

Paolino, M., Brindisi, M., Vallone, A., Butini, S., Campiani, G., Nannicini, C., et al. (2018). Development of potent inhibitors of the *Mycobacterium tuberculosis* virulence factor Zmp1 and evaluation of their effect on mycobacterial survival inside macrophages. *ChemMedChem* 13, 422–430. doi: 10.1002.201700759

Partridge, S. R., Kwong, S. M., Firth, N., and Jensen, S. O. (2018). Mobile genetic elements associated with antimicrobial resistance. *Clin. Microbiol. Rev.* 31:e00088-17. doi: 10.1128/CMR.00088-17

Paruzynski, A., Arens, A., Gabriel, R., Bartholomae, C. C., Scholz, S., Wang, W., et al. (2010). Genome-wide high-throughput integrome analyses by nrLAM-PCR and next-generation sequencing. *Nat. Protoc.* 5, 1379–1395. doi: 10.1038/nprot.2010.87

Perez-Casal, J., Prysliak, T., Maina, T., Suleman, M., and Jimbo, S. (2017). Status of the development of a vaccine against *Mycoplasma bovis*. *Vaccine* 35, 2902–2907. doi: 10.1016/j.vaccine.2017.03.095

Rasheed, M. A., Qi, J., Zhu, X., Chenfei, H., Menghwar, H., Khan, F. A., et al. (2017). Comparative genomics of *Mycoplasma bovis* strains reveals that decreased virulence with increasing passages might correlate with potential virulence-related factors. *Front. Cell. Infect. Microbiol.* 7:177. doi: 10.3389/fcimb.2017.00177

Razin, S., Yogev, D., and Naot, Y. (1998). Molecular biology and pathogenicity of mycoplasmas. *Microbiol. Mol. Biol. Rev.* 62, 1094–1156.

Renaudin, J., Breton, M., and Citti, C. (2014). “Molecular genetic tools for mollicutes,” in *Mollicutes: Molecular Biology and Pathogenesis*, eds G. F. Browning, and C. Citti, (Wymondham: Caister Academic Press), 55–76.

Robinson, M. D., McCarthy, D. J., and Smyth, G. K. (2010). edgeR: a bioconductor package for differential expression analysis of digital gene expression data. *Bioinformatics* 26, 139–140. doi: 10.1093/bioinformatics/btp616

Rocha, E. P. C., Cornet, E., and Michel, B. (2005). Comparative and evolutionary analysis of the bacterial homologous recombination systems. *PLoS Genet.* 1:e15. doi: 10.1371/journal.pgen.0010015

Rosati, S., Robino, P., Fadda, M., Pozzi, S., Mannelli, A., and Pittau, M. (2000). Expression and antigenic characterization of recombinant *Mycoplasma agalactiae* P48 major surface protein. *Vet. Microbiol.* 71, 201–210. doi: 10.1016/S0378-1135(99)00164-9

Rosengarten, R., Behrens, A., Stetefeld, A., Heller, M., Ahrens, M., Sachse, K., et al. (1994). Antigen heterogeneity among isolates of *Mycoplasma bovis* is generated by high-frequency variation of diverse membrane surface proteins. *Infect. Immun.* 62, 5066–5074.

Rossetti, B. C., Frey, J., and Pilo, P. (2010). Direct detection of *Mycoplasma bovis* in milk and tissue samples by real-time PCR. *Mol. Cell. Probes* 24, 321–323. doi: 10.1016/j.mcp.2010.05.001

Rouquette, C., De Chastellier, C., Nair, S., and Berche, P. (1998). The ClpC ATPase of *Listeria monocytogenes* is a general stress protein required for virulence and promoting early bacterial escape from the phagosome of macrophages. *Mol. Microbiol.* 27, 1235–1245. doi: 10.1046/j.1365-2958.1998.00775.x

Sachdeva, G., Kumar, K., Jain, P., and Ramachandran, S. (2005). SPAAN: a software program for prediction of adhesins and adhesin-like proteins using neural networks. *Bioinformatics* 21, 483–491. doi: 10.1093/bioinformatics/bti028

Sachse, K. (1998). Detection and analysis of mycoplasma adhesins. *Methods Mol. Biol.* 104, 299–307. doi: 10.1385/0-89603-525-5:299

Sachse, K., Grajetzki, C., Rosengarten, R., Hanel, I., Heller, M., and Pfutzner, H. (1996). Mechanisms and factors involved in *Mycoplasma bovis* adhesion to host cells. *Zentralbl. Bakteriol.* 284, 80–92. doi: 10.1016/S0934-8840(96)80157-5

Sachse, K., Helbig, J. H., Lysnyansky, I., Grajetzki, C., Muller, W., Jacobs, E., et al. (2000). Epitope mapping of immunogenic and adhesive structures in repetitive domains of *Mycoplasma bovis* variable surface lipoproteins. *Infect. Immun.* 68, 680–687. doi: 10.1128/IAI.68.2.680-687.2000

Sachse, K., Pfutzner, H., Heller, M., and Hanel, I. (1993). Inhibition of *Mycoplasma bovis* cytidherence by a monoclonal antibody and various carbohydrate substances. *Vet. Microbiol.* 36, 307–316. doi: 10.1016/0378-1135(93)90097-Q

Sharma, S., Markham, P. F., and Browning, G. F. (2014). Genes found essential in other mycoplasmas are dispensable in *Mycoplasma bovis*. *PLoS One* 9:e97100. doi: 10.1371/journal.pone.0097100

Sharma, S., Tivendale, K. A., Markham, P. F., and Browning, G. F. (2015). Disruption of the membrane nuclease gene (MBOVPG45_0215) of *Mycoplasma bovis* greatly reduces cellular nuclease activity. *J. Bacteriol.* 197, 1549–1558. doi: 10.1128/JB.00034-15

Sirand-Pugnet, P., Citti, C., Barre, A., and Blanchard, A. (2007a). Evolution of mollicutes: down a bumpy road with twists and turns. *Res. Microbiol.* 158, 754–766. doi: 10.1016/j.resmic.2007.09.007

Sirand-Pugnet, P., Lartigue, C., Marenda, M., Jacob, D., Barre, A., Barbe, V., et al. (2007b). Being pathogenic, plastic, and sexual while living with a nearly minimal bacterial genome. *PLoS Genet.* 3:e75. doi: 10.1371/journal.pgen.0030075

Solaimanpour, S., Sarmiento, F., and Mrázek, J. (2015). Tn-Seq explorer: a tool for analysis of high-throughput sequencing data of transposon mutant libraries. *PLoS One* 10:e0126070. doi: 10.1371/journal.pone.0126070

Song, Z., Li, Y., Liu, Y., Xin, J., Zou, X., and Sun, W. (2012). α-Enolase, an adhesion-related factor of *Mycoplasma bovis*. *PLoS One* 7:e38836. doi: 10.1371/journal.pone.0038836

Thomas, A., Sachse, K., Farnir, F., Dizier, I., Mainil, J., and Linden, A. (2003). Adherence of *Mycoplasma bovis* to bovine bronchial epithelial cells. *Microb. Pathog.* 34, 141–148. doi: 10.1016/S0882-4010(03)00003-2

Tseng, C. W., Kanci, A., Citti, C., Rosengarten, R., Chiu, C. J., Chen, Z. H., et al. (2013). Malf is essential for persistence of *Mycoplasma gallisepticum* in vivo. *Microbiology* 159, 1459–1470. doi: 10.1099/mic.0.067553-0

Vandecraen, J., Chandler, M., Aertsen, A., and Van Houdt, R. (2017). The impact of insertion sequences on bacterial genome plasticity and adaptability. *Critic. Rev. Microbiol.* 43, 709–730. doi: 10.1080/1040841X.2017.1303661

Wellnitz, O., and Kerr, D. E. (2004). Cryopreserved bovine mammary cells to model epithelial response to infection. *Vet. Immunol. Immunopathol.* 101, 191–202. doi: 10.1016/j.vetimm.2004.04.019
Wise, K. S., Calcutt, M. J., Foecking, M. F., Roske, K., Madupu, R., and Methe, B. A. (2011). Complete genome sequence of *Mycoplasma bovis* type strain PG45 (ATCC 25523). *Infect. Immun.* **79**, 982–983. doi: 10.1128/IAI.00726-10

Xie, C., Mao, X., Huang, J., Ding, Y., Wu, J., Dong, S., et al. (2011). KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. *Nucleic Acids Res.* **39**, W316–W322. doi: 10.1093/nar/gkr483

Zbinden, C., Pilo, P., Frey, J., Bruckmaier, R. M., and Wellnitz, O. (2015). The immune response of bovine mammary epithelial cells to live or heat-inactivated *Mycoplasma bovis*. *Vet. Microbiol.* **179**, 336–340. doi: 10.1016/j.vetmic.2015.07.007

Zhang, H., Zhao, G., Guo, Y., Menghwar, H., Chen, Y., Chen, H., et al. (2016). *Mycoplasma bovis* MBOV_RS02825 Encodes a secretory nuclease associated with cytotoxicity. *Int. J. Mol. Sci.* **17**:E628. doi: 10.3390/ijms17050628

Zhao, G., Zhang, H., Chen, X., Zhu, X., Guo, Y., He, C., et al. (2017). *Mycoplasma bovis* NADH oxidase functions as both a NADH oxidizing and O2 reducing enzyme and an adhesin. *Sci. Rep.* **7**:44. doi: 10.1038/s41598-017-00121-y

Zimmerman, C. U., and Herrmann, R. (2005). Synthesis of a small, cysteine-rich, 29 amino acids long peptide in *Mycoplasma pneumoniae*. *FEMS Microbiol. Lett.* **253**, 315–321. doi: 10.1016/j.femsle.2005.09.054

Zomer, A., Burghout, P., Bootsma, H. J., Hermans, P. W. M., and Van Hijum, S. A. F. T. (2012). ESSENTIALS: software for rapid analysis of high throughput transposon insertion sequencing data. *PLoS One* **7**:e43012. doi: 10.1371/journal.pone.0043012

Zou, X., Li, Y., Wang, Y., Zhou, Y., Liu, Y., and Xin, J. (2013). Molecular cloning and characterization of a surface-localized adhesion protein in *Mycoplasma bovis* Hubei-1 strain. *PLoS One* **8**:e69644. doi: 10.1371/journal.pone.0069644

**Conflict of Interest Statement:** 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.

Copyright © 2019 Josi, Bürki, Vidal, Dordet-Frisoni, Citti, Falquet and Pilo. 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.