Search NASA⌕ Search

SEARCH · Search NASA

Results for “DNA, 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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Microbial Metagenomes Across a Complete Phytoplankton Bloom Cycle: High-Resolution Sampling Every 4 Hours Over 22 Days

In May and June of 2021, marine microbial samples were collected for DNA sequencing in East Sound, WA, USA every 4 hours for 22 days. This high temporal resolution sampling effort captured the last 3 days of a Rhizosolenia sp. bloom, the initiation and complete bloom cycle of Chaetoceros socialis (8 days), and the following bacterial bloom (2 days). Metagenomes were completed on the time series, and the dataset includes 128 size-fractionated microbial samples (0.22–1.2 µm), providing gene abundances for the dominant members of bacteria, archaea, and viruses. This dataset also has time-matched nutrient analyses, flow cytometry data, and physical parameters of the environment at a single point of sampling within a coastal ecosystem that experiences regular bloom events, facilitating a range of modeling efforts that can be leveraged to understand microbial community structure and their influences on the growth, maintenance, and senescence of phytoplankton blooms.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

DL-TODA: A Deep Learning Tool for Omics Data Analysis

Metagenomics is a technique for genome-wide profiling of microbiomes; this technique generates billions of DNA sequences called reads. Given the multiplication of metagenomic projects, computational tools are necessary to enable the efficient and accurate classification of metagenomic reads without needing to construct a reference database. The program DL-TODA presented here aims to classify metagenomic reads using a deep learning model trained on over 3000 bacterial species. A convolutional neural network architecture originally designed for computer vision was applied for the modeling of species-specific features. Using synthetic testing data simulated with 2454 genomes from 639 species, DL-TODA was shown to classify nearly 75% of the reads with high confidence. The classification accuracy of DL-TODA was over 0.98 at taxonomic ranks above the genus level, making it comparable with Kraken2 and Centrifuge, two state-of-the-art taxonomic classification tools. DL-TODA also achieved an accuracy of 0.97 at the species level, which is higher than 0.93 by Kraken2 and 0.85 by Centrifuge on the same test set. Application of DL-TODA to the human oral and cropland soil metagenomes further demonstrated its use in analyzing microbiomes from diverse environments. Compared to Centrifuge and Kraken2, DL-TODA predicted distinct relative abundance rankings and is less biased toward a single taxon.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

High-throughput, single-microbe genomics with strain resolution, applied to a human gut microbiome

We present Microbe-seq, a high-throughput single-microbe method that yields strain-resolved genomes from complex microbial communities. We encapsulate individual microbes into droplets with microfluidics and liberate their DNA, which we amplify, tag with droplet-specific barcodes, and sequence. We use Microbe-seq to explore the human gut microbiome; we collect stool samples from a single individual, sequence over 20,000 microbes, and reconstruct nearly-complete genomes of almost 100 bacterial species, including several with multiple subspecies strains. We use these genomes to probe genomic signatures of microbial interactions: we reconstruct the horizontal gene transfer (HGT) network within the individual and observe far greater exchange within the same bacterial phylum than between different phyla. We probe bacteria-virus interactions; unexpectedly, we identify a significant in vivo association between crAssphage, an abundant bacteriophage, and a single strain of Bacteroides vulgatus. Microbe-seq contributes high-throughput culture-free capabilities to investigate genomic blueprints of complex microbial communities with single-microbe resolution.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

XA21-mediated resistance to Xanthomonas oryzae pv. oryzae is dose dependent

The rice receptor kinase XA21 confers broad-spectrum resistance to Xanthomonas oryzae pv. oryzae ( Xoo ), the causal agent of rice bacterial blight disease. To investigate the relationship between the expression level of XA21 and resulting resistance, we generated independent HA-XA21 transgenic rice lines accumulating the XA21 immune receptor fused with an HA epitope tag. Whole-genome sequence analysis identified the T-DNA insertion sites in sixteen independent T0 events. Further, through quantification of the HA-XA21 protein and assessment of the resistance to Xoo strain PXO99 in six independent transgenic lines, we observed that XA21-mediated resistance is dose dependent. In contrast, based on the four agronomic traits quantified in these experiments, yield is unlikely to be affected by the expression level of HA-XA21 . These findings extend our knowledge of XA21-mediated defense and contribute to the growing number of well-defined genomic landing pads in the rice genome that can be targeted for gene insertion without compromising yield.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Interrelationships among methods of estimating microbial biomass across multiple soil orders and biomes: Supporting data

This dataset contains environmental and soil measurements from 18 different locations across the globe including the SPRUCE experiment site and multiple sampling depths, with 17 of these locations having samples processed between 2012-2013 and one location (SPRUCE) collected in 2021 and processed in 2022. Environmental measurements include: mean annual temperature, mean annual precipitation, and 30-day presampling temperature. Soil physicochemical measurements include: particle size analysis (PSA), pH, gravimetric moisture content (GMC), bulk soil carbon (C) and nitrogen (N), total organic C and N, C:N ratio, and dissolved organic carbon (DOC). Soil biological measurements include: microbial biomass carbon (MBC) measured through chloroform fumigation extraction (CFE), gene copy numbers (GCN) of bacteria, fungi, and archaea measured through quantitative polymerase chain reaction (qPCR), DNA yield measured through Nanodrop spectrophotometry, and phospholipid fatty acids (PLFA) of bacteria and fungi measured through PLFA analysis. This data set contains one file in comma separate (*.csv) format.

archaea gene copy number↗
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 OSTI2020

Improved bacterial recombineering by parallelized protein discovery

Exploiting bacteriophage-derived homologous recombination processes has offered precise, multiplex editing of microbial genomes and the construction of billions of customized genetic variants in a single day. The techniques that enable this, multiplex automated genome engineering (MAGE) and directed evolution with random genomic mutations (DIvERGE), are however, currently limited to a handful of microorganisms for which single-stranded DNA-annealing proteins (SSAPs) that promote efficient recombineering have been identified. Thus, to enable genome-scale engineering in new hosts, efficient SSAPs must first be found. Here we present a high-throughput method for SSAP discovery that we call “serial enrichment for efficient recombineering” (SEER). By performing SEER in Escherichia coli to screen hundreds of putative SSAPs, we identify highly active variants PapRecT and CspRecT. CspRecT increases the efficiency of single-locus editing to as high as 50% and improves multiplex editing by 5- to 10-fold in E. coli , while PapRecT enables efficient recombineering in Pseudomonas aeruginosa , a concerning human pathogen. CspRecT and PapRecT are also active in other, clinically and biotechnologically relevant enterobacteria. We envision that the deployment of SEER in new species will pave the way toward pooled interrogation of genotype-to-phenotype relationships in previously intractable bacteria.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Longitudinal Multi-omics Reveal Phase-Dependent Viral Adaptive Strategies and Functional Potential During Formation of Algal-bacterial Granular Sludge

Virus-host interactions within microbial aggregates critically influence microbiome function and stability, yet how physicochemical stresses shape the interactive dynamics remains largely unexplored. Here, we investigated virus–host dynamics during the transition of algal-bacterial granular sludge (ABGS) from activated sludge under continuous hydraulic shear using integrated metagenomics and metatranscriptomics. Hydraulic stress initially reduced host a-diversity, which coincided with a marked increase in viral lysogenicity. During this host diversity bottleneck, viral microdiversity increased, and genes related to virion structure and DNA packaging were under positive selection (pN/pS >1). As host diversity recovered, viral microdiversity declined, while viral anti-defense systems (ADS) significantly increased in abundance. Lagged correlation analysis revealed a significant positive correlation between viral ADS and host defense systems (DS), suggesting an evolutionary arms race. Furthermore, active lysogenic infections were accompanied by enrichment of DS and auxiliary viral genes (AVGs) involved in genetic information processing and amino acid metabolism, potentially enhancing host fitness. Overall, our study unveils a phase-dependent co-evolutionary interplay between viruses and hosts during ABGS formation, providing insights into the development and maintenance of microbial structural and functional resilience in engineered ecosystems.

Qi, Huiyuan↗
DOE OSTI2024

Discovery of GuaB inhibitors with efficacy against Acinetobacter baumannii infection

ABSTRACT Guanine nucleotides are required for growth and viability of cells due to their structural role in DNA and RNA, and their regulatory roles in translation, signal transduction, and cell division. The natural antibiotic mycophenolic acid (MPA) targets the rate-limiting step inde novoguanine nucleotide biosynthesis executed by inosine-5´-monophosphate dehydrogenase (IMPDH). MPA is used clinically as an immunosuppressant, but whetherin vivoinhibition of bacterial IMPDH (GuaB) is a valid antibacterial strategy is controversial. Here, we describe the discovery of extremely potent small molecule GuaB inhibitors (GuaBi) specific to pathogenic bacteria with a low frequency of on-target spontaneous resistance and bactericidal efficacyin vivoagainstAcinetobacter baumanniimouse models of infection. The spectrum of GuaBi activity includes multidrug-resistant pathogens that are a critical priority of new antibiotic development. Co-crystal structures ofA. baumannii, Staphylococcus aureus, andEscherichia coliGuaB proteins bound to inhibitors show comparable binding modes of GuaBi across species and identifies key binding site residues that are predictive of whole-cell activity across both Gram-positive and Gram-negative clades of Bacteria. The clearin vivoefficacy of these small molecule GuaB inhibitors in a model ofA. baumanniiinfection validates GuaB as an essential antibiotic target. IMPORTANCE The emergence of multidrug-resistant bacteria worldwide has renewed interest in discovering antibiotics with novel mechanism of action. For the first time ever, we demonstrate that pharmacological inhibition ofde novoguanine biosynthesis is bactericidal in a mouse model ofAcinetobacter baumanniiinfection. Structural analyses of novel inhibitors explain differences in biochemical and whole-cell activity across bacterial clades and underscore why this discovery may have broad translational impact on treatment of the most recalcitrant bacterial infections.

Microbiology↗
DOE OSTI2021

In Vivo Entombment of Bacteria and Fungi during Calcium Oxalate, Brushite, and Struvite Urolithiasis

Human kidney stones form via repeated events of mineral precipitation, partial dissolution, and reprecipitation, which are directly analogous to similar processes in other natural and manmade environments, where resident microbiomes strongly influence biomineralization. High-resolution microscopy and high-fidelity metagenomic (microscopy-to-omics) analyses, applicable to all forms of biomineralization, have been applied to assemble definitive evidence of in vivo microbiome entombment during urolithiasis. Stone fragments were collected from a randomly chosen cohort of 20 patients using standard percutaneous nephrolithotomy (PCNL). Fourier transform infrared (FTIR) spectroscopy indicated that 18 of these patients were calcium oxalate (CaOx) stone formers, whereas one patient formed each formed brushite and struvite stones. This apportionment is consistent with global stone mineralogy distributions. Stone fragments from seven of these 20 patients (five CaOx, one brushite, and one struvite) were thin sectioned and analyzed using brightfield (BF), polarization (POL), confocal, super-resolution autofluorescence (SRAF), and Raman techniques. DNA from remaining fragments, grouped according to each of the 20 patients, were analyzed with amplicon sequencing of 16S rRNA gene sequences (V1–V3, V3–V5) and internal transcribed spacer (ITS1, ITS2) regions. Bulk-entombed DNA was sequenced from stone fragments in 11 of the 18 patients who formed CaOx stones, and the patients who formed brushite and struvite stones. These analyses confirmed the presence of an entombed low-diversity community of bacteria and fungi, including Actinobacteria, Bacteroidetes, Firmicutes, Proteobacteria, and Aspergillus niger. Bacterial cells approximately 1 μm in diameter were also optically observed to be entombed and well preserved in amorphous hydroxyapatite spherules and fans of needle-like crystals of brushite and struvite. In conclusion, these results indicate a microbiome is entombed during in vivo CaOx stone formation. Similar processes are implied for brushite and struvite stones. This evidence lays the groundwork for future in vitro and in vivo experimentation to determine how the microbiome may actively and/or passively influence kidney stone biomineralization.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Mechanistic insights from structure of Mycobacterium smegmatis topoisomerase I with ssDNA bound to both N- and C-terminal domains

Abstract Type IA topoisomerases interact with G-strand and T-strand ssDNA to regulate DNA topology. However, simultaneous binding of two ssDNA segments to a type IA topoisomerase has not been observed previously. We report here the crystal structure of a type IA topoisomerase with ssDNA segments bound in opposite polarity to the N- and C-terminal domains. Titration of small ssDNA oligonucleotides to Mycobacterium smegmatis topoisomerase I with progressive C-terminal deletions showed that the C-terminal region has higher affinity for ssDNA than the N-terminal active site. This allows the C-terminal domains to capture one strand of underwound negatively supercoiled DNA substrate first and position the N-terminal domains to bind and cleave the opposite strand in the relaxation reaction. Efficiency of negative supercoiling relaxation increases with the number of domains that bind ssDNA primarily with conserved aromatic residues and possibly with assistance from polar/basic residues. A comparison of bacterial topoisomerase I structures showed that a conserved transesterification unit (N-terminal toroid structure) for cutting and rejoining of a ssDNA strand can be combined with two different types of C-terminal ssDNA binding domains to form diverse bacterial topoisomerase I enzymes that are highly efficient in their physiological role of preventing excess negative supercoiling in the genome.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Automated model-predictive design of synthetic promoters to control transcriptional profiles in bacteria

Abstract Transcription rates are regulated by the interactions between RNA polymerase, sigma factor, and promoter DNA sequences in bacteria. However, it remains unclear how non-canonical sequence motifs collectively control transcription rates. Here, we combine massively parallel assays, biophysics, and machine learning to develop a 346-parameter model that predicts site-specific transcription initiation rates for any σ 70 promoter sequence, validated across 22132 bacterial promoters with diverse sequences. We apply the model to predict genetic context effects, design σ 70 promoters with desired transcription rates, and identify undesired promoters inside engineered genetic systems. The model provides a biophysical basis for understanding gene regulation in natural genetic systems and precise transcriptional control for engineering synthetic genetic systems.

42 ENGINEERING↗
DOE OSTI2023

Structural basis of Gabija anti-phage defence and viral immune evasion

Bacteria encode hundreds of diverse defence systems that protect them from viral infection and inhibit phage propagation. Gabija is one of the most prevalent anti-phage defence systems, occurring in more than 15% of all sequenced bacterial and archaeal genomes, but the molecular basis of how Gabija defends cells from viral infection remains poorly understood. Here we use X-ray crystallography and cryo-electron microscopy (cryo-EM) to define how Gabija proteins assemble into a supramolecular complex of around 500 kDa that degrades phage DNA. Gabija protein A (GajA) is a DNA endonuclease that tetramerizes to form the core of the anti-phage defence complex. Two sets of Gabija protein B (GajB) dimers dock at opposite sides of the complex and create a 4:4 GajA–GajB assembly (hereafter, GajAB) that is essential for phage resistance in vivo. We show that a phage-encoded protein, Gabija anti-defence 1 (Gad1), directly binds to the Gabija GajAB complex and inactivates defence. A cryo-EM structure of the virally inhibited state shows that Gad1 forms an octameric web that encases the GajAB complex and inhibits DNA recognition and cleavage. Our results reveal the structural basis of assembly of the Gabija anti-phage defence complex and define a unique mechanism of viral immune evasion.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Bovine Rectoanal Junction In Vitro Organ Culture Model System to Study Shiga Toxin-Producing Escherichia coli Adherence

Studies evaluating the interactions between Shiga toxin-producing Escherichia coli O157:H7 (O157) and the bovine recto–anal junction (RAJ) have been limited to either in vitro analyses of bacteria, cells, or nucleic acids at the RAJ, providing limited information. Alternatively, expensive in vivo studies in animals have been conducted. Therefore, our objective was to develop a comprehensive in vitro organ culture system of the RAJ (RAJ-IVOC) that accurately represents all cell types present in the RAJ. This system would enable studies that yield results similar to those observed in vivo. Pieces of RAJ tissue, obtained from unrelated cattle necropsies, were assembled and subjected to various tests in order to determine the optimal conditions for assaying bacterial adherence in a viable IVOC. O157 strain EDL933 and E. coli K12 with known adherence differences were used to standardize the RAJ-IVOC adherence assay. Tissue integrity was assessed using cell viability, structural cell markers, and histopathology, while the adherence of bacteria was evaluated via microscopy and culture methods. DNA fingerprinting verified the recovered bacteria against the inoculum. When the RAJ-IVOC was assembled in Dulbecco’s Modified Eagle Medium, maintained at a temperature of 39 °C with 5% CO 2 and gentle shaking for a duration of 3–4 h, it successfully preserved tissue integrity and reproduced the expected adherence phenotype of the bacteria being tested. The RAJ-IVOC model system provides a convenient method to pre-screen multiple bacteria-RAJ interactions prior to in vivo experiments, thereby reducing animal usage.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Structural visualization of transcription activated by a multidrug-sensing MerR family regulator

Bacterial RNA polymerase (RNAP) holoenzyme initiates transcription by recognizing the conserved –35 and –10 promoter elements that are optimally separated by a 17-bp spacer. The MerR family of transcriptional regulators activate suboptimal 19–20 bp spacer promoters in response to myriad cellular signals, ranging from heavy metals to drug-like compounds. The regulation of transcription by MerR family regulators is not fully understood. Here we report one crystal structure of a multidrug-sensing MerR family regulator EcmrR and nine cryo-electron microscopy structures that capture the EcmrR-dependent transcription process from promoter opening to initial transcription to RNA elongation. These structures reveal that EcmrR is a dual ligand-binding factor that reshapes the suboptimal 19-bp spacer DNA to enable optimal promoter recognition, sustains promoter remodeling to stabilize initial transcribing complexes, and finally dissociates from the promoter to reverse DNA remodeling and facilitate the transition to elongation. Our findings yield a comprehensive model for transcription regulation by MerR family factors and provide insights into the transition from transcription initiation to elongation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Site and Bioenergy Cropping System Similarly Affect Distinct Live and Total Soil Microbial Communities

Bioenergy crops are a promising energy alternative to fossil fuels. During bioenergy feedstock production, crop inputs shape the composition of soil microbial communities, which in turn influences nutrient cycling and plant productivity. In addition to cropping inputs, site characteristics (e.g., soil texture, climate) influence bacterial and fungal communities. We explored the response of soil microorganisms to bioenergy cropping system (switchgrass vs. maize) and site (sandy loam vs. silty loam) within two long-term experimental research stations. The live and total microbial community membership was investigated using 16S and ITS amplicon sequencing of soil RNA and DNA. For both nucleic acid types, we expected fungi and prokaryotes to be differentially impacted by crop and site due their dissimilar life strategies. We also expected live communities to be more strongly affected by site and crop than the total communities due to a sensitivity to recent stimuli. Instead, we found that prokaryotic and fungal community composition was primarily driven by site with a secondary crop effect, highlighting the importance of soil texture and fertility in shaping both communities. Specific highly abundant prokaryotic and fungal taxa within live communities were indicative of site and cropping systems, providing insight into treatment-specific, agriculturally relevant microbial taxa that were obscured within total community profiles. Within live prokaryote communities, predatory Myxobacteria spp. were largely indicative of silty and switchgrass communities. Within live fungal communities, Glomeromycota spp. were solely indicative of switchgrass soils, while a few very abundant Mortierellomycota spp. were indicative of silty soils. Site and cropping system had distinct effects on the live and total communities reflecting selection forces of plant inputs and environmental conditions over time. Comparisons between RNA and DNA communities uncovered live members obscured within the total community as well as members of the relic DNA pool. The associations between live communities and relic DNA are a product of the intimate relationship between the ephemeral responses of the live community and the accumulation of DNA within necromass that contributes to soil organic matter, and in turn shapes soil microbial dynamics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Horizontal Gene Transfer to a Defensive Symbiont with a Reduced Genome in a Multipartite Beetle Microbiome

Symbiotic mutualisms of bacteria and animals are ubiquitous in nature, running a continuum from facultative to obligate from the perspectives of both partners. The loss of functions required for living independently but not within a host gives rise to reduced genomes in many symbionts. Although the phenomenon of genome reduction can be explained by existing evolutionary models, the initiation of the process is not well understood. Here, we describe the microbiome associated with the eggs of the beetle Lagria villosa, consisting of multiple bacterial symbionts related to Burkholderia gladioli, including a reduced-genome symbiont thought to be the exclusive producer of the defensive compound lagriamide. We show that the putative lagriamide-producing symbiont is the only member of the microbiome undergoing genome reduction and that it has already lost the majority of its primary metabolism and DNA repair pathways. The key step preceding genome reduction in the symbiont was likely the horizontal acquisition of the putative lagriamide lga biosynthetic gene cluster. Unexpectedly, we uncovered evidence of additional horizontal transfers to the symbiont’s genome while genome reduction was occurring and despite a current lack of genes needed for homologous recombination. These gene gains may have given the genome-reduced symbiont a selective advantage in the microbiome, especially given the maintenance of the large lga gene cluster despite ongoing genome reduction.

59 BASIC BIOLOGICAL SCIENCES↗
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↗