Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Genome, Bacterial”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Standardized phylogenetic and molecular evolutionary analysis applied to species across the microbial tree of life

There is growing interest in reconstructing phylogenies from the copious amounts of genome sequencing projects that target related viral, bacterial or eukaryotic organisms. To facilitate the construction of standardized and robust phylogenies for disparate types of projects, we have developed a complete bioinformatic workflow, with a web-based component to perform phylogenetic and molecular evolutionary (PhaME) analysis from sequencing reads, draft assemblies or completed genomes of closely related organisms. Furthermore, the ability to incorporate raw data, including some metagenomic samples containing a target organism (e.g. from clinical samples with suspected infectious agents), shows promise for the rapid phylogenetic characterization of organisms within complex samples without the need for prior assembly.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Mapping the soil microbiome functions shaping wetland methane emissions

Accounting for only 8% of Earth’s land cover, freshwater wetlands remain the foremost contributors to global methane emissions. Yet the microorganisms and processes underlying methane emissions from wetland soils remain poorly understood. Over a five-year period, we surveyed the microbial membership and in situ methane measurements from over 700 samples in one of the most prolific methane-emitting wetlands in the United States. We constructed a catalog of 2,502 metagenome-assembled genomes (MAGs), with more than half of the 70 bacterial and archaeal phyla sampled containing novel lineages. Integration of these data with 133 soil metatranscriptomes provided a genome-resolved view of the biogeochemical specialization and versatility expressed over wetland soil spatial and temporal gradients. Centimeter-scale depth differences best explained patterns of microbial community structure and transcribed functionalities, even more than land cover or temporal information. Moreover, while extended flooding restructured soil redox, this perturbation failed to reconfigure the transcriptional profiles of methane-cycling microorganisms, contrasting with theoretically expected responses to hydrological perturbations. Co-expression analyses, coupled with depth-resolved methane measurements, revealed the metabolisms and trophic structures most predictive of methane hotspots. Mapping the spatiotemporal transcriptional patterns on this compendium of biogeochemically classified soil-derived genomes begins to untangle the microbial carbon, energy, and nutrient processing contributing to wetland methane production.

MAG↗
DOE OSTI2025

Infection and Genomic Properties of Single- and Double-Stranded DNA Cellulophaga Phages

Bacterial viruses (phages) are abundant and ecologically impactful, but laboratory-based experimental model systems vastly under-represent known phage diversity, particularly for ssDNA phages. Here, we characterize the genomes and infection properties of two unrelated marine flavophages—ssDNA generalist phage phi18:4 (6.5 Kbp) and dsDNA specialist phage phi18:1 (39.2 Kbp)—when infecting the same Cellulophaga baltica strain #18 (Cba18), of the class Flavobacteriia. Phage phi18:4 belongs to a new family of ssDNA phages, has an internal lipid membrane, and its genome encodes primarily structural proteins, as well as a DNA replication protein common to ssDNA phages and a unique lysis protein. Phage phi18:1 is a siphovirus that encodes several virulence genes, despite not having a known temperate lifestyle, a CAZy enzyme likely for regulatory purposes, and four DNA methyltransferases dispersed throughout the genome that suggest both host modulation and phage DNA protection against host restriction. Physiologically, ssDNA phage phi18:4 has a shorter latent period and smaller burst size than dsDNA phage phi18:1, and both phages efficiently infect this host. These results help augment the diversity of characterized environmental phage–host model systems by studying infections of genomically diverse phages (ssDNA vs. dsDNA) on the same host.

Howard-Varona, Cristina (ORCID:0000000241495818)↗
DOE OSTI2021

Searching for fat tails in CRISPR-Cas systems: Data analysis and mathematical modeling

Understanding CRISPR-Cas systems—the adaptive defence mechanism that about half of bacterial species and most of archaea use to neutralise viral attacks—is important for explaining the biodiversity observed in the microbial world as well as for editing animal and plant genomes effectively. The CRISPR-Cas system learns from previous viral infections and integrates small pieces from phage genomes called spacers into the microbial genome. The resulting library of spacers collected in CRISPR arrays is then compared with the DNA of potential invaders. One of the most intriguing and least well understood questions about CRISPR-Cas systems is the distribution of spacers across the microbial population. Here, using empirical data, we show that the global distribution of spacer numbers in CRISPR arrays across multiple biomes worldwide typically exhibits scale-invariant power law behaviour, and the standard deviation is greater than the sample mean. We develop a mathematical model of spacer loss and acquisition dynamics which fits observed data from almost four thousand metagenomes well. In analogy to the classical ‘rich-get-richer’ mechanism of power law emergence, the rate of spacer acquisition is proportional to the CRISPR array size, which allows a small proportion of CRISPRs within the population to possess a significant number of spacers. Our study provides an alternative explanation for the rarity of all-resistant super microbes in nature and why proliferation of phages can be highly successful despite the effectiveness of CRISPR-Cas systems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Increasing Bacterial Tolerance and Metabolism of the Biofuel, N-Butanol Using Community-Level Evolution and Functional Genomics

Bioremediation capability should be developed along with biofuel technology to mitigate the potential damage of future spills. One biofuel that is being developed is biobutanol, since n-butanol is more energy dense and less volatile than ethanol. A bottleneck for industrial production of biobutanol is its toxicity; most microbes cannot survive about 1.5% v/v. Thus, microbial bioremediation of n-butanol would need microbes that can both tolerate and metabolize butanol. We used ecological and evolutionary biology approaches to find bacteria that could tolerate and metabolize butanol. We then tried to increase the metabolism of butanol by promising bacterial strains and communities. Ecological community-level assays and screenings were conducted followed by 16S amplicon sequencing to identify butanol-tolerant artificial bacterial communities. Promising communities were then tested for growth with butanol as the sole carbon source. Secondly, we looked for bacteria with alcohol dehydrogenase enzymes and looked to increase butanol metabolism. We found that the tolerance for n-butanol may be improved with repeated exposure, but it was difficult to switch from tolerance to metabolism. Bacterial community v dynamics may be influenced by n-butanol concentration, and there was putative butanol metabolism found with both research approaches.

09 BIOMASS FUELS↗
DOE OSTI2026

Isolation of genome-predicted Caldatribacterium ( Atribacterota ) reveals pervasive microbial cultivation problem due to folate precipitation

Most bacterial phyla have few or no pure cultures, including Atribacterota , comprised of ubiquitous anaerobes. Here, we report genome-guided enrichment and isolation of two Atribacterota species representing a new family, Caldatribacterium saccharofermentans from a hot spring, and Caldatribacterium inferamans from a deep aquifer. Both were co-enriched with sulfate-reducing bacteria and initially resisted isolation, which we link to inadvertent removal of precipitated folic acid by filter-sterilization of unbuffered Wolin’s vitamin solution. We then predict folate auxotrophy across the Atribacterota and ~29% of all bacteria, with extensive auxotrophy in 27% of phyla. Since ≥604 of 791 ( ≥ 76%) media with folic acid additions in the MediaDive database use unbuffered vitamin solutions in which folic acid is likely removed during filter-sterilization, we propose that folate auxotrophy limits culturability in defined media en masse. We also uncover unusual features of Caldatribacterium , including three lipid membrane-like layers (LMLs), with the inner LML surrounding the nucleoid, and a high percentage of secreted proteins, supporting a unique cell biology of Atribacterota .

Biological and medical sciences↗
DOE OSTI2021

Fungal-Bacterial Cooccurrence Patterns Differ between Arbuscular Mycorrhizal Fungi and Nonmycorrhizal Fungi across Soil Niches

Soil bacteria and fungi are known to form niche-specific communities that differ between actively growing and decaying roots. Yet almost nothing is known about the cross-kingdom interactions that frame these communities and the environmental filtering that defines these potentially friendly or competing neighbors. We explored the temporal and spatial patterns of soil fungal (mycorrhizal and nonmycorrhizal) and bacterial cooccurrence near roots of wild oat grass, Avena fatua, growing in its naturalized soil in a greenhouse experiment. Amplicon sequences of the fungal internal transcribed spacer (ITS) and bacterial 16S rRNA genes from rhizosphere and bulk soils collected at multiple plant growth stages were used to construct covariation-based networks as a step toward identifying fungal-bacterial associations. Corresponding stable-isotope-enabled metagenome-assembled genomes (MAGs) of bacteria identified in cooccurrence networks were used to inform potential mechanisms underlying the observed links. Bacterial-fungal networks were significantly different in rhizosphere versus bulk soils and between arbuscular mycorrhizal fungi (AMF) and nonmycorrhizal fungi. Over 12 weeks of plant growth, nonmycorrhizal fungi formed increasingly complex networks with bacteria in rhizosphere soils, while AMF more frequently formed networks with bacteria in bulk soils. Analysis of network-associated bacterial MAGs suggests that some of the fungal-bacterial links that we identified are potential indicators of bacterial breakdown and consumption of fungal biomass, while others intimate shared ecological niches.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Genome-Resolved Metaproteomics Decodes the Microbial and Viral Contributions to Coupled Carbon and Nitrogen Cycling in River Sediments

Rivers have a significant role in global carbon and nitrogen cycles, serving as a nexus for nutrient transport between terrestrial and marine ecosystems. Although rivers have a small global surface area, they contribute substantially to worldwide greenhouse gas emissions through microbially mediated processes within the river hyporheic zone. Despite this importance, research linking microbial and viral communities to specific biogeochemical reactions is still nascent in these sediment environments. To survey the metabolic potential and gene expression underpinning carbon and nitrogen biogeochemical cycling in river sediments, we collected an integrated data set of 33 metagenomes, metaproteomes, and paired metabolomes. We reconstructed over 500 microbial metagenome-assembled genomes (MAGs), which we dereplicated into 55 unique, nearly complete medium- and high-quality MAGs spanning 12 bacterial and archaeal phyla. We also reconstructed 2,482 viral genomic contigs, which were dereplicated into 111 viral MAGs (vMAGs) of >10 kb in size. As a result of integrating gene expression data with geochemical and metabolite data, we created a conceptual model that uncovered new roles for microorganisms in organic matter decomposition, carbon sequestration, nitrogen mineralization, nitrification, and denitrification. We show how these metabolic pathways, integrated through shared resource pools of ammonium, carbon dioxide, and inorganic nitrogen, could ultimately contribute to carbon dioxide and nitrous oxide fluxes from hyporheic sediments. Further, by linking viral MAGs to these active microbial hosts, we provide some of the first insights into viral modulation of river sediment carbon and nitrogen cycling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

CRISPR-Cas “Non-Target” Sites Inhibit On-Target Cutting Rates

CRISPR-Cas systems have become ubiquitous for genome editing in eukaryotic as well as bacterial systems. Cas9 forms a complex with a guide RNA (gRNA) and searches DNA for a matching sequence (target site) next to a protospacer adjacent motif (PAM). Once found, Cas9 cuts the DNA. Cas9 is revolutionary for the ability to change the RNA sequence and target a new site easily. However, while algorithms have been developed to predict gRNA-specific Cas9 activity, a fundamental biological understanding of gRNA-specific activity is lacking. Here, the number of PAM sites in the genome is effectively a large pool of inhibitory substrates, competing with the target site for the Cas9/gRNA complex. We demonstrate that increasing the number of non-target sites for a given gRNA reduces on-target activity in a dose-dependent manner. Furthermore, we show that the use of Cas9 mutants with increased PAM specificity toward a smaller subset of PAMs (or smaller pool of competitive substrates) improves cutting rates, while increased PAM promiscuity decreases cutting rates. Decreasing the potential search space by increasing PAM specificity provides a path toward improving on-target activity for slower high-fidelity Cas9 variants. Engineering improved PAM specificity to reduce the competitive search space offers an alternative strategy to engineer Cas9 variants with increased specificity and maintained on-target activity.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Bacterial Origin and Reductive Evolution of the CPR Group

The candidate phyla radiation (CPR) is a proposed subdivision within the bacterial domain comprising several candidate phyla. CPR organisms are united by small genome and physical sizes, lack several metabolic enzymes, and populate deep branches within the bacterial subtree of life. Here, these features raise intriguing questions regarding their origin and mode of evolution. In this study, we performed a comparative and phylogenomic analysis to investigate CPR origin and evolution. Unlike previous gene/protein sequence-based reports of CPR evolution, we used protein domain superfamilies classified by protein structure databases to resolve the evolutionary relationships of CPR with non-CPR bacteria, Archaea, Eukarya, and viruses. Across all supergroups, CPR shared maximum superfamilies with non-CPR bacteria and were placed as deep branching bacteria in most phylogenomic trees. CPR contributed 1.22% of new superfamilies to bacteria including the ribosomal protein L19e and encoded four core superfamilies that are likely involved in cell-to-cell interaction and establishing episymbiotic lifestyles. Although CPR and non-CPR bacterial proteomes gained common superfamilies over the course of evolution, CPR and Archaea had more common losses. These losses mostly involved metabolic superfamilies. In fact, phylogenies built from only metabolic protein superfamilies separated CPR and non-CPR bacteria. These findings indicate that CPR are bacterial organisms that have probably evolved in an Archaea-like manner via the early loss of metabolic functions. We also discovered that phylogenies built from metabolic and informational superfamilies gave contrasting views of the groupings among Archaea, Bacteria, and Eukarya, which add to the current debate on the evolutionary relationships among superkingdoms.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

The Genetics of Aerotolerant Growth in Zymomonas mobilis

Reduced genome bacteria are genetically simplified systems that facilitate biological study and industrial use. The free-living Alphaproteobacterium, Zymomonas mobilis, has a naturally reduced genome containing fewer than 2000 protein coding genes. Despite its small genome, Z. mobilis thrives in diverse conditions including the presence or absence of atmospheric oxygen. However, insufficient characterization of essential and conditionally essential genes has limited broader adoption of Z. mobilis as a model Alphaproteobacterium. Here, we use genome-scale CRISPRi-seq to systematically identify and characterize Z. mobilis genes that are conditionally essential for aerotolerant or anaerobic growth, or are generally essential across both conditions. Comparative genomics revealed that the essentiality of most "generally essential" genes was shared between Z. mobilis and other Alphaproteobacteria, validating Z. mobilis as reduced genome model. Among conditionally essential genes, we found that the DNA repair gene, recJ, was critical only for aerobic growth but reduced the mutation rate under both conditions. Further, we show that genes encoding the F1FO ATP synthase and Rnf respiratory complex are required for anaerobic growth of Z. mobilis. Combining CRISPRi partial knockdowns with metabolomics and membrane potential measurements, we determined that the ATP synthase generates membrane potential that is consumed by Rnf to power downstream processes. Rnf knockdown strains accumulated isoprenoid biosynthesis intermediates, suggesting a key role for Rnf in powering essential biosynthetic reactions. Our work establishes Z. mobilis as a streamlined model for alphaproteobacterial genetics, has broad implications in bacterial energy coupling, and informs Z. mobilis genome manipulation for optimized production of valuable isoprenoid-based bioproducts

Alphaproteobacteria↗
DOE OSTI2023

The genetics of aerotolerant growth in an alphaproteobacterium with a naturally reduced genome

Reduced genome bacteria are genetically simplified systems that facilitate biological study and industrial use. The free-living alphaproteobacterium Zymomonas mobilis has a naturally reduced genome containing fewer than 2,000 protein-coding genes. Despite its small genome, Z. mobilis thrives in diverse conditions including the presence or absence of atmospheric oxygen. However, insufficient characterization of essential and conditionally essential genes has limited broader adoption of Z. mobilis as a model alphaproteobacterium. Here, we use genome-scale CRISPRi-seq (clustered regularly interspaced short palindromic repeats interference sequencing) to systematically identify and characterize Z. mobilis genes that are conditionally essential for aerotolerant or anaerobic growth or are generally essential across both conditions. Comparative genomics revealed that the essentiality of most “generally essential” genes was shared between Z. mobilis and other Alphaproteobacteria, validating Z. mobilis as a reduced genome model. Among conditionally essential genes, we found that the DNA repair gene, recJ, was critical only for aerobic growth but reduced the mutation rate under both conditions. Further, we show that genes encoding the F1FO ATP synthase and Rhodobacter nitrogen fixation (Rnf) respiratory complex are required for the anaerobic growth of Z. mobilis. Combining CRISPRi partial knockdowns with metabolomics and membrane potential measurements, we determined that the ATP synthase generates membrane potential that is consumed by Rnf to power downstream processes. Rnf knockdown strains accumulated isoprenoid biosynthesis intermediates, suggesting a key role for Rnf in powering essential biosynthetic reactions. Our work establishes Z. mobilis as a streamlined model for alphaproteobacterial genetics, has broad implications in bacterial energy coupling, and informs Z. mobilis genome manipulation for optimized production of valuable isoprenoid-based bioproducts.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Comparative genomics of the Liberibacter genus reveals widespread diversity in genomic content and positive selection history

‘Candidatus Liberibacter’ is a group of bacterial species that are obligate intracellular plant pathogens and cause Huanglongbing disease of citrus trees and Zebra Chip in potatoes. Here, we examined the extent of intra- and interspecific genetic diversity across the genus using comparative genomics. Our approach examined a wide set of Liberibacter genome sequences including five pathogenic species and one species not known to cause disease. By performing comparative genomics analyses, we sought to understand the evolutionary history of this genus and to identify genes or genome regions that may affect pathogenicity. With a set of 52 genomes, we performed comparative genomics, measured genome rearrangement, and completed statistical tests of positive selection. We explored markers of genetic diversity across the genus, such as average nucleotide identity across the whole genome. These analyses revealed the highest intraspecific diversity amongst the ‘Ca. Liberibacter solanacearum’ species, which also has the largest plant host range. We identified sets of core and accessory genes across the genus and within each species and measured the ratio of nonsynonymous to synonymous mutations (dN/dS) across genes. We identified ten genes with evidence of a history of positive selection in the Liberibacter genus, including genes in the Tad complex, which have been previously implicated as being highly divergent in the ‘Ca. L. capsica’ species based on high values of dN.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Bleach Rescues Nannochloropsis from an Obligate Parasite and Alters Microbial and Metabolite Signatures of Outdoor Cultures

Chemical agents are commonly used to protect algal crops. Yet, few studies have characterized the effects of these agents on associated microbial communities to understand effects on microbial functions relevant to algal crop production and protection. Here, we used shotgun metagenomic sequencing and untargeted exometabolite profiling to link the application of bleach, a -cidal agent used to protect algae from pests, to changes in community composition, metabolic pathways, and exometabolies - at a whole community level. Bleach protected the algal crop from crashing but altered bacterial diversity. Analysis of metagenome-assembled genomes (MAGs) revealed a classic predator-prey cycle between Oligoflexus and our target alga Nannochloropsis. Olifoflexus genomes from our study were notably similar to a previously identified BALO (Bdellovibrio and like organism), FD111, known to kill Nannochloropsis cultures, providing strong evidence that an FD111-like organism was responsible for the crash. Metabolic pathway composition differed between bleached and unbleached ponds, with abundance of twelve pathways related to stress tolerance, including the superpathway of methylglyoxal degradation, lipid IVA biosynthesis, and ectoine biosynthesis, greater in bleached ponds compared to unbleached ponds. Virulence factors related to adherence, biofilm formation, motility, and pathogenicity increased dramatically in bleached ponds with time, although this increase was not coupled with an increase in pathogens - algal or otherwise - or a decline in algal health. Our study highlights the importance of coupling 16S rRNA gene sequencing with whole genome data and other -omics tools to sketch a larger picture of community structure and function in crop systems. Moreover, our results highlight that continued long-term bleaching may lead to negative effects to crop health or downstream adverse health effects to humans or animals, depending on the algal product (i.e. human supplements or animal feedstocks). Future work on alternative treatment methods that would reduce resistance is necessary in the field.

09 BIOMASS FUELS↗
DOE OSTI2024

Transport of metformin metabolites by guanidinium exporters of the small multidrug resistance family

Proteins from the small multidrug resistance (SMR) family are frequently associated with horizontally transferred multidrug resistance gene arrays found in bacteria from wastewater and the human-adjacent biosphere. Recent studies suggest that a subset of SMR transporters might participate in the metabolism of the common pharmaceutical metformin by bacterial consortia. Here, we show that both genomic and plasmid-associated transporters of the SMRGdx functional subtype export byproducts of microbial metformin metabolism, with particularly high export efficiency for guanylurea. We use solid-supported membrane electrophysiology to evaluate the transport kinetics for guanylurea and native substrate guanidinium by four representative SMRGdx homologs. Using an internal reference to normalize independent electrophysiology experiments, we show that transport rates are comparable for genomic and plasmid-associated SMRGdx homologs, and using a proteoliposome-based transport assay, we show that 2 proton:1 substrate transport stoichiometry is maintained. Additional characterization of guanidinium and guanylurea export properties focuses on the structurally characterized homolog, Gdx-Clo, for which we examined the pH dependence and thermodynamics of substrate binding and solved an x-ray crystal structure with guanylurea bound. Together, these experiments contribute in two main ways. By providing the first detailed kinetic examination of the structurally characterized SMRGdx homolog Gdx-Clo, they provide a functional framework that will inform future mechanistic studies of this model transport protein. Second, this study casts light on a potential role for SMRGdx transporters in microbial handling of metformin and its microbial metabolic byproducts, providing insight into how native transport physiologies are co-opted to contend with new selective pressures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2019

Methylation Pattern Detection in the Genome of Bacillus Pumilus Strain SAFR-032

Bacillus pumilus SAFR-032, an endospore-forming bacterial strain that was isolated from a spacecraft assembly facility (SAFR), was investigated to determine its methylation pattern (methylome) across the genome in comparison to the previously sequenced reference genome. In addition, a version of SAFR-032 that was flown as spores for 18 months on the International Space Station (ISS) was also investigated for possible genomic changes due to long-duration ISS-flight and to determine if methylation patterns may have changed. Both the genomics and methylomics were conducted using a Nanopore MinION sequencing device. In addition to the omics investigation, the two SAFR-032 strains, ISS flown and non-ISS flown, were compared phenotypically in chamber experiments testing individual environmental insults: ionizing radiation, UV exposure, and cold desiccation (i.e. freeze drying). Results from this study inform on Planetary Protection concerns and will reveal potential DNA damage associated with long-term spaceflight and how such damage may influence survivors after being transported to an extraterrestrial environment, such as Mars.

Serda, Bianca M.↗
DOE OSTI2021

Comparative Virulence and Genomic Analysis of Streptococcus suis Isolates

Streptococcus suis is a zoonotic bacterial swine pathogen causing substantial economic and health burdens to the pork industry. Mechanisms used by S. suis to colonize and cause disease remain unknown and vaccines and/or intervention strategies currently do not exist. Studies addressing virulence mechanisms used by S. suis have been complicated because different isolates can cause a spectrum of disease outcomes ranging from lethal systemic disease to asymptomatic carriage. The objectives of this study were to evaluate the virulence capacity of nine United States S. suis isolates following intranasal challenge in swine and then perform comparative genomic analyses to identify genomic attributes associated with swine-virulent phenotypes. No correlation was found between the capacity to cause disease in swine and the functional characteristics of genome size, serotype, sequence type (ST), or in vitro virulence-associated phenotypes. A search for orthologs found in highly virulent isolates and not found in non-virulent isolates revealed numerous predicted protein coding sequences specific to each category. While none of these predicted protein coding sequences have been previously characterized as potential virulence factors, this analysis does provide a reliable one-to-one assignment of specific genes of interest that could prove useful in future allelic replacement and/or functional genomic studies. Collectively, this report provides a framework for future allelic replacement and/or functional genomic studies investigating genetic characteristics underlying the spectrum of disease outcomes caused by S. suis isolates.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Dynamics of phage-host interactions in Bacteroides fragilis resolved by single-cell transcriptomics

The interactions between lytic phages and their hosts are typically studied in bulk culture, which obscures cell-cell differences in infection susceptibility or expression of protective factors. Here, we use bacterial single-cell RNA sequencing to profile the transcriptomes of ~50,000 cells from cultures of a human pathobiont, Bacteroides fragilis, infected with a lytic bacteriophage. From a single sampling, we quantified the asynchronous progression of phage infection in individual bacterial cells and reconstructed the infection timeline, characterizing both host and phage transcriptomic changes as infection unfolded. Further, we discovered phenotypic subpopulations of bacteria that remained uninfected. Each cell’s vulnerability to phage infection was influenced by expression of multiple genetic loci, most prominently phase-variable capsular polysaccharide (CPS) biosynthesis pathways and an operon predicted to encode fimbrial genes. These findings uncovered genome-wide phase variation and stochasticity that enable bacterial survival and re-growth without acquiring additional mutations. Overall, we establish bacterial single-cell RNA sequencing as a powerful platform for investigating the dynamics of host-phage interactions and revealing the roles of phase variation and stochasticity in bacterial defenses.

Bacteria↗