I am happy to announce that the opinion/review piece I wrote for Current Opinion in Food Science has been published. The paper (1) extends on some of my thoughts on how high-throughput sequencing and metagenomics can aid in risk assessment of antibiotic resistant bacteria that I outlined in my PhD thesis (2), but specifically focuses on the food supply chain and its role in resistance dissemination and selection.
In the paper, I argue for that the food supply chain is a special type of setting in the resistance puzzle, as it not only serves as a connection between environmental habitats for bacteria and humans, but also sometimes presents a substantial selection for resistance, due to use of antibiotics in agri- and aquaculture. International food standards are clear that both selection and dissemination of foodborne resistance should be considered in the risk analysis of food production (3). However, the current main use of DNA sequencing in food safety is whole genome sequencing to delineate which specific strains that are involved in foodborne disease outbreaks, including the resistance factors they may carry (4,5). Further, I argue that while shotgun metagenomics could be used to screen samples for a large number of genes involved in resistance and virulence in the food supply chain, it would at present be very costly and therefore of doubtful benefit to employ in routine screening programs. Still, metagenomics can contribute knowledge that can be used in quantitative risk assessment of antibiotic resistance in the food supply chain.
The entire paper can be read here.
- Bengtsson-Palme J: Antibiotic resistance in the food supply chain: Where can sequencing and metagenomics aid risk assessment? Current Opinion in Food Science, in press (2017). doi: 10.1016/j.cofs.2017.01.010 [Paper link]
- Bengtsson-Palme J: Antibiotic resistance in the environment: a contribution from metagenomic studies. Doctoral thesis (medicine), Department of Infectious Diseases, Institute of Biomedicine, Sahlgrenska Academy, University of Gothenburg, 2016. [Link]
- Codex Alimentarius Commission: Guidelines for risk analysis of foodborne antimicrobial resistance. Food and Agriculture Organization of the United Nations & World Health Organization2011. [Link]
- Franz E, Gras LM, Dallman T: Significance of whole genome sequencing for surveillance, source attribution and microbial risk assessment of foodborne pathogens. Current Opinion in Food Science, 8, 74-79 (2016). doi: 10.1016/j.cofs.2016.04.004
- Stasiewicz MJ, Bakker den HC, Wiedmann M: Genomics tools in microbial food safety. Current Opinion in Food Science, 4, 105-110 (2015). doi: 10.1016/j.cofs.2015.06.002
Late yesterday, Microbiome put online our most recent work, covering the resistomes to antibiotics, biocides and metals across a vast range of environments. In the paper (1), we perform the largest characterization of resistance genes, mobile genetic elements (MGEs) and bacterial taxonomic compositions to date, covering 864 different metagenomes from humans (350), animals (145) and external environments such as soil, water, sewage, and air (369 in total). All the investigated metagenomes were sequenced to at least 10 million reads each, using Illumina technology, making the results more comparable across environments than in previous studies (2-4).
We found that the environment types had clear differences both in terms of resistance profiles and bacterial community composition. Humans and animals hosted microbial communities with limited taxonomic diversity as well as low abundance and diversity of biocide/metal resistance genes and MGEs. On the contrary, the abundance of ARGs was relatively high in humans and animals. External environments, on the other hand, showed high taxonomic diversity and high diversity of biocide/metal resistance genes and MGEs. Water, sediment and soil generally carried low relative abundance and few varieties of known ARGs, whereas wastewater and sludge were on par with the human gut. The environments with the largest relative abundance and diversity of ARGs, including genes encoding resistance to last resort antibiotics, were those subjected to industrial antibiotic pollution and air samples from a Beijing smog event.
A paper investigating this vast amount of data is of course hard to describe in a blog post, so I strongly suggest the interested reader to head over to Microbiome’s page and read the full paper (1). However, here’s a ver short summary of the findings:
- The median relative abundance of ARGs across all environments was 0.035 copies per bacterial 16S rRNA
- Antibiotic-polluted environments have (by far) the highest abundances of ARGs
- Urban air samples carried high abundance and diversity of ARGs
- Human microbiota has high abundance and diversity of known ARGs, but low taxonomic diversity compared to the external environment
- The human and animal resistomes are dominated by tetracycline resistance genes
- Over half of the ARGs were only detected in external environments, while 20.5 % were found in human, animal and at least one of the external environments
- Biocide and metal resistance genes are more common in external environments than in the human microbiota
- Human microbiota carries low abundance and richness of MGEs compared to most external environments
Importantly, less than 1.5 % of all detected ARGs were found exclusively in the human microbiome. At the same time, 57.5 % of the known ARGs were only detected in metagenomes from environmental samples, despite that the majority of the investigated ARGs were initially encountered in pathogens. Still, our analysis suggests that most of these genes are relatively rare in the human microbiota. Environmental samples generally contained a wider distribution of resistance genes to a more diverse set of antibiotics classes. For example, the relative abundance of beta-lactam resistance genes was much larger in external environments than in human and animal microbiomes. This suggests that the external environment harbours many more varieties of resistance genes than the ones currently known from the clinic. Indeed, functional metagenomics has resulted in the discovery of many novel ARGs in external environments (e.g. 5). This all fits well with an overall much higher taxonomic diversity of environmental microbial communities. In terms of consequences associated with the potential transfer of ARGs to human pathogens, we argue that unknown resistance genes are of greater concern than those already known to circulate among human-associated bacteria (6).
This study describes the potential for many external environments, including those subjected to pharmaceutical pollution, air and wastewater/sludge, to serve as hotspots for resistance development and/or transmission of ARGs. In addition, our results indicate that these environments may play important roles in the mobilization of yet unknown ARGs and their further transmission to human pathogens. To provide guidance for risk-reducing actions we – based on this study – suggest strict regulatory measures of waste discharges from pharmaceutical industries and encourage more attention to air in the transmission of antibiotic resistance (1).
- Pal C, Bengtsson-Palme J, Kristiansson E, Larsson DGJ: The structure and diversity of human, animal and environmental resistomes. Microbiome, 4, 54 (2016). doi: 10.1186/s40168-016-0199-5
- Durso LM, Miller DN, Wienhold BJ. Distribution and quantification of antibiotic resistant genes and bacteria across agricultural and non-agricultural metagenomes. PLoS One. 2012;7:e48325.
- Nesme J, Delmont TO, Monier J, Vogel TM. Large-scale metagenomic-based study of antibiotic resistance in the environment. Curr Biol. 2014;24:1096–100.
- Fitzpatrick D, Walsh F. Antibiotic resistance genes across a wide variety of metagenomes. FEMS Microbiol Ecol. 2016. doi:10.1093/femsec/fiv168.
- Allen HK, Moe LA, Rodbumrer J, Gaarder A, Handelsman J. Functional metagenomics reveals diverse β-lactamases in a remote Alaskan soil. ISME J. 2009;3:243–51.
- Bengtsson-Palme J, Larsson DGJ: Antibiotic resistance genes in the environment: prioritizing risks. Nature Reviews Microbiology, 13, 369 (2015). doi: 10.1038/nrmicro3399-c1
After a long wait (1), Science of the Total Environment has finally decided to make our paper on selection of antibiotic resistance genes in sewage treatment plants (STPs) available (2). STPs are often suggested to be “hotspots” for emergence and dissemination of antibiotic-resistant bacteria (3-6). However, we actually do not know if the selection pressures within STPs, that can be caused either by residual antibiotics or other co-selective agents, are sufficiently large to specifically promote resistance. To better understand this, we used shotgun metagenomic sequencing of samples from different steps of the treatment process (incoming water, treated water, primary sludge, recirculated sludge and digested sludge) in three Swedish STPs in the Stockholm area to characterize the frequencies of resistance genes to antibiotics, biocides and metal, as well as mobile genetic elements and taxonomic composition. In parallel, we also measured concentrations of antibiotics, biocides and metals.
We found that only the concentrations of tetracycline and ciprofloxacin in the influent water were above those that we predict to cause resistance selection (7). However, there was no consistent enrichment of resistance genes to any particular class of antibiotics in the STPs, neither for biocide and metal resistance genes. Instead, the most substantial change of the bacterial communities compared to human feces (sampled from Swedes in another study of ours (8)) occurred already in the sewage pipes, and was manifested by a strong shift from obligate to facultative anaerobes. Through the treatment process, resistance genes against antibiotics, biocides and metals were not reduced to the same extent as fecal bacteria were.
Worryingly, the OXA-48 beta-lactamase gene was consistently enriched in surplus and digested sludge. OXA-48 is still rare in Swedish clinical isolates (9), but provides resistance to carbapenems, one of our most critically important classes of antibiotics. However, taken together metagenomic sequencing did not provide clear support for any specific selection of antibiotic resistance. Rather, since stronger selective forces affect gross taxonomic composition, and thereby also resistance gene abundances, it is very hard to interpret the metagenomic data from a risk-for-selection perspective. We therefore think that comprehensive analyses of resistant vs. non-resistant strains within relevant species are warranted.
Taken together, the main take-home messages of the paper (2) are:
- There were no apparent evidence for direct selection of resistance genes by antibiotics or co-selection by biocides or metals
- Abiotic factors (mostly oxygen availability) strongly shape taxonomy and seems to be driving changes of resistance genes
- Metagenomic and/or PCR-based community studies may not be sufficiently sensitive to detect selection effects, as important shifts towards resistant may occur within species and not on the community level
- The concentrations of antibiotics, biocides and metals were overall reduced, but not removed in STPs. Incoming concentrations of antibiotics in Swedish STPs are generally low
- Resistance genes are overall reduced through the treatment process, but far from eliminated
References and notes
- Okay, those who takes notes know that I have already complained once before on Science of the Total Environment’s ridiculously long production handling times. But, seriously, how can a journal’s production team return the proofs for after three days of acceptance, and then wait seven weeks before putting the final proofs online? I still wonder what is going on beyond the scenes, which is totally obscure because the production office also refuses to respond to e-mails. Not a nice publishing experience this time either.
- Bengtsson-Palme J, Hammarén R, Pal C, Östman M, Björlenius B, Flach C-F, Kristiansson E, Fick J, Tysklind M, Larsson DGJ: Elucidating selection processes for antibiotic resistance in sewage treatment plants using metagenomics. Science of the Total Environment, in press (2016). doi: 10.1016/j.scitotenv.2016.06.228 [Paper link]
- Rizzo L, Manaia C, Merlin C, Schwartz T, Dagot C, Ploy MC, Michael I, Fatta-Kassinos D: Urban wastewater treatment plants as hotspots for antibiotic resistant bacteria and genes spread into the environment: a review. Science of the Total Environment, 447, 345–360 (2013). doi: 10.1016/j.scitotenv.2013.01.032
- Laht M, Karkman A, Voolaid V, Ritz C, Tenson T, Virta M, Kisand V: Abundances of Tetracycline, Sulphonamide and Beta-Lactam Antibiotic Resistance Genes in Conventional Wastewater Treatment Plants (WWTPs) with Different Waste Load. PLoS ONE, 9, e103705 (2014). doi: 10.1371/journal.pone.0103705
- Yang Y, Li B, Zou S, Fang HHP, Zhang T: Fate of antibiotic resistance genes in sewage treatment plant revealed by metagenomic approach. Water Research, 62, 97–106 (2014). doi: 10.1016/j.watres.2014.05.019
- Berendonk TU, Manaia CM, Merlin C, Fatta-Kassinos D, Cytryn E, Walsh F, et al.: Tackling antibiotic resistance: the environmental framework. Nature Reviews Microbiology, 13, 310–317 (2015). doi: 10.1038/nrmicro3439
- Bengtsson-Palme J, Larsson DGJ: Concentrations of antibiotics predicted to select for resistant bacteria: Proposed limits for environmental regulation. Environment International, 86, 140–149 (2016). doi: 10.1016/j.envint.2015.10.015
- 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. Antimicrobial Agents and Chemotherapy, 59, 10, 6551–6560 (2015). doi: 10.1128/AAC.00933-15
- Hellman J, Aspevall O, Bengtsson B, Pringle M: SWEDRES-SVARM 2014. Consumption of antimicrobials and occurrence of antimicrobial resistance in Sweden. Public Health Agency of Sweden and National Veterinary Institute, Solna/Uppsala, Sweden. Report No.: 14027. Available from: http://www.folkhalsomyndigheten.se/publicerat-material/ (2014)
I am happy to announce that our Viewpoint article on strategies for improving sequence databases has now been published in the journal Proteomics. The paper (1) defines some central problems hampering genomic, proteomic and metagenomic analyses and suggests five strategies to improve the situation:
- Clearly separate experimentally verified and unverified sequence entries
- Enable a system for tracing the origins of annotations
- Separate entries with high-quality, informative annotation from less useful ones
- Integrate automated quality-control software whenever such tools exist
- Facilitate post-submission editing of annotations and metadata associated with sequences
The paper is not long, so I encourage you to read it in its entirety. We believe that spreading this knowledge and pushing solutions to problems related to poor annotation metadata is vastly important in this era of big data. Although we specifically address protein-coding genes in this paper, the same logic also applies to other types of biological sequences. In this way the paper is related to my previous work with Henrik Nilsson on improving annotation data for taxonomic barcoding genes (2-4). This paper was one of the main end-results of the GoBiG network, and the backstory on the paper follows below the references…
- Bengtsson-Palme J, Boulund F, Edström R, Feizi A, Johnning A, Jonsson VA, Karlsson FH, Pal C, Pereira MB, Rehammar A, Sánchez J, Sanli K, Thorell K: Strategies to improve usability and preserve accuracy in biological sequence databases. Proteomics, Early view (2016). doi: 10.1002/pmic.201600034
- Kõljalg U, Nilsson RH, Abarenkov K, Tedersoo L, Taylor AFS, Bahram M, Bates ST, Bruns TT, Bengtsson-Palme J, Callaghan TM, Douglas B, Drenkhan T, Eberhardt U, Dueñas M, Grebenc T, Griffith GW, Hartmann M, Kirk PM, Kohout P, Larsson E, Lindahl BD, Lücking R, Martín MP, Matheny PB, Nguyen NH, Niskanen T, Oja J, Peay KG, Peintner U, Peterson M, Põldmaa K, Saag L, Saar I, Schüßler A, Senés C, Smith ME, Suija A, Taylor DE, Telleria MT, Weiß M, Larsson KH: Towards a unified paradigm for sequence-based identification of Fungi. Molecular Ecology, 22, 21, 5271–5277 (2013). doi: 10.1111/mec.12481
- Nilsson RH, Hyde KD, Pawlowska J, Ryberg M, Tedersoo L, Aas AB, Alias SA, Alves A, Anderson CL, Antonelli A, Arnold AE, Bahnmann B, Bahram M, Bengtsson-Palme J, Berlin A, Branco S, Chomnunti P, Dissanayake A, Drenkhan R, Friberg H, Frøslev TG, Halwachs B, Hartmann M, Henricot B, Jayawardena R, Jumpponen A, Kauserud H, Koskela S, Kulik T, Liimatainen K, Lindahl B, Lindner D, Liu J-K, Maharachchikumbura S, Manamgoda D, Martinsson S, Neves MA, Niskanen T, Nylinder S, Pereira OL, Pinho DB, Porter TM, Queloz V, Riit T, Sanchez-García M, de Sousa F, Stefaczyk E, Tadych M, Takamatsu S, Tian Q, Udayanga D, Unterseher M, Wang Z, Wikee S, Yan J, Larsson E, Larsson K-H, Kõljalg U, Abarenkov K: Improving ITS sequence data for identification of plant pathogenic fungi. Fungal Diversity, 67, 1, 11–19 (2014). doi: 10.1007/s13225-014-0291-8
- Nilsson RH, Tedersoo L, Ryberg M, Kristiansson E, Hartmann M, Unterseher M, Porter TM, Bengtsson-Palme J, Walker D, de Sousa F, Gamper HA, Larsson E, Larsson K-H, Kõljalg U, Edgar R, Abarenkov K: A comprehensive, automatically updated fungal ITS sequence dataset for reference-based chimera control in environmental sequencing efforts. Microbes and Environments, 30, 2, 145–150 (2015). doi: 10.1264/jsme2.ME14121
In June 2013, the Gothenburg Bioinformatics Group for junior scientists (GoBiG) arranged a workshop with two themes: “Parallelized quantification of genes in large metagenomic datasets” and “Assigning functional predictions to NGS data”. The following discussion on how to database quality influenced results and what could be done to improve the situation was rather intense, and several good ideas were thrown around. I took notes from the meeting, and in the evening I put them down during a warm summer night at the balcony. In fact, the notes were good enough to be an early embryo for a manuscript. So I sent it to some of the most active GoBiG members (Kaisa Thorell and Fredrik Boulund), who were positive regarding the idea to turn it into a manuscript. I wrote it together more properly and we decided that everyone who contributed with ideas at the meeting would be invited to become co-authors. We submitted the manuscript in early 2014, only to see it (rather brutally) rejected. At that point most of us were sucked up in their own projects, so nothing happened to this manuscript for over a year. Then we decided to give it another go, updated the manuscript heavily and changed a few parts to better reflect the current database situation (at this point, e.g., UniProt had already started implementing some of our suggested ideas). Still, some of the proposed strategies were more radical in 2013 than they would be now, more than three years later. We asked the Proteomics editors if they would be interested in the manuscript, and they turned out to be very positive. Indeed, the entire experience with the editors at Proteomics has been very pleasant. I am very thankful to the GoBiG team for this time, and to the editors at Proteomics who saw the value of this manuscript.
Late last year, we introduced FARAO – the Flexible All-Round Annotation Organizer – a software tool that allows visualization of annotated features on contigs. Today, the Applications Note describing the software was published as an advance access paper in Bioinformatics (1). As I have described before, storing and visualizing annotation and coverage information in FARAO has a number of advantages. FARAO is able to:
- Integrate annotation and coverage information for the same sequence set, enabling coverage estimates of annotated features
- Scale across millions of sequences and annotated features
- Filter sequences, such that only entries with annotations satisfying certain given criteria will be outputted
- Handle annotation and coverage data produced by a range of different bioinformatics tools
- Handle custom parsers through a flexible interface, allowing for adaption of the software to virtually any bioinformatic tool not supported out of the box
- Produce high-quality EPS output
- Integrate with MySQL databases
I have previously used FARAO to produce annotation figures in our paper on a polluted Indian lake (2), as well as in a paper on sewage treatment plants (which is in press and should be coming out any day now). We hope that the tool will find many more uses in other projects in the future!
- Hammarén R, Pal C, Bengtsson-Palme J: FARAO: The Flexible All-Round Annotation Organizer. Bioinformatics, advance access (2016). doi: 10.1093/bioinformatics/btw499 [Paper link]
- 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 [Paper link]
So, on Thursday (May 26th) I will defend my thesis, titled “Antibiotic resistance in the environment: a contribution from metagenomic studies”. I will not dwell into this by writing a novel text, but will instead shamelessly reproduce the press release, which should give a reasonable overview of what I have been doing:
More and more people are infected with antibiotic resistant bacteria. But how do bacteria become resistant? A doctoral thesis from the Centre for Antibiotic Resistance Research at University of Gothenburg has investigated the role of the environment in the development of antibiotic resistance.
“An important question we asked was how low concentrations of antibiotics that can favour the growth of resistant bacteria in the environment”, says Johan Bengtsson-Palme, author of the thesis.
“Based on our analyses, we propose emission limits for 111 antibiotics that should not be exceeded in order to avoid that environmental bacteria become more resistant.”
A starting point to regulate antibiotic pollution
A recent report, commissioned by the British Prime Minister David Cameron, proposes that the emission limits suggested in Johan’s thesis should be used as a starting point to regulate antibiotic pollution from, for example, pharmaceutical production – globally.
“Many people are surprised that such regulations are not already in place, but today it is actually not a crime to discharge wastewater contaminated with large amounts of antibiotics, not even in Europe”, says Johan Bengtsson-Palme.
In one of the studies in the thesis, the researchers show that resistance genes against a vast range of antibiotics are enriched in an Indian lake polluted by dumping of wastewater from pharmaceutical production.
“It’s scary. Not only do the bacteria carry a multitude of resistance genes. They are also unusually well adapted to share those genes with other bacteria. If a disease-causing bacterium ends up in the lake, it may quickly pick up the genes it needs to become resistant. Since the lake is located close to residential areas, such spread of resistant bacteria to humans is not hard to imagine”, says Johan Bengtsson-Palme.
Spreading by travelers
The thesis also shows that resistant bacteria spread in the intestines of travelers who have visited India or Central Africa, even if the travelers themselves have not become sick.
“That resistant bacteria spread so quickly across the planet highlights that we must adopt a global perspective on the resistance problem”, says Johan Bengtsson-Palme. “Furthermore, it is not enough to reduce the use of antibiotics in healthcare. We must also reduce the use of antibiotics for animals, and try to limit the releases of antibiotics into the environment to try to get control over the growing antibiotic resistance problem before it is too late”.
The thesis Antibiotic resistance in the environment: a contribution from metagenomic studies will be defended on a dissertation on May 26th.
Today marks the five year anniversary for the Metaxa software’s initial release. Much has happened to the software since; Metaxa started off as an rRNA extraction utility for metagenomic data (1), including coarse classification to organism/organelle type. Since it has gained full-scale taxonomic classification ability better or on par with other software packages (2), much greater speed, support for the LSU gene, gained a range of related software tools (3), and spurred development of other tools such as ITSx (4). I have also been involved in no less than four peer-reviewed publications directly related to the software (1-3,5).
But it does not end here; these five years were just the beginning. We are – in different constellations – working on further enhancements to Metaxa2, including support for more genes, an updated classification database, and better customization options. I am very much still devoted to keep Metaxa2 alive and relevant as a tool for taxonomic analysis of metagenomes, applicable whenever accuracy is a key parameter. Thanks for being part of the community for these five years!
- Bengtsson J, Eriksson KM, Hartmann M, Wang Z, Shenoy BD, Grelet G, Abarenkov K, Petri A, Alm Rosenblad M, Nilsson RH: Metaxa: A software tool for automated detection and discrimination among ribosomal small subunit (12S/16S/18S) sequences of archaea, bacteria, eukaryotes, mitochondria, and chloroplasts in metagenomes and environmental sequencing datasets. Antonie van Leeuwenhoek, 100, 3, 471–475 (2011). doi:10.1007/s10482-011-9598-6. [Paper link]
- 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, 15, 6, 1403–1414 (2015). doi: 10.1111/1755-0998.12399 [Paper link]
- Bengtsson-Palme J, Thorell K, Wurzbacher C, Sjöling Å, Nilsson RH: Metaxa2 Diversity Tools: Easing microbial community analysis with Metaxa2. Ecological Informatics, 33, 45–50 (2016). doi: 10.1016/j.ecoinf.2016.04.004 [Paper link]
- Bengtsson-Palme J, Ryberg M, Hartmann M, Branco S, Wang Z, Godhe A, De Wit P, Sánchez-García M, Ebersberger I, de Souza F, Amend AS, Jumpponen A, Unterseher M, Kristiansson E, Abarenkov K, Bertrand YJK, Sanli K, Eriksson KM, Vik U, Veldre V, Nilsson RH: Improved software detection and extraction of ITS1 and ITS2 from ribosomal ITS sequences of fungi and other eukaryotes for use in environmental sequencing. Methods in Ecology and Evolution, 4, 10, 914–919 (2013). doi: 10.1111/2041-210X.12073 [Paper link]
- Bengtsson-Palme J, Hartmann M, Eriksson KM, Nilsson RH: Metaxa, overview. In:Nelson K. (Ed.) Encyclopedia of Metagenomics: SpringerReference (www.springerreference.com). Springer-Verlag Berlin Heidelberg (2013). doi: 10.1007/978-1-4614-6418-1_239-6 [Link]
Yesterday, Ecological Informatics put our paper describing Metaxa2 Diversity Tools online (1). Metaxa2 Diversity Tools was introduced with Metaxa2 version 2.1 and consists of
- metaxa2_dc – a tool for collecting several .taxonomy.txt output files into one large abundance matrix, suitable for analysis in, e.g., R
- metaxa2_rf – generates resampling rarefaction curves (2) based on the .taxonomy.txt output
- metaxa2_si – species inference based on guessing species data from the other species present in the .taxonomy.txt output file
- metaxa2_uc – a tool for determining if the community composition of a sample is significantly different from others through resampling analysis
At the same time as I did this update to the web site, I also took the opportunity to update the Metaxa2 FAQ to better reflect recent updates to the Metaxa2 software.
Metaxa2 Diversity Tools
One often requested feature of Metaxa2 (3) has been the ability to make simple analyses from the data after classification. The Metaxa2 Diversity Tools included in Metaxa2 2.1 is a seed for such an effort (although not close to a full-fledged community analysis package comparable to QIIME (4) or Mothur (5)). It currently consist of four tools.
The Metaxa2 Data Collector (metaxa2_dc) is the simplest of them (but probably the most requested), designed to merge the output of several *.level_X.txt files from the Metaxa2 Taxonomic Traversal Tool into one large abundance matrix, suitable for further analysis in, for example, R. The Metaxa2 Species Inference tool (metaxa2_si) can be used to further infer taxon information on, for example, the species level at a lower reliability than what would be permitted by the Metaxa2 classifier, using a complementary algorithm. The idea is that is if only a single species is present in, e.g., a family and a read is assigned to this family, but not classified to the species level, that sequence will be inferred to the same species as the other reads, given that it has more than 97% sequence identity to its best reference match. This can be useful if the user really needs species or genus classifications but many organisms in the studied species group have similar rRNA sequences, making it hard for the Metaxa2 classifier to classify sequences to the species level.
The Metaxa2 Rarefaction analysis tool (metaxa2_rf) performs a resampling rarefaction analysis (2) based on the output from the Metaxa2 classifier, taking into account also the unclassified portion of rRNAs. The Metaxa2 Uniqueness of Community analyzer (metaxa2_uc), finally, allows analysis of whether the community composition of two or more samples or groups is significantly different. Using resampling of the community data, the null hypothesis that the taxonomic content of two communities is drawn from the same set of taxa (given certain abundances) is tested. All these tools are further described in the manual and the recent paper (1).
The latest version of Metaxa2, including the Metaxa2 Diversity Tools, can be downloaded here.
- Bengtsson-Palme J, Thorell K, Wurzbacher C, Sjöling Å, Nilsson RH: Metaxa2 Diversity Tools: Easing microbial community analysis with Metaxa2. Ecological Informatics, 33, 45–50 (2016). doi: 10.1016/j.ecoinf.2016.04.004 [Paper link]
- Gotelli NJ, Colwell RK: Quantifying biodiversity: procedures and pitfalls in the measurement and comparison of species richness. Ecology Letters, 4, 379–391 (2000). doi:10.1046/j.1461-0248.2001.00230.x
- 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).
After a long wait (1) Sara Lundström’s paper establishing minimal selective concentrations (MSCs) for the antibiotic tetracycline in complex microbial communities (2), of which I am a co-author, has gone online. Personally, I think this paper is among the finest work I have been involved in; a lot of good science have gone into this publication. Risk assessment and management of antibiotics pollution is in great need of scientific data to underpin mitigation efforts (3). This paper describes a method to determine the minimal selective concentrations of antibiotics, and investigates different endpoints for measuring those MSCs. The method involves a testing system highly relevant for aquatic communities, in which bacteria are allowed to form biofilms in aquaria under controlled antibiotic exposure. Using the system, we find that 1 μg/L tetracycline selects for the resistance genes tetA and tetG, while 10 μg/L tetracycline is required to detect changes of phenotypic resistance. In short, the different endpoints studied (and their corresponding MSCs) were:
- CFU counts on R2A plates with 20 μg/mL tetracycline – MSC = 10 μg/L
- MIC range – MSC ~ 10-100 μg/L
- PICT, leucine uptake after short-term TC challenge – MSC ~ 100 μg/L
- Increased resistance gene abundances, metagenomics – MSC range: 0.1-10 μg/L
- Increased resistance gene abundances, qPCR (tetA and tetG) – MSC ≤ 1 μg/L
- Changes to taxonomic diversity – no significant changes detected
- Changes to taxonomic community composition – MSC ~ 1-10 μg/L
This study confirms that the estimated PNECs we reported recently (4) correspond well to experimentally determined MSCs, at least for tetracycline. Importantly, the selective concentrations we report for tetracycline overlap with those that have been reported in sewage treatment plants (5). We also see that tetracycline not only selects for tetracycline resistance genes, but also resistance genes against other classes of antibiotics, including sulfonamides, beta-lactams and aminoglycosides. Finally, the approach we describe can be used for improved in risk assessment for (also other) antibiotics, and to refine the emission limits we suggested in a recent paper based on theoretical calculations (4).
References and notes
- Okay, seriously: how can a journal’s production team return the proofs for a paper within 24 hours of acceptance, and then wait literally five weeks before putting the final proofs online? Nothing against STOTEN, but I honestly wonder what was going on beyond the scenes here.
- Lundström SV, Östman M, Bengtsson-Palme J, Rutgersson C, Thoudal M, Sircar T, Blanck H, Eriksson KM, Tysklind M, Flach C-F, Larsson DGJ: Minimal selective concentrations of tetracycline in complex aquatic bacterial biofilms. Science of the Total Environment, 553, 587–595 (2016). doi: 10.1016/j.scitotenv.2016.02.103 [Paper link]
- Ågerstrand M, Berg C, Björlenius B, Breitholtz M, Brunstrom B, Fick J, Gunnarsson L, Larsson DGJ, Sumpter JP, Tysklind M, Rudén C: Improving environmental risk assessment of human pharmaceuticals. Environmental Science and Technology (2015). doi:10.1021/acs.est.5b00302
- Bengtsson-Palme J, Larsson DGJ: Concentrations of antibiotics predicted to select for resistant bacteria: Proposed limits for environmental regulation. Environment International, 86, 140-149 (2016). doi: 10.1016/j.envint.2015.10.015
- Michael I, Rizzo L, McArdell CS, Manaia CM, Merlin C, Schwartz T, Dagot C, Fatta-Kassinos D: Urban wastewater treatment plants as hotspots for the release of antibiotics in the environment: a review. Water Research, 47, 957–995 (2013). doi:10.1016/j.watres.2012.11.027
The main problem with the previous version of TriMetAss was that the Trinity developers had changed many options in the Trinity software, which rendered more recent versions of Trinity incompatible with TriMetAss. TriMetAss was not the only external software using Trinity that was affected by these changes. As far as my testing goes, these incompatibilities should now be fixed, by improved Trinity version determination in TriMetAss. This is still not a guarantee for future changes though, so just to make sure, use one of the Trinity versions tested with TriMetAss (versions v2.1.1 or trinityrnaseq_r2013_08_14).
This time I would like to thank Artemis Louyakis at the Univesity of Florida and Tatsuya Unno at the Jeju National University (Korea) for their input on TriMetAss.