Mariculture waters as yet another hotbed for the creation and transfer of new antibiotic-resistant pathogenome

With the rapid growth of aquaculture globally, large amounts of antibiotics have been used to treat aquatic disease, which may accelerate induction and spread of antibiotic-resistant bacteria (ARB) and antibiotic resistance genes (ARGs) in aquaculture environments. Herein, metagenomic and 16S rRNA analyses were used to analyze the potentials and co-occurrence patterns of pathogenome (culturable and unculturable pathogens), antibiotic resistome (ARGs), and mobilome (mobile genetic elements (MGEs)) from mariculture waters near 5000 km coast of South China. Total 207 species of pathogens were identified, with only 10 culturable species. Furthermore, more pathogen species were detected in mariculture waters than those in coastal waters, and mariculture waters were prone to become reservoirs of unculturable pathogens. In addition, 913 subtypes of 21 ARG types were also identified, with multidrug resistance genes as the majority. MGEs including plasmids, integrons, transposons, and insertion sequences were abundantly present in mariculture waters. The co-occurrence network pattern between pathogenome, antibiotic resistome, and mobilome suggested that most of pathogens may be potential multidrug resistant hosts, possibly due to high frequency of horizontal gene transfer. These findings increase our understanding of mariculture waters as reservoirs of antibiotic resistome and mobilome, and as yet another hotbed for creation and transfer of new antibiotic-resistant pathogenome.

Exploring the resistome, virulome, and mobilome of multidrug-resistant Klebsiella pneumoniae isolates: deciphering the molecular basis of carbapenem resistance

BACKGROUND

Klebsiella pneumoniae, a notorious pathogen for causing nosocomial infections has become a major cause of neonatal septicemia, leading to high morbidity and mortality worldwide. This opportunistic bacterium has become highly resistant to antibiotics due to the widespread acquisition of genes encoding a variety of enzymes such as extended-spectrum beta-lactamases (ESBLs) and carbapenemases. We collected Klebsiella pneumoniae isolates from a local tertiary care hospital from February 2019-February 2021. To gain molecular insight into the resistome, virulome, and genetic environment of significant genes of multidrug-resistant K. pneumoniae isolates, we performed the short-read whole-genome sequencing of 10 K. pneumoniae isolates recovered from adult patients, neonates, and hospital tap water samples.

RESULTS

The draft genomes of the isolates varied in size, ranging from 5.48 to 5.96 Mbp suggesting the genome plasticity of this pathogen. Various genes conferring resistance to different classes of antibiotics e.g., aminoglycosides, quinolones, sulfonamides, tetracycline, and trimethoprim were identified in all sequenced isolates. The highest resistance was observed towards carbapenems, which has been putatively linked to the presence of both class B and class D carbapenemases, blaNDM, and blaOXA, respectively. Moreover, the biocide resistance gene qacEdelta1 was found in 6/10 of the sequenced strains. The sequenced isolates exhibited a broad range of sequence types and capsular types. The significant antibiotic resistance genes (ARGs) were bracketed by a variety of mobile genetic elements (MGEs). Various spontaneous mutations in genes other than the acquired antibiotic-resistance genes were observed, which play an indirect role in making these bugs resistant to antibiotics. Loss or deficiency of outer membrane porins, combined with ESBL production, played a significant role in carbapenem resistance in our sequenced isolates. Phylogenetic analysis revealed that the study isolates exhibited evolutionary relationships with strains from China, India, and the USA suggesting a shared evolutionary history and potential dissemination of similar genes amongst the isolates of different origins.

CONCLUSIONS

This study provides valuable insight into the presence of multiple mechanisms of carbapenem resistance in K. pneumoniae strains including the acquisition of multiple antibiotic-resistance genes through mobile genetic elements. Identification of rich mobilome yielded insightful information regarding the crucial role of insertion sequences, transposons, and integrons in shaping the genome of bacteria for the transmission of various resistance-associated genes. Multi-drug resistant isolates that had the fewest resistance genes exhibited a significant number of mutations. K. pneumoniae isolate from water source displayed comparable antibiotic resistance determinants to clinical isolates and the highest number of virulence-associated genes suggesting the possible interplay of ARGs amongst bacteria from different sources.

Long-read sequencing of extrachromosomal circular DNA and genome assembly of a Solanum lycopersicum breeding line revealed active LTR retrotransposons originating from S. Peruvianum L. introgressions

Transposable elements (TEs) are a major force in the evolution of plant genomes. Differences in the transposition activities and landscapes of TEs can vary substantially, even in closely related species. Interspecific hybridization, a widely employed technique in tomato breeding, results in the creation of novel combinations of TEs from distinct species. The implications of this process for TE transposition activity have not been studied in modern cultivars. In this study, we used nanopore sequencing of extrachromosomal circular DNA (eccDNA) and identified two highly active Ty1/Copia LTR retrotransposon families of tomato (Solanum lycopersicum), called Salsa and Ketchup. Elements of these families produce thousands of eccDNAs under controlled conditions and epigenetic stress. EccDNA sequence analysis revealed that the major parts of eccDNA produced by Ketchup and Salsa exhibited low similarity to the S. lycopersicum genomic sequence. To trace the origin of these TEs, whole-genome nanopore sequencing and de novo genome assembly were performed. We found that these TEs occurred in a tomato breeding line via interspecific introgression from S. peruvianum. Our findings collectively show that interspecific introgressions can contribute to both genetic and phenotypic diversity not only by introducing novel genetic variants, but also by importing active transposable elements from other species.

Comparative genomics of Staphylococcus aureus strains from wild birds and pig farms elucidates levels of mobilomes, antibiotic pressure and host adaptation

OBJECTIVES

This study characterized the resistome, mobilome and phylogenomic relatedness of Staphylococcus aureus strains previously obtained from healthy nestling storks (HNS), pigs (HP) and pig farmers (HPF) to analyse possible transmission pathways of S. aureus with implications for the spread of antimicrobial resistance.

METHODS

The genomic contents of 52 S. aureus strains obtained from the nasal cavity of HNS, HP and HPF in Spain were sequenced using the Illumina NextSeq platform to characterize their resistome, virulome and mobile genetic elements. The relatedness of strains was assessed by core-genome single nucleotide polymorphisms (SNPs).

RESULTS

The frequencies of multidrug-resistance phenotype and transposons were significantly lower in strains from HNS than in those from HP and HPF (P < 0.005). However, the presence of human immune evasion cluster genes in S. aureus strains from HNS was significantly higher than in those from HP and HPF (P < 0.005). Interestingly, the frequencies of plasmids and phages were not significantly associated with the host (P > 0.05). The phylogenetic analysis identified a cluster of all the MSSA-CC398 strains carrying φSa3 and ermT on rep13 separately from the two MRSA-CC398 strains (carrying ermT on repUS18). Highly related MRSA-CC398 strains were detected in some pigs and related farmers (<10 SNPs).

CONCLUSION

This study confirms high-level antibiotic selection in S. aureus in HP and HPF in comparison to HNS. Furthermore, our findings highlight the continuous transmission of MRSA-CC398 in the pig-to-human interface and MSSA-CC398 with human adaptation markers in HNS. Molecular surveillance of S. aureus using the One Health model is required to establish appropriate control strategies.

Mobile genetic elements affect the dissemination of antibiotic resistance genes (ARGs) of clinical importance in the environment

The prevalence of antibiotic resistance genes (ARGs) in the environment is a quintessential One Health issue that threats both human and ecosystem health; however, the source and transmission of ARGs, especially clinically important ARGs (CLIARGs), in the environment have not yet been well studied. In the present study, shotgun metagenomic approaches were used to characterize the microbiome, resistome, and mobilome composition in human feces and six different environment sample types in South China. Overall, the resistome harbored 157 CLIARGs, with specific ARG hotspots (e.g., human feces, wastewater treatment plants, livestock manure and wastewater) excreting significantly higher abundance of CLIARGs compared with the natural environment. A redundancy analysis (RDA) was performed and revealed that the bacterial community compositions and mobile genetic elements (MGEs) explained 55.08% and 34.68% of the variations in ARG abundance, respectively, indicating that both bacterial community and MGEs are key contributors to the maintenance and dissemination of CLIARGs in the environment. The network analysis revealed non-random co-occurrence patterns between 200 bacterial genera and 147 CLIARGs, as well as between 135 MGEs and 123 CLIARGs. In addition to numerous co-shared CLIARGs among different sample types, the source tracking program based on the FEAST probabilistic model was used to estimate the relative contributions of the CLIARGs from potential sources to the natural environment. The source tracking analysis results delineated that mobilome, more than microbiome, contributed CLIARG transmission from those ARG hotspots into natural environment, and the MGEs in WWTPs seem to play the most significant role in the spread of CLIARGs to the natural environment (average contribution 32.9%-46.4%). Overall, this study demonstrated the distribution and dissemination of CLIARGs in the environment, and aimed to better inform strategies to control the spread of CLIARGs into the natural environment.

Genome-resolved metagenomics of Venice Lagoon surface sediment bacteria reveals high biosynthetic potential and metabolic plasticity as successful strategies in an impacted environment

Bacteria living in sediments play essential roles in marine ecosystems and deeper insights into the ecology and biogeochemistry of these largely unexplored organisms can be obtained from 'omics' approaches. Here, we characterized metagenome-assembled-genomes (MAGs) from the surface sediment microbes of the Venice Lagoon (northern Adriatic Sea) in distinct sub-basins exposed to various natural and anthropogenic pressures. MAGs were explored for biodiversity, major marine metabolic processes, anthropogenic activity-related functions, adaptations at the microscale, and biosynthetic gene clusters. Starting from 126 MAGs, a non-redundant dataset of 58 was compiled, the majority of which (35) belonged to (Alpha- and Gamma-) Proteobacteria. Within the broad microbial metabolic repertoire (including C, N, and S metabolisms) the potential to live without oxygen emerged as one of the most important features. Mixotrophy was also found as a successful lifestyle. Cluster analysis showed that different MAGs encoded the same metabolic patterns (e.g., C fixation, sulfate oxidation) thus suggesting metabolic redundancy. Antibiotic and toxic compounds resistance genes were coupled, a condition that could promote the spreading of these genetic traits. MAGs showed a high biosynthetic potential related to antimicrobial and biotechnological classes and to organism defense and interactions as well as adaptive strategies for micronutrient uptake and cellular detoxification. Our results highlighted that bacteria living in an impacted environment, such as the surface sediments of the Venice Lagoon, may benefit from metabolic plasticity as well as from the synthesis of a wide array of secondary metabolites, promoting ecosystem resilience and stability toward environmental pressures.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00192-z.

Snapshot of resistome, virulome and mobilome in aquaculture

Aquaculture environments can be hotspots for resistance genes through the surrounding environment. Our objective was to study the resistome, virulome and mobilome of Gram-negative bacteria isolated in seabream and bivalve molluscs, using a WGS approach. Sixty-six Gram-negative strains (Aeromonadaceae, Enterobacteriaceae, Hafniaceae, Morganellaceae, Pseudomonadaceae, Shewanellaceae, Vibrionaceae, and Yersiniaceae families) were selected for genomic characterization. The species and MLST were determined, and antibiotic/disinfectants/heavy metals resistance genes, virulence determinants, MGE, and pathogenicity to humans were investigated. Our study revealed new sequence-types (e.g. Aeromonas spp. ST879, ST880, ST881, ST882, ST883, ST887, ST888; Shewanella spp. ST40, ST57, ST58, ST60, ST61, ST62; Vibrio spp. ST206, ST205). >140 different genes were identified in the resistome of seabream and bivalve molluscs, encompassing genes associated with β-lactams, tetracyclines, aminoglycosides, quinolones, sulfonamides, trimethoprim, phenicols, macrolides and fosfomycin resistance. Disinfectant resistance genes qacE-type, sitABCD-type and formA-type were found. Heavy metals resistance genes mdt, acr and sil stood out as the most frequent. Most resistance genes were associated with antibiotics/disinfectants/heavy metals commonly used in aquaculture settings. We also identified 25 different genes related with increased virulence, namely associated with adherence, colonization, toxins production, red blood cell lysis, iron metabolism, escape from the immune system of the host. Furthermore, 74.2 % of the strains analysed were considered pathogenic to humans. We investigated the genetic environment of several antibiotic resistance genes, including blaTEM-1B, blaFOX-18, aph(3″)-Ib, dfrA-type, aadA1, catA1-type, tet(A)/(E), qnrB19 and sul1/2. Our analysis also focused on identifying MGE in proximity to these genes (e.g. IntI1, plasmids and TnAs), which could potentially facilitate the spread of resistance among bacteria across different environments. This study provides a comprehensive examination of the diversity of resistance genes that can be transferred to both humans and the environment, with the recognition that aquaculture and the broader environment play crucial roles as intermediaries within this complex transmission network.

Key taxa and mobilome-mediated responses co-reshape the soil antibiotic resistome under dazomet fumigation stress

Agrochemicals are emergingly being implicated in the widespread dissemination of antibiotic resistance genes (ARGs) in agroecosystems. However, minimal research exists on the disturbance of fumigant on soil ARGs. Focusing on a typical fumigant dazomet in a simulated soil microcosm, we characterized the dazomet-triggered timely response and longstanding dynamic of ARGs at one-fold and two-fold field recommended doses using metagenome and quantitative PCR. Dazomet treatments reduced 13.17%-69.98% of absolute abundance of 16S rRNA gene and targeted ARGs, but, awfully, boosted diversity and relative abundance of ARGs up to 1.33-1.60 and 1.62-1.90 folds, respectively. Approximately 77.28% of changes in relative abundance of ARGs could be explained by bacterial community and mobile genetic elements (MGEs). Mechanistically, primary hosts of ARGs shifted from Proteobacteria (control) to Firmicutes and Actinobacteria (treatments) accompanied with corresponding changes in their abundance by combining community analysis, host tracking analysis and antibiotic resistant bacteria assay. Meanwhile, dazomet exposure significantly increased the incidence of MGEs and stimulated the conjugation of antibiotic-resistant plasmid. In addition, absolute abundance of targeted ARGs gradually recovered in the post-fumigation stage. Collectively, our results elucidate the dazomet-triggered emergence and spread of soil ARGs and highlight the importance of navigating toward rational use of fumigant in agricultural fields.

Pangenome analysis reveals the genetic basis for taxonomic classification of the Lactobacillaceae family

Lactobacillaceae represent a large family of important microbes that are foundational to the food industry. Many genome sequences of Lactobacillaceae strains are now available, enabling us to conduct a comprehensive pangenome analysis of this family. We collected 3591 high-quality genomes from public sources and found that: 1) they contained enough genomes for 26 species to perform a pangenomic analysis, 2) the normalized Heap's coefficient λ (a measure of pangenome openness) was found to have an average value of 0.27 (ranging from 0.07 to 0.37), 3) the pangenome openness was correlated with the abundance and genomic location of transposons and mobilomes, 4) the pangenome for each species was divided into core, accessory, and rare genomes, that highlight the species-specific properties (such as motility and restriction-modification systems), 5) the pangenome of Lactiplantibacillus plantarum (which contained the highest number of genomes found amongst the 26 species studied) contained nine distinct phylogroups, and 6) genome mining revealed a richness of detected biosynthetic gene clusters, with functions ranging from antimicrobial and probiotic to food preservation, but ∼93% were of unknown function. This study provides the first in-depth comparative pangenomics analysis of the Lactobacillaceae family.

Genomic characterization of Staphylococcus aureus from Canastra Minas Artisanal Cheeses

Canastra Minas Artisanal Cheese is produced in the Brazilian State of Minas Gerais using raw milk, rennet, and pingo, a natural endogenous starter culture (fermented whey) collected from the previous day's production. Due to the use of raw milk, the product can carry microorganisms that may cause foodborne diseases (FBD), including Staphylococcus aureus. Genomic characterization of S. aureus is an important tool to assess diversity, virulence, antimicrobial resistance, and the potential for causing food poisoning due to enterotoxin production. This study is aimed at exploring the genomic features of S. aureus strains isolated from Canastra Minas Artisanal Cheeses. Multilocus sequence typing (MLST) classified these strains as ST1, ST5, and a new profile ST7849 (assigned to the clonal complex CC97). These strains belonged to four spa types: t008, t127, t359, and t992. We identified antimicrobial resistance genes with phenotypic correlation against methicillin (MRSA) and tetracycline. Virulome analysis revealed genes associated with iron uptake, immune evasion, and potential capacity for adherence and biofilm formation. The toxigenic potential included cyto- and exotoxins genes, and all strains presented the genes that encode for Panton-Valentine toxin and hemolysin, and two strains encoded 4 and 8 Staphylococcal enterotoxin (SE) genes. The results revealed the pathogenic potential of the evaluated S. aureus strains circulating in the Canastra region, representing a potential risk to public health. This study also provides useful information to monitor and guide the application of control measures to the artisanal dairy food production chain.

Metagenomic evidence for antibiotics-driven co-evolution of microbial community, resistome and mobilome in hospital sewage

Overconsumption of antibiotics is an immediate cause for the emergence of antimicrobial resistance (AMR) and antibiotic resistant bacteria (ARB), though its environmental impact remains inadequately clarified. There is an urgent need to dissect the complex links underpinning the dynamic co-evolution of ARB and their resistome and mobilome in hospital sewage. Metagenomic and bioinformatic methods were employed to analyze the microbial community, resistome and mobilome in hospital sewage, in relation to data on clinical antibiotic use collected from a tertiary-care hospital. In this study, resistome (1,568 antibiotic resistance genes, ARGs, corresponding to 29 antibiotic types/subtypes) and mobilome (247 types of mobile genetic elements, MGEs) were identified. Networks connecting co-occurring ARGs with MGEs encompass 176 nodes and 578 edges, in which over 19 types of ARGs had significant correlations with MGEs. Prescribed dosage and time-dependent antibiotic consumption were associated with the abundance and distributions of ARGs, and conjugative transfer of ARGs via MGEs. Variation partitioning analyses show that effects of conjugative transfer were most likely the main contributors to transient propagation and persistence of AMR. We have presented the first evidence supporting idea that use of clinical antibiotics is a potent driving force for the development of co-evolving resistome and mobilome, which in turn supports the growth and evolution of ARB in hospital sewage. The use of clinical antibiotics calls for greater attention in antibiotic stewardship and management.

Microbiome, resistome and mobilome of chlorine-free drinking water treatment systems

Drinking water treatment plants (DWTPs) are designed to remove physical, chemical, and biological contaminants. However, until recently, the role of DWTPs in minimizing the cycling of antibiotic resistance determinants has got limited attention. In particular, the risk of selecting antibiotic-resistant bacteria (ARB) is largely overlooked in chlorine-free DWTPs where biological processes are applied. Here, we combined high-throughput quantitative PCR and metagenomics to analyze the abundance and dynamics of microbial communities, antibiotic resistance genes (ARGs), and mobile genetic elements (MGEs) across the treatment trains of two chlorine-free DWTPs involving dune-based and reservoir-based systems. The microbial diversity of the water increased after all biological unit operations, namely rapid and slow sand filtration (SSF), and granular activated carbon filtration. Both DWTPs reduced the concentration of ARGs and MGEs in the water by circa 2.5 log gene copies mL^-1, despite their relative increase in the disinfection sub-units (SSF in dune-based and UV treatment in reservoir-based DWTPs). The total microbial concentration was also reduced (2.5 log units), and none of the DWTPs enriched for bacteria containing genes linked to antibiotic resistance. Our findings highlight the effectiveness of chlorine-free DWTPs in supplying safe drinking water while reducing the concentration of antibiotic resistance determinants. To the best of our knowledge, this is the first study that monitors the presence and dynamics of antibiotic resistance determinants in chlorine-free DWTPs.

Development of early life gut resistome and mobilome across gestational ages and microbiota-modifying treatments

BACKGROUND

Gestational age (GA) and associated level of gastrointestinal tract maturation are major factors driving the initial gut microbiota composition in preterm infants. Besides, compared to term infants, premature infants often receive antibiotics to treat infections and probiotics to restore optimal gut microbiota. How GA, antibiotics, and probiotics modulate the microbiota's core characteristics, gut resistome and mobilome, remains nascent.

METHODS

We analysed metagenomic data from a longitudinal observational study in six Norwegian neonatal intensive care units to describe the bacterial microbiota of infants of varying GA and receiving different treatments. The cohort consisted of probiotic-supplemented and antibiotic-exposed extremely preterm infants (n = 29), antibiotic-exposed very preterm (n = 25), antibiotic-unexposed very preterm (n = 8), and antibiotic-unexposed full-term (n = 10) infants. The stool samples were collected on days of life 7, 28, 120, and 365, and DNA extraction was followed by shotgun metagenome sequencing and bioinformatical analysis.

FINDINGS

The top predictors of microbiota maturation were hospitalisation length and GA. Probiotic administration rendered the gut microbiota and resistome of extremely preterm infants more alike to term infants on day 7 and ameliorated GA-driven loss of microbiota interconnectivity and stability. GA, hospitalisation, and both microbiota-modifying treatments (antibiotics and probiotics) contributed to an elevated carriage of mobile genetic elements in preterm infants compared to term controls. Finally, Escherichia coli was associated with the highest number of antibiotic-resistance genes, followed by Klebsiella pneumoniae and Klebsiella aerogenes.

INTERPRETATION

Prolonged hospitalisation, antibiotics, and probiotic intervention contribute to dynamic alterations in resistome and mobilome, gut microbiota characteristics relevant to infection risk.

FUNDING

Odd-Berg Group, Northern Norway Regional Health Authority.

Assessment of plasmids for relating the 2020 Salmonella enterica serovar Newport onion outbreak to farms implicated by the outbreak investigation

BACKGROUND

The Salmonella enterica serovar Newport red onion outbreak of 2020 was the largest foodborne outbreak of Salmonella in over a decade. The epidemiological investigation suggested two farms as the likely source of contamination. However, single nucleotide polymorphism (SNP) analysis of the whole genome sequencing data showed that none of the Salmonella isolates collected from the farm regions were linked to the clinical isolates-preventing the use of phylogenetics in source identification. Here, we explored an alternative method for analyzing the whole genome sequencing data driven by the hypothesis that if the outbreak strain had come from the farm regions, then the clinical isolates would disproportionately contain plasmids found in isolates from the farm regions due to horizontal transfer.

RESULTS

SNP analysis confirmed that the clinical isolates formed a single, nearly-clonal clade with evidence for ancestry in California going back a decade. The clinical clade had a large core genome (4,399 genes) and a large and sparsely distributed accessory genome (2,577 genes, at least 64% on plasmids). At least 20 plasmid types occurred in the clinical clade, more than were found in the literature for Salmonella Newport. A small number of plasmids, 14 from 13 clinical isolates and 17 from 8 farm isolates, were found to be highly similar (> 95% identical)-indicating they might be related by horizontal transfer. Phylogenetic analysis was unable to determine the geographic origin, isolation source, or time of transfer of the plasmids, likely due to their promiscuous and transient nature. However, our resampling analysis suggested that observing a similar number and combination of highly similar plasmids in random samples of environmental Salmonella enterica within the NCBI Pathogen Detection database was unlikely, supporting a connection between the outbreak strain and the farms implicated by the epidemiological investigation.

CONCLUSION

Horizontally transferred plasmids provided evidence for a connection between clinical isolates and the farms implicated as the source of the outbreak. Our case study suggests that such analyses might add a new dimension to source tracking investigations, but highlights the need for detailed and accurate metadata, more extensive environmental sampling, and a better understanding of plasmid molecular evolution.

Metagenomic binning and assembled genome analysis revealed the distinct composition of resistome and mobilome in the Ili River

Rivers play an important role in receiving and transporting the resistome among different environmental compartments. However, the difference in resistome and mobilome between the water and sediment and their underlying mechanisms were still poorly understood. In this study, the Ili River, an important water source in the arid area of Central Asia, was selected as the studied target. The comprehensive profile of resistome and mobilome and their host in water and sediment were studied based on metagenomic binning and assembled genome (MAG) analysis. The relative abundance of resistome and mobilome in sediment were 28.0 - 67.8 × /Gb and 46.5 - 121.1 × /Gb, respectively, which were significantly higher than those in water (23.1 - 52.8 ×/Gb and 25.3 - 67.7 ×/Gb). Multidrug and macrolides-lincosamides-streptogramin (MLS) resistance genes were the main ARG types in both water and sediment from relative abundance. Transposases dominated the relative abundance of mobilome, followed by insert elements and integrases. Strong correlations were found between the relative abundance of resistome and mobilome (r > 0.6 and p < 0.01) in both water and sediment, indicating the mobilome played an important role in the propagation of resistome in the Ili River. The main hosts for multidrug resistance genes via MAG analysis differed in water (Alphaproteobacteria and Gammaproteobacteria) and sediment (Gammaproteobacteria). Distinct compositions of resistome and mobilome existed between water and sediment in the Ili River. Specificity-occupancy analysis of the differential resistome and mobilome showed that occurrence frequencies and habitat selections of the differential ARGs shaped the resistome of water and sediment. In contrast, habitat was the main driver that shaped the mobilome in the Ili River.

Genomic and proteomic characterization of two strains of Shigella flexneri 2 isolated from infants' stool samples in Argentina

Background

Shigella specie is a globally important intestinal pathogen disseminated all over the world. In this study we analyzed the genome and the proteomic component of two Shigella flexneri 2a clinical isolates, collected from pediatric patients with gastroenteritis of the Northwest region of Argentina (NWA) in two periods of time, with four years of difference. Our goal was to determine putative changes at molecular levels occurred during these four years, that could explain the presence of this Shigella`s serovar as the prevalent pathogen in the population under study.

Results

As previously reported, our findings support the idea of Shigella has a conserved “core” genome, since comparative studies of CI133 and CI172 genomes performed against 80 genomes obtained from the NCBI database, showed that there is a large number of genes shared among all of them. However, we observed that CI133 and CI172 harbors a small number of strain-specific genes, several of them present in mobile genetic elements, supporting the hypothesis that these isolates were established in the population by horizontal acquisition of genes. These differences were also observed at proteomic level, where it was possible to detect the presence of certain secreted proteins in a culture medium that simulates the host environment.

Conclusion

Great similarities were observed between the CI133 and CI172 strains, confirming the high percentage of genes constituting the “core” genome of S. flexneri 2. However, numerous strain specific genes were also determined. The presence of the here identified molecular elements into other strain of our culture collation, is currently used to develop characteristic markers of local pathogens. In addition, the most outstanding result of this study was the first description of a S. flexneri 2 producing Colicin E, as one of the characteristics that allows S. flexneri 2 to persist in the microbial community. These findings could also contribute to clarify the mechanism and the evolution strategy used by this pathogen to specifically colonize, survive, and cause infection within the NWA population.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08711-5.

Spatiotemporal dynamics of the resistome and virulome of riverine microbiomes disturbed by a mining mud tsunami

Aquatic ecosystems are highly vulnerable to anthropogenic activities. However, it remains unclear how the microbiome responds to press disturbance events in these ecosystems. We examined the impact of the world's largest mining disaster (Brazil, 2015) on sediment microbiomes in two disturbed rivers compared to an undisturbed river during 390 days post-disturbance. The diversity and structure of the virulome and microbiome, and of antibiotic and metal resistomes, consistently differed between the disturbed and undisturbed rivers, particularly at day 7 post-disturbance. 684 different ARGs were predicted, 38% were exclusive to the disturbed rivers. Critical antibiotic resistance genes (ARGs), e.g., mcr and ereA2, were significantly more common in the disturbed microbiomes. 401 different ARGs were associated with mobile genetic elements (MGEs), 95% occurred in the disturbed rivers. While plasmids were the most common MGEs with a broad spectrum of ARGs, spanning 16 antibiotic classes, integrative conjugative elements (ICEs) and integrons disseminated ARGs associated with aminoglycoside and tetracycline, and aminoglycoside and beta-lactam, respectively. A significant increase in the relative abundance of class 1 integrons, ICEs, and pathogens was identified at day 7 in the disturbed microbiomes, 72-, 14- and 3- fold higher, respectively, compared with the undisturbed river. Mobile ARGs associated with ESKAPEE group pathogens, while metal resistance genes and virulence factor genes in nonpathogenic hosts predominated in all microbiomes. Network analysis showed highly interconnected ARGs in the disturbed communities, including genes targeting antibiotics of last resort. Interactions between copper and beta-lactam/aminoglycoside/macrolide resistance genes, mostly mobile and critical, were also uncovered. We conclude that the mud tsunami resulted in resistome expansion, enrichment of pathogens, and increases in promiscuous and mobile ARGs. From a One Health perspective, mining companies need to move toward more environmentally friendly and socially responsible mining practices to reduce risks associated with pathogens and critical and mobile ARGs.

Nanopore-based metagenomics analysis reveals prevalence of mobile antibiotic and heavy metal resistome in wastewater

In-depth studies of the microbiome and mobile resistome profile of different environments is central to understanding the role of the environment in antimicrobial resistance (AMR), which is one of the urgent threats to global public health. In this study, we demonstrated the use of a rapid (and easily portable) sequencing approach coupled with user-friendly bioinformatics tools, the MinION (Oxford Nanopore Technologies), on the evaluation of the microbial as well as mobile metal and antibiotic resistome profile of semi-rural wastewater. A total of 20 unique phyla, 43 classes, 227 genera, and 469 species were identified in samples collected from the Amherst Wastewater Treatment Plant, both from primary and secondary treated wastewater. Alpha diversity indices indicated that primary samples were significantly richer and more microbially diverse than secondary samples. A total of 1041 ARGs, 68 MRGs, and 17 MGEs were detected in this study. There were more classes of AMR genes in primary than secondary wastewater, but in both cases multidrug, beta-lactam and peptide AMR predominated. Of note, OXA β-lactamases, some of which are also carbapenemases, were enriched in secondary samples. Metal resistance genes against arsenic, copper, zinc and molybdenum were the dominant MRGs in the majority of the samples. A larger proportion of resistome genes were located in chromosome-derived sequences except for mobilome genes, which were predominantly located in plasmid-derived sequences. Genetic elements related to transposase were the most common MGEs in all samples. Mobile or MGE/plasmid-associated resistome genes that confer resistance to last resort antimicrobials such as carbapenems and colistin were detected in most samples. Worryingly, several of these potentially transferable genes were found to be carried by clinically-relevant hosts including pathogenic bacterial species in the orders Aeromonadales, Clostridiales, Enterobacterales and Pseudomonadales. This study demonstrated that the MinION can be used as a metagenomics approach to evaluate the microbiome, resistome, and mobilome profile of primary and secondary wastewater.

Metagenomic exploration of antibiotic resistome in treated wastewater effluents and their receiving water

Environmental dissemination of antimicrobial resistance is a global health problem. Antimicrobial-resistant bacteria and antibiotic-resistant genes (ARGs) are constantly released into the environment through effluents (EFs) from wastewater treatment plants (WWTPs). Thus, requiring a better understanding of the selection and fate of ARGs in wastewater treatment processes. Therefore, we investigated the impacts of urban WWTP EFs on receiving water in the context of their resistomes and mobilomes. We used a HiSeq-based short read metagenomic approach to address the dynamics and diversity of ARGs in WWTP EF as well as the upstream (UP) and downstream (DN) river waters, followed by an investigation of plasmid-mediated ARGs. The abundance of ARGs at each site varied from 7.2 × 10^-2 to 7.4 × 10^-1 ARG copies per 16S rRNA gene copy, and EF samples showed the highest abundance, followed by DN and UP water samples. ARG diversity ranged from 121 to 686 types per site, and EF had the most diverse ARGs. Commonly identified ARGs in the EF and DN samples were clinically important and were absent in UP samples. The abundance of ARGs, mobile genetic elements (MGEs), and plasmid contigs found only in EF and DN were positively correlated with each other, indicating the importance of mobilomes in the dissemination of ARGs in the environment. Moreover, the proportions of plasmid-mediated ARGs was highest in the EF samples, followed by the DN and UP samples. These findings suggest that WWTP EF may act as a driving factor shaping the resistomes and mobilomes of receiving waters. In particular, a higher abundance of plasmid-mediated ARGs in WWTP EF suggests higher transmissibility in the DN environment.

Comparative genomics of the ADA clade within the Nostocales

The ADA clade of Nostocales cyanobacteria, a group that is prominent in current harmful algal bloom events, now includes over 40 genome sequences with the recent addition of sixteen novel sequenced genomes (Dreher et al., Harmful Algae, 2021). Fourteen genomes are complete (closed), enabling highly detailed assessments of gene content and genome architecture. ADA genomes contain 5 rRNA operons, genes expected to support a photoautotrophic and diazotrophic lifestyle, and a varied array of genes for the synthesis of bioactive secondary metabolites. Genes for the production of the taste-and-odor compound geosmin and the four major classes of cyanotoxins - anatoxin-a, cylindrospermopsin, microcystin and saxitoxin - are represented in members of the ADA clade. Notably, the gene array for the synthesis of cylindrospermopsin by Dolichospermum sp. DET69 was located on a plasmid, raising the possibility of facile horizontal transmission. However, genes supporting independent conjugative transfer of this plasmid are lacking. Further, analysis of genomic loci containing this and other cyanotoxin gene arrays shows evidence that these arrays have long-term stability and do not appear to be genomic islands easily capable of horizontal transmission to other cells. There is considerable diversity in the gene complements of individual ADA genomes, including the variable presence of physiologically important genes: genomes in three species-level subclades lack the gas vesicle genes that facilitate a planktonic lifestyle, and, surprisingly, the genome of Cuspidothrix issatschenkoi CHARLIE-1, a reported diazotroph, lacks the genes for nitrogen fixation. Notably, phylogenetically related genomes possess limited synteny, indicating a prominent role for chromosome rearrangements during ADA strain evolution. The genomes contain abundant insertion sequences and repetitive transposase genes, which could be the main drivers of genome rearrangement through active transposition and homologous recombination. No prophages were found, and no evidence of viral infection was observed in the bloom population samples from which the genomes discussed here were derived. Phages thus seem to have a limited influence on ADA evolution.

Approaches for characterizing and tracking hospital-associated multidrug-resistant bacteria

Hospital-associated infections are a major concern for global public health. Infections with antibiotic-resistant pathogens can cause empiric treatment failure, and for infections with multidrug-resistant bacteria which can overcome antibiotics of "last resort" there exists no alternative treatments. Despite extensive sanitization protocols, the hospital environment is a potent reservoir and vector of antibiotic-resistant organisms. Pathogens can persist on hospital surfaces and plumbing for months to years, acquire new antibiotic resistance genes by horizontal gene transfer, and initiate outbreaks of hospital-associated infections by spreading to patients via healthcare workers and visitors. Advancements in next-generation sequencing of bacterial genomes and metagenomes have expanded our ability to (1) identify species and track distinct strains, (2) comprehensively profile antibiotic resistance genes, and (3) resolve the mobile elements that facilitate intra- and intercellular gene transfer. This information can, in turn, be used to characterize the population dynamics of hospital-associated microbiota, track outbreaks to their environmental reservoirs, and inform future interventions. This review provides a detailed overview of the approaches and bioinformatic tools available to study isolates and metagenomes of hospital-associated bacteria, and their multi-layered networks of transmission.

Complete Genomes of the Anaerobic Degradation Specialists Aromatoleum petrolei ToN1T and Aromatoleum bremense PbN1T

The betaproteobacterial genus Aromatoleum comprises facultative denitrifiers specialized in the anaerobic degradation of recalcitrant organic compounds (aromatic and terpenoid). This study reports on the complete and manually annotated genomes of Ar. petrolei ToN1T (5.41 Mbp) and Ar. bremense PbN1T (4.38 Mbp), which cover the phylogenetic breadth of the genus Aromatoleum together with previously genome sequenced Ar. aromaticum EbN1T [Rabus et al., Arch Microbiol. 2005 Jan;183(1):27-36]. The gene clusters for the anaerobic degradation of aromatic and terpenoid (strain ToN1T only) compounds are scattered across the genomes of strains ToN1T and PbN1T. The richness in mobile genetic elements is shared with other Aromatoleum spp., substantiating that horizontal gene transfer should have been a major driver in shaping the genomes of this genus. The composite catabolic network of strains ToN1T and PbN1T comprises 88 proteins, the coding genes of which occupy 86.1 and 76.4 kbp (1.59 and 1.75%) of the respective genome. The strain-specific gene clusters for anaerobic degradation of ethyl-/propylbenzene (strain PbN1T) and toluene/monoterpenes (strain ToN1T) share high similarity with their counterparts in Ar. aromaticum strains EbN1T and pCyN1, respectively. Glucose is degraded via the ED-pathway in strain ToN1T, while gluconeogenesis proceeds via the reverse EMP-pathway in strains ToN1T, PbN1T, and EbN1T. The diazotrophic, endophytic lifestyle of closest related genus Azoarcus is known to be associated with nitrogenase and type-6 secretion system (T6SS). By contrast, strains ToN1T, PbN1T, and EbN1T lack nif genes for nitrogenase (including cofactor synthesis and enzyme maturation). Moreover, strains PbN1T and EbN1T do not possess tss genes for T6SS, while strain ToN1T does and facultative endophytic "Aromatoleum" sp. CIB is known to even have both. These findings underpin the functional heterogeneity among Aromatoleum members, correlating with the high plasticity of their genomes.

Antimicrobial use and production system shape the fecal, environmental, and slurry resistomes of pig farms

BACKGROUND

The global threat of antimicrobial resistance (AMR) is a One Health problem impacted by antimicrobial use (AMU) for human and livestock applications. Extensive Iberian swine production is based on a more sustainable and eco-friendly management system, providing an excellent opportunity to evaluate how sustained differences in AMU impact the resistome, not only in the animals but also on the farm environment. Here, we evaluate the resistome footprint of an extensive pig farming system, maintained for decades, as compared to that of industrialized intensive pig farming by analyzing 105 fecal, environmental and slurry metagenomes from 38 farms.

RESULTS

Our results evidence a significantly higher abundance of antimicrobial resistance genes (ARGs) on intensive farms and a link between AMU and AMR to certain antimicrobial classes. We observed differences in the resistome across sample types, with a higher richness and dispersion of ARGs within environmental samples than on those from feces or slurry. Indeed, a deeper analysis revealed that differences among the three sample types were defined by taxa-ARGs associations. Interestingly, mobilome analyses revealed that the observed AMR differences between intensive and extensive farms could be linked to differences in the abundance of mobile genetic elements (MGEs). Thus, while there were no differences in the abundance of chromosomal-associated ARGs between intensive and extensive herds, a significantly higher abundance of integrons in the environment and plasmids, regardless of the sample type, was detected on intensive farms.

CONCLUSIONS

Overall, this study shows how AMU, production system, and sample type influence, mainly through MGEs, the profile and dispersion of ARGs in pig production. Video Abstract.

Dynamics of class 1 integrons in aerobic biofilm reactors spiked with antibiotics

Class 1 integrons are strongly associated with the dissemination of antibiotic resistance in bacteria. However, little is known about whether the presence of antibiotics affects the abundance of integrons and antibiotic resistance genes during biological wastewater treatment. To explore the roles of class 1 integrons in spreading antibiotic resistance genes in environmental compartments, the dynamics of integrons were followed in biofilm reactors treating synthetic wastewater respectively spiked with streptomycin (STM) and oxytetracycline (OTC). The relative abundance of the integron-integrase gene (intI1) increased 12 or 29-fold respectively when treated with STM or OTC, under incrementally increasing dosage regimes from 0 to 50 mg L^-1. Significant increases in intI1 abundance initially occurred at an antibiotic dose of 0.1 mg L^-1. At the beginning of the experiment, 51% to 64% of integrons carried no gene cassettes. In STM and OTC spiked systems, there was a significant increase in the proportion of integrons that contained resistance gene cassettes, particularly at intermediate and higher antibiotic concentrations. Gene cassettes encoding resistance to aminoglycosides, trimethoprim, beta-lactam, erythromycin, and quaternary ammonium compounds were all detected in the treated systems. Three tetracycline resistance genes (tetA, tetC, tetG) were significantly correlated with the abundance of intI1 (p < 0.01), despite no tet resistance being present as a gene cassette. Genome sequencing of isolates showed synteny between the tet resistance genes and intI1, mediated through linkage to transposable elements including Tn3, IS26 and ISCR3. Class 1 integrons appeared to be under positive selection in the presence of antibiotics, and might have actively acquired new gene cassettes during the experiment.

Diverse mobilome of Dichotomius (Luederwaldtinia) schiffleri (Coleoptera: Scarabaeidae) reveals long-range horizontal transfer events of DNA transposons

Transposable elements (TEs) are mobile DNA sequences that are able to move from one genomic location to another. These selfish elements are known as genomic parasites, since they hijack the host molecular machinery to generate new copies of themselves. The mobilization of TEs can be seen as a natural mutagen because new TE copies can insert into different loci and impact host genomic structure through different mechanisms. Although our knowledge about TEs is improving with new genomes available, there is still very limited data about the mobilome of species from the Coleoptera order, the most diverse order of insects, including species from the Scarabaeidae family. Therefore, the main goal of this study was to characterize the mobilome of D. (Luederwaldtinia) schiffleri, based on low-coverage genome sequencing, and reconstruct their evolutionary history. We used a combination of four different approaches for TE characterization and maximum likelihood phylogenetic analysis to study their evolution. We found a large and diverse mobilome composed of 38 TE superfamilies, 20 DNA transposon and 18 retrotransposons, accounting for 21% of the genome. Moreover, we found a number of incongruences between the TE and host phylogenetic trees in three DNA transposon TE superfamilies, which represents five TE families, suggesting possible horizontal transfer events between highly divergent taxa. In summary, we found an abundant and diverse mobilome and a number of horizontal transfer events that have shaped the evolutionary history of this species.

Pathogenomic Analysis of a Novel Extensively Drug-Resistant Citrobacter freundii Isolate Carrying a blaNDM-1 Carbapenemase in South Africa

Pathogenomic analysis was performed on a novel carbapenem-resistant Citrobacter freundii isolate (H2730R) from a rectal swab of an adult male patient admitted to a tertiary hospital, Durban, South Africa. H2730R was identified using selective media and API 20e kit. Confirmatory identification and antibiotic susceptibility testing were performed using the VITEK II. H2730R was whole-genome sequenced on the Illumina MiSeq platform. H2730R was resistant to all tested antibiotics except tigecycline and was defined as ST498 by the C. freundii multilocus sequence typing (MLST) database. The estimated pathogenic potential predicted a higher probability (P_score ≈ 0.875), supporting H2730R as a human pathogen. H2730R harbored 25 putative acquired resistance genes, 4 plasmid replicons, 4 intact prophages, a class 1 integron (IntI1), 2 predominant insertion sequences (IS3 and IS5), numerous efflux genes, and virulome. BLASTn analysis of the bla_NDM-1 encoding contig (00022) and its flanking sequences revealed the bla_NDM-1 was located on a plasmid similar to the multireplicon p18-43_01 plasmid reported for the spread of carbapenem resistance in South Africa. Phylogenomic analysis showed clustering of H2730R with CF003/CF004 strains in the same clade, suggesting a possible association between C. freundii strains/clones. Acquiring the p18-43_01 plasmid containing bla_NDM-1, the diversity, and complex resistome, virulome, and mobilome of this pathogen makes its incidence very worrying regarding mobilized resistance. This study presents the background genomic information for future surveillance and tracking of the spread of carbapenem-resistant Enterobacteriaceae in South Africa.

Comparative genomics of multidrug-resistant Enterococcus spp. isolated from wastewater treatment plants

BACKGROUND

Wastewater treatment plants (WWTPs) are considered hotspots for the environmental dissemination of antimicrobial resistance (AMR) determinants. Vancomycin-Resistant Enterococcus (VRE) are candidates for gauging the degree of AMR bacteria in wastewater. Enterococcus faecalis and Enterococcus faecium are recognized indicators of fecal contamination in water. Comparative genomics of enterococci isolated from conventional activated sludge (CAS) and biological aerated filter (BAF) WWTPs was conducted.

RESULTS

VRE isolates, including E. faecalis (n = 24), E. faecium (n = 11), E. casseliflavus (n = 2) and E. gallinarum (n = 2) were selected for sequencing based on WWTP source, species and AMR phenotype. The pangenomes of E. faecium and E. faecalis were both open. The genomic fraction related to the mobilome was positively correlated with genome size in E. faecium (p < 0.001) and E. faecalis (p < 0.001) and with the number of AMR genes in E. faecium (p = 0.005). Genes conferring vancomycin resistance, including vanA and vanM (E. faecium), vanG (E. faecalis), and vanC (E. casseliflavus/E. gallinarum), were detected in 20 genomes. The most prominent functional AMR genes were efflux pumps and transporters. A minimum of 16, 6, 5 and 3 virulence genes were detected in E. faecium, E. faecalis, E. casseliflavus and E. gallinarum, respectively. Virulence genes were more common in E. faecalis and E. faecium, than E. casseliflavus and E. gallinarum. A number of mobile genetic elements were shared among species. Functional CRISPR/Cas arrays were detected in 13 E. faecalis genomes, with all but one also containing a prophage. The lack of a functional CRISPR/Cas arrays was associated with multi-drug resistance in E. faecium. Phylogenetic analysis demonstrated differential clustering of isolates based on original source but not WWTP. Genes related to phage and CRISPR/Cas arrays could potentially serve as environmental biomarkers.

CONCLUSIONS

There was no discernible difference between enterococcal genomes from the CAS and BAF WWTPs. E. faecalis and E. faecium have smaller genomes and harbor more virulence, AMR, and mobile genetic elements than other Enterococcus spp.

The Mobilome: Metagenomic Analysis of Circular Plasmids, Viruses, and Other Extrachromosomal Elements

Isolation, sequencing, and analysis of circular genetic elements bring new insights to mobile genetic elements related to microbial ecology. One method used to study circular plasmids, viruses, and other elements is called the mobilome method. The mobilome method presented here is an unamplified mobilome approach allowing fast isolation of circular DNA elements from a variety of samples followed by directly building unamplified Illumina-compatible sequencing libraries using enzymatic tagging and fragmentation. Several methods for bioinformatic analysis of mobilome data are also suggested.

LTR-TEs abundance, timing and mobility in Solanum commersonii and S. tuberosum genomes following cold-stress conditions

Copia/Ale is the youngest lineage in both Solanum tuberosum and S. commersonii. Within it, we identified nightshade, a new LTR element active in the cultivated potato. From an evolutionary perspective, long-terminal repeat retrotransposons (LTR-RT) activity during stress may be viewed as a mean by which organisms can keep up rates of genetic adaptation to changing conditions. Potato is one of the most important crop consumed worldwide, but studies on LTR-RT characterization are still lacking. Here, we assessed the abundance, insertion time and activity of LTR-RTs in both cultivated Solanum tuberosum and its cold-tolerant wild relative S. commersonii genomes. Gypsy elements were more abundant than Copia ones, suggesting that the former was somehow more successful in colonizing potato genomes. However, Copia elements, and in particular, the Ale lineage, are younger than Gypsy ones, since their insertion time was in average ~ 2 Mya. Due to the ability of LTR-RTs to be circularized by the host DNA repair mechanisms, we identified via mobilome-seq a Copia/Ale element (called nightshade, informal name used for potato family) active in S. tuberosum genome. Our analyses represent a valuable resource for comparative genomics within the Solanaceae, transposon-tagging and for the design of cultivar-specific molecular markers in potato.

The complete and fully assembled genome sequence of Aeromonas salmonicida subsp. pectinolytica and its comparative analysis with other Aeromonas species: investigation of the mobilome in environmental and pathogenic strains

BACKGROUND

Due to the predominant usage of short-read sequencing to date, most bacterial genome sequences reported in the last years remain at the draft level. This precludes certain types of analyses, such as the in-depth analysis of genome plasticity.

RESULTS

Here we report the finalized genome sequence of the environmental strain Aeromonas salmonicida subsp. pectinolytica 34mel, for which only a draft genome with 253 contigs is currently available. Successful completion of the transposon-rich genome critically depended on the PacBio long read sequencing technology. Using finalized genome sequences of A. salmonicida subsp. pectinolytica and other Aeromonads, we report the detailed analysis of the transposon composition of these bacterial species. Mobilome evolution is exemplified by a complex transposon, which has shifted from pathogenicity-related to environmental-related gene content in A. salmonicida subsp. pectinolytica 34mel.

CONCLUSION

Obtaining the complete, circular genome of A. salmonicida subsp. pectinolytica allowed us to perform an in-depth analysis of its mobilome. We demonstrate the mobilome-dependent evolution of this strain's genetic profile from pathogenic to environmental.

Discovery of rare, diagnostic AluYb8/9 elements in diverse human populations

BACKGROUND

Polymorphic human Alu elements are excellent tools for assessing population structure, and new retrotransposition events can contribute to disease. Next-generation sequencing has greatly increased the potential to discover Alu elements in human populations, and various sequencing and bioinformatics methods have been designed to tackle the problem of detecting these highly repetitive elements. However, current techniques for Alu discovery may miss rare, polymorphic Alu elements. Combining multiple discovery approaches may provide a better profile of the polymorphic Alu mobilome. AluYb8/9 elements have been a focus of our recent studies as they are young subfamilies (~2.3 million years old) that contribute ~30% of recent polymorphic Alu retrotransposition events. Here, we update our ME-Scan methods for detecting Alu elements and apply these methods to discover new insertions in a large set of individuals with diverse ancestral backgrounds.

RESULTS

We identified 5,288 putative Alu insertion events, including several hundred novel AluYb8/9 elements from 213 individuals from 18 diverse human populations. Hundreds of these loci were specific to continental populations, and 23 non-reference population-specific loci were validated by PCR. We provide high-quality sequence information for 68 rare AluYb8/9 elements, of which 11 have hallmarks of an active source element. Our subfamily distribution of rare AluYb8/9 elements is consistent with previous datasets, and may be representative of rare loci. We also find that while ME-Scan and low-coverage, whole-genome sequencing (WGS) detect different Alu elements in 41 1000 Genomes individuals, the two methods yield similar population structure results.

CONCLUSION

Current in-silico methods for Alu discovery may miss rare, polymorphic Alu elements. Therefore, using multiple techniques can provide a more accurate profile of Alu elements in individuals and populations. We improved our false-negative rate as an indicator of sample quality for future ME-Scan experiments. In conclusion, we demonstrate that ME-Scan is a good supplement for next-generation sequencing methods and is well-suited for population-level analyses.

Insertion and deletion polymorphisms of the ancient AluS family in the human genome

Background

Polymorphic Alu elements account for 17% of structural variants in the human genome. The majority of these belong to the youngest AluY subfamilies, and most structural variant discovery efforts have focused on identifying Alu polymorphisms from these currently retrotranspositionally active subfamilies. In this report we analyze polymorphisms from the evolutionarily older AluS subfamily, whose peak activity was tens of millions of years ago. We annotate the AluS polymorphisms, assess their likely mechanism of origin, and evaluate their contribution to structural variation in the human genome.

Results

Of 52 previously reported polymorphic AluS elements ascertained for this study, 48 were confirmed to belong to the AluS subfamily using high stringency subfamily classification criteria. Of these, the majority (77%, 37/48) appear to be deletion polymorphisms. Two polymorphic AluS elements (4%) have features of non-classical Alu insertions and one polymorphic AluS element (2%) likely inserted by a mechanism involving internal priming. Seven AluS polymorphisms (15%) appear to have arisen by the classical target-primed reverse transcription (TPRT) retrotransposition mechanism. These seven TPRT products are 3′ intact with 3′ poly-A tails, and are flanked by target site duplications; L1 ORF2p endonuclease cleavage sites were also observed, providing additional evidence that these are L1 ORF2p endonuclease-mediated TPRT insertions. Further sequence analysis showed strong conservation of both the RNA polymerase III promoter and SRP9/14 binding sites, important for mediating transcription and interaction with retrotransposition machinery, respectively. This conservation of functional features implies that some of these are fairly recent insertions since they have not diverged significantly from their respective retrotranspositionally competent source elements.

Conclusions

Of the polymorphic AluS elements evaluated in this report, 15% (7/48) have features consistent with TPRT-mediated insertion, thus suggesting that some AluS elements have been more active recently than previously thought, or that fixation of AluS insertion alleles remains incomplete. These data expand the potential significance of polymorphic AluS elements in contributing to structural variation in the human genome. Future discovery efforts focusing on polymorphic AluS elements are likely to identify more such polymorphisms, and approaches tailored to identify deletion alleles may be warranted.

Electronic supplementary material

The online version of this article (doi:10.1186/s13100-017-0089-9) contains supplementary material, which is available to authorized users.

Changes of resistome, mobilome and potential hosts of antibiotic resistance genes during the transformation of anaerobic digestion from mesophilic to thermophilic

This study aimed to reveal how antibiotic resistance genes (ARGs) and their horizontal and vertical transfer-related items (mobilome and bacterial hosts) respond to the transformation of anaerobic digestion (AD) from mesophilic to thermophilic using one-step temperature increase. The resistomes and mobilomes of mesophilic and thermophilic sludge were investigated using metagenome sequencing, and the changes in 24 representative ARGs belonging to three categories, class 1 integron and bacterial genera during the transition period were further followed using quantitative PCR and 454-pyrosequencing. After the temperature increase, resistome abundance in the digested sludge decreased from 125.97 ppm (day 0, mesophilic) to 50.65 ppm (day 57, thermophilic) with the reduction of most ARG types except for the aminoglycoside resistance genes. Thermophilic sludge also had a smaller mobilome, including plasmids, insertion sequences and integrons, than that of mesophilic sludge, suggesting the lower horizontal transfer potential of ARGs under thermophilic conditions. On the other hand, the total abundance of 18 bacterial genera, which were suggested as the possible hosts for 13 ARGs through network analysis, decreased from 23.27% in mesophilic sludge to 11.92% in thermophilic sludge, indicating fewer hosts for the vertical expansion of ARGs after the increase in temperature. These results indicate that the better reduction of resistome abundance by thermophilic AD might be associated with the decrease of both the horizontal and vertical transferability of ARGs.

Structural and functional analysis of the finished genome of the recently isolated toxic Anabaena sp. WA102

Background

Very few closed genomes of the cyanobacteria that commonly produce toxic blooms in lakes and reservoirs are available, limiting our understanding of the properties of these organisms. A new anatoxin-a-producing member of the Nostocaceae, Anabaena sp. WA102, was isolated from a freshwater lake in Washington State, USA, in 2013 and maintained in non-axenic culture.

Results

The Anabaena sp. WA102 5.7 Mbp genome assembly has been closed with long-read, single-molecule sequencing and separately a draft genome assembly has been produced with short-read sequencing technology. The closed and draft genome assemblies are compared, showing a correlation between long repeats in the genome and the many gaps in the short-read assembly. Anabaena sp. WA102 encodes anatoxin-a biosynthetic genes, as does its close relative Anabaena sp. AL93 (also introduced in this study). These strains are distinguished by differences in the genes for light-harvesting phycobilins, with Anabaena sp. AL93 possessing a phycoerythrocyanin operon. Biologically relevant structural variants in the Anabaena sp. WA102 genome were detected only by long-read sequencing: a tandem triplication of the anaBCD promoter region in the anatoxin-a synthase gene cluster (not triplicated in Anabaena sp. AL93) and a 5-kbp deletion variant present in two-thirds of the population. The genome has a large number of mobile elements (160). Strikingly, there was no synteny with the genome of its nearest fully assembled relative, Anabaena sp. 90.

Conclusion

Structural and functional genome analyses indicate that Anabaena sp. WA102 has a flexible genome. Genome closure, which can be readily achieved with long-read sequencing, reveals large scale (e.g., gene order) and local structural features that should be considered in understanding genome evolution and function.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2738-7) contains supplementary material, which is available to authorized users.

Current perspectives on the dynamics of antibiotic resistance in different reservoirs

Antibiotic resistance consists of a dynamic web. In this review, we describe the path by which different antibiotic residues and antibiotic resistance genes disseminate among relevant reservoirs (human, animal, and environmental settings), evaluating how these events contribute to the current scenario of antibiotic resistance. The relationship between the spread of resistance and the contribution of different genetic elements and events is revisited, exploring examples of the processes by which successful mobile resistance genes spread across different niches. The importance of classic and next generation molecular approaches, as well as action plans and policies which might aid in the fight against antibiotic resistance, are also reviewed.

Genetic variability of psychrotolerant Acidithiobacillus ferrivorans revealed by (meta)genomic analysis

Acidophilic microorganisms inhabit low pH environments such as acid mine drainage that is generated when sulfide minerals are exposed to air. The genome sequence of the psychrotolerant Acidithiobacillus ferrivorans SS3 was compared to a metagenome from a low temperature acidic stream dominated by an A. ferrivorans-like strain. Stretches of genomic DNA characterized by few matches to the metagenome, termed 'metagenomic islands', encoded genes associated with metal efflux and pH homeostasis. The metagenomic islands were enriched in mobile elements such as phage proteins, transposases, integrases and in one case, predicted to be flanked by truncated tRNAs. Cus gene clusters predicted to be involved in copper efflux and further Cus-like RND systems were predicted to be located in metagenomic islands and therefore, constitute part of the flexible gene complement of the species. Phylogenetic analysis of Cus clusters showed both lineage specificity within the Acidithiobacillus genus as well as niche specificity associated with an acidic environment. The metagenomic islands also contained a predicted copper efflux P-type ATPase system and a polyphosphate kinase potentially involved in polyphosphate mediated copper resistance. This study identifies genetic variability of low temperature acidophiles that likely reflects metal resistance selective pressures in the copper rich environment.