Global Distribution and Evolution of a Toxinogenic Burkholderia-Rhizopus Symbiosis

Toxinogenic endobacteria were isolated from a collection of Rhizopus spp. representing highly diverse geographic origins and ecological niches. All endosymbionts belonged to the Burkholderia rhizoxinica complex according to matrix-assisted laser desorption ionization-time of flight biotyping and multilocus sequence typing, suggesting a common ancestor. Comparison of host and symbiont phylogenies provides insights into possible cospeciation and horizontal-transmission events.Bacterial symbionts and their metabolic potential play essential roles for many organisms. They may benefit from improved fitness, survival, and even acquired virulence (7, 12, 22). In the course of our studies of the biosynthesis of rhizoxin, the causative agent of rice seedling blight (10), we found that the phytotoxin is produced not by the fungus Rhizopus microsporus but by symbiotic bacteria (Burkholderia rhizoxinica) that reside within the fungus cytosol (13, 15, 23). Furthermore, cloning and sequencing of the rhizoxin biosynthesis gene cluster revealed the molecular basis of bacterial toxin production (14). In sum, this represents an unparalleled example for a symbiosis in which a fungus harbors bacteria for the production of a virulence factor. In analogy, we found that the first reported “mycotoxins” from lower fungi, the highly toxic cyclopeptides rhizonin A and B (25, 28), are also produced by symbiotic bacteria (Burkholderia endofungorum) and not by the fungus (16). While both rhizoxins and rhizonins have been believed to promote zygomycoses (21), there is no indication for toxin-producing endosymbiotic bacteria in clinical isolates (18).In nature, toxin production plays a pivotal role in the development of the fungus-bacterium association. Studies of the evolution of host resistance indicate that the association resulted from a pathogenicity mutualism shift in insensitive zygomycetes (24). The fungus lost its ability to sporulate independently and became totally dependent on endobacteria for reproduction through spores, thus warranting the persistence of the symbiosis and its efficient distribution through vegetative spores (17).To gain a broader view of the occurrence, biosynthetic potential, and relationship of toxinogenic endofungal bacteria, we investigated a collection of Rhizopus spp. consisting of 20 isolates classified as R. microsporus (of which 13 belong to R. microsporus var. microsporus, four to R. microsporus var. chinensis, two to R. microsporus var. oligosporus, and one to R. microsporus var. rhizopodiformis), one isolate classified as Rhizopus sp., and one Rhizopus oryzae strain. We initially monitored the presence of bacterial symbionts by PCR using universal primers (16S rRNA genes) and rhizoxin production in all available Rhizopus strains. Liquid cultivation of fungi in production medium with and without antibiotic followed by organic solvent extraction yielded crude extracts that were analyzed by high-performance liquid chromatography (HPLC) and mass spectrometry (MS). In total, eight fungal strains were identified or confirmed as rhizoxin positive and thus expected to harbor endosymbionts. In all cases, this assumption was verified by PCR and confocal scanning microscopy. By means of an optimized protocol, we finally succeeded in the isolation and cultivation of all eight bacterial symbiont strains in pure cultures (isolates B1 to B8) (Table ​(Table11).

TABLE 1.

Fungal strains and their bacterial endosymbionts

Taxon	Strain designationa	Origin	Bacterial endosymbiont (isolate)
Rhizopus microsporus van Tieghem	ATCC 62417	Rice seedlings, Japan	Burkholderia rhizoxinica HKI-0454 (B1)
Rhizopus sp. strain F-1360	ATCC 20577	Soil, Japan	Burkholderia sp. strain HKI-0512 (B2)
Rhizopus microsporus Tieghem var. microsporus	CBS 111563	Sufu starter culture, rice wine tablet, Vietnam	Burkholderia sp. strain HKI-0455 (B3)
Rhizopus microsporus Tieghem var. microsporus	CBS 699.68	Soil, Ukraine	Burkholderia sp. strain HKI-402 (B4)
Rhizopus microsporus Tieghem	CBS 112285	Ground nuts, Mozambique	Burkholderia endofungorum HKI-0456 (B5)
Rhizopus microsporus var. chinensis (Saito) Schipper & Stalpers	CBS 261.28	Not specified, United States of America	Burkholderia sp. strain HKI-0513 (B6)
Rhizopus microsporus Tieghem var. microsporus	CBS 700.68	Forest soil, Georgia	Burkholderia sp. strain HKI-0403 (B7)
Rhizopus microsporus Tieghem var. microsporus	CBS 308.87	Man, from deep necrotic tissue within the hand following a spider bite, Australia	Burkholderia sp. strain HKI-0404 (B8)

Open in a separate window^aATCC, American Type Culture Collection, Manassas, VA; CBS, Centraalbureau voor Schimmelcultures, Utrecht, The Netherlands.Notably, the eight Rhizopus isolates are from geographically distinct collection sites, covering all five continents (Africa, America, Asia, Australia, and Europe) and representing diverse ecological niches of the host (plants, soil, food, and necrotic tissue) (Fig. ​(Fig.1;1; Table ​Table1).1). HPLC and MS analyses of the metabolic profiles and comparison with authentic references revealed that all endofungal bacterial strains are capable of producing considerable amounts of rhizoxin derivatives 1 to 7 (23) (Fig. ​(Fig.2).2). Among the rhizoxin derivatives, rhizoxin S2 (derivative 3) is the main product formed by all isolates, followed by compounds WF-1360F (derivative 6) (11), rhizoxin Z1 (derivative 4), and rhizoxin S1 (derivative 1) (23), while derivatives 2, 5, and 7 are formed only in minor amounts. Significant differences in production of rhizoxins were not found among the isolates (see Fig. S1 in the supplemental material). Only one isolate, Burkholderia endofungorum HKI-0456 (isolate B5), also produces the hepatotoxic cyclopeptides rhizonin A (derivative 8) and B (derivative 9) under laboratory conditions (16).Open in a separate windowFIG. 1.Survey of collection sites of toxinogenic R. microsporus strains used in this study.Open in a separate windowFIG. 2.Structures of the main rhizoxin derivatives (derivatives 1 to 7) produced by all eight fungal endosymbionts (isolates B1 to B8) and structures of rhizonin A and B (derivatives 8 and 9), produced by the symbiont B. endofungorum HKI-0456 (isolate B5).A preceding phylogenetic analysis of the 16S rRNA gene of the type strains B. rhizoxinica HKI-0454 (B1) and B. endofungorum HKI-0456 (B5) showed that both isolates belong to the genus Burkholderia (13). Although the two strains resemble each other in terms of endofungal lifestyle and physiology, DNA-DNA hybridization experiments enforced the division of the two isolates into two distinct species. To establish the taxonomic positions of all eight bacterial symbionts, we isolated genomic DNA from the recovered strains and obtained full-length 16S rRNA gene sequences by PCR using 16S universal primers (15). Sequence comparisons revealed that all isolated endosymbiotic bacteria are closely related to species of the genus Burkholderia.However, the close relationship of the symbionts is particularly intriguing considering the highly diverse collection localities of the host strains (Table ​(Table1).1). Despite the clear grouping of the bacteria associated with Rhizopus, the phylogenetic relationship within the endofungal symbiont complex could not be resolved by 16S rRNA gene data alone (see Fig. S2 in the supplemental material). Several computational methods failed to infer a statistically meaningful phylogeny. To overcome uncertainties in the 16S rRNA gene and biotyping analyses and to further characterize the genotypes of the eight isolates of the B. rhizoxinica complex, we performed a multilocus sequence typing (MLST) analysis. Seven conserved gene loci from all isolates were amplified by PCR, sequenced, and phylogenetically analyzed. To facilitate the comparison with an MLST study of the related bacterium Burkholderia pseudomallei (8), fragments of the following genes were chosen: ace (acetoacetyl-coenzyme A reductase), gltB (glutamate synthase, large subunit), gmhD (ADP-l-glycero-d-manno-heptose-6-epimerase), lepA (GTP-binding protein), lipA (lipoate synthase), and ndh (NADH:ubiquinone oxidoreductase). Shotgun sequencing of the genomes of two symbiont isolates indicated that the narK locus used in the B. pseudomallei study is obviously not present in the endofungal Burkholderia strains (G. Lackner, L. P. Partida-Martinez, and C. Hertweck, unpublished results). Therefore, as a characteristic feature of the ecotype, a locus from the rhizoxin biosynthesis gene cluster, rhiE (14), was sequenced in all isolates. The rhiE locus codes for a part of the polyketide synthase involved in rhizoxin biosynthesis in endofungal bacteria (14). It should be mentioned that a homologous rhizoxin biosynthesis gene cluster has been identified in the phylogenetically distant strain Pseudomonas fluorescens Pf-5 (3, 19).All loci were analyzed independently to test for incongruence between the data sets. The majority of the single-locus trees yielded a topology similar to that shown in Fig. S3 in the supplemental material. Only the ace locus resulted in a different tree: it showed a split (Eurasian, B8) (B5, others) different from the split (Eurasian, B5) (B8, others) found in the majority of single-locus trees (gltB, gmhD, and lipA) (see Fig. S3 in the supplemental material). The remaining loci (lepA, ndh, and rhiE) failed to infer statistically supported clades containing the strain B5 or B8. The phi test for recombination implemented in the program SplitsTree4 (5, 9a) indicated evidence for recombination, if the ace locus was included in the data set (P = 0.044). The removal of the ace locus abolished the signal. We conducted concatenated analyses including gltB, gmhD, lipA, lepA, ndh, and rhiE in the presence and absence of the ace locus. The resulting phylogenetic trees were recovered using distance matrix, maximum-parsimony, and Bayesian methods (see Fig. S3 and S4 in the supplemental material). We found that the tree topology is independent from the presence or absence of the ace locus. An alternative way to illustrate the phylogenetic groupings is in a network (9) (Fig. ​(Fig.3).3). The type strain B. rhizoxinica HKI-0454 (isolate B1) and isolate B6 share identical alleles in all sequenced loci and thus could be considered the same species. This high degree of similarity is supported by matrix-assisted laser desorption ionization protein profiling (see Fig. S5 in the supplemental material). However, it is remarkable that the geographic origins of isolates B1 and B6 (Japan and the United States, respectively) are different. Another member of this “Pacific group,” isolate B2 from Japan, is the closest relative. This observation strongly suggests that the Japanese and U.S. isolates have a common ancestor. Another highly supported clade, the “Eurasian group,” consists of the isolates B3 (Vietnam), B4 (Ukraine), and B7 (Georgia). Again, strain B3 and B7 are highly similar at the nucleotide level despite their geographic distance. Although related to this clade, the B. endofungorum type strain HKI-0456 (isolate B5), isolated from ground nuts in Mozambique, is unique in both genotypic and phenotypic aspects. Isolate B8 from Australia is related even more remotely to all other strains. The phylogenetic data obtained in this study suggest that all Burkholderia symbiont strains found in Rhizopus have a common ancestor.Open in a separate windowFIG. 3.Phylogenetic network of the endofungal symbiont complex (isolates B1 to B8) and B. pseudomallei (BPS), based on MLST. The graph was obtained by the neighbor-net method, implemented with the SplitsTree4 program. Uncertainties in the data are visualized by the network structure in the center.To test whether the phylogenetic data obtained from the MLST analysis contain further information about the evolution of the endofungal symbiosis, we analyzed the extent of detectable recombination in the data set. The presence of recombination could be explained by the exchange of genetic material between bacterial lineages that might have occurred during horizontal transmission of endosymbionts. Less likely, recombination could mean that the mutualistic association was established several times in some of the lineages.The number of incompatible splits in the center of the phylogenetic-network structure (Fig. ​(Fig.3)3) leaves some uncertainty about the correct placement of the strains B8 and B5, and the phi test for recombination indicated evidence for recombination if the ace locus was included. However, recombination does not appear to be a dominant factor in the evolution of the core genome of the known endosymbionts. Notably, the strong congruence between the trees retrieved from the symbiont data is in stark contrast to the recently reported high rate of recombination in free-living Burkholderia spp. (2, 6). Presumably, this is not due to a lack of recombination machinery, as homologous recombination works fine in at least three of the strains under laboratory conditions. Rather, we assume that the mainly vertically transmitted and geographically separated symbionts have evolved primarily separately from each other. Remaining traces of recombination might be hints of coinfection events in the early history of the symbiosis.Another genetic feature of the endosymbionts presented here is their relatively low GC content compared to that for related, but free-living Burkholderia species. The bacterial endosymbiont of aphids, Buchnera aphidicola, is known to have some mutational bias toward low GC content (26). Indeed, all of the endosymbionts have significantly lower GC contents than their sequenced free-living relatives, e.g., B. pseudomallei, Burkholderia thailandensis, and Burkholderia cenocepacia, in all conserved loci (see Fig. S6 in the supplemental material). Although the data obtained in this study are only preliminary evidence and future studies at the whole-genome level could certainly provide more insights into the nucleotide evolution of fungal endosymbionts, it is possible that mechanisms similar to those for Buchnera species are responsible for the reduced GC content in the fungal endosymbionts. The observation that the reduction in GC content is not as striking as that in Buchnera species is then well in accordance with the expectation that the fungus-bacterium endosymbiosis is young compared to the Buchnera-aphid mutualistic relationship.To obtain hints about possible cospeciation or horizontal-transmission events, we compared the phylogenetic relationships between the endobacteria and their fungal hosts. Nucleotide sequences of the 18S ribosomal DNA (rDNA), 28S rDNA, and internal transcribed spacer (ITS) regions were chosen to elucidate phylogenetic relationships between fungal hosts (1). These attempts were hampered since 18S rDNA sequences were highly conserved among the fungal isolates. The only variable site distinguished symbiotic from nonsymbiotic R. microsporus strains (see Fig. S7 in the supplemental material). 28S rDNA sequences are known to accumulate single nucleotide changes at a relatively low rate as well (27). Our 28S rDNA data set, with a total length of 604 nucleotides, contained only eight variable sites. With three of them being parsimony informative, no meaningful phylogram could be inferred from the 28S rDNA sequences. The ITS region is known to evolve more rapidly and is used to provide discrimination within species (27). The curated ITS alignments consisted of 621 sites, 14 being variable and 10 being parsimony informative. Four short insertions or deletions, which are mostly ignored by phylogeny inference software, were found. The phylogram based on both ITS data sets was juxtaposed with the endosymbiont tree determined by MLST (Fig. ​(Fig.4).4). The host strains of the Burkholderia sp. isolates B1, B2, and B6 (ATCC 62417, ATCC 20577, and CBS 261.28, respectively) were identical and reproduced the Pacific group of endobacteria. Strikingly, in contrast to their bacterial partners, the fungal host strains CBS 308.87 (Burkholderia sp. isolate B8) from Australia and CBS 111563 (Burkholderia sp. isolate B3) from Vietnam are members of the Pacific group as well. The strains CBS 700.08 (Burkholderia sp. isolate B7) from Georgia and ATCC 699.68 (Burkholderia sp. isolate B4) from Ukraine appeared to represent the Eurasian branch. This group was known from the bacterial phylogeny, but the fungal clade missed the close relationship to CBS 111563 (Burkholderia sp. isolate B3). Again, the African branch, consisting exclusively of CBS 112285 (Burkholderia sp. isolate B5), shared a common ancestor with the Eurasian group. These results are in accordance with cospeciation of some fungal hosts and their endosymbionts, especially for Burkholderia sp. isolates B1, B2, and B6 and B5, B4, and B7 (Fig. ​(Fig.4).4). However, there might be first evidence for some host switching events in the history of the endofungal bacteria (Burkholderia sp. isolates B3 and B8). Although this hypothesis is based mainly on a few informative sites within the ITS data set, three insertion or deletion events within the alignment support the extended Pacific group (Burkholderia sp. isolates B1, B2, B6, B3, and B8), indicating the horizontal transfer of symbionts and/or genetic material between strains. While it is possible that multiple events led to this unusual symbiosis, a scenario in which all symbiont strains are derived from an ancestral association seems to be more likely. The endosymbiont-dependent sporulation of the host strain indicates that the fungus-bacterium interaction is highly specialized. Furthermore, vertical transmission of the symbionts through spores is an efficient strategy for rapid distribution (4, 20). Nonetheless, our data suggest that the horizontal transmission of symbionts might also have played a role during the evolution of the endofungal bacteria.Open in a separate windowFIG. 4.Juxtaposition of phylogenetic trees derived from the MLST data of the endofungal symbiont complex (isolates B1 to B8) and the ITS sequence data of the fungal host (strain designations of fungal isolates of the genus Rhizopus are given). Dashed lines are representative of a symbiotic relationship. The numbers on top of the branches indicate the clade probability values. Shading designates similar clade affiliations for the bacterial symbiont and the fungal host.In conclusion, we have investigated eight bacterial endosymbiont strains isolated from toxinogenic R. microsporus strains in pure culture. All isolates are representatives of the same unique “endofungal” ecotype, albeit the hosts'' origins cover all five continents and occur in highly diverse niches. The bacterial endosymbionts share characteristic phenotypic traits, like secondary metabolite production and protein profile, as demonstrated by HPLC-MS and matrix-assisted laser desorption ionization-time of flight biotyping, respectively. Phylogenetic analyses (16S rRNA genes) provide strong evidence that all symbiont strains originate from a common ancestor and form a new complex within the genus Burkholderia. This observation is strongly supported by MLST, according to which all eight symbiont isolates can be grouped into continental branches. Results revealing both similar and deviating geographical groupings of fungal isolates in comparison to bacterial endosymbionts allow hypothesizing about the possible cospeciation of fungal and bacterial symbionts and some extent of horizontal-transmission events. All bacterial strains investigated seem to have evolved mainly separately from each other, not showing extensive recombination. In addition, we present preliminary evidence that there might be a mutational bias toward high AT contents, as is known for other endosymbiotic bacteria.