Draft Genome Sequence of Aeromonas veronii Hm21, a Symbiotic Isolate from the Medicinal Leech Digestive Tract

Acquisition of Type I methyltransferase via horizontal gene transfer increases the drug resistance of Aeromonas veronii

Aeromonas veronii is an opportunistic pathogen that affects both fish and mammals, including humans, leading to bacteraemia, sepsis, meningitis and even death. The increasing virulence and drug resistance of A. veronii are of significant concern and pose a severe risk to public safety. The Type I restriction-modification (RM) system, which functions as a bacterial defence mechanism, can influence gene expression through DNA methylation. However, little research has been conducted to explore its origin, evolutionary path, and relationship to virulence and drug resistance in A. veronii. In this study, we analysed the pan-genome of 233 A. veronii strains, and the results indicated that it was 'open', meaning that A. veronii has acquired additional genes from other species. This suggested that A. veronii had the potential to adapt and evolve rapidly, which might have contributed to its drug resistance. One Type I methyltransferase (MTase) and two complete Type I RM systems were identified, namely AveC4I, AveC4II and AveC4III in A. veronii strain C4, respectively. Notably, AveC4I was exclusive to A. veronii C4. Phylogenetic analysis revealed that AveC4I was derived from horizontal gene transfer from Thiocystis violascens and exchanged genes with the human pathogen Comamonas kerstersii. Single molecule real-time sequencing was applied to identify the motif methylated by AveC4I, which was unique and not recognized by any reported MTases in the REBASE database. We also annotated the functions and pathways of the genes containing the motif, revealing that AveC4I may control drug resistance in A. veronii C4. Our findings provide new insight on the mechanisms underlying drug resistance in pathogenic bacteria. By identifying the specific genes and pathways affected by AveC4I, this study may aid in the development of new therapeutic approaches to combat A. veronii infections.

Polar Flagella Glycosylation in Aeromonas: Genomic Characterization and Involvement of a Specific Glycosyltransferase (Fgi-1) in Heterogeneous Flagella Glycosylation

Polar flagella from mesophilic Aeromonas strains have previously been shown to be modified with a range of glycans. Mass spectrometry studies of purified polar flagellins suggested the glycan typically includes a putative pseudaminic acid like derivative; while some strains are modified with this single monosaccharide, others modified with a heterologous glycan. In the current study, we demonstrate that genes involved in polar flagella glycosylation are clustered in highly polymorphic genomic islands flanked by pseudaminic acid biosynthetic genes (pse). Bioinformatic analysis of mesophilic Aeromonas genomes identified three types of polar flagella glycosylation islands (FGIs), denoted Group I, II and III. FGI Groups I and III are small genomic islands present in Aeromonas strains with flagellins modified with a single monosaccharide pseudaminic acid derivative. Group II were large genomic islands, present in strains found to modify polar flagellins with heterogeneous glycan moieties. Group II, in addition to pse genes, contained numerous glycosyltransferases and other biosynthetic enzymes. All Group II strains shared a common glycosyltransferase downstream of luxC that we named flagella glycosylation island 1, fgi-1, in A. piscicola AH-3. We demonstrate that Fgi-1 transfers the first sugar of the heterogeneous glycan to the pseudaminic acid derivative linked to polar flagellins and could be used as marker for polysaccharidic glycosylation of Aeromonas polar flagella.

Closed Genome Sequence of Aeromonas veronii Strain Hm21, an Isolate from the Medicinal Leech Hirudo verbana

Aeromonas veronii strain Hm21 was isolated from the medicinal leech and is used for genetic studies. We present here the 4.71-Mbp genome with a 56-kb plasmid and identify the mutations present in strains commonly used for genetic engineering.

Aeromonas veronii strain Hm21 was isolated from the medicinal leech Hirudo verbana and is used for genetic studies. We present here the 4.71-Mbp genome with a 56-kb plasmid and identify the mutations present in strains commonly used for genetic engineering.

Identification of Antimicrobial Resistance Determinants in Aeromonas veronii Strain MS-17-88 Recovered From Channel Catfish (Ictalurus punctatus)

Aeromonas veronii is a Gram-negative species ubiquitous in different aquatic environments and capable of causing a variety of diseases to a broad host range. Aeromonas species have the capability to carry and acquire antimicrobial resistance (AMR) elements, and currently multi-drug resistant (MDR) Aeromonas isolates are commonly found across the world. A. veronii strain MS-17-88 is a MDR strain isolated from catfish in the southeastern United States. The present study was undertaken to uncover the mechanism of resistance in MDR A. veronii strain MS-17-88 through the detection of genomic features. To achieve this, genomic DNA was extracted, sequenced, and assembled. The A. veronii strain MS-17-88 genome comprised 5,178,226-bp with 58.6% G+C, and it encoded several AMR elements, including imiS, ampS, mcr-7.1, mcr-3, catB2, catB7, catB1, floR, vat(F), tet(34), tet(35), tet(E), dfrA3, and tetR. The phylogeny and resistance profile of a large collection of A. veronii strains, including MS-17-88, were evaluated. Phylogenetic analysis showed a close relationship between MS-17-88 and strain Ae5 isolated from fish in China and ARB3 strain isolated from pond water in Japan, indicating a common ancestor of these strains. Analysis of phage elements revealed 58 intact, 63 incomplete, and 15 questionable phage elements among the 53 A. veronii genomes. The average phage element number is 2.56 per genome, and strain MS-17-88 is one of two strains having the maximum number of identified prophage elements (6 elements each). The profile of resistance against various antibiotics across the 53 A. veronii genomes revealed the presence of tet(34), mcr-7.1, mcr-3, and dfrA3 in all genomes (100%). By comparison, sul1 and sul2 were detected in 7.5% and 1.8% of A. veronii genomes. Nearly 77% of strains carried tet(E), and 7.5% of strains carried floR. This result suggested a low abundance and prevalence of sulfonamide and florfenicol resistance genes compared with tetracycline resistance among A. veronii strains. Overall, the present study provides insights into the resistance patterns among 53 A. veronii genomes, which can inform therapeutic options for fish affected by A. veronii.

Virulence characteristics and antimicrobial resistance of Aeromonas veronii biovar sobria 312M, a clinical isolate

Aeromonas are bacteria widely distributed in the environment, and some species are able to cause infections in humans, of which diarrhea is the most common. The objective of this study was to evaluate the presence of virulence and antimicrobial resistance associated characteristics in A. veronii biovar sobria strain 312M isolated from diarrheal stools. For this, the genome sequencing and phenotypical tests were performed. The draft genome annotation revealed several complete pathways associated with carbon metabolism and a mucin-desulfating sulfatase which may contribute to intestine colonization, and a large number of virulence-associated genes encoding structures associated with adhesion, toxins, and secretion systems. The strain exhibited swimming and swarming motility, biofilm formation, and hemolytic activity. It was resistant to ampicillin, ampicillin/sulbactam, and amoxicillin-clavulanic acid. Although a cphA gene encoding a narrow-spectrum carbapenase was identified in the strain genome, no carbapenemase activity was detected in the antimicrobial susceptibility test. When compared with other A. veronii with complete genomes, the main differences in virulence characteristics are related to lateral flagella and type III and VI secretion systems; the antimicrobial resistance spectrum also varied among strains. The results indicated that A. veronii biovar sobria 312M presents high virulence potential and resistance to limited classes of antimicrobials.

Comparative genomics of Aeromonas veronii: Identification of a pathotype impacting aquaculture globally

Aeromonas veronii is a gram-negative species abundant in aquatic environments that causes disease in humans as well as terrestrial and aquatic animals. In the current study, 41 publicly available A. veronii genomes were compared to investigate distribution of putative virulence genes, global dissemination of pathotypes, and potential mechanisms of virulence. The complete genome of A. veronii strain ML09-123 from an outbreak of motile aeromonas septicemia in farm-raised catfish in the southeastern United States was included. Dissemination of A. veronii strain types was discovered in dispersed geographical locations. Isolate ML09-123 is highly similar to Chinese isolate TH0426, suggesting the two strains have a common origin and may represent a pathotype impacting aquaculture in both countries. Virulence of strain ML09-123 in catfish in a dose-dependent manner was confirmed experimentally. Subsystem category disposition showed the majority of genomes exhibit similar distribution of genomic elements. The type I secretion system (T1SS), type II secretion system (T2SS), type 4 pilus (T4P), and flagellum core elements are conserved in all A. veronii genomes, whereas the type III secretion system (T3SS), type V secretion system (T5SS), type VI secretion system (T6SS), and tight adherence (TAD) system demonstrate variable dispersal. Distribution of mobile elements is dependent on host and geographic origin, suggesting this species has undergone considerable genetic exchange. The data presented here lends insight into the genomic variation of A. veronii and identifies a pathotype impacting aquaculture globally.

Identification and characterization of putative Aeromonas spp. T3SS effectors

The genetic determinants of bacterial pathogenicity are highly variable between species and strains. However, a factor that is commonly associated with virulent Gram-negative bacteria, including many Aeromonas spp., is the type 3 secretion system (T3SS), which is used to inject effector proteins into target eukaryotic cells. In this study, we developed a bioinformatics pipeline to identify T3SS effector proteins, applied this approach to the genomes of 105 Aeromonas strains isolated from environmental, mutualistic, or pathogenic contexts and evaluated the cytotoxicity of the identified effectors through their heterologous expression in yeast. The developed pipeline uses a two-step approach, where candidate Aeromonas gene families are initially selected using Hidden Markov Model (HMM) profile searches against the Virulence Factors DataBase (VFDB), followed by strict comparisons against positive and negative control datasets, greatly reducing the number of false positives. This approach identified 21 Aeromonas T3SS likely effector families, of which 8 represent known or characterized effectors, while the remaining 13 have not previously been described in Aeromonas. We experimentally validated our in silico findings by assessing the cytotoxicity of representative effectors in Saccharomyces cerevisiae BY4741, with 15 out of 21 assayed proteins eliciting a cytotoxic effect in yeast. The results of this study demonstrate the utility of our approach, combining a novel in silico search method with in vivo experimental validation, and will be useful in future research aimed at identifying and authenticating bacterial effector proteins from other genera.

Modernized Tools for Streamlined Genetic Manipulation and Comparative Study of Wild and Diverse Proteobacterial Lineages

A great challenge in microbiota research is the immense diversity of symbiotic bacteria with the capacity to impact the lives of plants and animals. Moving beyond correlative DNA sequencing-based studies to define the cellular and molecular mechanisms by which symbiotic bacteria influence the biology of their hosts is stalling because genetic manipulation of new and uncharacterized bacterial isolates remains slow and difficult with current genetic tools. Moreover, developing tools de novo is an arduous and time-consuming task and thus represents a significant barrier to progress. To address this problem, we developed a suite of engineering vectors that streamline conventional genetic techniques by improving postconjugation counterselection, modularity, and allelic exchange. Our modernized tools and step-by-step protocols will empower researchers to investigate the inner workings of both established and newly emerging models of bacterial symbiosis.

Correlating the presence of bacteria and the genes they carry with aspects of plant and animal biology is rapidly outpacing the functional characterization of naturally occurring symbioses. A major barrier to mechanistic studies is the lack of tools for the efficient genetic manipulation of wild and diverse bacterial isolates. To address the need for improved molecular tools, we used a collection of proteobacterial isolates native to the zebrafish intestinal microbiota as a testbed to construct a series of modernized vectors that expedite genetic knock-in and knockout procedures across lineages. The innovations that we introduce enhance the flexibility of conventional genetic techniques, making it easier to manipulate many different bacterial isolates with a single set of tools. We developed alternative strategies for domestication-free conjugation, designed plasmids with customizable features, and streamlined allelic exchange using visual markers of homologous recombination. We demonstrate the potential of these tools through a comparative study of bacterial behavior within the zebrafish intestine. Live imaging of fluorescently tagged isolates revealed a spectrum of distinct population structures that differ in their biogeography and dominant growth mode (i.e., planktonic versus aggregated). Most striking, we observed divergent genotype-phenotype relationships: several isolates that are predicted by genomic analysis and in vitro assays to be capable of flagellar motility do not display this trait within living hosts. Together, the tools generated in this work provide a new resource for the functional characterization of wild and diverse bacterial lineages that will help speed the research pipeline from sequencing-based correlations to mechanistic underpinnings.

A conserved bacterial protein induces pancreatic beta cell expansion during zebrafish development

Resident microbes play important roles in the development of the gastrointestinal tract, but their influence on other digestive organs is less well explored. Using the gnotobiotic zebrafish, we discovered that the normal expansion of the pancreatic β cell population during early larval development requires the intestinal microbiota and that specific bacterial members can restore normal β cell numbers. These bacteria share a gene that encodes a previously undescribed protein, named herein BefA (β Cell Expansion Factor A), which is sufficient to induce β cell proliferation in developing zebrafish larvae. Homologs of BefA are present in several human-associated bacterial species, and we show that they have conserved capacity to stimulate β cell proliferation in larval zebrafish. Our findings highlight a role for the microbiota in early pancreatic β cell development and suggest a possible basis for the association between low diversity childhood fecal microbiota and increased diabetes risk.

Host Matters: Medicinal Leech Digestive-Tract Symbionts and Their Pathogenic Potential

Digestive-tract microbiota exert tremendous influence over host health. Host-symbiont model systems are studied to investigate how symbioses are initiated and maintained, as well as to identify host processes affected by resident microbiota. The medicinal leech, Hirudo verbana, is an excellent model to address such questions owing to a microbiome that is consistently dominated by two species, Aeromonas veronii and Mucinivorans hirudinis, both of which are cultivable and have sequenced genomes. This review outlines current knowledge about the dynamics of the H. verbana microbiome. We discuss in depth the factors required for A. veronii colonization and proliferation in the leech crop and summarize the current understanding of interactions between A. veronii and its annelid host. Lastly, we discuss leech usage in modern medicine and highlight how leech-therapy associated infections, often attributable to Aeromonas spp., are of growing clinical concern due in part to an increased prevalence of fluoroquinolone resistant strains.

Draft Genome Sequence of Aeromonas sp. Strain EERV15

We report here the draft genome sequence of Aeromonas sp. strain EERV15 isolated from sand filter. The organism most closely related to Aeromonas sp. EERV15 is Aeromonas veronii B565, with an average 83% amino acid sequence similarity of putatively encoded protein open reading frames.

Lessons from Digestive-Tract Symbioses Between Bacteria and Invertebrates

In most animals, digestive tracts harbor the greatest number of bacteria in the animal that contribute to its health: by aiding in the digestion of nutrients, provisioning essential nutrients and protecting against colonization by pathogens. Invertebrates have been used to enhance our understanding of metabolic processes and microbe-host interactions owing to experimental advantages. This review describes how advances in DNA sequencing technologies have dramatically altered how researchers investigate microbe-host interactions, including 16S rRNA gene surveys, metagenome experiments, and metatranscriptome studies. Advantages and challenges of each of these approaches are described herein. Hypotheses generated through omics studies can be directly tested using site-directed mutagenesis, and findings from transposon studies and site-directed experiments are presented. Finally, unique structural aspects of invertebrate digestive tracts that contribute to symbiont specificity are presented. The combination of omics approaches with genetics and microscopy allows researchers to move beyond correlations to identify conserved mechanisms of microbe-host interactions.

A Tale of Transmission: Aeromonas veronii Activity within Leech-Exuded Mucus

Transmission, critical to the establishment and persistence of host-associated microbiotas, also exposes symbionts to new environmental conditions. With horizontal transmission, these different conditions represent major lifestyle shifts. Yet genome-wide analyses of how microbes adjust their transcriptomes toward these dramatic shifts remain understudied. Here, we provide a comprehensive and comparative analysis of the global transcriptional profiles of a symbiont as it shifts between lifestyles during transmission. The gammaproteobacterium Aeromonas veronii is transmitted from the gut of the medicinal leech to other hosts via host mucosal castings, yet A. veronii can also transition from mucosal habitancy to a free-living lifestyle. These three lifestyles are characterized by distinct physiological constraints and consequently lifestyle-specific changes in the expression of stress-response genes. Mucus-bound A. veronii had the greatest expression in terms of both the number of loci and levels of transcription of stress-response mechanisms. However, these bacteria are still capable of proliferating within the mucus, suggesting the availability of nutrients within this environment. We found that A. veronii alters transcription of loci in a synthetic pathway that obtains and incorporates N-acetylglucosamine (NAG; a major component of mucus) into the bacterial cell wall, enabling proliferation. Our results demonstrate that symbionts undergo dramatic local adaptation, demonstrated by widespread transcriptional changes, throughout the process of transmission that allows them to thrive while they encounter new environments which further shape their ecology and evolution.

Identification of Population Bottlenecks and Colonization Factors during Assembly of Bacterial Communities within the Zebrafish Intestine

The zebrafish, Danio rerio, is a powerful model for studying bacterial colonization of the vertebrate intestine, but the genes required by commensal bacteria to colonize the zebrafish gut have not yet been interrogated on a genome-wide level. Here we apply a high-throughput transposon mutagenesis screen to Aeromonas veronii Hm21 and Vibrio sp. strain ZWU0020 during their colonization of the zebrafish intestine alone and in competition with each other, as well as in different colonization orders. We use these transposon-tagged libraries to track bacterial population sizes in different colonization regimes and to identify gene functions required during these processes. We show that intraspecific, but not interspecific, competition with a previously established bacterial population greatly reduces the ability of these two bacterial species to colonize. Further, using a simple binomial sampling model, we show that under conditions of interspecific competition, genes required for colonization cannot be identified because of the population bottleneck experienced by the second colonizer. When bacteria colonize the intestine alone or at the same time as the other species, we find shared suites of functional requirements for colonization by the two species, including a prominent role for chemotaxis and motility, regardless of the presence of another species.

Zebrafish larvae, which are amenable to large-scale gnotobiotic studies, comprehensive sampling of their intestinal microbiota, and live imaging, are an excellent model for investigations of vertebrate intestinal colonization dynamics. We sought to develop a mutagenesis and tagging system in order to understand bacterial population dynamics and functional requirements during colonization of the larval zebrafish intestine. We explored changes in bacterial colonization dynamics and functional requirements when bacteria colonize a bacterium-free intestine, one previously colonized by their own species, or one colonized previously or simultaneously with a different species. This work provides a framework for rapid identification of colonization factors important under different colonization conditions. Furthermore, we demonstrate that when colonizing bacterial populations are very small, this approach is not accurate because random sampling of the input pool is sufficient to explain the distribution of inserts recovered from bacteria that colonized the intestines.

Identification of iron and heme utilization genes in Aeromonas and their role in the colonization of the leech digestive tract

It is known that many pathogens produce high-affinity iron uptake systems like siderophores and/or proteins for utilizing iron bound to heme-containing molecules, which facilitate iron-acquisition inside a host. In mutualistic digestive-tract associations, iron uptake systems have not been as well studied. We investigated the importance of two iron utilization systems within the beneficial digestive-tract association Aeromonas veronii and the medicinal leech, Hirudo verbana. Siderophores were detected in A. veronii using chrome azurol S. Using a mini Tn5, a transposon insertion in viuB generated a mutant unable to utilize iron using siderophores. The A. veronii genome was then searched for genes potentially involved in iron utilization bound to heme-containing molecules. A putative outer membrane heme receptor (hgpB) was identified with a transcriptional activator, termed hgpR, downstream. The hgpB gene was interrupted with an antibiotic resistance cassette in both the parent strain and the viuB mutant, yielding an hgpB mutant and a mutant with both iron uptake systems inactivated. In vitro assays indicated that hgpB is involved in utilizing iron bound to heme and that both iron utilization systems are important for A. veronii to grow in blood. In vivo colonization assays revealed that the ability to acquire iron from heme-containing molecules is critical for A. veronii to colonize the leech gut. Since iron and specifically heme utilization is important in this mutualistic relationship and has a potential role in virulence factor of other organisms, genomes from different Aeromonas strains (both clinical and environmental) were queried with iron utilization genes of A. veronii. This analysis revealed that in contrast to the siderophore utilization genes heme utilization genes are widely distributed among aeromonads. The importance of heme utilization in the colonization of the leech further confirms that symbiotic and pathogenic relationships possess similar mechanisms for interacting with animal hosts.

Spatial and temporal features of the growth of a bacterial species colonizing the zebrafish gut

The vertebrate intestine is home to microbial ecosystems that play key roles in host development and health. Little is known about the spatial and temporal dynamics of these microbial communities, limiting our understanding of fundamental properties, such as their mechanisms of growth, propagation, and persistence. To address this, we inoculated initially germ-free zebrafish larvae with fluorescently labeled strains of an Aeromonas species, representing an abundant genus in the zebrafish gut. Using light sheet fluorescence microscopy to obtain three-dimensional images spanning the gut, we quantified the entire bacterial load, as founding populations grew from tens to tens of thousands of cells over several hours. The data yield the first ever measurements of the growth kinetics of a microbial species inside a live vertebrate intestine and show dynamics that robustly fit a logistic growth model. Intriguingly, bacteria were nonuniformly distributed throughout the gut, and bacterial aggregates showed considerably higher growth rates than did discrete individuals. The form of aggregate growth indicates intrinsically higher division rates for clustered bacteria, rather than surface-mediated agglomeration onto clusters. Thus, the spatial organization of gut bacteria both relative to the host and to each other impacts overall growth kinetics, suggesting that spatial characterizations will be an important input to predictive models of host-associated microbial community assembly.

Our intestines are home to vast numbers of microbes that influence many aspects of health and disease. Though we now know a great deal about the constituents of the gut microbiota, we understand very little about their spatial structure and temporal dynamics in humans or in any animal: how microbial populations establish themselves, grow, fluctuate, and persist. To address this, we made use of a model organism, the zebrafish, and a new optical imaging technique, light sheet fluorescence microscopy, to visualize for the first time the colonization of a live, vertebrate gut by specific bacteria with sufficient resolution to quantify the population over a range from a few individuals to tens of thousands of bacterial cells. Our results provide unprecedented measures of bacterial growth kinetics and also show the influence of spatial structure on bacterial populations, which can be revealed only by direct imaging.

Bioinformatic genome comparisons for taxonomic and phylogenetic assignments using Aeromonas as a test case

Prokaryotic taxonomy is the underpinning of microbiology, as it provides a framework for the proper identification and naming of organisms. The “gold standard” of bacterial species delineation is the overall genome similarity determined by DNA-DNA hybridization (DDH), a technically rigorous yet sometimes variable method that may produce inconsistent results. Improvements in next-generation sequencing have resulted in an upsurge of bacterial genome sequences and bioinformatic tools that compare genomic data, such as average nucleotide identity (ANI), correlation of tetranucleotide frequencies, and the genome-to-genome distance calculator, or in silico DDH (isDDH). Here, we evaluate ANI and isDDH in combination with phylogenetic studies using Aeromonas, a taxonomically challenging genus with many described species and several strains that were reassigned to different species as a test case. We generated improved, high-quality draft genome sequences for 33 Aeromonas strains and combined them with 23 publicly available genomes. ANI and isDDH distances were determined and compared to phylogenies from multilocus sequence analysis of housekeeping genes, ribosomal proteins, and expanded core genes. The expanded core phylogenetic analysis suggested relationships between distant Aeromonas clades that were inconsistent with studies using fewer genes. ANI values of ≥96% and isDDH values of ≥70% consistently grouped genomes originating from strains of the same species together. Our study confirmed known misidentifications, validated the recent revisions in the nomenclature, and revealed that a number of genomes deposited in GenBank are misnamed. In addition, two strains were identified that may represent novel Aeromonas species.

Improvements in DNA sequencing technologies have resulted in the ability to generate large numbers of high-quality draft genomes and led to a dramatic increase in the number of publically available genomes. This has allowed researchers to characterize microorganisms using genome data. Advantages of genome sequence-based classification include data and computing programs that can be readily shared, facilitating the standardization of taxonomic methodology and resolving conflicting identifications by providing greater uniformity in an overall analysis. Using Aeromonas as a test case, we compared and validated different approaches. Based on our analyses, we recommend cutoff values for distance measures for identifying species. Accurate species classification is critical not only to obviate the perpetuation of errors in public databases but also to ensure the validity of inferences made on the relationships among species within a genus and proper identification in clinical and veterinary diagnostic laboratories.

Draft Genome Sequences of Amoeba-Resistant Aeromonas spp. Isolated from Aquatic Environments

Amoeba-resistant Aeromonas veronii ARB3 and Aeromonas media ARB13 and ARB20, which may be important intracellular pathogens of eukaryotic hosts, were isolated from pond and river waters. The draft genome sequences indicate that the strains harbor multiple protein secretion systems and toxins that induce disruption of the actin cytoskeleton.

Hitchhiking of host biology by beneficial symbionts enhances transmission

Transmission plays a key role in the evolution of symbiosis. Mixed mode transmission combines horizontal and vertical mechanisms for symbiont acquisition. However, features that enable mixed transmission are poorly understood. Here, we determine the mechanistic basis for the recruitment of the beneficial bacterium, Aeromonas veronii by the leech, Hirudo verbana. We demonstrate that host mucosal secretions complement imperfect symbiont vertical transmission. First, we show that the A. veronii population within secretions originates from the host digestive tract and proliferates synchronously with shedding frequency, demonstrating the coupling of partner biology. Furthermore, leeches are attracted to these castings with oral contact proving sufficient for symbiont transmission. Leech attraction to mucus is not affected by the symbiont state of either the host or mucus, suggesting that A. veronii exploits preexisting host behavior and physiological traits. A dual transmission mode, integrating multiple layers of host contributions, may prove evolutionarily advantageous for a wide range of symbioses. Using such a strategy, host infection is ensured, while also providing access to a higher genetic diversity of symbionts. Countless host-associated microbes exhibit mixed mode transmission, supporting the use of the leech symbiosis as a model for enhancing our understanding of the specificity, establishment and persistence of microbiotas.

Metagenomic analysis of the medicinal leech gut microbiota

There are trillions of microbes found throughout the human body and they exceed the number of eukaryotic cells by 10-fold. Metagenomic studies have revealed that the majority of these microbes are found within the gut, playing an important role in the host's digestion and nutrition. The complexity of the animal digestive tract, unculturable microbes, and the lack of genetic tools for most culturable microbes make it challenging to explore the nature of these microbial interactions within this niche. The medicinal leech, Hirudo verbana, has been shown to be a useful tool in overcoming these challenges, due to the simplicity of the microbiome and the availability of genetic tools for one of the two dominant gut symbionts, Aeromonas veronii. In this study, we utilize 16S rRNA gene pyrosequencing to further explore the microbial composition of the leech digestive tract, confirming the dominance of two taxa, the Rikenella-like bacterium and A. veronii. The deep sequencing approach revealed the presence of additional members of the microbial community that suggests the presence of a moderately complex microbial community with a richness of 36 taxa. The presence of a Proteus strain as a newly identified resident in the leech crop was confirmed using fluorescence in situ hybridization (FISH). The metagenome of this community was also pyrosequenced and the contigs were binned into the following taxonomic groups: Rikenella-like (3.1 MB), Aeromonas (4.5 MB), Proteus (2.9 MB), Clostridium (1.8 MB), Eryspelothrix (0.96 MB), Desulfovibrio (0.14 MB), and Fusobacterium (0.27 MB). Functional analyses on the leech gut symbionts were explored using the metagenomic data and MG-RAST. A comparison of the COG and KEGG categories of the leech gut metagenome to that of other animal digestive-tract microbiomes revealed that the leech digestive tract had a similar metabolic potential to the human digestive tract, supporting the usefulness of this system as a model for studying digestive-tract microbiomes. This study lays the foundation for more detailed metatranscriptomic studies and the investigation of symbiont population dynamics.