IS Families/IS481 family

Historical

IS481 was first identified in the late 1980s as a repeated sequence in the agent of whooping cough, Bordetella pertussis. It was estimated, from Southern blots of restriction digested B. pertusis genomic DNA, that IS481 was present in multiple copies of about 20, with some variation between strains. They were not observed as repeated sequences in the single strain of B. parapertussis or in B. bronchiseptica which were analyzed, although homology was observed as a single band in B. parapertussis strain BPAH1 genomic digests and two bands in digests of B. bronchiseptica strain BBRH1.

The DNA sequence later identified IS481 as a 1052bp element with 28bp terminal inverted repeats and confirmed that it was not present in multiple copies in either theB. parapertussis or B. bronchiseptica genomes. This was later reassessed to be 1046bp in length (Fig. IS481.1; ) and to massively hybridize with pertussis genomic DNA of 100 B. pertusis clinical isolates. There were a number of small sequence variations observed in the results from different authors.

It was also observed to transpose. When the bvgA gene, a global virulence regulator, was present as both a plasmid and chromosomal copy, host growth was inhibited when the plasmid-borne gene carries a small deletion . Those derivatives in which growth inhibition was relieved had acquired IS481 insertions into the plasmid-borne gene or IS481/IS1002 insertions into the chromosomal gene and to generate 6bp target DR .

Fig. IS481.1. Organization of IS481. Top: A map of IS481. The insertion sequence is shown as a yellow horizontal box with its transposase gene as a lilac horizontal arrow to indicate the direction of expression. The 28 base pair inverted repeated are shown as horizontal gray arrows. Bottom: Transcriptional organization at the left IS end. The DNA sequence including the left inverted repeat (gray shading) together with the transposase promoter (Pin) oriented to the right with its -35 and -10 components labelled, the transcription start shown by a red arrow and the transposase initiation codon, ATG, shown in red. An outward directed promoter Pout, overlapping with Pin and oriented to the left with its -35 and -10 components labelled, the transcription start shown by a blue arrow is also included .

Relationship to other IS.

Initially, IS481 appeared to be an IS3 family derivative which had been truncated for the N-terminal end of the Tpase and includes a C-terminal extension. The DDE active site domain and the IR (ending in 5’ TGT 3’) are similar to those of IS3 family members. Although, at that time IS481 had not been demonstrated to transpose, their presence in high copy number in some species and the identification of at least 130 distinct but related IS from over 90 species strongly suggested that these represent a distinct transpositionally active family. Different members generate DR of between 4 and 15 bp. Moreover, certain members (e.g. ISSav7) insert specifically into the tetranucleotide CTAG which becomes the flanking DR and provides the UAG termination codon for the Tpase. In contrast to the vast majority of IS3 family members, the IS481 Tpase is not produced by frameshifting. There is no evidence for a leucine zipper as in IS3. This can be clearly seen from the alignment and remarkable predicted structural similarities in the C-terminal catalytic domain and small N-terminal domain (Fig. IS481.2).

They also show a conserved 5’TG 3’ tip to the IR which is typical of this and a number of other types of prokaryotic mobile elements.

IS481 is also distantly related to the IS1202 family .

Fig. IS481.2 Comparison of the IS481 transposase with that of the IS3 family IS911. Predicted structures were obtained from Alphafold and compared using FATCAT tool .

Transporter (t)IS derivatives

Some members carry passenger genes (Fig. IS481.3) including antibiotic resistance (Cm^R, cmx, for IS5564 and the closely related ISCgl1), or potential transcriptional regulators among other open reading frames (ISKrh1, ISPfr21- which contains a termination codon within the transposase, and ISSav7).

Fig. IS481.3. IS481 tIS Organization. Examples of transporter (t) IS481 family members are shown. The tIS are represented by horizontal yellow boxes and their terminal inverted repeats are shown as gray horizontal arrows. Transposase genes are indicated by horizontal lilac arrows indicating their direction of expression. The chloramphenicol antibiotic resistance gene, cmx, is shown in red, while other genes are shown in peach color. Note that the function of most of these genes has not been experimentally determined.

Use of IS481 for Diagnostic Purposes

There are PCR diagnostic clinical protocols for Bordetella based on IS481 although several instances of false positive results have been reported.

Relationship to Eukaryotic IS

IS481 family IS are distantly related to the eukaryotic Banshee transposon which at present is restricted to the anaerobic flagellated protozoan Trichomonas vaginalis . They share the highly conserved Pfam integrase core domain identified in the IS3 family and retroviruses . Kojima and Bao , have identified a number of related IS from Trichomonas vaginalis which they have called IS481EU.

Interestingly, these were not well separated from a group of IS3 family-related eukaryotic sequences they called IS3EU/GingerRoot (similarly to the prokaryotic relatives) (Fig. IS481.4). It would be interesting to determine whether IS481EU transpose using a dsDNA circular intermediate as do IS3 family members.

Fig.IS481.4 Phylogenetic tree based on the transposase domains (DDE) of related prokaryotic insertion sequences compared to examples of different eukaryotic transposons. IS were those in the June 2022 version of ISfinder. Lineages are those with over 50% bootstrap support.

IS481 in the Evolution of Bordetella

IS481 itself has played a fundamental role in the evolution of the genomes of the Bordetellae. The genus Bordetella includes 16 named species (see Wigland et al ). Three of these, B. bronchiseptica, B. parapertussis and B. pertussis are closely related and colonize the mammalian respiratory tract. Of the three, B. bronchiseptica is thought to be ancestral to B. parapertussis and B. pertussis although B. parapertussis is more closely related to B. bronchiseptica than to B. pertussis . While B. bronchiseptica can infect a number of mammals, B. pertussis is a strict human pathogen and B. parapertussis infects both humans and sheep.

More detailed analysis revealed a relative of IS481, IS1002, with identical IR but only 61.5% sequence identity was present in B. pertussis and B. parapertussis strains from humans but absent in B. parapertussis isolated from sheep, suggesting that little or no transmission has occurred between sheep and humans.

It was postulated that, once in the human host, B. parapertussis probably acquired IS1002 from B. pertussis. In contrast to human B. parapertussis, isolates, B. pertussis strains produced polymorphic IS1002-related DNA fingerprint patterns .

Genomic sequences, IS expansion, genome streamlining and decay

A remarkable result came from a comparison of the whole genome sequences of the three Bordatellae: B. bronchiseptica RB50; B. parapertussis 12822); and B. pertussis Tohama I (Parkhill et al., 2003 PMID: 12910271)(Fig. IS481.4) showing a large expansion in one particular IS, IS481 accompanied by generation of an extensive collection of pseudogenes. (see……..????)

This analysis confirmed that B. parapertussis and B. pertussis are independent derivatives of B. bronchiseptica-like ancestors (Musser et al., 1986 PMID: 3957867). Indeed, it has even been proposed that B. pertussis has evolved from a human-associated B. bronchiseptica lineage (Diavatopoulos et al, 2005 PMID: 16389302). Moreover, the process was driven by prophage loss from B. bronchiseptica and a massive amplification of IS481/IS1002 copies (261 in B. pertussis and 112 in B. parapertussis; Preston et al 2004 PMID: 15100691) accompanied by significant genome reduction and extensive gene inactivation and loss (Fig. IS481.XXX) and in the two host-restricted species. Although there are a number of collinear regions between B. bronchiseptica and B. parapertussis, most of the rearrangements involve either IS1001 or the IS481-related IS1002.

The rearrangements and deletions in B. pertussis were similar to those in B. parapertussis but more extreme. Short, highly conserved, segments were observed to be fragmented into ~150 individual rearrangements, with the vast majority (88%) bounded primarily by copies of IS481. Thus, the authors point out, adaptation to specific hosts results in loss rather than gain of function. A study which also included B. holmesii found that inversions in the genomes of B. holmesii, B. pertussis, and B. parapertussis were generally flanked by the multicopy IS which occur in their genomes (Wigland et al PMID: 31744907).

Changes at the gene level

There is little net gain of genes in the host restrictive species: very few genes were present in B. pertussis and B. parapertussis but missing in B. bronchiseptica RB50. On the other hand, B. bronchiseptica had over 600 genes not present in either of the other two species: over 1,000 were shared only with B. parapertussis while just over 100 were shared only with B. pertussis. Moreover, IS expansion was accompanied by a large increase in the number of pseudogenes (Fig. IS481.3).

General

Initially, IS481 appeared to be an IS3 family derivative which had been truncated for the N-terminal end of the Tpase and includes a C-terminal extension. The DDE active site domain and the IR (ending in 5’ TGT 3’) are similar to those of IS3 family members. Their presence in high copy number in some species and the identification of at least 130 distinct but related IS from over 90 species strongly suggests that these represent a distinct transpositionally active family.

Different members generate DR of between 4 and 15 bp. Moreover, certain members (e.g. ISSav7) insert specifically into the tetranucleotide CTAG which becomes the flanking DR and provides the UAG termination codon for the Tpase. In contrast to the vast majority of IS3 family members, the IS481 Tpase is not produced by frameshifting. There is no evidence for a leucine zipper as in IS3.

Some members include passenger genes including antibiotic resistance (Chloramphenicol - CmR gene for IS5564 and ISCgl1), or potential transcriptional regulators (ISKrh1, ISPfr21, ISSav7). IS481 itself has played a fundamental role in the evolution of the genomes of the Bordetellae where, in B. pertusis it has undergone extensive amplification to several hundred copies with accompanying genome decay (Fig.IS481.1),

Fig. IS481.1. Colinearity between the genomes of B. pertussis (Top), B. bronchiseptica (Center) and B. parapertussis (Bottom). Top: B. pertussis. Black triangles represent IS481 elements. Center: B. bronchiseptica. Pink boxes represent prophage. Bottom: B. parapertussis. Black triangles represent IS481 elements. The red lines between the genomes represent BLASTN matches between the two sequences. (extracted from: Parkhill J, Sebaihia M, Preston A, Murphy LD, Thomson N, Harris DE, et al. Comparative analysis of the genome sequences of Bordetella pertussis, Bordetella parapertussis and Bordetella bronchiseptica. Nat Genet 2003;35:32-40.)

These IS are distantly related to the eukaryotic Banshee transposon which at present is restricted to the anaerobic flagellated protozoan Trichomonas vaginalis (Pritham per. comm.). They share the highly conserved Pfam integrase core domain identified initially in the IS3 family and retroviruses. They also show a conserved 5’ TG 3’ tip to the IR which is typical of this and other types of mobile elements. It would be interesting to determine whether Banshee transposes using a dsDNA circular intermediate as do IS3 family members.

IS1202 group (ISNCY)

A small group including IS1202, which had been included in the ISNCY (Not Classified Yet) group appears distantly related to IS481. Members are between 1400 and 1700 bp (except for ISKpn21 which includes a passenger gene annotated as "hypothetical protein”) with a Tpase orf of between 400 and 500 aa in a single reading frame. Their IR begin with TGT as do those of the IS3 and IS481 families. They generate DR of between 5 and, unusually, 27 bp.

They appear to have similarities at the level of their Tpases particularly in their DDE domains (e.g. IS1202 is 39% aa similar to ISPfr5 of the IS481 family). Furthermore, they include glutamine (Q) seven residues C-terminal to the conserved E instead of the characteristic K/R. Identification of additional IS will be necessary to clearly define this group.

---