An in silico pipeline to filter the Toxoplasma gondii proteome for proteins that could traffic to the host cell nucleus and influence host cell epigenetic regulation.

Toxoplasma gondii uses epigenetic mechanisms to regulate both endogenous and host cell gene expression. To identify genes with putative epigenetic functions, we developed an in silico pipeline to interrogate the T. gondii proteome of 8313 proteins. Step 1 employs PredictNLS and NucPred to identify genes predicted to target eukaryotic nuclei. Step 2 uses GOLink to identify proteins of epigenetic function based on Gene Ontology terms. This resulted in 611 putative nuclear localised proteins with predicted epigenetic functions. Step 3 filtered for secretory proteins using SignalP, SecretomeP, and experimental data. This identified 57 of the 611 putative epigenetic proteins as likely to be secreted. The pipeline is freely available online, uses open access tools and software with user-friendly Perl scripts to automate and manage the results, and is readily adaptable to undertake any such in silico search for genes contributing to particular functions.

An in silico pipeline to filter the Toxoplasma gondii proteome for proteins that could traffic to the host cell nucleus and influence host cell epigenetic regulation

Toxoplasma gondii uses epigenetic mechanisms to regulate both endogenous and host cell gene expression. To identify genes with putative epigenetic functions, we developed an in silico pipeline to interrogate the T. gondii proteome of 8313 proteins. Step 1 employs PredictNLS and NucPred to identify genes predicted to target eukaryotic nuclei. Step 2 uses GOLink to identify proteins of epigenetic function based on Gene Ontology terms. This resulted in 611 putative nuclear localised proteins with predicted epigenetic functions. Step 3 filtered for secretory proteins using SignalP, SecretomeP, and experimental data. This identified 57 of the 611 putative epigenetic proteins as likely to be secreted. The pipeline is freely available online, uses open access tools and software with user-friendly Perl scripts to automate and manage the results, and is readily adaptable to undertake any such in silico search for genes contributing to particular functions.

Analysis of expression of FLI1 and MMP1 in American cutaneous leishmaniasis caused by Leishmania braziliensis infection.

FLI1 (Friend leukemia virus integration 1) and IL6 (interleukin 6; IL-6) are associated with Leishmania braziliensis susceptibility. Cutaneous lesions show exaggerated matrix metalloproteinase 1 (MMP1). In other skin diseases, FLI1 promoter methylation reduces FLI1 expression, and low FLI1 down-regulates MMP1. IL-6 increases FLI1 expression. We hypothesized that epigenetic regulation of FLI1 in cutaneous leishmaniasis, together with IL-6, might determine MMP1 expression. While generally low (<10%), percent FLI1 promoter methylation was lower (P=0.001) in lesion biopsies than normal skin. Contrary to expectation, a strong positive correlation occurred between FLI1 methylation and gene expression in lesions (r=0.98, P=0.0005) and in IL-6-treated L. braziliensis-infected macrophages (r=0.99, P=0.0004). In silico analysis of the FLI1 promoter revealed co-occurring active H3K27ac and repressive DNA methylation marks to enhance gene expression. FLI1 expression was enhanced between 3 and 24hour post infection in untreated (P=0.0002) and IL-6-treated (P=0.028) macrophages. MMP1 was enhanced in lesion biopsies (P=0.0002), induced (P=0.007) in infected macrophages, but strongly inhibited by IL-6. No correlations occurred between FLI1 and MMP1 expression in lesions or infected macrophages (with/without IL-6). We conclude that MMP1 is regulated by factors other than FLI1, and that the influence of IL-6 on MMP1 was independent of its effect on FLI1.

The -308 bp TNF gene polymorphism influences tumor necrosis factor expression in leprosy patients in Bahia State, Brazil.

Leprosy or Hansen's disease is a debilitating chronic granulomatous disease caused by Mycobacterium leprae, with high incidence and prevalence in Brazil. The -308 bp G/A single nucleotide polymorphism (SNP rs1800629) in the tumor necrosis factor (TNF) gene promoter is a proposed risk factor for leprosy. In Brazil, Northern India, Egypt and Nepal, the common G allele was associated with leprosy. In Eastern India, Thailand and Malawi the minor A allele was the risk factor. Allele A was previously associated with high TNF. We genotyped rs1800629 in 326 leprosy cases from Bahia State, Brazil, including 72 paucibacillary (PB) and 47 multibacillary (MB) without reactions, and 69 reversal reaction (RR) and 78 erythema nodosum leprosum (ENL) with reactions. Logistic regression was used to compare patient groups with 331 healthy controls. Relative TNF mRNA was determined in peripheral blood leukocytes by QRTPCR, and serum TNF levels measured by ELISA. We found that TNF mRNA expression was higher (P=0.03) in leprosy patients compared to endemic controls, but did not differ significantly between clinical subgroups. Carriage of the minor A allele was associated (P=0.003) with low TNF mRNA across leprosy patients. Nevertheless, we found no evidence for either allele at this SNP as a risk factor for leprosy per se (OR=1.12, 95% CI 0.79-1.60, P=0.52), PB (OR=0.99, 95% CI 0.54-1.81, P=0.97), MB (OR=0.86, 95% CI 0.40-1.83, P=0.70), RR (OR=1.37, 95% CI 0.79-2.38, P=0.27) or ENL (OR=0.76, 95% CI 0.40-1.45, P=0.42) when compared to endemic controls. Further studies are required to determine whether the influence of the minor A allele on TNF mRNA levels determines response to treatment, particularly in the context of ENL reaction treatment with anti-TNF therapies and RR reactions where treatment with prednisolone is known to reduce TNF levels. Our findings contribute to understanding TNF as an important determinant of leprosy immunopathology in Brazil.

Wound healing genes and susceptibility to cutaneous leishmaniasis in Brazil.

Leishmania braziliensis causes cutaneous (CL) and mucosal (ML) leishmaniasis. In the mouse, Fli1 was identified as a gene influencing enhanced wound healing and resistance to CL caused by Leishmania major. Polymorphism at FLI1 is associated with CL caused by L. braziliensis in humans, with an inverse association observed for ML disease. Here we extend the analysis to look at other wound healing genes, including CTGF, TGFB1, TGFBR1/2, SMADS 2/3/4/7 and FLII, all functionally linked along with FLI1 in the TGF beta pathway. Haplotype tagging single nucleotide polymorphisms (tag-SNPs) were genotyped using Taqman technology in 325 nuclear families (652 CL cases; 126 ML cases) from Brazil. Robust case-pseudocontrol (CPC) conditional logistic regression analysis showed associations between CL and SNPs at CTGF (SNP rs6918698; CC genotype; OR 1.67; 95%CI 1.10-2.54; P=0.016), TGFBR2 (rs1962859; OR 1.50; 95%CI 1.12-1.99; P=0.005), SMAD2 (rs1792658; OR 1.57; 95%CI 1.04-2.38; P=0.03), SMAD7 (rs4464148; AA genotype; OR 2.80; 95%CI 1.00-7.87; P=0.05) and FLII (rs2071242; OR 1.60; 95%CI 1.14-2.24; P=0.005), and between ML and SNPs at SMAD3 (rs1465841; OR 2.15; 95%CI 1.13-4.07; P=0.018) and SMAD7 (rs2337107; TT genotype; OR 3.70; 95%CI 1.27-10.7; P=0.016). Stepwise logistic regression analysis showed that all SNPs associated with CL at FLI1, CTGF, TGFBR2, and FLII showed independent effects from each other, but SNPs at SMAD2 and SMAD7 did not add independent effects to SNPs from other genes. These results suggest that TGFß signalling via SMAD2 is important in directing events that contribute to CL, whereas signalling via SMAD3 is important in ML. Both are modulated by the inhibitory SMAD7 that acts upstream of SMAD2 and SMAD3 in this signalling pathway. Along with the published FLI1 association, these data further contribute to the hypothesis that wound healing processes are important determinants of pathology associated with cutaneous forms of leishmaniasis.

The -2518bp promoter polymorphism at CCL2/MCP1 influences susceptibility to mucosal but not localized cutaneous leishmaniasis in Brazil.

Mucosal leishmaniasis (ML) follows localized cutaneous leishmaniasis (CL) caused by Leishmania braziliensis. Proinflammatory responses mediate CL self-healing but are exaggerated in ML. Proinflammatory monocyte chemoattractant protein 1 (MCP-1; encoded by CCL2) is associated with CL. We explore its role in CL/ML through analysis of the regulatory CCL2 -2518bp promoter polymorphism in CL/ML population samples and families from Brazil. Genotype frequencies were compared among ML/CL cases and control groups using logistic regression and the family-based association test (FBAT). MCP-1 was measured in plasma and macrophages. The GG recessive genotype at CCL2 -2518bp was more common in patients with ML (N=67) than in neighborhood control (NC; N=60) subjects (OR 1.78; 95% CI 1.01-3.14; P=0.045), than in NC combined with leishmanin skin-test positive (N=60) controls (OR 4.40; 95% CI 1.42-13.65; P=0.010), and than in controls combined with CL (N=60) patients (OR 2.78; 95% CI 1.13-6.85; P=0.045). No associations were observed for CL compared to any groups. FBAT (91 ML and 223 CL cases in families) confirmed recessive association of ML with allele G (Z=2.679; P=0.007). Higher levels of MCP-1 occurred in plasma (P=0.03) and macrophages (P<0.0001) from GG compared to AA individuals. These results suggest that high MCP-1 increases risk of ML.

CXCR1 and SLC11A1 polymorphisms affect susceptibility to cutaneous leishmaniasis in Brazil: a case-control and family-based study.

BACKGROUND: L. braziliensis causes cutaneous (CL) and mucosal (ML) leishmaniasis. Wound healing neutrophil (PMN) and macrophage responses made following the bite of the vector sand fly contribute to disease progression in mice. To look at the interplay between PMN and macrophages in disease progression in humans we asked whether polymorphisms at genes that regulate their infiltration or function are associated with different clinical phenotypes. Specifically, CXCR1 (IL8RA) and CXCR2 (IL8RB) are receptors for chemokines that attract PMN to inflammatory sites. They lie 30-260 kb upstream of SLC11A1, a gene known primarily for its role in regulating macrophage activation, resistance to leishmaniasis, and wound healing responses in mice, but also known to be expressed in PMN, macrophages and dendritic cells. METHODS: Polymorphic variants at CXCR1, CXCR2 and SLC11A1 were analysed using Taqman or ABI fragment separation technologies in cases (60 CL; 60 ML), unrelated controls (n = 120), and multicase families (104 nuclear families; 88 ML, 250 CL cases) from Brazil. Logistic regression analysis, family-based association testing (FBAT) and haplotype analysis (TRANSMIT) were performed. RESULTS: Case-control analysis showed association between the common C allele (OR 2.38; 95% CI 1.23-4.57; P = 0.009) of CXCR1_rs2854386 and CL, supported by family-based (FBAT; Z score 2.002; P = 0.045) analysis (104 nuclear families; 88 ML, 250 CL cases). ML associated with the rarer G allele (Z score 1.999; P = 0.046). CL associated with a 3' insertion/deletion polymorphism at SLC11A1 (Z score 2.549; P = 0.011). CONCLUSIONS: The study supports roles for CXCR1 and SLC11A1 in the outcome of L. braziliensis infection in humans. Slc11a1 does not influence cutaneous lesion development following needle injection of Leishmania in mice, suggesting that its role here might relate to the action of PMN, macrophage and/or dendritic cells in the wound healing response to the sand fly bite. Together with the CXCR1 association, the data are consistent with hypotheses relating to the possible role of PMN in initiation of a lesion following the delivery of parasites via the sand fly bite. Association of ML with the rare derived G allele suggests that PMN also have an important positive role to play in preventing this form of the disease.

Candidate gene analysis of ocular toxoplasmosis in Brazil: evidence for a role for toll-like receptor 9 (TLR9)

Toxoplasma gondii infection is an important mediator of ocular disease in Brazil more frequently than reported from elsewhere. Infection and pathology are characterized by a strong proinflammatory response which in mice is triggered by interaction of the parasite with the toll-like receptor (TLR)/MyD88 pathway. A powerful way to identify the role of TLRs in humans is to determine whether polymorphisms at these loci influence susceptibility to T. gondii-mediated pathologies. Here we report on a small family-based study (60 families; 68 affected offspring) undertaken in Brazil which was powered for large effect sizes using single nucleotide polymorphisms with minor alleles frequencies > 0.3. Of markers in TLR2, TLR5 and TLR9 that met these criteria, we found an association Family Based Association Tests [(FBAT) Z score = 4.232; p = 1.5 x 10-5; p corrected = 1.2 x 10-4] between the C allele (frequency = 0.424; odds ratio = 7; 95 percent confidence interval 1.6-30.8) of rs352140 at TLR9 and toxoplasmic retinochoroiditis in Brazil. This supports the hypothesis that direct interaction between T. gondii and TLR9 may trigger proinflammatory responses that lead to severe pathologies such as the ocular disease that is associated with this infection in Brazil.

Host genetic and epigenetic factors in toxoplasmosis.

Analysing human genetic variation provides a powerful tool in understanding risk factors for disease. Toxoplasma gondii acquired by the mother can be transmitted to the fetus. Infants with the most severe clinical signs in brain and eye are those infected early in pregnancy when fetal immunity is least well developed. Genetic analysis could provide unique insight into events in utero that are otherwise difficult to determine. We tested the hypothesis that propensity for T. gondii to cause eye disease is associated with genes previously implicated in congenital or juvenile onset ocular disease. Using mother-child pairs from Europe (EMSCOT) and child/parent trios from North America (NCCCTS), we demonstrated that ocular and brain disease in congenital toxoplasmosis associate with polymorphisms in ABCA4 encoding ATP-binding cassette transporter, subfamily A, member 4 previously associated with juvenile onset retinal dystrophies including Stargardt's disease. Polymorphisms at COL2A1 encoding type II collagen, previously associated with Stickler syndrome, associated only with ocular disease in congenital toxoplasmosis. Experimental studies showed that both ABCA4 and COL2A1 show isoform-specific epigenetic modifications consistent with imprinting, which provided an explanation for the patterns of inheritance observed. These genetic and epigenetic risk factors provide unique insight into molecular pathways in the pathogenesis of disease.

Host genetic and epigenetic factors in toxoplasmosis

Analysing human genetic variation provides a powerful tool in understanding risk factors for disease. Toxoplasma gondii acquired by the mother can be transmitted to the fetus. Infants with the most severe clinical signs in brain and eye are those infected early in pregnancy when fetal immunity is least well developed. Genetic analysis could provide unique insight into events in utero that are otherwise difficult to determine. We tested the hypothesis that propensity for T. gondii to cause eye disease is associated with genes previously implicated in congenital or juvenile onset ocular disease. Using mother-child pairs from Europe (EMSCOT) and child/parent trios from North America (NCCCTS), we demonstrated that ocular and brain disease in congenital toxoplasmosis associate with polymorphisms in ABCA4 encoding ATP-binding cassette transporter, subfamily A, member 4 previously associated with juvenile onset retinal dystrophies including Stargardt's disease. Polymorphisms at COL2A1 encoding type II collagen, previously associated with Stickler syndrome, associated only with ocular disease in congenital toxoplasmosis. Experimental studies showed that both ABCA4 and COL2A1 show isoform-specific epigenetic modifications consistent with imprinting, which provided an explanation for the patterns of inheritance observed. These genetic and epigenetic risk factors provide unique insight into molecular pathways in the pathogenesis of disease.

Candidate gene analysis of ocular toxoplasmosis in Brazil: evidence for a role for toll-like receptor 9 (TLR9).

Toxoplasma gondii infection is an important mediator of ocular disease in Brazil more frequently than reported from elsewhere. Infection and pathology are characterized by a strong proinflammatory response which in mice is triggered by interaction of the parasite with the toll-like receptor (TLR)/MyD88 pathway. A powerful way to identify the role of TLRs in humans is to determine whether polymorphisms at these loci influence susceptibility to T. gondii-mediated pathologies. Here we report on a small family-based study (60 families; 68 affected offspring) undertaken in Brazil which was powered for large effect sizes using single nucleotide polymorphisms with minor alleles frequencies > 0.3. Of markers in TLR2, TLR5 and TLR9 that met these criteria, we found an association Family Based Association Tests [(FBAT) Z score = 4.232; p = 1.5 x 10-5; p corrected = 1.2 x 10-4] between the C allele (frequency = 0.424; odds ratio = 7; 95% confidence interval 1.6-30.8) of rs352140 at TLR9 and toxoplasmic retinochoroiditis in Brazil. This supports the hypothesis that direct interaction between T. gondii and TLR9 may trigger proinflammatory responses that lead to severe pathologies such as the ocular disease that is associated with this infection in Brazil.

IL6 -174 G/C promoter polymorphism influences susceptibility to mucosal but not localized cutaneous leishmaniasis in Brazil.

BACKGROUND: Mucosal leishmaniasis (ML) is associated with exaggerated tumor necrosis factor- alpha and interferon- gamma responses and tissue destruction. ML follows localized cutaneous leishmaniasis (CL) caused by Leishmania braziliensis infection. Interleukin (IL)-6 down-regulates T helper (Th) cell type 1 differentiation and drives Th2 cell differentiation. The IL6 -174 G/C polymorphism is associated with proinflammatory diseases and IL-6 regulation. METHODS: The -174 G/C polymorphism was genotyped in population samples and families with CL and ML from Brazil. Genotype frequencies were compared among patients with ML, patients with CL, and 2 control groups by logistic regression and family-based association test (FBAT) analysis. IL-6 levels were measured in macrophages. RESULTS: The C allele was more common in patients with ML than in patients with CL (odds ratio [OR], 2.55 [95% confidence interval {CI}, 1.32-4.91]; P=.005), than in patients who were leishmanin skin-test positive (OR, 2.23 [95% CI, 1.23-4.05]; P=.009), and than in neighborhood control subjects (OR, 2.47 [95% CI, 1.24-4.90]; P=.01). FBAT analysis confirmed an association between allele C and ML under both additive (z=4.295; P=.000017) and dominant (z=4.325; P=.000015) models. Significantly lower levels of IL-6 were measured in unstimulated macrophages from CC individuals than from GG individuals (P=.003) as well as after stimulation with soluble leishmania antigen (P=.009). CONCLUSIONS: IL-6 may regulate type 1 proinflammatory responses, putting individuals with low macrophage IL-6 levels at increased risk for ML.