Horizontally acquired genes in early-diverging pathogenic fungi enable the use of host nucleosides and nucleotides

William G. Alexander, Jennifer H. Wisecaver, Antonis Rokas, and Chris Todd Hittinger
PNAS April 12, 2016 113 (15) 4116-4121; first published March 29, 2016 https://doi.org/10.1073/pnas.1517242113

  1. Edited by Wen-Hsiung Li, Academia Sinica, Taipei, Taiwan, and approved March 2, 2016 (received for review August 28, 2015)

Significance

Early-diverging fungi from the phylum Microsporidia are opportunistic pathogens of humans and other animals. A genome-wide search was conducted in four species of Microsporidia and one species of Cryptomycota for genes acquired by horizontal transfer from other organisms. Up to 2% of their genes had strong evidence for horizontal transfer. We showed that many transferred nucleic acid metabolism genes lie on the boundary of the host and pathogen metabolic networks. Finally, we functionally characterized the gene encoding thymidine kinase, whose multiple parallel transfers from three different sources strongly suggest an integral adaptive role in the lifestyles of these intracellular parasites. Microsporidian thymidine kinases activate a prodrug, suggesting a treatment route for microsporidian infections.

Abstract

Horizontal gene transfer (HGT) among bacteria, archaea, and viruses is widespread, but the extent of transfers from these lineages into eukaryotic organisms is contentious. Here we systematically identify hundreds of genes that were likely acquired horizontally from a variety of sources by the early-diverging fungal phyla Microsporidia and Cryptomycota. Interestingly, the Microsporidia have acquired via HGT several genes involved in nucleic acid synthesis and salvage, such as those encoding thymidine kinase (TK), cytidylate kinase, and purine nucleotide phosphorylase. We show that these HGT-derived nucleic acid synthesis genes tend to function at the interface between the metabolic networks of the host and pathogen. Thus, these genes likely play vital roles in diversifying the useable nucleic acid components available to the intracellular parasite, often through the direct capture of resources from the host. Using an in vivo viability assay, we also demonstrate that one of these genes, TK, encodes an enzyme that is capable of activating known prodrugs to their active form, which suggests a possible treatment route for microsporidiosis. We further argue that interfacial genes with well-understood activities, especially those horizontally transferred from bacteria or viruses, could provide medical treatments for microsporidian infections.

Horizontal gene transfer (HGT), or the nonvertical transmission of genetic information between distantly related organisms, is common in bacteria, archaea, and viruses (1⇓–3). The importance and scale of HGT in eukaryotes, however, is a matter of debate (4⇓⇓–7). In particular, HGT into fungi was thought to be rare, but several examples from bacteria into fungi or between fungi have recently been described (5, 8, 9). For example, one phylum of early-diverging fungi, the Microsporidia, has had a handful of HGT events documented (5, 10⇓–12), but these fungi are thought to have acquired relatively few genes through HGT (13).

Microsporidia are obligate intracellular parasites of animals. Their genomes are highly compact, and they have eliminated many core metabolic processes in favor of relying on their host for the synthesis of essential molecules (14). They are opportunistic pathogens that primarily infect immunocompromised individuals, such as AIDS patients and organ-transplant recipients (15). They also infect a number of economically important animals as potent zoonotic pathogens (16). Members of a related early-diverging fungal phylum, Cryptomycota, are also obligate intracellular parasites that infect algae (17), amoeboids (18), and other fungi (19), highlighting the similarity between their lifestyles and their evolutionary affinity.

Prior analysis of microsporidian genomes demonstrated that, like other microbial eukaryotes, this phylum of fungi has obtained multiple genes through HGT (20), which may have provided a portion of the raw material required for adaptation, as seen in other organisms (21). Known HGT events into Microsporidia include the following: (i) the ADP/ATP translocase gene family, originating from an HGT event that transferred the founding gene from a member of the bacterial phylum Chlamydia (10), which are known to steal energy-bearing molecules from their host; (ii) a six-gene folate synthesis pathway transferred into Encephalitozoon hellem from multiple donors, a transfer hypothesized to reduce host metabolic stress (11); and (iii) the acquisition of a glutamate-ammonia ligase from an unknown prokaryotic source by Spraguea lophii, which is thought to provide spores a mechanism for defense against the ammonia generated by the decomposing flesh in which they are embedded (12). These case studies, the first of which is shared with the Cryptomycota (19), suggest a role for HGT in the evolution of their unusual pathogenic metabolisms and indicate the need for a thorough, systematic analysis of HGT events in these early-diverging fungi.

Results

To quantify the number of HGT-derived genes in these intracellular parasites, we deployed a high-throughput analytical pipeline to analyze the mode of inheritance of 14,914 genes (Table 1) from sequenced genomes of Encephalitozoon cuniculi (22, 23), E. hellem (11), Nosema ceranae (24, 25), Nematocida parisii (26), and Rozella allomycis (19), the latter being the sole Cryptomycota genome available (Table S1). Specifically, for each gene, we first generated Alien Index (AI) scores (27) that compared the similarity of the gene between specified ingroup and outgroup taxa (e.g., fungi and bacteria, respectively; Fig. S1). We chose a more relaxed AI score cutoff than the original study because the relatively small number of genes in these organisms made manual curation of gene trees feasible. Using this AI cutoff, we were able to recover 8 of the 10 previously noted HGT events that fell within the scope of the species surveyed (5, 10, 11, 13, 22, 26, 28, 29); even using our relaxed AI threshold, two previously described HGTs were not identified, in one case narrowly (Table S2).

View this table:
Table 1.

HGT events among early-diverging fungi

View this table:
Table S1.

Genomes analyzed

Fig. S1.
Fig. S1.

Workflow for genome-wide AI score calculation. See Methods for specific details on genome-wide AI calculation.

View this table:
Table S2.

Calibration of HGT calling methodology using prior works

Detection of most, but not all, known HGT events suggests that our methodology is sensitive but conservative. Between 44 and 64 microsporidian genes per genome and 327 R. allomycis genes were designated as high-positive AI genes (genes with AI ≥ 20 or with 10 ≤ AI < 20 and no significant BLAST hit to other fungi; Dataset S1). The high-positive AI genes were cross-referenced with Eukaryotic Orthologous Group (KOG), Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Ontology (GO), and InterPro annotations provided by the Joint Genome Institute MycoCosm Portal (30). Between 39 and 56 microsporidian genes per genome and 268 R. allomycis genes had at least one functional annotation and were designated as genes of interest. These genes of interest were analyzed phylogenetically, and 10–22 microsporidian genes per genome and 57 R. allomycis genes had strong phylogenetic support for HGT (Table 1, Fig. 1, and Dataset S2). In addition, the majority of these likely transferred genes were nested within well-established higher-order taxa, which further supports their HGT origin. Depending on the stringency of criteria applied, the range of HGT events varied between 0.34% and 2.11% of the total genes present in a genome.

Fig. 1.
Fig. 1.

Plots of the AI for five early-diverging fungal eukaryotes. (A) AI scores were calculated for every gene in the five genomes analyzed, and they were ordered by decreasing AI score. (B) Genes with an AI score ≥ 10 were plotted with colors corresponding to whether they passed or failed the AI score or phylogenetic filters (see main text).

Many of the putatively HGT-derived genes identified in this screen are annotated as aminoacyl-tRNA synthetases, DNA repair enzymes, and nucleic acid anabolic enzymes (Dataset S2). The latter group of genes is particularly interesting because it includes previously characterized HGTs. For example, purine nucleotide phosphorylases (PNPs) and phosphoribosyltransferases (PRTs) are hypothesized to play roles in guanine metabolism in E. hellem (11), whereas ADP:ATP translocases transferred from Chlamydia are found in all Microsporidia (10). In addition, cytidylate kinase (CK) was previously noted as a possible HGT (5, 26), but its role in Microsporidian biology and its connection to other HGTs is uncharacterized. These, together with the newly identified nucleic acid salvage gene thymidine kinase (TK), may enable the harvest of host nucleosides and nucleotides.

Closer examination of the nucleic acid subpathway in the Microsporidia E. hellem revealed it to be completely disconnected from the core metabolic network, an observation conserved in the other microsporidians examined (Fig. 2). Thus, many, if not all, Microsporidia are unable to synthesize nucleic acid components de novo and must rely on scavenging free nucleic acid components from the host. Remarkably, each of the HGT-derived genes in the E. hellem metabolic subpathway lies on or near the perimeter. The only HGT-derived gene (PNP) that is removed by one step coexists with an HGT-derived gene (PRT) that directly connects it to the perimeter of the subpathway (11, 29). Thus, each HGT provides a new interface that connects the pathogen's nucleic acid subpathway to the host metabolic network. To test the hypothesis that HGT-derived metabolic nodes tend to lie at the host–pathogen interface, we permuted the E. hellem nucleic acid metabolic subpathway, independently evaluating enzymatic steps and genes for the number of steps remaining to the perimeter. We found that the observed subnetwork was among the most extreme permutations (P = 0.0027 and P = 0.019, respectively; Dataset S3), leading us to focus on the genes at the perimeter of this network.

Fig. 2.
Fig. 2.

Early-diverging fungal HGT events and nucleic acid metabolism. A schematic of the E. hellem nucleic acid anabolic subpathway indicates that HGT-derived genes are found near the host–parasite interface, which provides the parasite with additional sources of nucleosides and nucleotides. HGT-derived enzymatic functions are denoted by red arrows (10, 11, 29), whereas the dCMP deaminase present in some other Microsporidia (but absent in E. hellem) is denoted by a dashed black arrow. HGTs are also contained within the E. hellem folate biosynthesis pathway (orange), which donates a methyl group to convert dUMP to dTMP. CK, cytidylate kinase; PNP, purine nucleotide phosphorylase; PRT, phosphoribosyltransferase; TK, thymidine kinase.

The transfer of TK into microsporidians would enable the direct phosphorylation of thymidine from the host into dTMP for further synthesis by the pathogen, relieving shortages that might be encountered at various stages of the host’s cell cycle. In particular, TK reduces the need for folate by microsporidians, the ultimate source of which is the diet of the animal host. Folate donates a methyl group to convert dUMP to dTMP, and this particular stress has been invoked previously to explain the horizontal transfer of several genes encoding folate synthesis enzymes into E. hellem (11). Interestingly, not all microsporidians possess the bacterial TK, such as Vavraia culicis and Trachipleistophora hominis. Instead, these organisms have retained an ancestral gene encoding a deoxycytidine monophosphate deaminase, which produces dUMP and, ultimately, dTMP from dCMP (Fig. 2). Thus, microsporidians have evolved three distinct, but sometimes overlapping, strategies for the production of dTMP, suggesting that production of this nucleotide is central to their intracellular parasitic lifestyle.

Of the three strategies for making dTMP, the HGT of TK stands out because the gene is thought to have been lost in the fungal lineage shortly after it diverged from animals (31). The phylogeny of TK sequences supports two independent horizontal transfers of bacterial TK genes into several Microsporidia taxa and into R. allomycis (Fig. 3). Likelihood ratio tests rejected alternate tree topologies where the Microsporidia TKs were sister to either Animalia or Eukarya TKs; likewise, alternate topologies in which the R. allomycis TK was sister to either Animalia or Eukarya TKs were rejected. Finally, an alternate topology in which R. allomycis and the Microsporidia TKs were grouped together was also rejected, further supporting the independent transfer of TK into these phyla (Table 2). In addition to these two HGT events, we found an independent HGT of a putative, viral-like TK into N. parisii (Fig. 4 and Fig. S2), documenting a third parallel transfer of TK into these fungal intracellular parasites.

Fig. 3.
Fig. 3.

TK has been transferred from bacteria into Microsporidia and Cryptomycota by independent HGT events. Top BLASTP hits for E. cuniculi, No. ceranae, and R. allomycis TK were aligned along with TK sequences from representative eukaryotes and other bacterial phyla, and the resulting alignment was used to infer a maximum-likelihood phylogeny. Red stars indicate the location on the tree where likelihood-ratio tests were performed to evaluate the relationship of the sequences in question with Animalia or other Eukarya representatives. Bootstrap values for key nodes are noted. The N. parisii TK results from a third independent HGT, but is absent from the phylogenetic tree because of its distinctive viral domain structure and amino acid sequence (Fig. S2). Note that the No. ceranae TK has also been excluded due its long branch length (Dataset S2), but Nosema bombycis is shown.

View this table:
Table 2.

Likelihood ratio tests support two independent transfers into Microsporidia and Cryptomycota

Fig. 4.
Fig. 4.

Schematic of nucleic acid metabolism HGT events in the Microsporidia and Cryptomycota. A cladogram of representative members of early-diverging fungi highlighting HGT-derived genes that function in the nucleic acid subpathway. Green boxes denote three independent horizontal transfers of TK, while orange boxes denote other horizontal transfers of nucleic acid metabolism genes (10, 11, 27).

Fig. S2.
Fig. S2.

N. parisii has a putative viral TK. (A) Domain analysis of the N. parisii gene suggests the presence of a viral-like TK gene similar to those found in herpesviruses from various animals. (B) Phylogenetic analysis suggests a loose relationship between this hypothetical N. parisii gene and herpesvirus TK.

The absolute dependence of microsporidians on their hosts for nucleic acid components suggests that this subpathway could be a viable target for antimicrosporidian drugs. TK is an especially promising target for disruption because of the apparent importance placed on dTMP synthesis during evolution and the wide range of TK-activated prodrugs currently used as antiviral therapies. First, we confirmed that TK is expressed by analyzing a previously published E. cuniculi transcriptome dataset (32); TK was in the top quartile of expressed genes at all three postinfection time points (Table S3 and Dataset S4). We then evaluated the predicted TK proteins of E. cuniculi, N. parisii, and R. allomycis in an assay in which the nucleoside kinase activity of TK would activate the prodrug 5-fluoro-2-deoxyuridine (FUdR), converting it into the toxic compound fluorodeoxyuridine monophosphate and killing the tester strain of Saccharomyces cerevisiae (31). The growth of strains expressing Homo sapiens, E. cuniculi, and N. parisii TK were completely inhibited when FUdR was present, whereas strains expressing either predicted version of the R. allomycis TK behaved identically to the blank vector control strain, perhaps due to genome annotation or heterologous expression issues. Indeed, we also noted that the N. parisii TK strain did not grow as robustly as the other four TK-expressing strains or the blank vector control strain in any media tested, suggesting that N. parisii TK is detrimental to S. cerevisiae growth (Table 3). Thus, we conclude that at least two of the three TK proteins that underwent HGT—one from viruses and one from bacteria—are active and capable of activating known prodrugs.

View this table:
Table 3.

Endpoint analysis of optical density of TK-expressing strains relative to control cultures

View this table:
Table S3.

Expression analysis of E. cuniculi horizontally transferred and control genes

Discussion

We have shown that members of the early-diverging fungal phyla Microsporidia and Cryptomycota have acquired dozens of genes by HGT, the vast majority of which are, to our knowledge, reported here for the first time. These findings are inconsistent with the hypothesis that intracellular parasites possess few HGT-derived genes (13). Despite their prevalence, the mechanism by which these genes are transferred into early-diverging fungi is still unclear, and previously proposed eukaryotic HGT mechanisms, such as through viral intermediates (7, 33), seem unlikely for intracellular parasites. One possible mechanism is that phagocytosis of spores by host cells into intracellular vesicles (34) may bring nucleic acids from lysed environmental bacteria into contact with receptive spores. Because microsporidians reside in the host cytoplasm, another possibility is that host mRNA molecules are reverse-transcribed and incorporated into microsporidian genomes. Whatever the molecular mechanism of the HGTs, their impact on the microsporidian genome and metabolism has been widespread.

Nucleic acid anabolic genes, in particular, have been transferred multiple times into the Cryptomycota and Microsporidia. These HGTs provide the intracellular parasites with novel sources of nucleosides and nucleotides, while simultaneously reducing metabolic stress on their hosts. This evolutionary strategy has been described in the intracellular Apicomplexan parasite Cryptosporidium parvum (35), which, like many microsporidians, has acquired a bacterial TK to diversify its nucleotide sources. Herpes viruses also rely on TK to increase dTMP levels, and TK strongly affects its pathogenicity (36).

Previous work has shown that N. parisii replication can be slowed by FUdR treatment, although whether this effect was a function of host or parasite metabolism was unclear (37). Here we show that some Microsporidia TK genes, including N. parisii TK, indeed encode active nucleoside kinases, providing a likely drug target with a known mechanism of action (38, 39). Thus, antiviral prodrugs that are activated by TK, such as FUdR or acyclovir, may offer a path to a viable treatment option for microsporidiosis. These drugs are already used to treat diseases, such as herpes simplex virus (40), HIV (41), and certain types of cancers (42). In the latter cases, cancer cells are preferentially targeted because a viral exogenous TK has been inserted into the cell that has lower specificity than the endogenous TK (38), or native human TK expression has been increased due to uncontrolled cellular division (43). Although we did not detect any activity for R. allomycis TK, a negative result in a heterologous gene expression assay is not definitive. R. allomycis TK may require a chaperone that S. cerevisiae does not have; it may require a posttranslational modification that S. cerevisiae is unable to perform; it may function in a narrow range of cellular conditions; or the intron–exon boundaries may have been misannotated.

The treatment of fungal pathogens of animals has traditionally proven challenging due to the shared physiological traits of these sister kingdoms. Because of its inability to synthesize nucleic acid components, nucleotide analogs that are activated by the existing microsporidian nucleic acid subnetwork could provide new treatments for microsporidiosis in humans and livestock. For example, bacterial CK prefers CMP or dCMP as a substrate, whereas the UMP/CMP kinase found in animals and fungi prefers CMP or UMP as a substrate (44, 45). Because Microsporidia possess the bacterial CK, they may be much more susceptible to dCMP analogs than their host cells, a treatment strategy that has been used against HIV (46). Many other pathogenic fungi that infect animals, such as Aspergillus fumigatus, have also acquired an appreciable number of genes through HGT (47), including entire secondary metabolism gene clusters (48). Using HGT events to identify novel candidate drug targets may thus prove a promising general strategy across diverse clades of pathogenic fungi and other recalcitrant pathogens whose physiology is similar to humans.

Methods

AI Analysis of Early-Diverging Fungal Genomes.

All predicted proteins from five early-diverging fungal genomes (Table S1) were queried by using BLAST (E-value = 0.001, max target seqs = 1,000) against a custom database consisting of the National Center for Biotechnology Information’s (NCBI’s) nonredundant protein database (last updated November 24, 2014), as well as additional protein sequences from 411 fungal and plant genome assemblies (Dataset S5). These additional genomes served to increase the representation of close relatives to microsporidians in the database, which is critical when searching for HGTs based on relative BLAST scores. They also allowed the taxonomic extent of HGTs into the Microsporidia to be assessed without risking false positives that could be caused by genomes of variable quality (e.g., bacterial contamination). An AI approach was used to screen these genomes for genes with significantly better BLAST hits to distantly related organisms (e.g., metazoans or bacteria) than to closely related ones (e.g., other fungi) and thereby identify HGT candidates. Two taxonomic lineages were first specified: the recipient lineage into which possible HGT events may have occurred (i.e., Microsporidia or Cryptomycota, depending on the query) and a larger group of related taxa (i.e., Fungi). When parsing the BLAST output, all hits to the recipient lineage were skipped. The AI score is given by the formula:AI=(ln(bbhG+1×10−200)−ln(bbhO+1×10−200)),

where bbhG is the E-value of the best BLAST hit to a species within the group lineage (i.e., the best nonmicrosporidian or, in the case of R. allomycis, noncryptomycotan fungal match) and bbhO is the E-value of the best blast hit to a species outside of the group lineage (i.e., the best nonfungal match). In cases where there were no significant BLAST hits, the corresponding bbhG or bbhO was set to 1. AI can range from 461 to −461, and AI is greater than zero if the gene has a better BLAST hit to a species outside of the group lineage. Contamination in the genome assembly could also result in a positive AI score; therefore, any assembly contigs containing only genes with positive AI scores were flagged as potential contamination and eliminated from the analysis.

Systematic Identification of Likely HGT Events.

Genes with AI ≥ 20 or genes with 20 > AI ≥ 10 and with no significant BLAST hit to other Fungi, were classified as high-positive AI genes. These genes were cross-referenced with KOG, KEGG, GO, and InterPro annotations provided by the Joint Genome Institute (30), and all genes with at least one functional annotation were classified as genes of interest and analyzed by using a phylogenetic pipeline (Fig. S1). For each HGT candidate, a Perl script based publicly available software (49) extracted up to 400 homologs from the custom database (referenced above) based on BLAST similarity, allowing up to five orders of magnitude difference between query and subject lengths, and extracting up to five sequences per species. To reduce the number of sequences per gene tree, highly similar sequences were collapsed with CD-HIT using default parameters (50). Sequences were aligned with MAFFT by using the E-INS-i strategy (51). The resulting alignment was trimmed with TRIMAL by using the automated1 strategy (52). All genes with trimmed alignments <100 amino acids were discarded. Phylogenetic trees were constructed by using RAxML (53) with the PROTGAMMAAUTOF amino acid model of substitution and 100 bootstrap replicates. Trees were midpoint-rooted, and branches with <50% bootstrap support were collapsed by using TreeCollapseCL4 (emmahodcroft.com/TreeCollapseCL.htm). The resulting phylogenies were manually inspected to assess each gene’s mode of transmission. First, the presence or absence of fungal BLASTP hits was determined. Second, the nearest neighbor node to the query was identified as prokaryotic, eukaryotic, opisthokont, and/or fungal (multiple nearest neighbors are possible because poorly supported nodes were collapsed). Next, if the nearest node was opisthokont, the likelihood of parallel loss to explain the topology was assessed. Finally, we noted whether the query was nested within a well-established alien clade. A gene’s mode of transmission was labeled as a likely HGT event if and only if its phylogeny was inconsistent with the organismal phylogeny (19) and any parallel-loss scenario (e.g., a protein sequence from E. cuniculi being nested within a clade of Proteobacteria instead of sister to the Fungi). If a phylogeny could be explained by simple parallel loss or any mechanism other than horizontal transfer, its mode of transmission was labeled as unlikely HGT. If the phylogeny provided no clear support for either vertical or horizontal transmission, its mode was labeled as ambiguous. Finally, some phylogenies were labeled as “no call” because they lacked an interpretable topology (e.g., “star phylogeny”). All phylogenetic trees and calls are available in Dataset S2.

Permutation Analyses of the E. hellem Nucleic Acid Subpathway.

The enzymatic steps at the periphery of the complete E. hellem nucleic acid subpathway interface directly with the metabolism of the host (Fig. 2). These steps had a distance to the host of zero, whereas other steps were assigned a distance corresponding to the smallest number of steps from their nearest peripheral metabolites. The observed average distance to the host for enzymatic steps was 0.167 for steps experiencing HGTs, whereas the average for the complete subpathway was 1.394. Because CK (and some non-HGT enzymes) can perform multiple steps, we also considered genes in a separate analysis by averaging the values of each step they encode, yielding an observed average distance of 0.2 for genes experiencing HGTs vs. an average of 1.163 for the complete subpathway. To assess significance, 1 million random permutations were simulated, and the proportion of permuted networks that were at least as extreme as the observed network was determined. Complete documentation is provided in Dataset S3.

Phylogenetic Analysis of TK.

The E. cuniculi and R. allomycis TK sequences were used as BLASTP queries via the NCBI BLAST website (www.ncbi.nlm.nih.gov/BLAST), and the top hits were downloaded along with the TK sequences of representatives from eukaryotic and prokaryotic clades. These protein sequences were aligned in MUSCLE (54) by using default settings, and maximum-likelihood analysis was performed by using the PROTGAMMALG model of evolution in multithreaded RAxML (53). To evaluate alternative possibilities for the source of these sequences, five different alternate trees were drawn, also by using the PROTGAMMALG model in RAxML, and the likelihood ratio was calculated for each tree compared with the best tree (55). P values were calculated from this ratio in the Python Scipy library χ2 probability module (56).

Plasmid and Strain Construction.

The protein sequence for soluble TK of H. sapiens was downloaded from NCBI (BAG70082.1), and the TK protein sequences for E. cuniculi (ECU01_0740im.01) (22), N. parisii (Nempa11372) (26), and R. allomycis (Rozal12947) (19) were downloaded from JGI; note that the R. allomycis constructs were designed twice, once without and once with predicted intron sequences present, resulting in R. allomycis TK versions A and B, respectively. These protein sequences were codon-optimized for expression in S. cerevisiae by using the Integrated DNA Technologies (IDT) webtool (www.idtdna.com/CodonOpt). Although critical to ensuring proper heterologous protein expression (57), this approach necessarily eliminates any RNA regulation that may occur in E. cuniculi, N. parisii, and R. allomycis. gBlocks for the H. sapiens and E. cuniculi sequences, as well as plasmids containing synthesized sequence for the N. parisii and both versions of R. allomycis, were ordered from IDT. Primers (Table S4) were designed to amplify each gBlock or synthesized gene, while adding ∼40 bp of sequence on both ends to enable homology repair (58) of the gBlock sequence into the plasmid pRS426-HygMX under the control of the TDH3 promoter and the CYC1 terminator. This cloning was performed in the wild diploid strain M22 of S. cerevisiae (59) by transforming it using linearized pRS426-HygMX and PCR product via the lithium acetate/PEG method (60). Plasmid repair and yeast transformation was selected for on yeast extract/peptone/dextrose (YPD) +Hyg (10 g/L yeast extract, 20 g/L peptone, 20 g/L glucose, 18 g/L agar, and 200 mg/L hygromycin B) after 3 h of recovery in liquid YPD. Resulting yeast colonies had their total DNA extracted, which was used to transform Escherichia coli by electroporation, and the recovered plasmids had their insertions Sanger-sequenced. Once the gene-insertion sequences were confirmed, plasmids were again transformed into M22 by selection on YPD +Hyg plates as described above.

View this table:
Table S4.

Primers used in this study

Functional Assay of TK.

Single colonies of the six plasmid-bearing strains of S. cerevisiae (Table S5) were grown overnight in SC-MSG +Hyg (1.72 g/L yeast nitrogen base, 2 g/L Complete Drop-Out Mix, 1 g/L monosodium glutamate, 20 g/L glucose, and 200 mg/L hygromycin B) liquid medium. A total of 30,000 cells of each strain were used to inoculate 3 mL of SC-MSG +Hyg +100 mg/L FUdR in quadruplicate; control medium tubes containing SC-MSG +Hyg without FUdR were also inoculated with 30,000 cells for each strain in duplicate. All cultures, except for the N. parisii TK-expressing strain, were grown for 60 h, and their optical density was measured with a 600-nm photometer (Implen GmbH). The N. parisii cultures were measured after 168 h.

View this table:
Table S5.

Strains used in this study

Acknowledgments

We thank Issac Knoflicek, Michael G. Mangiore, and Dana Opulente for computational support and advice and Amanda B. Hulfachor for artwork. This work was conducted in part using the resources of the Advanced Computing Center for Research and Education at Vanderbilt University. This material is based upon work supported by National Science Foundation Grants IOS-1401682 (to J.H.W.), DEB-1442113 (to A.R.), and DEB-1253634 and DEB-1442148 (to C.T.H.); in part by the DOE Great Lakes Bioenergy Research Center through DOE Office of Science Grant BER DE-FC02-07ER64494; and USDA (United States Department of Agriculture) National Institute of Food and Agriculture Hatch Project 1003258. C.T.H. is an Alfred Toepfer Faculty Fellow and Pew Scholar in the Biomedical Sciences, supported by the Alexander von Humboldt Foundation and Pew Charitable Trusts, respectively.

Footnotes

  • ↵1To whom correspondence should be addressed. Email: cthittinger{at}wisc.edu.

  • Author contributions: W.G.A. and C.T.H. conceived the project; W.G.A., J.H.W., A.R., and C.T.H. designed research; W.G.A., J.H.W., and C.T.H. performed research; W.G.A., J.H.W., and C.T.H. contributed new reagents/analytic tools; W.G.A., J.H.W., and C.T.H. analyzed data; and W.G.A., J.H.W., A.R., and C.T.H. wrote the paper.

  • The authors declare no conflict of interest.

  • This article is a PNAS Direct Submission.

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1517242113/-/DCSupplemental.

References

  1. ↵
    1. Pál C,
    2. Papp B,
    3. Lercher MJ
    (2005) Adaptive evolution of bacterial metabolic networks by horizontal gene transfer. Nat Genet 37(12):1372–1375
    .
    OpenUrlCrossRefPubMed
  2. ↵
    1. Treangen TJ,
    2. Rocha EPC
    (2011) Horizontal transfer, not duplication, drives the expansion of protein families in prokaryotes. PLoS Genet 7(1):e1001284
    .
    OpenUrlCrossRefPubMed
  3. ↵
    1. Moreira D,
    2. López-García P
    (2009) Ten reasons to exclude viruses from the tree of life. Nat Rev Microbiol 7(4):306–311
    .
    OpenUrlCrossRefPubMed
  4. ↵
    1. Keeling PJ,
    2. Palmer JD
    (2008) Horizontal gene transfer in eukaryotic evolution. Nat Rev Genet 9(8):605–618
    .
    OpenUrlCrossRefPubMed
  5. ↵
    1. Marcet-Houben M,
    2. Gabaldón T
    (2010) Acquisition of prokaryotic genes by fungal genomes. Trends Genet 26(1):5–8
    .
    OpenUrlCrossRefPubMed
  6. ↵
    1. Huang J
    (2013) Horizontal gene transfer in eukaryotes: The weak-link model. BioEssays 35(10):868–875
    .
    OpenUrlPubMed
  7. ↵
    1. Schönknecht G,
    2. Weber APM,
    3. Lercher MJ
    (2014) Horizontal gene acquisitions by eukaryotes as drivers of adaptive evolution. BioEssays 36(1):9–20
    .
    OpenUrlCrossRefPubMed
  8. ↵
    1. Wisecaver JH,
    2. Slot JC,
    3. Rokas A
    (2014) The evolution of fungal metabolic pathways. PLoS Genet 10(12):e1004816
    .
    OpenUrlCrossRefPubMed
  9. ↵
    1. Wisecaver JH,
    2. Rokas A
    (2015) Fungal metabolic gene clusters-caravans traveling across genomes and environments. Front Microbiol 6:161
    .
    OpenUrlPubMed
  10. ↵
    1. Tsaousis AD, et al.
    (2008) A novel route for ATP acquisition by the remnant mitochondria of Encephalitozoon cuniculi. Nature 453(7194):553–556
    .
    OpenUrlCrossRefPubMed
  11. ↵
    1. Pombert J-F, et al.
    (2012) Gain and loss of multiple functionally related, horizontally transferred genes in the reduced genomes of two microsporidian parasites. Proc Natl Acad Sci USA 109(31):12638–12643
    .
    OpenUrlAbstract/FREE Full Text
  12. ↵
    1. Campbell SE, et al.
    (2013) The genome of Spraguea lophii and the basis of host-microsporidian interactions. PLoS Genet 9(8):e1003676
    .
    OpenUrlCrossRefPubMed
  13. ↵
    1. Nakjang S, et al.
    (2013) Reduction and expansion in microsporidian genome evolution: New insights from comparative genomics. Genome Biol Evol 5(12):2285–2303
    .
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Keeling PJ,
    2. Fast NM
    (2002) Microsporidia: Biology and evolution of highly reduced intracellular parasites. Annu Rev Microbiol 56:93–116
    .
    OpenUrlCrossRefPubMed
  15. ↵
    1. Didier ES
    (2005) Microsporidiosis: An emerging and opportunistic infection in humans and animals. Acta Trop 94(1):61–76
    .
    OpenUrlCrossRefPubMed
  16. ↵
    1. Santín M,
    2. Fayer R
    (2011) Microsporidiosis: Enterocytozoon bieneusi in domesticated and wild animals. Res Vet Sci 90(3):363–371
    .
    OpenUrlCrossRefPubMed
  17. ↵
    1. Ishida S,
    2. Nozaki D,
    3. Grossart H-P,
    4. Kagami M
    (2015) Novel basal, fungal lineages from freshwater phytoplankton and lake samples. Environ Microbiol Rep 7(3):435–441
    .
    OpenUrlCrossRefPubMed
  18. ↵
    1. Corsaro D, et al.
    (2014) Rediscovery of Nucleophaga amoebae, a novel member of the Rozellomycota. Parasitol Res 113(12):4491–4498
    .
    OpenUrlCrossRefPubMed
  19. ↵
    1. James TY, et al.
    (2013) Shared signatures of parasitism and phylogenomics unite Cryptomycota and microsporidia. Curr Biol 23(16):1548–1553
    .
    OpenUrlCrossRefPubMed
  20. ↵
    1. Alsmark C, et al.
    (2013) Patterns of prokaryotic lateral gene transfers affecting parasitic microbial eukaryotes. Genome Biol 14(2):R19
    .
    OpenUrlCrossRefPubMed
  21. ↵
    1. Hall C,
    2. Brachat S,
    3. Dietrich FS
    (2005) Contribution of horizontal gene transfer to the evolution of Saccharomyces cerevisiae. Eukaryot Cell 4(6):1102–1115
    .
    OpenUrlAbstract/FREE Full Text
  22. ↵
    1. Katinka MD, et al.
    (2001) Genome sequence and gene compaction of the eukaryote parasite Encephalitozoon cuniculi. Nature 414(6862):450–453
    .
    OpenUrlCrossRefPubMed
  23. ↵
    1. Peyretaillade E, et al.
    (2009) Identification of transcriptional signals in Encephalitozoon cuniculi widespread among Microsporidia phylum: Support for accurate structural genome annotation. BMC Genomics 10:607–619
    .
    OpenUrlCrossRefPubMed
  24. ↵
    1. Chen YP, et al.
    (2009) Morphological, molecular, and phylogenetic characterization of Nosema ceranae, a microsporidian parasite isolated from the European honey bee, Apis mellifera. J Eukaryot Microbiol 56(2):142–147
    .
    OpenUrlCrossRefPubMed
  25. ↵
    1. Cornman RS, et al.
    (2009) Genomic analyses of the microsporidian Nosema ceranae, an emergent pathogen of honey bees. PLoS Pathog 5(6):e1000466
    .
    OpenUrlCrossRefPubMed
  26. ↵
    1. Cuomo CA, et al.
    (2012) Microsporidian genome analysis reveals evolutionary strategies for obligate intracellular growth. Genome Res 22(12):2478–2488
    .
    OpenUrlAbstract/FREE Full Text
  27. ↵
    1. Gladyshev EA,
    2. Meselson M,
    3. Arkhipova IR
    (2008) Massive horizontal gene transfer in bdelloid rotifers. Science 320(5880):1210–1213
    .
    OpenUrlAbstract/FREE Full Text
  28. ↵
    1. Heinz E, et al.
    (2012) The genome of the obligate intracellular parasite Trachipleistophora hominis: New insights into microsporidian genome dynamics and reductive evolution. PLoS Pathog 8(10):e1002979
    .
    OpenUrlCrossRefPubMed
  29. ↵
    1. Selman M, et al.
    (2011) Acquisition of an animal gene by microsporidian intracellular parasites. Curr Biol 21(15):R576–R577
    .
    OpenUrlCrossRefPubMed
  30. ↵
    1. Grigoriev IV, et al.
    (2014) MycoCosm portal: Gearing up for 1000 fungal genomes. Nucleic Acids Res 42(Database issue):D699–D704
    .
    OpenUrlAbstract/FREE Full Text
  31. ↵
    1. Alexander WG,
    2. Doering DT,
    3. Hittinger CT
    (2014) High-efficiency genome editing and allele replacement in prototrophic and wild strains of Saccharomyces. Genetics 198(3):859–866
    .
    OpenUrlAbstract/FREE Full Text
  32. ↵
    1. Grisdale CJ,
    2. Bowers LC,
    3. Didier ES,
    4. Fast NM
    (2013) Transcriptome analysis of the parasite Encephalitozoon cuniculi: An in-depth examination of pre-mRNA splicing in a reduced eukaryote. BMC Genomics 14:207
    .
    OpenUrlCrossRefPubMed
  33. ↵
    1. Liu H, et al.
    (2010) Widespread horizontal gene transfer from double-stranded RNA viruses to eukaryotic nuclear genomes. J Virol 84(22):11876–11887
    .
    OpenUrlAbstract/FREE Full Text
  34. ↵
    1. Sibley LD
    (2011) Invasion and intracellular survival by protozoan parasites. Immunol Rev 240(1):72–91
    .
    OpenUrlCrossRefPubMed
  35. ↵
    1. Striepen B, et al.
    (2004) Gene transfer in the evolution of parasite nucleotide biosynthesis. Proc Natl Acad Sci USA 101(9):3154–3159
    .
    OpenUrlAbstract/FREE Full Text
  36. ↵
    1. Efstathiou S,
    2. Kemp S,
    3. Darby G,
    4. Minson AC
    (1989) The role of herpes simplex virus type 1 thymidine kinase in pathogenesis. J Gen Virol 70(Pt 4):869–879
    .
    OpenUrlCrossRefPubMed
  37. ↵
    1. Bakowski MA, et al.
    (2014) Ubiquitin-mediated response to microsporidia and virus infection in C. elegans. PLoS Pathog 10(6):e1004200
    .
    OpenUrlCrossRefPubMed
  38. ↵
    1. Kokoris MS,
    2. Black ME
    (2002) Characterization of herpes simplex virus type 1 thymidine kinase mutants engineered for improved ganciclovir or acyclovir activity. Protein Sci 11(9):2267–2272
    .
    OpenUrlCrossRefPubMed
  39. ↵
    1. Longley DB,
    2. Harkin DP,
    3. Johnston PG
    (2003) 5-fluorouracil: Mechanisms of action and clinical strategies. Nat Rev Cancer 3(5):330–338
    .
    OpenUrlCrossRefPubMed
  40. ↵
    1. Simmons A
    (2002) Clinical manifestations and treatment considerations of herpes simplex virus infection. J Infect Dis 186(Suppl 1):S71–S77
    .
    OpenUrlAbstract/FREE Full Text
  41. ↵
    1. Furman PA,
    2. Barry DW
    (1988) Spectrum of antiviral activity and mechanism of action of zidovudine. An overview. Am J Med 85(2A):176–181
    .
    OpenUrlPubMed
  42. ↵
    1. Power DG,
    2. Kemeny NE
    (2009) The role of floxuridine in metastatic liver disease. Mol Cancer Ther 8(5):1015–1025
    .
    OpenUrlAbstract/FREE Full Text
  43. ↵
    1. Munch-Petersen B,
    2. Cloos L,
    3. Jensen HK,
    4. Tyrsted G
    (1995) Human thymidine kinase 1. Regulation in normal and malignant cells. Adv Enzyme Regul 35:69–89
    .
    OpenUrlCrossRefPubMed
  44. ↵
    1. Hsu C-H,
    2. Liou JY,
    3. Dutschman GE,
    4. Cheng YC
    (2005) Phosphorylation of cytidine, deoxycytidine, and their analog monophosphates by human UMP/CMP kinase is differentially regulated by ATP and magnesium. Mol Pharmacol 67(3):806–814
    .
    OpenUrlAbstract/FREE Full Text
  45. ↵
    1. Briozzo P, et al.
    (1998) Structures of Escherichia coli CMP kinase alone and in complex with CDP: A new fold of the nucleoside monophosphate binding domain and insights into cytosine nucleotide specificity. Structure 6(12):1517–1527
    .
    OpenUrlCrossRefPubMed
  46. ↵
    1. Ray AS, et al.
    (2003) Probing the mechanistic consequences of 5-fluorine substitution on cytidine nucleotide analogue incorporation by HIV-1 reverse transcriptase. Antivir Chem Chemother 14(3):115–125
    .
    OpenUrlAbstract/FREE Full Text
  47. ↵
    1. Mallet LV,
    2. Becq J,
    3. Deschavanne P
    (2010) Whole genome evaluation of horizontal transfers in the pathogenic fungus Aspergillus fumigatus. BMC Genomics 11:171
    .
    OpenUrlCrossRefPubMed
  48. ↵
    1. Campbell MA,
    2. Rokas A,
    3. Slot JC
    (2012) Horizontal transfer and death of a fungal secondary metabolic gene cluster. Genome Biol Evol 4(3):289–293
    .
    OpenUrlAbstract/FREE Full Text
  49. ↵
    1. DeBlasio D,
    2. Wiscaver J
    (2013) SICLE: A high-throughput tool for extracting evolutionary relationships from phylogenetic trees. arXiv:1303.5785
    .
  50. ↵
    1. Li W,
    2. Godzik A
    (2006) Cd-hit: A fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 22(13):1658–1659
    .
    OpenUrlAbstract/FREE Full Text
  51. ↵
    1. Katoh K,
    2. Standley DM
    (2013) MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol Biol Evol 30(4):772–780
    .
    OpenUrlAbstract/FREE Full Text
  52. ↵
    1. Capella-Gutiérrez S,
    2. Silla-Martínez JM,
    3. Gabaldón T
    (2009) trimAl: A tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 25(15):1972–1973
    .
    OpenUrlAbstract/FREE Full Text
  53. ↵
    1. Stamatakis A
    (2006) RAxML-VI-HPC: Maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 22(21):2688–2690
    .
    OpenUrlAbstract/FREE Full Text
  54. ↵
    1. Edgar RC
    (2004) MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 32(5):1792–1797
    .
    OpenUrlAbstract/FREE Full Text
  55. ↵
    1. Shimodaira H,
    2. Hasegawa M
    (1999) Multiple comparisons of log-likelihoods with applications to phylogenetic inference. Mol Biol Evol 16:1114–1116
    .
    OpenUrlCrossRef
  56. ↵
    1. Oliphant T
    (2007) Python for scientific computing. Comput Sci Eng 13:10–20
    .
    OpenUrl
  57. ↵
    1. Presnyak V, et al.
    (2015) Codon optimality is a major determinant of mRNA stability. Cell 160(6):1111–1124
    .
    OpenUrlCrossRefPubMed
  58. ↵
    1. Weir M,
    2. Keeney JB
    (2014) PCR mutagenesis and gap repair in yeast. Methods Mol Biol 1205:29–35
    .
    OpenUrlCrossRefPubMed
  59. ↵
    1. Mortimer RK
    (2000) Evolution and variation of the yeast (Saccharomyces) genome. Genome Res 10(4):403–409
    .
    OpenUrlAbstract/FREE Full Text
  60. ↵
    1. Gietz RD
    (2014) Yeast transformation by the LiAc/SS carrier DNA/PEG method. Methods Mol Biol 1205:1–12
    .
    OpenUrlCrossRefPubMed
View Abstract
Back to top
Article Alerts
Email Article
Citation Tools
Request Permissions
Horizontal transfer in early-diverging fungi
William G. Alexander, Jennifer H. Wisecaver, Antonis Rokas, Chris Todd Hittinger
Proceedings of the National Academy of Sciences Apr 2016, 113 (15) 4116-4121; DOI: 10.1073/pnas.1517242113
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
Proceedings of the National Academy of Sciences: 116 (30)
Current Issue

Submit

Article Classifications

Jump to section

You May Also be Interested in

Chemical evidence suggests ritual use of psychoactive plants used to brew ayahuasca in pre-Columbian Bolivia. Image courtesy of Juan V. Albarracin-Jordan and José M. Capriles.
Hints of ayahuasca in pre-Columbian rituals
Chemical evidence suggests ritual use of psychoactive plants used to brew ayahuasca in pre-Columbian Bolivia.
Image courtesy of Juan Albarracin-Jordan and José M. Capriles.
Changes in subarachnoid and intracerebral spaces, including superior sagittal sinus (1) and ventricles (2), of cosmonauts. Image courtesy of R. Maxine Rühl.
Brain fluid shifts following spaceflight
Volume of cosmonauts’ brain ventricles increased by an average of 12% after spaceflight—a potential mechanism to cope with increased brain fluid volume, a study suggests.
Image courtesy of R. Maxine Rühl.

Similar Articles