The holobiont transcriptome of teneral tsetse fly species of varying vector competence.

BACKGROUND: Tsetse flies are the obligate vectors of African trypanosomes, which cause Human and Animal African Trypanosomiasis. Teneral flies (newly eclosed adults) are especially susceptible to parasite establishment and development, yet our understanding of why remains fragmentary. The tsetse gut microbiome is dominated by two Gammaproteobacteria, an essential and ancient mutualist Wigglesworthia glossinidia and a commensal Sodalis glossinidius. Here, we characterize and compare the metatranscriptome of teneral Glossina morsitans to that of G. brevipalpis and describe unique immunological, physiological, and metabolic landscapes that may impact vector competence differences between these two species. RESULTS: An active expression profile was observed for Wigglesworthia immediately following host adult metamorphosis. Specifically, 'translation, ribosomal structure and biogenesis' followed by 'coenzyme transport and metabolism' were the most enriched clusters of orthologous genes (COGs), highlighting the importance of nutrient transport and metabolism even following host species diversification. Despite the significantly smaller Wigglesworthia genome more differentially expressed genes (DEGs) were identified between interspecific isolates (n = 326, ~ 55% of protein coding genes) than between the corresponding Sodalis isolates (n = 235, ~ 5% of protein coding genes) likely reflecting distinctions in host co-evolution and adaptation. DEGs between Sodalis isolates included genes involved in chitin degradation that may contribute towards trypanosome susceptibility by compromising the immunological protection provided by the peritrophic matrix. Lastly, G. brevipalpis tenerals demonstrate a more immunologically robust background with significant upregulation of IMD and melanization pathways. CONCLUSIONS: These transcriptomic differences may collectively contribute to vector competence differences between tsetse species and offers translational relevance towards the design of novel vector control strategies.

Different laboratory populations similar bacterial profile? The case of Glossina palpalis gambiensis.

BACKGROUND: Microbiota plays an important role in the biology, ecology and evolution of insects including tsetse flies. The bacterial profile of 3 Glossina palpalis gambiensis laboratory colonies was examined using 16S rRNA gene amplicon sequencing to evaluate the dynamics of the bacterial diversity within and between each G. p. gambiensis colony. RESULTS: The three G. p. gambiensis laboratory colonies displayed similar bacterial diversity indices and OTU distribution. Larval guts displayed a higher diversity when compared with the gastrointestinal tract of adults while no statistically significant differences were observed between testes and ovaries. Wigglesworthia and Sodalis were the most dominant taxa. In more detail, the gastrointestinal tract of adults was more enriched by Wigglesworthia while Sodalis were prominent in gonads. Interestingly, in larval guts a balanced co-existence between Wigglesworthia and Sodalis was observed. Sequences assigned to Wolbachia, Propionibacterium, and Providencia were also detected but to a much lesser degree. Clustering analysis indicated that the bacterial profile in G. p. gambiensis exhibits tissue tropism, hence distinguishing the gut bacterial profile from that present in reproductive organs. CONCLUSIONS: Our results indicated that age, gender and the origin of the laboratory colonies did not significantly influence the formation of the bacterial profile, once these populations were kept under the same rearing conditions. Within the laboratory populations a tissue tropism was observed between the gut and gonadal bacterial profile.

Into the Wild: Parallel Transcriptomics of the Tsetse-Wigglesworthia Mutualism within Kenyan Populations.

Tsetse flies (Diptera: Glossinidae) have medical significance as the obligate vectors of African trypanosomes. In addition, tsetse harbor a simple gut microbiota. A predominant gut microbiota member, the Gammaproteobacterium Wigglesworthia spp., has coevolved with tsetse for a significant portion of Glossina radiation proving critical to tsetse fitness. Although multiple roles have been described for Wigglesworthia within colony flies, little research has been dedicated towards functional characterization within wild tsetse. Here, dual RNA-Seq was performed to characterize the tsetse-Wigglesworthia symbiosis within flies captured in Nguruman, Kenya. A significant correlation in Gene Ontology (GO) distribution between tsetse and Wigglesworthia was observed, with homogeneous enrichment in metabolic and transport categories, likely supporting a hallmark of the symbiosis-bidirectional metabolic exchange. Within field flies, highly transcribed Wigglesworthia loci included those involved in B vitamin synthesis and in substrate translocation, including amino acid transporters and multidrug efflux pumps, providing a molecular means for interaction. The universal expression of several Wigglesworthia and G. pallidipes orthologs, putatively involved in nutrient provisioning and resource allocation, was confirmed in sister tsetse species. These transcriptional profiles varied through host age and mating status likely addressing varying symbiont demands and also confirming their global importance within Glossina. This study, not only supports symbiont nutrient provisioning roles, but also serves as a foundation for insight into novel roles and molecular mechanisms associated with vector-microbiota interactions. The role of symbiont B vitamin provisioning towards impacting host epigenetics is discussed. Knowledge of vector-microbiota interactions may lead to the discovery of novel targets in pest control.

Challenging the Wigglesworthia, Sodalis, Wolbachia symbiosis dogma in tsetse flies: Spiroplasma is present in both laboratory and natural populations.

Profiling of wild and laboratory tsetse populations using 16S rRNA gene amplicon sequencing allowed us to examine whether the "Wigglesworthia-Sodalis-Wolbachia dogma" operates across species and populations. The most abundant taxa, in wild and laboratory populations, were Wigglesworthia (the primary endosymbiont), Sodalis and Wolbachia as previously characterized. The species richness of the microbiota was greater in wild than laboratory populations. Spiroplasma was identified as a new symbiont exclusively in Glossina fuscipes fuscipes and G. tachinoides, members of the palpalis sub-group, and the infection prevalence in several laboratory and natural populations was surveyed. Multi locus sequencing typing (MLST) analysis identified two strains of tsetse-associated Spiroplasma, present in G. f. fuscipes and G. tachinoides. Spiroplasma density in G. f. fuscipes larva guts was significantly higher than in guts from teneral and 15-day old male and female adults. In gonads of teneral and 15-day old insects, Spiroplasma density was higher in testes than ovaries, and was significantly higher density in live versus prematurely deceased females indicating a potentially mutualistic association. Higher Spiroplasma density in testes than in ovaries was also detected by fluorescent in situ hybridization in G. f. fuscipes.

Analysis of multiple tsetse fly populations in Uganda reveals limited diversity and species-specific gut microbiota.

The invertebrate microbiome contributes to multiple aspects of host physiology, including nutrient supplementation and immune maturation processes. We identified and compared gut microbial abundance and diversity in natural tsetse flies from Uganda using five genetically distinct populations of Glossina fuscipes fuscipes and multiple tsetse species (Glossina morsitans morsitans, G. f. fuscipes, and Glossina pallidipes) that occur in sympatry in one location. We used multiple approaches, including deep sequencing of the V4 hypervariable region of the 16S rRNA gene, 16S rRNA gene clone libraries, and bacterium-specific quantitative PCR (qPCR), to investigate the levels and patterns of gut microbial diversity from a total of 151 individuals. Our results show extremely limited diversity in field flies of different tsetse species. The obligate endosymbiont Wigglesworthia dominated all samples (>99%), but we also observed wide prevalence of low-density Sodalis (tsetse's commensal endosymbiont) infections (<0.05%). There were also several individuals (22%) with high Sodalis density, which also carried coinfections with Serratia. Albeit in low density, we noted differences in microbiota composition among the genetically distinct G. f. fuscipes flies and between different sympatric species. Interestingly, Wigglesworthia density varied in different species (10(4) to 10(6) normalized genomes), with G. f. fuscipes having the highest levels. We describe the factors that may be responsible for the reduced diversity of tsetse's gut microbiota compared to those of other insects. Additionally, we discuss the implications of Wigglesworthia and Sodalis density variations as they relate to trypanosome transmission dynamics and vector competence variations associated with different tsetse species.

Population dynamics of Glossina palpalis gambiensis symbionts, Sodalis glossinidius, and Wigglesworthia glossinidia, throughout host-fly development.

The tsetse fly (Diptera: Glossinidae), the vector of trypanosomes causing human and animal trypanosomiasis, harbors symbiotic microorganisms including the primary symbiont Wigglesworthia glossinidia, involved in the fly's nutrition and fertility, and the secondary symbiont Sodalis glossinidius, involved in the trypanosome establishment in the fly's midgut. Both symbionts are maternally transmitted to the intrauterine progeny through the fly's milk gland secretions. In this study, we investigated the population dynamics of these symbionts during fly development. Wigglesworthia and Sodalis densities were estimated using quantitative PCR performed on Glossina palpalis gambiensis at different developmental stages. The results showed that the density of the primary Wigglesworthia symbiont was higher than that of Sodalis for all host developmental stages. Sodalis densities remained constant in pupae, but increased significantly in adult flies. The opposite situation was observed for Wigglesworthia, whose density increased in pupae and remained constant during the female adult stage. Moreover, Wigglesworthia density increased significantly during the transition from the pupal to the teneral stage, while mating had a contradictory effect depending on the age of the fly. Finally, tsetse fly colonization by both symbionts appears as a continuous and adaptive process throughout the insect's development. Last, the study demonstrated both symbionts of G. p. gambiensis, the vector of the chronic form of human African trypanosomiasis, to be permanent inhabitants of the colony flies throughout their life span. This was expected for the primary symbiont, Wigglesworthia, but not necessarily for the secondary symbiont, S. glossinidius, whose permanent presence is not required for the fly's survival. This result is of importance as Sodalis could be involved in the tsetse fly vector competence and may constitute a target in the frame of sleeping sickness fighting strategies.

The tsetse fly obligate mutualist Wigglesworthia morsitans alters gene expression and population density via exogenous nutrient provisioning.

The obligate mutualist Wigglesworthia morsitans provisions nutrients to tsetse flies. The symbiont's response to thiamine (B(1)) supplementation of blood meals, specifically towards the regulation of thiamine biosynthesis and population density, is described. Despite an ancient symbiosis and associated genome tailoring, Wigglesworthia responds to nutrient availability, potentially accommodating a decreased need.

Insight into the transmission biology and species-specific functional capabilities of tsetse (Diptera: glossinidae) obligate symbiont Wigglesworthia.

UNLABELLED: Ancient endosymbionts have been associated with extreme genome structural stability with little differentiation in gene inventory between sister species. Tsetse flies (Diptera: Glossinidae) harbor an obligate endosymbiont, Wigglesworthia, which has coevolved with the Glossina radiation. We report on the ~720-kb Wigglesworthia genome and its associated plasmid from Glossina morsitans morsitans and compare them to those of the symbiont from Glossina brevipalpis. While there was overall high synteny between the two genomes, a large inversion was noted. Furthermore, symbiont transcriptional analyses demonstrated host tissue and development-specific gene expression supporting robust transcriptional regulation in Wigglesworthia, an unprecedented observation in other obligate mutualist endosymbionts. Expression and immunohistochemistry confirmed the role of flagella during the vertical transmission process from mother to intrauterine progeny. The expression of nutrient provisioning genes (thiC and hemH) suggests that Wigglesworthia may function in dietary supplementation tailored toward host development. Furthermore, despite extensive conservation, unique genes were identified within both symbiont genomes that may result in distinct metabolomes impacting host physiology. One of these differences involves the chorismate, phenylalanine, and folate biosynthetic pathways, which are uniquely present in Wigglesworthia morsitans. Interestingly, African trypanosomes are auxotrophs for phenylalanine and folate and salvage both exogenously. It is possible that W. morsitans contributes to the higher parasite susceptibility of its host species. IMPORTANCE: Genomic stasis has historically been associated with obligate endosymbionts and their sister species. Here we characterize the Wigglesworthia genome of the tsetse fly species Glossina morsitans and compare it to its sister genome within G. brevipalpis. The similarity and variation between the genomes enabled specific hypotheses regarding functional biology. Expression analyses indicate significant levels of transcriptional regulation and support development- and tissue-specific functional roles for the symbiosis previously not observed in obligate mutualist symbionts. Retention of the genetically expensive flagella within these small genomes was demonstrated to be significant in symbiont transmission and tailored to the unique tsetse fly reproductive biology. Distinctions in metabolomes were also observed. We speculate an additional role for Wigglesworthia symbiosis where infections with pathogenic trypanosomes may depend upon symbiont species-specific metabolic products and thus influence the vector competence traits of different tsetse fly host species.

Multiple origins of endosymbiosis within the Enterobacteriaceae (γ-Proteobacteria): convergence of complex phylogenetic approaches.

BACKGROUND: The bacterial family Enterobacteriaceae gave rise to a variety of symbiotic forms, from the loosely associated commensals, often designated as secondary (S) symbionts, to obligate mutualists, called primary (P) symbionts. Determination of the evolutionary processes behind this phenomenon has long been hampered by the unreliability of phylogenetic reconstructions within this group of bacteria. The main reasons have been the absence of sufficient data, the highly derived nature of the symbiont genomes and lack of appropriate phylogenetic methods. Due to the extremely aberrant nature of their DNA, the symbiotic lineages within Enterobacteriaceae form long branches and tend to cluster as a monophyletic group. This state of phylogenetic uncertainty is now improving with an increasing number of complete bacterial genomes and development of new methods. In this study, we address the monophyly versus polyphyly of enterobacterial symbionts by exploring a multigene matrix within a complex phylogenetic framework. RESULTS: We assembled the richest taxon sampling of Enterobacteriaceae to date (50 taxa, 69 orthologous genes with no missing data) and analyzed both nucleic and amino acid data sets using several probabilistic methods. We particularly focused on the long-branch attraction-reducing methods, such as a nucleotide and amino acid data recoding and exclusion (including our new approach and slow-fast analysis), taxa exclusion and usage of complex evolutionary models, such as nonhomogeneous model and models accounting for site-specific features of protein evolution (CAT and CAT+GTR). Our data strongly suggest independent origins of four symbiotic clusters; the first is formed by Hamiltonella and Regiella (S-symbionts) placed as a sister clade to Yersinia, the second comprises Arsenophonus and Riesia (S- and P-symbionts) as a sister clade to Proteus, the third Sodalis, Baumannia, Blochmannia and Wigglesworthia (S- and P-symbionts) as a sister or paraphyletic clade to the Pectobacterium and Dickeya clade and, finally, Buchnera species and Ishikawaella (P-symbionts) clustering with the Erwinia and Pantoea clade. CONCLUSIONS: The results of this study confirm the efficiency of several artifact-reducing methods and strongly point towards the polyphyly of P-symbionts within Enterobacteriaceae. Interestingly, the model species of symbiotic bacteria research, Buchnera and Wigglesworthia, originated from closely related, but different, ancestors. The possible origins of intracellular symbiotic bacteria from gut-associated or pathogenic bacteria are suggested, as well as the role of facultative secondary symbionts as a source of bacteria that can gradually become obligate maternally transferred symbionts.

Influence of host phylogeographic patterns and incomplete lineage sorting on within-species genetic variability in Wigglesworthia species, obligate symbionts of tsetse flies.

Vertical transmission of obligate symbionts generates a predictable evolutionary history of symbionts that reflects that of their hosts. In insects, evolutionary associations between symbionts and their hosts have been investigated primarily among species, leaving population-level processes largely unknown. In this study, we investigated the tsetse (Diptera: Glossinidae) bacterial symbiont, Wigglesworthia glossinidia, to determine whether observed codiversification of symbiont and tsetse host species extends to a single host species (Glossina fuscipes fuscipes) in Uganda. To explore symbiont genetic variation in G. f. fuscipes populations, we screened two variable loci (lon and lepA) from the Wigglesworthia glossinidia bacterium in the host species Glossina fuscipes fuscipes (W. g. fuscipes) and examined phylogeographic and demographic characteristics in multiple host populations. Symbiont genetic variation was apparent within and among populations. We identified two distinct symbiont lineages, in northern and southern Uganda. Incongruence length difference (ILD) tests indicated that the two lineages corresponded exactly to northern and southern G. f. fuscipes mitochondrial DNA (mtDNA) haplogroups (P = 1.0). Analysis of molecular variance (AMOVA) confirmed that most variation was partitioned between the northern and southern lineages defined by host mtDNA (85.44%). However, ILD tests rejected finer-scale congruence within the northern and southern populations (P = 0.009). This incongruence was potentially due to incomplete lineage sorting that resulted in novel combinations of symbiont genetic variants and host background. Identifying these novel combinations may have public health significance, since tsetse is the sole vector of sleeping sickness and Wigglesworthia is known to influence host vector competence. Thus, understanding the adaptive value of these host-symbiont combinations may afford opportunities to develop vector control methods.

Reductive evolution of bacterial genome in insect gut environment.

Obligate endocellular symbiotic bacteria of insects and other organisms generally exhibit drastic genome reduction. Recently, it was shown that symbiotic gut bacteria of some stinkbugs also have remarkably reduced genomes. Here, we report the complete genome sequence of such a gut bacterium Ishikawaella capsulata of the plataspid stinkbug Megacopta punctatissima. Gene repertoire and evolutionary patterns, including AT richness and elevated evolutionary rate, of the 745,590 bp genome were strikingly similar to those of obligate γ-proteobacterial endocellular insect symbionts like Buchnera in aphids and Wigglesworthia in tsetse flies. Ishikawaella was suggested to supply essential amino acids for the plant-sucking stinkbug as Buchnera does for the host aphid. Although Buchnera is phylogenetically closer to Wigglesworthia than to Ishikawaella, in terms of gene repertoire Buchnera was similar to Ishikawaella rather than to Wigglesworthia, providing a possible case of genome-level convergence of gene content. Meanwhile, several notable differences were identified between the genomes of Ishikawaella and Buchnera, including retention of TCA cycle genes and lack of flagellum-related genes in Ishikawaella, which may reflect their adaptation to distinct symbiotic habitats. Unexpectedly, Ishikawaella retained fewer genes related to cell wall synthesis and lipid metabolism than many endocellular insect symbionts. The plasmid of Ishikawaella encoded genes for arginine metabolism and oxalate detoxification, suggesting the possibility of additional Ishikawaella roles similar to those of human gut bacteria. Our data highlight strikingly similar evolutionary patterns that are shared between the extracellular and endocellular insect symbiont genomes.

The tsetse fly Glossina fuscipes fuscipes (Diptera: Glossina) harbours a surprising diversity of bacteria other than symbionts.

Three different bacterial species are regularly described from tsetse flies. However, no broad screens have been performed to investigate the existence of other bacteria in this medically and agriculturally important vector insect. Utilising both culture dependent and independent methods we show that Kenyan populations of Glossina fuscipes fuscipes harbour a surprising diversity of bacteria. Bacteria were isolated from 72% of flies with 23 different bacterial species identified. The Firmicutes phylum dominated with 16 species of which seven belong to the genus Bacillus. The tsetse fly primary symbiont, Wigglesworthia glossinidia, was identified by the culture independent pathway. However, neither the secondary symbiont Sodalis nor Wolbachia was detected with either of the methods used. Two other bacterial species were identified with the DNA based method, Bacillus subtilis and Serratia marcescens. Further studies are needed to determine how tsetse flies, which only ever feed on vertebrate blood, pick up bacteria and to investigate the possible impact of these bacteria on Glossina longevity and vector competence.

Chance and necessity in the evolution of minimal metabolic networks.

It is possible to infer aspects of an organism's lifestyle from its gene content. Can the reverse also be done? Here we consider this issue by modelling evolution of the reduced genomes of endosymbiotic bacteria. The diversity of gene content in these bacteria may reflect both variation in selective forces and contingency-dependent loss of alternative pathways. Using an in silico representation of the metabolic network of Escherichia coli, we examine the role of contingency by repeatedly simulating the successive loss of genes while controlling for the environment. The minimal networks that result are variable in both gene content and number. Partially different metabolisms can thus evolve owing to contingency alone. The simulation outcomes do preserve a core metabolism, however, which is over-represented in strict intracellular bacteria. Moreover, differences between minimal networks based on lifestyle are predictable: by simulating their respective environmental conditions, we can model evolution of the gene content in Buchnera aphidicola and Wigglesworthia glossinidia with over 80% accuracy. We conclude that, at least for the particular cases considered here, gene content of an organism can be predicted with knowledge of its distant ancestors and its current lifestyle.

Gene expression levels influence amino acid usage and evolutionary rates in endosymbiotic bacteria.

Most endosymbiotic bacteria have extremely reduced genomes, accelerated evolutionary rates, and strong AT base compositional bias thought to reflect reduced efficacy of selection and increased mutational pressure. Here, we present a comparative study of evolutionary forces shaping five fully sequenced bacterial endosymbionts of insects. The results of this study were three-fold: (i) Stronger conservation of high expression genes at not just nonsynonymous, but also synonymous, sites. (ii) Variation in amino acid usage strongly correlates with GC content and expression level of genes. This pattern is largely explained by greater conservation of high expression genes, leading to their higher GC content. However, we also found indication of selection favoring GC-rich amino acids that contrasts with former studies. (iii) Although the specific nutritional requirements of the insect host are known to affect gene content of endosymbionts, we found no detectable influence on substitution rates, amino acid usage, or codon usage of bacterial genes involved in host nutrition.

The relationship between palindrome avoidance and intragenic codon usage variations: a Monte Carlo study.

Several studies have shown that codon usage within genes varies, as it seems dependent on both codon context and codon position within the gene. Given that palindromes in addition often are avoided in genomes, this study aimed at finding out if intragenic variations in codon usage may be a way to control the amount and location of palindromes. A Monte Carlo algorithm was written which resampled the codons in genes while keeping the amino acid sequence of the translation product constant. On the resampled sequences, palindromes were counted and their intragenic positions mapped. Escherichia coli K12 uses type II restriction-modification systems and displays pronounced codon usage phenomena. Using this as a reference organism it was clearly shown that the number of palindromes in genes is generally lower than the amount of palindromes in resampled genes; thus, the succession of codons seems to be a way to decrease the number of palindromes. The intragenic position of palindromes in resampled sequences, however, was largely equal to the position in the native genes, so codon usage phenomena are unlikely to be a way to control the intragenic position of palindromes. The analysis was repeated on two bacteriophages and gave similar same results, even though the virus genomes are much smaller. Studies on the endosymbionts Buchnera sp. APS and Wigglesworthia sp., which seemingly have no type II restriction-modification systems, showed that in these species there is only weak evidence for codon usage acting to control the number of palindromes.

Gene expression level influences amino acid usage, but not codon usage, in the tsetse fly endosymbiont Wigglesworthia.

Wigglesworthia glossinidia brevipalpis, the obligate bacterial endosymbiont of the tsetse fly Glossina brevipalpis, is characterized by extreme genome reduction and AT nucleotide composition bias. Here, multivariate statistical analyses are used to test the hypothesis that mutational bias and genetic drift shape synonymous codon usage and amino acid usage of Wigglesworthia. The results show that synonymous codon usage patterns vary little across the genome and do not distinguish genes of putative high and low expression levels, thus indicating a lack of translational selection. Extreme AT composition bias across the genome also drives relative amino acid usage, but predicted high-expression genes (ribosomal proteins and chaperonins) use GC-rich amino acids more frequently than do low-expression genes. The levels and configuration of amino acid differences between Wigglesworthia and Escherichia coli were compared to test the hypothesis that the relatively GC-rich amino acid profiles of high-expression genes reflect greater amino acid conservation at these loci. This hypothesis is supported by reduced levels of protein divergence at predicted high-expression Wigglesworthia genes and similar configurations of amino acid changes across expression categories. Combined, the results suggest that codon and amino acid usage in the Wigglesworthia genome reflect a strong AT mutational bias and elevated levels of genetic drift, consistent with expected effects of an endosymbiotic lifestyle and repeated population bottlenecks. However, these impacts of mutation and drift are apparently attenuated by selection on amino acid composition at high-expression genes.

Horizontal gene transfer of "prototype" Nramp in bacteria.

Eukaryotic Nramp genes encode divalent metal ion permeases important for nutrition and resistance to microbial infection. Bacterial homologs encode proton-dependent transporters of manganese (MntH), and other divalent metal ions. Bacterial MntH were classified in three homology groups (A, B, C) and MntH C further subdivided in Calpha, Cbeta, Cgamma. The proteins from C. tepidum (MntH B) and E. faecalis (MntH Cbeta1, 2), divergent in sequence and hydropathy profile, conferred increased metal sensitivity when expressed in E. coli, suggesting conservation of divalent metal transport function in MntH B and C. Several genomic evidence suggest horizontal gene transfer (HGT) of mntH C genes: (i) The enterobacteria Wigglesworthia mntH Cbeta gene is linked to an Asn t-RNA, and its sequence most conserved with Gram positive bacteria homologs; (ii) all the Cbeta genes identified in oral streptococcaceae are associated with different potentially mobile DNA elements; (iii) Lactococcus lactis and Burkholderia mallei genomes contain an mntH gene prematurely terminated and a novel full-length mntH C gene; (iv) remarkable sequence relatedness between the unicellular alga C. reinhardtii "prototype" Nramp and some MntH Calpha (e.g., Nostoc spp., Listeria spp.) suggests HGT between Eukarya and Bacteria. Other "prototype" Nramp genes (intronless, encoding proteins strongly conserved with MntH A and B proteins) identified in invertebrates represent a possible source for transfer of Nramp genes toward opportunistic bacteria. This study demonstrates complex evolution of MntH in Bacteria. It is proposed that "prototype" Nramp are ancestors of bacterial MntH C proteins, which could facilitate bacterial infection.