Gene duplication drives genome expansion in a major lineage of Thaumarchaeota

Acknowledgements This work and P.O.S. were financially supported by UKRI through the NERC grant NE/R001529/1. In addition, C.G.-R. and T.A.W. were both supported by a Royal Society University Research Fellowship (URF150571 and UF140626). C.Q. was funded through an MRC fellowship (MR/M50161X/1) as part of the CLoud Infrastructure for Microbial Genomics (CLIMB) consortium (MR/L015080/1). S.R. was funded through the BBSRC grant BB/R015171/1. The Thames Metagenome Database was funded through the NERC grants NE/M011674/1, NE/M011259/1 and NE/M01133X/1. We thank Dr Tony Travis for his support with Biolinux and acknowledge Prof Jim Prosser for his critical reading of the manuscript. The authors would also like to acknowledge the support of the Maxwell computer cluster funded by the University of Aberdeen.Peer reviewedPublisher PD

Gene duplication drives genome expansion in a major lineage of Thaumarchaeota

Ammonia-oxidising archaea of the phylum Thaumarchaeota are important organisms in the nitrogen cycle, but the mechanisms driving their radiation into diverse ecosystems remain underexplored. Here, existing thaumarchaeotal genomes are complemented with 12 genomes belonging to the previously under-sampled Nitrososphaerales to investigate the impact of lateral gene transfer (LGT), gene duplication and loss across thaumarchaeotal evolution. We reveal a major role for gene duplication in driving genome expansion subsequent to early LGT. In particular, two large LGT events are identified into Nitrososphaerales and the fate of these gene families is highly lineage-specific, being lost in some descendant lineages, but undergoing extensive duplication in others, suggesting niche-specific roles. Notably, some genes involved in carbohydrate transport or coenzyme metabolism were duplicated, likely facilitating niche specialisation in soils and sediments. Overall, our results suggest that LGT followed by gene duplication drives Nitrososphaerales evolution, highlighting a previously under-appreciated mechanism of genome expansion in archaea

Impact of sulfamethoxazole on a riverine microbiome

The continued emergence of bacterial pathogens presenting antimicrobial resistance is widely recognised as a global health threat and recent attention focused on potential environmental reservoirs of antibiotic resistance genes (ARGs). Freshwater environments such as rivers represent a potential hotspot for ARGs and antibiotic resistant bacteria as they are receiving systems for effluent discharges from wastewater treatment plants (WWTPs). Effluent also contains low levels of different antimicrobials including antibiotics and biocides. Sulfonamides are antibacterial chemicals widely used in clinical, veterinary and agricultural settings and are frequently detected in sewage sludge and manure in addition to riverine ecosystems. The impact of such exposure on ARG prevalence and diversity is unknown, so the aim of this study was to investigate the release of a sub-lethal concentration of the sulfonamide compound sulfamethoxazole (SMX) on the river bacterial microbiome using a flume system. This system was a semi-natural in vitro flume using river water (30 L) and sediment (6 kg) with circulation to mimic river flow. A combination of ‘omics’ approaches were conducted to study the impact of SMX exposure on the microbiomes within the flumes. Metagenomic analysis showed that the addition of low concentrations of SMX (<4 μg L−1) had a limited effect on the bacterial resistome in the water fraction only, with no impact observed in the sediment. Metaproteomics did not show differences in ARGs expression with SMX exposure in water. Overall, the river bacterial community was resilient to short term exposure to sub-lethal concentrations of SMX which mimics the exposure such communities experience downstream of WWTPs throughout the year

De novo extraction of microbial strains from metagenomes reveals intra-species niche partitioning

Background: We introduce DESMAN for De novo Extraction of Strains from MetAgeNomes. Metagenome sequencing generates short reads from throughout the genomes of a microbial community. Increasingly large, multi-sample metagenomes, stratified in space and time are being generated from communities with thousands of species. Repeats result in fragmentary co-assemblies with potentially millions of contigs. Contigs can be binned into metagenome assembled genomes (MAGs) but strain level variation will remain. DESMAN identifies variants on core genes, then uses co-occurrence across samples to link variants into strain sequences and abundance profiles. These strain profiles are then searched for on non-core genes to determine the accessory genes present in each strain. Results: We validated DESMAN on a synthetic twenty genome community with 64 samples. We could resolve the five E. coli strains present with 99.58% accuracy across core gene variable sites and their gene complement with 95.7% accuracy. Similarly, on real fecal metagenomes from the 2011 E. coli (STEC) O104:H4 outbreak, the outbreak strain was reconstructed with 99.8% core sequence accuracy. Application to an anaerobic digester metagenome time series reveals that strain level variation is endemic with 16 out of 26 MAGs (61.5%) examined exhibiting two strains. In almost all cases the strain proportions were not statistically different between replicate reactors, suggesting intra-species niche partitioning. The only exception being when the two strains had almost identical gene complement and, hence, functional capability. Conclusions: DESMAN will provide a provide a powerful tool for de novo resolution of fine-scale variation in microbial communities. It is available as open source software from https://github.com/chrisquince/DESMAN