I have had the pleasure to be chosen as a speaker for next week’s (ten days from now) Swedish Bioinformatics Workshop. My talk is entitled “Turn up the signal – wipe out the noise: Gaining insights into bacterial community functions using metagenomic data“, and will largely deal with the same questions as my talk on EDAR3 in May this year. As then, the talk will highlight the some particular pitfalls related to interpretation of data, and exemplify how flawed analysis practices can result in misleading conclusions regarding community function, and use examples from our studies of environments subjected to pharmaceutical pollution in India, the effect of travel on the human resistome, and modern municipal wastewater treatment processes.
The talk will take place on Thursday, September 24, 2015 at 16:30. The full program for the conference can be found here. And also, if you want a sneak peak of the talk, you can drop by on Friday 13.00 at Chemistry and Molecular Biology, where I will give a seminar on the same topic in the Monthly Bioinformatic Practical Meetings series.
Earlier today, my most recent paper (1) became available online, describing resistance gene patterns in the gut microbiota of Swedes before and after travel to the Indian peninsula and central Africa. In this work, we have used metagenomic sequencing of the intestinal microbiome of Swedish students returning from exchange programs to show that the abundance of antibiotic resistance genes in several classes are increased after travel. This work reiterates the findings of several papers describing uptake of resistant bacteria (2-8) or resistance genes (9-11) after travel to destinations with worse resistance situation.
Our paper is important because it:
- Addresses the abundance of a vast range of resistance genes (more than 300).
- Finds evidence for that the overall relative abundance of antibiotic resistance genes increased after travel, without any intake of antibiotics.
- Shows that the sensitivity of metagenomics was, despite very deep sequencing efforts, not sufficient to detect acquisition of the low-abundant (CTX-M) resistance genes responsible for observed ESBL phenotypes.
- Reveals a “core resistome” of resistance genes that are more or less omnipresent, and remain relatively stable regardless of travel, while changes seem to occur in the more variable part of the resistome.
- Hints at increased abundance of Proteobacteria after travel, although this increase could not specifically be linked to resistance gene increases.
- Uses de novo metagenomic assembly to physically link resistance genes in the same sample, giving hints of co-resistance patterns in the gut microbiome.
The paper was a collaboration with Martin Angelin, Helena Palmgren and Anders Johansson at Umeå University, and was made possible by bioinformatics support from SciLifeLab in Stockholm. I highly recommend reading it as a complement to e.g. the Forslund et al. paper (12) describing country-specific antibiotic resistance patterns in the gut microbiota.
Taken together, this study offers a broadened perspective on how the antibiotic resistance potential of the human gut microbiome changes after travel, providing an independent complement to previous studies targeting a limited number of bacterial species or antibiotic resistance genes. Understanding how resistance genes travels the globe is hugely important, since resistance in principle only need to appear in a pathogen once; improper hygiene and travel may then spread novel resistance genes across continents rapidly (13,14).
- Bengtsson-Palme J, Angelin M, Huss M, Kjellqvist S, Kristiansson E, Palmgren H, Larsson DGJ, Johansson A: The human gut microbiome as a transporter of antibiotic resistance genes between continents. Antimicrob Agents Chemother Accepted manuscript posted online (2015). doi: 10.1128/AAC.00933-15 [Paper link]
- Gaarslev K, Stenderup J: Changes during travel in the composition and antibiotic resistance pattern of the intestinal Enterobacteriaceae flora: results from a study of mecillinam prophylaxis against travellers’ diarrhoea. Curr Med Res Opin 9:384–387 (1985).
- Paltansing S, Vlot JA, Kraakman MEM, Mesman R, Bruijning ML, Bernards AT, Visser LG, Veldkamp KE: Extended-spectrum β-lactamase-producing enterobacteriaceae among travelers from the Netherlands. Emerging Infect. Dis. 19:1206–1213 (2013).
- Ruppé E, Armand-Lefèvre L, Estellat C, El-Mniai A, Boussadia Y, Consigny PH, Girard PM, Vittecoq D, Bouchaud O, Pialoux G, Esposito-Farèse M, Coignard B, Lucet JC, Andremont A, Matheron S: Acquisition of carbapenemase-producing Enterobacteriaceae by healthy travellers to India, France, February 2012 to March 2013. Euro Surveill. 19 (2014).
- Kennedy K, Collignon P: Colonisation with Escherichia coli resistant to “critically important” antibiotics: a high risk for international travellers. Eur J Clin Microbiol Infect Dis 29:1501–1506 (2010).
- Tham J, Odenholt I, Walder M, Brolund A, Ahl J, Melander E: Extended-spectrum beta-lactamase-producing Escherichia coli in patients with travellers’ diarrhoea. Scand. J. Infect. Dis. 42:275–280 (2010).
- Östholm-Balkhed Å, Tärnberg M, Nilsson M, Nilsson LE, Hanberger H, Hällgren A, Travel Study Group of Southeast Sweden: Travel-associated faecal colonization with ESBL-producing Enterobacteriaceae: incidence and risk factors. J Antimicrob Chemother 68:2144–2153 (2013).
- Kantele A, Lääveri T, Mero S, Vilkman K, Pakkanen SH, Ollgren J, Antikainen J, Kirveskari J: Antimicrobials increase travelers’ risk of colonization by extended-spectrum betalactamase-producing enterobacteriaceae. Clin Infect Dis 60:837–846 (2015).
- von Wintersdorff CJH, Penders J, Stobberingh EE, Oude Lashof AML, Hoebe CJPA, Savelkoul PHM, Wolffs PFG: High rates of antimicrobial drug resistance gene acquisition after international travel, The Netherlands. Emerging Infect. Dis. 20:649–657 (2014).
- Tängdén T, Cars O, Melhus A, Löwdin E: Foreign travel is a major risk factor for colonization with Escherichia coli producing CTX-M-type extended-spectrum beta-lactamases: a prospective study with Swedish volunteers. Antimicrob Agents Chemother 54:3564–3568 (2010).
- Dhanji H, Patel R, Wall R, Doumith M, Patel B, Hope R, Livermore DM, Woodford N: Variation in the genetic environments of bla(CTX-M-15) in Escherichia coli from the faeces of travellers returning to the United Kingdom. J Antimicrob Chemother 66:1005–1012 (2011).
- Forslund K, Sunagawa S, Kultima JR, Mende DR, Arumugam M, Typas A, Bork P: Country-specific antibiotic use practices impact the human gut resistome. Genome Res 23:1163–1169 (2013).
- Bengtsson-Palme J, Larsson DGJ: Antibiotic resistance genes in the environment: prioritizing risks. Nat Rev Microbiol 13:396 (2015).
- Larsson DGJ: Antibiotics in the environment. Ups J Med Sci 119:108–112 (2014).
- Multiple input files can now be specified by adding several -1 and -2 options.
- TriMetAss now automatically stops if the candidate reads are the same for two iterations in a row.
- Support for recent versions of Trinity that no longer contain the Trinity.pl script.
- A minor bug causing TriMetAss to use more memory than necessary has been fixed.
- Fixed the
--stop_totaloption so that TriMetAss actually uses this option (rather than
- Allowed complicated paths to be supplied for the output directory.
I would like to thank users Rickard Hammarén, Dr. Tatsuya Unno, Dr. Gisle Vestergaard and Dr. Joseph Nesme for providing me with the underlying information to provide these fixes. Thanks a lot!
Late last year, an opinion paper by José Martínez, Teresa Coque and Fernando Baquero was published in Nature Reviews Microbiology (1). In this paper, the authors present a system – resistance readiness conditions (RESCon) – for ranking the risks associated with the detection of antibiotic resistance genes. They also outline the obstacles associated with determining risks presented by antibiotic resistance genes in environmental microbial communities in terms of their potential to transfer to human pathogens. Generally, I am very positive about this paper, which I think is a must-read for anyone who works with antibiotic resistance genes in metagenomes, regardless of it they stem from the human gut or the external environment.
There is, however, one very important aspect that struck me and many other members of our research group as curious: the proposed system assign antibiotic resistance genes already present on mobile genetic elements in human pathogens to the highest risk category (RESCon 1), while resistance genes encoding novel resistance mechanisms not yet been found on mobile elements in a pathogen are considered to be part of lower risk categories. We believe that this system will overestimate the risks associated with well-known resistance factors that are already circulating among human pathogens and under-appreciate the potentially disastrous consequences that the transfer of previously unknown resistance determinants from the environmental resistome could have (exemplified by the rapid clinical spread of the NDM-1 metallo-beta-lactamase gene (2,3)).
With this in mind me and Joakim Larsson wrote a response letter to Nature Reviews Microbiology that went online last monday (4), together with the authors’ reply to us (5). (I strongly suggest that you read the entire original paper (1) before you read the reply (5) to our response letter (4), since Martinez et al. changes the scope slightly from the original paper in their response letter, and these clarifications may (or may not) have been in response to our arguments.)
In our response, we also stress that the abundances of resistance genes, and not only their presence, should be accounted for when estimating risks (although that last point might have been slightly obscured due to the very low word limit). In other words, we think that identifying environmental hotspots for antibiotic resistance genes, where novel resistance genes could be selected for (6,7,8), is of great importance for mitigating public health risks related to environmental antibiotic resistance. Please read our full thoughts on the matter in Nature Reviews Microbiology.
- Martinez JL, Coque TM, Baquero F: What is a resistance gene? Ranking risk in resistomes. Nat Rev Microbiol 2015, 13:116–123.
- Kumarasamy KK, et al.: Emergence of a new antibiotic resistance mechanism in India, Pakistan, and the UK: a molecular, biological, and epidemiological study. Lancet Infect Dis 2010, 10:597–602.
- Walsh TR, Weeks J, Livermore DM, Toleman MA: Dissemination of NDM‐1 positive bacteria in the New Delhi environment and its implications for human health: an environmental point prevalence study. Lancet Infect Dis 2011, 11:355–362.
- Bengtsson-Palme J, Larsson DGJ: Antibiotic resistance genes in the environment: prioritizing risks. Nat Rev Microbiol 2015, Advance online publication. doi:10.1038/nrmicro3399‐c1
- Martinez JL, Coque TM, Baquero F: Prioritizing risks of antibiotic resistance genes in all metagenomes. Nat Rev Microbiol 2015, Advance online publication. doi:10.1038/nrmicro3399‐c2
- Kristiansson E, et al.: Pyrosequencing of antibiotic‐contaminated river sediments reveals high levels of resistance and gene transfer elements. PLoS ONE 2011, 6:e17038.
- Bengtsson‐Palme J, Boulund F, Fick J, Kristiansson E, Larsson DGJ: Shotgun metagenomics reveals a wide array of antibiotic resistance genes and mobile elements in a polluted lake in India. Front Microbiol 2014, 5:648.
- Marathe NP, et al.: A treatment plant receiving waste water from multiple bulk drug manufacturers is a reservoir for highly multi‐drug resistant integron‐bearing bacteria. PLoS ONE 2013, 8:e77310.
I will be giving a talk at the Third International symposium on the environmental dimension of antibiotic resistance (EDAR2015) next month (five weeks from now. The talk is entitled “Turn up the signal – wipe out the noise: Gaining insights into antibiotic resistance of bacterial communities using metagenomic data“, and will deal with handling of metagenomic data in antibiotic resistance gene research. The talk will highlight the some particular pitfalls related to interpretation of data, and exemplify how flawed analysis practices can result in misleading conclusions regarding antibiotic resistance risks. I will particularly address how taxonomic composition influences the frequencies of resistance genes, the importance of knowledge of the functions of the genes in the databases used, and how normalization strategies influence the results. Furthermore, we will show how the context of resistance genes can allow inference of their potential to spread to human pathogens from environmental or commensal bacteria. All these aspects will be exemplified by data from our studies of environments subjected to pharmaceutical pollution in India, the effect of travel on the human resistome, and modern municipal wastewater treatment processes.
The talk will take place on Monday, May 18, 2015 at 13:20. The full scientific program for the conference can be found here. Registration for the conference is still possible, although not for the early-bird price. I look forward to see a lot of the people who will attend the conference, and hopefully also you!
Metaxa2 has been updated to version 2.0.2 and can be downloaded from the Metaxa2 web site. The 2.0.2 update fixes two minor bugs; one causing the “.graph” file to display incorrect or no names for the regions of the LSU regions, and one causing misreporting of the number of sequences in single-end FASTQ files (paired-end files were reported correctly). The update also brings a slightly improved classifier. Thanks to Marco Severgnini for reporting the FASTQ file issue! The update is available here.
A couple of days ago a paper was published in Environmental Sciences Europe summarizing the EU report on effect-based tools for use in toxicology in the aquatic environment I have been involved in (1). This report was officially published last spring (2), and can be found here, with the annex available on the European Commission document website. My contribution to the paper was, as with the report, in the genomics and metagenomics section. The paper briefly presents modern bioassays, biomarkers and ecological methods that can be used for aquatic monitoring of the environment.
- Wernersson A-S, Carere M, Maggi C, Tusil P, Soldan P, James A, Sanchez W, Dulio V, Broeg K, Reifferscheid G, Buchinger S, Maas H, Van Der Grinten E, O’Toole S, Ausili A, Manfra L, Marziali L, Polesello S, Lacchetti I, Mancini L, Lilja K, Linderoth M, Lundeberg T, Fjällborg B, Porsbring T, Larsson DGJ, Bengtsson-Palme J, Förlin L, Kienle C, Kunz P, Vermeirssen E, Werner I, Robinson CD, Lyons B, Katsiadaki I, Whalley C, den Haan K, Messiaen M, Clayton H, Lettieri T, Negrão Carvalho R, Gawlik BM, Hollert H, Di Paolo C, Brack W. Kammann U, Kase R: The European technical report on aquatic effect-based monitoring tools under the water framework directive. Environmental Sciences Europe, 27, 7 (2015). doi: 10.1186/s12302-015-0039-4 [Paper link]
- Wernersson A-S, Carere M, Maggi C, Tusil P, Soldan P, James A, Sanchez W, Broeg K, Kammann U, Reifferscheid G, Buchinger S, Maas H, Van Der Grinten E, Ausili A, Manfra L, Marziali L, Polesello S, Lacchetti I, Mancini L, Lilja K, Linderoth M, Lundeberg T, Fjällborg B, Porsbring T, Larsson DGJ, Bengtsson-Palme J, Förlin L, Kase R, Kienle C, Kunz P, Vermeirssen E, Werner I, Robinson CD, Lyons B, Katsiadaki I, Whalley C, den Haan K, Messiaen M, Clayton H, Lettieri T, Negrão Carvalho R, Gawlik BM, Dulio V, Hollert H, Di Paolo C, Brack W (2014). Technical Report on Aquatic Effect-Based Monitoring Tools. European Commission. Technical Report 2014-077, Office for Official Publications of European Communities, ISBN: 978-92-79-35787-9. doi:10.2779/7260
After almost a year in different stages of review and revision, in which the paper (but not the software) saw a total transformation, I am happy to announce that the paper describing Metaxa2 has been accepted in Molecular Ecology Resources and is available in a rudimentary online early form. The figures in this version are not that pretty, but those who wants to read the paper asap, you have the possibility to do so.
This means that if you have been using Metaxa2 for a publication, there is now a new preferred way of citing this, namely:
Bengtsson-Palme J, Hartmann M, Eriksson KM, Pal C, Thorell K, Larsson DGJ, Nilsson RH: Metaxa2: Improved Identification and Taxonomic Classification of Small and Large Subunit rRNA in Metagenomic Data. Molecular Ecology Resources (2015). doi: 10.1111/1755-0998.12399
The paper (1), apart from describing the new Metaxa version, also brings a very thorough evaluation of the software, compared to other tools for taxonomic classification implemented in QIIME (2). In short, we show that:
- Metaxa2 can make trustworthy taxonomic classifications even with reads as short as 100 bp
- Generally, the performance is reliable across the entire SSU rRNA gene, regardless of which V-region a read is derived from
- Metaxa2 can reliably recapture species composition from short-read metagenomic data, comparable with results of amplicon sequencing
- Metaxa2 outperforms other popular tools such as Mothur (3), the RDP Classifier (4), Rtax (5) and the QIIME implementation of Uclust (6) in terms of proportion of correctly classified reads from metagenomic data
- The false positive rate of Metaxa2 is very close to zero; far superior to many of the above mentioned tools, many of which assume that reads must derive from the rRNA gene
Metaxa2 can be downloaded here. We have already used it for around two years internally, and it forms the base of the taxonomic classifications in e.g. our recently published paper on antibiotic resistance in a polluted Indian lake (7).
- Bengtsson-Palme J, Hartmann M, Eriksson KM, Pal C, Thorell K, Larsson DGJ, Nilsson RH: Metaxa2: Improved Identification and Taxonomic Classification of Small and Large Subunit rRNA in Metagenomic Data. Molecular Ecology Resources (2015). doi: 10.1111/1755-0998.12399 [Paper link]
- Caporaso JG, Kuczynski J, Stombaugh J et al.: QIIME allows analysis of high-throughput community sequencing data. Nature Methods, 7, 335–336 (2010).
- Schloss PD, Westcott SL, Ryabin T et al.: Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Applied and Environmental Microbiology, 75, 7537–7541 (2009).
- Wang Q, Garrity GM, Tiedje JM, Cole JR: Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Applied and Environmental Microbiology, 73, 5261–5267 (2007).
- Soergel DAW, Dey N, Knight R, Brenner SE: Selection of primers for optimal taxonomic classification of environmental 16S rRNA gene sequences. The ISME Journal, 6, 1440–1444 (2012).
- Edgar RC: Search and clustering orders of magnitude faster than BLAST. Bioinformatics, 26, 2460–2461 (2010).
- Bengtsson-Palme J, Boulund F, Fick J, Kristiansson E, Larsson DGJ: Shotgun metagenomics reveals a wide array of antibiotic resistance genes and mobile elements in a polluted lake in India. Frontiers in Microbiology, 5, 648 (2014).
Our paper describing the bacterial community of a polluted lake in India has now been typeset and appears in its final form in Frontiers in Microbiology. If I may say so, I think that the paper turned out to be very goodlooking and it is indeed nice to finally see it in print. The paper describes an unprecedented diversity and abundance of antibiotic resistance genes and genes enabling transfer of DNA between bacteria. We also describe a range of potential novel plasmids from the lake. Finally, the paper briefly describes a new approach to targeted assembly of metagenomic data — TriMetAss — which can be downloaded here.
Bengtsson-Palme J, Boulund F, Fick J, Kristiansson E, Larsson DGJ: Shotgun metagenomics reveals a wide array of antibiotic resistance genes and mobile elements in a polluted lake in India. Frontiers in Microbiology, 5, 648 (2014). doi: 10.3389/fmicb.2014.00648
With the publication of my latest paper last week (1), I also would like to highlight some of the software underpinning the findings a bit. To get around the problem that extremely common resistance genes could be present in multiple contexts and variants, causing assembler such as Velvet (2) to perform sub-optimally, we have written a software tool that utilizes Vmatch (3) and Trinity (4) to iteratively construct contigs from reads associated with resistance genes. This could of course be used in many other situations as well, when you want to specifically assemble a certain portion of a metagenome, but suspect that that portion might be found in multiple contexts.
TriMetAss is a Perl program, employing Vmatch and Trinity to construct multi-context contigs. TriMetAss uses extracted reads associated with, e.g., resistance genes as seeds for a Vmatch search against the complete set of read pairs, extracting reads matching with at least 49 bp (by default) to any of the seed reads. These reads are then assembled using Trinity. The resulting contigs are then used as seeds for another search using Vmatch to the complete set of reads, as above. All matches (including the previously matching read pairs) are again then used for a Trinity assembly. This iterative process is repeated until a stop criteria is met, e.g. when the total number of assembled nucleotides starts to drop rather than increase. The software can be downloaded here.
- Bengtsson-Palme J, Boulund F, Fick J, Kristiansson E, Larsson DGJ: Shotgun metagenomics reveals a wide array of antibiotic resistance genes and mobile elements in a polluted lake in India. Frontiers in Microbiology, 5, 648 (2014). doi: 10.3389/fmicb.2014.00648
- Zerbino DR, Birney E: Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res 18, 821–829 (2008). doi:10.1101/gr.074492.107
- Kurtz S: The Vmatch large scale sequence analysis software (2010). http://vmatch.de/
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, et al.: Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol 29, 644–652 (2011). doi:10.1038/nbt.1883