Discovery of antimicrobial peptides in the global microbiome with machine learning.

Novel antibiotics are urgently needed to combat the antibiotic-resistance crisis. We present a machine-learning-based approach to predict antimicrobial peptides (AMPs) within the global microbiome and leverage a vast dataset of 63,410 metagenomes and 87,920 prokaryotic genomes from environmental and host-associated habitats to create the AMPSphere, a comprehensive catalog comprising 863,498 non-redundant peptides, few of which match existing databases. AMPSphere provides insights into the evolutionary origins of peptides, including by duplication or gene truncation of longer sequences, and we observed that AMP production varies by habitat. To validate our predictions, we synthesized and tested 100 AMPs against clinically relevant drug-resistant pathogens and human gut commensals both in vitro and in vivo. A total of 79 peptides were active, with 63 targeting pathogens. These active AMPs exhibited antibacterial activity by disrupting bacterial membranes. In conclusion, our approach identified nearly one million prokaryotic AMP sequences, an open-access resource for antibiotic discovery.

Translocated microbiome composition determines immunological outcome in treated HIV infection.

The impact of the microbiome on HIV disease is widely acknowledged although the mechanisms downstream of fluctuations in microbial composition remain speculative. We detected rapid, dynamic changes in translocated microbial constituents during two years after cART initiation. An unbiased systems biology approach revealed two distinct pathways driven by changes in the abundance ratio of Serratia to other bacterial genera. Increased CD4 T cell numbers over the first year were associated with high Serratia abundance, pro-inflammatory innate cytokines, and metabolites that drive Th17 gene expression signatures and restoration of mucosal integrity. Subsequently, decreased Serratia abundance and downregulation of innate cytokines allowed re-establishment of systemic T cell homeostasis promoting restoration of Th1 and Th2 gene expression signatures. Analyses of three other geographically distinct cohorts of treated HIV infection established a more generalized principle that changes in diversity and composition of translocated microbial species influence systemic inflammation and consequently CD4 T cell recovery.

The Human Gut Microbiome: From Association to Modulation.

Our understanding of the human gut microbiome continues to evolve at a rapid pace, but practical application of thisknowledge is still in its infancy. This review discusses the type of studies that will be essential for translating microbiome research into targeted modulations with dedicated benefits for the human host.

Ambient oxygen levels regulate intestinal dysbiosis and GVHD severity after allogeneic stem cell transplantation.

The severity of T cell-mediated gastrointestinal (GI) diseases such as graft-versus-host disease (GVHD) and inflammatory bowel diseases correlates with a decrease in the diversity of the host gut microbiome composition characterized by loss of obligate anaerobic commensals. The mechanisms underpinning these changes in the microbial structure remain unknown. Here, we show in multiple specific pathogen-free (SPF), gnotobiotic, and germ-free murine models of GI GVHD that the initiation of the intestinal damage by the pathogenic T cells altered ambient oxygen levels in the GI tract and caused dysbiosis. The change in oxygen levels contributed to the severity of intestinal pathology in a host intestinal HIF-1α- and a microbiome-dependent manner. Regulation of intestinal ambient oxygen levels with oral iron chelation mitigated dysbiosis and reduced the severity of the GI GVHD. Thus, targeting ambient intestinal oxygen levels may represent a novel, non-immunosuppressive strategy to mitigate T cell-driven intestinal diseases.

Paternal microbiome perturbations impact offspring fitness.

The gut microbiota operates at the interface of host-environment interactions to influence human homoeostasis and metabolic networks1-4. Environmental factors that unbalance gut microbial ecosystems can therefore shape physiological and disease-associated responses across somatic tissues5-9. However, the systemic impact of the gut microbiome on the germline-and consequently on the F1 offspring it gives rise to-is unexplored10. Here we show that the gut microbiota act as a key interface between paternal preconception environment and intergenerational health in mice. Perturbations to the gut microbiota of prospective fathers increase the probability of their offspring presenting with low birth weight, severe growth restriction and premature mortality. Transmission of disease risk occurs via the germline and is provoked by pervasive gut microbiome perturbations, including non-absorbable antibiotics or osmotic laxatives, but is rescued by restoring the paternal microbiota before conception. This effect is linked with a dynamic response to induced dysbiosis in the male reproductive system, including impaired leptin signalling, altered testicular metabolite profiles and remapped small RNA payloads in sperm. As a result, dysbiotic fathers trigger an elevated risk of in utero placental insufficiency, revealing a placental origin of mammalian intergenerational effects. Our study defines a regulatory 'gut-germline axis' in males, which is sensitive to environmental exposures and programmes offspring fitness through impacting placenta function.

Functional and evolutionary significance of unknown genes from uncultivated taxa.

Many of the Earth's microbes remain uncultured and understudied, limiting our understanding of the functional and evolutionary aspects of their genetic material, which remain largely overlooked in most metagenomic studies1. Here we analysed 149,842 environmental genomes from multiple habitats2-6 and compiled a curated catalogue of 404,085 functionally and evolutionarily significant novel (FESNov) gene families exclusive to uncultivated prokaryotic taxa. All FESNov families span multiple species, exhibit strong signals of purifying selection and qualify as new orthologous groups, thus nearly tripling the number of bacterial and archaeal gene families described to date. The FESNov catalogue is enriched in clade-specific traits, including 1,034 novel families that can distinguish entire uncultivated phyla, classes and orders, probably representing synapomorphies that facilitated their evolutionary divergence. Using genomic context analysis and structural alignments we predicted functional associations for 32.4% of FESNov families, including 4,349 high-confidence associations with important biological processes. These predictions provide a valuable hypothesis-driven framework that we used for experimental validatation of a new gene family involved in cell motility and a novel set of antimicrobial peptides. We also demonstrate that the relative abundance profiles of novel families can discriminate between environments and clinical conditions, leading to the discovery of potentially new biomarkers associated with colorectal cancer. We expect this work to enhance future metagenomics studies and expand our knowledge of the genetic repertory of uncultivated organisms.

Gut microbiome-derived metabolites modulate intestinal epithelial cell damage and mitigate graft-versus-host disease.

The effect of alterations in intestinal microbiota on microbial metabolites and on disease processes such as graft-versus-host disease (GVHD) is not known. Here we carried out an unbiased analysis to identify previously unidentified alterations in gastrointestinal microbiota-derived short-chain fatty acids (SCFAs) after allogeneic bone marrow transplant (allo-BMT). Alterations in the amount of only one SCFA, butyrate, were observed only in the intestinal tissue. The reduced butyrate in CD326(+) intestinal epithelial cells (IECs) after allo-BMT resulted in decreased histone acetylation, which was restored after local administration of exogenous butyrate. Butyrate restoration improved IEC junctional integrity, decreased apoptosis and mitigated GVHD. Furthermore, alteration of the indigenous microbiota with 17 rationally selected strains of high butyrate-producing Clostridia also decreased GVHD. These data demonstrate a heretofore unrecognized role of microbial metabolites and suggest that local and specific alteration of microbial metabolites has direct salutary effects on GVHD target tissues and can mitigate disease severity.

Targeting placental growth factor/neuropilin 1 pathway inhibits growth and spread of medulloblastoma.

Medulloblastoma is the most common pediatric malignant brain tumor. Although current therapies improve survival, these regimens are highly toxic and are associated with significant morbidity. Here, we report that placental growth factor (PlGF) is expressed in the majority of medulloblastomas, independent of their subtype. Moreover, high expression of PlGF receptor neuropilin 1 (Nrp1) correlates with poor overall survival in patients. We demonstrate that PlGF and Nrp1 are required for the growth and spread of medulloblastoma: PlGF/Nrp1 blockade results in direct antitumor effects in vivo, resulting in medulloblastoma regression, decreased metastasis, and increased mouse survival. We reveal that PlGF is produced in the cerebellar stroma via tumor-derived Sonic hedgehog (Shh) and show that PlGF acts through Nrp1-and not vascular endothelial growth factor receptor 1-to promote tumor cell survival. This critical tumor-stroma interaction-mediated by Shh, PlGF, and Nrp1 across medulloblastoma subtypes-supports the development of therapies targeting PlGF/Nrp1 pathway.

Towards the biogeography of prokaryotic genes.

Microbial genes encode the majority of the functional repertoire of life on earth. However, despite increasing efforts in metagenomic sequencing of various habitats1-3, little is known about the distribution of genes across the global biosphere, with implications for human and planetary health. Here we constructed a non-redundant gene catalogue of 303 million species-level genes (clustered at 95% nucleotide identity) from 13,174 publicly available metagenomes across 14 major habitats and use it to show that most genes are specific to a single habitat. The small fraction of genes found in multiple habitats is enriched in antibiotic-resistance genes and markers for mobile genetic elements. By further clustering these species-level genes into 32 million protein families, we observed that a small fraction of these families contain the majority of the genes (0.6% of families account for 50% of the genes). The majority of species-level genes and protein families are rare. Furthermore, species-level genes, and in particular the rare ones, show low rates of positive (adaptive) selection, supporting a model in which most genetic variability observed within each protein family is neutral or nearly neutral.

Combinatorial, additive and dose-dependent drug-microbiome associations.

During the transition from a healthy state to cardiometabolic disease, patients become heavily medicated, which leads to an increasingly aberrant gut microbiome and serum metabolome, and complicates biomarker discovery1-5. Here, through integrated multi-omics analyses of 2,173 European residents from the MetaCardis cohort, we show that the explanatory power of drugs for the variability in both host and gut microbiome features exceeds that of disease. We quantify inferred effects of single medications, their combinations as well as additive effects, and show that the latter shift the metabolome and microbiome towards a healthier state, exemplified in synergistic reduction in serum atherogenic lipoproteins by statins combined with aspirin, or enrichment of intestinal Roseburia by diuretic agents combined with beta-blockers. Several antibiotics exhibit a quantitative relationship between the number of courses prescribed and progression towards a microbiome state that is associated with the severity of cardiometabolic disease. We also report a relationship between cardiometabolic drug dosage, improvement in clinical markers and microbiome composition, supporting direct drug effects. Taken together, our computational framework and resulting resources enable the disentanglement of the effects of drugs and disease on host and microbiome features in multimedicated individuals. Furthermore, the robust signatures identified using our framework provide new hypotheses for drug-host-microbiome interactions in cardiometabolic disease.

A male-killing Wolbachia endosymbiont is concealed by another endosymbiont and a nuclear suppressor.

Bacteria that live inside the cells of insect hosts (endosymbionts) can alter the reproduction of their hosts, including the killing of male offspring (male killing, MK). MK has only been described in a few insects, but this may reflect challenges in detecting MK rather than its rarity. Here, we identify MK Wolbachia at a low frequency (around 4%) in natural populations of Drosophila pseudotakahashii. MK Wolbachia had a stable density and maternal transmission during laboratory culture, but the MK phenotype which manifested mainly at the larval stage was lost rapidly. MK Wolbachia occurred alongside a second Wolbachia strain expressing a different reproductive manipulation, cytoplasmic incompatibility (CI). A genomic analysis highlighted Wolbachia regions diverged between the 2 strains involving 17 genes, and homologs of the wmk and cif genes implicated in MK and CI were identified in the Wolbachia assembly. Doubly infected males induced CI with uninfected females but not females singly infected with CI-causing Wolbachia. A rapidly spreading dominant nuclear suppressor genetic element affecting MK was identified through backcrossing and subsequent analysis with ddRAD SNPs of the D. pseudotakahashii genome. These findings highlight the complexity of nuclear and microbial components affecting MK endosymbiont detection and dynamics in populations and the challenges of making connections between endosymbionts and the host phenotypes affected by them.

SPIRE: a Searchable, Planetary-scale mIcrobiome REsource.

Meta'omic data on microbial diversity and function accrue exponentially in public repositories, but derived information is often siloed according to data type, study or sampled microbial environment. Here we present SPIRE, a Searchable Planetary-scale mIcrobiome REsource that integrates various consistently processed metagenome-derived microbial data modalities across habitats, geography and phylogeny. SPIRE encompasses 99 146 metagenomic samples from 739 studies covering a wide array of microbial environments and augmented with manually-curated contextual data. Across a total metagenomic assembly of 16 Tbp, SPIRE comprises 35 billion predicted protein sequences and 1.16 million newly constructed metagenome-assembled genomes (MAGs) of medium or high quality. Beyond mapping to the high-quality genome reference provided by proGenomes3 (http://progenomes.embl.de), these novel MAGs form 92 134 novel species-level clusters, the majority of which are unclassified at species level using current tools. SPIRE enables taxonomic profiling of these species clusters via an updated, custom mOTUs database (https://motu-tool.org/) and includes several layers of functional annotation, as well as crosslinks to several (micro-)biological databases. The resource is accessible, searchable and browsable via http://spire.embl.de.

Quantitative analysis of sterol-modulated monomer-dimer equilibrium of the ß1-adrenergic receptor by DEER spectroscopy.

G protein-coupled receptors (GPCR) activate numerous intracellular signaling pathways. The oligomerization properties of GPCRs, and hence their cellular functions, may be modulated by various components within the cell membrane (such as the presence of cholesterol). Modulation may occur directly via specific interaction with the GPCR or indirectly by affecting the physical properties of the membrane. Here, we use pulsed Q-band double electron-electron resonance (DEER) spectroscopy to probe distances between R1 nitroxide spin labels attached to Cys163 and Cys344 of the ß1-adrenergic receptor (ß1AR) in n-dodecyl-ß-D-maltoside micelles upon titration with two soluble cholesterol analogs, cholesteryl hemisuccinate (CHS) and sodium cholate. The former, like cholesterol, inserts itself into the lipid membrane, parallel to the phospholipid chains; the latter is aligned parallel to the surface of membranes. Global quantitative analysis of DEER echo curves upon titration of spin-labeled ß1AR with CHS and sodium cholate reveal the following: CHS binds specifically to the ß1AR monomer at a site close to the Cys163-R1 spin label with an equilibrium dissociation constant [Formula: see text] ~1.4 ± 0.4 mM. While no direct binding of sodium cholate to the ß1AR receptor was observed by DEER, sodium cholate induces specific ß1AR dimerization ([Formula: see text] ~35 ± 6 mM and a Hill coefficient n ~ 2.5 ± 0.4) with intersubunit contacts between transmembrane helices 1 and 2 and helix 8. Analysis of the DEER data obtained upon the addition of CHS to the ß1AR dimer in the presence of excess cholate results in dimer dissociation with species occupancies as predicted from the individual KD values.

Further pharmacological and genetic evidence for the efficacy of PlGF inhibition in cancer and eye disease.

Our findings that PlGF is a cancer target and anti-PlGF is useful for anticancer treatment have been challenged by Bais et al. Here we take advantage of carcinogen-induced and transgenic tumor models as well as ocular neovascularization to report further evidence in support of our original findings of PlGF as a promising target for anticancer therapies. We present evidence for the efficacy of additional anti-PlGF antibodies and their ability to phenocopy genetic deficiency or silencing of PlGF in cancer and ocular disease but also show that not all anti-PlGF antibodies are effective. We also provide additional evidence for the specificity of our anti-PlGF antibody and experiments to suggest that anti-PlGF treatment will not be effective for all tumors and why. Further, we show that PlGF blockage inhibits vessel abnormalization rather than density in certain tumors while enhancing VEGF-targeted inhibition in ocular disease. Our findings warrant further testing of anti-PlGF therapies.

proGenomes3: approaching one million accurately and consistently annotated high-quality prokaryotic genomes.

The interpretation of genomic, transcriptomic and other microbial 'omics data is highly dependent on the availability of well-annotated genomes. As the number of publicly available microbial genomes continues to increase exponentially, the need for quality control and consistent annotation is becoming critical. We present proGenomes3, a database of 907 388 high-quality genomes containing 4 billion genes that passed stringent criteria and have been consistently annotated using multiple functional and taxonomic databases including mobile genetic elements and biosynthetic gene clusters. proGenomes3 encompasses 41 171 species-level clusters, defined based on universal single copy marker genes, for which pan-genomes and contextual habitat annotations are provided. The database is available at http://progenomes.embl.de/.

Time-resolved DEER EPR and solid-state NMR afford kinetic and structural elucidation of substrate binding to Ca2+-ligated calmodulin.

Recent advances in rapid mixing and freeze quenching have opened the path for time-resolved electron paramagnetic resonance (EPR)-based double electron-electron resonance (DEER) and solid-state NMR of protein-substrate interactions. DEER, in conjunction with phase memory time filtering to quantitatively extract species populations, permits monitoring time-dependent probability distance distributions between pairs of spin labels, while solid-state NMR provides quantitative residue-specific information on the appearance of structural order and the development of intermolecular contacts between substrate and protein. Here, we demonstrate the power of these combined approaches to unravel the kinetic and structural pathways in the binding of the intrinsically disordered peptide substrate (M13) derived from myosin light-chain kinase to the universal eukaryotic calcium regulator, calmodulin. Global kinetic analysis of the data reveals coupled folding and binding of the peptide associated with large spatial rearrangements of the two domains of calmodulin. The initial binding events involve a bifurcating pathway in which the M13 peptide associates via either its N- or C-terminal regions with the C- or N-terminal domains, respectively, of calmodulin/4Ca2+ to yield two extended "encounter" complexes, states A and A*, without conformational ordering of M13. State A is immediately converted to the final compact complex, state C, on a timescale τ ≤ 600 µs. State A*, however, only reaches the final complex via a collapsed intermediate B (τ ∼ 1.5 to 2.5 ms), in which the peptide is only partially ordered and not all intermolecular contacts are formed. State B then undergoes a relatively slow (τ ∼ 7 to 18 ms) conformational rearrangement to state C.

Deconvoluting Monomer- and Dimer-Specific Distance Distributions between Spin Labels in a Monomer/Dimer Mixture Using T1-Edited DEER EPR Spectroscopy.

Double electron-electron resonance (DEER) EPR is a powerful tool in structural biology, providing distances between pairs of spin labels. When the sample consists of a mixture of oligomeric species (e.g., monomer and dimer), the question arises as to how to assign the peaks in the DEER-derived probability distance distribution to the individual species. Here, we propose incorporating an EPR longitudinal electron relaxation (T1) inversion recovery experiment within a DEER pulse sequence to resolve this problem. The apparent T1 between dipolar coupled electron spins measured from the inversion recovery time (τinv) dependence of the peak intensities in the T1-edited DEER-derived probability P(r) distance distribution will be affected by the number of nitroxide labels attached to the biomolecule of interest, for example, two for a monomer and four for a dimer. We show that global fitting of all the T1-edited DEER echo curves, recorded over a range of τinv values, permits the deconvolution of distances between spin labels originating from monomeric (longer T1) and dimeric (shorter T1) species. This is especially useful when the trapping of spin labels in different conformational states during freezing gives rise to complex P(r) distance distributions. The utility of this approach is demonstrated for two systems, the ß1 adrenergic receptor and a construct of the huntingtin exon-1 protein fused to the immunoglobulin domain of protein G, both of which exist in a monomer-dimer equilibrium.

Stratifying esophago-gastric cancer treatment using a patient-derived organoid-based threshold.

BACKGROUND AND AIMS: This study sought to determine the value of patient-derived organoids (PDOs) from esophago-gastric adenocarcinoma (EGC) for response prediction to neoadjuvant chemotherapy (neoCTx). METHODS: Endoscopic biopsies of patients with locally advanced EGC (n = 120) were taken into culture and PDOs expanded. PDOs' response towards the single substances of the FLOT regimen and the combination treatment were correlated to patients' pathological response using tumor regression grading. A classifier based on FLOT response of PDOs was established in an exploratory cohort (n = 13) and subsequently confirmed in an independent validation cohort (n = 13). RESULTS: EGC PDOs reflected patients' diverse responses to single chemotherapeutics and the combination regimen FLOT. In the exploratory cohort, PDOs response to single 5-FU and FLOT combination treatment correlated with the patients' pathological response (5-FU: Kendall's τ = 0.411, P = 0.001; FLOT: Kendall's τ = 0.694, P = 2.541e-08). For FLOT testing, a high diagnostic precision in receiver operating characteristic (ROC) analysis was reached with an AUCROC of 0.994 (CI 0.980 to 1.000). The discriminative ability of PDO-based FLOT testing allowed the definition of a threshold, which classified in an independent validation cohort FLOT responders from non-responders with high sensitivity (90%), specificity (100%) and accuracy (92%). CONCLUSION: In vitro drug testing of EGC PDOs has a high predictive accuracy in classifying patients' histological response to neoadjuvant FLOT treatment. Taking into account the high rate of successful PDO expansion from biopsies, the definition of a threshold that allows treatment stratification paves the way for an interventional trial exploring PDO-guided treatment of EGC patients.

A decade of stability for wMel Wolbachia in natural Aedes aegypti populations.

Mosquitoes carrying Wolbachia endosymbionts are being released in many countries for arbovirus control. The wMel strain of Wolbachia blocks Aedes-borne virus transmission and can spread throughout mosquito populations by inducing cytoplasmic incompatibility. Aedes aegypti mosquitoes carrying wMel were first released into the field in Cairns, Australia, over a decade ago, and with wider releases have resulted in the near elimination of local dengue transmission. The long-term stability of Wolbachia effects is critical for ongoing disease suppression, requiring tracking of phenotypic and genomic changes in Wolbachia infections following releases. We used a combination of field surveys, phenotypic assessments, and Wolbachia genome sequencing to show that wMel has remained stable in its effects for up to a decade in Australian Ae. aegypti populations. Phenotypic comparisons of wMel-infected and uninfected mosquitoes from near-field and long-term laboratory populations suggest limited changes in the effects of wMel on mosquito fitness. Treating mosquitoes with antibiotics used to cure the wMel infection had limited effects on fitness in the next generation, supporting the use of tetracycline for generating uninfected mosquitoes without off-target effects. wMel has a temporally stable within-host density and continues to induce complete cytoplasmic incompatibility. A comparison of wMel genomes from pre-release (2010) and nine years post-release (2020) populations show few genomic differences and little divergence between release locations, consistent with the lack of phenotypic changes. These results indicate that releases of Wolbachia-infected mosquitoes for population replacement are likely to be effective for many years, but ongoing monitoring remains important to track potential evolutionary changes.

Consistency across multi-omics layers in a drug-perturbed gut microbial community.

Multi-omics analyses are used in microbiome studies to understand molecular changes in microbial communities exposed to different conditions. However, it is not always clear how much each omics data type contributes to our understanding and whether they are concordant with each other. Here, we map the molecular response of a synthetic community of 32 human gut bacteria to three non-antibiotic drugs by using five omics layers (16S rRNA gene profiling, metagenomics, metatranscriptomics, metaproteomics and metabolomics). We find that all the omics methods with species resolution are highly consistent in estimating relative species abundances. Furthermore, different omics methods complement each other for capturing functional changes. For example, while nearly all the omics data types captured that the antipsychotic drug chlorpromazine selectively inhibits Bacteroidota representatives in the community, the metatranscriptome and metaproteome suggested that the drug induces stress responses related to protein quality control. Metabolomics revealed a decrease in oligosaccharide uptake, likely caused by Bacteroidota depletion. Our study highlights how multi-omics datasets can be utilized to reveal complex molecular responses to external perturbations in microbial communities.