Recommendations for Uniform Variant Calling of SARS-CoV-2 Genome Sequence across Bioinformatic Workflows.

Genomic sequencing of clinical samples to identify emerging variants of SARS-CoV-2 has been a key public health tool for curbing the spread of the virus. As a result, an unprecedented number of SARS-CoV-2 genomes were sequenced during the COVID-19 pandemic, which allowed for rapid identification of genetic variants, enabling the timely design and testing of therapies and deployment of new vaccine formulations to combat the new variants. However, despite the technological advances of deep sequencing, the analysis of the raw sequence data generated globally is neither standardized nor consistent, leading to vastly disparate sequences that may impact identification of variants. Here, we show that for both Illumina and Oxford Nanopore sequencing platforms, downstream bioinformatic protocols used by industry, government, and academic groups resulted in different virus sequences from same sample. These bioinformatic workflows produced consensus genomes with differences in single nucleotide polymorphisms, inclusion and exclusion of insertions, and/or deletions, despite using the same raw sequence as input datasets. Here, we compared and characterized such discrepancies and propose a specific suite of parameters and protocols that should be adopted across the field. Consistent results from bioinformatic workflows are fundamental to SARS-CoV-2 and future pathogen surveillance efforts, including pandemic preparation, to allow for a data-driven and timely public health response.

Host and nonhost bacteria support bacteriophage dissemination along mycelia and abiotic dispersal networks.

Bacteriophages play a crucial role in shaping bacterial communities, yet the mechanisms by which nonmotile bacteriophages interact with their hosts remain poorly understood. This knowledge gap is especially pronounced in structured environments like soil, where spatial constraints and air-filled zones hinder aqueous diffusion. In soil, hyphae of filamentous microorganisms form a network of 'fungal highways' (FHs) that facilitate the dispersal of other microorganisms. We propose that FHs also promote bacteriophage dissemination. Viral particles can diffuse in liquid films surrounding hyphae or be transported by infectable (host) or uninfectable (nonhost) bacterial carriers coexisting on FH networks. To test this, two bacteriophages that infect Pseudomonas putida DSM291 (host) but not KT2440 (nonhost) were used. In the absence of carriers, bacteriophages showed limited diffusion on 3D-printed abiotic networks, but diffusion was significantly improved in Pythium ultimum-formed FHs when the number of connecting hyphae exceeded 20. Transport by both host and nonhost carriers enhanced bacteriophage dissemination. Host carriers were five times more effective in transporting bacteriophages, particularly in FHs with over 30 connecting hyphae. This study enhances our understanding of bacteriophage dissemination in nonsaturated environments like soils, highlighting the importance of biotic networks and bacterial hosts in facilitating this process.

Bridging Place-Based Astrobiology Education with Genomics, Including Descriptions of Three Novel Bacterial Species Isolated from Mars Analog Sites of Cultural Relevance.

Democratizing genomic data science, including bioinformatics, can diversify the STEM workforce and may, in turn, bring new perspectives into the space sciences. In this respect, the development of education and research programs that bridge genome science with "place" and world-views specific to a given region are valuable for Indigenous students and educators. Through a multi-institutional collaboration, we developed an ongoing education program and model that includes Illumina and Oxford Nanopore sequencing, free bioinformatic platforms, and teacher training workshops to address our research and education goals through a place-based science education lens. High school students and researchers cultivated, sequenced, assembled, and annotated the genomes of 13 bacteria from Mars analog sites with cultural relevance, 10 of which were novel species. Students, teachers, and community members assisted with the discovery of new, potentially chemolithotrophic bacteria relevant to astrobiology. This joint education-research program also led to the discovery of species from Mars analog sites capable of producing N-acyl homoserine lactones, which are quorum-sensing molecules used in bacterial communication. Whole genome sequencing was completed in high school classrooms, and connected students to funded space research, increased research output, and provided culturally relevant, place-based science education, with participants naming three novel species described here. Students at St. Andrew's School (Honolulu, Hawai'i) proposed the name Bradyrhizobium prioritasuperba for the type strain, BL16AT, of the new species (DSM 112479T = NCTC 14602T). The nonprofit organization Kauluakalana proposed the name Brenneria ulupoensis for the type strain, K61T, of the new species (DSM 116657T = LMG = 33184T), and Hawai'i Baptist Academy students proposed the name Paraflavitalea speifideiaquila for the type strain, BL16ET, of the new species (DSM 112478T = NCTC 14603T).

Associated bacterial communities, confrontation studies, and comparative genomics reveal important interactions between Morchella with Pseudomonas spp.

Members of the fungal genus Morchella are widely known for their important ecological roles and significant economic value. In this study, we used amplicon and genome sequencing to characterize bacterial communities associated with sexual fruiting bodies from wild specimens, as well as vegetative mycelium and sclerotia obtained from Morchella isolates grown in vitro. These investigations included diverse representatives from both Elata and Esculenta Morchella clades. Unique bacterial community compositions were observed across the various structures examined, both within and across individual Morchella isolates or specimens. However, specific bacterial taxa were frequently detected in association with certain structures, providing support for an associated core bacterial community. Bacteria from the genus Pseudomonas and Ralstonia constituted the core bacterial associates of Morchella mycelia and sclerotia, while other genera (e.g., Pedobacter spp., Deviosa spp., and Bradyrhizobium spp.) constituted the core bacterial community of fruiting bodies. Furthermore, the importance of Pseudomonas as a key member of the bacteriome was supported by the isolation of several Pseudomonas strains from mycelia during in vitro cultivation. Four of the six mycelial-derived Pseudomonas isolates shared 16S rDNA sequence identity with amplicon sequences recovered directly from the examined fungal structures. Distinct interaction phenotypes (antagonistic or neutral) were observed in confrontation assays between these bacteria and various Morchella isolates. Genome sequences obtained from these Pseudomonas isolates revealed intriguing differences in gene content and annotated functions, specifically with respect to toxin-antitoxin systems, cell adhesion, chitinases, and insecticidal toxins. These genetic differences correlated with the interaction phenotypes. This study provides evidence that Pseudomonas spp. are frequently associated with Morchella and these associations may greatly impact fungal physiology.

Genomic and Phenotypic Characterization of Shiga Toxin-Producing Escherichia albertii Strains Isolated from Wild Birds in a Major Agricultural Region in California.

Escherichia albertii is an emerging foodborne pathogen. To better understand the pathogenesis and health risk of this pathogen, comparative genomics and phenotypic characterization were applied to assess the pathogenicity potential of E. albertii strains isolated from wild birds in a major agricultural region in California. Shiga toxin genes stx2f were present in all avian strains. Pangenome analyses of 20 complete genomes revealed a total of 11,249 genes, of which nearly 80% were accessory genes. Both core gene-based phylogenetic and accessory gene-based relatedness analyses consistently grouped the three stx2f-positive clinical strains with the five avian strains carrying ST7971. Among the three Stx2f-converting prophage integration sites identified, ssrA was the most common one. Besides the locus of enterocyte effacement and type three secretion system, the high pathogenicity island, OI-122, and type six secretion systems were identified. Substantial strain variation in virulence gene repertoire, Shiga toxin production, and cytotoxicity were revealed. Six avian strains exhibited significantly higher cytotoxicity than that of stx2f-positive E. coli, and three of them exhibited a comparable level of cytotoxicity with that of enterohemorrhagic E. coli outbreak strains, suggesting that wild birds could serve as a reservoir of E. albertii strains with great potential to cause severe diseases in humans.

Fungal drops: a novel approach for macro- and microscopic analyses of fungal mycelial growth.

This study presents an inexpensive approach for the macro- and microscopic observation of fungal mycelial growth. The 'fungal drops' method allows to investigate the development of a mycelial network in filamentous microorganisms at the colony and hyphal scales. A heterogeneous environment is created by depositing 15-20 µl drops on a hydrophobic surface at a fixed distance. This system is akin to a two-dimensional (2D) soil-like structure in which aqueous-pockets are intermixed with air-filled pores. The fungus (spores or mycelia) is inoculated into one of the drops, from which hyphal growth and exploration take place. Hyphal structures are assessed at different scales using stereoscopic and microscopic imaging. The former allows to evaluate the local response of regions within the colony (modular behaviour), while the latter can be used for fractal dimension analyses to describe the hyphal network architecture. The method was tested with several species to underpin the transferability to multiple species. In addition, two sets of experiments were carried out to demonstrate its use in fungal biology. First, mycelial reorganization of Fusarium oxysporum was assessed as a response to patches containing different nutrient concentrations. Second, the effect of interactions with the soil bacterium Pseudomonas putida on habitat colonization by the same fungus was assessed. This method appeared as fast and accessible, allowed for a high level of replication, and complements more complex experimental platforms. Coupled with image analysis, the fungal drops method provides new insights into the study of fungal modularity both macroscopically and at a single-hypha level.

Metabolic versatility of Caldarchaeales from geothermal features of Hawai'i and Chile as revealed by five metagenome-assembled genomes.

Members of the archaeal order Caldarchaeales (previously the phylum Aigarchaeota) are poorly sampled and are represented in public databases by relatively few genomes. Additional representative genomes will help resolve their placement among all known members of Archaea and provide insights into their roles in the environment. In this study, we analyzed 16S rRNA gene amplicons belonging to the Caldarchaeales that are available in public databases, which demonstrated that archaea of the order Caldarchaeales are diverse, widespread, and most abundant in geothermal habitats. We also constructed five metagenome-assembled genomes (MAGs) of Caldarchaeales from two geothermal features to investigate their metabolic potential and phylogenomic position in the domain Archaea. Two of the MAGs were assembled from microbial community DNA extracted from fumarolic lava rocks from Mauna Ulu, Hawai'i, and three were assembled from DNA obtained from hot spring sinters from the El Tatio geothermal field in Chile. MAGs from Hawai'i are high quality bins with completeness >95% and contamination <1%, and one likely belongs to a novel species in a new genus recently discovered at a submarine volcano off New Zealand. MAGs from Chile have lower completeness levels ranging from 27 to 70%. Gene content of the MAGs revealed that these members of Caldarchaeales are likely metabolically versatile and exhibit the potential for both chemoorganotrophic and chemolithotrophic lifestyles. The wide array of metabolic capabilities exhibited by these members of Caldarchaeales might help them thrive under diverse harsh environmental conditions. All the MAGs except one from Chile harbor putative prophage regions encoding several auxiliary metabolic genes (AMGs) that may confer a fitness advantage on their Caldarchaeales hosts by increasing their metabolic potential and make them better adapted to new environmental conditions. Phylogenomic analysis of the five MAGs and over 3,000 representative archaeal genomes showed the order Caldarchaeales forms a monophyletic group that is sister to the clade comprising the orders Geothermarchaeales (previously Candidatus Geothermarchaeota), Conexivisphaerales and Nitrososphaerales (formerly known as Thaumarchaeota), supporting the status of Caldarchaeales members as a clade distinct from the Thaumarchaeota.

Identification of mobile genetic elements with geNomad.

Identifying and characterizing mobile genetic elements in sequencing data is essential for understanding their diversity, ecology, biotechnological applications and impact on public health. Here we introduce geNomad, a classification and annotation framework that combines information from gene content and a deep neural network to identify sequences of plasmids and viruses. geNomad uses a dataset of more than 200,000 marker protein profiles to provide functional gene annotation and taxonomic assignment of viral genomes. Using a conditional random field model, geNomad also detects proviruses integrated into host genomes with high precision. In benchmarks, geNomad achieved high classification performance for diverse plasmids and viruses (Matthews correlation coefficient of 77.8% and 95.3%, respectively), substantially outperforming other tools. Leveraging geNomad's speed and scalability, we processed over 2.7 trillion base pairs of sequencing data, leading to the discovery of millions of viruses and plasmids that are available through the IMG/VR and IMG/PR databases. geNomad is available at https://portal.nersc.gov/genomad .

Comparative genomics of Mollicutes-related endobacteria supports a late invasion into Mucoromycota fungi.

Diverse members of early-diverging Mucoromycota, including mycorrhizal taxa and soil-associated Mortierellaceae, are known to harbor Mollicutes-related endobacteria (MRE). It has been hypothesized that MRE were acquired by a common ancestor and transmitted vertically. Alternatively, MRE endosymbionts could have invaded after the divergence of Mucoromycota lineages and subsequently spread to new hosts horizontally. To better understand the evolutionary history of MRE symbionts, we generated and analyzed four complete MRE genomes from two Mortierellaceae genera: Linnemannia (MRE-L) and Benniella (MRE-B). These genomes include the smallest known of fungal endosymbionts and showed signals of a tight relationship with hosts including a reduced functional capacity and genes transferred from fungal hosts to MRE. Phylogenetic reconstruction including nine MRE from mycorrhizal fungi revealed that MRE-B genomes are more closely related to MRE from Glomeromycotina than MRE-L from the same host family. We posit that reductions in genome size, GC content, pseudogene content, and repeat content in MRE-L may reflect a longer-term relationship with their fungal hosts. These data indicate Linnemannia and Benniella MRE were likely acquired independently after their fungal hosts diverged from a common ancestor. This work expands upon foundational knowledge on minimal genomes and provides insights into the evolution of bacterial endosymbionts.

The endohyphal microbiome: current progress and challenges for scaling down integrative multi-omic microbiome research.

As microbiome research has progressed, it has become clear that most, if not all, eukaryotic organisms are hosts to microbiomes composed of prokaryotes, other eukaryotes, and viruses. Fungi have only recently been considered holobionts with their own microbiomes, as filamentous fungi have been found to harbor bacteria (including cyanobacteria), mycoviruses, other fungi, and whole algal cells within their hyphae. Constituents of this complex endohyphal microbiome have been interrogated using multi-omic approaches. However, a lack of tools, techniques, and standardization for integrative multi-omics for small-scale microbiomes (e.g., intracellular microbiomes) has limited progress towards investigating and understanding the total diversity of the endohyphal microbiome and its functional impacts on fungal hosts. Understanding microbiome impacts on fungal hosts will advance explorations of how "microbiomes within microbiomes" affect broader microbial community dynamics and ecological functions. Progress to date as well as ongoing challenges of performing integrative multi-omics on the endohyphal microbiome is discussed herein. Addressing the challenges associated with the sample extraction, sample preparation, multi-omic data generation, and multi-omic data analysis and integration will help advance current knowledge of the endohyphal microbiome and provide a road map for shrinking microbiome investigations to smaller scales. Video Abstract.

Diversity of fungal microbiome obtained from plant rhizoplanes.

Microbial communities, and their ecological importance, have been investigated in several habitats. However, so far, most studies could not describe the closest microbial interactions and their functionalities. This study investigates the co-occurring interactions between fungi and bacteria in plant rhizoplanes and their potential functions. The partnerships were obtained using fungal-highway columns with four plant-based media. The fungi and associated microbiomes isolated from the columns were identified by sequencing the ITS (fungi) and 16S rRNA genes (bacteria). Statistical analyses including Exploratory Graph and Network Analysis were used to visualize the presence of underlying clusters in the microbial communities and evaluate the metabolic functions associated with the fungal microbiome (PICRUSt2). Our findings characterize the presence of both unique and complex bacterial communities associated with different fungi. The results showed that Bacillus was associated as exo-bacteria in 80 % of the fungi but occurred as putative endo-bacteria in 15 %. A shared core of putative endo-bacterial genera, potentially involved in the nitrogen cycle was found in 80 % of the isolated fungi. The comparison of potential metabolic functions of the putative endo- and exo-communities highlighted the potential essential factors to establish an endosymbiotic relationship, such as the loss of pathways associated with metabolites obtained from the host while maintaining pathways responsible for bacterial survival within the hypha.

Cryptosporulation in Kurthia spp. forces a rethinking of asporogenesis in Firmicutes.

Endosporulation is a complex morphophysiological process resulting in a more resistant cellular structure that is produced within the mother cell and is called endospore. Endosporulation evolved in the common ancestor of Firmicutes, but it is lost in descendant lineages classified as asporogenic. While Kurthia spp. is considered to comprise only asporogenic species, we show here that strain 11kri321, which was isolated from an oligotrophic geothermal reservoir, produces phase-bright spore-like structures. Phylogenomics of strain 11kri321 and other Kurthia strains reveals little similarity to genetic determinants of sporulation known from endosporulating Bacilli. However, morphological hallmarks of endosporulation were observed in two of the four Kurthia strains tested, resulting in spore-like structures (cryptospores). In contrast to classic endospores, these cryptospores did not protect against heat or UV damage and successive sub-culturing led to the loss of the cryptosporulating phenotype. Our findings imply that a cryptosporulation phenotype may have been prevalent and subsequently lost by laboratory culturing in other Firmicutes currently considered as asporogenic. Cryptosporulation might thus represent an ancestral but unstable and adaptive developmental state in Firmicutes that is under selection under harsh environmental conditions.

Intrahost SARS-CoV-2 k-mer Identification Method (iSKIM) for Rapid Detection of Mutations of Concern Reveals Emergence of Global Mutation Patterns.

Despite unprecedented global sequencing and surveillance of SARS-CoV-2, timely identification of the emergence and spread of novel variants of concern (VoCs) remains a challenge. Several million raw genome sequencing runs are now publicly available. We sought to survey these datasets for intrahost variation to study emerging mutations of concern. We developed iSKIM ("intrahost SARS-CoV-2 k-mer identification method") to relatively quickly and efficiently screen the many SARS-CoV-2 datasets to identify intrahost mutations belonging to lineages of concern. Certain mutations surged in frequency as intrahost minor variants just prior to, or while lineages of concern arose. The Spike N501Y change common to several VoCs was found as a minor variant in 834 samples as early as October 2020. This coincides with the timing of the first detected samples with this mutation in the Alpha/B.1.1.7 and Beta/B.1.351 lineages. Using iSKIM, we also found that Spike L452R was detected as an intrahost minor variant as early as September 2020, prior to the observed rise of the Epsilon/B.1.429/B.1.427 lineages in late 2020. iSKIM rapidly screens for mutations of interest in raw data, prior to genome assembly, and can be used to detect increases in intrahost variants, potentially providing an early indication of novel variant spread.

Design and construction of 3D printed devices to investigate active and passive bacterial dispersal on hydrated surfaces.

BACKGROUND: To disperse in water-unsaturated environments, such as the soil, bacteria rely on the availability and structure of water films forming on biotic and abiotic surfaces, and, especially, along fungal mycelia. Dispersal along such "fungal highways" may be driven both by mycelial physical properties and by interactions between bacteria and fungi. However, we still do not have a way to disentangle the biotic and abiotic elements. RESULTS: We designed and 3D printed two devices establishing stable liquid films that support bacteria dispersal in the absence of biotic interactions. The thickness of the liquid film determined the presence of hydraulic flow capable of transporting non-motile cells. In the absence of flow, only motile cells can disperse in the presence of an energy source. Non-motile cells could not disperse autonomously without flow but dispersed as "hitchhikers" when co-inoculated with motile cells. CONCLUSIONS: The 3D printed devices can be used as an abiotic control to study bacterial dispersal on hydrated surfaces, such as plant roots and fungal hyphae networks in the soil. By teasing apart the abiotic and biotic dimensions, these 3D printed devices will stimulate further research on microbial dispersal in soil and other water-unsaturated environments.

Complete Genome Sequences of the Soil Oxalotrophic Bacterium Cupriavidus oxalaticus Strain Ox1 and Its Derived mCherry-Tagged Strain.

Here, we report the complete genome sequences of the soil oxalotrophic bacterium Cupriavidus oxalaticus Ox1 and a derived mCherry-tagged strain. The genome size is approximately 6.69 Mb, with a GC content of 66.9%. The genome sequence of C. oxalaticus Ox1 contains a complete operon for the degradation and assimilation of oxalate.

Islands Within Islands: Bacterial Phylogenetic Structure and Consortia in Hawaiian Lava Caves and Fumaroles.

Lava caves, tubes, and fumaroles in Hawai'i present a range of volcanic, oligotrophic environments from different lava flows and host unexpectedly high levels of bacterial diversity. These features provide an opportunity to study the ecological drivers that structure bacterial community diversity and assemblies in volcanic ecosystems and compare the older, more stable environments of lava tubes, to the more variable and extreme conditions of younger, geothermally active caves and fumaroles. Using 16S rRNA amplicon-based sequencing methods, we investigated the phylogenetic distinctness and diversity and identified microbial interactions and consortia through co-occurrence networks in 70 samples from lava tubes, geothermal lava caves, and fumaroles on the island of Hawai'i. Our data illustrate that lava caves and geothermal sites harbor unique microbial communities, with very little overlap between caves or sites. We also found that older lava tubes (500-800 yrs old) hosted greater phylogenetic diversity (Faith's PD) than sites that were either geothermally active or younger (<400 yrs old). Geothermally active sites had a greater number of interactions and complexity than lava tubes. Average phylogenetic distinctness, a measure of the phylogenetic relatedness of a community, was higher than would be expected if communities were structured at random. This suggests that bacterial communities of Hawaiian volcanic environments are phylogenetically over-dispersed and that competitive exclusion is the main driver in structuring these communities. This was supported by network analyses that found that taxa (Class level) co-occurred with more distantly related organisms than close relatives, particularly in geothermal sites. Network "hubs" (taxa of potentially higher ecological importance) were not the most abundant taxa in either geothermal sites or lava tubes and were identified as unknown families or genera of the phyla, Chloroflexi and Acidobacteria. These results highlight the need for further study on the ecological role of microbes in caves through targeted culturing methods, metagenomics, and long-read sequence technologies.

Importance of appropriate genome information for the design of mating type primers in black and yellow morel populations.

Morels are highly prized edible fungi where sexual reproduction is essential for fruiting-body production. As a result, a comprehensive understanding of their sexual reproduction is of great interest. Central to this is the identification of the reproductive strategies used by morels. Sexual reproduction in fungi is controlled by mating-type (MAT) genes and morels are thought to be mainly heterothallic with two idiomorphs, MAT1-1 and MAT1-2. Genomic sequencing of black (Elata clade) and yellow (Esculenta clade) morel species has led to the development of PCR primers designed to amplify genes from the two idiomorphs for rapid genotyping of isolates from these two clades. To evaluate the design and theoretical performance of these primers we performed a thorough bioinformatic investigation, including the detection of the MAT region in publicly available Morchella genomes and in-silico PCR analyses. All examined genomes, including those used for primer design, appeared to be heterothallic. This indicates an inherent fault in the original primer design which utilized a single Morchella genome, as the use of two genomes with complementary mating types would be required to design accurate primers for both idiomorphs. Furthermore, potential off-targets were identified for some of the previously published primer sets, but verification was challenging due to lack of adequate genomic information and detailed methodologies for primer design. Examinations of the black morel specific primer pairs (MAT11L/R and MAT22L/R) indicated the MAT22 primers would correctly target and amplify the MAT1-2 idiomorph, but the MAT11 primers appear to be capable of amplifying incorrect off-targets within the genome. The yellow morel primer pairs (EMAT1-1 L/R and EMAT1-2 L/R) appear to have reporting errors, as the published primer sequences are dissimilar with reported amplicon sequences and the EMAT1-2 primers appear to amplify the RNA polymerase II subunit (RPB2) gene. The lack of the reference genome used in primer design and descriptive methodology made it challenging to fully assess the apparent issues with the primers for this clade. In conclusion, additional work is still required for the generation of reliable primers to investigate mating types in morels and to assess their performance on different clades and across multiple geographical regions.

Intrahost SARS-CoV-2 k-mer identification method (iSKIM) for rapid detection of mutations of concern reveals emergence of global mutation patterns.

Despite unprecedented global sequencing and surveillance of SARS-CoV-2, timely identification of the emergence and spread of novel variants of concern (VoCs) remains a challenge. Several million raw genome sequencing runs are now publicly available. We sought to survey these datasets for intrahost variation to study emerging mutations of concern. We developed iSKIM ("intrahost SARS-CoV-2 k-mer identification method") to relatively quickly and efficiently screen the many SARS-CoV-2 datasets to identify intrahost mutations belonging to lineages of concern. Certain mutations surged in frequency as intrahost minor variants just prior to, or while lineages of concern arose. The Spike N501Y change common to several VoCs was found as a minor variant in 834 samples as early as October 2020. This coincides with the timing of the first detected samples with this mutation in the Alpha/B.1.1.7 and Beta/B.1.351 lineages. Using iSKIM, we also found that Spike L452R was detected as an intrahost minor variant as early as September 2020, prior to the observed rise of the Epsilon/B.1.429/B.1.427 lineages in late 2020. iSKIM rapidly screens for mutations of interest in raw data, prior to genome assembly, and can be used to detect increases in intrahost variants, potentially providing an early indication of novel variant spread.

Advances and Challenges in Fluorescence in situ Hybridization for Visualizing Fungal Endobacteria.

Several bacteria have long been known to interact intimately with fungi, but molecular approaches have only recently uncovered how cosmopolitan these interactions are in nature. Currently, bacterial-fungal interactions (BFI) are inferred based on patterns of co-occurrence in amplicon sequencing investigations. However, determining the nature of these interactions, whether the bacteria are internally or externally associated, remains a grand challenge in BFI research. Fluorescence in situ hybridization (FISH) is a robust method that targets unique sequences of interest which can be employed for visualizing intra-hyphal targets, such as mitochondrial organelles or, as in this study, bacteria. We evaluate the challenges and employable strategies to resolve intra-hyphal BFI to address pertinent criteria in BFI research, such as culturing media, spatial distribution of bacteria, and abundance of bacterial 16S rRNA copies for fluorescent labeling. While these experimental factors influence labeling and detection of endobacteria, we demonstrate how to overcome these challenges thorough permeabilization, appropriate media choice, and targeted amplification using hybridization chain reaction FISH. Such microscopy imaging approaches can now be utilized by the broader research community to complement sequence-based investigations and provide more conclusive evidence on the nature of specific bacterial-fungal relationships.