Meeting Report of the Second Symposium of the Belgian Society for Viruses of Microbes and Launch of the Phage Valley.

The second symposium of the Belgian Society for Viruses of Microbes (BSVoM) took place on 8 September 2023 at the University of Liège with 141 participants from 10 countries. The meeting program covered three thematic sessions opened by international keynote speakers: two sessions were devoted to "Fundamental research in phage ecology and biology" and the third one to the "Present and future applications of phages". During this one day symposium, four invited keynote lectures, nine selected talks and eight student pitches were given along with thirty presented posters. The president of the Belgian Society for Viruses of Microbes, Prof. Yves Briers, took advantage of this symposium to launch the Phage Valley concept that will put the spotlight on the exceptionally high density of researchers investigating viruses of microbes as well as the successful triple helix approach between academia, industry and government in Belgium.

Foundation of the Belgian Society for Viruses of Microbes and Meeting Report of Its Inaugural Symposium.

The Belgian Society for Viruses of Microbes (BSVoM) was founded on 9 June 2022 to capture and enhance the collaborative spirit among the expanding community of microbial virus researchers in Belgium. The sixteen founders are affiliated to fourteen different research entities across academia, industry and government. Its inaugural symposium was held on 23 September 2022 in the Thermotechnical Institute at KU Leuven. The meeting program covered three thematic sessions launched by international keynote speakers: (1) virus-host interactions, (2) viral ecology, evolution and diversity and (3) present and future applications. During the one-day symposium, four invited keynote lectures, ten selected talks and eight student pitches were given along with 41 presented posters. The meeting hosted 155 participants from twelve countries.

Structure of SpoT reveals evolutionary tuning of catalysis via conformational constraint.

Stringent factors orchestrate bacterial cell reprogramming through increasing the level of the alarmones (p)ppGpp. In Beta- and Gammaproteobacteria, SpoT hydrolyzes (p)ppGpp to counteract the synthetase activity of RelA. However, structural information about how SpoT controls the levels of (p)ppGpp is missing. Here we present the crystal structure of the hydrolase-only SpoT from Acinetobacter baumannii and uncover the mechanism of intramolecular regulation of 'long'-stringent factors. In contrast to ribosome-associated Rel/RelA that adopt an elongated structure, SpoT assumes a compact τ-shaped structure in which the regulatory domains wrap around a Core subdomain that controls the conformational state of the enzyme. The Core is key to the specialization of long RelA-SpoT homologs toward either synthesis or hydrolysis: the short and structured Core of SpoT stabilizes the τ-state priming the hydrolase domain for (p)ppGpp hydrolysis, whereas the longer, more dynamic Core domain of RelA destabilizes the τ-state priming the monofunctional RelA for efficient (p)ppGpp synthesis.

Rhizomania: Hide and Seek of Polymyxa betae and the Beet Necrotic Yellow Vein Virus with Beta vulgaris.

The molecular interactions between Polymyxa betae, the protist vector of sugar beet viruses, beet necrotic yellow vein virus (BNYVV), the causal agent of rhizomania, and Beta vulgaris have not been extensively studied. Here, the transmission of BNYVV to sugar beet by P. betae zoospores was optimized using genetically characterized organisms. Molecular interactions of aviruliferous and viruliferous protist infection on sugar beet were highlighted by transcriptomic analysis. P. betae alone induced limited gene expression changes in sugar beet, as a biotrophic asymptomatic parasite. Most differentially expressed plant genes were down-regulated and included resistance gene analogs and cell wall peroxidases. Several enzymes involved in stress regulation, such as the glutathione-S-transferases, were significantly induced. With BNYVV, the first stages of the P. betae life cycle on sugar beet were accelerated with a faster increase of relative protist DNA level and an earlier appearance of sporangia and sporosori in plants roots. A clear activation of plant defenses and the modulation of genes involved in plant cell wall metabolism were observed. The P. betae transcriptome in the presence of BNYVV revealed induction of genes possibly involved in the switch to the survival stage. The interactions were different depending on the presence or absence of the virus. P. betae alone alleviates plant defense response, playing hide-and-seek with sugar beet and allowing for their mutual development. Conversely, BNYVV manipulates plant defense and promotes the rapid invasion of plant roots by P. betae. This accelerated colonization is accompanied by the development of thick-walled resting spores, supporting the virus survival. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Metagenomics approach for Polymyxa betae genome assembly enables comparative analysis towards deciphering the intracellular parasitic lifestyle of the plasmodiophorids.

Genomic knowledge of the tree of life is biased to specific groups of organisms. For example, only six full genomes are currently available in the rhizaria clade. Here, we have applied metagenomic techniques enabling the assembly of the genome of Polymyxa betae (Rhizaria, Plasmodiophorida) RES F41 isolate from unpurified zoospore holobiont and comparison with the A26-41 isolate. Furthermore, the first P. betae mitochondrial genome was assembled. The two P. betae nuclear genomes were highly similar, each with just ~10.2 k predicted protein coding genes, ~3% of which were unique to each isolate. Extending genomic comparisons revealed a greater overlap with Spongospora subterranea than with Plasmodiophora brassicae, including orthologs of the mammalian cation channel sperm-associated proteins, raising some intriguing questions about zoospore physiology. This work validates our metagenomics pipeline for eukaryote genome assembly from unpurified samples and enriches plasmodiophorid genomics; providing the first full annotation of the P. betae genome.

Disentangling environmental effects in microbial association networks.

BACKGROUND: Ecological interactions among microorganisms are fundamental for ecosystem function, yet they are mostly unknown or poorly understood. High-throughput-omics can indicate microbial interactions through associations across time and space, which can be represented as association networks. Associations could result from either ecological interactions between microorganisms, or from environmental selection, where the association is environmentally driven. Therefore, before downstream analysis and interpretation, we need to distinguish the nature of the association, particularly if it is due to environmental selection or not. RESULTS: We present EnDED (environmentally driven edge detection), an implementation of four approaches as well as their combination to predict which links between microorganisms in an association network are environmentally driven. The four approaches are sign pattern, overlap, interaction information, and data processing inequality. We tested EnDED on networks from simulated data of 50 microorganisms. The networks contained on average 50 nodes and 1087 edges, of which 60 were true interactions but 1026 false associations (i.e., environmentally driven or due to chance). Applying each method individually, we detected a moderate to high number of environmentally driven edges-87% sign pattern and overlap, 67% interaction information, and 44% data processing inequality. Combining these methods in an intersection approach resulted in retaining more interactions, both true and false (32% of environmentally driven associations). After validation with the simulated datasets, we applied EnDED on a marine microbial network inferred from 10 years of monthly observations of microbial-plankton abundance. The intersection combination predicted that 8.3% of the associations were environmentally driven, while individual methods predicted 24.8% (data processing inequality), 25.7% (interaction information), and up to 84.6% (sign pattern as well as overlap). The fraction of environmentally driven edges among negative microbial associations in the real network increased rapidly with the number of environmental factors. CONCLUSIONS: To reach accurate hypotheses about ecological interactions, it is important to determine, quantify, and remove environmentally driven associations in marine microbial association networks. For that, EnDED offers up to four individual methods as well as their combination. However, especially for the intersection combination, we suggest using EnDED with other strategies to reduce the number of false associations and consequently the number of potential interaction hypotheses. Video abstract.

HOXA1 Is an Antagonist of ERα in Breast Cancer.

Breast cancer is a heterogeneous disease and the leading cause of female cancer mortality worldwide. About 70% of breast cancers express ERα. HOX proteins are master regulators of embryo development which have emerged as being important players in oncogenesis. HOXA1 is one of them. Here, we present bioinformatic analyses of genome-wide mRNA expression profiles available in large public datasets of human breast cancer samples. We reveal an extremely strong opposite correlation between HOXA1 versus ER expression and that of 2,486 genes, thereby supporting a functional antagonism between HOXA1 and ERα. We also demonstrate in vitro that HOXA1 can inhibit ERα activity. This inhibition is at least bimodal, requiring an intact HOXA1 DNA-binding homeodomain and involving the DNA-binding independent capacity of HOXA1 to activate NF-κB. We provide evidence that the HOXA1-PBX interaction known to be critical for the transcriptional activity of HOXA1 is not involved in the ERα inhibition. Finally, we reveal that HOXA1 and ERα can physically interact but that this interaction is not essential for the HOXA1-mediated inhibition of ERα. Like other HOX oncoproteins interacting with ERα, HOXA1 could be involved in endocrine therapy resistance.

Deep ocean metagenomes provide insight into the metabolic architecture of bathypelagic microbial communities.

The deep sea, the largest ocean's compartment, drives planetary-scale biogeochemical cycling. Yet, the functional exploration of its microbial communities lags far behind other environments. Here we analyze 58 metagenomes from tropical and subtropical deep oceans to generate the Malaspina Gene Database. Free-living or particle-attached lifestyles drive functional differences in bathypelagic prokaryotic communities, regardless of their biogeography. Ammonia and CO oxidation pathways are enriched in the free-living microbial communities and dissimilatory nitrate reduction to ammonium and H2 oxidation pathways in the particle-attached, while the Calvin Benson-Bassham cycle is the most prevalent inorganic carbon fixation pathway in both size fractions. Reconstruction of the Malaspina Deep Metagenome-Assembled Genomes reveals unique non-cyanobacterial diazotrophic bacteria and chemolithoautotrophic prokaryotes. The widespread potential to grow both autotrophically and heterotrophically suggests that mixotrophy is an ecologically relevant trait in the deep ocean. These results expand our understanding of the functional microbial structure and metabolic capabilities of the largest Earth aquatic ecosystem.

Toxin-Antitoxin Gene Pairs Found in Tn3 Family Transposons Appear To Be an Integral Part of the Transposition Module.

Much of the diversity of prokaryotic genomes is contributed by the tightly controlled recombination activity of transposons (Tns). The Tn3 family is arguably one of the most widespread transposon families. Members carry a large range of passenger genes incorporated into their structures. Family members undergo replicative transposition using a DDE transposase to generate a cointegrate structure which is then resolved by site-specific recombination between specific DNA sequences (res) on each of the two Tn copies in the cointegrate. These sites also carry promoters controlling expression of the recombinase and transposase. We report here that a number of Tn3 members encode a type II toxin-antitoxin (TA) system, typically composed of a stable toxin and a labile antitoxin that binds the toxin and inhibits its lethal activity. This system serves to improve plasmid maintenance in a bacterial population and, until recently, was believed to be associated with bacterial persistence. At least six different TA gene pairs are associated with various Tn3 members. Our data suggest that several independent acquisition events have occurred. In contrast to most Tn3 family passenger genes, which are generally located away from the transposition module, the TA gene pairs abut the res site upstream of the resolvase genes. Although their role when part of Tn3 family transposons is unclear, this finding suggests a potential role for the embedded TA in stabilizing the associated transposon with the possibility that TA expression is coupled to expression of transposase and resolvase during the transposition process itself.IMPORTANCE Transposable elements (TEs) are important in genetic diversification due to their recombination properties and their ability to promote horizontal gene transfer. Over the last decades, much effort has been made to understand TE transposition mechanisms and their impact on prokaryotic genomes. For example, the Tn3 family is ubiquitous in bacteria, molding their host genomes by the paste-and-copy mechanism. In addition to the transposition module, Tn3 members often carry additional passenger genes (e.g., conferring antibiotic or heavy metal resistance and virulence), and three were previously known to carry a toxin-antitoxin (TA) system often associated with plasmid maintenance; however, the role of TA systems within the Tn3 family is unknown. The genetic context of TA systems in Tn3 members suggests that they may play a regulatory role in ensuring stable invasion of these Tns during transposition.

Population-level analysis of Blastocystis subtype prevalence and variation in the human gut microbiota.

OBJECTIVE: Human gut microbiome studies are mainly bacteria- and archaea-oriented, overlooking the presence of single-cell eukaryotes such as Blastocystis, an enteric stramenopiles with worldwide distribution. Here, we surveyed the prevalence and subtype variation of Blastocystis in faecal samples collected as part of the Flemish Gut Flora Project (FGFP), a Western population cohort. We assessed potential links between Blastocystis subtypes and identified microbiota-host covariates and quantified microbiota differentiation relative to subtype abundances. DESIGN: We profiled stool samples from 616 healthy individuals from the FGFP cohort as well as 107 patients with IBD using amplicon sequencing targeting the V4 variable region of the 16S rRNA and 18S rRNA genes. We evaluated associations of Blastocystis, and their subtypes, with host parameters, diversity and composition of bacterial and archaeal communities. RESULTS: Blastocystis prevalence in the non-clinical population cohort was 30% compared with 4% among Flemish patients with IBD. Within the FGFP cohort, out of 69 previously identified gut microbiota covariates, only age was associated with Blastocystis subtype carrier status. In contrast, a strong association between microbiota community composition and Blastocystis subtypes was observed, with effect sizes larger than that of host covariates. Microbial richness and diversity were linked to both Blastocystis prevalence and subtype variation. All Blastocystis subtypes detected in this cohort were found to be less prevalent in Bacteroides enterotyped samples. Interestingly, Blastocystis subtypes 3 and 4 were inversely correlated with Akkermansia, suggesting differential associations of subtypes with host health. CONCLUSIONS: These results emphasise the role of Blastocystis as a common constituent of the healthy gut microbiota. We show its prevalence is reduced in patients with active IBD and demonstrate that subtype characterisation is essential for assessing the relationship between Blastocystis, microbiota profile and host health. These findings have direct clinical applications, especially in donor selection for faecal transplantation.

A global ocean atlas of eukaryotic genes.

While our knowledge about the roles of microbes and viruses in the ocean has increased tremendously due to recent advances in genomics and metagenomics, research on marine microbial eukaryotes and zooplankton has benefited much less from these new technologies because of their larger genomes, their enormous diversity, and largely unexplored physiologies. Here, we use a metatranscriptomics approach to capture expressed genes in open ocean Tara Oceans stations across four organismal size fractions. The individual sequence reads cluster into 116 million unigenes representing the largest reference collection of eukaryotic transcripts from any single biome. The catalog is used to unveil functions expressed by eukaryotic marine plankton, and to assess their functional biogeography. Almost half of the sequences have no similarity with known proteins, and a great number belong to new gene families with a restricted distribution in the ocean. Overall, the resource provides the foundations for exploring the roles of marine eukaryotes in ocean ecology and biogeochemistry.

Species-function relationships shape ecological properties of the human gut microbiome.

Despite recent progress, the organization and ecological properties of the intestinal microbial ecosystem remain under-investigated. Here, using a manually curated metabolic module framework for (meta-)genomic data analysis, we studied species-function relationships in gut microbial genomes and microbiomes. Half of gut-associated species were found to be generalists regarding overall substrate preference, but we observed significant genus-level metabolic diversification linked to bacterial life strategies. Within each genus, metabolic consistency varied significantly, being low in Firmicutes genera and higher in Bacteroides. Differentiation of fermentable substrate degradation potential contributed to metagenomic functional repertoire variation between individuals, with different enterotypes showing distinct saccharolytic/proteolytic/lipolytic profiles. Finally, we found that module-derived functional redundancy was reduced in the low-richness Bacteroides enterotype, potentially indicating a decreased resilience to perturbation, in line with its frequent association to dysbiosis. These results provide insights into the complex structure of gut microbiome-encoded metabolic properties and emphasize the importance of functional and ecological assessment of gut microbiome variation in clinical studies.

Population-level analysis of gut microbiome variation.

Fecal microbiome variation in the average, healthy population has remained under-investigated. Here, we analyzed two independent, extensively phenotyped cohorts: the Belgian Flemish Gut Flora Project (FGFP; discovery cohort; N = 1106) and the Dutch LifeLines-DEEP study (LLDeep; replication; N = 1135). Integration with global data sets (N combined = 3948) revealed a 14-genera core microbiota, but the 664 identified genera still underexplore total gut diversity. Sixty-nine clinical and questionnaire-based covariates were found associated to microbiota compositional variation with a 92% replication rate. Stool consistency showed the largest effect size, whereas medication explained largest total variance and interacted with other covariate-microbiota associations. Early-life events such as birth mode were not reflected in adult microbiota composition. Finally, we found that proposed disease marker genera associated to host covariates, urging inclusion of the latter in study design.

Cyanobacterial symbionts diverged in the late Cretaceous towards lineage-specific nitrogen fixation factories in single-celled phytoplankton.

The unicellular cyanobacterium UCYN-A, one of the major contributors to nitrogen fixation in the open ocean, lives in symbiosis with single-celled phytoplankton. UCYN-A includes several closely related lineages whose partner fidelity, genome-wide expression and time of evolutionary divergence remain to be resolved. Here we detect and distinguish UCYN-A1 and UCYN-A2 lineages in symbiosis with two distinct prymnesiophyte partners in the South Atlantic Ocean. Both symbiotic systems are lineage specific and differ in the number of UCYN-A cells involved. Our analyses infer a streamlined genome expression towards nitrogen fixation in both UCYN-A lineages. Comparative genomics reveal a strong purifying selection in UCYN-A1 and UCYN-A2 with a diversification process ∼91 Myr ago, in the late Cretaceous, after the low-nutrient regime period occurred during the Jurassic. These findings suggest that UCYN-A diversified in a co-evolutionary process, wherein their prymnesiophyte partners acted as a barrier driving an allopatric speciation of extant UCYN-A lineages.

Cross-biome comparison of microbial association networks.

Clinical and environmental meta-omics studies are accumulating an ever-growing amount of microbial abundance data over a wide range of ecosystems. With a sufficiently large sample number, these microbial communities can be explored by constructing and analyzing co-occurrence networks, which detect taxon associations from abundance data and can give insights into community structure. Here, we investigate how co-occurrence networks differ across biomes and which other factors influence their properties. For this, we inferred microbial association networks from 20 different 16S rDNA sequencing data sets and observed that soil microbial networks harbor proportionally fewer positive associations and are less densely interconnected than host-associated networks. After excluding sample number, sequencing depth and beta-diversity as possible drivers, we found a negative correlation between community evenness and positive edge percentage. This correlation likely results from a skewed distribution of negative interactions, which take place preferentially between less prevalent taxa. Overall, our results suggest an under-appreciated role of evenness in shaping microbial association networks.

Ocean plankton. Structure and function of the global ocean microbiome.

Microbes are dominant drivers of biogeochemical processes, yet drawing a global picture of functional diversity, microbial community structure, and their ecological determinants remains a grand challenge. We analyzed 7.2 terabases of metagenomic data from 243 Tara Oceans samples from 68 locations in epipelagic and mesopelagic waters across the globe to generate an ocean microbial reference gene catalog with >40 million nonredundant, mostly novel sequences from viruses, prokaryotes, and picoeukaryotes. Using 139 prokaryote-enriched samples, containing >35,000 species, we show vertical stratification with epipelagic community composition mostly driven by temperature rather than other environmental factors or geography. We identify ocean microbial core functionality and reveal that >73% of its abundance is shared with the human gut microbiome despite the physicochemical differences between these two ecosystems.

Ocean plankton. Eukaryotic plankton diversity in the sunlit ocean.

Marine plankton support global biological and geochemical processes. Surveys of their biodiversity have hitherto been geographically restricted and have not accounted for the full range of plankton size. We assessed eukaryotic diversity from 334 size-fractionated photic-zone plankton communities collected across tropical and temperate oceans during the circumglobal Tara Oceans expedition. We analyzed 18S ribosomal DNA sequences across the intermediate plankton-size spectrum from the smallest unicellular eukaryotes (protists, >0.8 micrometers) to small animals of a few millimeters. Eukaryotic ribosomal diversity saturated at ~150,000 operational taxonomic units, about one-third of which could not be assigned to known eukaryotic groups. Diversity emerged at all taxonomic levels, both within the groups comprising the ~11,200 cataloged morphospecies of eukaryotic plankton and among twice as many other deep-branching lineages of unappreciated importance in plankton ecology studies. Most eukaryotic plankton biodiversity belonged to heterotrophic protistan groups, particularly those known to be parasites or symbiotic hosts.

Ocean plankton. Determinants of community structure in the global plankton interactome.

Species interaction networks are shaped by abiotic and biotic factors. Here, as part of the Tara Oceans project, we studied the photic zone interactome using environmental factors and organismal abundance profiles and found that environmental factors are incomplete predictors of community structure. We found associations across plankton functional types and phylogenetic groups to be nonrandomly distributed on the network and driven by both local and global patterns. We identified interactions among grazers, primary producers, viruses, and (mainly parasitic) symbionts and validated network-generated hypotheses using microscopy to confirm symbiotic relationships. We have thus provided a resource to support further research on ocean food webs and integrating biological components into ocean models.

Metagenomic 16S rDNA Illumina tags are a powerful alternative to amplicon sequencing to explore diversity and structure of microbial communities.

Sequencing of 16S rDNA polymerase chain reaction (PCR) amplicons is the most common approach for investigating environmental prokaryotic diversity, despite the known biases introduced during PCR. Here we show that 16S rDNA fragments derived from Illumina-sequenced environmental metagenomes (mi tags) are a powerful alternative to 16S rDNA amplicons for investigating the taxonomic diversity and structure of prokaryotic communities. As part of the Tara Oceans global expedition, marine plankton was sampled in three locations, resulting in 29 subsamples for which metagenomes were produced by shotgun Illumina sequencing (ca. 700 Gb). For comparative analyses, a subset of samples was also selected for Roche-454 sequencing using both shotgun (m454 tags; 13 metagenomes, ca. 2.4 Gb) and 16S rDNA amplicon (454 tags; ca. 0.075 Gb) approaches. Our results indicate that by overcoming PCR biases related to amplification and primer mismatch, mi tags may provide more realistic estimates of community richness and evenness than amplicon 454 tags. In addition, mi tags can capture expected beta diversity patterns. Using mi tags is now economically feasible given the dramatic reduction in high-throughput sequencing costs, having the advantage of retrieving simultaneously both taxonomic (Bacteria, Archaea and Eukarya) and functional information from the same microbial community.