Speeding genomic island discovery through systematic design of reference database composition.

BACKGROUND: Genomic islands (GIs) are mobile genetic elements that integrate site-specifically into bacterial chromosomes, bearing genes that affect phenotypes such as pathogenicity and metabolism. GIs typically occur sporadically among related bacterial strains, enabling comparative genomic approaches to GI identification. For a candidate GI in a query genome, the number of reference genomes with a precise deletion of the GI serves as a support value for the GI. Our comparative software for GI identification was slowed by our original use of large reference genome databases (DBs). Here we explore smaller species-focused DBs. RESULTS: With increasing DB size, recovery of our reliable prophage GI calls reached a plateau, while recovery of less reliable GI calls (FPs) increased rapidly as DB sizes exceeded ~500 genomes; i.e., overlarge DBs can increase FP rates. Paradoxically, relative to prophages, FPs were both more frequently supported only by genomes outside the species and more frequently supported only by genomes inside the species; this may be due to their generally lower support values. Setting a DB size limit for our SMAll Ranked Tailored (SMART) DB design speeded runtime ~65-fold. Strictly intra-species DBs would tend to lower yields of prophages for small species (with few genomes available); simulations with large species showed that this could be partially overcome by reaching outside the species to closely related taxa, without an FP burden. Employing such taxonomic outreach in DB design generated redundancy in the DB set; as few as 2984 DBs were needed to cover all 47894 prokaryotic species. CONCLUSIONS: Runtime decreased dramatically with SMART DB design, with only minor losses of prophages. We also describe potential utility in other comparative genomics projects.

Genome sequence and characterization of a novel Pseudomonas putida phage, MiCath.

Pseudomonads are ubiquitous bacteria with importance in medicine, soil, agriculture, and biomanufacturing. We report a novel Pseudomonas putida phage, MiCath, which is the first known phage infecting P. putida S12, a strain increasingly used as a synthetic biology chassis. MiCath was isolated from garden soil under a tomato plant using P. putida S12 as a host and was also found to infect four other P. putida strains. MiCath has a ~ 61 kbp double-stranded DNA genome which encodes 97 predicted open reading frames (ORFs); functions could only be predicted for 48 ORFs using comparative genomics. Functions include structural phage proteins, other common phage proteins (e.g., terminase), a queuosine gene cassette, a cas4 exonuclease, and an endosialidase. Restriction digestion analysis suggests the queuosine gene cassette encodes a pathway capable of modification of guanine residues. When compared to other phage genomes, MiCath shares at most 74% nucleotide identity over 2% of the genome with any sequenced phage. Overall, MiCath is a novel phage with no close relatives, encoding many unique gene products.

Characterization of novel recombinant mycobacteriophages derived from homologous recombination between two temperate phages.

Comparative analyses of mycobacteriophage genomes reveals extensive genetic diversity in genome organization and gene content, contributing to widespread mosaicism. We previously reported that the prophage of mycobacteriophage Butters (cluster N) provides defense against infection by Island3 (subcluster I1). To explore the anti-Island3 defense mechanism, we attempted to isolate Island3 defense escape mutants on a Butters lysogen, but only uncovered phages with recombinant genomes comprised of regions of Butters and Island3 arranged from left arm to right arm as Butters-Island3-Butters (BIBs). Recombination occurs within two distinct homologous regions that encompass lysin A, lysin B, and holin genes in one segment, and RecE and RecT genes in the other. Structural genes of mosaic BIB genomes are contributed by Butters while the immunity cassette is derived from Island3. Consequently, BIBs are morphologically identical to Butters (as shown by transmission electron microscopy) but are homoimmune with Island3. Recombinant phages overcome antiphage defense and silencing of the lytic cycle. We leverage this observation to propose a stratagem to generate novel phages for potential therapeutic use.

Helper-embedded satellites from an integrase clade that repeatedly targets prophage late genes.

Satellites such as phage-induced chromosomal islands (PICIs) are mobile genetic elements relying on helper phages for their mobilization, through trans-regulatory interactions. We discovered a PICI with a more intimate cis-regulatory configuration, integrated within a late gene of its helper prophage. This helper-embedded PICI (HE-PICI) configuration delays expression of the interrupted helper late gene until the satellite excises and provides passive helper-driven components to both HE-PICI replication and late transcription. Upon induction of a helper-satellite composite, precise excision of the entire composite was observed, followed by composite replication, then satellite excision. We mapped 491 additional HE-PICIs to one of 14 sites in cognates of phage lambda late genes. Associated integrases form a single phylogenetic clade with subclades respecting the 14 site groups, exhibiting repeated tropism for prophage late genes as new integration sites evolve. Four ordered zones in a general gram-negative PICI genome organization are: an integration zone encoding integrase and AlpA, a dynamic zone encoding members of the Bro-N network of domain-swapping DNA-interactive proteins and immunity repressor RNAs, a replication zone, and a dynamic late zone in which clusters as large as 17 consecutive helper prophage late genes have been captured. Helper-embedded satellites present new dimensions in satellite/helper relationships.

Genome Sequence of Mycobacteriophage Bassalto.

Bassalto is a newly isolated phage of Mycobacterium smegmatis mc2155 from the campus grounds of Norfolk State University in Norfolk, VA. Bassalto belongs to the cluster B and subcluster B3 mycobacteriophages, based on the nucleotide composition and comparison to known mycobacteriophages.

Whole-Cell MALDI-ToF MS Coupled with Untargeted Metabolomics Facilitates Investigations of Microbial Chemical Interactions.

The emergence of drug-resistant pathogens necessitates the development of new countermeasures. In this regard, the introduction of probiotics to directly attack or competitively exclude pathogens presents a useful strategy. Application of this approach requires an understanding of how a probiotic and its target pathogen interact. A key means of probiotic-pathogen interaction involves the production of small molecules called natural products (NPs). Here, we report the use of whole-cell matrix-assisted laser desorption/ionization time-of-flight (MALDI-ToF) mass spectrometry to characterize NP production by candidate probiotics (mouse airway microbiome isolates) when co-cultured with the respiratory pathogen Burkholderia. We found that a Bacillus velezensis strain inhibits growth of and elicits NP production by Burkholderia thailandensis. Dereplication of known NPs detected in the metabolome of this B. velezensis strain suggests that a previously unannotated bioactive compound is involved. Thus, we present the use of whole-cell MALDI as a broadly applicable method for screening the NP composition of microbial co-cultures; this can be combined with other -omics methods to characterize probiotic-pathogen and other microbe-microbe interactions.

Improved Mobilome Delineation in Fragmented Genomes.

The mobilome of a microbe, i.e., its set of mobile elements, has major effects on its ecology, and is important to delineate properly in each genome. This becomes more challenging for incomplete genomes, and even more so for metagenome-assembled genomes (MAGs), where misbinning of scaffolds and other losses can occur. Genomic islands (GIs), which integrate into the host chromosome, are a major component of the mobilome. Our GI-detection software TIGER, unique in its precise mapping of GI termini, was applied to 74,561 genomes from 2,473 microbial species, each species containing at least one MAG and one isolate genome. A species-normalized deficit of ∼1.6 GIs/genome was measured for MAGs relative to isolates. To test whether this undercount was due to the higher fragmentation of MAG genomes, TIGER was updated to enable detection of split GIs whose termini are on separate scaffolds or that wrap around the origin of a circular replicon. This doubled GI yields, and the new split GIs matched the quality of single-scaffold GIs, except that highly fragmented GIs may lack central portions. Cross-scaffold search is an important upgrade to GI detection as fragmented genomes increasingly dominate public databases. TIGER2 better captures MAG microdiversity, recovering niche-defining GIs and supporting microbiome research aims such as virus-host linking and ecological assessment.

BigDNA: Primer Design Software for Overlap-Based Assembly of Phage Genomes and Larger DNAs.

Background: Gibson assembly and assembly-in-yeast are strategies to create long synthetic DNAs from diverse fragments, for example, when engineering bacteriophage genomes. Design for these methods requires terminal sequence overlaps in the fragments, determining the order of assembly. Design to rebuild a genomic fragment that is too long for a single PCR presents a puzzle since some candidate joint regions cannot yield satisfactory primers for the overlap. No existing overlap assembly design software is open-source, and none explicitly supports rebuilding. Methods: We describe here bigDNA software that solves the rebuilding puzzle by recursive backtracking, with options to remove or introduce genes; it also tests for mispriming on the template DNA. BigDNA was tested with 3082 prophages and other genomic islands (GIs), from 20 to 100 kb, and the synthetic Mycoplasma genitalium genome. Results: Rebuilding assembly design succeeded for all but ∼1% of GIs. Conclusion: BigDNA will speed and standardize assembly design.

Development of potent and effective synthetic SARS-CoV-2 neutralizing nanobodies.

The respiratory virus responsible for coronavirus disease 2019 (COVID-19), severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has affected nearly every aspect of life worldwide, claiming the lives of over 3.9 million people globally, at the time of this publication. Neutralizing humanized nanobody (VHH)-based antibodies (VHH-huFc) represent a promising therapeutic intervention strategy to address the current SARS-CoV-2 pandemic and provide a powerful toolkit to address future virus outbreaks. Using a synthetic, high-diversity VHH bacteriophage library, several potent neutralizing VHH-huFc antibodies were identified and evaluated for their capacity to tightly bind to the SARS-CoV-2 receptor-binding domain, to prevent binding of SARS-CoV-2 spike (S) to the cellular receptor angiotensin-converting enzyme 2, and to neutralize viral infection. Preliminary preclinical evaluation of multiple VHH-huFc antibody candidates demonstrate that they are prophylactically and therapeutically effective in vivo against wildtype SARS-CoV-2. The identified and characterized VHH-huFc antibodies described herein represent viable candidates for further preclinical evaluation and another tool to add to our therapeutic arsenal to address the COVID-19 pandemic.

Mycobacterium Phage Butters-Encoded Proteins Contribute to Host Defense against Viral Attack.

A diverse set of prophage-mediated mechanisms protecting bacterial hosts from infection has been recently uncovered within cluster N mycobacteriophages isolated on the host, Mycobacterium smegmatis mc2155. In that context, we unveil a novel defense mechanism in cluster N prophage Butters. By using bioinformatics analyses, phage plating efficiency experiments, microscopy, and immunoprecipitation assays, we show that Butters genes located in the central region of the genome play a key role in the defense against heterotypic viral attack. Our study suggests that a two-component system, articulated by interactions between protein products of genes 30 and 31, confers defense against heterotypic phage infection by PurpleHaze (cluster A/subcluster A3) or Alma (cluster A/subcluster A9) but is insufficient to confer defense against attack by the heterotypic phage Island3 (cluster I/subcluster I1). Therefore, based on heterotypic phage plating efficiencies on the Butters lysogen, additional prophage genes required for defense are implicated and further show specificity of prophage-encoded defense systems.IMPORTANCE Many sequenced bacterial genomes, including those of pathogenic bacteria, contain prophages. Some prophages encode defense systems that protect their bacterial host against heterotypic viral attack. Understanding the mechanisms undergirding these defense systems is crucial to appreciate the scope of bacterial immunity against viral infections and will be critical for better implementation of phage therapy that would require evasion of these defenses. Furthermore, such knowledge of prophage-encoded defense mechanisms may be useful for developing novel genetic tools for engineering phage-resistant bacteria of industrial importance.

Computational Basis for On-Demand Production of Diversified Therapeutic Phage Cocktails.

New therapies are necessary to combat increasingly antibiotic-resistant bacterial pathogens. We have developed a technology platform of computational, molecular biology, and microbiology tools which together enable on-demand production of phages that target virtually any given bacterial isolate. Two complementary computational tools that identify and precisely map prophages and other integrative genetic elements in bacterial genomes are used to identify prophage-laden bacteria that are close relatives of the target strain. Phage genomes are engineered to disable lysogeny, through use of long amplicon PCR and Gibson assembly. Finally, the engineered phage genomes are introduced into host bacteria for phage production. As an initial demonstration, we used this approach to produce a phage cocktail against the opportunistic pathogen Pseudomonas aeruginosa PAO1. Two prophage-laden P. aeruginosa strains closely related to PAO1 were identified, ATCC 39324 and ATCC 27853. Deep sequencing revealed that mitomycin C treatment of these strains induced seven phages that grow on P. aeruginosa PAO1. The most diverse five phages were engineered for nonlysogeny by deleting the integrase gene (int), which is readily identifiable and typically conveniently located at one end of the prophage. The Δint phages, individually and in cocktails, killed P. aeruginosa PAO1 in liquid culture as well as in a waxworm (Galleria mellonella) model of infection.IMPORTANCE The antibiotic resistance crisis has led to renewed interest in phage therapy as an alternative means of treating infection. However, conventional methods for isolating pathogen-specific phage are slow, labor-intensive, and frequently unsuccessful. We have demonstrated that computationally identified prophages carried by near-neighbor bacteria can serve as starting material for production of engineered phages that kill the target pathogen. Our approach and technology platform offer new opportunity for rapid development of phage therapies against most, if not all, bacterial pathogens, a foundational advance for use of phage in treating infectious disease.

Genome Sequences of Burkholderia thailandensis Strains E421, E426, and DW503.

We present the draft genome sequences of three Burkholderia thailandensis strains, E421, E426, and DW503. E421 consists of 90 contigs of 6,639,935 bp and 67.73% GC content. E426 consists of 106 contigs of 6,587,853 bp and 67.73% GC content. DW503 consists of 102 contigs of 6,458,767 bp and 67.64% GC content.

New candidates for regulated gene integrity revealed through precise mapping of integrative genetic elements.

Integrative genetic elements (IGEs) are mobile multigene DNA units that integrate into and excise from host bacterial genomes. Each IGE usually targets a specific site within a conserved host gene, integrating in a manner that preserves target gene function. However, a small number of bacterial genes are known to be inactivated upon IGE integration and reactivated upon excision, regulating phenotypes of virulence, mutation rate, and terminal differentiation in multicellular bacteria. The list of regulated gene integrity (RGI) cases has been slow-growing because IGEs have been challenging to precisely and comprehensively locate in genomes. We present software (TIGER) that maps IGEs with unprecedented precision and without attB site bias. TIGER uses a comparative genomic, ping-pong BLAST approach, based on the principle that the IGE integration module (i.e. its int-attP region) is cohesive. The resultant IGEs from 2168 genomes, along with integrase phylogenetic analysis and gene inactivation tests, revealed 19 new cases of genes whose integrity is regulated by IGEs (including dut, eccCa1, gntT, hrpB, merA, ompN, prkA, tqsA, traG, yifB, yfaT and ynfE), as well as recovering previously known cases (in sigK, spsM, comK, mlrA and hlb genes). It also recovered known clades of site-promiscuous integrases and identified possible new ones.

Genome Sequences of Six Cluster N Mycobacteriophages, Kevin1, Nenae, Parmesanjohn, ShrimpFriedEgg, Smurph, and SpongeBob, Isolated on Mycobacterium smegmatis mc2155.

The annotation of six cluster N Mycobacterium smegmatis phages (Kevin1, Nenae, Parmesanjohn, ShrimpFriedEgg, Smurph, and SpongeBob) reveals regions of genomic diversity, particularly within the central region of the genome. The genome of Kevin1 includes two orphams (genes with no similarity to other phage genes), with one predicted to encode an AAA-ATPase.

Specialized eRpL22 paralogue-specific ribosomes regulate specific mRNA translation in spermatogenesis in Drosophila melanogaster.

The functional significance of ribosome heterogeneity in development and differentiation is relatively unexplored. We present the first in vivo evidence of ribosome heterogeneity playing a role in specific mRNA translation in a multicellular eukaryote. Eukaryotic-specific ribosomal protein paralogues eRpL22 and eRpL22-like are essential in development and required for sperm maturation and fertility in Drosophila. eRpL22 and eRpL22-like roles in spermatogenesis are not completely interchangeable. Flies depleted of eRpL22 and rescued by eRpL22-like overexpression have reduced fertility, confirming that eRpL22-like cannot substitute fully for eRpL22 function, and that paralogues have functionally distinct roles, not yet defined. We investigated the hypothesis that specific RNAs differentially associate with eRpL22 or eRpL22-like ribosomes, thereby establishing distinct ribosomal roles. RNA-seq identified 12,051 transcripts (mRNAs/noncoding RNAs) with 50% being enriched on specific polysome types. Analysis of ∼10% of the most abundant mRNAs suggests ribosome specialization for translating groups of mRNAs expressed at specific stages of spermatogenesis. Further, we show enrichment of "model" eRpL22-like polysome-associated testis mRNAs can occur outside the germline within S2 cells transfected with eRpL22-like, indicating that germline-specific factors are not required for selective translation. This study reveals specialized roles in translation for eRpL22 and eRpL22-like ribosomes in germline differentiation.

Functional interplay between ribosomal protein paralogues in the eRpL22 family in Drosophila melanogaster.

Duplicated ribosomal protein (RP) genes in the Drosophila melanogaster eRpL22 family encode structurally-divergent and differentially-expressed rRNA-binding RPs. eRpL22 is expressed ubiquitously and eRpL22-like expression is tissue-restricted with highest levels in the adult male germline. We explored paralogue functional equivalence using the GAL4-UAS system for paralogue knockdown or overexpression and a conditional eRpL22-like knockout in a heat- shock flippase/FRT line. Ubiquitous eRpL22 knockdown with Actin-GAL4 resulted in embryonic lethality, confirming eRpL22 essentiality. eRpL22-like knockdown (60%) was insufficient to cause lethality; yet, conditional eRpL22-like knockout at one hour following egg deposition caused lethality within each developmental stage. Therefore, each paralogue is essential. Variation in timing of heat-shock-induced eRpL22-like knockout highlighted early embryogenesis as the critical period where eRpL22-like expression (not compensated for by eRpL22) is required for normal development of several organ systems, including testis development and subsequent sperm production. To determine if eRpL22-like can substitute for eRpL22, we used Actin-GAL4 for ubiquitous eRpL22 knockdown and eRpL22-like-FLAG (or FLAG-eRpL22: control) overexpression. Emergence of adults demonstrated that ubiquitous eRpL22-like-FLAG or FLAG-eRpL22 expression eliminates embryonic lethality resulting from eRpL22 depletion. Adults rescued by eRpL22-like-FLAG (but not by FLAG-eRpL22) overexpression had reduced fertility and longevity. We conclude that eRpL22 paralogue roles are not completely interchangeable and include functionally-diverse roles in development and spermatogenesis. Testis-specific paralogue knockdown revealed molecular phenotypes, including increases in eRpL22 protein and mRNA levels following eRpL22-like depletion, implicating a negative crosstalk mechanism regulating eRpL22 expression. Paralogue depletion unmasked mechanisms, yet to be defined that impact paralogue co-expression within germ cells.

Eight Genome Sequences of Cluster BE1 Phages That Infect Streptomyces Species.

Cluster BE1 Streptomyces bacteriophages belong to the Siphoviridae, with genome sizes over 130 kbp, and they contain direct terminal repeats of approximately 11 kbp. Eight newly isolated closely related cluster BE1 phages contain 43 to 48 tRNAs, one transfer-messenger RNA (tmRNA), and 216 to 236 predicted open reading frames (ORFs), but few of their genes are shared with other phages, including those infecting Streptomyces species.

Genome Sequences of Chancellor, Mitti, and Wintermute, Three Subcluster K4 Phages Isolated Using Mycobacterium smegmatis mc2155.

Mycobacteriophages Chancellor, Mitti, and Wintermute infect Mycobacterium smegmatis mc2155 and are closely related to phages Cheetobro and Fionnbharth in subcluster K4. Genome sizes range from 57,697 bp to 58,046 bp. Phages are predicted to be temperate and to infect the pathogen Mycobacterium tuberculosis.

Genome Sequences of Mycobacteriophages Jane and Sneeze, New Members of Cluster G.

Jane and Sneeze are newly isolated phages of Mycobacterium smegmatis mc2155 from Hillsborough, NJ, and Palo Verde, Costa Rica, respectively. Both are cluster G, subcluster G1 mycobacteriophages. Notable nucleotide differences exist between genomes in the right half, including the presence of mycobacteriophage mobile element 1 (MPME1) in Jane.

Genome Sequence of Cluster W Mycobacteriophage Taptic.

The Taptic genome is the first to be annotated from the W cluster of mycobacteriophages infecting Mycobacterium smegmatis mc2155. All 92 predicted open reading frames (ORFs) and a single tRNA specifying glycine (tRNA-gly) are transcribed rightward. Many functionally uncharacterized ORFs appear to be W cluster specific, as nucleotide similarity is shared only with other W cluster genomes.