# | Rank | Similarity | Title + Abs. | Year | PMID |
|---|---|---|---|---|---|
| 0 | 1 | 2 | 3 | 4 | 5 |
| 8670 | 0 | 0.9920 | Complete Genome Analysis of Subtercola sp. PAMC28395: Genomic Insights into Its Potential Role for Cold Adaptation and Biotechnological Applications. This study reports the complete genome sequence of Subtercola sp. PAMC28395, a strain isolated from cryoconite in Uganda. This strain possesses several active carbohydrate-active enzyme (CAZyme) genes involved in glycogen and trehalose metabolism. Additionally, two specific genes associated with α-galactosidase (GH36) and bacterial alpha-1,2-mannosidase (GH92) were identified in this strain. The presence of these genes indicates the likelihood that they can be expressed, enabling the strain to break down specific polysaccharides derived from plants or the shells of nearby crabs. The authors performed a comparative analysis of CAZyme patterns and biosynthetic gene clusters (BGCs) in several Subtercola strains and provided annotations describing the unique characteristics of these strains. The comparative analysis of BGCs revealed that four strains, including PAMC28395, have oligosaccharide BGCs, and we confirmed that the pentose phosphate pathway was configured perfectly in the genome of PAMC28395, which may be associated with adaptation to low temperatures. Additionally, all strains contained antibiotic resistance genes, indicating a complex self-resistance system. These results suggest that PAMC28395 can adapt quickly to the cold environment and produce energy autonomously. This study provides valuable information on novel functional enzymes, particularly CAZymes, that operate at low temperatures and can be used for biotechnological applications and fundamental research purposes. | 2023 | 37374983 |
| 9066 | 1 | 0.9917 | VRprofile: gene-cluster-detection-based profiling of virulence and antibiotic resistance traits encoded within genome sequences of pathogenic bacteria. VRprofile is a Web server that facilitates rapid investigation of virulence and antibiotic resistance genes, as well as extends these trait transfer-related genetic contexts, in newly sequenced pathogenic bacterial genomes. The used backend database MobilomeDB was firstly built on sets of known gene cluster loci of bacterial type III/IV/VI/VII secretion systems and mobile genetic elements, including integrative and conjugative elements, prophages, class I integrons, IS elements and pathogenicity/antibiotic resistance islands. VRprofile is thus able to co-localize the homologs of these conserved gene clusters using HMMer or BLASTp searches. With the integration of the homologous gene cluster search module with a sequence composition module, VRprofile has exhibited better performance for island-like region predictions than the other widely used methods. In addition, VRprofile also provides an integrated Web interface for aligning and visualizing identified gene clusters with MobilomeDB-archived gene clusters, or a variety set of bacterial genomes. VRprofile might contribute to meet the increasing demands of re-annotations of bacterial variable regions, and aid in the real-time definitions of disease-relevant gene clusters in pathogenic bacteria of interest. VRprofile is freely available at http://bioinfo-mml.sjtu.edu.cn/VRprofile. | 2018 | 28077405 |
| 9074 | 2 | 0.9915 | BacAnt: A Combination Annotation Server for Bacterial DNA Sequences to Identify Antibiotic Resistance Genes, Integrons, and Transposable Elements. Whole genome sequencing (WGS) of bacteria has become a routine method in diagnostic laboratories. One of the clinically most useful advantages of WGS is the ability to predict antimicrobial resistance genes (ARGs) and mobile genetic elements (MGEs) in bacterial sequences. This allows comprehensive investigations of such genetic features but can also be used for epidemiological studies. A plethora of software programs have been developed for the detailed annotation of bacterial DNA sequences, such as rapid annotation using subsystem technology (RAST), Resfinder, ISfinder, INTEGRALL and The Transposon Registry. Unfortunately, to this day, a reliable annotation tool of the combination of ARGs and MGEs is not available, and the generation of genbank files requires much manual input. Here, we present a new webserver which allows the annotation of ARGs, integrons and transposable elements at the same time. The pipeline generates genbank files automatically, which are compatible with Easyfig for comparative genomic analysis. Our BacAnt code and standalone software package are available at https://github.com/xthua/bacant with an accompanying web application at http://bacant.net. | 2021 | 34367079 |
| 9072 | 3 | 0.9915 | PanGeT: Pan-genomics tool. A decade after the concept of Pan-genome was first introduced; research in this field has spread its tentacles to areas such as pathogenesis of diseases, bacterial evolutionary studies and drug resistance. Gene content-based differentiation of virulent and a virulent strains of bacteria and identification of pathogen specific genes is imperative to understand their physiology and gain insights into the mechanism of genome evolution. Subsequently, this will aid in identifying diagnostic targets and in developing and selecting vaccines. The root of pan-genomic studies, however, is to identify the core genes, dispensable genes and strain specific genes across the genomes belonging to a clade. To this end, we have developed a tool, "PanGeT - Pan-genomics Tool" to compute the 'pan-genome' based on comparisons at the genome as well as the proteome levels. This automated tool is implemented using LaTeX libraries for effective visualization of overall pan-genome through graphical plots. Links to retrieve sequence information and functional annotations have also been provided. PanGeT can be downloaded from http://pranag.physics.iisc.ernet.in/PanGeT/ or https://github.com/PanGeTv1/PanGeT. | 2017 | 27851981 |
| 9067 | 4 | 0.9914 | PIPdb: a comprehensive plasmid sequence resource for tracking the horizontal transfer of pathogenic factors and antimicrobial resistance genes. Plasmids, as independent genetic elements, carrying resistance or virulence genes and transfer them among different pathogens, posing a significant threat to human health. Under the 'One Health' approach, it is crucial to control the spread of plasmids carrying such genes. To achieve this, a comprehensive characterization of plasmids in pathogens is essential. Here we present the Plasmids in Pathogens Database (PIPdb), a pioneering resource that includes 792 964 plasmid segment clusters (PSCs) derived from 1 009 571 assembled genomes across 450 pathogenic species from 110 genera. To our knowledge, PIPdb is the first database specifically dedicated to plasmids in pathogenic bacteria, offering detailed multi-dimensional metadata such as collection date, geographical origin, ecosystem, host taxonomy, and habitat. PIPdb also provides extensive functional annotations, including plasmid type, insertion sequences, integron, oriT, relaxase, T4CP, virulence factors genes, heavy metal resistance genes and antibiotic resistance genes. The database features a user-friendly interface that facilitates studies on plasmids across diverse host taxa, habitats, and ecosystems, with a focus on those carrying antimicrobial resistance genes (ARGs). We have integrated online tools for plasmid identification and annotation from assembled genomes. Additionally, PIPdb includes a risk-scoring system for identifying potentially high-risk plasmids. The PIPdb web interface is accessible at https://nmdc.cn/pipdb. | 2025 | 39460620 |
| 8259 | 5 | 0.9913 | Secondary Metabolite Transcriptomic Pipeline (SeMa-Trap), an expression-based exploration tool for increased secondary metabolite production in bacteria. For decades, natural products have been used as a primary resource in drug discovery pipelines to find new antibiotics, which are mainly produced as secondary metabolites by bacteria. The biosynthesis of these compounds is encoded in co-localized genes termed biosynthetic gene clusters (BGCs). However, BGCs are often not expressed under laboratory conditions. Several genetic manipulation strategies have been developed in order to activate or overexpress silent BGCs. Significant increases in production levels of secondary metabolites were indeed achieved by modifying the expression of genes encoding regulators and transporters, as well as genes involved in resistance or precursor biosynthesis. However, the abundance of genes encoding such functions within bacterial genomes requires prioritization of the most promising ones for genetic manipulation strategies. Here, we introduce the 'Secondary Metabolite Transcriptomic Pipeline' (SeMa-Trap), a user-friendly web-server, available at https://sema-trap.ziemertlab.com. SeMa-Trap facilitates RNA-Seq based transcriptome analyses, finds co-expression patterns between certain genes and BGCs of interest, and helps optimize the design of comparative transcriptomic analyses. Finally, SeMa-Trap provides interactive result pages for each BGC, allowing the easy exploration and comparison of expression patterns. In summary, SeMa-Trap allows a straightforward prioritization of genes that could be targeted via genetic engineering approaches to (over)express BGCs of interest. | 2022 | 35580059 |
| 9219 | 6 | 0.9910 | Knowing and Naming: Phage Annotation and Nomenclature for Phage Therapy. Bacteriophages, or phages, are viruses that infect bacteria shaping microbial communities and ecosystems. They have gained attention as potential agents against antibiotic resistance. In phage therapy, lytic phages are preferred for their bacteria killing ability, while temperate phages, which can transfer antibiotic resistance or toxin genes, are avoided. Selection relies on plaque morphology and genome sequencing. This review outlines annotating genomes, identifying critical genomic features, and assigning functional labels to protein-coding sequences. These annotations prevent the transfer of unwanted genes, such as antimicrobial resistance or toxin genes, during phage therapy. Additionally, it covers International Committee on Taxonomy of Viruses (ICTV)-an established phage nomenclature system for simplified classification and communication. Accurate phage genome annotation and nomenclature provide insights into phage-host interactions, replication strategies, and evolution, accelerating our understanding of the diversity and evolution of phages and facilitating the development of phage-based therapies. | 2023 | 37932119 |
| 4354 | 7 | 0.9909 | ARDB--Antibiotic Resistance Genes Database. The treatment of infections is increasingly compromised by the ability of bacteria to develop resistance to antibiotics through mutations or through the acquisition of resistance genes. Antibiotic resistance genes also have the potential to be used for bio-terror purposes through genetically modified organisms. In order to facilitate the identification and characterization of these genes, we have created a manually curated database--the Antibiotic Resistance Genes Database (ARDB)--unifying most of the publicly available information on antibiotic resistance. Each gene and resistance type is annotated with rich information, including resistance profile, mechanism of action, ontology, COG and CDD annotations, as well as external links to sequence and protein databases. Our database also supports sequence similarity searches and implements an initial version of a tool for characterizing common mutations that confer antibiotic resistance. The information we provide can be used as compendium of antibiotic resistance factors as well as to identify the resistance genes of newly sequenced genes, genomes, or metagenomes. Currently, ARDB contains resistance information for 13,293 genes, 377 types, 257 antibiotics, 632 genomes, 933 species and 124 genera. ARDB is available at http://ardb.cbcb.umd.edu/. | 2009 | 18832362 |
| 6125 | 8 | 0.9909 | Complete Genome Sequence Analysis of Brevibacillus laterosporus Bl-zj Reflects its Potential Algicidal Response. We analyzed the complete genome of the bacteria Brevibacillus laterosporus Bl-zj. Its genome has a total length of 5,202,546 bp with 4594 annotated genes. The functional groups included transporters, pathogen-host interaction factors, antibiotic resistance genes, virulence factor, and secreted proteins were predicted, and carbon and nitrogen metabolism and transporters were mapped. A total of 34 genes possibly involved in algae-lysing processes were further screened, including 8 virulence factors, 18 secreted proteases, and 8 antibiotic-resistant genes, which could be playing important roles in host identification, invasion, and the destruction of algal cells. This study will provide a theoretical framework for the algicidal mechanism of algae-lysing bacteria and possible application to algal control. | 2021 | 33649996 |
| 8470 | 9 | 0.9908 | Genomic Insights into Vaccinium spp. Endophytes B. halotolerans and B. velezensis and Their Antimicrobial Potential. Plant microbiota contributes to nutrient absorption, and the production of hormones and vitamins, and plays a crucial role in responding to environmental stress. We hypothesized that Vaccinium spp. harbour a unique microbiota that enables them to coexist in extreme environments such as saline, nutrient-poor, and waterlogged conditions. Upon examining Bacillus spp. endophytes isolated from blueberries, cranberries and lingonberries in vitro, we identified B. halotolerans (Bil-LT1_1, Bil-LT1_2) and B. velezensis (Cran-LT1_8, Ling-NOR4_15) strains that inhibit the growth of five pathogenic fungi and five foodborne bacteria. Whole-genome sequencing provided insights into genome organization and plasticity, helping identify mobile elements and genes potentially acquired through horizontal gene transfer. Functional annotation identified genes associated with plant colonization, stress tolerance, biocontrol activity, and plant growth promotion. Comparative genomic analyses revealed key biosynthetic gene clusters (BGCs) responsible for producing antifungal metabolites, including lipopeptides and polyketides. Genes supporting plant nutrition, growth, and environmental adaptation were present also in these strains. Notably, isolated endophytes exhibited particularly high levels of genomic plasticity, likely due to horizontal gene transfer involving gene ontology (GO) pathways related to survival in polymicrobial and foreign environments. | 2025 | 40724928 |
| 4358 | 10 | 0.9908 | Genomic profiling of pediococcus acidilactici BCB1H and identification of its key features for Biotechnological innovation, food technology and medicine. Lactic acid bacteria has been extensively used in food industry because of widespread properties and Pediococcus is among one of them. This study aims to conduct a comprehensive genomic analysis of Pediococcus acidilactici strain BCB1H to elucidate its genetic composition, functional elements, and potential biotechnological applications. The objectives include identifying key genomic features such as coding sequences, tRNA and rRNA genes, antibiotic resistance genes, and secondary metabolite biosynthetic gene clusters, which will highlight the adaptability and potential of P. acidilactici strain BCB1H for use in a variety of industrial and therapeutic applications. P. acidilactici strain BCB1H was analyzed using whole-genome sequencing, which used advanced sequencing technologies to obtain comprehensive genomic data. Key genomic features, such as coding sequences, tRNA and rRNA genes, antibiotic resistance genes, and secondary metabolite biosynthetic gene clusters, were identified through bioinformatics analyses. The genomic analysis of P. acidilactici strain BCB1H revealed a genome size of approximately 1.92 million base pairs with a GC content of 42.4%. The annotation identified 1,895 genes across 192 subsystems, highlighting the metabolic pathways and functional categories. Notably, specialty genes associated with carbohydrate metabolism, stress response, pathogenicity, and amino acid synthesis were identified, underscoring the versatility and potential applications in food technology and medicine. These findings shed light on the genetic makeup and functional potential of P. acidilactici strain BCB1H, highlighting its flexibility and industrial importance. The genetic traits discovered suggest its prospective use in probiotics, food preservation, and biotechnological advancements. | 2025 | 39971970 |
| 9077 | 11 | 0.9908 | The PLSDB 2025 update: enhanced annotations and improved functionality for comprehensive plasmid research. Plasmids are extrachromosomal DNA molecules in bacteria and archaea, playing critical roles in horizontal gene transfer, antibiotic resistance, and pathogenicity. Since its first release in 2018, our database on plasmids, PLSDB, has significantly grown and enhanced its content and scope. From 34 513 records contained in the 2021 version, PLSDB now hosts 72 360 entries. Designed to provide life scientists with convenient access to extensive plasmid data and to support computer scientists by offering curated datasets for artificial intelligence (AI) development, this latest update brings more comprehensive and accurate information for plasmid research, with interactive visualization options. We enriched PLSDB by refining the identification and classification of plasmid host ecosystems and host diseases. Additionally, we incorporated annotations for new functional structures, including protein-coding genes and biosynthetic gene clusters. Further, we enhanced existing annotations, such as antimicrobial resistance genes and mobility typing. To accommodate these improvements and to host the increase plasmid sets, the webserver architecture and underlying data structures of PLSDB have been re-reconstructed, resulting in decreased response times and enhanced visualization of features while ensuring that users have access to a more efficient and user-friendly interface. The latest release of PLSDB is freely accessible at https://www.ccb.uni-saarland.de/plsdb2025. | 2025 | 39565221 |
| 9079 | 12 | 0.9908 | Review, Evaluation, and Directions for Gene-Targeted Assembly for Ecological Analyses of Metagenomes. Shotgun metagenomics has greatly advanced our understanding of microbial communities over the last decade. Metagenomic analyses often include assembly and genome binning, computationally daunting tasks especially for big data from complex environments such as soil and sediments. In many studies, however, only a subset of genes and pathways involved in specific functions are of interest; thus, it is not necessary to attempt global assembly. In addition, methods that target genes can be computationally more efficient and produce more accurate assembly by leveraging rich databases, especially for those genes that are of broad interest such as those involved in biogeochemical cycles, biodegradation, and antibiotic resistance or used as phylogenetic markers. Here, we review six gene-targeted assemblers with unique algorithms for extracting and/or assembling targeted genes: Xander, MegaGTA, SAT-Assembler, HMM-GRASPx, GenSeed-HMM, and MEGAN. We tested these tools using two datasets with known genomes, a synthetic community of artificial reads derived from the genomes of 17 bacteria, shotgun sequence data from a mock community with 48 bacteria and 16 archaea genomes, and a large soil shotgun metagenomic dataset. We compared assemblies of a universal single copy gene (rplB) and two N cycle genes (nifH and nirK). We measured their computational efficiency, sensitivity, specificity, and chimera rate and found Xander and MegaGTA, which both use a probabilistic graph structure to model the genes, have the best overall performance with all three datasets, although MEGAN, a reference matching assembler, had better sensitivity with synthetic and mock community members chosen from its reference collection. Also, Xander and MegaGTA are the only tools that include post-assembly scripts tuned for common molecular ecology and diversity analyses. Additionally, we provide a mathematical model for estimating the probability of assembling targeted genes in a metagenome for estimating required sequencing depth. | 2019 | 31749830 |
| 8421 | 13 | 0.9907 | Dynamic stepwise opening of integron attC DNA hairpins by SSB prevents toxicity and ensures functionality. Biologically functional DNA hairpins are found in archaea, prokaryotes and eukaryotes, playing essential roles in various DNA transactions. However, during DNA replication, hairpin formation can stall the polymerase and is therefore prevented by the single-stranded DNA binding protein (SSB). Here, we address the question how hairpins maintain their functional secondary structure despite SSB's presence. As a model hairpin, we used the recombinogenic form of the attC site, essential for capturing antibiotic-resistance genes in the integrons of bacteria. We found that attC hairpins have a conserved high GC-content near their apical loop that creates a dynamic equilibrium between attC fully opened by SSB and a partially structured attC-6-SSB complex. This complex is recognized by the recombinase IntI, which extrudes the hairpin upon binding while displacing SSB. We anticipate that this intriguing regulation mechanism using a base pair distribution to balance hairpin structure formation and genetic stability is key to the dissemination of antibiotic resistance genes among bacteria and might be conserved among other functional hairpins. | 2017 | 28985409 |
| 9071 | 14 | 0.9907 | RAC: Repository of Antibiotic resistance Cassettes. Antibiotic resistance in bacteria is often due to acquisition of resistance genes associated with different mobile genetic elements. In Gram-negative bacteria, many resistance genes are found as part of small mobile genetic elements called gene cassettes, generally found integrated into larger elements called integrons. Integrons carrying antibiotic resistance gene cassettes are often associated with mobile elements and here are designated 'mobile resistance integrons' (MRIs). More than one cassette can be inserted in the same integron to create arrays that contribute to the spread of multi-resistance. In many sequences in databases such as GenBank, only the genes within cassettes, rather than whole cassettes, are annotated and the same gene/cassette may be given different names in different entries, hampering analysis. We have developed the Repository of Antibiotic resistance Cassettes (RAC) website to provide an archive of gene cassettes that includes alternative gene names from multiple nomenclature systems and allows the community to contribute new cassettes. RAC also offers an additional function that allows users to submit sequences containing cassettes or arrays for annotation using the automatic annotation system Attacca. Attacca recognizes features (gene cassettes, integron regions) and identifies cassette arrays as patterns of features and can also distinguish minor cassette variants that may encode different resistance phenotypes (aacA4 cassettes and bla cassettes-encoding β-lactamases). Gaps in annotations are manually reviewed and those found to correspond to novel cassettes are assigned unique names. While there are other websites dedicated to integrons or antibiotic resistance genes, none includes a complete list of antibiotic resistance gene cassettes in MRI or offers consistent annotation and appropriate naming of all of these cassettes in submitted sequences. RAC thus provides a unique resource for researchers, which should reduce confusion and improve the quality of annotations of gene cassettes in integrons associated with antibiotic resistance. DATABASE URL: http://www2.chi.unsw.edu.au/rac. | 2011 | 22140215 |
| 6144 | 15 | 0.9907 | Efficient arsenate reduction by As-resistant bacterium Bacillus sp. strain PVR-YHB1-1: Characterization and genome analysis. Arsenate (AsV) reduction in bacteria is essential to alleviate their arsenic (As) toxicity. We isolated a Bacillus strain PVR-YHB1-1 from the roots of As-hyperaccumulator Pteris vittata. The strain was efficient in reducing AsV to arsenite (AsIII), but the associated mechanisms were unclear. Here, we investigated its As resistance and reduction behaviors and associated genes at genome level. Results showed that the strain tolerated up to 20 mM AsV. When grown in 1 mM AsV, 96% AsV was reduced to AsIII in 48 h, with its AsV reduction ability being positively correlated to bacterial biomass. Two ars operons arsRacr3arsCDA and arsRKacr3arsC for As metabolisms were identified based on draft genome sequencing and gene annotations. Our data suggested that both operons might have attributed to efficient As resistance and AsV reduction in PVR-YHB1-1, providing clues to better understand As transformation in bacteria and their roles in As transformation in the environment. | 2019 | 30609485 |
| 6139 | 16 | 0.9906 | Complete genome and two plasmids sequences of Lactiplantibacillus plantarum L55 for probiotic potentials. In this study, we report the complete genome sequence of Lactiplantibacillus plantarum L55, a probiotic strain of lactic acid bacteria isolated from kimchi. The genome consists of one circular chromosome (2,077,416 base pair [bp]) with a guanine cytosine (GC) content of 44.5%, and two circular plasmid sequences (54,267 and 19,592 bp, respectively). We also conducted a comprehensive analysis of the genome, which identified the presence of functional genes, genomic islands, and antibiotic-resistance genes. The genome sequence data presented in this study provide insights into the genetic basis of L. plantarum L55, which could be beneficial for the future development of probiotic applications. | 2023 | 38616876 |
| 9666 | 17 | 0.9905 | The comprehensive antibiotic resistance database. The field of antibiotic drug discovery and the monitoring of new antibiotic resistance elements have yet to fully exploit the power of the genome revolution. Despite the fact that the first genomes sequenced of free living organisms were those of bacteria, there have been few specialized bioinformatic tools developed to mine the growing amount of genomic data associated with pathogens. In particular, there are few tools to study the genetics and genomics of antibiotic resistance and how it impacts bacterial populations, ecology, and the clinic. We have initiated development of such tools in the form of the Comprehensive Antibiotic Research Database (CARD; http://arpcard.mcmaster.ca). The CARD integrates disparate molecular and sequence data, provides a unique organizing principle in the form of the Antibiotic Resistance Ontology (ARO), and can quickly identify putative antibiotic resistance genes in new unannotated genome sequences. This unique platform provides an informatic tool that bridges antibiotic resistance concerns in health care, agriculture, and the environment. | 2013 | 23650175 |
| 8452 | 18 | 0.9905 | Functional differentiation and spatial-temporal co-expression networks of the NBS-encoding gene family in Jilin ginseng, Panax ginseng C.A. Meyer. Ginseng, Panax ginseng C.A. Meyer, is one of the most important medicinal plants for human health and medicine. It has been documented that over 80% of genes conferring resistance to bacteria, viruses, fungi and nematodes are contributed by the nucleotide binding site (NBS)-encoding gene family. Therefore, identification and characterization of NBS genes expressed in ginseng are paramount to its genetic improvement and breeding. However, little is known about the NBS-encoding genes in ginseng. Here we report genome-wide identification and systems analysis of the NBS genes actively expressed in ginseng (PgNBS genes). Four hundred twelve PgNBS gene transcripts, derived from 284 gene models, were identified from the transcriptomes of 14 ginseng tissues. These genes were classified into eight types, including TNL, TN, CNL, CN, NL, N, RPW8-NL and RPW8-N. Seven conserved motifs were identified in both the Toll/interleukine-1 receptor (TIR) and coiled-coil (CC) typed genes whereas six were identified in the RPW8 typed genes. Phylogenetic analysis showed that the PgNBS gene family is an ancient family, with a vast majority of its genes originated before ginseng originated. In spite of their belonging to a family, the PgNBS genes have functionally dramatically differentiated and been categorized into numerous functional categories. The expressions of the across tissues, different aged roots and the roots of different genotypes. However, they are coordinating in expression, forming a single co-expression network. These results provide a deeper understanding of the origin, evolution and functional differentiation and expression dynamics of the NBS-encoding gene family in plants in general and in ginseng particularly, and a NBS gene toolkit useful for isolation and characterization of disease resistance genes and for enhanced disease resistance breeding in ginseng and related species. | 2017 | 28727829 |
| 8408 | 19 | 0.9905 | The defensome of prokaryotes in aquifers. Groundwater harbors a pristine biosphere where microbes co-evolve with less human interference, yet the ancient and ongoing arms race between prokaryotes and viruses remains largely unknown in such ecosystems. Based on our recent nationwide groundwater monitoring campaign across China, we construct a metagenomic groundwater prokaryotic defensome catalogue (GPDC), encompassing 190,810 defense genes, 90,824 defense systems, 139 defense families, and 669 defense islands from 141 prokaryotic phyla. Over 94% of the defense genes in GPDC are novel and contribute vast microbial immune resources in groundwater. We find that candidate phyla radiation (CPR) bacteria possess higher defense system density and diversity against intense phage infection, while microbes as a whole exhibit an inverse relationship between defense systems and adaptive traits like resistance genes in groundwater. We further identify five first-line defense families covering 69.2% of the total defense systems, and high-turnover accessory immune genes are mostly conveyed to defense islands by mobile genetic elements. Our study also reveals viral resistance to microbial defense through co-localized anti-defense genes and interactions between CRISPR-Cas9 and anti-CRISPR protein. These findings expand our understanding of microbial immunity in pristine ecosystems and offer valuable immune resources for potential biotechnological applications. | 2025 | 40659683 |