Lack of validation of genetic variants associated with anti–tumor necrosis factor therapy response in rheumatoid arthritis: a genome-wide association study replication and meta-analysis

Introduction

In this study, our aim was to elucidate the role of four polymorphisms identified in a prior large genome-wide association study (GWAS) in which the investigators analyzed the responses of patients with rheumatoid arthritis (RA) to treatment with tumor necrosis factor inhibitors (TNFi). The authors of that study reported that the four genetic variants were significantly associated. However, none of the associations reached GWAS significance, and two subsequent studies failed to replicate these associations.

Methods

The four polymorphisms (rs12081765, rs1532269, rs17301249 and rs7305646) were genotyped in a total of 634 TNFi-treated RA patients of Spanish Caucasian origin. Four outcomes were evaluated: changes in the Disease Activity Score in 28 joints (DAS28) after 6 and 12 months of treatment and classification according to the European League Against Rheumatism (EULAR) response criteria at the same time points. Association with DAS28 changes was assessed by linear regression using an additive genetic model. Contingency tables of genotype and allele frequencies between EULAR responder and nonresponder patients were compared. In addition, we combined our data with those of previously reported studies in a meta-analysis including 2,998 RA patients.

Results

None of the four genetic variants showed an association with response to TNFi in any of the four outcomes analyzed in our Spanish patients. In addition, only rs1532269 yielded a suggestive association (P = 0.0033) with the response to TNFi when available data from previous studies were combined in the meta-analysis.

Conclusion

Our data suggest that the rs12081765, rs1532269, rs17301249 and rs7305646 genetic variants do not have a role as genetic predictors of TNFi treatment outcomes.     

A systems genetics approach provides a bridge from discovered genetic variants to biological pathways in rheumatoid arthritis

Genome-wide association studies (GWAS) have yielded novel genetic loci underlying common diseases. We propose a systems genetics approach to utilize these discoveries for better understanding of the genetic architecture of rheumatoid arthritis (RA). Current evidence of genetic associations with RA was sought through PubMed and the NHGRI GWAS catalog. The associations of 15 single nucleotide polymorphisms and HLA-DRB1 alleles were confirmed in 1,287 cases and 1,500 controls of Japanese subjects. Among these, HLA-DRB1 alleles and eight SNPs showed significant associations and all but one of the variants had the same direction of effect as identified in the previous studies, indicating that the genetic risk factors underlying RA are shared across populations. By receiver operating characteristic curve analysis, the area under the curve (AUC) for the genetic risk score based on the selected variants was 68.4%. For seropositive RA patients only, the AUC improved to 70.9%, indicating good but suboptimal predictive ability. A simulation study shows that more than 200 additional loci with similar effect size as recent GWAS findings or 20 rare variants with intermediate effects are needed to achieve AUC = 80.0%. We performed the random walk with restart (RWR) algorithm to prioritize genes for future mapping studies. The performance of the algorithm was confirmed by leave-one-out cross-validation. The RWR algorithm pointed to ZAP70 in the first rank, in which mutation causes RA-like autoimmune arthritis in mice. By applying the hierarchical clustering method to a subnetwork comprising RA-associated genes and top-ranked genes by the RWR, we found three functional modules relevant to RA etiology: “leukocyte activation and differentiation”, “pattern-recognition receptor signaling pathway”, and “chemokines and their receptors”.These results suggest that the systems genetics approach is useful to find directions of future mapping strategies to illuminate biological pathways.     

Prophages integrating into prophages: A mechanism to accumulate type III secretion effector genes and duplicate Shiga toxin-encoding prophages in Escherichia coli

Bacteriophages (or phages) play major roles in the evolution of bacterial pathogens via horizontal gene transfer. Multiple phages are often integrated in a host chromosome as prophages, not only carrying various novel virulence-related genetic determinants into host bacteria but also providing various possibilities for prophage-prophage interactions in bacterial cells. In particular, Escherichia coli strains such as Shiga toxin (Stx)-producing E. coli (STEC) and enteropathogenic E. coli (EPEC) strains have acquired more than 10 prophages (up to 21 prophages), many of which encode type III secretion system (T3SS) effector gene clusters. In these strains, some prophages are present at a single locus in tandem, which is usually interpreted as the integration of phages that use the same attachment (att) sequence. Here, we present phages integrating into T3SS effector gene cluster-associated loci in prophages, which are widely distributed in STEC and EPEC. Some of the phages integrated into prophages are Stx-encoding phages (Stx phages) and have induced the duplication of Stx phages in a single cell. The identified attB sequences in prophage genomes are apparently derived from host chromosomes. In addition, two or three different attB sequences are present in some prophages, which results in the generation of prophage clusters in various complex configurations. These phages integrating into prophages represent a medically and biologically important type of inter-phage interaction that promotes the accumulation of T3SS effector genes in STEC and EPEC, the duplication of Stx phages in STEC, and the conversion of EPEC to STEC and that may be distributed in other types of E. coli strains as well as other prophage-rich bacterial species.     

Combinations of Susceptibility Genes Are Associated with Higher Risk for Multiple Sclerosis and Imply Disease Course Specificity

Multiple sclerosis (MS) is a chronic autoimmune disease of the central nervous system that predominantly affects young adults. The genetic contributions to this multifactorial disease were underscored by a genome wide association study (GWAS) conducted by the International Multiple Sclerosis Genetic Consortium in a multinational cohort prompting the discovery of 57 non-MHC MS-associated common genetic variants. Hitherto, few of these newly reported variants have been replicated in larger independent patient cohorts. We genotyped a cohort of 1033 MS patients and 644 healthy controls with a consistent genetic background for the 57 non-MHC variants reported to be associated with MS by the first large GWAS as well as the HLA DRB1*1501 tagging SNP rs3135388. We robustly replicated three of the 57 non-MHC reported MS-associated single nucleotide polymorphisms (SNPs). In addition, our study revealed several genotype-genotype combinations with an evidently higher degree of disease association than the genotypes of the single SNPs. We further correlated well-defined clinical phenotypes, i.e. ataxia, visual impairment due to optic neuritis and paresis with single SNPs and genotype combinations, and identified several associations. The results may open new avenues for clinical implications of the MS associated genetic variants reported from large GWAS.     

Isolation of an Aptamer that Binds Specifically to E. coli

Escherichia coli is a bacterial species found ubiquitously in the intestinal flora of animals, although pathogenic variants cause major public health problems. Aptamers are short oligonucleotides that bind to targets with high affinity and specificity, and have great potential for use in diagnostics and therapy. We used cell-based Systematic Evolution of Ligands by EXponential enrichment (cell-SELEX) to isolate four single stranded DNA (ssDNA) aptamers that bind strongly to E. coli cells (ATCC generic strain 25922), with Kd values in the nanomolar range. Fluorescently labeled aptamers label the surface of E. coli cells, as viewed by fluorescent microscopy. Specificity tests with twelve different bacterial species showed that one of the aptamers–called P12-31—is highly specific for E. coli. Importantly, this aptamer binds to Meningitis/sepsis associated E. coli (MNEC) clinical isolates, and is the first aptamer described with potential for use in the diagnosis of MNEC-borne pathologies.     

Epidemiology and Pathogenesis of Staphylococcus Bloodstream Infections in Humans: a Review

Staphylococci are among the most frequent human microbiota components associated with the high level of bloodstream infection (BSI) episodes. In predisposed patients, there is a high risk of transformation of BSI episodes to sepsis. Both bacterial and host factors are crucial for the outcomes of BSI and sepsis. The highest rates of BSI episodes were reported in Africa, where these infections were up to twice as high as the European rates. However, there remains a great need to analyze African data for comprehensive quantification of staphylococcal BSI prevalence. The lowest rates of BSI exist in Australia. Asian, European, and North American data showed similar frequency values. Worldwide analysis indicated that both Staphylococcus aureus and coagulase-negative staphylococci (CoNS) are the most frequent BSI agents. In the second group, the most prevalent species was Staphylococcus epidermidis, although CoNS were not identified at the species level in many studies. The lack of a significant worldwide decrease in BSI episodes indicates a great need to implement standardized diagnostic methods and research etiological factors using advanced genetic methods.     

Deep resequencing identifies candidate functional genes in leprosy GWAS loci

Leprosy is the second most prevalent mycobacterial disease globally. Despite the existence of an effective therapy, leprosy incidence has consistently remained above 200,000 cases per year since 2010. Numerous host genetic factors have been identified for leprosy that contribute to the persistently high case numbers. In the past decade, genetic epidemiology approaches, including genome-wide association studies (GWAS), identified more than 30 loci contributing to leprosy susceptibility. However, GWAS loci commonly encompass multiple genes, which poses a challenge to define causal candidates for each locus. To address this problem, we hypothesized that genes contributing to leprosy susceptibility differ in their frequencies of rare protein-altering variants between cases and controls. Using deep resequencing we assessed protein-coding variants for 34 genes located in GWAS or linkage loci in 555 Vietnamese leprosy cases and 500 healthy controls. We observed 234 nonsynonymous mutations in the targeted genes. A significant depletion of protein-altering variants was detected for the IL18R1 and BCL10 genes in leprosy cases. The IL18R1 gene is clustered with IL18RAP and IL1RL1 in the leprosy GWAS locus on chromosome 2q12.1. Moreover, in a recent GWAS we identified an HLA-independent signal of association with leprosy on chromosome 6p21. Here, we report amino acid changes in the CDSN and PSORS1C2 genes depleted in leprosy cases, indicating them as candidate genes in the chromosome 6p21 locus. Our results show that deep resequencing can identify leprosy candidate susceptibility genes that had been missed by classic linkage and association approaches.     

Insight into mechanisms of pig lncRNA FUT3-AS1 regulating E. coli F18-bacterial diarrhea

Escherichia coli F18 is a common conditional pathogen that is associated with a variety of infections in humans and animals. LncRNAs have emerged as critical players in pathogen infection, but their role in the resistance of the host to bacterial diarrhea remains unknown. Here, we used piglets as animal model and identified an antisense lncRNA termed FUT3-AS1 as a host regulator related to E. coli F18 infection by RNA sequencing. Downregulation of FUT3-AS1 expression contributed to the enhancement of E. coli F18 resistance in IPEC-J2 cells. FUT3-AS1 knockdown reduced FUT3 expression via decreasing the H4K16ac level of FUT3 promoter. Besides, the FUT3-AS1/miR-212 axis could act as a competing endogenous RNA to regulate FUT3 expression. Functional analysis demonstrated that target FUT3 plays a vital role in the resistance of IPEC-J2 cells to E. coli F18 invasion. A Fut3-knockout mice model was established and Fut3-knockout mice obviously improved the ability of resistance to bacterial diarrhea. Interestingly, FUT3 could enhance E. coli F18 susceptibility by activating glycosphingolipid biosynthesis and toll-like receptor signaling which are related to receptor formation and immune response, respectively. In summary, we have identified a novel biomarker FUT3-AS1 that modulates E. coli F18 susceptibility via histone H4 modifications or miR-212/FUT3 axis, which will provide theoretical guidance to develop novel strategies for combating bacterial diarrhea in piglets.     

Connecting the Dots: Potential of Data Integration to Identify Regulatory SNPs in Late-Onset Alzheimer's Disease GWAS Findings

Late-onset Alzheimer''s disease (LOAD) is a multifactorial disorder with over twenty loci associated with disease risk. Given the number of genome-wide significant variants that fall outside of coding regions, it is possible that some of these variants alter some function of gene expression rather than tagging coding variants that alter protein structure and/or function. RegulomeDB is a database that annotates regulatory functions of genetic variants. In this study, we utilized RegulomeDB to investigate potential regulatory functions of lead single nucleotide polymorphisms (SNPs) identified in five genome-wide association studies (GWAS) of risk and age-at onset (AAO) of LOAD, as well as SNPs in LD (r²≥0.80) with the lead GWAS SNPs. Of a total 614 SNPs examined, 394 returned RegulomeDB scores of 1–6. Of those 394 variants, 34 showed strong evidence of regulatory function (RegulomeDB score <3), and only 3 of them were genome-wide significant SNPs (ZCWPW1/rs1476679, CLU/rs1532278 and ABCA7/rs3764650). This study further supports the assumption that some of the non-coding GWAS SNPs are true associations rather than tagged associations and demonstrates the application of RegulomeDB to GWAS data.     

Pleiotropy informed adaptive association test of multiple traits using genome‐wide association study summary data

Genetic variants associated with disease outcomes can be used to develop personalized treatment. To reach this precision medicine goal, hundreds of large‐scale genome‐wide association studies (GWAS) have been conducted in the past decade to search for promising genetic variants associated with various traits. They have successfully identified tens of thousands of disease‐related variants. However, in total these identified variants explain only part of the variation for most complex traits. There remain many genetic variants with small effect sizes to be discovered, which calls for the development of (a) GWAS with more samples and more comprehensively genotyped variants, for example, the NHLBI Trans‐Omics for Precision Medicine (TOPMed) Program is planning to conduct whole genome sequencing on over 100 000 individuals; and (b) novel and more powerful statistical analysis methods. The current dominating GWAS analysis approach is the “single trait” association test, despite the fact that many GWAS are conducted in deeply phenotyped cohorts including many correlated and well‐characterized outcomes, which can help improve the power to detect novel variants if properly analyzed, as suggested by increasing evidence that pleiotropy, where a genetic variant affects multiple traits, is the norm in genome‐phenome associations. We aim to develop pleiotropy informed powerful association test methods across multiple traits for GWAS. Since it is generally very hard to access individual‐level GWAS phenotype and genotype data for those existing GWAS, due to privacy concerns and various logistical considerations, we develop rigorous statistical methods for pleiotropy informed adaptive multitrait association test methods that need only summary association statistics publicly available from most GWAS. We first develop a pleiotropy test, which has powerful performance for truly pleiotropic variants but is sensitive to the pleiotropy assumption. We then develop a pleiotropy informed adaptive test that has robust and powerful performance under various genetic models. We develop accurate and efficient numerical algorithms to compute the analytical P‐value for the proposed adaptive test without the need of resampling or permutation. We illustrate the performance of proposed methods through application to joint association test of GWAS meta‐analysis summary data for several glycemic traits. Our proposed adaptive test identified several novel loci missed by individual trait based GWAS meta‐analysis. All the proposed methods are implemented in a publicly available R package.     

Complete Genome Sequence and Comparative Analysis of the Wild-type Commensal Escherichia coli Strain SE11 Isolated from a Healthy Adult

We sequenced and analyzed the genome of a commensal Escherichia coli (E. coli) strain SE11 (O152:H28) recently isolated from feces of a healthy adult and classified into E. coli phylogenetic group B1. SE11 harbored a 4.8 Mb chromosome encoding 4679 protein-coding genes and six plasmids encoding 323 protein-coding genes. None of the SE11 genes had sequence similarity to known genes encoding phage- and plasmid-borne virulence factors found in pathogenic E. coli strains. The comparative genome analysis with the laboratory strain K-12 MG1655 identified 62 poorly conserved genes between these two non-pathogenic strains and 1186 genes absent in MG1655. These genes in SE11 were mostly encoded in large insertion regions on the chromosome or in the plasmids, and were notably abundant in genes of fimbriae and autotransporters, which are cell surface appendages that largely contribute to the adherence ability of bacteria to host cells and bacterial conjugation. These data suggest that SE11 may have evolved to acquire and accumulate the functions advantageous for stable colonization of intestinal cells, and that the adhesion-associated functions are important for the commensality of E. coli in human gut habitat.Key words: Escherichia coli, commensal, human gut, genome sequencing     

Tuberculosis infection and lung adenocarcinoma: Mendelian randomization and pathway analysis of genome-wide association study data from never-smoking Asian women

We investigated whether genetic susceptibility to tuberculosis (TB) influences lung adenocarcinoma development among never-smokers using TB genome-wide association study (GWAS) results within the Female Lung Cancer Consortium in Asia. Pathway analysis with the adaptive rank truncated product method was used to assess the association between a TB-related gene-set and lung adenocarcinoma using GWAS data from 5512 lung adenocarcinoma cases and 6277 controls. The gene-set consisted of 31 genes containing known/suggestive associations with genetic variants from previous TB-GWAS. Subsequently, we followed-up with Mendelian Randomization to evaluate the association between TB and lung adenocarcinoma using three genome-wide significant variants from previous TB-GWAS in East Asians. The TB-related gene-set was associated with lung adenocarcinoma (p = 0.016). Additionally, the Mendelian Randomization showed an association between TB and lung adenocarcinoma (OR = 1.31, 95% CI: 1.03, 1.66, p = 0.027). Our findings support TB as a causal risk factor for lung cancer development among never-smoking Asian women.     

Evaluation of European Schizophrenia GWAS Loci in Asian Populations via Comprehensive Meta-Analyses

Schizophrenia is a severe and highly heritable neuropsychiatric disorder. Recent genetic analyses including genome-wide association studies (GWAS) have implicated multiple genome-wide significant variants for schizophrenia among European populations. However, many of these risk variants were not largely validated in other populations of different ancestry such as Asians. To validate whether these European GWAS significant loci are associated with schizophrenia in Asian populations, we conducted a systematic literature search and meta-analyses on 19 single nucleotide polymorphisms (SNPs) in Asian populations by combining all available case-control and family-based samples, including up to 30,000 individuals. We employed classical fixed (or random) effects inverse variance weighted methods to calculate summary odds ratios (ORs) and 95 % confidence intervals (CIs). Among the 19 GWAS loci, we replicated the risk associations of nine markers (e.g., SNPs at VRK2, ITIH3/4, NDST3, NOTCH4) surpassing significance level (two-tailed P?<?0.05), and three additional SNPs in MIR137 and ZNF804A also showed trend associations (one-tailed P?<?0.05). These risk associations are in the same directions of allelic effects between Asian replication samples and initial European GWAS findings, and the successful replications of these GWAS loci in a different ethnic group provide stronger evidence for their clinical associations with schizophrenia. Further studies, focusing on the molecular mechanisms of these GWAS significant loci, will become increasingly important for understanding of the pathogenesis to schizophrenia.     

A Phylogeny-aware GWAS Framework to Correct for Heritable Pathogen Effects on Infectious Disease Traits

Infectious diseases are particularly challenging for genome-wide association studies (GWAS) because genetic effects from two organisms (pathogen and host) can influence a trait. Traditional GWAS assume individual samples are independent observations. However, pathogen effects on a trait can be heritable from donor to recipient in transmission chains. Thus, residuals in GWAS association tests for host genetic effects may not be independent due to shared pathogen ancestry. We propose a new method to estimate and remove heritable pathogen effects on a trait based on the pathogen phylogeny prior to host GWAS, thus restoring independence of samples. In simulations, we show this additional step can increase GWAS power to detect truly associated host variants when pathogen effects are highly heritable, with strong phylogenetic correlations. We applied our framework to data from two different host–pathogen systems, HIV in humans and X. arboricola in A. thaliana. In both systems, the heritability and thus phylogenetic correlations turn out to be low enough such that qualitative results of GWAS do not change when accounting for the pathogen shared ancestry through a correction step. This means that previous GWAS results applied to these two systems should not be biased due to shared pathogen ancestry. In summary, our framework provides additional information on the evolutionary dynamics of traits in pathogen populations and may improve GWAS if pathogen effects are highly phylogenetically correlated amongst individuals in a cohort.     

Small Heat-Shock Proteins,IbpAB, Protect Non-Pathogenic Escherichia coli from Killing by Macrophage-Derived Reactive Oxygen Species

Many intracellular bacterial pathogens possess virulence factors that prevent detection and killing by macrophages. However, similar virulence factors in non-pathogenic bacteria are less well-characterized and may contribute to the pathogenesis of chronic inflammatory conditions such as Crohn’s disease. We hypothesize that the small heat shock proteins IbpAB, which have previously been shown to reduce oxidative damage to proteins in vitro and be upregulated in luminal non-pathogenic Escherichia strain NC101 during experimental colitis in vivo, protect commensal E. coli from killing by macrophage-derived reactive oxygen species (ROS). Using real-time PCR, we measured ibpAB expression in commensal E. coli NC101 within wild-type (wt) and ROS-deficient (gp91phox^-/-) macrophages and in NC101 treated with the ROS generator paraquat. We also quantified survival of NC101 and isogenic mutants in wt and gp91phox^-/- macrophages using gentamicin protection assays. Similar assays were performed using a pathogenic E. coli strain O157:H7. We show that non-pathogenic E. coli NC101inside macrophages upregulate ibpAB within 2 hrs of phagocytosis in a ROS-dependent manner and that ibpAB protect E. coli from killing by macrophage-derived ROS. Moreover, we demonstrate that ROS-induced ibpAB expression is mediated by the small E. coli regulatory RNA, oxyS. IbpAB are not upregulated in pathogenic E. coli O157:H7 and do not affect its survival within macrophages. Together, these findings indicate that ibpAB may be novel virulence factors for certain non-pathogenic E. coli strains.     

Genome-wide association studies reveal the role of polymorphisms affecting factor H binding protein expression in host invasion by Neisseria meningitidis

Many invasive bacterial diseases are caused by organisms that are ordinarily harmless components of the human microbiome. Effective interventions against these microbes require an understanding of the processes whereby symbiotic or commensal relationships transition into pathology. Here, we describe bacterial genome-wide association studies (GWAS) of Neisseria meningitidis, a common commensal of the human respiratory tract that is nevertheless a leading cause of meningitis and sepsis. An initial GWAS discovered bacterial genetic variants, including single nucleotide polymorphisms (SNPs), associated with invasive meningococcal disease (IMD) versus carriage in several loci across the meningococcal genome, encoding antigens and other extracellular components, confirming the polygenic nature of the invasive phenotype. In particular, there was a significant peak of association around the fHbp locus, encoding factor H binding protein (fHbp), which promotes bacterial immune evasion of human complement by recruiting complement factor H (CFH) to the meningococcal surface. The association around fHbp with IMD was confirmed by a validation GWAS, and we found that the SNPs identified in the validation affected the 5’ region of fHbp mRNA, altering secondary RNA structures, thereby increasing fHbp expression and enhancing bacterial escape from complement-mediated killing. This finding is consistent with the known link between complement deficiencies and CFH variation with human susceptibility to IMD. These observations demonstrate the importance of human and bacterial genetic variation across the fHbp:CFH interface in determining IMD susceptibility, the transition from carriage to disease.     

Repositioning of Memantine as a Potential Novel Therapeutic Agent against Meningitic E. coli–Induced Pathogenicities through Disease-Associated Alpha7 Cholinergic Pathway and RNA Sequencing-Based Transcriptome Analysis of Host Inflammatory Responses

Neonatal sepsis and meningitis (NSM) remains a leading cause worldwide of mortality and morbidity in newborn infants despite the availability of antibiotics over the last several decades. E. coli is the most common gram-negative pathogen causing NSM. Our previous studies show that α7 nicotinic receptor (α7 nAChR), an essential regulator of inflammation, plays a detrimental role in the host defense against NSM. Despite notable successes, there still exists an unmet need for new effective therapeutic approaches to treat this disease. Using the in vitro/in vivo models of the blood-brain barrier (BBB) and RNA-seq, we undertook a drug repositioning study to identify unknown antimicrobial activities for known drugs. We have demonstrated for the first time that memantine (MEM), a FDA-approved drug for treatment of Alzheimer’s disease, could very efficiently block E. coli-caused bacteremia and meningitis in a mouse model of NSM in a manner dependent on α7 nAChR. MEM was able to synergistically enhance the antibacterial activity of ampicillin in HBMEC infected with E. coli K1 (E44) and in neonatal mice with E44-caused bacteremia and meningitis. Differential gene expression analysis of RNA-Seq data from mouse BMEC infected with E. coli K1 showed that several E44-increased inflammatory factors, including IL33, IL18rap, MMP10 and Irs1, were significantly reduced by MEM compared to the infected cells without drug treatment. MEM could also significantly up-regulate anti-inflammatory factors, including Tnfaip3, CISH, Ptgds and Zfp36. Most interestingly, these factors may positively and negatively contribute to regulation of NF-κB, which is a hallmark feature of bacterial meningitis. Furthermore, we have demonstrated that circulating BMEC (cBMEC) are the potential novel biomarkers for NSM. MEM could significantly reduce E44-increased blood level of cBMEC in mice. Taken together, our data suggest that memantine can efficiently block host inflammatory responses to bacterial infection through modulation of both inflammatory and anti-inflammatory pathways.     

Characterisation of novel functionality within the Blastocystis tryptophanase gene

In recent years, the human gut microbiome has been recognised to play a pivotal role in the health of the host. Intestinal homeostasis relies on this intricate and complex relationship between the gut microbiota and the human host. While much effort and attention has been placed on the characterization of the organisms that inhabit the gut microbiome, the complex molecular cross-talk between the microbiota could also exert an effect on gastrointestinal conditions. Blastocystis is a single-cell eukaryotic parasite of emerging interest, as its beneficial or pathogenic role in the microbiota has been a subject of contention even to-date. In this study, we assessed the function of the Blastocystis tryptophanase gene (BhTnaA), which was acquired by horizontal gene transfer and likely to be of bacterial origin within Blastocystis. Bioinformatic analysis and phylogenetic reconstruction revealed distinct divergence of BhTnaA versus known bacterial homologs. Despite sharing high homology with the E. coli tryptophanase gene, we show that Blastocystis does not readily convert tryptophan into indole. Instead, BhTnaA preferentially catalyzes the conversion of indole to tryptophan. We also show a direct link between E. coli and Blastocystis tryptophan metabolism: In the presence of E. coli, Blastocystis ST7 is less able to metabolise indole to tryptophan. This study examines the potential for functional variation in horizontally-acquired genes relative to their canonical counterparts, and identifies Blastocystis as a possible producer of tryptophan within the gut.