mSystems

Introduction. Polycyclic aromatic hydrocarbon (PAH) toxicology has been primarily defined by host metabolic bioactivation to reactive intermediates. Yet ingested PAHs and biliary-delivered metabolites also enter the intestinal lumen, where they encounter the gut microbiota, an additional but understudied metabolic compartment that may alter toxicant fate. Although PAH interactions are well characterized in environmental microbes, resulting in defined enzymatic pathways for their catabolism, whether gut commensal organisms express related machinery in vivo remains poorly understood. To address this gap, we used an oral murine-exposure model of naphthalene and integrated compartment-resolved 16S rRNA gene profiling, global metaproteomics, and chemoproteomics to resolve the PAH-processing repertoire of the gut microbiota and its spatial distribution. Results. Naphthalene, a common environmental toxicant and representative PAH, elicited a compartment-dependent functional response that was most pronounced in the large intestine. Distally, exposure was associated with 2,386 differential proteins and accounted for 48.4% of metaproteomic variation, compared with 7.9% of 16S rRNA gene compositional variation. The small intestine showed a reduced proteomic response, with 994 differentially abundant proteins. Hydrocarbon catabolic markers showed a 4.8-fold distal expansion (36 to 172 proteins), including the naphthalene-catabolic nahE, proteins associated with phenylacetate, {beta}-ketoadipate, catechol/TOL, and anaerobic CoA-linked redox machinery. The small-intestinal repertoire decreased from 19 to 14 proteins, though aromatic-catabolic markers doubled, indicating a selective expansion. Critically, these homologs concentrated in Lachnospiraceae, a gut family that is prevalent in humans, primarily known for anaerobic fermentation, with no priorly described PAH-processing capacity. An uncultured Acetatifactor showed the largest hydrocarbon hits (14 to 82 protein increase), with additional contributions from Kineothrix, Eisenbergiella, and Petralouisia. Chemoproteomics, using a naphthalene-based probe, revealed enrichment of an orthogonal proteome of proteins interacting with naphthalene, capturing 171 small-intestinal and 67 large-intestinal proteins that were undetected in the global metaproteomes. Conclusions. Gut commensals express homologs of aromatic-hydrocarbon catabolic enzymes characterized in environmental microorganisms, and naphthalene amplifies this repertoire most strongly in the large intestine. The decoupling of taxonomic structure, expressed-protein representation, and probe engagement identifies a functional dimension of PAH exposure that composition-based profiling does not resolve and positions microbiome variation as a candidate modifier of PAH metabolic fate.

Schistosoma mansoni infection is associated with alterations of the host gut microbial community structure, yet its impact on microbial genome dynamics and horizontal gene transfer remains poorly understood. In this study, we investigated the activity of Mobile Genetic Elements (MGEs), and of Insertion Sequences (IS) in particular, in the gut microbiome of wild-type (WT) and human-microbiota-associated (HMA) mice. Taxonomic profiling revealed significant differences in bacterial diversity between WT and HMA mice, with reduced richness and alpha diversity observed in HMA gut microbiomes regardless of the infection status. Principal Coordinates Analysis demonstrated clear separation between microbial communities, reflecting differences in host background and bacterial composition. In contrast, differences were observed between plasmid and viral community profiles and between these and bacterial diversity, suggesting that different MGEs may follow distinct ecological dynamics within the gut microbiome. IS elements within bacterial genomes were widespread across all microbiome communities analysed in this study, with higher insertion counts observed not only in HMA compared to WT mice, but also in infected compared to uninfected mice. IS elements were predominantly associated with dominant gut bacterial taxa, particularly members of Bacteroidia and Clostridia. At the genomic level, insertion events occurred more frequently in intergenic regions than in coding sequences, suggesting a preferential selective pressure towards alteration of gene expression rather than gene disruption. Furthermore, IS activity was dominated by a subset of IS families, with consistent patterns observed across experimental groups. Notably, infection was associated with increased insertion frequencies across multiple taxa and IS families, suggesting enhanced mobilisation under conditions of host-associated stress (the physiological and molecular pressure experienced by microorganisms, pathogens, or the host itself during their interactions within a symbiotic or infected environment) and microbiome perturbation. Together, these findings demonstrate that S. mansoni infection may promote IS mobilisation in the mouse gut microbiome, potentially contributing to microbial genome plasticity and facilitating horizontal gene transfer within host-associated microbial communities.

Mesoplasma florum is a fast-growing, near-minimal bacterium and an emerging model for systems and synthetic biology. However, its dependence on complex serum-containing media limits experimental control and complicates the interpretation of cellular phenotypes. Here, we developed CMRL-AT, a serum-free, quasi-defined medium that supports rapid growth comparable to the commonly used ATCC 1161 medium. Despite supporting similar biomass, CMRL-AT profoundly reshaped the M. florum transcriptome, with approximately one-third of the annotated protein-encoding genes being differentially expressed relative to ATCC 1161. These changes revealed distinct physiological programs associated with rapid growth in complex medium and higher nutrient acquisition in CMRL-AT, illustrating how medium composition alters the functional priorities of a near-minimal cell. Transcriptome profiling across six energy sources further uncovered distinct sugar-responsive expression programs. Combining these responses with transcription-unit organization, protein-domain predictions, and metabolic context resolved fructose- and sucrose-responsive modules, and allowed the assignment of previously ambiguous phosphotransferase system components to specific sugar-utilization pathways. CMRL-AT provides an experimental framework to help resolving gene functions, refining metabolic models, and designing reduced genomes adapted to defined environments.

Spontaneous vegetable fermentation is driven by the indigenous microbiota associated with the raw ingredients. The process often involves an initial heterofermentative phase followed by a predominantly homofermentative phase, both ensured by lactic acid bacteria (LAB). While microbial dynamics during fermentation are well characterized, phage communities remain poorly understood, despite evidence linking phage activity to LAB succession. We monitored viral and bacterial succession during sauerkraut and carrot fermentations using a metagenomics approach to characterize community dynamics. To further assess the impact of phage PLE001 on bacterial community dynamics, we developed a simplified carrot fermentation model comprising representative strains of Rahnella aquatilis, Leuconostoc mesenteroides and Lactiplantibacillus plantarum. Two consortia were assembled: one containing the PLE001-susceptible strain L. mesenteroides CIM3004, and the other containing also its phage-resistant derivative, CIM3004 BIM. We observed major shifts in the viral community during the first week of fermentation, coinciding with changes in bacterial structure during hetero- and homofermentative phases transition. In the simplified model, phage PLE001 reduced the population of the host during the heterofermentative phase but did not affect overall fermentation dynamics nor metabolite production. These results suggest that bacteriophages contribute to bacterial succession during fermentation, while communities remain resilient to phage infection, preserving fermentation stability.

Uropathogenic Escherichia coli (UPEC) is the predominant causative agent of urinary tract infections. The UPEC strain, CFT073, has a mosaic genome containing 60 genomic islands enriched in virulence genes. In this study, we explore the small proteome of CFT073 using ribosome profiling (Ribo-seq) coupled with RNA sequencing (RNA-seq), and reveal 338 translated, previously unannotated, putative small open reading frames (sORFs) with less than 100 codons. Among them, 78 sORFs are located in intergenic regions. Fifteen intergenic sORFs encode a predicted transmembrane helix (smemORFs), of which seven are hypothetical or have not been characterized in other E. coli strains. The majority of the remaining putative sORFs encodes uncharacterized or novel soluble proteins (ssolORFs). Considering the evolution of many sORFs is likely to be recent, functional sORFs are more likely to have conserved gene neighborhood in other genomes. Searching into the genomes of 10 other E. coli strains and four enterobacteria species, synteny analysis reveals that 10 of 15 putative smemORFs, and 27 of 63 putative ssolORFs, have a homologous sequence in a similar genomic context in at least one other genome, supporting their likelihood of encoding functional proteins. Further validation of putative sORFs using proteomics, cell-free protein synthesis, and mNeonGreen tagging provide strong evidence for the translation of nine predicted intergenic sORFs. Notably, four previously uncharacterized small membrane proteins have increased abundance under oxidative stress, suggesting a potential role in the stress response. Together, our findings reveal the small proteome of an extraintestinal pathogenic E. coli and provide a catalog of putative sORFs for further investigation.

Bacterial genomes are dynamic and harbour accessory regions that can be gained, lost, or rearranged in response to selective pressures. Members of the genus Komagataeibacter are used for bacterial cellulose production and researched as platforms for engineered living materials. However, their genomes are rich in mobile genetic elements (MGEs), including insertion sequences (ISs), genomic islands and prophage remnants, which drive genome rearrangements under agitated cultivation conditions. Whether this genome architecture contributes to genomic changes during routine strain maintenance and distribution remains poorly understood. Comparative genome sequencing of independent Komagataeibacter rhaeticus iGEM stocks revealed recurrent loss of a 210 kbp chromosomal region flanked by homologous IS5 sequences. This deletion was identified in a stock received as an agar stab and reproduced from an intact-genome population subjected to simulated agar-stab handling. Phylogenomic analysis of the affected locus revealed gene-tree/species-tree discordance, heterogeneous GC content and multiple ISs, supporting its characterization as a dynamic accessory region. Populations carrying the larger genome reduction showed improved growth and cellulose production relative to the intact genome reference, whereas populations carrying a smaller nested deletion showed reduced performance. These findings demonstrate that routine strain handling is not necessarily genomically neutral and can expose structural genomic variabilities that generate population-level genomic heterogeneity with phenotypic consequences. For bacteria with MGE-rich genomes, handling history may represent an underappreciated source of genomic and phenotypic variation between laboratories. More broadly, identifying unstable accessory regions may highlight genomic content that is dispensable under laboratory conditions, providing candidates for genome streamlining and strain improvement.

Vietnam's rapid socioeconomic and dietary transition offers a rare window into how industrialization is reflected in the gut microbiome, yet Vietnamese populations remain underrepresented in global reference datasets. We investigated whether healthy Vietnamese adults in Hanoi harbor microbial community configurations associated with different stages of industrialization. Metagenomic data from stool samples of 183 clinically screened participants in Hanoi were analyzed alongside 629 publicly available metagenomic datasets from populations in low to middle income (LMI) and high income (HI) countries. Species-level taxonomic composition and gene-content-based metabolic pathways were characterized using community clustering and multivariable association analyses. Rather than forming a single population profile, the Vietnamese cohort separated into two enterotypes that mirrored the global industrialization gradient. Enterotype 1 more closely resembled populations from LMICs and showed greater microbial richness and diversity, whereas Enterotype 2 aligned with taxonomic and functional profiles observed in HICs. Both taxonomic and functional potential shows similar trajectories. Diet and lifestyle patterns describe some of the variance between enterotypes, supporting an association between contemporary exposures and microbiome structure without establishing causation. Antibiotic resistance profiles also varied between enterotypes, suggesting that the transition may encompass functions relevant to antimicrobial exposure. Together, these findings position urban Vietnam at a microbiome transition, where community states associated with different levels of industrialization coexist within a single population. This Vietnamese reference cohort expands representation in global metagenomic research and provides a baseline for tracking how rapid socioeconomic change may reshape microbial ecology, metabolic capacity, and antimicrobial resistance.

The genomic and epigenomic complexity of microbiomes with low species richness remains poorly characterized. Here, we integrated 111 non-redundant isolate genomes with deep long-read metagenomics to investigate an undefined mesophilic dairy starter. Closely related isolates harboured distinct combinations of features associated with fermentation traits and phage defence, including surface polysaccharide genotypes, plasmids, prophages, defence systems, and methylation profiles. Metagenomics revealed additional elements and inferred host associations. However, the limited overlap between approaches across core-genome lineages and accessory repertoires underscored complementary but incomplete recovery, to the extent that the dominant metagenomic cell wall polysaccharide genotype (29.6% abundance) was absent from isolates. Metagenomic methylation signatures supported 36 species-level plasmid-host assignments, corroborated by isolate genomes in nine of ten testable cases. Our findings revealed extensive combinatorial diversity hidden from either approach alone, providing a plausible molecular basis for phage resilience, while likely representing only a lower bound on population complexity.

Fecal-derived in vitro microbial communities represent promising models for investigating gut microbiome function and facilitating translational applications. However, the factors governing community establishment and the underlying mechanisms remain incompletely understood. We employed large-scale quantitative metaproteomics to investigate the functional and community dynamics of human fecal microbiomes cultured in three media over 25 serial passages. Our results demonstrated that both culture medium and donor-specific microbiome characteristics significantly influence community establishment. Taxon-specific functional analyses revealed nutrient availability in the medium and the microbial stress responses as key determinants of in vitro community establishment. Notably, the metabolism of mucin glycans in RapidAIM medium contributed to the maintenance of Bacteroidaceae species that were poorly supported in other media, resulting in an in vitro community that most closely resembled the uncultured fecal inoculum. Overall, these findings provide mechanistic insights into in vitro gut microbial community establishment, supporting further optimization and application of fecal-derived in vitro microbiome models.

Background: Microbiome research has linked microbial community structure to host phenotype, but translating microbiome-based discovery into defined interventions with reproducible long-term host effects remains challenging. We evaluated whether a defined four-member consortium of rumen native microorganisms identified through associations with host phenotype could influence herd retention. Cull outcomes (sold or died vs. retained in herd) were analyzed across 20 controlled commercial dairy trials encompassing more than 40,000 cattle housed across 110 experimental pens, with each trial conducted for at least 5 months. Treatment pens received a defined daily dose of Galaxis Frontier (GF, Native Microbials, San Diego, CA), containing Pichia kudriavzevii, Clostridium beijerinckii, Ruminococcus bovis, and Butyrivibrio fibrisolvens, incorporated into the ration; control and treatment groups were otherwise managed equivalently. Cow-level time to culling was modeled within each trial, accounting for pen-level clustering, and trial estimates were combined using random-effects meta-analysis. Results: Daily treatment with rumen native microbes was associated with a 9% lower hazard of culling across 20 trials (hazard ratio = 0.91, 95% confidence interval: 0.86-0.96; P < 0.001). Between trial heterogeneity was negligible, and treatment effects were not significantly moderated by breed or geographic region. Detailed health records from one trial were examined to further characterize health outcomes associated with GF; treatment was associated with fewer transition-related disorders, particularly ketosis and metritis, whereas mastitis and lameness were largely unaffected. These findings suggest that the health effects associated with GF may be more pronounced for disorders linked to physiological and metabolic adaptation during the transition period than for health events more broadly. Conclusion: Across diverse commercial dairy environments, treatment with a defined consortium of rumen native microbes was associated with improved dairy cow retention. The consistency of this association across independent trials provides large scale evidence that targeted application of native microorganisms can translate into reproducible, long-term outcomes for the host. Complementary health event data suggest that improved transition period health and metabolic stability may contribute to this effect. These findings support the potential for microbiome targeted interventions to influence productive longevity under commercial conditions.

Abstract: Background Despite the ecological importance of viral anti-defense systems (ADS), yet how viral ADS repertoires vary across microbial ecosystems remains poorly understood. Here, we integrated global virus datasets with metagenomes from heavy metal-contaminated soils and chromium-stress microcosms to investigate ecological determinants of viral ADS. Results We found that ADS repertoires were constrained by viral genome capacity, with most ADS-carrying viruses encoded a single system, whereas multi-ADS occurred mainly in larger genomes. Viral lifestyle partitioned ADS, with temperate viruses enriched in ADS targeting early-stage defenses and virulent viruses preferentially encoding ADS against mid- and late-stage defense. ADS composition also varied by host lineage and across habitats, linking viral ADS repertoires to host immune landscapes and ecosystem-specific viral strategies. Moreover, along the heavy metal contamination gradient, changes in the prevalence of ADS-carrying viruses occurred alongside shifts in prokaryotic defense profiles. In addition, expanded ADS repertoires may further increase dissemination potential of MGE-borne adaptive genes, highlighting viral ADS as a potential contributor to the adaptive evolution of microbial community. Conclusions These findings reveal that viral ADS are shaped by genomic traits, host immunity and environmental stress, with implications for virus-host interaction and microbial adaptive evolution. More broadly, this study extends anti-defense research beyond individual virus-host systems, offering an integrative, community-level perspective to guide future studies testing whether and how virus-prokaryote immune interactions contribute to microbial community dynamics and evolution in changing environments.

Background: Plasmids are key drivers of horizontal gene transfer (HGT), enabling the dissemination of accessory traits that shape microbial adaptation and ecological interactions. In Haloarchaea-dominant members of hypersaline environments-characterizing plasmidomes remains particularly challenging because most available genomes are incomplete, leaving many contigs unassigned to either chromosomal or plasmidic origin. This limitation hampers our ability to reconstruct archaeal plasmid diversity and to evaluate how plasmid-borne genes circulate within Haloarchaea and, potentially, outside this class. Results: Here, we present HaloClassifier, the first machine-learning tool specifically designed to distinguish plasmid- and chromosome-derived contigs in haloarchaeal genomes. Unlike existing bacterial-focused approaches, HaloClassifier uniquely integrates haloarchaeal-specific genomic signatures with coding-derived features-an underused but highly informative class of predictors-and incorporates a variable-selection strategy that substantially reduces model complexity and computational cost while maintaining high predictive accuracy. When trained on simulated contigs generated from complete haloarchaeal genomes, the model achieves 96.61% accuracy on extralarge contigs (40-100 kb), increasing to 98.1% at a 0.6 classification threshold, with only 3.55% of contigs left unclassified. For small contigs (1-5 kb), the model achieves 79.51% accuracy overall and 84.42% at the 0.6 classification threshold, with 17.07% of contigs left unclassified. The model also maintains robust performance in metagenomic datasets. Conclusions: By enabling reliable plasmid identification in fragmented assemblies, HaloClassifier fills a major gap and establishes the foundation for large-scale haloarchaeal plasmidome reconstruction. This framework will support future investigations into gene mobility, ecological adaptation, and the evolutionary dynamics of plasmids and HGT in Haloarchaea.

The growth and performance of cattle is dependent on efficient feed conversion, facilitated by the gut microbiome, including both prokaryotic and eukaryotic organisms. Extrinsic factors such as management practices and intrinsic factors like age and genetics can shape these communities. To characterize the relative effects of such factors on the prokaryotic and eukaryotic microbiomes, fresh fecal samples from dairy heifers (n = 8), beef heifers (n = 8), bred beef heifers (n = 8), and adult beef cows (n = 8) were subjected to paired 16S rRNA (V4) and a novel 18S rRNA (V4-V5) metabarcoding assay, and 18S rRNA amplicons were annotated using two databases in parallel. Taxonomic resolution of 18S rRNA amplicons was excellent, with SILVA and EUKARYOME databases each providing better resolution of different groups of organisms. While no differences in eukaryotic richness or alpha diversity were detected, beta diversity differed in all pairwise comparisons except between beef heifers and bred beef heifers. Prokaryotic richness was lower in adult beef cows than the other groups, and beta diversity differed in all pairwise comparisons. Adult beef cows were the most distinct, characterized by replacement of fungal organisms with protozoa and loss of thermophilic bacteria in the eukaryotic and prokaryotic compartments, respectively. Among numerous differentially abundant features, dairy heifers harbored greater abundance of Cooperia oncophora while beef heifers contained greater abundance of Buxtonella sulcata and Mucor. Correlation network analysis revealed multiple significant positive and negative associations between bacterial, fungal, and protozoal features during maturation of the fecal microbiome. Collectively, these data highlight the richness and diversity of both the prokaryotic and eukaryotic compartments of the bovine fecal microbiome, and the influence of management practices and age on both communities. From a technical standpoint, these data demonstrate the ability to resolve fungal, protozoal, and metazoan organisms, including potential parasites and pathogens, with a single metabarcoding assay. They also highlight the nuanced differences between databases and the value of parallel analyses using more than one database.

Microbiome has co-evolved with insects over millions of years, establishing symbiotic relationships that influence their fitness and shape evolutionary trajectories. These long-term associations are shaped by the hosts ecology, diet, and physiology which determines microbial diversity. In turn, external stresses such as pathogenic infections can adversely affect the host microbiome, further compounding the effects of pathogenic infections on the host. Using a computational pipeline and publicly available transcriptomics datasets, we analyzed the effect of viral infections on transcript derived microbiome diversity and composition in various insects including Drosophila melanogaster, Acyrthosiphon pisum, Antheraea pernyi, and Culex pipiens. The meta-transcriptomic reanalysis provides an indirect but comprehensive comparison of infection-associated transcriptional microbiome changes across hosts with distinct life histories. Our study reveals that viral infections influence the host-associated microbial transcripts, leading to shifts in transcript-derived microbial profiles and offering insights into host-microbe dynamics. In Cx. pipiens, microbial transcript profiles remained relatively stable following RVFV exposure and were dominated by the bacterial endosymbiont Wolbachia. Our study highlights viral infections as key drivers of infection-associated shifts in microbial profiles and suggests a potential association between Wolbachia dominance and transcriptional microbiome stability in insects, with age as a contributing factor in shaping microbial transcript profiles in Culex. This cross-species framework uncovers conserved and host-specific patterns of virus-microbiome interactions, expanding our understanding of microbial resilience mechanisms in insects.

During phage infection, phages reprogram their host to secure energy, building blocks for their genome, and proteins, as well as the transcription and translation machinery required to replicate their genome and produce phage particles. This also includes reprogramming the host transcriptome to modify host gene expression and align it with the needs of the phage during infection. In this study, we present a putative early-expressed protein, Gp204 of Mosigvirus AV110, renamed Host transcription modulator (Htm), that has a broad impact on the extended-spectrum beta-lactamase (ESBL)- associated E. coli transcriptome. Overexpression of Htm decreased the growth rate of ESBL-E. coli and caused DNA condensation. Htm is structurally similar to bacterial enhancer-binding proteins, binds ssDNA, and alters the expression levels of over 700 E. coli genes when expressed. The majority of these genes were upregulated and are involved in multiple pathways, such as galactose metabolism and pentose and glucuronate interconversions, which may promote the synthesis of purines, pyrimidines, and histidine, all of which are required for the replication, transcription, and translation of the phage genome. Furthermore, energy production is likely to be upregulated by increased availability of acetyl-CoA, resulting from, among other factors, fatty acid degradation. Htm also downregulated the expression of the structural genes involved in flagellar assembly and multiple two-component systems. In conclusion, Htm could promote phage progeny production by increasing the availability of nucleotides, amino acids, and energy, while inhibiting non-essential processes, such as flagella formation.

Microbial metabolism in extreme environments has long been shrouded in profound mystery. Distinct from the direct assimilation of exogenous organic carbon, microorganisms predominantly exploit inorganic nitrogen substrates as their primary nitrogen sources. As a critical nitrogen assimilation pathway, assimilatory nitrite reduction plays a vital role in microbial nitrogen acquisition, especially in extreme environments with poor nutrient availability. In this study, we identified the coexistence of two isoforms of NADH-dependent nitrite reductase in single genomes. NIR- and NIR-{beta}, named based on their phylogenetic relationships, exhibit significant differences in protein sequence, genomic context and evolutionary relationship, a phenomenon that has rarely been reported before. They also show distinct expression patterns, which may allow bacteria to achieve efficient nitrogen conversion in nutrient-poor environments. This finding calls for a revision of the conventional model of nitrogen cycling gene combinations inferred from common model microorganisms and holds numerous implications for how we understand and study the environmental adaptation mechanisms of extremophiles.

The gut-microbiome-brain axis (GMBA) links the gut and brain through metabolic, immunological, and neuroendocrine pathways, and can become dysregulated during psychological stress (PS). While previous studies primarily infer microbial function and use physical stressors that incompletely model human PS, here, functional multi-omics was applied to a chronic social-isolation mouse model. Although isolation caused no measurable behavioral changes, brain proteomics revealed stress-associated molecular responses. Taxonomic and short-chain fatty acid analyses suggested limited microbiome changes; however, metaproteomics uncovered sex-specific restructuring of bacterial ecology and carbohydrate-active enzymes. In females, isolation increased Petralouisia muris proteins involved in bacterial microcompartment metabolism and {beta}-glucuronidases across Lachnospiraceae taxa. Enriched {beta}-glucuronidase abundance coincided with reduced circulating glucuronides and altered brain estrogen-receptor activity, suggesting a microbial deconjugation-host estrogen-signaling interaction. These findings demonstrate that functional GMBA alterations may precede or occur independently of overt behavioral phenotypes and highlight the value of direct functional profiling in PS.

Bacteriophages exhibit staggering genomic diversity, yet phenotypic characterization remains bottlenecked by assays that generally yield simple binary infection outcomes. Quantitative traits can be extracted from standard liquid cultures, but many phages fail to produce a measurable effect in these settings, necessitating an improved quantitative approach on semi-solid media. While end-point plaque size is commonly used for quantitative phenotyping, it only provides a static snapshot that fails to disentangle phage performance across changing bacterial growth phases. Here, we report a cost-effective, scalable method to quantify the bacterial killing rates of diverse phages by continuously tracking plaque development using a consumer-grade flatbed scanner and custom open-source computer vision software. By fitting a minimal, five-parameter phenomenological model to the diverse BASEL collection of Escherichia coli phages, we accurately captured plaque expansion dynamics across all bacterial growth phases. We discovered that plaque expansion parameters vary significantly across individual phages and taxonomic families, with the initial expansion rate demonstrating an inverse correlation with viral genome size. Notably, these plaque expansion dynamics did not correlate with bacterial collapse times in liquid culture, indicating that they represent distinct physiological and kinetic processes. Finally, we demonstrated that phage taxonomic family can be predicted with 72.1% accuracy solely from dynamic plaque parameters using a random forest classifier, whereas conventional endpoint plaque sizes held no predictive power. This time-lapse phenotyping approach yields robust kinetic data superior to end-point measurements, establishing a scalable platform to accelerate phage discovery and genotype-phenotype mapping.

The gut-airway axis is increasingly recognised as an important determinant of respiratory health in livestock. However, the interactions among respiratory infections, microbial communities, mucosal barrier function, and microbial metabolism that underpin this axis remain poorly understood. Here, we combined longitudinal immunological, mucosal, microbiological, and metabolic analyses to characterise host-microbiota interactions during antibiotic treatment and experimental bovine respiratory disease (BRD). We found that the bovine respiratory tract is organised as a hierarchically structured mucosal ecosystem, with compartment-specific mucin expression patterns associated with distinct microbial communities. Within this framework, increasing respiratory pathogen loads were associated with greater clinical severity, oxidative stress, acute inflammation, and targeted alterations in mucosal regulation and microbial composition, while serum amyloid A emerged as the host biomarker most closely reflecting infection intensity. Although infection altered host and microbial states, antimicrobial treatment represented the dominant ecological perturbation. Antibiotic exposure reshaped both respiratory and intestinal microbial communities, reducing diversity, increasing community drift, and driving coordinated ecological transitions across the gut-airway axis that exceeded the effects of infection alone. Importantly, these compositional changes were accompanied by profound reductions in acetate, propionate, butyrate, and valerate, identifying short-chain fatty acids as a major functional consequence of gut microbiota disruption. Integrative multi-omics analyses revealed that mucosal organisation, microbial diversity, microbial metabolism, and immune activity were embedded within a highly structured gut-airway network. Together, these findings support a holobiont view of BRD in which disease outcomes emerge not only from pathogen-host interactions but also from ecological disturbances that propagate across interconnected mucosal compartments. Notably, antimicrobial treatment affected ecosystem structure and function more than infection itself, highlighting potential trade-offs between pathogen control and microbiota stability. Preserving microbial resilience and its metabolic functions may therefore represent an important component of future strategies to improve respiratory health and disease resilience in cattle.

Long-term microbial population dynamics can reveal how ocean microbiomes respond to environmental variation, but resolving within-species change requires suitable genomic references. We integrated 30 PacBio HiFi long-read metagenome-assembled genomes (MAGs) from the northwestern Mediterranean Sea with monthly Illumina metagenomes from neighboring time series spanning 12 and 7 years. Population differentiation varied widely among genomes and often differed in magnitude between sites, while a subset showed increasing differentiation with temporal distance. Genetic diversity and population structure were frequently seasonal, with populations associated with warm or cold waters and others spanning both thermal regimes. Candidate adaptive genes accounted for 1.6-3.0% of coding sequences, and many were detected in both time series. Their pN/pS trajectories included regular seasonal fluctuations consistent with changing strain or ecotype contributions, alongside irregular temporal patterns. These results show that combining long-read MAG references with dense short-read time series can resolve seasonal population structure, site-specific differentiation, and candidate adaptive variation within marine microbial species, providing a framework for tracking contemporary population-genomic change in the ocean.