Abstract
Background:
Macrolide resistance genes (erm and mef families) and heavy metal resistance genes (cadA) are increasingly disseminated among streptococci via diverse mobile genetic elements.
Methods:
Through whole-genome sequencing of 16 Streptococcus agalactiae isolates resistant to both erythromycin and clindamycin, we identified 19 integrative and conjugative elements (ICEs), a type of self-transfer genetic elements, conferring antibiotic resistance. Among these, a novel composite ICE designated ICESag39 was identified in S. agalactiae Sag39 through comparative analysis with the NCBI database.
Results:
ICESag39 measured 113,125 bp in length, and it featured a nested “Russian doll” structure comprising an ICESa2603 family backbone integrated with an internal Tn1806-like ICE. The embedded Tn1806-like ICE contained four variable regions (VR1-VR4) that serve as insertion hotspots; among these, VR3 and VR4 carry erm(TR) from ICESp2907 and the cadmium resistance gene cadA from an uncharacterized genetic element, respectively. Conjugation and excision assays confirmed that ICESag39 transfers at a frequency of 8.2 × 10−9 and co-transfers both resistance determinants. Under cadmium stress, transconjugants carrying ICESag39 displayed enhanced growth relative to the recipient. Although the internal Tn1806-like ICE was also capable of independent transfer, its efficiency was significantly lower (< 10−9), and its circular form is undetectable by PCR. Database screening identified 199 structurally similar ICEs (ICESag39-like ICE), 62.8% (125/199) of which co-carried erm(TR) and cadA, underscoring the prevalence of this ICE and its associated resistance traits.
Conclusions:
This study characterizes a new composite ICE and elucidates a modular mechanism that facilitates the assembly and dissemination of resistance genes, thereby promoting bacterial genome diversification.
1 Introduction
Macrolide resistance in streptococci is primarily mediated through two distinct mechanisms: active drug efflux and ribosomal target site modification (; ). The first mechanism involves efflux pumps encoded by mef genes [such as mef (E), mef (A), and mef (I)], which reduce intracellular concentrations by expelling the antibiotic from the cell. This efflux-based strategy confers low-level resistance to 14- and 15-membered macrolides, and it is frequently observed in species such as Streptococcus pneumoniae and Streptococcus pyogenes. The second mechanism entails target site modification mediated by methylases encoded by erm family genes [including erm(A), erm(B), and erm(TR)]. These enzymes methylate the 23S rRNA of the bacterial ribosome, thereby inhibiting antibiotic binding. This modification confers high-level resistance to macrolides, lincosamides, and streptogramin B (the MLSB phenotype), representing the most prevalent and clinically significant macrolide resistance determinants in streptococci and enterococci ().
The widespread dissemination of erm family genes in S. agalactiae and other streptococci is largely driven by integrative and conjugative elements (ICEs) (; ; ). These mobile genetic elements excise from the chromosome, form a circular intermediate, and transfer between bacterial cells via conjugation machinery such as type IV secretion systems (T4SSs) and type IV coupling proteins (T4CPs) (). Following transfer, a site-specific integrase catalyses ICE integration into specific genomic sites, facilitating the efficient spread of antibiotic resistance genes (ARGs) (). Beyond antibiotic resistance, some ICEs can also transmit heavy metal resistance determinants, conferring related tolerance (; ).
Two groups of ICEs are currently recognized as common vehicles for the dissemination of erm genes in streptococci, namely the ICESa2603 and Tn1806 families. The ICESa2603 family, first identified in S. agalactiae 2603V/R, is the most frequently reported group (; ; ). These elements often carry multiple antibiotic resistance genes and integrate site-specifically 3′ rplL gene, which encodes the 50S ribosomal subunit protein L7/L12. Representative members include ICESsu32457 from Streptococcus suis, carrying aadE-aphA (conferring aminoglycoside resistance), erm(B) and tet(40)-tet(O/W/32/O) (tetracycline resistance) (), as well as ICESpn529IQ from S. pneumoniae, which harbors catQ (chloramphenicol resistance), tetM (tetracycline resistance), mef (I), and erm(TR) (). In addition, some ICESa2603-like family ICEs (characterized by distinct integrase modules and integration into alternative sites) also carry erm genes, such as ICESsuYZDH1 from S. suis () and ICESag37 from S. agalactiae (). The second group, Tn1806, was initially identified in S. pneumoniae, and it is characterized by its low-integration-site specificity and high transfer efficiency, promoting its dissemination across streptococci and enterococci (; ). However, ICEs are dynamic and modular. Genetic exchange and interaction can occur between different ICE groups, potentially giving rise to chimeric elements with novel characteristics and resistance profiles (). For instance, one ICE or mobile genetic element can integrate into a specific site within another, forming larger, more complex composite structures reminiscent of a “Russian doll” (; ). Although such composite elements/ICEs have rarely been documented, they can accumulate advantageous genes from multiple ICEs or other mobile genetic elements, thereby enhancing bacterial survival under diverse antibiotic pressures.
In this study, we analyzed the genetic context of antibiotic resistance genes in 16 macrolide-resistant isolates of S. agalactiae and identified a novel composite ICE. This ICE incorporates modules homologous to several previously characterized genetic elements, including an ICESa2603 family backbone, a Tn1806-like ICE, the erm(TR)-carrying element ICESp2907, and an uncharacterized cadA-carrying segment. These genetic elements are assembled at appropriate genomic locations according to their functional properties, resulting in a chimeric ICE that combines structural and functional features from distinct ancestral elements.
2 Materials and methods
2.1 Bacterial strains and susceptibility tests
This study included 47 non-duplicate clinical isolates of S. agalactiae collected and stored by the Shanghai Center for Clinical Laboratory (Shanghai, China) between 2021 and 2025. Among these, 16 isolates exhibited resistance to both erythromycin and clindamycin [minimum inhibitor concentration (MIC) > 1 mg/L]. Recipient strains for conjugation experiments were selected from the 47 S. agalactiae isolates. All bacterial strains were cultured in Todd–Hewitt (TH) broth at 37 °C for 8 h. Antimicrobial susceptibility testing was performed using the broth microdilution method, with breakpoints interpreted according to CLSI (the Clinical and Laboratory Standards Institute) guidelines. The reference strain S. agalactiae ATCC 13813 was obtained from the Huashan Hospital strain repository (Shanghai, China). The clinical isolates were partially collected through routine hospital laboratory procedures, and no additional personal information was involved in this study.
2.2 Genome comparative analysis and ICE identification
Genomic DNA was extracted from S. agalactiae using the QIAGEN Midi Kit (Qiagen, Hilden, Germany) and sequenced by next-generation sequencing (NGS) at BGI Genomics (Shenzhen, China) using HiSeq X (Illumina, San Diego, CA, USA). High-quality reads were assembled with SPAdes (), and the assembled sequences were annotated through the RAST online server (rast.nmpdr.org) (). Contigs were ordered against the complete reference genome S. agalactiae 2603V/R using the “Order Contigs” function in Mauve v2.4.0 (; ). The reordered genomes were then aligned and visualized against the reference using the “Comparison” function in Mauve, enabling preliminary localization of large genetic elements. Putative mobile genetic elements were identified by screening for genes encoding T4SSs, T4CPs, relaxase, and integrase. These candidate regions were further analyzed using ICEfinder to confirm their classification as ICEs.
2.3 Amplification experiments
PCR assays were performed for three objectives (Supplementary Table 1): to screen candidate recipient strains by detecting the presence of an unoccupied rplL gene; to amplify conserved ICE genes in trans-conjugants for preliminary confirmation of ICE integration; and to examine ICE integration, excision, and circularization using the primer pairs P1/P2 and P3/P4 to detect integration, the primer pair P1/P4 to detect excision, and the inverse primer pair P2/P3 to detect circular intermediates.
2.4 Conjugative transfer assays
Recipient strains were selected from clinical isolates based on two criteria: a resistance profile distinct from that of the donor strain and the presence of at least one unoccupied ICE integration site. These characteristics were initially assessed using antimicrobial susceptibility testing and PCR and subsequently confirmed by whole-genome sequencing.
Conjugation was performed using a filter mating method as previously described (). Fifty mL Donor and 250 mL recipient strains were grown to the late-log phase, normalized to the same CFU, and mixed at a donor-to-recipient ratio of 1:5. Prior to mixing, 200 μL of the donor culture was sampled for CFU enumeration to calculate the conjugation frequency. A nitrocellulose membrane was cut to an appropriate size, treated with 10 μg/mL DNase to prevent transformation-mediated DNA transfer, and placed on a TH agar plate. The mixed culture was pelleted by centrifugation, resuspended, and evenly spread onto the DNase-treated membrane, followed by incubation at 37 °C for 4 h. After incubation, cells were harvested from the membrane, serially diluted, and plated onto TH agar containing 50 μg/mL erythromycin, 50 μg/mL clindamycin, and 20 μg/mL levofloxacin. Transconjugant colonies were counted after 2–3 days of incubation and randomly selected for validation.
The conjugation frequency was calculated as the number of transconjugants divided by the number of donor cells. Donor count: 100 μL culture was taken from a 50 mL donor culture and serially diluted 10-fold in 1 × PBS for 6-10 times. After thorough mixing, 100 μL dilution was spread evenly onto counting plates. Colony counts between 30 and 300 were considered valid. Transconjugant count: After mating, cells were collected from the nitrocellulose membrane using a sterile swab and resuspended. The suspension was then diluted 10-fold, 100-fold, or used undiluted, and plated onto transconjugant selective medium. Plates with 30-300 colonies were considered valid. If the undiluted suspension yielded transconjugants but fewer than 30 colonies, the transfer was considered successful, but the frequency was too low to be calculated accurately.
2.5 Growth curve analysis under cadmium stress
To evaluate the growth dynamics of the recipient and transconjugants under cadmium stress, growth curves were generated in TH broth with or without cadmium as previous described with minor modification (). The recipient and transconjugant strains were inoculated into two culture conditions: unsupplemented TH broth and TH broth supplemented with 15 μM CdCl2 to induce metal stress. Cultures were incubated at 37 °C with shaking at 160 rpm, and optical density at 600 nm was measured hourly using a microplate reader until the stationary phase was reached. Growth curves were plotted using GraphPad Prism (GraphPad, Boston, MA, USA), and the growth of transconjugants was compared with that of the recipient strain under cadmium stress. All experiments were performed in at least two independent biological replicates to ensure reproducibility.
2.6 Bioinformatics identification of ICEs similar to composite ICESag39
The core genes of ICESag39, including the tandem serine integrase, relaxase, and conjugation modules (including T4SS and T4CP cluster) of the Tn1806 family, along with the tyrosine integrase, relaxase, and conjugation modules of the ICESa2603 family, were used as templates for BlastN analysis in Nucleotide Collection and Whole-genome shotgun contigs databases from NCBI. According to the criteria defined by ICEFinder (db-mml.sjtu.edu.cn/ICEberg) and previous reports (; ), a genomic region was considered to harbor an ICE similar to ICESag39 if it contained genes homologous to these core genes (>60% nucleotide identity); the maximum distance between these genes was < 100 kb; Minimum sequence coverage of BLAST hits >25%. In contigs database, a sequenced strain was considered to contain the composite ICE only when the complete element was identified on one scaffold/contig. Boundaries and attachment sites of ICEs were manually checked by terminal integrase and rplL genes. ICE-carried AMR genes were identified as previously described.
2.7 GeneBank accession number
Sequence data generated in this study have been deposited in the NCBI GenBank, including the NGS drafts of 16 erythromycin- and clindamycin-resistant S. agalactiae isolates, the complete genomes of recipient S. agalactiae SagR272 and SagR31, and the NGS draft of transconjugants SagR272_TC and SagR31_TC. Public data used in this study included Tn1806, ICESp2907, ICESpy009, ICESpy009, and the complete genome of S. agalactiae 2603V/R.
3 Results
3.1 Characterization of erythromycin- and clindamycin-resistant isolates
The NGS data of 16 erythromycin- and clindamycin-resistant S. agalactiae isolates revealed diverse phenotypic and genotypic profiles (Table 1). Ten carried the erm(B) gene, and one contained the erm(TR) gene. Three isolates harbored the macrolide resistance genes mef (E) and/or mel, whereas tetracycline resistance, associated with the tet(M) gene, was observed in seven isolates.
Table 1
| Strains | Year | MIC (mg/L) | Resistance genotype | ARGs-carrying genetic elements | Serotype | Accession number | ||||
|---|---|---|---|---|---|---|---|---|---|---|
| ERY | CLI | TET | LEV | CHL | ||||||
| Sag191 | 2021 | >256 | >256 | 128 | ≤0.5 | ≤0.5 | erm(B), tetM | Tn1806-like ICE | ST1 | - |
| Sag199 | 2021 | >256 | >256 | ≤0.5 | ≤0.5 | ≤0.5 | erm(B) | Tn1806-like ICE | ST10 | - |
| Sag228 | 2022 | >256 | >256 | ≤0.5 | ≤0.5 | ≤0.5 | erm(B) | Tn1806-like ICE | ST10 | - |
| Sag257 | 2022 | 64 | >256 | 64 | 32 | 1 | mef(E)-mel, tetM | ICESa2603, Tn916 | ST19 | JBAPEL000000000 |
| Sag273 | 2022 | >256 | >256 | ≤0.5 | ≤0.5 | ≤0.5 | erm(B), mel | Tn1806-like ICE | ST10 | - |
| Sag2111 | 2022 | 4 | 2 | 32 | ≤0.5 | ≤0.5 | tetM | Tn916 | ST23 | JBAPEO000000000 |
| Sag2142 | 2022 | >256 | >256 | ≤0.5 | ≤0.5 | ≤0.5 | erm(B) | Tn1806-like ICE | ST10 | - |
| Sag2157 | 2022 | >256 | >256 | ≤0.5 | ≤0.5 | ≤0.5 | erm(B) | Tn1806-like ICE | ST10 | - |
| Sag2167 | 2022 | >256 | >256 | ≤0.5 | ≤0.5 | ≤0.5 | erm(B) | Tn1806-like ICE | ST10 | JBAPEQ000000000 |
| Sag35 | 2023 | >256 | >256 | ≤0.5 | ≤0.5 | ≤0.5 | erm(B) | Tn1806-like ICE | ST10 | - |
| Sag39 | 2023 | >256 | >256 | >256 | ≤0.5 | ≤0.5 | erm(TR), tetM | ICESag39, Tn916 | ST19 | JBAPEP000000000 |
| Sag3155 | 2023 | >256 | >256 | 1 | ≤0.5 | ≤0.5 | erm(B) | Tn916 | ST10 | - |
| Sag3220 | 2023 | >256 | 256 | >256 | ≤0.5 | ≤0.5 | mef(E)-mel, tetM | ICESa2603 | ST17 | JBAPER000000000 |
| Sag4240 | 2024 | >256 | >256 | ≤0.5 | 2 | ≤0.5 | erm(B) | Tn1806-like ICE | ST10 | JBAPES000000000 |
| Sag4272 | 2024 | 1 | 1 | >256 | ≤0.5 | ≤0.5 | tetM | Tn916 | ST10 | JBAPEB000000000 |
| Sag516 | 2025 | 1 | 2 | 32 | 1 | 1 | tetM | Tn916 | ST10 | JBAPEH000000000 |
Phenotypic, genotypic, and other characteristics of the 16 erythromycin- and clindamycin-resistant Streptococcus agalactiae isolates.
MIC, minimum inhibitory concentration; ERY, erythromycin; CLI, clindamycin; TET, tetracycline; LEV, levofloxacin; CHL, chloramphenicol; ARG, antibiotic resistance gene. The “-” indicate that the corresponding sequence was not essential for this study and therefore has not been uploaded to the NCBI database. Bolded entries represent the primary subjects, host of ICESag39, of this investigation.
Multi-locus sequence typing identified ST10 as the most predominant sequence type (11 isolates), followed by ST19 (two isolates). ST1, ST17, and ST23 (one isolate each) were also detected.
3.2 Identification of antibiotic resistance gene-carrying genetic elements
To localize antibiotic resistance genes and identify their associated genetic elements, we compared the genomes of the 16 S. agalactiae isolates with that of the reference strain 2603V/R (). This analysis identified 35 major genomic regions of difference, 19 of which contained antibiotic resistance genes along with integrase genes, suggesting integration of mobile genetic elements into the genome (Figure 1). Comparative analysis using GenBank classified the antibiotic resistance gene-carrying elements into several known families: the Tn1806 family carrying erm(B), the Tn916 family carrying tet(M), and the ICESa2603 family carrying mef (E)-mel (Table 1). A novel composite genetic structure carrying erm(TR) was identified in S. agalactiae Sag39, and it was formed by the integration of Tn1806- and ICESa2603-like ICE 3′ rplL gene. This element also carried the cadmium resistance gene cadA, and because of the presence of multiple integrase and conjugation-related genes, it was classified as a composite ICE and designated ICESag39.
Figure 1
3.3 Characterization of ICESag39 and the internal Tn1806-like ICE
ICESag39 measured 112,852 bp in length and featured 102 open-reading frames and 35% GC content. The entire ICE was flanked by direct repeats (TTATTTAAGAGTAAC–TTATTTAAGAGTAAC) 3′ rplL gene. Comparative genomic analysis revealed that ICESag39 comprises four modules (Figure 2A). The first module shared high similarity with ICESpy009, a member of the ICESa2603 family, but lacked the mef (E)-mel gene cluster present in S. pyogenes (). A Tn1806-like ICE inserted within the snf2 (methyltransferase encoding the sucrose non-fermenting 2 protein) gene of the ICESpy009 backbone constituted the second module. Notably, this internal Tn1806-like ICE itself carried two additional mobile genetic units: ICESp2907, which lacks a conjugation system but carries erm(TR) (), and an uncharacterized element harboring cadA.
Figure 2
To explore the structural organization and integration hotspots for exogenous DNA, we analyzed the relatively variable Tn1806-like ICE. This analysis involved 41 additional ICEs with high sequence similarity to Tn1806 from NCBI databases, and 25 of these ICEs carried at least one antibiotic resistance genes (Supplementary Table 2). These elements integrate with low specificity into a 4-bp direct repeat sequence (GGGG) within snf2 or diverse methyltransferase encoding genes, truncating the target gene upon insertion. Beyond the conserved integrase and conjugation modules, four variable regions (VR1–VR4) were identified as hotspots for the integration of foreign genetic elements (Figure 2B, Supplementary Table 3). VR1 and VR2 frequently harbored erm(A) and erm(B), along with the aminoglycoside resistance gene aac(6′)-aph(2″) in some cases. Entire mobile elements such as Tn5253 (), Tn916 (), and Tn6002 () were also identified within these regions. VR3 served as a specific integration site for ICESp2907, which carries erm(TR). VR4, uniquely identified in ICESag39, contained the cadA-carrying genetic element.
3.4 Excision and circulation of ICESag39 and the internal Tn1806-like ICE
Chromosomal excision and subsequent circularization are essential for the rolling circle replication and conjugative transfer of ICEs. To preliminarily assess the transfer potential of ICESag39 and the embedded Tn1806-like ICE, we designed PCR primers to detect these processes. Amplification with the primer pairs P1/P2 and P3/P4 as well as Pi1/Pi2 and Pi3/Pi4 (targeting ICE-flanking regions) yielded positive products for both ICESag39 and the internal Tn1806-like ICE (Figure 3), confirming their chromosomal integration. However, only ICESag39 produced faint PCR products with the primer pairs P1/P4 (excision) and P2/P3 (circular intermediate), whereas corresponding amplicons Pi1/Pi4 and Pi2/Pi3 were undetectable for the Tn1806-like ICE (Figure 3). These results suggest that ICESag39 is capable of excising and forming circular intermediates, whereas the internal Tn1806-like ICE appears incapable of autonomous excision, or its transfer frequency is too low to be detected by PCR.
Figure 3
3.5 Transferability of ICESag39
S. agalactiae Sag39 carrying ICESag39 was used as a donor (resistant to erythromycin, clindamycin, but sensitive to levofloxacin). The recipient strain, S. agalactiae SagR272 (sensitive to erythromycin, clindamycin, but resistant to levofloxacin), which possesses an unoccupied rplL gene for ICESag39, was selected (Figure 4A). In mating experiments, the complete ICESag39 was transferred from the donor to the recipient at a frequency of 8.2 × 10−9 (Table 2). Three randomly selected transconjugants acquired high-level resistance to erythromycin and clindamycin while retaining the recipient's resistance to levofloxacin. Under cadmium stress (15 μM CdCl2 in TH broth), transconjugants exhibited significantly better growth than the recipient (Figure 4D), suggesting that transfer of cadmium resistance. Whole-genome sequencing confirmed integration of ICESag39 downstream of the rplL gene in the transconjugants SagR272_TC comparing with recipient SagR272, as expected (Figure 4C).
Figure 4
Table 2
| Strains | Purpose | MIC (μg/mL) | Transfer frequency | Putative integration sites | Transfer ICE | Accession No. | Objective ICE location | |||
|---|---|---|---|---|---|---|---|---|---|---|
| ERY | CLI | TET | LEV | |||||||
| Sag39 | Donor | >256 | 256 | 32 | ≤0.5 | JBAPEP000000000 | Sequence02, 12144-62118 | |||
| SagR272 | Recipient | ≤0.5 | 1 | 128 | 32 | 8.2 × 10−9 | rplL | ICESag39 | CM130719 | |
| SagR272_TC | Transconjugant | >256 | >256 | 128 | 32 | JBRZQM000000000 | ||||
| SagR31 | Recipient | 4 | ≤0.5 | 4 | 64 | <10−9 | hsdM | Tn1806-like ICE | CP138369 | |
| SagR31_TC | Transconjugant | >256 | 256 | 4 | 64 | JBRZQN000000000 | ||||
Characteristics of the donor, recipient, and related transconjugants.
The bold values were obtained by antibiotic resistance screening via a conjugation transfer experiment. MIC, minimum inhibitory concentration; ERY, erythromycin; CLI, clindamycin; TET, tetracycline; LEV, levofloxacin; ICE, integrative and conjugative element.
To assess whether the internal Tn1806-like ICE could excise and transfer independently, another isolate, S. agalactiae SagR31 (sensitive to clindamycin, but resistant to erythromycin and levofloxacin), was used as the recipient. This strain contains an hsdM gene (encoding type I methyltransferase), a known integration site for Tn1806 family ICEs, and the rplL site has been occupied by a resident ICE (Figure 4B). After three independent mating experiments, the Tn1806-like ICE transferred only once. The transconjugant counts obtained were 0, 0, and 7, respectively, none of which reached the minimum countable threshold of 30 and were therefore considered as < 30; the donor counts in the mating experiments were 3.4 × 1010 CFU, 4.1 × 1010 CFU, and 3.2 × 1010 CFU. Therefore, calculated conjugation transfer frequency was < 10−9. Whole-genome sequencing confirmed Tn1806-like ICE were compensatorily integrated into hsdM genes in recipient SagR31 (Figure 4C).
3.6 Prevalence of composite genomic elements resembling ICESag39
We screened the NCBI GenBank database using BlastN to identify composite ICEs structurally similar to ICESag39 in S. agalactiae. Details of the identified ICE, including the host strain, accession number, genomic location, length, and carried antibiotic resistance genes, are summarized in Supplementary Table 4. In total, 199 ICESag39-like ICEs were identified, including seven from complete genomes and 192 from whole-genome shotgun contigs, all integrated 3′ rplL gene. These strains were unloaded from different years (2012–2025), with 74.4% (148/199) from the United States, 16.6% (33/199) from China, 4.0% (8/199) from South Korea, 2.0% (4/199) from Brazil, 1.5% (3/199) from the United Kingdom, 1.0% (2/199) from Italy, and 0.5% (1/199) from Japan. Their sizes ranged from 90,473 to 138,541 bp. Analysis of resistance gene carriage revealed that 62.8% (125/199) of these ICEs harbored both erm(TR) and cadA, 22.1% (44/199) carried only cadA, and 4.5% (9/199) carried only erm(TR). Additionally, 1.5% (3/199) carried aac(6′)-aph(2″) together with erm(TR), whereas 7.0% (14/199) possessed erm(TR) along with catQ–mef (I). Four ICEs lacked any known resistance genes. Notably, nearly all resistance genes were located within the embedded Tn1806-like ICE module of these composite elements.
4 Discussion
Macrolide resistance determinants, such as erm family genes and the heavy metal resistance gene cadA, are widely recognized as common resistance factors. The erm(B) gene in particular exhibits a detection rate as high as 89.3% in S. agalactiae, posing a challenge of global significance (; ). In 2019, the US Centers for Disease Control and Prevention identified 21 antibiotic resistance threats, classifying clindamycin-resistant S. agalactiae and erythromycin-resistant S. pyogenes at the “concerning” level (). Moreover, cadA has frequently been detected on plasmids from erm(B)- or mef (A)-positive isolates, indicating a potential linkage between these commonly co-occurring resistance genes ().
This study identified a novel composite ICE characterized by a nested architecture, formed by an ICESa2603 family backbone integrated with a Tn1806-like ICE, spanning 112,852 bp in length. The embedded Tn1806-like ICE constitutes an autonomous genetic module featuring four variable regions (VR1–VR4) that act as insertion hotspots for additional resistance determinants. Specifically, ICESp2907 carrying erm(TR) and an uncharacterized cadA-carrying element were identified within VR3 and VR4, respectively (). Comparative genomic analysis further revealed that VR1 and VR2 might harbor resistance genes such as aac(6′)-aph(2″) (aminoglycoside), erm(B) (MLSB), and tet(M) and potentially accommodate entire mobile elements such as Tn6002, Tn916, and Tn5253. The backbone of ICESag39 represents a typical member of the ICESa2603 family, as evidenced by its tyrosine integrase, specific integration site rplL gene, and flanking 15-bp direct repeat sequences (; ). This backbone provides a stable genomic platform for the acquisition of exogenous DNA and rplL as the defined integration locus. By contrast, the embedded Tn1806-like ICE substantially expands the accessory gene content available to the ICESa2603 family by serving as an insertion site for smaller genetic elements. This Russian doll-like nesting of mobile genetic elements produced the formation of the composite ICE ICESag39, facilitating the co-localization and potential co-transfer of diverse resistance genes.
Through a systematic comparative analysis of Tn1806-like ICEs, this study first expanded the recognized diversity of this ICE group and elucidated the molecular basis for the formation of the composite ICE ICESag39. The Tn1806-like ICE identified in this study was integrated within the 6-kb (GGG–GGG) snf2 gene of ICESag39, an integration site that differs from the previously characterized hsdM (encoding a type I methyltransferase) target observed in Tn1806 and ICESpy10750 (; ; ). An expanded screening of Tn1806-like ICEs identified four additional Tn1806-like ICEs from Gen Bank also integrated within the snf2 gene, indicating that snf2 represents a previously overlooked integration site for this ICE family. Notably, snf2 has been reported as a conserved gene in the ICESa2603 family, and it is also present in enterococcal pheromone-responsive plasmids such as pBET_5 (; ). The targeting of snf2 as an integration site strongly suggests that the composite structure of ICESag39, arising from the combination of ICESa2603 and Tn1806-like ICEs, is not an occasional event, but reflects a high-affinity interaction. This interpretation is further supported by the identification of 199 similar ICESag39-like ICEs in public databases. However, no composite plasmid-ICE structures associated with Tn1806-like ICEs or enterococcal pheromone-responsive plasmids were detected, despite the presence of the snf2 gene in such plasmids and prior evidence supporting the interspecies transferability of Tn1806-like ICEs. This absence might reflect specific genetic or physiological barriers that limit the formation of such composite, a possibility that remains to be investigated.
ICESag39 represents a novel class of composite ICEs that was demonstrated for the first time in this study to possess a modular architecture and to be transferable to other Sagalactiae strains. The embedded Tn1806-like ICE, which acts as a primary carrier of antibiotic and heavy metal resistance genes, is not a fixed component of the structure. Conjugation experiments confirmed that individual internal ICEs within ICESag39 can be excised from the composite and transferred independently. However, detection of circular intermediates and conjugation frequency indicated that the transfer efficiency of internal ICEs is remarkably low, suggesting that ICESag39 generally maintains its structural integrity under most conditions. In situations unfavorable for the transfer of the composite ICE, such as when its length imposes a high fitness cost or when suitable integration sites are absent in the recipient, the ability of internal Tn1806-like ICEs to transfer independently might provide an alternative dissemination pathway. This structural flexibility highlights a distinct adaptive advantage of composite ICEs over simple genetic elements. The assembly of erm(TR) and cadA within a single transferable ICE provides important insights into the evolution of ICE diversity and the horizontal dissemination of resistance determinants. Integrating these advantageous genes into a large, mobile composite ICE both increases the risk of multidrug-resistant “superbug” emergence and accelerates the diversification of bacterial genomes.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.
Author contributions
SY: Validation, Resources, Writing – review & editing, Funding acquisition, Formal analysis, Project administration, Conceptualization, Methodology, Writing – original draft, Data curation, Visualization, Investigation. CS: Methodology, Investigation, Writing – original draft, Formal analysis. LY: Software, Writing – original draft, Formal analysis. XX: Formal analysis, Writing – original draft, Software. XW: Writing – original draft, Software, Formal analysis.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This work was supported by the Clinical Research Special Program of Shanghai Municipal Health Commission [Grant Number 2025Y0054].
Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2026.1763839/full#supplementary-material
References
1
AmbrosetC.ColuzziC.GuédonG.DevignesM. D.LouxV.LacroixT.et al. (2016). New insights into the classification and integration specificity of Streptococcus integrative conjugative elements through extensive genome exploration. Front. Microbiol.6:1483doi: 10.3389/fmicb.2015.01483
2
AzizR. K.BartelsD.BestA. A.DejonghM.DiszT.EdwardsR. A.et al. (2008). The RAST Server: rapid annotations using subsystems technology. BMC Genom.9:75. doi: 10.1186/1471-2164-9-75
3
BankevichA.NurkS.AntipovD.GurevichA. A.DvorkinM.KulikovA. S.et al. (2012). SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol.19, 455–477. doi: 10.1089/cmb.2012.0021
4
BeresS. B.MusserJ. M. (2007). Contribution of exogenous genetic elements to the group A Streptococcus metagenome. PLoS ONE2:e800. doi: 10.1371/journal.pone.0000800
5
BrencianiA.BacciagliaA.VecchiM.VitaliL. A.VaraldoP. E.GiovanettiE. (2007). Genetic elements carrying erm(B) in Streptococcus pyogenes and association with tet(M) tetracycline resistance gene. Antimicrob. Agents Chemother.51, 1209–1216. doi: 10.1128/AAC.01484-06
6
CamilliR.Del GrossoM.IannelliF.PantostiA. (2008). New genetic element carrying the erythromycin resistance determinant erm(TR) in Streptococcus pneumoniae. Antimicrob. Agents Chemother.52, 619–625. doi: 10.1128/AAC.01081-07
7
CDC (2019). Centers for Disease Control and Prevention. Avaliable online at: http://www.cdc.gov/drugresistance/threat-report-2019. (Accessed December 01, 2023)
8
ChenL.HuangJ.HuangX.HeY.SunJ.DaiX.et al. (2021). Horizontal transfer of different erm(B)-carrying mobile elements among Streptococcus suis strains with different serotypes. Front. Microbiol.12:628740. doi: 10.3389/fmicb.2021.628740
9
CostaS. S.PalmaC.KadlecK.FesslerA. T.ViveirosM.Melo-CristinoJ.et al. (2016). Plasmid-borne antimicrobial resistance of Staphylococcus aureus isolated in a hospital in Lisbon, Portugal. Microb. Drug Resist.22, 617–626. doi: 10.1089/mdr.2015.0352
10
DarlingA. C.MauB.BlattnerF. R.PernaN. T. (2004). Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res.14, 1394–1403. doi: 10.1101/gr.2289704
11
Del GrossoM.CamilliR.RizziE.PietrelliA.De BellisG.PantostiA. (2016). ICESpy009, a conjugative genetic element carrying mef (E) in Streptococcus pyogenes. Antimicrob. Agents Chemother.60, 3906–3912. doi: 10.1128/AAC.03082-15
12
GaworJ.ZuchniewiczK.Ojeda SaavedraM.BeresS. B.KiedrowskaM.Wrobel-PawelczykI.et al. (2025). ICESp1109, a Novel hybrid integrative conjugative element of macrolide-resistant Streptococcus pyogenes serotype M77 collected between 2003 and 2017 in Poland. J. Infect. Dis.231, e521–e530. doi: 10.1093/infdis/jiae473
13
Gawron-BurkeC.ClewellD. B. (1982). A transposon in Streptococcus faecalis with fertility properties. Nature300, 281–284. doi: 10.1038/300281a0
14
GiovanettiE.BrencianiA.TiberiE.BacciagliaA.VaraldoP. E. (2012). ICESp2905, the erm(TR)-tet(O) element of Streptococcus pyogenes, is formed by two independent integrative and conjugative elements. Antimicrob. Agents Chemother.56, 591–594. doi: 10.1128/AAC.05352-11
15
HatrongjitR.BoueroyP.ChopjittP.WongsurawatT.JenjaroenpunP.WankaewN.et al. (2025). Genomic characterization and virulence of Streptococcus suis serotype 7 sequence type 373 of clonal complex 94. Vet. Res.56:125. doi: 10.1186/s13567-025-01562-4
16
HuangJ.LiangY.GuoD.ShangK.GeL.KashifJ.et al. (2016a). Comparative genomic analysis of the ICESa2603 family ICEs and spread of erm(B)- and tet(O)-carrying transferable 89K-subtype ICEs in swine and bovine isolates in China. Front. Microbiol.7:55. doi: 10.3389/fmicb.2016.00055
17
HuangJ.MaJ.ShangK.HuX.LiangY.LiD.et al. (2016b). Evolution and diversity of the antimicrobial resistance associated mobilome in Streptococcus suis: a probable mobile genetic elements reservoir for other streptococci. Front. Cell. Infect. Microbiol.6:118. doi: 10.3389/fcimb.2016.00118
18
HuangJ.SunJ.WuY.ChenL.DuanD.LvX.et al. (2019). Identification and pathogenicity of an XDR Streptococcus suis isolate that harbours the phenicol-oxazolidinone resistance genes optrA and cfr, and the bacitracin resistance locus bcrABDR. Int. J. Antimicrob. Agents54, 43–48. doi: 10.1016/j.ijantimicag.2019.04.003
19
LeclercqR. (2002). Mechanisms of resistance to macrolides and lincosamides: nature of the resistance elements and their clinical implications. Clin. Infect. Dis.34, 482–492. doi: 10.1086/324626
20
MariniE.PalmieriC.MagiG.FacinelliB. (2015). Recombination between Streptococcus suis ICESsu32457 and Streptococcus agalactiae ICESa2603 yields a hybrid ICE transferable to Streptococcus pyogenes. Vet. Microbiol.178, 99–104. doi: 10.1016/j.vetmic.2015.04.013
21
MingoiaM.MoriciE.MariniE.BrencianiA.GiovanettiE.VaraldoP. E. (2016). Macrolide resistance gene erm(TR) and erm(TR)-carrying genetic elements in Streptococcus agalactiae: characterization of ICESagTR7, a new composite element containing IMESp2907. J. Antimicrob. Chemother.71, 593–600. doi: 10.1093/jac/dkv408
22
MingoiaM.MoriciE.MorroniG.GiovanettiE.Del GrossoM.PantostiA.et al. (2014). Tn5253 family integrative and conjugative elements carrying mef (I) and catQ determinants in Streptococcus pneumoniae and Streptococcus pyogenes. Antimicrob. Agents Chemother.58, 5886–5893. doi: 10.1128/AAC.03638-14
23
MoriciE.SimoniS.BrencianiA.GiovanettiE.VaraldoP. E.MingoiaM. (2017). A new mosaic integrative and conjugative element from Streptococcus agalactiae carrying resistance genes for chloramphenicol (catQ) and macrolides [mef (I) and erm(TR)]. J. Antimicrob. Chemother.72, 64–67. doi: 10.1093/jac/dkw367
24
MoroiH.KimuraK.IdoA.BannoH.JinW.WachinoJ. I.et al. (2019). Erythromycin-susceptible but clindamycin-resistant phenotype of clinical ermB-PCR-positive group B streptococci isolates with is1216e-inserted ermB. Jpn. J. Infect. Dis.72, 420–422. doi: 10.7883/yoken.JJID.2019.015
25
PalmieriC.MagiG.MingoiaM.BagnarelliP.RipaS.VaraldoP. E.et al. (2012). Characterization of a Streptococcus suis tet(O/W/32/O)-carrying element transferable to major streptococcal pathogens. Antimicrob. Agents Chemother.56, 4697–4702. doi: 10.1128/AAC.00629-12
26
ParsonsC.LeeS.JayeolaV.KathariouS. (2017). Novel Cadmium resistance determinant in listeria monocytogenes. Appl. Environ. Microbiol.83:e02580. doi: 10.1128/AEM.02580-16
27
PartridgeS. R.KwongS. M.FirthN.JensenS. O. (2018). Mobile genetic elements associated with antimicrobial resistance. Clin. Microbiol. Rev.31, e00088-17. doi: 10.1128/CMR.00088-17
28
RiceL. B. (1998). Tn916 family conjugative transposons and dissemination of antimicrobial resistance determinants. Antimicrob. Agents Chemother.42, 1871–1877. doi: 10.1128/AAC.42.8.1871
29
Rodriguez-BlancoA.LemosM. L.OsorioC. R. (2012). Integrating conjugative elements as vectors of antibiotic, mercury, and quaternary ammonium compound resistance in marine aquaculture environments. Antimicrob. Agents Chemother.56, 2619–2626. doi: 10.1128/AAC.05997-11
30
ShanX.LiX. S.SchwarzS.ChenY.XuC.DuX. D. (2022). Plasmid-assisted horizontal transfer of a large lsa(e)-carrying genomic island in enterococcus faecalis. Microbiol. Spectr.10:e0015422. doi: 10.1128/spectrum.00154-22
31
VaraldoP. E.MontanariM. P.GiovanettiE. (2009). Genetic elements responsible for erythromycin resistance in streptococci. Antimicrob. Agents Chemother.53, 343–353. doi: 10.1128/AAC.00781-08
32
WangR.ChenH.LiuY. (2025a). Metagenomic insights into the characteristics and co-migration of antibiotic resistome and metal(loid) resistance genes in urban landfill soil and groundwater. Environ. Res.285:122285. doi: 10.1016/j.envres.2025.122285
33
WangY.LiuT.SidaY.ZhuY. (2025b). Diversity and evolution of the mobilome associated with antibiotic resistance genes in Streptococcus anginosus. Microb. Drug Resist.31, 52–63. doi: 10.1089/mdr.2024.0229
34
WozniakR. A.WaldorM. K. (2010). Integrative and conjugative elements: mosaic mobile genetic elements enabling dynamic lateral gene flow. Nat. Rev. Microbiol.8, 552–563. doi: 10.1038/nrmicro2382
35
YangY.YangX.ZhouH.NiuY.LiJ.FuX.et al. (2022). Bisphenols promote the pheromone-responsive plasmid-mediated conjugative transfer of antibiotic resistance genes in Enterococcus faecalis. Environ. Sci. Technol.56, 17653–17662. doi: 10.1021/acs.est.2c05349
36
YiS.HuangJ.HuX.ChenL.DaiX.SunJ.et al. (2021). Nonconservative integration and diversity of a new family of integrative and conjugative elements associated with antibiotic resistance in zoonotic pathogen Streptococcus suis. Vet. Microbiol.254:109009. doi: 10.1016/j.vetmic.2021.109009
37
ZhouK.XieL.HanL.GuoX.WangY.SunJ. (2017). ICESag37, a Novel integrative and conjugative element carrying antimicrobial resistance genes and potential virulence factors in Streptococcus agalactiae. Front. Microbiol.8:1921 doi: 10.3389/fmicb.2017.01921
Summary
Keywords
antibiotic resistance genes, antimicrobial susceptibility, cadmium resistance, integrative and conjugative elements (ICEs), Streptococcus agalactiae
Citation
Yi S, Shi C, Yin L, Xu X and Wang X (2026) Characterization of a novel composite ICE in Streptococcus agalactiae conferring resistance to macrolides [erm(TR)] and cadmium (cadA). Front. Microbiol. 17:1763839. doi: 10.3389/fmicb.2026.1763839
Received
09 December 2025
Revised
06 January 2026
Accepted
12 January 2026
Published
09 March 2026
Volume
17 - 2026
Edited by
Fangyu Ren, Nankai University, China
Reviewed by
Xingyang Dai, Nanjing Agricultural University, China
Zhoumiao Chen, First Affiliated Hospital of Wenzhou Medical University, China
Xiao-Qing Zhu, South China Agricultural University, China
Updates
Copyright
© 2026 Yi, Shi, Yin, Xu and Wang.
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.
*Correspondence: Xueliang Wang, xlwang12@126.com
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.