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.

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.

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.

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.

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.

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.

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.

Arylsulphatase A Pseudodeficiency (ARSA-PD), hypertension and chronic renal disease in Aboriginal Australians.

Chronic renal disease (CRD) associated with cardiovascular disease (CVD) and/or type 2 diabetes (T2D) is a significant problem in Aboriginal Australians. Whole exome sequencing data (N = 72) showed enrichment for ClinVar pathogenic variants in gene sets/pathways linking lipoprotein, lipid and glucose metabolism. The top Ingenuity Pathway Analysis canonical pathways were Farsenoid X Receptor and Retinoid Receptor (FXR/RXR; (P = 1.86 × 10-7), Liver X Receptor and Retinoid Receptor (LXR/RXR; P = 2.88 × 10-6), and atherosclerosis signalling (P = 3.80 × 10-6). Top pathways/processes identified using Enrichr included: Reactome 2016 chylomicron-mediated lipid transport (P = 3.55 × 10-7); Wiki 2016 statin (P = 8.29 × 10-8); GO Biological Processes 2017 chylomicron remodelling (P = 1.92 × 10-8). ClinVar arylsulfatase A pseudodeficiency (ARSA-PD) pathogenic variants were common, including the missense variant c.511 G > A (p.Asp171Asn; rs74315466; frequency 0.44) only reported in Polynesians. This variant is in cis with known ARSA-PD 3' regulatory c.*96 A > G (rs6151429; frequency 0.47) and missense c.1055 A > G (p.Asn352Ser; rs2071421; frequency 0.47) variants. These latter two variants are associated with T2D (risk haplotype GG; odds ratio 2.67; 95% CI 2.32-3.08; P = 2.43 × 10-4) in genome-wide association data (N = 402), but are more strongly associated with quantitative traits (DBP, SBP, ACR, eGFR) for hypertension and renal function in non-diabetic than diabetic subgroups. Traits associated with CVD, CRD and T2D in Aboriginal Australians provide novel insight into function of ARSA-PD variants.

Functional Alleles of Chicken BG Genes, Members of the Butyrophilin Gene Family, in Peripheral T Cells.

γδ T cells recognize a wide variety of ligands in mammals, among them members of the butyrophilin (BTN) family. Nothing is known about γδ T cell ligands in chickens, despite there being many such cells in blood and lymphoid tissues, as well as in mucosal surfaces. The major histocompatibility complex (MHC) of chickens was discovered because of polymorphic BG genes, part of the BTN family. All but two BG genes are located in the BG region, oriented head-to-tail so that unequal crossing-over has led to copy number variation (CNV) as well as hybrid (chimeric) genes, making it difficult to identify true alleles. One approach is to examine BG genes expressed in particular cell types, which likely have the same functions in different BG haplotypes and thus can be considered "functional alleles." We cloned nearly full-length BG transcripts from peripheral T cells of four haplotypes (B2, B15, B19, and B21), and compared them to the BG genes of the B12 haplotype that previously were studied in detail. A dominant BG gene was found in each haplotype, but with significant levels of subdominant transcripts in three haplotypes (B2, B15, and B19). For three haplotypes (B15, B19, and B21), most sequences are closely-related to BG8, BG9, and BG12 from the B12 haplotype. We found that variation in the extracellular immunoglobulin-variable-like (Ig-V) domain is mostly localized to the membrane distal loops but without evidence for selection. However, variation in the cytoplasmic tail composed of many amino acid heptad repeats does appear to be selected (although not obviously localized), consistent with an intriguing clustering of charged and polar residues in an apparent α-helical coiled-coil. By contrast, the dominantly-expressed BG gene in the B2 haplotype is identical to BG13 from the B12 haplotype, and most of the subdominant sequences are from the BG5-BG7-BG11 clade. Moreover, alternative splicing leading to intron read-through results in dramatically truncated cytoplasmic tails, particularly for the dominantly-expressed BG gene of the B2 haplotype. The approach of examining "functional alleles" has yielded interesting data for closely-related genes, but also thrown up unexpected findings for at least one haplotype.

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.

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.

Fine mapping under linkage peaks for symptomatic or asymptomatic outcomes of Leishmania infantum infection in Brazil.

Infection with the protozoan Leishmania infantum can lead to asymptomatic infection and protective immunity, or to the progressive and potentially fatal disease visceral leishmaniasis (VL). Published studies show host genetic background determines in part whether infected individuals will develop a symptomatic or asymptomatic outcome. The purpose of the current study was to fine map chromosome regions previously linked with risk for symptomatic (chromosome 9) or asymptomatic (chromosomes 15 and 19) manifestations of L. infantum infection. We conducted a family-based genetic study of VL and asymptomatic infection (detected by a DTH skin test) with a final post quality control sample of 961 individuals with full genotype and phenotype information from highly endemic neighborhoods of northeast Brazil. A total of 5485 SNPs under the linkage peaks on chromosomes 9, 15 and 19 were genotyped. No strong SNP associations were observed for the DTH phenotype. The most significant associations with the VL phenotype were with SNP rs1470217 (p=5.9e-05; pcorrected=0.057) on chromosome 9, and with SNP rs8107014 (p=1.4e-05; pcorrected=0.013) on chromosome 19. SNP rs1470217 is situated in a 180kb intergenic region between TMEM215 (Transmembrane protein 215) and APTX (Aprataxin). SNP rs8107014 lies in the intron between exons 26 and 27 of a 34 exon transcript (ENST00000204005) of LTBP4, (Latent transforming growth factor-beta-binding protein 4a). The latter supports growing evidence that the transforming growth factor-beta pathway is important in the immunopathogenesis of VL.

First genome-wide association study in an Australian aboriginal population provides insights into genetic risk factors for body mass index and type 2 diabetes.

A body mass index (BMI) >22kg/m2 is a risk factor for type 2 diabetes (T2D) in Aboriginal Australians. To identify loci associated with BMI and T2D we undertook a genome-wide association study using 1,075,436 quality-controlled single nucleotide polymorphisms (SNPs) genotyped (Illumina 2.5M Duo Beadchip) in 402 individuals in extended pedigrees from a Western Australian Aboriginal community. Imputation using the thousand genomes (1000G) reference panel extended the analysis to 6,724,284 post quality-control autosomal SNPs. No associations achieved genome-wide significance, commonly accepted as P<5x10-8. Nevertheless, genes/pathways in common with other ethnicities were identified despite the arrival of Aboriginal people in Australia >45,000 years ago. The top hit (rs10868204 Pgenotyped = 1.50x10-6; rs11140653 Pimputed_1000G = 2.90x10-7) for BMI lies 5' of NTRK2, the type 2 neurotrophic tyrosine kinase receptor for brain-derived neurotrophic factor (BDNF) that regulates energy balance downstream of melanocortin-4 receptor (MC4R). PIK3C2G (rs12816270 Pgenotyped = 8.06x10-6; rs10841048 Pimputed_1000G = 6.28x10-7) was associated with BMI, but not with T2D as reported elsewhere. BMI also associated with CNTNAP2 (rs6960319 Pgenotyped = 4.65x10-5; rs13225016 Pimputed_1000G = 6.57x10-5), previously identified as the strongest gene-by-environment interaction for BMI in African-Americans. The top hit (rs11240074 Pgenotyped = 5.59x10-6, Pimputed_1000G = 5.73x10-6) for T2D lies 5' of BCL9 that, along with TCF7L2, promotes beta-catenin's transcriptional activity in the WNT signaling pathway. Additional hits occurred in genes affecting pancreatic (KCNJ6, KCNA1) and/or GABA (GABRR1, KCNA1) functions. Notable associations observed for genes previously identified at genome-wide significance in other populations included MC4R (Pgenotyped = 4.49x10-4) for BMI and IGF2BP2 Pimputed_1000G = 2.55x10-6) for T2D. Our results may provide novel functional leads in understanding disease pathogenesis in this Australian Aboriginal population.

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.

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.

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.

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.

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.