Dispersion of antimicrobial resistant bacteria in pig farms and in the surrounding environment.

BACKGROUND: Antimicrobial resistance has been identified as a major threat to global health. The pig food chain is considered an important source of antimicrobial resistance genes (ARGs). However, there is still a lack of knowledge on the dispersion of ARGs in pig production system, including the external environment. RESULTS: In the present study, we longitudinally followed one swine farm located in Italy from the weaning phase to the slaughterhouse to comprehensively assess the diversity of ARGs, their diffusion, and the bacteria associated with them. We obtained shotgun metagenomic sequences from 294 samples, including pig feces, farm environment, soil around the farm, wastewater, and slaughterhouse environment. We identified a total of 530 species-level genome bins (SGBs), which allowed us to assess the dispersion of microorganisms and their associated ARGs in the farm system. We identified 309 SGBs being shared between the animals gut microbiome, the internal and external farm environments. Specifically, these SGBs were characterized by a diverse and complex resistome, with ARGs active against 18 different classes of antibiotic compounds, well matching antibiotic use in the pig food chain in Europe. CONCLUSIONS: Collectively, our results highlight the urgency to implement more effective countermeasures to limit the dispersion of ARGs in the pig food systems and the relevance of metagenomics-based approaches to monitor the spread of ARGs for the safety of the farm working environment and the surrounding ecosystems.

Routes of dispersion of antibiotic resistance genes from the poultry farm system.

Poultry farms are hotspots for the development and spread of antibiotic resistance genes (ARGs), due to high stocking densities and extensive use of antibiotics, posing a threat of spread and contagion to workers and the external environment. Here, we applied shotgun metagenome sequencing to characterize the gut microbiome and resistome of poultry, workers and their households - also including microbiomes from the internal and external farm environment - in three different farms in Italy during a complete rearing cycle. Our results highlighted a relevant overlap among the microbiomes of poultry, workers, and their families (gut and skin), with clinically relevant ARGs and associated mobile elements shared in both poultry and human samples. On a finer scale, the reconstruction of species-level genome bins (SGBs) allowed us to delineate the dynamics of microorganism and ARGs dispersion from farm systems. We found the associations with worker microbiomes representing the main route of ARGs dispersion from poultry to human populations. Collectively, our findings clearly demonstrate the urgent need to implement more effective procedures to counteract ARGs dispersion from poultry food systems and the relevance of metagenomics-based metacommunity approaches to monitor the ARGs dispersion process for the safety of the working environment on farms.

Protocol for a multicentre cross-sectional, longitudinal ambulatory clinical trial in rheumatoid arthritis and Parkinson's disease patients analysing the relation between the gut microbiome, fasting and immune status in Germany (ExpoBiome).

INTRODUCTION: Chronic inflammatory diseases like rheumatoid arthritis (RA) and neurodegenerative disorders like Parkinson's disease (PD) have recently been associated with a decreased diversity in the gut microbiome, emerging as key driver of various diseases. The specific interactions between gut-borne microorganisms and host pathophysiology remain largely unclear. The microbiome can be modulated by interventions comprising nutrition.The aim of our clinical study is to (1) examine effects of prolonged fasting (PF) and time-restricted eating (TRE) on the outcome parameters and the immunophenotypes of RA and PD with (2) special consideration of microbial taxa and molecules associated with changes expected in (1), and (3) identify factors impacting the disease course and treatment by in-depth screening of microorganisms and molecules in personalised HuMiX gut-on-chip models, to identify novel targets for anti-inflammatory therapy. METHODS AND ANALYSIS: This trial is an open-label, multicentre, controlled clinical trial consisting of a cross-sectional and a longitudinal study. A total of 180 patients is recruited. For the cross-sectional study, 60 patients with PD, 60 patients with RA and 60 healthy controls are recruited at two different, specialised clinical sites. For the longitudinal part, 30 patients with PD and 30 patients with RA undergo 5-7 days of PF followed by TRE (16:8) for a period of 12 months. One baseline visit takes place before the PF intervention and 10 follow-up visits will follow over a period of 12 months (April 2021 to November 2023). ETHICS AND DISSEMINATION: Ethical approval was obtained to plan and conduct the trial from the institutional review board of the Charité-Universitätsmedizin Berlin (EA1/204/19), the ethics committee of the state medical association (Landesärztekammer) of Hessen (2021-2230-zvBO) and the Ethics Review Panel (ERP) of the University of Luxembourg (ERP 21-001 A ExpoBiome). The results of this study will be disseminated through peer-reviewed publications, scientific presentations and social media. TRIAL REGISTRATION NUMBER: NCT04847011.

Altered infective competence of the human gut microbiome in COVID-19.

BACKGROUND: Infections with SARS-CoV-2 have a pronounced impact on the gastrointestinal tract and its resident microbiome. Clear differences between severe cases of infection and healthy individuals have been reported, including the loss of commensal taxa. We aimed to understand if microbiome alterations including functional shifts are unique to severe cases or a common effect of COVID-19. We used high-resolution systematic multi-omic analyses to profile the gut microbiome in asymptomatic-to-moderate COVID-19 individuals compared to a control group. RESULTS: We found a striking increase in the overall abundance and expression of both virulence factors and antimicrobial resistance genes in COVID-19. Importantly, these genes are encoded and expressed by commensal taxa from families such as Acidaminococcaceae and Erysipelatoclostridiaceae, which we found to be enriched in COVID-19-positive individuals. We also found an enrichment in the expression of a betaherpesvirus and rotavirus C genes in COVID-19-positive individuals compared to healthy controls. CONCLUSIONS: Our analyses identified an altered and increased infective competence of the gut microbiome in COVID-19 patients. Video Abstract.

The gut microbiome molecular complex in human health and disease.

The human gut microbiome produces a functional complex of biomolecules, including nucleic acids, (poly)peptides, structural molecules, and metabolites. This impacts human physiology in multiple ways, especially by triggering inflammatory pathways in disease. At present, much remains to be learned about the identity of key effectors and their causal roles.

SnapShot: The Expobiome Map.

The human gut microbiome is intricately connected to health and disease. Microbiome-derived molecules are implicated in many chronic diseases involving inflammation. Herein, we summarize the diverse complex of such immunogenic molecules, including nucleic acids, (poly)peptides, structural molecules, and metabolites. The interactions between this "expobiome" and human immune pathways are specifically illustrated in the context of chronic diseases. To view this SnapShot, open or download the PDF.

Functional meta-omics provide critical insights into long- and short-read assemblies.

Real-world evaluations of metagenomic reconstructions are challenged by distinguishing reconstruction artifacts from genes and proteins present in situ. Here, we evaluate short-read-only, long-read-only and hybrid assembly approaches on four different metagenomic samples of varying complexity. We demonstrate how different assembly approaches affect gene and protein inference, which is particularly relevant for downstream functional analyses. For a human gut microbiome sample, we use complementary metatranscriptomic and metaproteomic data to assess the metagenomic data-based protein predictions. Our findings pave the way for critical assessments of metagenomic reconstructions. We propose a reference-independent solution, which exploits the synergistic effects of multi-omic data integration for the in situ study of microbiomes using long-read sequencing data.

PathoFact: a pipeline for the prediction of virulence factors and antimicrobial resistance genes in metagenomic data.

BACKGROUND: Pathogenic microorganisms cause disease by invading, colonizing, and damaging their host. Virulence factors including bacterial toxins contribute to pathogenicity. Additionally, antimicrobial resistance genes allow pathogens to evade otherwise curative treatments. To understand causal relationships between microbiome compositions, functioning, and disease, it is essential to identify virulence factors and antimicrobial resistance genes in situ. At present, there is a clear lack of computational approaches to simultaneously identify these factors in metagenomic datasets. RESULTS: Here, we present PathoFact, a tool for the contextualized prediction of virulence factors, bacterial toxins, and antimicrobial resistance genes with high accuracy (0.921, 0.832 and 0.979, respectively) and specificity (0.957, 0.989 and 0.994). We evaluate the performance of PathoFact on simulated metagenomic datasets and perform a comparison to two other general workflows for the analysis of metagenomic data. PathoFact outperforms all existing workflows in predicting virulence factors and toxin genes. It performs comparably to one pipeline regarding the prediction of antimicrobial resistance while outperforming the others. We further demonstrate the performance of PathoFact on three publicly available case-control metagenomic datasets representing an actual infection as well as chronic diseases in which either pathogenic potential or bacterial toxins are hypothesized to play a role. In each case, we identify virulence factors and AMR genes which differentiated between the case and control groups, thereby revealing novel gene associations with the studied diseases. CONCLUSION: PathoFact is an easy-to-use, modular, and reproducible pipeline for the identification of virulence factors, bacterial toxins, and antimicrobial resistance genes in metagenomic data. Additionally, our tool combines the prediction of these pathogenicity factors with the identification of mobile genetic elements. This provides further depth to the analysis by considering the genomic context of the pertinent genes. Furthermore, PathoFact's modules for virulence factors, toxins, and antimicrobial resistance genes can be applied independently, thereby making it a flexible and versatile tool. PathoFact, its models, and databases are freely available at https://pathofact.lcsb.uni.lu . Video abstract.

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.

Systematic characterization of human gut microbiome-secreted molecules by integrated multi-omics.

The human gut microbiome produces a complex mixture of biomolecules that interact with human physiology and play essential roles in health and disease. Crosstalk between micro-organisms and host cells is enabled by different direct contacts, but also by the export of molecules through secretion systems and extracellular vesicles. The resulting molecular network, comprised of various biomolecular moieties, has so far eluded systematic study. Here we present a methodological framework, optimized for the extraction of the microbiome-derived, extracellular biomolecular complement, including nucleic acids, (poly)peptides, and metabolites, from flash-frozen stool samples of healthy human individuals. Our method allows simultaneous isolation of individual biomolecular fractions from the same original stool sample, followed by specialized omic analyses. The resulting multi-omics data enable coherent data integration for the systematic characterization of this molecular complex. Our results demonstrate the distinctiveness of the different extracellular biomolecular fractions, both in terms of their taxonomic and functional composition. This highlights the challenge of inferring the extracellular biomolecular complement of the gut microbiome based on single-omic data. The developed methodological framework provides the foundation for systematically investigating mechanistic links between microbiome-secreted molecules, including those that are typically vesicle-associated, and their impact on host physiology in health and disease.

Dichloromethane Degradation Pathway from Unsequenced Hyphomicrobium sp. MC8b Rapidly Explored by Pan-Proteomics.

Several bacteria are able to degrade the major industrial solvent dichloromethane (DCM) by using the conserved dehalogenase DcmA, the only system for DCM degradation characterised at the sequence level so far. Using differential proteomics, we rapidly identified key determinants of DCM degradation for Hyphomicrobium sp. MC8b, an unsequenced facultative methylotrophic DCM-degrading strain. For this, we designed a pan-proteomics database comprising the annotated genome sequences of 13 distinct Hyphomicrobium strains. Compared to growth with methanol, growth with DCM induces drastic changes in the proteome of strain MC8b. Dichloromethane dehalogenase DcmA was detected by differential pan-proteomics, but only with poor sequence coverage, suggesting atypical characteristics of the DCM dehalogenation system in this strain. More peptides were assigned to DcmA by error-tolerant search, warranting subsequent sequencing of the genome of strain MC8b, which revealed a highly divergent set of dcm genes in this strain. This suggests that the dcm enzymatic system is less strongly conserved than previously believed, and that substantial molecular evolution of dcm genes has occurred beyond their horizontal transfer in the bacterial domain. Our study showed the power of pan-proteomics for quick characterization of new strains belonging to branches of the Tree of Life that are densely genome-sequenced.

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.

Integrating Culture-based Antibiotic Resistance Profiles with Whole-genome Sequencing Data for 11,087 Clinical Isolates.

Emerging antibiotic resistance is a major global health threat. The analysis of nucleic acid sequences linked to susceptibility phenotypes facilitates the study of genetic antibiotic resistance determinants to inform molecular diagnostics and drug development. We collected genetic data (11,087 newly-sequenced whole genomes) and culture-based resistance profiles (10,991 out of the 11,087 isolates comprehensively tested against 22 antibiotics in total) of clinical isolates including 18 main species spanning a time period of 30 years. Species and drug specific resistance patterns were observed including increased resistance rates for Acinetobacter baumannii to carbapenems and for Escherichia coli to fluoroquinolones. Species-level pan-genomes were constructed to reflect the genetic repertoire of the respective species, including conserved essential genes and known resistance factors. Integrating phenotypes and genotypes through species-level pan-genomes allowed to infer gene-drug resistance associations using statistical testing. The isolate collection and the analysis results have been integrated into GEAR-base, a resource available for academic research use free of charge at https://gear-base.com.

Assessing the heterogeneity of in silico plasmid predictions based on whole-genome-sequenced clinical isolates.

High-throughput next-generation shotgun sequencing of pathogenic bacteria is growing in clinical relevance, especially for chromosomal DNA-based taxonomic identification and for antibiotic resistance prediction. Genetic exchange is facilitated for extrachromosomal DNA, e.g. plasmid-borne antibiotic resistance genes. Consequently, accurate identification of plasmids from whole-genome sequencing (WGS) data remains one of the major challenges for sequencing-based precision medicine in infectious diseases. Here, we assess the heterogeneity of four state-of-the-art tools (cBar, PlasmidFinder, plasmidSPAdes and Recycler) for the in silico prediction of plasmid-derived sequences from WGS data. Heterogeneity, sensitivity and precision were evaluated by reference-independent and reference-dependent benchmarking using 846 Gram-negative clinical isolates. Interestingly, the majority of predicted sequences were tool-specific, resulting in a pronounced heterogeneity across tools for the reference-independent assessment. In the reference-dependent assessment, sensitivity and precision values were found to substantially vary between tools and across taxa, with cBar exhibiting the highest median sensitivity (87.45%) but a low median precision (27.05%). Furthermore, integrating the individual tools into an ensemble approach showed increased sensitivity (95.55%) while reducing the precision (25.62%). CBar and plasmidSPAdes exhibited the strongest concordance with respect to identified antibiotic resistance factors. Moreover, false-positive plasmid predictions typically contained only few antibiotic resistance factors. In conclusion, while high degrees of heterogeneity and variation in sensitivity and precision were observed across the different tools and taxa, existing tools are valuable for investigating the plasmid-borne resistome. Nevertheless, additional studies on representative clinical data sets will be necessary to translate in silico plasmid prediction approaches from research to clinical application.

Comparison of initial oral microbiomes of young adults with and without cavitated dentin caries lesions using an in situ biofilm model.

Dental caries is caused by acids released from bacterial biofilms. However, the in vivo formation of initial biofilms in relation to caries remains largely unexplored. The aim of this study was to compare the oral microbiome during the initial phase of bacterial colonization for individuals with (CC) and without (NC) cavitated dentin caries lesions. Bovine enamel slabs on acrylic splints were worn by the volunteers (CC: 14, NC: 13) for in situ biofilm formation (2 h, 4 h, 8 h, 1 ml saliva as reference). Sequencing of the V1/V2 regions of the 16S rRNA gene was performed (MiSeq). The relative abundances of individual operational taxonomic units (OTUs) were compared between samples from the CC group and the NC group. Random forests models were furthermore trained to separate the groups. While the overall heterogeneity did not differ substantially between CC and NC individuals, several individual OTUs were found to have significantly different relative abundances. For the 8 h samples, most of the significant OTUs showed higher relative abundances in the CC group, while the majority of significant OTUs in the saliva samples were more abundant in the NC group. Furthermore, using OTU signatures enabled a separation between both groups, with area-under-the-curve (AUC) values of ~0.8. In summary, the results suggest that initial oral biofilms provide the potential to differentiate between CC and NC individuals.

Comparing genome versus proteome-based identification of clinical bacterial isolates.

Whole-genome sequencing (WGS) is gaining importance in the analysis of bacterial cultures derived from patients with infectious diseases. Existing computational tools for WGS-based identification have, however, been evaluated on previously defined data relying thereby unwarily on the available taxonomic information.Here, we newly sequenced 846 clinical gram-negative bacterial isolates representing multiple distinct genera and compared the performance of five tools (CLARK, Kaiju, Kraken, DIAMOND/MEGAN and TUIT). To establish a faithful 'gold standard', the expert-driven taxonomy was compared with identifications based on matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) mass spectrometry (MS) analysis. Additionally, the tools were also evaluated using a data set of 200 Staphylococcus aureus isolates.CLARK and Kraken (with k =31) performed best with 626 (100%) and 193 (99.5%) correct species classifications for the gram-negative and S. aureus isolates, respectively. Moreover, CLARK and Kraken demonstrated highest mean F-measure values (85.5/87.9% and 94.4/94.7% for the two data sets, respectively) in comparison with DIAMOND/MEGAN (71 and 85.3%), Kaiju (41.8 and 18.9%) and TUIT (34.5 and 86.5%). Finally, CLARK, Kaiju and Kraken outperformed the other tools by a factor of 30 to 170 fold in terms of runtime.We conclude that the application of nucleotide-based tools using k-mers-e.g. CLARK or Kraken-allows for accurate and fast taxonomic characterization of bacterial isolates from WGS data. Hence, our results suggest WGS-based genotyping to be a promising alternative to the MS-based biotyping in clinical settings. Moreover, we suggest that complementary information should be used for the evaluation of taxonomic classification tools, as public databases may suffer from suboptimal annotations.

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.

Integrated meta-omic analyses of the gastrointestinal tract microbiome in patients undergoing allogeneic hematopoietic stem cell transplantation.

In patients undergoing allogeneic hematopoietic stem cell transplantation (allo-HSCT), treatment-induced changes to the gastrointestinal tract (GIT) microbiome have been linked to adverse outcomes, most notably graft-versus-host disease (GvHD). However, it is presently unknown whether this relationship is causal or consequential. Here, we performed an integrated meta-omic analysis to probe deeper into the GIT microbiome changes during allo-HSCT and its accompanying treatments. We used 16S and 18S rRNA gene amplicon sequencing to resolve archaea, bacteria, and eukaryotes within the GIT microbiomes of 16 patients undergoing allo-HSCT for the treatment of hematologic malignancies. These results revealed a major shift in the GIT microbiome after allo-HSCT including a marked reduction in bacterial diversity, accompanied by only limited changes in eukaryotes and archaea. An integrated analysis of metagenomic and metatranscriptomic data was performed on samples collected from a patient before and after allo-HSCT for acute myeloid leukemia. This patient developed severe GvHD, leading to death 9 months after allo-HSCT. In addition to drastically decreased bacterial diversity, the post-treatment microbiome showed a higher overall number and higher expression levels of antibiotic resistance genes (ARGs). One specific Escherichia coli strain causing a paravertebral abscess was linked to GIT dysbiosis, suggesting loss of intestinal barrier integrity. The apparent selection for bacteria expressing ARGs suggests that prophylactic antibiotic administration may adversely affect the overall treatment outcome. We therefore assert that such analyses including information about the selection of pathogenic bacteria expressing ARGs may assist clinicians in "personalizing" regimens for individual patients to improve overall outcomes.

BusyBee Web: metagenomic data analysis by bootstrapped supervised binning and annotation.

Metagenomics-based studies of mixed microbial communities are impacting biotechnology, life sciences and medicine. Computational binning of metagenomic data is a powerful approach for the culture-independent recovery of population-resolved genomic sequences, i.e. from individual or closely related, constituent microorganisms. Existing binning solutions often require a priori characterized reference genomes and/or dedicated compute resources. Extending currently available reference-independent binning tools, we developed the BusyBee Web server for the automated deconvolution of metagenomic data into population-level genomic bins using assembled contigs (Illumina) or long reads (Pacific Biosciences, Oxford Nanopore Technologies). A reversible compression step as well as bootstrapped supervised binning enable quick turnaround times. The binning results are represented in interactive 2D scatterplots. Moreover, bin quality estimates, taxonomic annotations and annotations of antibiotic resistance genes are computed and visualized. Ground truth-based benchmarks of BusyBee Web demonstrate comparably high performance to state-of-the-art binning solutions for assembled contigs and markedly improved performance for long reads (median F1 scores: 70.02-95.21%). Furthermore, the applicability to real-world metagenomic datasets is shown. In conclusion, our reference-independent approach automatically bins assembled contigs or long reads, exhibits high sensitivity and precision, enables intuitive inspection of the results, and only requires FASTA-formatted input. The web-based application is freely accessible at: https://ccb-microbe.cs.uni-saarland.de/busybee.

Colonization and Succession within the Human Gut Microbiome by Archaea, Bacteria, and Microeukaryotes during the First Year of Life.

Perturbations to the colonization process of the human gastrointestinal tract have been suggested to result in adverse health effects later in life. Although much research has been performed on bacterial colonization and succession, much less is known about the other two domains of life, archaea, and eukaryotes. Here we describe colonization and succession by bacteria, archaea and microeukaryotes during the first year of life (samples collected around days 1, 3, 5, 28, 150, and 365) within the gastrointestinal tract of infants delivered either vaginally or by cesarean section and using a combination of quantitative real-time PCR as well as 16S and 18S rRNA gene amplicon sequencing. Sequences from organisms belonging to all three domains of life were detectable in all of the collected meconium samples. The microeukaryotic community composition fluctuated strongly over time and early diversification was delayed in infants receiving formula milk. Cesarean section-delivered (CSD) infants experienced a delay in colonization and succession, which was observed for all three domains of life. Shifts in prokaryotic succession in CSD infants compared to vaginally delivered (VD) infants were apparent as early as days 3 and 5, which were characterized by increased relative abundances of the genera Streptococcus and Staphylococcus, and a decrease in relative abundance for the genera Bifidobacterium and Bacteroides. Generally, a depletion in Bacteroidetes was detected as early as day 5 postpartum in CSD infants, causing a significantly increased Firmicutes/Bacteroidetes ratio between days 5 and 150 when compared to VD infants. Although the delivery mode appeared to have the strongest influence on differences between the infants, other factors such as a younger gestational age or maternal antibiotics intake likely contributed to the observed patterns as well. Our findings complement previous observations of a delay in colonization and succession of CSD infants, which affects not only bacteria but also archaea and microeukaryotes. This further highlights the need for resolving bacterial, archaeal, and microeukaryotic dynamics in future longitudinal studies of microbial colonization and succession within the neonatal gastrointestinal tract.