CYP8B1 downregulation mediates the metabolic effects of vertical sleeve gastrectomy in mice.

BACKGROUND AND AIMS: Although the benefits of vertical sleeve gastrectomy (VSG) surgery are well known, the molecular mechanisms by which VSG alleviates obesity and its complications remain unclear. We aim to determine the role of CYP8B1 (cytochrome P450, family 8, subfamily B, polypeptide 1) in mediating the metabolic benefits of VSG. APPROACH AND RESULTS: We found that expression of CYP8B1, a key enzyme in controlling the 12α-hydroxylated (12α-OH) bile acid (BA) to non-12α-OH BA ratio, was strongly downregulated after VSG. Using genetic mouse models of CYP8B1 overexpression, knockdown, and knockout, we demonstrated that overexpression of CYP8B1 dampened the metabolic improvements associated with VSG. In contrast, short hairpin RNA-mediated CYP8B1 knockdown improved metabolism similar to those observed after VSG. Cyp8b1 deficiency diminished the metabolic effects of VSG. Further, VSG-induced alterations to the 12α-OH/non-12α-OH BA ratio in the BA pool depended on CYP8B1 expression level. Consequently, intestinal lipid absorption was restricted, and the gut microbiota (GM) profile was altered. Fecal microbiota transplantation from wild type-VSG mice (vs. fecal microbiota transplantation from wild-type-sham mice) improved metabolism in recipient mice, while there were no differences between mice that received fecal microbiota transplantation from knockout-sham and knockout-VSG mice. CONCLUSIONS: CYP8B1 is a critical downstream target of VSG. Modulation of BA composition and gut microbiota profile by targeting CYP8B1 may provide novel insight into the development of therapies that noninvasively mimic bariatric surgery to treat obesity and its complications.

Comparison of Brucella abortus population structure based on genotyping methods with different levels of resolution.

Numerous genotyping techniques based on different principles and with different costs and levels of resolution are currently available for understanding the transmission dynamics of brucellosis worldwide. We aimed to compare the population structure of the genomes of 53 Brazilian Brucella abortus isolates using eight different genotyping methods: multiple-locus variable-number tandem-repeat analysis (MLVA8, MLVA11, MLVA16), multilocus sequence typing (MLST9, MLST21), core genome MLST (cgMLST) and two techniques based on single nucleotide polymorphism (SNP) detection (parSNP and NASP) from whole genomes. The strains were isolated from six different Brazilian states between 1977 and 2008 and had previously been analyzed using MLVA8, MLVA11, and MLVA16. Their whole genomes were sequenced, assembled, and subjected to MLST9 MLST21, cgMLST, and SNP analyses. All the genotypes were compared by hierarchical grouping method based on the average distances between the correlation matrices of each technique. MLST9 and MLST21 had the lowest level of resolution, both revealing only four genotypes. MLVA8, MLVA11, and MLVA16 had progressively increasing levels of resolution as more loci were analyzed, identifying 6, 16, and 44 genotypes, respectively. cgMLST showed the highest level of resolution, identifying 45 genotypes, followed by the SNP-based methods, both of which had 44 genotypes. In the assessed population, MLVA was more discriminatory than MLST and was easier and cheaper to perform. SNP techniques and cgMLST provided the highest levels of resolution and the results from the two methods were in close agreement. In conclusion, the choice of genotyping technique can strongly affect one's ability to make meaningful epidemiological conclusions but is dependent on available resources: while the VNTR based techniques are more indicated to high prevalence scenarios, the WGS methods are the ones with the best discriminative power and therefore recommended for outbreaks investigation.

Multi-faceted metagenomic analysis of spacecraft associated surfaces reveal planetary protection relevant microbial composition.

The National Aeronautics and Space Administration (NASA) has been monitoring the microbial burden of spacecraft since the 1970's Viking missions. Originally culture-based and then focused 16S sequencing techniques were used, but we have now applied whole metagenomic sequencing to a variety of cleanroom samples at the Jet Propulsion Lab (JPL), including the Spacecraft Assembly Facility (SAF) with the goals of taxonomic identification and for functional assignment. Our samples included facility pre-filters, cleanroom vacuum debris, and surface wipes. The taxonomic composition was carried out by three different analysis tools to contrast marker, k-mer, and true alignment approaches. Hierarchical clustering analysis of the data separated vacuum particles from other SAF DNA samples. Vacuum particle samples were the most diverse while DNA samples from the ISO (International Standards Organization) compliant facilities and the SAF were the least diverse; all three were dominated by Proteobacteria. Wipe samples had higher diversity and were predominated by Actinobacteria, including human commensals Cutibacterium acnes and Corynebacterium spp. Taxa identified by the three methods were not identical, supporting the use of multiple methods for metagenome characterization. Likewise, functional annotation was performed using multiple methods. Vacuum particles and SAF samples contained strong signals of the tricarboxylic acid cycle and of amino acid biosynthesis, suggesting that many of the identified microorganisms have the ability to grow in nutrient-limited environments. In total, 18 samples generated high quality metagenome assembled genomes (MAG), which were dominated by Moraxella osloensis or Malassezia restricta. One M. osloensis MAG was assembled into a single circular scaffold and gene annotated. This study includes a rigorous quantitative determination of microbial loads and a qualitative dissection of microbial composition. Assembly of multiple specimens led to greater confidence for the identification of particular species and their predicted functional roles.

Analysis of Gut Microbiome Using Explainable Machine Learning Predicts Risk of Diarrhea Associated With Tyrosine Kinase Inhibitor Neratinib: A Pilot Study.

Neratinib has great efficacy in treating HER2+ breast cancer but is associated with significant gastrointestinal toxicity. The objective of this pilot study was to understand the association of gut microbiome and neratinib-induced diarrhea. Twenty-five patients (age ≥ 60) were enrolled in a phase II trial evaluating safety and tolerability of neratinib in older adults with HER2+ breast cancer (NCT02673398). Fifty stool samples were collected from 11 patients at baseline and during treatment. 16S rRNA analysis was performed and relative abundance data were generated. Shannon's diversity was calculated to examine gut microbiome dysbiosis. An explainable tree-based approach was utilized to classify patients who might experience neratinib-related diarrhea (grade ≥ 1) based on pre-treatment baseline microbial relative abundance data. The hold-out Area Under Receiver Operating Characteristic and Area Under Precision-Recall Curves of the model were 0.88 and 0.95, respectively. Model explanations showed that patients with a larger relative abundance of Ruminiclostridium 9 and Bacteroides sp. HPS0048 may have reduced risk of neratinib-related diarrhea and was confirmed by Kruskal-Wallis test (p ≤ 0.05, uncorrected). Our machine learning model identified microbiota associated with reduced risk of neratinib-induced diarrhea and the result from this pilot study will be further verified in a larger study. CLINICAL TRIAL REGISTRATION: ClinicalTrials.gov, identifier NCT02673398.

Randomized trial assessing impact of probiotic supplementation on gut microbiome and clinical outcome from targeted therapy in metastatic renal cell carcinoma.

Studies suggest a link between the gut microbiome and metastatic renal cell carcinoma (mRCC) outcomes, including evidence that mRCC patients possess a lower abundance of Bifidobacterium spp. compared to healthy adults. We sought to assess if a Bifidobacterium-containing yogurt product could modulate the gut microbiome and clinical outcome from vascular endothelial growth factor-tyrosine kinase inhibitors (VEGF-TKIs). mRCC patients initiating VEGF-TKIs, regardless of the line of therapy, were randomized to probiotic-supplemented (two 4 oz. servings of the probiotic yogurt product daily) or probiotic-restricted arms. Stool samples were collected prior to therapy and at weeks 2, 3, 4, and 12. Microbiome composition was assessed using whole-metagenome sequencing. A total of 20 patients were randomized. Bifidobacterium animalis, the active ingredient of the probiotic supplement, reached detectable levels in all patients in the probiotic-supplemented arm versus two patients in the probiotic-restricted arm. Clinical benefit rate was similar in probiotic-supplemented versus probiotic-restricted arms (70% vs. 80%, p = 0.606). Linear discriminant analysis (LDA) effect size analysis of MetaPhIAn2 abundance data predicted 25 enriched species demonstrating an LDA score >3 in either clinical benefit or no clinical benefit. In patients with clinical benefit (vs. no clinical benefit), Barnesiella intestinihominis and Akkermansia muciniphila were significantly more abundant (p = 7.4 × 10-6 and p = 5.6 × 10-3 , respectively). This is the first prospective randomized study demonstrating modulation of the gut microbiome with a probiotic in mRCC. Probiotic supplementation successfully increased the Bifidobacterium spp. levels. Analysis of longitudinal stool specimens identified an association between B. intestinihominis, A. muciniphila, and clinical benefit with therapy. Trial Registration: NCT02944617.

Roles of bacteriophages, plasmids and CRISPR immunity in microbial community dynamics revealed using time-series integrated meta-omics.

Viruses and plasmids (invasive mobile genetic elements (iMGEs)) have important roles in shaping microbial communities, but their dynamic interactions with CRISPR-based immunity remain unresolved. We analysed generation-resolved iMGE-host dynamics spanning one and a half years in a microbial consortium from a biological wastewater treatment plant using integrated meta-omics. We identified 31 bacterial metagenome-assembled genomes encoding complete CRISPR-Cas systems and their corresponding iMGEs. CRISPR-targeted plasmids outnumbered their bacteriophage counterparts by at least fivefold, highlighting the importance of CRISPR-mediated defence against plasmids. Linear modelling of our time-series data revealed that the variation in plasmid abundance over time explained more of the observed community dynamics than phages. Community-scale CRISPR-based plasmid-host and phage-host interaction networks revealed an increase in CRISPR-mediated interactions coinciding with a decrease in the dominant 'Candidatus Microthrix parvicella' population. Protospacers were enriched in sequences targeting genes involved in the transmission of iMGEs. Understanding the factors shaping the fitness of specific populations is necessary to devise control strategies for undesirable species and to predict or explain community-wide phenotypes.

A Phase II Clinical Trial of Pembrolizumab and Enobosarm in Patients with Androgen Receptor-Positive Metastatic Triple-Negative Breast Cancer.

LESSONS LEARNED: The combination of enobosarm and pembrolizumab was well tolerated and showed a modest clinical benefit rate of 25% at 16 weeks. Future trials investigating androgen receptor-targeted therapy in combination with immune checkpoint inhibitors are warranted. BACKGROUND: Luminal androgen receptor is a distinct molecular subtype of triple-negative breast cancer (TNBC) defined by overexpression of androgen receptor (AR). AR-targeted therapy has shown modest activity in AR-positive (AR+) TNBC. Enobosarm (GTx-024) is a nonsteroidal selective androgen receptor modulator (SARM) that demonstrates preclinical and clinical activity in AR+ breast cancer. The current study was designed to explore the safety and efficacy of the combination of enobosarm and pembrolizumab in patients with AR+ metastatic TNBC (mTNBC). METHODS: This study was an open-label phase II study for AR+ (≥10%, 1+ by immunohistochemistry [IHC]) mTNBC. Eligible patients received pembrolizumab 200 mg intravenous (IV) every 3 weeks and enobosarm 18 mg oral daily. The primary objective was to evaluate the safety of enobosarm plus pembrolizumab and determine the response rate. Peripheral blood, tumor biopsies, and stool samples were collected for correlative analysis. RESULTS: The trial was stopped early because of the withdrawal of GTx-024 drug supply. Eighteen patients were enrolled, and 16 were evaluable for responses. Median age was 64 (range 36-81) years. The combination was well tolerated, with only a few grade 3 adverse events: one dry skin, one diarrhea, and one musculoskeletal ache. The responses were 1 of 16 (6%) complete response (CR), 1 of 16 (6%) partial response (PR), 2 of 16 (13%) stable disease (SD), and 12 of 16 (75%) progressive disease (PD). Response rate (RR) was 2 of 16 (13%). Clinical benefit rate (CBR) at 16 weeks was 4 of 16 (25%). Median follow-up was 24.9 months (95% confidence interval [CI], 17.5-30.9). Progression-free survival (PFS) was 2.6 months (95% CI, 1.9-3.1) and overall survival (OS) was 25.5 months (95% CI, 10.4-not reached [NR]). CONCLUSION: The combination of enobosarm and pembrolizumab was well tolerated, with a modest clinical benefit rate of 25% at 16 weeks in heavily pretreated AR+ TNBC without preselected programmed death ligand-1 (PD-L1). Future clinical trials combining AR-targeted therapy with immune checkpoint inhibitor (ICI) for AR+ TNBC warrant investigation.

Integration of time-series meta-omics data reveals how microbial ecosystems respond to disturbance.

The development of reliable, mixed-culture biotechnological processes hinges on understanding how microbial ecosystems respond to disturbances. Here we reveal extensive phenotypic plasticity and niche complementarity in oleaginous microbial populations from a biological wastewater treatment plant. We perform meta-omics analyses (metagenomics, metatranscriptomics, metaproteomics and metabolomics) on in situ samples over 14 months at weekly intervals. Based on 1,364 de novo metagenome-assembled genomes, we uncover four distinct fundamental niche types. Throughout the time-series, we observe a major, transient shift in community structure, coinciding with substrate availability changes. Functional omics data reveals extensive variation in gene expression and substrate usage amongst community members. Ex situ bioreactor experiments confirm that responses occur within five hours of a pulse disturbance, demonstrating rapid adaptation by specific populations. Our results show that community resistance and resilience are a function of phenotypic plasticity and niche complementarity, and set the foundation for future ecological engineering efforts.

Stool Microbiome Profiling of Patients with Metastatic Renal Cell Carcinoma Receiving Anti-PD-1 Immune Checkpoint Inhibitors.

Preclinical models and early clinical data suggest an interplay between the gut microbiome and response to immunotherapy in solid tumors including metastatic renal cell carcinoma (mRCC). We sought to characterize the stool microbiome of mRCC patients receiving a checkpoint inhibitor (CPI) and to assess treatment-related changes in microbiome composition over the course of CPI therapy. Stool was collected from 31 patients before initiation of nivolumab (77%) or nivolumab plus ipilimumab (23%) therapy, of whom 58% experienced clinical benefit. Greater microbial diversity was associated with clinical benefit from CPI therapy (p = 0.001), and multiple species were associated with clinical benefit or lack thereof. Temporal profiling of the microbiome indicated that the relative abundance of Akkermansia muciniphila increased in patients deriving clinical benefit from CPIs. This study substantiates results from previous CPI-related microbiome profiling studies in mRCC. Temporal changes in microbiome composition suggest potential utility in modulating the microbiome for more successful CPI outcomes. PATIENT SUMMARY: We compared the composition and diversity of the gut microbiome in patients receiving immunotherapy for renal cell carcinoma. We found that higher microbial diversity is associated with better treatment outcomes. Treatment response is characterized by changes in microbial species over the course of treatment.

First draft genome sequence of a strain belonging to the Zoogloea genus and its gene expression in situ.

The Gram-negative beta-proteobacterium Zoogloea sp. LCSB751 (LMG 29444) was newly isolated from foaming activated sludge of a municipal wastewater treatment plant. Here, we describe its draft genome sequence and annotation together with a general physiological and genomic analysis, as the first sequenced representative of the Zoogloea genus. Moreover, Zoogloea sp. gene expression in its environment is described using metatranscriptomic data obtained from the same treatment plant. The presented genomic and transcriptomic information demonstrate a pronounced capacity of this genus to synthesize poly-ß-hydroxyalkanoate within wastewater.

Phylogeographic, genomic, and meropenem susceptibility analysis of Burkholderia ubonensis.

The bacterium Burkholderia ubonensis is commonly co-isolated from environmental specimens harbouring the melioidosis pathogen, Burkholderia pseudomallei. B. ubonensis has been reported in northern Australia and Thailand but not North America, suggesting similar geographic distribution to B. pseudomallei. Unlike most other Burkholderia cepacia complex (Bcc) species, B. ubonensis is considered non-pathogenic, although its virulence potential has not been tested. Antibiotic resistance in B. ubonensis, particularly towards drugs used to treat the most severe B. pseudomallei infections, has also been poorly characterised. This study examined the population biology of B. ubonensis, and includes the first reported isolates from the Caribbean. Phylogenomic analysis of 264 B. ubonensis genomes identified distinct clades that corresponded with geographic origin, similar to B. pseudomallei. A small proportion (4%) of strains lacked the 920kb chromosome III replicon, with discordance of presence/absence amongst genetically highly related strains, demonstrating that the third chromosome of B. ubonensis, like other Bcc species, probably encodes for a nonessential pC3 megaplasmid. Multilocus sequence typing using the B. pseudomallei scheme revealed that one-third of strains lack the "housekeeping" narK locus. In comparison, all strains could be genotyped using the Bcc scheme. Several strains possessed high-level meropenem resistance (≥32 µg/mL), a concern due to potential transmission of this phenotype to B. pseudomallei. In silico analysis uncovered a high degree of heterogeneity among the lipopolysaccharide O-antigen cluster loci, with at least 35 different variants identified. Finally, we show that Asian B. ubonensis isolate RF23-BP41 is avirulent in the BALB/c mouse model via a subcutaneous route of infection. Our results provide several new insights into the biology of this understudied species.

Mapping and comparing bacterial microbiota in the sinonasal cavity of healthy, allergic rhinitis, and chronic rhinosinusitis subjects.

BACKGROUND: The role of microbiota in sinonasal inflammation can be further understood by targeted sampling of healthy and diseased subjects. We compared the microbiota of the middle meatus (MM) and inferior meatus (IM) in healthy, allergic rhinitis (AR), and chronic rhinosinusitis (CRS) subjects to characterize intrasubject, intersubject, and intergroup differences. METHODS: Subjects were recruited in the office, and characterized into healthy, AR, and CRS groups. Endoscopically-guided swab samples were obtained from the MM and IM bilaterally. Bacterial microbiota were characterized by sequencing the V3-V4 region of the 16S ribosomal RNA (rRNA) gene. RESULTS: Intersubject microbiome analyses were conducted in 65 subjects: 8 healthy, 11 AR, and 46 CRS (25 CRS with nasal polyps [CRSwNP]; 21 CRS without nasal polyps [CRSsNP]). Intrasubject analyses were conducted for 48 individuals (4 controls, 11 AR, 8 CRSwNP, and 15 CRSwNP). There was considerable intersubject microbiota variability, but intrasubject profiles were similar (p = 0.001, nonparametric t test). Intrasubject bacterial diversity was significantly reduced in MM of CRSsNP subjects compared to IM samples (p = 0.022, nonparametric t test). CRSsNP MM samples were enriched in Streptococcus, Haemophilus, and Fusobacterium spp. but exhibited loss of diversity compared to healthy, CRSwNP, and AR subject-samples (p < 0.05; nonparametric t test). CRSwNP patients were enriched in Staphylococcus, Alloiococcus, and Corynebacterium spp. CONCLUSION: This study presents the sinonasal microbiome profile in one of the larger populations of non-CRS and CRS subjects, and is the first office-based cohort in the literature. In contrast to healthy, AR, and CRSwNP subjects, CRSsNP MM samples exhibited decreased microbiome diversity and anaerobic enrichment. CRSsNP MM samples had reduced diversity compared to same-subject IM samples, a novel finding.

More than 50% of Clostridium difficile Isolates from Pet Dogs in Flagstaff, USA, Carry Toxigenic Genotypes.

Nosocomial acquisition of Clostridium difficile is well documented, yet recent studies have highlighted the importance of community acquired infections and identified community associated reservoirs for this pathogen. Multiple studies have implicated companion pets and farm animals as possible sources of community acquired C. difficile infections in humans. To explore the potential role of pet dogs in human C. difficile infections we systematically collected canine fecal samples (n = 197) in Flagstaff, AZ. Additionally, nineteen fecal samples were collected at a local veterinary clinic from diarrheic dogs. We used these combined samples to investigate important questions regarding C. difficile colonization in pet canines: 1) What is the prevalence and diversity of C. difficile in this companion pet population, and 2) Do C. difficile isolates collected from canines genetically overlap with isolates that cause disease in humans? We used a two-step sequence typing approach, including multilocus sequence typing to determine the overall genetic diversity of C. difficile present in Flagstaff canines, and whole-genome sequencing to assess the fine-scale diversity patterns within identical multilocus sequence types from isolates obtained within and among multiple canine hosts. We detected C. difficile in 17% of the canine fecal samples with 10% containing toxigenic strains that are known to cause human disease. Sequencing analyses revealed similar genotypes in dogs and humans. These findings suggest that companion pets are a potential source of community acquired C. difficile infections in humans.

A Bacillus anthracis Genome Sequence from the Sverdlovsk 1979 Autopsy Specimens.

Anthrax is a zoonotic disease that occurs naturally in wild and domestic animals but has been used by both state-sponsored programs and terrorists as a biological weapon. A Soviet industrial production facility in Sverdlovsk, USSR, proved deficient in 1979 when a plume of spores was accidentally released and resulted in one of the largest known human anthrax outbreaks. In order to understand this outbreak and others, we generated a Bacillus anthracis population genetic database based upon whole-genome analysis to identify all single-nucleotide polymorphisms (SNPs) across a reference genome. Phylogenetic analysis has defined three major clades (A, B, and C), B and C being relatively rare compared to A. The A clade has numerous subclades, including a major polytomy named the trans-Eurasian (TEA) group. The TEA radiation is a dominant evolutionary feature of B. anthracis, with many contemporary populations having resulted from a large spatial dispersal of spores from a single source. Two autopsy specimens from the Sverdlovsk outbreak were deep sequenced to produce draft B. anthracis genomes. This allowed the phylogenetic placement of the Sverdlovsk strain into a clade with two Asian live vaccine strains, including the Russian Tsiankovskii strain. The genome was examined for evidence of drug resistance manipulation or other genetic engineering, but none was found. The Soviet Sverdlovsk strain genome is consistent with a wild-type strain from Russia that had no evidence of genetic manipulation during its industrial production. This work provides insights into the world's largest biological weapons program and provides an extensive B. anthracis phylogenetic reference. IMPORTANCE: The 1979 Russian anthrax outbreak resulted from an industrial accident at the Soviet anthrax spore production facility in the city of Sverdlovsk. Deep genomic sequencing of two autopsy specimens generated a draft genome and phylogenetic placement of the Soviet Sverdlovsk anthrax strain. While it is known that Soviet scientists had genetically manipulated Bacillus anthracis with the potential to evade vaccine prophylaxis and antibiotic therapeutics, there was no genomic evidence of this from the Sverdlovsk production strain genome. The whole-genome SNP genotype of the Sverdlovsk strain was used to precisely identify it and its close relatives in the context of an extensive global B. anthracis strain collection. This genomic identity can now be used for forensic tracking of this weapons material on a global scale and for future anthrax investigations.

Rapid Drug Susceptibility Testing of Drug-Resistant Mycobacterium tuberculosis Isolates Directly from Clinical Samples by Use of Amplicon Sequencing: a Proof-of-Concept Study.

Increasingly complex drug-resistant tuberculosis (DR-TB) is a major global health concern and one of the primary reasons why TB is now the leading infectious cause of death worldwide. Rapid characterization of a DR-TB patient's complete drug resistance profile would facilitate individualized treatment in place of empirical treatment, improve treatment outcomes, prevent amplification of resistance, and reduce the transmission of DR-TB. The use of targeted next-generation sequencing (NGS) to obtain drug resistance profiles directly from patient sputum samples has the potential to enable comprehensive evidence-based treatment plans to be implemented quickly, rather than in weeks to months, which is currently needed for phenotypic drug susceptibility testing (DST) results. In this pilot study, we evaluated the performance of amplicon sequencing of Mycobacterium tuberculosis DNA from patient sputum samples using a tabletop NGS technology and automated data analysis to provide a rapid DST solution (the Next Gen-RDST assay). One hundred sixty-six out of 176 (94.3%) sputum samples from the Republic of Moldova yielded complete Next Gen-RDST assay profiles for 7 drugs of interest. We found a high level of concordance of our Next Gen-RDST assay results with phenotypic DST (97.0%) and pyrosequencing (97.8%) results from the same clinical samples. Our Next Gen-RDST assay was also able to estimate the proportion of resistant-to-wild-type alleles down to mixtures of ≤1%, which demonstrates the ability to detect very low levels of resistant variants not detected by pyrosequencing and possibly below the threshold for phenotypic growth methods. The assay as described here could be used as a clinical or surveillance tool.

NASP: an accurate, rapid method for the identification of SNPs in WGS datasets that supports flexible input and output formats.

Whole-genome sequencing (WGS) of bacterial isolates has become standard practice in many laboratories. Applications for WGS analysis include phylogeography and molecular epidemiology, using single nucleotide polymorphisms (SNPs) as the unit of evolution. NASP was developed as a reproducible method that scales well with the hundreds to thousands of WGS data typically used in comparative genomics applications. In this study, we demonstrate how NASP compares with other tools in the analysis of two real bacterial genomics datasets and one simulated dataset. Our results demonstrate that NASP produces similar, and often better, results in comparison with other pipelines, but is much more flexible in terms of data input types, job management systems, diversity of supported tools and output formats. We also demonstrate differences in results based on the choice of the reference genome and choice of inferring phylogenies from concatenated SNPs or alignments including monomorphic positions. NASP represents a source-available, version-controlled, unit-tested method and can be obtained from tgennorth.github.io/NASP.

Using Whole Genome Analysis to Examine Recombination across Diverse Sequence Types of Staphylococcus aureus.

Staphylococcus aureus is an important clinical pathogen worldwide and understanding this organism's phylogeny and, in particular, the role of recombination, is important both to understand the overall spread of virulent lineages and to characterize outbreaks. To further elucidate the phylogeny of S. aureus, 35 diverse strains were sequenced using whole genome sequencing. In addition, 29 publicly available whole genome sequences were included to create a single nucleotide polymorphism (SNP)-based phylogenetic tree encompassing 11 distinct lineages. All strains of a particular sequence type fell into the same clade with clear groupings of the major clonal complexes of CC8, CC5, CC30, CC45 and CC1. Using a novel analysis method, we plotted the homoplasy density and SNP density across the whole genome and found evidence of recombination throughout the entire chromosome, but when we examined individual clonal lineages we found very little recombination. However, when we analyzed three branches of multiple lineages, we saw intermediate and differing levels of recombination between them. These data demonstrate that in S. aureus, recombination occurs across major lineages that subsequently expand in a clonal manner. Estimated mutation rates for the CC8 and CC5 lineages were different from each other. While the CC8 lineage rate was similar to previous studies, the CC5 lineage was 100-fold greater. Fifty known virulence genes were screened in all genomes in silico to determine their distribution across major clades. Thirty-three genes were present variably across clades, most of which were not constrained by ancestry, indicating horizontal gene transfer or gene loss.

Comparative integrated omics: identification of key functionalities in microbial community-wide metabolic networks.

BACKGROUND: Mixed microbial communities underpin important biotechnological processes such as biological wastewater treatment (BWWT). A detailed knowledge of community structure and function relationships is essential for ultimately driving these systems towards desired outcomes, e.g., the enrichment in organisms capable of accumulating valuable resources during BWWT. METHODS: A comparative integrated omic analysis including metagenomics, metatranscriptomics and metaproteomics was carried out to elucidate functional differences between seasonally distinct oleaginous mixed microbial communities (OMMCs) sampled from an anoxic BWWT tank. A computational framework for the reconstruction of community-wide metabolic networks from multi-omic data was developed. These provide an overview of the functional capabilities by incorporating gene copy, transcript and protein abundances. To identify functional genes, which have a disproportionately important role in community function, we define a high relative gene expression and a high betweenness centrality relative to node degree as gene-centric and network topological features, respectively. RESULTS: Genes exhibiting high expression relative to gene copy abundance include genes involved in glycerolipid metabolism, particularly triacylglycerol lipase, encoded by known lipid accumulating populations, e.g., CandidatusMicrothrix parvicella. Genes with a high relative gene expression and topologically important positions in the network include genes involved in nitrogen metabolism and fatty acid biosynthesis, encoded by Nitrosomonas spp. and Rhodococcus spp. Such genes may be regarded as 'keystone genes' as they are likely to be encoded by keystone species. CONCLUSION: The linking of key functionalities to community members through integrated omics opens up exciting possibilities for devising prediction and control strategies for microbial communities in the future.

Community-integrated omics links dominance of a microbial generalist to fine-tuned resource usage.

Microbial communities are complex and dynamic systems that are primarily structured according to their members' ecological niches. To investigate how niche breadth (generalist versus specialist lifestyle strategies) relates to ecological success, we develop and apply an integrative workflow for the multi-omic analysis of oleaginous mixed microbial communities from a biological wastewater treatment plant. Time- and space-resolved coupled metabolomic and taxonomic analyses demonstrate that the community-wide lipid accumulation phenotype is associated with the dominance of the generalist bacterium Candidatus Microthrix spp. By integrating population-level genomic reconstructions (reflecting fundamental niches) with transcriptomic and proteomic data (realised niches), we identify finely tuned gene expression governing resource usage by Candidatus Microthrix parvicella over time. Moreover, our results indicate that the fluctuating environmental conditions constrain the accumulation of genetic variation in Candidatus Microthrix parvicella likely due to fitness trade-offs. Based on our observations, niche breadth has to be considered as an important factor for understanding the evolutionary processes governing (microbial) population sizes and structures in situ.

Cryptococcus gattii in North American Pacific Northwest: whole-population genome analysis provides insights into species evolution and dispersal.

The emergence of distinct populations of Cryptococcus gattii in the temperate North American Pacific Northwest (PNW) was surprising, as this species was previously thought to be confined to tropical and semitropical regions. Beyond a new habitat niche, the dominant emergent population displayed increased virulence and caused primary pulmonary disease, as opposed to the predominantly neurologic disease seen previously elsewhere. Whole-genome sequencing was performed on 118 C. gattii isolates, including the PNW subtypes and the global diversity of molecular type VGII, to better ascertain the natural source and genomic adaptations leading to the emergence of infection in the PNW. Overall, the VGII population was highly diverse, demonstrating large numbers of mutational and recombinational events; however, the three dominant subtypes from the PNW were of low diversity and were completely clonal. Although strains of VGII were found on at least five continents, all genetic subpopulations were represented or were most closely related to strains from South America. The phylogenetic data are consistent with multiple dispersal events from South America to North America and elsewhere. Numerous gene content differences were identified between the emergent clones and other VGII lineages, including genes potentially related to habitat adaptation, virulence, and pathology. Evidence was also found for possible gene introgression from Cryptococcus neoformans var. grubii that is rarely seen in global C. gattii but that was present in all PNW populations. These findings provide greater understanding of C. gattii evolution in North America and support extensive evolution in, and dispersal from, South America. Importance: Cryptococcus gattii emerged in the temperate North American Pacific Northwest (PNW) in the late 1990s. Beyond a new environmental niche, these emergent populations displayed increased virulence and resulted in a different pattern of clinical disease. In particular, severe pulmonary infections predominated in contrast to presentation with neurologic disease as seen previously elsewhere. We employed population-level whole-genome sequencing and analysis to explore the genetic relationships and gene content of the PNW C. gattii populations. We provide evidence that the PNW strains originated from South America and identified numerous genes potentially related to habitat adaptation, virulence expression, and clinical presentation. Characterization of these genetic features may lead to improved diagnostics and therapies for such fungal infections. The data indicate that there were multiple recent introductions of C. gattii into the PNW. Public health vigilance is warranted for emergence in regions where C. gattii is not thought to be endemic.