# Deciphering the role of insertion sequences in the evolution of bacterial epidemic pathogens with panISa software

Charlotte Couchoud, Xavier Bertrand, Benoît Valot, Didier Hocquet

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

## To cite this version:

Charlotte Couchoud, Xavier Bertrand, Benoît Valot, Didier Hocquet. Deciphering the role of insertion sequences in the evolution of bacterial epidemic pathogens with panISa software. Microbial Genomics, 2020, 6 (6), pp.e000356. 10.1099/mgen.0.000356. hal-02539297

HAL Id: hal-02539297

https://hal.science/hal-02539297v1

Submitted on 26 Sep 2024

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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_2.png" /></figure>

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

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

# Deciphering the role of insertion sequences in the evolution of bacterial epidemic pathogens with *panISa* software

Charlotte Couchoud¹,²*, Xavier Bertrand¹,², Benoit Valot²,³ and Didier Hocquet¹,²,³

## Abstract

Next-generation sequencing (NGS) is now widely used in microbiology to explore genome evolution and the structure of patho- gen outbreaks. Bioinformatics pipelines readily detect single-nucleotide polymorphisms or short indels. However, bacterial genomes also evolve through the action of small transposable elements called insertion sequences (ISs), which are difficult to detect due to their short length and multiple repetitions throughout the genome. We designed *panISa* software for the *ab initio* detection of IS insertions in the genomes of prokaryotes. *PanISa* has been released as open source software (GPL3) available from https://github.com/bvalot/panISa. In this study, we assessed the utility of this software for evolutionary studies, by reana- lysing five published datasets for outbreaks of human major pathogens in which ISs had not been specifically investigated. We reanalysed the raw data from each study, by aligning the reads against reference genomes and running *panISa* on the align- ments. Each hit was automatically curated and IS-related events were validated on the basis of nucleotide sequence similar- ity, by comparison with the ISFinder database. In *Acinetobacter baumannii*, the *panISa* pipeline identified *ISAba1* or *ISAba125* upstream from the *ampC* gene, which encodes a cephalosporinase in all third-generation cephalosporin-resistant isolates. In the genomes of *Vibrio cholerae* isolates, we found that early Haitian isolates had the same ISs as Nepalese isolates, confirm- ing the inferred history of the contamination of this island. In *Enterococcus faecalis*, *panISa* identified regions of high plasticity, including a pathogenicity island enriched in IS-related events. The overall distribution of ISs deduced with *panISa* was con- sistent with SNP-based phylogenetic trees, for all species considered. The role of ISs in pathogen evolution has probably been underestimated due to difficulties detecting these transposable elements. We show here that *panISa* is a useful addition to the bioinformatics toolbox for analyses of the evolution of bacterial genomes. *PanISa* will facilitate explorations of the functional impact of ISs and improve our understanding of prokaryote evolution.

## DATA SUMMARY

We confirm that all supporting data, code and protocols have been provided within the article or through Supplementary Material.

## INTRODUCTION

Whole-genome sequencing (WGS) is becoming the gold- standard technique for investigating the evolution of bacte- rial pathogen genomes during their spread. Application of the appropriate pipelines to sequencing data results in the

detection of single-nucleotide polymorphisms (SNPs) or small insertion/deletion (indels) after the alignment of reads with a reference genome sequence. However, bacterial genomes also evolve through the insertion of insertion sequences (ISs), which are widespread and occur in all domains of life [1]. ISs are mobile autonomous elements formed by (i) one or two transposase-encoding genes, (ii) two terminal inverted repeats (IRs), and (iii) two direct repeated sequences (DRs) [2]. ISs are sorted into families using the amino acid simi- larity of their transposase [3]. In 2019, the ISFinder database reported more than 4000 ISs belonging to 29 families [2, 3].

Received 28 August 2019; Accepted 02 March 2020; Published 26 March 2020

**Author affiliations:** ¹Laboratoire d'Hygiène Hospitalière, Centre Hospitalier Régional Universitaire, Besançon, France; ²UMR CNRS 6249 Chrono-environnement, Université de Bourgogne Franche-Comté, Besançon, France; ³Bioinformatique et big data au service de la santé, UFR Santé, Université de Bourgogne Franche-Comté, Besançon, France.

*Correspondence: Charlotte Couchoud, ccouchoud@chu-besancon.fr

**Keywords:** whole-genome sequencing; outbreak; insertion sequence; bacterial evolution.

**Abbreviations:** GC1, global clone 1; GFF, general feature format; ISs, insertion sequences; NGS, next-generation sequencing; PAI, pathogenicity island; SNPs, single-nucleotide polymorphisms; SRA, sequence read archive; WGS, Whole-genome sequencing.

Sequence Read Archive accession numbers for all re-analysed datasets are available in Supplementary Data S1.

**Data statement:** All supporting data, code and protocols have been provided within the article or through Supplementary Material. Supplementary Material is available with the online version of this article.

000356 © 2020 The Authors

<figure><img src="image_5.png" /></figure>
Genes are inactivated by the insertion of an IS into their coding sequences. ISs can also modulate the expression of a gene if they disrupt its promoter or create an alternative promoter [4]. Most known examples of IS transposition are linked to antibiotic resistance, because the resulting phenotypes are easy to detect. For example, the insertion of IS1 or IS10 upstream from the efflux pump gene *acrEF* increases the resistance of *Salmonella enterica* to fluoroquinolones [5]. Similarly, the insertion of IS*Aba1* or IS*Aba125* upstream from *ampC* increases resistance to third-generation cephalosporin in *Acinetobacter baumannii* [6, 7]. This ability of IS insertion to affect bacterial resistance to antibiotics and virulence can help bacterial pathogens to adapt to new niches [4]. IS dynamics are rarely investigated, but an understanding of these dynamics during outbreaks of pathogenic prokaryotes could be highly informative in evolutionary studies.

The detection of ISs is challenging, because read lengths are usually shorter (<300 bp) than ISs and the same IS may be repeated in the genome. We have developed the *panISa* program to detect new and unknown insertions *ab initio* (i.e.

with a database-free approach) in bacterial genomes, based on the detection of structural variants in short-read data. *PanISa* requires only short reads and a reference genome as input. The presence of ISs and their repeated nature renders the IS localization from the assembly very challenging, which is why we used a software that localizes ISs from the raw reads. The program has been validated on simulated data and compared with existing tools, but the benefits of IS detection for epidemiological studies remain to be evaluated [8].

We therefore assessed the dynamics of IS insertion during the spread of bacterial pathogens, by reanalysing the WGS data from five published studies describing genome evolution for major epidemic bacterial pathogens.

# METHODS

## Selection of the datasets

We reanalysed five published datasets from studies aiming to decipher the evolution of a specific clone of a bacterial pathogen during its spread (Table 1). All the species concerned were major human pathogens for which ISs have already been described.

## Impact Statement

Insertion sequences (ISs) are small transposable elements playing a key role in bacterial genome organization and evolution. They are difficult to detect in sequencing data. We therefore designed *panISa* software for the *ab initio* detection of IS insertion in prokaryotic genomes. Here, to evaluate the potential of this new tool for use in evolutionary studies, we selected five published studies describing genome evolution in five major human epidemic pathogens. None of these studies had used bioinformatics pipelines to retrieve ISs. The ISs retrieved by *panISa* had a genomic distribution consistent with SNP-based phylogenetic analysis. Our pipeline rapidly detected IS-related mechanisms of resistance to antibiotics and identified genomic regions of high plasticity with a high concentration of IS-related events. The proportion of genomes displaying IS insertions varied considerably between bacterial species, but was at least 86% for *Acinetobacter baumannii*, *Vibrio cholerae*, and *Enterococcus faecalis*. *PanISa* is a useful addition to the bioinformatics toolbox for analyses of prokaryote evolution. It will help us to determine the role of ISs in pathogen evolution, which has probably been underestimated.

The first dataset was from a study describing the global epidemiology of the seventh cholera pandemic that aimed to identify the geographic origin of the contamination of Haiti [9]. We also selected a dataset for isolates of the transcontinental epidemic strain of *Vibrio parahaemolyticus* ST36 [10]. The authors reconstructed the evolution of this clone over the last 25 years, by genome-wide analysis. We selected a dataset for genomes from a global collection of isolates of *Enterococcus faecalis* retrieved between 1958 and 2012 from bloodstream infections in Ireland, the UK and the USA [11]. Another collection of genomes of *Salmonella enterica* subspecies *enterica* serotype Tennessee (S. Tennessee) isolates from an outbreak involving transmission in peanut butter was also selected [12]. In this study, SNP-based analysis revealed

<table><caption>**Table 1.** Collections of genomes of bacterial pathogens reanalysed with *panISa* software</caption><thead><tr><th rowspan="2">Reference</th><th rowspan="2">Reference genome (NCBI accession numbers)</th><th colspan="2">Isolates (n)</th></tr><tr><th>Total no.</th><th>Isolates with available data</th></tr></thead><tbody><tr><td>Eppinger et al. [9]</td><td>Vibrio cholerae O1 biovar El tor str. N16961 (AE003852.1; AE003583.1)</td><td>116</td><td>110</td></tr><tr><td>Martinez-Urtaza et al. [10]</td><td>Vibrio parahaemolyticus RIMD 2210633 (BA000031.2; BA000032.2)</td><td>48</td><td>16</td></tr><tr><td>Raven et al. [11]</td><td>Enterococcus faecalis V583 (AE016830.1)</td><td>168</td><td>168</td></tr><tr><td>Wilson et al. [12]</td><td>Salmonella Tennessee str. TXSC_TXSC08-19 (CP007505.1)</td><td>69</td><td>68</td></tr><tr><td>Holt et al. [13]</td><td>Acinetobacter baumannii strain A1 (CP010781.1)</td><td>44</td><td>35</td></tr></tbody></table>
that the contamination was of environmental origin. Finally, we selected a dataset of genomes of isolates of *Acinetobacter baumannii* global clone 1 (GC1) collected between 1960 and 2011, the analysis of which provided insight into the evolution of these genomes and the phylodynamics of GC1 [13].

## Downloading and pretreatment of short-read data

For the five datasets, we first determined the number of isolates for which short-read data were available in the sequence read archive (SRA) database (1. Eppinger et al. [9], 2. Martinez-Urtaza et al. [10], 3. Raven et al. [11], 4. Wilson et al. [12], 5. Holt et al. [13]) (Table 1; Data S1, available in the online version of this article). We used prefetch and fastq-dump command-line tools from NCBI toolkits to extract and convert the data into an optimized input format for *panISa* [14]. Reads were subsampled to a final coverage of 60× and mapped against the same reference genome as in the original studies with the Burrows-Wheeler aligner [15].

## PanISa search

*PanISa* identifies IS insertions through comparison with the pysam library. Briefly, each read that maps partially on the reference genome was detected by *panISa* as a clipped read. When clipped reads were in opposite directions on two close positions of the genome, the program identified a potential IS insertion (also called a hit). The flanking parts of the clipped reads defined the boundaries of the insertion, with start clipped reads defining the IRR and end clipped reads defining the IRL (see *Treepong* et al. [8] for a detailed description of the functioning of the software). *PanISa* has been released as open-source software (GPL3) available from https://github.com/bvalot/panISa. Each alignment file (.bam) was used as input for *panISa* with the minimum clipped reads option set to 10 (default settings for all other options). As *panISa* detects all insertion events, manual curation of the list of potential ISs is required. We checked for sequence similarity (identity >90% over >80% of the length of the sequence) between the boundaries of the IS (IRR and IRL), with the ISFinder database, to confirm IS insertions [16]. To clarify the terms used thereafter, *panISa* produces a 'hit' for each detection of whichever insertion, while an 'IS-related event' describes the

insertion of an IS. In other words, if one specific IS inserted at two different positions, we described it as one IS and two IS-related events.

## Estimation of the biological impact of IS insertion

We downloaded the annotations (i.e. GFF general feature format) of each reference genome to identify the function of the proteins encoded by the genes disrupted by ISs. In addition, for ISs inserted into intergenic regions, we identified the function of the proteins encoded by genes with translation start sites <100 bp away from the insertion site of an IS.

## RESULTS AND DISCUSSION

In this study, we used *panISa* to reanalyse published WGS datasets. Unsurprisingly, none of these published studies could identify new ISs or new IS insertion sites. However, Holt *et al.* [13] reported the detection of *ISAba1* and *ISAba125* upstream from the cephalosporinase-encoding *ampC* gene by a PCR approach for the identification of genetic events leading to antibiotic resistance [13]. ISs can contribute to antibiotic resistance and virulence, depending on the nature of the genes that are disrupted or imported. They can also, more widely, contribute to bacterial genome rearrangements. The identification of IS-related events with an accurate tool would therefore improve our understanding of genome evolution in epidemic pathogens.

## IS detection in the genomes of epidemic bacterial pathogens

*PanISa* retrieved 692 to 15878 hits from the five datasets, the smallest number of hits being obtained for the study of *V. parahaemolyticus* genomes and the largest number for the study of *E. faecalis* (Table 2). Only hits matching sequences in the ISFinder database were considered to be IS-related events. Thus, 207 of the 2913 hits were identified as IS-related events in the *V. cholerae* dataset [9]. *PanISa* identified 1348 IS-related events among the 15878 hits for the *E. faecalis* dataset and 345 IS-related events among the 1371 hits for the *A. baumannii* dataset [11, 13]. None of the 692 hits for the *Vibrio parahaemolyticus* dataset was associated with an IS

<table><caption>**Table 2.** Result of the reanalysis of five genome collections with *panISa*.</caption><thead><tr><th>Reference</th><th>Species</th><th>panISa hits (n)</th><th>IS-related events (n)</th><th>ISs (n)</th><th>Proportion (%) of isolates with ≥1 IS</th></tr></thead><tbody><tr><td>Eppinger et al. [9]</td><td>Vibrio cholerae</td><td>2913</td><td>207</td><td>5</td><td>91</td></tr><tr><td>Martinez-Urtaza et al. [10]</td><td>Vibrio parahaemolyticus</td><td>692</td><td>0</td><td>0</td><td>0</td></tr><tr><td>Raven et al. [11]</td><td>Enterococcus faecalis</td><td>15878</td><td>1348</td><td>29</td><td>100</td></tr><tr><td>Wilson et al. [12]</td><td>Salmonella Tennessee</td><td>727</td><td>4</td><td>1</td><td>1.4</td></tr><tr><td>Holt et al. [13]</td><td>Acinetobacter baumannii</td><td>1371</td><td>345</td><td>19</td><td>86</td></tr></tbody></table>

The third column shows the number of insertion events identified by *panISa* from WGS datasets of bacterial pathogens, the fourth column gives the number of hits matching sequences in the ISFinder database, considered to correspond to IS-related events, the fifth column gives the number of different ISs found among the IS-related events, and the last column gives the proportion of isolates concerned.
[10]. We identified four IS-related events among the 727 hits for S. Tennessee, all involving the same IS in a single isolate [12]. The small numbers of IS-related events detected in these last two datasets may be due to either the short length of the reads (<100 bp) creating artefactual repeat sequences or the low coverage (<40×) for several genomes [8]. Furthermore, as Wilson *et al.* [12] studied genomes of S. Tennessee isolates collected between 2006 and 2008, the study period may have been too short for IS insertion to happen [12] (Data S2). In addition, the short length of the reads or the presence of repeated regions may account for the high level of background noise.

## Phylogeny of the Vibrio cholerae outbreak

Haiti had been free from *V. cholerae* for 100 years before being hit, in 2010, by an outbreak that lasted 2 years. Phylogenetic analysis revealed that Haiti had been contaminated by a Nepalese strain brought to the island by United Nations peacekeepers [17]. We reanalysed datasets of genome sequences for *V. cholerae* isolates from Haiti, Nepal and Bangladesh

[9]. Phylogenetic analysis of these sequences performed by Eppinger *et al.* [9] revealed these isolates to be highly homogeneous. *PanISa* identified 207 insertions of five different ISs. The IS256-like element, inserted at position 2440662, was common to all isolates of the dataset. The IS256-like element in position 327875 was specific to the genomes of isolates from cluster 'a', which grouped together isolates obtained in Haiti in 2010 and isolates from Nepal (Fig. 1a). It is unlikely that independent insertions of the same IS occurred at the same genomic position, so the insertions of the IS256-like element probably occurred in a common ancestor of the Haitian isolates and the Nepalese isolates of the 'Nepal-4' subcluster (Fig. 1a). This IS insertion pattern is consistent with the close relationship showed by the phylogeny between the Haitian and Nepalese isolates.

Similarly, the IS1634-like element inserted at position 91094 was specific to isolates clustering in the 'Bangladesh-1' group, and the three IS1634-like elements inserted at positions 452943–878317 and 911076 were specific to isolates

<figure><img src="image_6.png" /><figcaption>**Fig. 1.** Comparison of (a) phylogenetic analysis from Eppinger *et al.* [9] and (b) IS-related events identified by *panISa*. Red boxes represent insertions of IS1634-like elements and green boxes represent insertions of IS256-like elements. Each column represents an insertion site on a specific chromosome. The two colours refer to the two different ISs detected in the dataset, and the different shades of colours represent the different positions of insertion of each IS [at the chromosomal positions indicated at the bottom of (b)].</figcaption></figure>
clustering in the 'Nepal-1' group, as defined by Eppinger *et al.* [9]. These insertions undoubtedly occurred in a common ancestor and were then transmitted to all isolates derived from that ancestor (Fig. 1b).

The localization, by *panISa*, of the IS insertions in the genomes made it possible to identify the genes affected and to predict the biological consequences of the gene disruption or the modification of gene expression. Unfortunately, the disrupted genes encoded proteins did not seem to be involved in known virulence or resistance (Data S3), making it difficult to implicate any of these changes in the development of antibiotic resistance or virulence during the spread of *V. cholerae*. However, the identification of IS insertions with

*panISa* provided information consolidating the phylogenetic trees built from SNP data.

## IS insertions in the antibiotic-resistant *Acinetobacter baumannii* global clone 1

We reanalysed the *A. baumannii* genome dataset with *panISa* and retrieved 345 IS-related events in 35 genomes, with 86% of the genomes displaying at least one IS-related event [13] (Table 2; Data S4). The wide temporal coverage of the collection (~50 years) probably accounts for the large number of IS events retrieved (Fig. 2a). Likewise, the isolates were selected to represent the maximum diversity among the global clone 1, the high variability of sources and locations

<figure><img src="image_7.png" /><figcaption>**Fig. 2.** Phylogenetic analysis of *A. baumannii* GC1. (a) Temporal and phylogenetic analysis performed by Holt *et al.* [13]. The dot colours represent the capsule type of each isolate. (b) Hierarchical clustering based on the presence/absence of 345 IS-related events implicating 19 ISs in 35 isolates of *A. baumannii* with available SRA data. Each of the 345 columns represents one IS-related event, with grey indicating no event, black indicating an event and white if the sequence data could not be obtained. (c) Representation of the insertion of *ISAba1* and *ISAba125* upstream from *ampC*.</figcaption></figure>
presumably accounted for the higher number of IS-related events detected. We performed hierarchical clustering based on the presence/absence of the 345 IS-related events, which revealed similarities between the isolates (Fig. 2b). Similar IS-related event profiles were observed in some clusters of two or three isolates, such as Canada-BC1 and Canada-BC5, and isolates D78 and D81. The global results for IS-based clustering were consistent with those based on SNP data reported by Holt *et al.* [13] (Fig. 2a).

The original study aimed to identify the determinants of resistance to third-generation cephalosporins. The insertion of an IS upstream from the cephalosporinase-encoding *ampC* gene can lead to the overproduction of AmpC, increasing resistance to third-generation cephalosporins [6, 7]. The authors of the original study used PCR and sequencing to search for ISs upstream from *ampC*, but *panISa* rapidly and correctly identified all the *ISAba1* and *ISAba125* insertions upstream from *ampC* in all 13 AmpC-overproducing isolates (Fig. 2c) [13]. Thirteen of the 19 different ISs retrieved by *panISa* had already been described in this pathogen [16]. The proportion of isolates displaying at least one IS-related event was consistent with previous reports [18]. The most active IS in the global clone was *ISAba1*, with 7% of its insertions disrupting the genome upstream from *ampC* [18]. We found that the most common position for insertions in these genomes was upstream from an IS256 (after nucleotide 287 817) element already present in the reference genome.

Thus, *panISa* is a useful tool for identifying IS insertions with a well-described biological impact, such as antibiotic resistance. *PanISa* can not only retrieve IS insertions in well-known location (i.e. upstream from *ampC*) but also everywhere in an Illumina-sequenced genome. *PanISa* can accelerate the laborious task of targeted searches for ISs in genomes. Our genomic data also indicate that this pipeline can detect unexpected insertion sites, thereby improving our understanding of the genomic events leading to particular phenotypes.

## Distribution of IS insertions in the genome of *Enterococcus faecalis*

We reanalysed the data for an international collection of isolates of *E. faecalis* from bloodstream infections collected over the last 50 years [11]. *PanISa* retrieved 1348 IS-related events at 472 different genomic sites. All of the 168 genomes for which SRA data were available presented at least one IS insertion (Table 2). More than a third of the IS insertions involved *ISEnfa4*, an *E. faecalis* IS from the IS256 family. The large number of ISs identified presumably reflects the long time period covered by the collection, the large phylogenetic distance between isolates and the large number of clones collected. Most (73%; 363 of 496) of the sites of IS-related events were unique to a single isolate (Fig. 3a). We investigated the 28 IS-related events common to at least ten isolates to identify events that had been selected through evolution. Six of these events occurred more than 100 bp away from a translation start site (Data S5). Eight of the 23 IS-related events potentially affecting gene function disrupted a gene,

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

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

**Fig. 3.** Distribution of the IS-related events in the genomes of a collection of 168 clinical isolates of *E. faecalis* [11]. (a) Distribution of the number of IS-related events as a function of the number of isolates affected. The y-axis is drawn to a logarithmic scale. Red bars represent IS-related events that occurred >100 bp from a coding sequence or occurred in less than ten isolates. Black bars represent IS-related events that occurred within 100 bp of a coding sequence and are reported in Table 3. (b) The outer green circle represents the number of IS-related events, the red central circle indicates the number of genomes from the collection of 168 isolates of *E. faecalis* containing the region, and the inner grey circle represents GC-content relative to the mean value. The 150 kb pathogenicity island and the ICE (ICEEEfaV583-1) of *E. faecalis* V583 are indicated by the green and red sectors, respectively.

ten occurred close (<100 bp) to the translation start site of a single gene, and five events occurred close to the transcription start sites of two genes (Table 3). One third of the IS events in our analysis of *E. faecalis*, disrupted genes but most insertions (20 of 28) occurred in intergenic regions, as already reported in *Shigella flexneri* strain 2457T [19]. We therefore tried to determine the functions of the proteins encoded by the genes with promoters or coding sequences disrupted by an IS. A third of the IS-disrupted genes encoded proteins of unknown function (Table 3). None of the IS insertions identified affected known antibiotic resistance or virulence genes.
<table><caption>**Table 3.** IS-related events common to at least ten of the 168 *E. faecalis* isolates from the Raven *et al.* dataset [11]</caption><thead><tr><th>Isolates (n)</th><th>IS</th><th>Position</th><th>Protein potentially affected by the IS insertion</th><th>GeneID</th></tr><tr><th></th><th></th><th>Genomic</th><th>In relation to the closest gene</th><th>Function</th><th></th></tr></thead><tbody><tr><td>10</td><td>ISEnf4</td><td>499863</td><td>In</td><td>Hypothetical protein</td><td>gene540</td></tr><tr><td rowspan="2">10</td><td rowspan="2">ISEfa10</td><td rowspan="2">2594127</td><td>Upstream</td><td>HAD superfamily hydrolase</td><td>gene2620</td></tr><tr><td>Downstream</td><td>Hypothetical protein</td><td>gene2621</td></tr><tr><td>10</td><td>IS6770</td><td>1627490</td><td>In</td><td>ABC transporter ATP-binding protein</td><td>gene1652</td></tr><tr><td>10</td><td>ISEnf4</td><td>705382</td><td>Downstream</td><td>cell wall surface anchor family protein</td><td>gene733</td></tr><tr><td>11</td><td>ISEnf4</td><td>491272</td><td>Upstream</td><td>Hypothetical protein</td><td>gene530</td></tr><tr><td>11</td><td>IS1062</td><td>991057</td><td>In</td><td>Phosphorylase</td><td>gene1006</td></tr><tr><td>11</td><td>ISEnf4</td><td>2443962</td><td>Downstream</td><td>Conjugal transfer protein</td><td>gene2462</td></tr><tr><td rowspan="2">12</td><td rowspan="2">ISChh1</td><td rowspan="2">1337504</td><td>Upstream</td><td>Hypothetical protein</td><td>gene1348</td></tr><tr><td>Downstream</td><td>hydroxymethylglutaryl-CoA synthase</td><td>gene1349</td></tr><tr><td>13</td><td>ISEnf3</td><td>1224046</td><td>Downstream</td><td>ABC transporter ATP-binding protein</td><td>gene1240</td></tr><tr><td>14</td><td>ISEfa10</td><td>218868</td><td>Upstream</td><td>Hypothetical protein</td><td>gene222</td></tr><tr><td>14</td><td>IS6770</td><td>832847</td><td>Upstream</td><td>Potassium uptake protein</td><td>gene850</td></tr><tr><td>19</td><td>ISEfa10</td><td>2650256</td><td>Upstream</td><td>Lipoate-protein ligase A</td><td>gene2676</td></tr><tr><td rowspan="2">21</td><td rowspan="2">ISEfa10</td><td rowspan="2">2594120</td><td>Downstream</td><td>HAD superfamily hydrolase</td><td>gene2620</td></tr><tr><td>Upstream</td><td>Hypothetical protein</td><td>gene2621</td></tr><tr><td rowspan="2">21</td><td rowspan="2">ISEfa10</td><td rowspan="2">2809601</td><td>Downstream</td><td>Hypothetical protein</td><td>gene2863</td></tr><tr><td>Upstream</td><td>Valyl-tRNA synthetase</td><td>gene2864</td></tr><tr><td>29</td><td>ISApl3</td><td>352544</td><td>In</td><td>DadA family oxidoreductase</td><td>gene418</td></tr><tr><td>31</td><td>ISEfa10</td><td>2404627</td><td>Upstream</td><td>Hypothetical protein</td><td>gene2430</td></tr><tr><td>32</td><td>ISEfm1</td><td>1317331</td><td>In</td><td>Hypothetical protein</td><td>gene1332</td></tr><tr><td rowspan="2">34</td><td rowspan="2">ISEfa10</td><td rowspan="2">2594130</td><td>Downstream</td><td>HAD superfamily hydrolase</td><td>gene2620</td></tr><tr><td>Upstream</td><td>Hypothetical protein</td><td>gene2621</td></tr><tr><td>37</td><td>ISEnf4</td><td>1954670</td><td>In</td><td>Hypothetical protein</td><td>gene1987</td></tr><tr><td>39</td><td>IS1485</td><td>1805579</td><td>In</td><td>3-methyl-2-oxobutanoate<br>hydroxymethyltransferase</td><td>gene1826</td></tr><tr><td>43</td><td>ISEfa4</td><td>608941</td><td>In</td><td>DeoR family transcriptional regulator</td><td>gene644</td></tr><tr><td>43</td><td>IS6770</td><td>641514</td><td>Upstream</td><td>Rotamase</td><td>gene672</td></tr></tbody></table>

The second and third columns give the name of the IS and the position of its insertion in the reference genome *E. faecalis* V583. The fourth column gives the position of the IS insertion in relation to the gene potentially affected, and the fifth column gives the function of the protein and the gene ID of the gene potentially affected by the IS. More detailed information are given in Data S5.

We then explored the distribution of IS-related events in the genomes of *E. faecalis*. Surprisingly, 21% of these events (281 of 1348) were clustered together in a 150 kb region in which the IS insertion frequency (1.8 insertions/1000 bp) was four times that elsewhere in the genome (0.4 insertions/1000 bp) (Fig. 3b, green sector). We also analysed GC-content in the genome and found that this region had a lower GC-content that the rest of the genome, indicating recent horizontal gene transfer [19]. This region is a pathogenicity island (PAI) encompassing genes encoding the cytolysin toxin, the

enterococcal surface protein Esp, Glس-24-like proteins, and proteins of unknown function [20–23] (Fig. 3b). Overall, 74% of the genomes studied (125 of 168) presented at least one IS-related event in this 150 kb region, despite the uneven coverage of this region between genomes. The accelerated genetic drift of this PAI might suggest a neutral or positive biological impact of gene disruption but needs further exploration. The gene content of this PAI is known to be highly variable, but the effect of IS insertion on the evolution of this region has yet to be explored [24]. However, ISs have
been implicated in the diversification of *E. faecium* [25]. Conversely, we also identified a region (between positions 2204066 and 2258320 of the reference genome *E. faecalis* V583) with a higher GC-content (Fig. 3b, red sector). This region corresponds to an integrative conjugative element (ICEEfaV583-1) encompassing the vancomycin resistance cassette *vanB*. The detection of IS in that particular region was impossible since absent in all genomes in this dataset but one.

## Limitations and benefits of *panISa*

The choice of reference genome impacts the number of ISs detected, with the background noise increasing with the distance between the studied sequences and the reference genome. Moreover, the detection of IS in a genomic region absent from the reference genome is impossible with *panISa*. *PanISa* requires manual curation of the annotation after finding/validating the IS in the genomes. As the *panISa* pipeline runs with raw reads as an input, it can expedite the reanalysis of sequenced data, avoiding the step of assembly. The validation of the IS-related events is based on the reconstructed boundaries of the IS (i.e. the IR), therefore IS fragments are also detected and reported. *PanISa* is easy to install and requires few dependencies, is lightweight and can be run on a laptop, making server implementation unnecessary.

The dynamics of IS insertion within bacterial genomes remains incompletely understood, but several studies have shown that the sites of transposable element insertion are not randomly distributed between coding and non-coding regions [19, 26]. We show here that IS insertions are not randomly distributed throughout the genome. The detection of these events with *panISa* software will shed light on the dynamics of IS insertion.

## Conclusion

*PanISa* is a software pipeline for detecting IS insertions in prokaryotic genomes from short-read data. It expands the toolkit available for exploring the evolution of prokaryotic lineages. ISs are difficult to detect in short-read sequencing data, and this has probably resulted in an underestimation of the impact of these mobile elements on the evolution of their bacterial hosts. Studies of IS dynamics in bacterial genomes have also been hindered by difficulties in genome assembly.

*PanISa* can increase our understanding of the evolution of bacterial pathogens during their spread. For example, we show here that *panISa* can consolidate phylogenetic analyses of large datasets. It can also accelerate the identification of IS events with a known biological impact, such as those triggering resistance to third-generation cephalosporins in *A. baumannii*. *PanISa* can also be used to identify new functional impacts of ISs during the spread of pathogens. Determinations of the pattern or frequency of IS insertion over the genome will undoubtedly help us to decipher the evolution of bacterial lineages and the dynamics of IS insertions.

## Funding information

This work received no specific grant from any funding agency.

## Conflicts of interest

The authors have no conflicts of interest to declare.

## Data Bibliography

1. Eppinger M, Pearson T, Koenig SSK, Pearson O, Hicks N, et al. Genomic epidemiology of the Haitian cholera outbreak: a single introduction followed by rapid, extensive, and continued spread characterized the onset of the epidemic. *mBio* 2014;5:e01721-14.

2. Martinez-Urtaza J, van Aerle R, Abanto M, Haendiges J, Myers RA, et al. Genomic variation and evolution of *Vibrio parahaemolyticus* ST36 over the course of a transcontinental epidemic expansion. *mBio* 2017;8:e01425-17.

3. Raven KE, Reuter S, Gouliouris T, Reynolds R, Russell JE, et al. Genome-based characterization of hospital-adapted *Enterococcus faecalis* lineages. *Nature Microbiology* 2016;1:15033.

4. Wilson MR, Brown E, Keys C, Strain E, Luo Y, et al. Whole genome DNA sequence analysis of *Salmonella* subspecies enterica serotype Tennessee obtained from related peanut butter foodborne outbreaks. *PLOS ONE* 2016;11:e0146929.

5. Holt K, Kenyon JJ, Hamidian M, Schultz MB, Pickard DJ, et al. Five decades of genome evolution in the globally distributed, extensively antibiotic-resistant *Acinetobacter baumannii* global clone 1. *Microbial Genomics*;2. Epub ahead of print 9 February 2016. DOI: 10.1099/mgen.0.000052.

## References

1. **Siguier P, Gourbeyre E, Chandler M.** Bacterial insertion sequences: their genomic impact and diversity. *FEMS Microbiol Rev* 2014;38:865-891.

2. **Siguier P, Gourbeyre E, Varani A, Ton-Hoang B, Chandler M.** Everyman's guide to bacterial insertion sequences. *Microbiol Spectr* 2015;3:MDNA3-0030-2014.

3. **Mahillon J, Chandler M.** Insertion sequences. *Microbiol Mol Biol Rev* 1998;62:725-774.

4. **Vandecraen J, Chandler M, Aertsen A, Van Houdt R, Houdt RV.** The impact of insertion sequences on bacterial genome plasticity and adaptability. *Crit Rev Microbiol* 2017;43:709-730.

5. **Olliver A, Vallé M, Chaslus-Dancla E, Clockaert A.** Overexpression of the multidrug efflux operon *acrEF* by insertional activation with IS1 or IS10 elements in *Salmonella enterica* serovar typhimurium DT204 *acrB* mutants selected with fluoroquinolones. *Antimicrob Agents Chemother* 2005;49:289-301.

6. **Segal H, Nelson EC, Elisha BG.** Genetic environment and transcription of *ampC* in an *Acinetobacter baumannii* clinical isolate. *Antimicrob Agents Chemother* 2004;48:612-614.

7. **Hamidian M, Hall RM.** ISAbal targets a specific position upstream of the intrinsic *ampC* gene of *Acinetobacter baumannii* leading to cephalosporin resistance. *J Antimicrob Chemother* 2013;68:2682-2683.

8. **Treepong P, Guyeux C, Meunier A, Couchoud C, Hocquet D** et al. *panISa: ab initio* detection of insertion sequences in bacterial genomes from short read sequence data. *Bioinformatics* 2018;34:3795-3800.

9. **Eppinger M, Pearson T, Koenig SSK, Pearson O, Hicks N** et al. Genomic epidemiology of the Haitian cholera outbreak: a single introduction followed by rapid, extensive, and continued spread characterized the onset of the epidemic. *mBio* 2014;5:e01721-14.

10. **Martinez-Urtaza J, van Aerle R, Abanto M, Haendiges J, Myers RA** et al. Genomic variation and evolution of *Vibrio parahaemolyticus* ST36 over the course of a transcontinental epidemic expansion. *mBio* 2017;8:e01425-17.

11. **Raven KE, Reuter S, Gouliouris T, Reynolds R, Russell JE** et al. Genome-based characterization of hospital-adapted *Enterococcus faecalis* lineages. *Nature Microbiology* 2016;1:15033.

12. **Wilson MR, Brown E, Keys C, Strain E, Luo Y** et al. Whole genome DNA sequence analysis of *Salmonella* subspecies *enterica*
serotype Tennessee obtained from related peanut butter food-
borne outbreaks. *PLoS One* 2016;11:e0146929.

13. Holt K, Kenyon JJ, Hamidian M, Schultz MB, Pickard DJ et al. Five decades of genome evolution in the globally distributed, extensively antibiotic-resistant *Acinetobacter baumannii* global clone 1. *Microb Genom* 2016;2:e000052.

14. . Download: Software: Sequence Read Archive: NCBI/NLM/NIH. Download: Software: Sequence Read Archive: NCBI/NLM/NIH.380 https://trace.ncbi.nlm.nih.gov/Traces/sra/sra.cgi?view=software.

15. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. *Bioinformatics* 2009;25:1754-1760.

16. Siguier P, Perochon J, Lestrade L, Mahillon J, Chandler M. ISfinder: the reference centre for bacterial insertion sequences. *Nucleic Acids Res* 2006;34:D32–D36.

17. Hendriksen RS, Price LB, Schupp JM, Gillece JD, Kaas RS et al. Population genetics of *Vibrio cholerae* from Nepal in 2010: evidence on the origin of the Haitian outbreak. *mBio* 2011;2:e00157-11.

18. Adams MD, Bishop B, Wright MS. Quantitative assessment of insertion sequence impact on bacterial genome architecture. *Microb Genom* 2016;2:e000062.

19. Zaghloul L, Tang C, Chin HY, Bek EJ, Lan R et al. The distribution of insertion sequences in the genome of *Shigella flexneri* strain 2457T. *FEMS Microbiol Lett* 2007;277:197-204.

20. Shankar N, Baghdayan AS, Gilmore MS. Modulation of virulence within a pathogenicity island in vancomycin-resistant *Enterococcus faecalis*. *Nature* 2002;417:746–750.

21. Chow JW, Thal LA, Perri MB, Vazquez JA, Donabedian SM et al. Plasmid-associated hemolysin and aggregation substance production contribute to virulence in experimental enterococcal endocarditis. *Antimicrob Agents Chemother* 1993;37:2474–2477.

22. Teng F, Nannini EC, Murray BE. Importance of gls24 in virulence and stress response of *Enterococcus faecalis* and use of the Gls24 protein as a possible immunotherapy target. *J Infect Dis* 2005;191:472–480.

23. Shankar V, Baghdayan AS, Huycke MM, Lindahl G, Gilmore MS. Infection-derived *Enterococcus faecalis* strains are enriched in *esp*, a gene encoding a novel surface protein. *Infect Immun* 1999;67:193–200.

24. McBride SM, Coburn PS, Baghdayan AS, Willems RJL, Grande MJ et al. Genetic variation and evolution of the pathogenicity island of *Enterococcus faecalis*. *J Bacteriol* 2009;191:3392–3402.

25. Leavis HL, Willems RJL, van Wamel WJB, Schuren FH, Caspers MPM et al. Insertion sequence-driven diversification creates a globally dispersed emerging multiresistant subspecies of *E. faecium*. *PLoS Pathog* 2007;3:e7.

26. Plague GR. Intergenic transposable elements are not randomly distributed in bacteria. *Genome Biol Evol* 2010;2:584–590.

Five reasons to publish your next article with a Microbiology Society journal

1. The Microbiology Society is a not-for-profit organization.

2. We offer fast and rigorous peer review – average time to first decision is 4–6 weeks.

3. Our journals have a global readership with subscriptions held in research institutions around the world.

4. 80% of our authors rate our submission process as 'excellent' or 'very good'.

5. Your article will be published on an interactive journal platform with advanced metrics.

Find out more and submit your article at microbiologyresearch.org.