Anti-Interleukin-10 Unleashes Transcriptional Response to Leishmanial Antigens in Visceral Leishmaniasis Patients.

Visceral leishmaniasis (VL; Leishmania donovani) cases produce interferon-γ and tumor necrosis factor in response to soluble leishmanial antigen (SLA) in whole-blood assays. Using transcriptional profiling, we demonstrate the impact of interleukin-10 (IL-10), a cytokine implicated in VL, on this response. SLA stimulation identified 28 differentially expressed genes (DEGs), 17/28 in a single network with TNF as hub. SLA plus anti-IL-10 produced 454 DEGs, 292 in a single network with TNF, IFNG, NFKBIA, IL6, and IL1B as hubs in concert with a remarkable chemokine/cytokine storm. Our data demonstrate the singular effect of IL-10 as a potent immune modulator in VL.

Common and Rare Genetic Variants That Could Contribute to Severe Otitis Media in an Australian Aboriginal Population.

BACKGROUND: Our goal was to identify genetic risk factors for severe otitis media (OM) in Aboriginal Australians. METHODS: Illumina® Omni2.5 BeadChip and imputed data were compared between 21 children with severe OM (multiple episodes chronic suppurative OM and/or perforations or tympanic sclerosis) and 370 individuals without this phenotype, followed by FUnctional Mapping and Annotation (FUMA). Exome data filtered for common (EXaC_all ≥ 0.1) putative deleterious variants influencing protein coding (CADD-scaled scores ≥15] were used to compare 15 severe OM cases with 9 mild cases (single episode of acute OM recorded over ≥3 consecutive years). Rare (ExAC_all ≤ 0.01) such variants were filtered for those present only in severe OM. Enrichr was used to determine enrichment of genes contributing to pathways/processes relevant to OM. RESULTS: FUMA analysis identified 2 plausible genetic risk loci for severe OM: NR3C1 (Pimputed_1000G = 3.62 × 10-6) encoding the glucocorticoid receptor, and NREP (Pimputed_1000G = 3.67 × 10-6) encoding neuronal regeneration-related protein. Exome analysis showed: (i) association of severe OM with variants influencing protein coding (CADD-scaled ≥ 15) in a gene-set (GRXCR1, CDH23, LRP2, FAT4, ARSA, EYA4) enriched for Mammalian Phenotype Level 4 abnormal hair cell stereociliary bundle morphology and related phenotypes; (ii) rare variants influencing protein coding only seen in severe OM provided gene-sets enriched for "abnormal ear" (LMNA, CDH23, LRP2, MYO7A, FGFR1), integrin interactions, transforming growth factor signaling, and cell projection phenotypes including hair cell stereociliary bundles and cilium assembly. CONCLUSIONS: This study highlights interacting genes and pathways related to cilium structure and function that may contribute to extreme susceptibility to OM in Aboriginal Australian children.

A Genome-wide Association Study Identifies SERPINB10, CRLF3, STX7, LAMP3, IFNG-AS1, and KRT80 As Risk Loci Contributing to Cutaneous Leishmaniasis in Brazil.

BACKGROUND: Our goal was to identify genetic risk factors for cutaneous leishmaniasis (CL) caused by Leishmania braziliensis. METHODS: Genotyping 2066 CL cases and 2046 controls using Illumina HumanCoreExomeBeadChips provided data for 4 498 586 imputed single-nucleotide variants (SNVs). A genome-wide association study (GWAS) using linear mixed models took account of genetic diversity/ethnicity/admixture. Post-GWAS positional, expression quantitative trait locus (eQTL) and chromatin interaction mapping was performed in Functional Mapping and Annotation (FUMA). Transcriptional data were compared between lesions and normal skin, and cytokines measured using flow cytometry and Bioplex assay. RESULTS: Positional mapping identified 32 genomic loci associated with CL, none achieving genome-wide significance (P < 5 × 10-8). Lead SNVs at 23 loci occurred at protein coding or noncoding RNA genes, 15 with eQTLs for functionally relevant cells/tissues and/or showing differential expression in lesions. Of these, the 6 most plausible genetic risk loci were SERPINB10 (Pimputed_1000G = 2.67 × 10-6), CRLF3 (Pimputed_1000G = 5.12 × 10-6), STX7 (Pimputed_1000G = 6.06 × 10-6), KRT80 (Pimputed_1000G = 6.58 × 10-6), LAMP3 (Pimputed_1000G = 6.54 × 10-6), and IFNG-AS1 (Pimputed_1000G = 1.32 × 10-5). LAMP3 (Padjusted = 9.25 × 10-12; +6-fold), STX7 (Padjusted = 7.62 × 10-3; +1.3-fold), and CRLF3 (Padjusted = 9.19 × 10-9; +1.97-fold) were expressed more highly in CL biopsies compared to normal skin; KRT80 (Padjusted = 3.07 × 10-8; -3-fold) was lower. Multiple cis-eQTLs across SERPINB10 mapped to chromatin interaction regions of transcriptional/enhancer activity in neutrophils, monocytes, B cells, and hematopoietic stem cells. Those at IFNG-AS1 mapped to transcriptional/enhancer regions in T, natural killer, and B cells. The percentage of peripheral blood CD3+ T cells making antigen-specific interferon-γ differed significantly by IFNG-AS1 genotype. CONCLUSIONS: This first GWAS for CL identified multiple genetic risk loci including a novel lead to understanding CL pathogenesis through regulation of interferon-γ by IFNG antisense RNA 1.

Human genetics of leishmania infections.

Identifying genetic risk factors for parasitic infections such as the leishmaniases could provide important leads for improved therapies and vaccines. Until recently most genetic studies of human leishmaniasis were underpowered and/or not replicated. Here, we focus on recent genome-wide association studies of visceral leishmaniasis (VL) and cutaneous leishmaniasis (CL). For VL, analysis across 2287 cases and 2692 controls from three cohorts identified a single major peak of genome-wide significance (Pcombined = 2.76 × 10-17) at HLA-DRB1-HLA-DQA1. HLA-DRB1*1501 and DRB1*1404/DRB1*1301 were the most significant protective versus risk alleles, respectively, with specific residues at amino acid positions 11 and 13 unique to protective alleles. Epitope-binding studies showed higher frequency of basic AAs in DRB1*1404-/*1301-specific epitopes compared to hydrophobic and polar AAs in DRB1*1501-specific epitopes at anchor residues P4 and P6 which interact with residues at DRB1 positions 11 and 13. For CL, genome-wide significance was not achieved in combined analysis of 2066 cases and 2046 controls across 2 cohorts. Rather, multiple top hits at P < 5 × 10-5 were observed, amongst which IFNG-AS1 was of specific interest as a non-coding anti-sense RNA known to influence responses to pathogens by increasing IFN-γ secretion. Association at LAMP3 encoding dendritic cell lysosomal associated membrane protein 3 was also interesting. LAMP3 increases markedly upon activation of dendritic cells, localizing to the MHC Class II compartment immediately prior to translocation of Class II to the cell surface. Together these GWAS results provide firm confirmation for the importance of antigen presentation and the regulation of IFNγ in determining the outcome of Leishmania infections.

Epitope-Binding Characteristics for Risk versus Protective DRB1 Alleles for Visceral Leishmaniasis.

HLA-DRB1 is the major genetic risk factor for visceral leishmaniasis (VL). We used SNP2HLA to impute HLA-DRB1 alleles and SNPTEST to carry out association analyses in 889 human cases and 977 controls from India. NetMHCIIpan 2.1 was used to map epitopes and binding affinities across 49 Leishmania vaccine candidates, as well as across peptide epitopes captured from dendritic cells treated with crude Leishmania Ag and identified using mass spectrometry and alignment to amino acid sequences of a reference Leishmania genome. Cytokines were measured in peptide-stimulated whole blood from 26 cured VL cases and eight endemic healthy controls. HLA-DRB1*1501 and DRB1*1404/DRB1*1301 were the most significant protective and risk alleles, respectively, with specific residues at aa positions 11 and 13 unique to protective alleles. We observed greater peptide promiscuity in sequence motifs for 9-mer core epitopes predicted to bind to risk (*1404/*1301) compared with protective (*1501) DRB1 alleles. There was a higher frequency of basic amino acids in DRB1*1404/*1301-specific epitopes compared with hydrophobic and polar amino acids in DRB1*1501-specific epitopes at anchor residues pocket 4 and pocket 6, which interact with residues at DRB1 positions 11 and 13. Cured VL patients made variable, but robust, IFN-γ, TNF, and IL-10 responses to 20-mer peptides based on captured epitopes, with peptides based on DRB1*1501-captured epitopes resulting in a higher proportion (odds ratio 2.23, 95% confidence interval 1.17-4.25, p = 0.017) of patients with IFN-γ/IL-10 ratios > 2-fold compared with peptides based on DRB1*1301-captured epitopes. Our data provide insight into the molecular mechanisms underpinning the association of HLA-DRB1 alleles with risk versus protection in VL in humans.

HLA-DR Class II expression on myeloid and lymphoid cells in relation to HLA-DRB1 as a genetic risk factor for visceral leishmaniasis.

Genetic variation at HLA-DRB1 is a risk factor for visceral leishmaniasis (VL) caused by Leishmania donovani. The single nucleotide polymorphism rs9271252 upstream of the DRB1 gene provides a perfect tag for protective versus risk HLA-DRB1 four-digit alleles. In addition to the traditional role of the membrane-distal region of HLA class II molecules in antigen presentation and CD4 T-cell activation, the membrane-proximal region mediates 'non-traditional' multi-functional activation, differentiation, or death signals, including in DR-expressing T cells. To understand how HLA-DR contributes to disease pathogenesis, we examined expression at the protein level in circulating myeloid (CD14+ , CD16+ ) and lymphoid (CD4+ , CD8+ , CD19+ ) cells of VL patients (pre- and post-treatment) compared with endemic healthy controls (EHC). Although DR expression is reduced in circulating myeloid cells in active disease relative to EHC and post-treatment groups, expression is enhanced on CD4+  DR+ and CD8+  DR+ T cells consistent with T-cell activation. Cells of all myeloid and lymphoid populations from active cases were refractory to stimulation of DR expression with interferon-γ (IFN-γ). In contrast, all populations except CD19+ B cells from healthy blood bank controls showed enhanced DR expression following IFN-γ stimulation. The rs9271252 genotype did not impact significantly on IFN-γ-activated DR expression in myeloid, B or CD8+ T cells, but CD4+ T cells from healthy individuals homozygous for the risk allele were particularly refractory to activated DR expression. Further analysis of DR expression on subsets of CD4+ T cells regulating VL disease could uncover additional ways in which pleiotropy at HLA DRB1 contributes to disease pathogenesis.

Comprehensive candidate gene analysis for symptomatic or asymptomatic outcomes of Leishmania infantum infection in Brazil.

Genetic risk factors contribute to asymptomatic versus symptomatic visceral leishmaniasis (VL) outcomes following infection with Leishmania infantum. We therefore carried out a family-based (n = 918 post-quality control fully genotyped and phenotyped individuals) candidate gene study for symptomatic VL or asymptomatic delayed-type hypersensitivity (DTH) skin test phenotypes in highly endemic neighborhoods of northeast Brazil. A total of 248 SNPs were genotyped in 42 genes selected as candidates on the basis of prior genetic, immunological, and transcriptional profiling studies. The most significant association with the VL phenotype was with SNP rs6785358 (P = 5.7e-04; pcorrected = 0.026) 3.8 kb upstream of TGFBR2, the gene encoding the type 2 receptor for transforming growth factor beta (TGFß). A second inhibitory member of the TGBß superfamily signaling pathway, SMAD7, was associated with the DTH phenotype (SNP rs7238442: P = 0.001; pcorrected = 0.051). The most significant association for the DTH phenotype was with SNP rs10800309 (P = -8.4e-06; pcorrected = 3.9e-04) situated 3.1 kb upstream of FCGR2A, the gene encoding the low-affinity IIa receptor for the Fc fragment of IgG. Overall, our results imply a role for IgG-mediated inflammation in determining DTH associated with asymptomatic infection and contribute to growing evidence that the TGFß pathway is important in the immunopathogenesis of VL.

Comparison of methods to account for relatedness in genome-wide association studies with family-based data.

Approaches based on linear mixed models (LMMs) have recently gained popularity for modelling population substructure and relatedness in genome-wide association studies. In the last few years, a bewildering variety of different LMM methods/software packages have been developed, but it is not always clear how (or indeed whether) any newly-proposed method differs from previously-proposed implementations. Here we compare the performance of several LMM approaches (and software implementations, including EMMAX, GenABEL, FaST-LMM, Mendel, GEMMA and MMM) via their application to a genome-wide association study of visceral leishmaniasis in 348 Brazilian families comprising 3626 individuals (1972 genotyped). The implementations differ in precise details of methodology implemented and through various user-chosen options such as the method and number of SNPs used to estimate the kinship (relatedness) matrix. We investigate sensitivity to these choices and the success (or otherwise) of the approaches in controlling the overall genome-wide error-rate for both real and simulated phenotypes. We compare the LMM results to those obtained using traditional family-based association tests (based on transmission of alleles within pedigrees) and to alternative approaches implemented in the software packages MQLS, ROADTRIPS and MASTOR. We find strong concordance between the results from different LMM approaches, and all are successful in controlling the genome-wide error rate (except for some approaches when applied naively to longitudinal data with many repeated measures). We also find high correlation between LMMs and alternative approaches (apart from transmission-based approaches when applied to SNPs with small or non-existent effects). We conclude that LMM approaches perform well in comparison to competing approaches. Given their strong concordance, in most applications, the choice of precise LMM implementation cannot be based on power/type I error considerations but must instead be based on considerations such as speed and ease-of-use.

Insights into the possible role of IFNG and IFNGR1 in Kala-azar and Post Kala-azar Dermal Leishmaniasis in Sudanese patients.

BACKGROUND: Little is known about the parasite/host factors that lead to Post Kala-azar Dermal Leishmaniasis (PKDL) in some visceral leishmaniasis (VL) patients after drug-cure. Studies in Sudan provide evidence for association between polymorphisms in the gene (IFNGR1) encoding the alpha chain of interferon-γ receptor type I and risk of PKDL. This study aimed to identify putative functional polymorphisms in the IFNGR1 gene, and to determine whether differences in expression of interferon-γ (IFNG) and IFNGR1 at the RNA level are associated with pathogenesis of VL and/or PKDL in Sudan. METHODS: Sanger sequencing was used to re-sequence 841 bp of upstream, exon1 and intron1 of the IFNGR1 gene in DNA from 30 PKDL patients. LAGAN and SYNPLOT bioinformatics tools were used to compare human, chimpanzee and dog sequences to identify conserved noncoding sequences carrying putative regulatory elements. The relative expression of IFNG and IFNGR1 in paired pre- and post-treatment RNA samples from the lymph nodes of 24 VL patients, and in RNA samples from skin biopsies of 19 PKDL patients, was measured using real time PCR. Pre- versus post-treatment expression was evaluated statistically using the nonparametric Wilcoxon matched pairs signed-rank test. RESULTS: Ten variants were identified in the 841 bp of sequence, four of which are novel polymorphisms at -77A/G, +10 C/T, +18C/T and +91G/T relative to the IFNGR1 initiation site. A cluster of conserved non-coding sequences with putative regulatory variants was identified in the distal promoter of IFNGR1. Variable expression of IFNG was detected in lymph node aspirates of VL patients before treatment, with a marked reduction (P = 0.006) in expression following treatment. IFNGR1 expression was also variable in lymph node aspirates from VL patients, with no significant reduction in expression with treatment. IFNG expression was undetectable in the skin biopsies of PKDL cases, while IFNGR1 expression was also uniformly low. CONCLUSIONS: Uniformly low expression of IFN and IFNGR1 in PKDL skin biopsies could explain parasite persistence and is consistent with prior demonstration of genetic association with IFNGR1 polymorphisms. Identification of novel potentially functional rare variants at IFNGR1 makes an important general contribution to knowledge of rare variants of potential relevance in this Sudanese population.

Host genetic factors in American cutaneous leishmaniasis: a critical appraisal of studies conducted in an endemic area of Brazil.

American cutaneous leishmaniasis (ACL) is a vector-transmitted infectious disease with an estimated 1.5 million new cases per year. In Brazil, ACL represents a significant public health problem, with approximately 30,000 new reported cases annually, representing an incidence of 18.5 cases per 100,000 inhabitants. Corte de Pedra is in a region endemic for ACL in the state of Bahia (BA), northeastern Brazil, with 500-1,300 patients treated annually. Over the last decade, population and family-based candidate gene studies were conducted in Corte de Pedra, founded on previous knowledge from studies on mice and humans. Notwithstanding limitations related to sample size and power, these studies contribute important genetic biomarkers that identify novel pathways of disease pathogenesis and possible new therapeutic targets. The present paper is a narrative review about ACL immunogenetics in BA, highlighting in particular the interacting roles of the wound healing gene FLI1 with interleukin-6 and genes SMAD2 and SMAD3 of the transforming growth factor beta signalling pathway. This research highlights the need for well-powered genetic and functional studies on Leishmania braziliensis infection as essential to define and validate the role of host genes in determining resistance/susceptibility regarding this disease.

Association of a NOD2 gene polymorphism and T-helper 17 cells with presumed ocular toxoplasmosis.

Retinochoroiditis manifests in patients infected with Toxoplasma gondii. Here, we assessed 30 sibships and 89 parent/case trios of presumed ocular toxoplasmosis (POT) to evaluate associations with polymorphisms in the NOD2 gene. Three haplotype-tagging single-nucleotide polymorphisms (tag-SNPs) within the NOD2 gene were genotyped. The family-based association test showed that the tag-SNP rs3135499 is associated with retinochoroiditis (P = .039). We then characterized the cellular immune response of 59 cases of POT and 4 cases of active ocular toxoplasmosis (AOT). We found no differences in levels of interferon γ (IFN-γ) and interleukin 2 produced by T-helper 1 cells when comparing patients with AOT or POT to asymptomatic individuals. Unexpectedly, we found an increased interleukin 17A (IL-17A) production in patients with POT or OAT. In patients with POT or AOT, the main cellular source of IL-17A was CD4(+)CD45RO(+)T-bet(-)IFN-γ(-) T-helper 17 cells. Altogether, our results suggest that NOD2 influences the production of IL-17A by CD4(+) T lymphocytes and might contribute to the development of ocular toxoplasmosis.

A genome-wide search for type 2 diabetes susceptibility genes in an extended Arab family.

Twenty percent of people aged 20 to 79 have type 2 diabetes (T2D) in the United Arab Emirates (UAE). Genome-wide association studies (GWAS) to identify genes for T2D have not been reported for Arab countries. We performed a discovery GWAS in an extended UAE family (N=178; 66 diabetic; 112 healthy) genotyped on the Illumina Human 660 Quad Beadchip, with independent replication of top hits in 116 cases and 199 controls. Power to achieve genome-wide significance (commonly P=5×10(-8)) was therefore limited. Nevertheless, transmission disequilibrium testing in FBAT identified top hits at Chromosome 4p12-p13 (KCTD8: rs4407541, P=9.70×10(-6); GABRB1: rs10517178/rs1372491, P=4.19×10(-6)) and 14q13 (PRKD1: rs10144903, 3.92×10(-6)), supported by analysis using a linear mixed model approximation in GenABEL (4p12-p13 GABRG1/GABRA2: rs7662743, Padj-agesex=2.06×10(-5); KCTD8: rs4407541, Padj-agesex=1.42×10(-4); GABRB1: rs10517178/rs1372491, Padj-agesex=0.027; 14q13 PRKD1: rs10144903, Padj-agesex=6.95×10(-5)). SNPs across GABRG1/GABRA2 did not replicate, whereas more proximal SNPs rs7679715 (Padj-agesex=0.030) and rs2055942 (Padj-agesex=0.022) at COX7B2/GABRA4 did, in addition to a trend distally at KCTD8 (rs4695718: Padj-agesex=0.096). Modelling of discovery and replication data support independent signals at GABRA4 (rs2055942: Padj-agesex-combined=3×10(-4)) and at KCTD8 (rs4695718: Padj-agesex-combined=2×10(-4)). Replication was observed for PRKD1 rs1953722 (proxy for rs10144903; Padj-agesex=0.031; Padj-agesex-combined=2×10(-4)). These genes may provide important functional leads in understanding disease pathogenesis in this population.

Genetic and functional evidence implicating DLL1 as the gene that influences susceptibility to visceral leishmaniasis at chromosome 6q27.

BACKGROUND: Visceral leishmaniasis (VL) is caused by Leishmania donovani and Leishmania infantum chagasi. Genome-wide linkage studies from Sudan and Brazil identified a putative susceptibility locus on chromosome 6q27. METHODS: Twenty-two single-nucleotide polymorphisms (SNPs) at genes PHF10, C6orf70, DLL1, FAM120B, PSMB1, and TBP were genotyped in 193 VL cases from 85 Sudanese families, and 8 SNPs at genes PHF10, C6orf70, DLL1, PSMB1, and TBP were genotyped in 194 VL cases from 80 Brazilian families. Family-based association, haplotype, and linkage disequilibrium analyses were performed. Multispecies comparative sequence analysis was used to identify conserved noncoding sequences carrying putative regulatory elements. Quantitative reverse-transcription polymerase chain reaction measured expression of candidate genes in splenic aspirates from Indian patients with VL compared with that in the control spleen sample. RESULTS: Positive associations were observed at PHF10, C6orf70, DLL1, PSMB1, and TBP in Sudan, but only at DLL1 in Brazil (combined P = 3 × 10(-4) at DLL1 across Sudan and Brazil). No functional coding region variants were observed in resequencing of 22 Sudanese VL cases. DLL1 expression was significantly (P = 2 × 10(-7)) reduced (mean fold change, 3.5 [SEM, 0.7]) in splenic aspirates from patients with VL, whereas other 6q27 genes showed higher levels (1.27 × 10(-6) < P < .01) than did the control spleen sample. A cluster of conserved noncoding sequences with putative regulatory variants was identified in the distal promoter of DLL1. CONCLUSIONS: DLL1, which encodes Delta-like 1, the ligand for Notch3, is strongly implicated as the chromosome 6q27 VL susceptibility gene.

No evidence for association between SLC11A1 and visceral leishmaniasis in India.

BACKGROUND: SLC11A1 has pleiotropic effects on macrophage function and remains a strong candidate for infectious disease susceptibility. 5' and/or 3' polymorphisms have been associated with tuberculosis, leprosy, and visceral leishmaniasis (VL). Most studies undertaken to date were under-powered, and none has been replicated within a population. Association with tuberculosis has replicated variably across populations. Here we investigate SLC11A1 and VL in India. METHODS: Nine polymorphisms (rs34448891, rs7573065, rs2276631, rs3731865, rs17221959, rs2279015, rs17235409, rs17235416, rs17229009) that tag linkage disequilibrium blocks across SLC11A1 were genotyped in primary family-based (313 cases; 176 families) and replication (941 cases; 992 controls) samples. Family- and population-based analyses were performed to look for association between SLC11A1 variants and VL. Quantitative RT/PCR was used to compare SLC11A1 expression in mRNA from paired splenic aspirates taken before and after treatment from 24 VL patients carrying different genotypes at the functional promoter GTn polymorphism (rs34448891). RESULTS: No associations were observed between VL and polymorphisms at SLC11A1 that were either robust to correction for multiple testing or replicated across primary and replication samples. No differences in expression of SLC11A1 were observed when comparing pre- and post-treatment samples, or between individuals carrying different genotypes at the GTn repeat. CONCLUSIONS: This is the first well-powered study of SLC11A1 as a candidate for VL, which we conclude does not have a major role in regulating VL susceptibility in India.

Genetic and functional evaluation of the role of CXCR1 and CXCR2 in susceptibility to visceral leishmaniasis in north-east India.

BACKGROUND: IL8RA and IL8RB, encoded by CXCR1 and CXCR2, are receptors for interleukin (IL)-8 and other CXC chemokines involved in chemotaxis and activation of polymorphonuclear neutrophils (PMN). Variants at CXCR1 and CXCR2 have been associated with susceptibility to cutaneous and mucocutaneous leishmaniasis in Brazil. Here we investigate the role of CXCR1/CXCR2 in visceral leishmaniasis (VL) in India. METHODS: Three single nucleotide polymorphisms (SNPs) (rs4674259, rs2234671, rs3138060) that tag linkage disequilibrium blocks across CXCR1/CXCR2 were genotyped in primary family-based (313 cases; 176 nuclear families; 836 individuals) and replication (941 cases; 992 controls) samples. Family- and population-based analyses were performed to look for association between CXCR1/CXCR2 variants and VL. Quantitative RT/PCR was used to compare CXCR1/CXCR2 expression in mRNA from paired splenic aspirates taken before and after treatment from 19 VL patients. RESULTS: Family-based analysis using FBAT showed association between VL and SNPs CXCR1_rs2234671 (Z-score = 2.935, P = 0.003) and CXCR1_rs3138060 (Z-score = 2.22, P = 0.026), but not with CXCR2_rs4674259. Logistic regression analysis of the case-control data under an additive model of inheritance showed association between VL and SNPs CXCR2_rs4674259 (OR = 1.15, 95%CI = 1.01-1.31, P = 0.027) and CXCR1_rs3138060 (OR = 1.25, 95%CI = 1.02-1.53, P = 0.028), but not with CXCR1_rs2234671. The 3-locus haplotype T_G_C across these SNPs was shown to be the risk haplotype in both family- (TRANSMIT; P = 0.014) and population- (OR = 1.16, P = 0.028) samples (combined P = 0.002). CXCR2, but not CXCR1, expression was down regulated in pre-treatment compared to post-treatment splenic aspirates (P = 0.021). CONCLUSIONS: This well-powered primary and replication genetic study, together with functional analysis of gene expression, implicate CXCR2 in determining outcome of VL in India.

Epigenomic landscape of human colorectal cancer unveils an aberrant core of pan-cancer enhancers orchestrated by YAP/TAZ.

Cancer is characterized by pervasive epigenetic alterations with enhancer dysfunction orchestrating the aberrant cancer transcriptional programs and transcriptional dependencies. Here, we epigenetically characterize human colorectal cancer (CRC) using de novo chromatin state discovery on a library of different patient-derived organoids. By exploring this resource, we unveil a tumor-specific deregulated enhancerome that is cancer cell-intrinsic and independent of interpatient heterogeneity. We show that the transcriptional coactivators YAP/TAZ act as key regulators of the conserved CRC gained enhancers. The same YAP/TAZ-bound enhancers display active chromatin profiles across diverse human tumors, highlighting a pan-cancer epigenetic rewiring which at single-cell level distinguishes malignant from normal cell populations. YAP/TAZ inhibition in established tumor organoids causes extensive cell death unveiling their essential role in tumor maintenance. This work indicates a common layer of YAP/TAZ-fueled enhancer reprogramming that is key for the cancer cell state and can be exploited for the development of improved therapeutic avenues.

Y chromosome lineage- and village-specific genes on chromosomes 1p22 and 6q27 control visceral leishmaniasis in Sudan.

Familial clustering and ethnic differences suggest that visceral leishmaniasis caused by Leishmania donovani is under genetic control. A recent genome scan provided evidence for a major susceptibility gene on Chromosome 22q12 in the Aringa ethnic group in Sudan. We now report a genome-wide scan using 69 families with 173 affected relatives from two villages occupied by the related Masalit ethnic group. A primary ten-centimorgan scan followed by refined mapping provided evidence for major loci at 1p22 (LOD score 5.65; nominal p = 1.72 x 10(-7); empirical p < 1 x 10(-5); lambdaS = 5.1) and 6q27 (LOD score 3.74; nominal p = 1.68 x 10(-5); empirical p < 1 x 10(-4); lambdaS = 2.3) that were Y chromosome-lineage and village-specific. Neither village supported a visceral leishmaniasis susceptibility gene on 22q12. The results suggest strong lineage-specific genes due to founder effect and consanguinity in these recently immigrant populations. These chance events in ethnically uniform African populations provide a powerful resource in the search for genes and mechanisms that regulate this complex disease.

Genetics, Transcriptomics and Meta-Taxonomics in Visceral Leishmaniasis.

Visceral leishmaniasis (VL) caused by parasites of the Leishmania donovani complex can be fatal in susceptible individuals. Understanding the interactions between host and pathogen is one way to obtain leads to develop better drugs and for vaccine development. In recent years multiple omics-based approaches have assisted researchers to gain a more global picture of this interaction in leishmaniasis. Here we review results from studies using three omics-based approaches to study VL caused by L. donovani in India: (i) chip-based analysis of single nucleotide variants in the first genome-wide association study of host genetic risk factors for VL, followed by analysis of epitope binding to HLA DRB1 risk versus protective alleles; (ii) transcriptional profiling demonstrating pathways important in Amphotericin B treated compared to active VL cases, including demonstration that anti-interleukin-10 unleashes a storm of chemokines and cytokines in whole blood responses to soluble leishmania antigen in active cases; and (iii) a meta-taxonomic approach based on sequencing amplicons derived from regions of 16S ribosomal RNA (16S rRNA) and 18S rRNA genes that allowed us to determine composition of both prokaryotic and eukaryotic gut microflora in VL cases compared to endemic controls. Overall, our omics-based approaches demonstrate that global analyses of genetic risk factors, host responses to infection, and the interaction between host, parasite and the microbiome can point to the most critical factors that determine the outcome of infection.

Transcriptional blood signatures for active and amphotericin B treated visceral leishmaniasis in India.

Amphotericin B provides improved therapy for visceral leishmaniasis (VL) caused by Leishmania donovani, with single dose liposomal-encapsulated Ambisome providing the best cure rates. The VL elimination program aims to reduce the incidence rate in the Indian subcontinent to <1/10,000 population/year. Ability to predict which asymptomatic individuals (e.g. anti-leishmanial IgG and/or Leishmania-specific modified Quantiferon positive) will progress to clinical VL would help in monitoring disease outbreaks. Here we examined whole blood transcriptional profiles associated with asymptomatic infection, active disease, and in treated cases. Two independent microarray experiments were performed, with analysis focussed primarily on differentially expressed genes (DEGs) concordant across both experiments. No DEGs were identified for IgG or Quantiferon positive asymptomatic groups compared to negative healthy endemic controls. We therefore concentrated on comparing concordant DEGs from active cases with all healthy controls, and in examining differences in the transcriptome following different regimens of drug treatment. In these comparisons 6 major themes emerged: (i) expression of genes and enrichment of gene sets associated with erythrocyte function in active cases; (ii) strong evidence for enrichment of gene sets involved in cell cycle in comparing active cases with healthy controls; (iii) identification of IFNG encoding interferon-γ as the major hub gene in concordant gene expression patterns across experiments comparing active cases with healthy controls or with treated cases; (iv) enrichment for interleukin signalling (IL-1/3/4/6/7/8) and a prominent role for CXCL10/9/11 and chemokine signalling pathways in comparing active cases with treated cases; (v) the novel identification of Aryl Hydrocarbon Receptor signalling as a significant canonical pathway when comparing active cases with healthy controls or with treated cases; and (vi) global expression profiling support for more effective cure at day 30 post-treatment with a single dose of liposomal encapsulated amphotericin B compared to multi-dose non-liposomal amphotericin B treatment over 30 days. (296 words; 300 words allowed).

Determinants for progression from asymptomatic infection to symptomatic visceral leishmaniasis: A cohort study.

BACKGROUND: Asymptomatic Leishmania donovani infections outnumber clinical presentations, however the predictors for development of active disease are not well known. We aimed to identify serological, immunological and genetic markers for progression from L. donovani infection to clinical Visceral Leishmaniasis (VL). METHODS: We enrolled all residents >2 years of age in 27 VL endemic villages in Bihar (India). Blood samples collected on filter paper on two occasions 6-12 months apart, were tested for antibodies against L. donovani with rK39-ELISA and DAT. Sero converters, (negative for both tests in the first round but positive on either of the two during the second round) and controls (negative on both tests on both occasions) were followed for three years. At the start of follow-up venous blood was collected for the following tests: DAT, rK39- ELISA, Quantiferon assay, SNP/HLA genotyping and L.donovani specific quantitative PCR. RESULTS: Among 1,606 subjects enrolled,17 (8/476 seroconverters and 9/1,130 controls) developed VL (OR 3.1; 95% CI 1.1-8.3). High DAT and rK39 ELISA antibody titers as well as positive qPCR were strongly and significantly associated with progression from seroconversion to VL with odds ratios of 19.1, 30.3 and 20.9 respectively. Most VL cases arose early (median 5 months) during follow-up. CONCLUSION: We confirmed the strong association between high DAT and/or rK39 titers and progression to disease among asymptomatic subjects and identified qPCR as an additional predictor. Low predictive values do not warrant prophylactic treatment but as most progressed to VL early during follow-up, careful oberservation of these subjects for at least 6 months is indicated.