- Volume 7, Issue 8, 2021

Many foundation species in chemosynthesis-based ecosystems rely on environmentally acquired symbiotic bacteria for their survival. Hence, understanding the biogeographic distributions of these symbionts at regional scales is key to understanding patterns of connectivity and predicting resilience of their host populations (and thus whole communities). However, such assessments are challenging because they necessitate measuring bacterial genetic diversity at fine resolutions. For this purpose, the recently discovered clustered regularly interspaced short palindromic repeats (CRISPR) constitutes a promising new genetic marker. These DNA sequences harboured by about half of bacteria hold their viral immune memory, and as such, might allow discrimination of different lineages or strains of otherwise indistinguishable bacteria. In this study, we assessed the potential of CRISPR as a hypervariable phylogenetic marker in the context of a population genetic study of an uncultured bacterial species. We used high-throughput CRISPR-based typing along with multi-locus sequence analysis (MLSA) to characterize the regional population structure of the obligate but environmentally acquired symbiont species Candidatus Endoriftia persephone on the Juan de Fuca Ridge. Mixed symbiont populations of Ca. Endoriftia persephone were sampled across individual Ridgeia piscesae hosts from contrasting habitats in order to determine if environmental conditions rather than barriers to connectivity are more important drivers of symbiont diversity. We showed that CRISPR revealed a much higher symbiont genetic diversity than the other housekeeping genes. Several lines of evidence imply this diversity is indicative of environmental strains. Finally, we found with both CRISPR and gene markers that local symbiont populations are strongly differentiated across sites known to be isolated by deep-sea circulation patterns. This research showed the high power of CRISPR to resolve the genetic structure of uncultured bacterial populations and represents a step towards making keystone microbial species an integral part of conservation policies for upcoming mining operations on the seafloor.

Using a previously described metagenomics dataset of 27 billion reads, we reconstructed over 50 000 metagenome-assembled genomes (MAGs) of organisms resident in the porcine gut, 46.5 % of which were classified as >70 % complete with a <10 % contamination rate, and 24.4 % were nearly complete genomes. Here, we describe the generation and analysis of those MAGs using time-series samples. The gut microbial communities of piglets appear to follow a highly structured developmental programme in the weeks following weaning, and this development is robust to treatments including an intramuscular antibiotic treatment and two probiotic treatments. The high resolution we obtained allowed us to identify specific taxonomic ‘signatures’ that characterize the gut microbial development immediately after weaning. Additionally, we characterized the carbohydrate repertoire of the organisms resident in the porcine gut. We tracked the abundance shifts of 294 carbohydrate active enzymes, and identified the species and higher-level taxonomic groups carrying each of these enzymes in their MAGs. This knowledge can contribute to the design of probiotics and prebiotic interventions as a means to modify the piglet gut microbiome.

Ethnicity is consistently reported as a strong determinant of human gut microbiota. However, the bulk of these studies are from Western countries, where microbiota variations are mainly driven by relatively recent migration events. Malaysia is a multicultural society, but differences in gut microbiota persist across ethnicities. We hypothesized that migrant ethnic groups continue to share fundamental gut traits with the population in the country of origin due to shared cultural practices despite subsequent geographical separation. To test this hypothesis, the 16S rRNA gene amplicons from 16 studies comprising three major ethnic groups in Malaysia were analysed, covering 636 Chinese, 248 Indian and 123 Malay individuals from four countries (China, India, Indonesia and Malaysia). A confounder-adjusted permutational multivariate analysis of variance (PERMANOVA) detected a significant association between ethnicity and the gut microbiota (PERMANOVA R 2=0.005, pseudo-F=2.643, P=0.001). A sparse partial least squares – discriminant analysis model trained using the gut microbiota of individuals from China, India and Indonesia (representation of Chinese, Indian and Malay ethnic group, respectively) showed a better-than-random performance in classifying Malaysian of Chinese descent, although the performance for Indian and Malay were modest (true prediction rate, Chinese=0.60, Indian=0.49, Malay=0.44). Separately, differential abundance analysis singled out Ligilactobacillus as being elevated in Indians. We postulate that despite the strong influence of geographical factors on the gut microbiota, cultural similarity due to a shared ethnic origin drives the presence of a shared gut microbiota composition. The interplay of these factors will likely depend on the circumstances of particular groups of migrants.

Despite that obesity is associated with many metabolic diseases, a significant proportion (10–30 %) of obese individuals is recognized as ‘metabolically healthy obeses’ (MHOs). The aim of the current study is to characterize the gut microbiome for MHOs as compared to ‘metabolically unhealthy obeses’ (MUOs). We compared the gut microbiome of 172 MHO and 138 MUO individuals from Chongqing (China) (inclined to eat red meat and food with a spicy taste), and performed validation with selected biomarkers in 40 MHOs and 33 MUOs from Quanzhou (China) (inclined to eat seafood and food with a light/bland taste). The genera Alistipes , Faecalibacterium and Odoribacter had increased abundance in both Chongqing and Quanzhou MHOs. We also observed different microbial functions in MUOs compared to MHOs, including an increased abundance of genes associated with glycan biosynthesis and metabolism. In addition, the microbial gene markers identified from the Chongqing cohort bear a moderate accuracy [AUC (area under the operating characteristic curve)=0.69] for classifying MHOs distinct from MUOs in the Quanzhou cohort. These findings indicate that gut microbiome is significantly distinct between MHOs and MUOs, implicating the potential of the gut microbiome in stratification and refined management of obesity.

The activity of transposable elements (TEs) can be an important driver of genetic diversity with TE-mediated mutations having a wide range of fitness consequences. To avoid deleterious effects of TE activity, some fungi have evolved highly sophisticated genomic defences to reduce TE proliferation across the genome. Repeat-induced point mutation (RIP) is a fungal-specific TE defence mechanism efficiently targeting duplicated sequences. The rapid accumulation of RIPs is expected to deactivate TEs over the course of a few generations. The evolutionary dynamics of TEs at the population level in a species with highly repressive genome defences is poorly understood. Here, we analyse 366 whole-genome sequences of Parastagonospora nodorum, a fungal pathogen of wheat with efficient RIP. A global population genomics analysis revealed high levels of genetic diversity and signs of frequent sexual recombination. Contrary to expectations for a species with RIP, we identified recent TE activity in multiple populations. The TE composition and copy numbers showed little divergence among global populations regardless of the demographic history. Miniature inverted-repeat transposable elements (MITEs) and terminal repeat retrotransposons in miniature (TRIMs) were largely underlying recent intra-species TE expansions. We inferred RIP footprints in individual TE families and found that recently active, high-copy TEs have possibly evaded genomic defences. We find no evidence that recent positive selection acted on TE-mediated mutations rather that purifying selection maintained new TE insertions at low insertion frequencies in populations. Our findings highlight the complex evolutionary equilibria established by the joint action of TE activity, selection and genomic repression.

Unlike most   Salmonella enterica   subsp.   diarizonae  , which are predominantly associated with cold-blooded animals such as reptiles, the serovar IIIb 61:k:1,5,(7) (termed SASd) is regarded as host-adapted to sheep. The bacterium is rarely associated with disease in humans but, nevertheless, SASd isolates are sporadically obtained from human clinical samples. It is unclear whether these transmissions are directly linked to sheep or whether transmissions may, for example, occur through other domestic animals also carrying SASd. For this reason, we utilized whole-genome sequencing to investigate a set of 119 diverse SASd isolates, including sheep-associated and human-associated isolates, as well as isolates obtained from other matrices. We discovered that serovar IIIb 61:k:1,5,(7) is composed of two distinct lineages defined by their sequence types ST432 and ST439. These two lineages are distinguished by a number of genetic features, as well as their prevalence and reservoir. ST432 appears to be the more prevalent sequence type, with the majority of isolates investigated in this study belonging to ST432. In contrast, only a small number of isolates were attributed to ST439. While ST432 isolates were of sheep, human or other origin, all ST439 isolates with source information available, were obtained from human clinical samples. Regarding their genetic features, lineage ST432 shows increased pseudogenization, harbours a virB/D4 plasmid that encodes a type IV secretion system (T4SS) and does not possess the iro gene cluster, which encodes a salmochelin siderophore for iron acquisition. These characteristics likely contribute to the ability of ST432 to persistently colonize the intestines of sheep. Furthermore, we found isolates of the lineage ST432 to be highly clonal, with little variation over the sampling period of almost 20 years. We conclude from the genomic comparisons that SASd underlies a microevolutionary process and that it is specifically lineage ST432 that should be considered as host-adapted to sheep.

Whole-genome sequencing is becoming the de facto standard for bacterial outbreak surveillance and infection prevention. This is accompanied by a variety of bioinformatic tools and needs bioinformatics expertise for implementation. However, little is known about the concordance of reported outbreaks when using different bioinformatic workflows. In this multi-centre proficiency testing among 13 major Dutch healthcare-affiliated centres, bacterial whole-genome outbreak analysis was assessed. Centres who participated obtained two randomized bacterial datasets of Illumina sequences, a Klebsiella pneumoniae and a Vancomycin-resistant Enterococcus faecium, and were asked to apply their bioinformatic workflows. Centres reported back on antimicrobial resistance, multi-locus sequence typing (MLST), and outbreak clusters. The reported clusters were analysed using a method to compare landscapes of phylogenetic trees and calculating Kendall–Colijn distances. Furthermore, fasta files were analysed by state-of-the-art single nucleotide polymorphism (SNP) analysis to mitigate the differences introduced by each centre and determine standardized SNP cut-offs. Thirteen centres participated in this study. The reported outbreak clusters revealed discrepancies between centres, even when almost identical bioinformatic workflows were used. Due to stringent filtering, some centres failed to detect extended-spectrum beta-lactamase genes and MLST loci. Applying a standardized method to determine outbreak clusters on the reported de novo assemblies, did not result in uniformity of outbreak-cluster composition among centres.

Klebsiella pneumoniae has been implicated in wide-ranging nosocomial outbreaks, causing severe infections without effective treatments due to antibiotic resistance. Here, we performed genome sequencing of 70 extensively drug resistant clinical isolates, collected from Brasília’s hospitals (Brazil) between 2010 and 2014. The majority of strains (60 out of 70) belonged to a single clonal complex (CC), CC258, which has become distributed worldwide in the last two decades. Of these CC258 strains, 44 strains were classified as sequence type 11 (ST11) and fell into two distinct clades, but no ST258 strains were found. These 70 strains had a pan-genome size of 10 366 genes, with a core-genome size of ~4476 genes found in 95 % of isolates. Analysis of sequences revealed diverse mechanisms of resistance, including production of multidrug efflux pumps, enzymes with the same target function but with reduced or no affinity to the drug, and proteins that protected the drug target or inactivated the drug. β-Lactamase production provided the most notable mechanism associated with K. pneumoniae . Each strain presented two or three different β-lactamase enzymes, including class A (SHV, CTX-M and KPC), class B and class C AmpC enzymes, although no class D β-lactamase was identified. Strains carrying the NDM enzyme involved three different ST types, suggesting that there was no common genetic origin.

Bacillus anthracis, the causative agent of anthrax disease, is a worldwide threat to livestock, wildlife and public health. While analyses of genetic data from across the globe have increased our understanding of this bacterium’s population genomic structure, the influence of selective pressures on this successful pathogen is not well understood. In this study, we investigate the effects of antimicrobial resistance, phage diversity, geography and isolation source in shaping population genomic structure. We also identify a suite of candidate genes potentially under selection, driving patterns of diversity across 356 globally extant   B. anthracis   genomes. We report ten antimicrobial resistance genes and 11 different prophage sequences, resulting in the first large-scale documentation of these genetic anomalies for this pathogen. Results of random forest classification suggest genomic structure may be driven by a combination of antimicrobial resistance, geography and isolation source, specific to the population cluster examined. We found strong evidence that a recombination event linked to a gene involved in protein synthesis may be responsible for phenotypic differences between comparatively disparate populations. We also offer a list of genes for further examination of   B. anthracis   evolution, based on high-impact single nucleotide polymorphisms (SNPs) and clustered mutations. The information presented here sheds new light on the factors driving genomic structure in this notorious pathogen and may act as a road map for future studies aimed at understanding functional differences in terms of   B. anthracis   biogeography, virulence and evolution.

The evolution and emergence of drug-resistant tuberculosis (TB) has been studied extensively in some contexts, but the ecological drivers of these two processes remain poorly understood. This study sought to describe the joint evolutionary and epidemiological histories of a novel multidrug-resistant Mycobacterium tuberculosis strain recently identified in the capital city of the Republic of Moldova (MDR Ural/4.2), where genomic surveillance of drug-resistant M. tuberculosis has been limited thus far. Using whole genome sequence data and Bayesian phylogenomic methods, we reconstruct the stepwise acquisition of drug resistance mutations in the MDR Ural/4.2 strain, estimate its historical bacterial population size over time, and infer the migration history of this strain between Eastern European countries. We infer that MDR Ural/4.2 likely evolved (via acquisition of rpoB S450L, which confers resistance to rifampin) in the early 1990s, during a period of social turmoil following Moldovan independence from the Soviet Union. This strain subsequently underwent substantial population size expansion in the early 2000s, at a time when national guidelines encouraged inpatient treatment of TB patients. We infer exportation of this strain and its isoniazid-resistant ancestral precursor from Moldova to neighbouring countries starting as early as 1985. Our findings suggest temporal and ecological associations between specific public health practices, including inpatient hospitalization of drug-resistant TB cases from the early 2000s until 2013, and the evolution of drug-resistant M. tuberculosis in Moldova. These findings underscore the need for regional coordination in TB control and expanded genomic surveillance efforts across Eastern Europe.

The emergence of antimicrobial resistance (AMR) to first- and second-line treatment regimens of enteric fever is a global public-health problem, and routine genomic surveillance to inform clinical and public-health management guidance is essential. Here, we present the prospective analysis of genomic data to monitor trends in incidence, AMR and travel, and assess hierarchical clustering (HierCC) methodology of 1742 isolates of typhoidal salmonellae. Trend analysis of Salmonella Typhi and S. Paratyphi A cases per year increased 48 and 17.3%, respectively, between 2016 and 2019 in England, mainly associated with travel to South Asia. S. Paratyphi B cases have remained stable and are mainly associated with travel to the Middle East and South America. There has been an increase in the number of S. Typhi exhibiting a multidrug-resistant (MDR) profile and the emergence of extensively drug resistant (XDR) profiles. HierCC was a robust method to categorize clonal groups into clades and clusters associated with travel and AMR profiles. The majority of cases that had XDR S. Typhi reported recent travel to Pakistan (94 %) and belonged to a subpopulation of the 4.3.1 (H58) clone (HC5_1452). The phenotypic and genotypic AMR results showed high concordance for S. Typhi and S. Paratyphi A, B and C, with 99.99 % concordance and only three (0.01 %) discordant results out of a possible 23 178 isolate/antibiotic combinations. Genomic surveillance of enteric fever has shown the recent emergence and increase of MDR and XDR S. Typhi strains, resulting in a review of clinical guidelines to improve management of imported infections.

Shigellosis in men who have sex with men (MSM) is caused by multidrug resistant Shigellae, exhibiting resistance to antimicrobials including azithromycin, ciprofloxacin and more recently the third-generation cephalosporins. We sequenced four bla CTX-M-27-positive MSM Shigella isolates (2018–20) using Oxford Nanopore Technologies; three S. sonnei (identified as two MSM clade 2, one MSM clade 5) and one S. flexneri 3a, to explore AMR context. All S. sonnei isolates harboured Tn7/Int2 chromosomal integrons, whereas S. flexneri 3a contained the Shigella Resistance Locus. All strains harboured IncFII pKSR100-like plasmids (67-83kbp); where present bla CTX-M-27 was located on these plasmids flanked by IS26 and IS903B, however bla CTX-M-27 was lost in S. flexneri 3a during storage between Illumina and Nanopore sequencing. IncFII AMR regions were mosaic and likely reorganised by IS26; three of the four plasmids contained azithromycin-resistance genes erm(B) and mph(A) and one harboured the pKSR100 integron. Additionally, all S. sonnei isolates possessed a large IncB/O/K/Z plasmid, two of which carried aph(3’)-Ib/aph(6)-Id/sul2 and tet(A). Monitoring the transmission of mobile genetic elements with co-located AMR determinants is necessary to inform empirical treatment guidance and clinical management of MSM-associated shigellosis.

Escherichia coli is a ubiquitous bacterium that has been widely exposed to antibiotics over the last 70 years. It has adapted by acquiring different antibiotic-resistance genes (ARGs), the census of which we aim to characterize here. To do so, we analysed 70 301 E. coli genomes obtained from the EnteroBase database and detected 1 027 651 ARGs using the AMRFinder, Mustard and ResfinderFG ARG databases. We observed a strong phylogroup and clonal lineage specific distribution of some ARGs, supporting the argument for epistasis between ARGs and the strain genetic background. However, each phylogroup had ARGs conferring a similar antibiotic class resistance pattern, indicating phenotypic adaptive convergence. The G+C content or the type of ARG was not associated with the frequency of the ARG in the database. In addition, we identified ARGs from anaerobic, non- Proteobacteria bacteria in four genomes of E. coli , supporting the hypothesis that the transfer between anaerobic bacteria and E. coli can spontaneously occur but remains exceptional. In conclusion, we showed that phylum barrier and intra-species phylogenetic history are major drivers of the acquisition of a resistome in E. coli .

During industrial processes, yeasts are exposed to harsh conditions, which eventually lead to adaptation of the strains. In the laboratory, it is possible to use experimental evolution to link the evolutionary biology response to these adaptation pressures for the industrial improvement of a specific yeast strain. In this work, we aimed to study the adaptation of a wine industrial yeast in stress conditions of the high ethanol concentrations present in stopped fermentations and secondary fermentations in the processes of champagne production. We used a commercial Saccharomyces cerevisiae × S. uvarum hybrid and assessed its adaptation in a modified synthetic must (M-SM) containing high ethanol, which also contained metabisulfite, a preservative that is used during wine fermentation as it converts to sulfite. After the adaptation process under these selected stressful environmental conditions, the tolerance of the adapted strain (H14A7-etoh) to sulfite and ethanol was investigated, revealing that the adapted hybrid is more resistant to sulfite compared to the original H14A7 strain, whereas ethanol tolerance improvement was slight. However, a trade-off in the adapted hybrid was found, as it had a lower capacity to ferment glucose and fructose in comparison with H14A7. Hybrid genomes are almost always unstable, and different signals of adaptation on H14A7-etoh genome were detected. Each subgenome present in the adapted strain had adapted differently. Chromosome aneuploidies were present in S. cerevisiae chromosome III and in S. uvarum chromosome VII–XVI, which had been duplicated. Moreover, S. uvarum chromosome I was not present in H14A7-etoh and a loss of heterozygosity (LOH) event arose on S. cerevisiae chromosome I. RNA-sequencing analysis showed differential gene expression between H14A7-etoh and H14A7, which can be easily correlated with the signals of adaptation that were found on the H14A7-etoh genome. Finally, we report alterations in the lipid composition of the membrane, consistent with conserved tolerance mechanisms.

Ochrobactrum pseudogrignonense (newly named Brucella pseudogrignonensis ) is an emerging pathogen in immunodeficient and immunocompetent patients. Most documented cases associated with Ochrobactrum are frequently catheter-related and exhibit wide-spectrum β-lactam resistance. Misidentification of this pathogen using commercial bacterial identification kits is common. We identified a case of O. pseudogrignonense infection associated with cholelithiasis. The O. pseudogrignonense genome was sequenced and reconstructed using a Nanopore and Illumina hybrid strategy. A novel bla OXA-919 divergent from existing OXA members was identified and subsequent analysis revealed its existence in all available O. pseudogrignonense genomes, which forms a new phylogenetic subgroup distinct from other OXA clusters. Further analysis demonstrated the presence of the novel bla OXA-919 in the chromosome of several other Ochrobactrum species. Our study indicated that Ochrobactrum chromosomes may be a reservoir of bla OXA-919 β-lactamases.

Trypanosoma brucei evansi and T. brucei equiperdum are animal infective trypanosomes conventionally classified by their clinical disease presentation, mode of transmission, host range, kinetoplast DNA (kDNA) composition and geographical distribution. Unlike other members of the subgenus Trypanozoon, they are non-tsetse transmitted and predominantly morphologically uniform (monomorphic) in their mammalian host. Their classification as independent species or subspecies has been long debated and genomic studies have found that isolates within T. brucei evansi and T. brucei equiperdum have polyphyletic origins. Since current taxonomy does not fully acknowledge these polyphyletic relationships, we re-analysed publicly available genomic data to carefully define each clade of monomorphic trypanosome. This allowed us to identify, and account for, lineage-specific variation. We included a recently published isolate, IVM-t1, which was originally isolated from the genital mucosa of a horse with dourine and typed as T. equiperdum. Our analyses corroborate previous studies in identifying at least four distinct monomorphic T. brucei clades. We also found clear lineage-specific variation in the selection efficacy and heterozygosity of the monomorphic lineages, supporting their distinct evolutionary histories. The inferred evolutionary position of IVM-t1 suggests its reassignment to the T. brucei evansi type B clade, challenging the relationship between the Trypanozoon species, the infected host, mode of transmission and the associated pathological phenotype. The analysis of IVM-t1 also provides, to our knowledge, the first evidence of the expansion of T. brucei evansi type B, or a fifth monomorphic lineage represented by IVM-t1, outside of Africa, with important possible implications for disease diagnosis.