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.

Zika Virus Changes Methylation of Genes Involved in Immune Response and Neural Development in Brazilian Babies Born With Congenital Microcephaly.

The recent increase in babies born with brain and eye malformations in Brazil is associated with Zika virus (ZIKV) infection in utero. ZIKV alters host DNA methylation in vitro. Using genome-wide DNA methylation profiling we compared 18 babies born with congenital ZIKV microcephaly with 20 controls. We found ZIKV-associated alteration of host methylation patterns, notably at RABGAP1L which is important in brain development, at viral host immunity genes MX1 and ISG15, and in an epigenetic module containing the causal microcephaly gene MCPH1. Our data support the hypothesis that clinical signs of congenital ZIKV are associated with changes in DNA methylation.

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.

A flexible computational pipeline for research analyses of unsolved clinical exome cases.

Exome sequencing has enabled molecular diagnoses for rare disease patients but often with initial diagnostic rates of ~25-30%. Here we develop a robust computational pipeline to rank variants for reassessment of unsolved rare disease patients. A comprehensive web-based patient report is generated in which all deleterious variants can be filtered by gene, variant characteristics, OMIM disease and Phenolyzer scores, and all are annotated with an ACMG classification and links to ClinVar. The pipeline ranked 21/34 previously diagnosed variants as top, with 26 in total ranked ≤7th, 3 ranked ≥13th; 5 failed the pipeline filters. Pathogenic/likely pathogenic variants by ACMG criteria were identified for 22/145 unsolved cases, and a previously undefined candidate disease variant for 27/145. This open access pipeline supports the partnership between clinical and research laboratories to improve the diagnosis of unsolved exomes. It provides a flexible framework for iterative developments to further improve diagnosis.

Reference exome data for a Northern Brazilian population.

Exome sequencing is widely used in the diagnosis of rare genetic diseases and provides useful variant data for analysis of complex diseases. There is not always adequate population-specific reference data to assist in assigning a diagnostic variant to a specific clinical condition. Here we provide a catalogue of variants called after sequencing the exomes of 45 babies from Rio Grande do Nord in Brazil. Sequence data were processed using an 'intersect-then-combine' (ITC) approach, using GATK and SAMtools to call variants. A total of 612,761 variants were identified in at least one individual in this Brazilian Cohort, including 559,448 single nucleotide variants (SNVs) and 53,313 insertion/deletions. Of these, 58,111 overlapped with nonsynonymous (nsSNVs) or splice site (ssSNVs) SNVs in dbNSFP. As an aid to clinical diagnosis of rare diseases, we used the American College of Medicine Genetics and Genomics (ACMG) guidelines to assign pathogenic/likely pathogenic status to 185 (0.32%) of the 58,111 nsSNVs and ssSNVs. Our data set provides a useful reference point for diagnosis of rare diseases in Brazil. (169 words).

Reference exome data for Australian Aboriginal populations to support health-based research.

Whole exome sequencing (WES) is a popular and successful technology which is widely used in both research and clinical settings. However, there is a paucity of reference data for Aboriginal Australians to underpin the translation of health-based genomic research. Here we provide a catalogue of variants called after sequencing the exomes of 50 Aboriginal individuals from the Northern Territory (NT) of Australia and compare these to 72 previously published exomes from a Western Australian (WA) population of Martu origin. Sequence data for both NT and WA samples were processed using an 'intersect-then-combine' (ITC) approach, using GATK and SAMtools to call variants. A total of 289,829 variants were identified in at least one individual in the NT cohort and 248,374 variants in at least one individual in the WA cohort. Of these, 166,719 variants were present in both cohorts, whilst 123,110 variants were private to the NT cohort and 81,655 were private to the WA cohort. Our data set provides a useful reference point for genomic studies on Aboriginal Australians.

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.

Personalised analytics for rare disease diagnostics.

Whole genome and exome sequencing is a standard tool for the diagnosis of patients suffering from rare and other genetic disorders. The interpretation of the tens of thousands of variants returned from such tests remains a major challenge. Here we focus on the problem of prioritising variants with respect to the observed disease phenotype. We hypothesise that linking patterns of gene expression across multiple tissues to the phenotypes will aid in discovering disease causing variants. To test this, we construct classifiers that learn associations between tissue-specific gene expression and disease phenotypes. We find that using Genotype-Tissue Expression project (GTEx) expression data in conjunction with disease agnostic variant prioritisation methods (CADD or MetaSVM) results in consistent improvements in classification accuracy. Our method represents a previously overlooked avenue of utilising existing expression data for clinical diagnostics, and also opens the door to use of other functional genomic data sets in the same manner.

Meta-taxonomic analysis of prokaryotic and eukaryotic gut flora in stool samples from visceral leishmaniasis cases and endemic controls in Bihar State India.

Visceral leishmaniasis (VL) caused by Leishmania donovani remains of public health concern in rural India. Those at risk of VL are also at risk of other neglected tropical diseases (NTDs) including soil transmitted helminths. Intestinal helminths are potent regulators of host immune responses sometimes mediated through cross-talk with gut microbiota. We evaluate a meta-taxonomic approach to determine the composition of prokaryotic and eukaryotic gut microflora using amplicon-based sequencing of 16S ribosomal RNA (16S rRNA) and 18S rRNA gene regions. The most abundant bacterial taxa identified in faecal samples from Bihar State India were Prevotella (37.1%), Faecalibacterium (11.3%), Escherichia-Shigella (9.1%), Alloprevotella (4.5%), Bacteroides (4.1%), Ruminococcaceae UCG-002 (1.6%), and Bifidobacterium (1.5%). Eukaryotic taxa identified (excluding plant genera) included Blastocystis (57.9%; Order: Stramenopiles), Dientamoeba (12.1%; Family: Tritrichomonadea), Pentatrichomonas (10.1%; Family: Trichomonodea), Entamoeba (3.5%; Family: Entamoebida), Ascaridida (0.8%; Family: Chromodorea; concordant with Ascaris by microscopy), Rhabditida (0.8%; Family: Chromodorea; concordant with Strongyloides), and Cyclophyllidea (0.2%; Order: Eucestoda; concordant with Hymenolepis). Overall alpha (Shannon's, Faith's and Pielou's indices) and beta (Bray-Curtis dissimilarity statistic; weighted UniFrac distances) diversity of taxa did not differ significantly by age, sex, geographic subdistrict, or VL case (N = 23) versus endemic control (EC; N = 23) status. However, taxon-specific associations occurred: (i) Ruminococcaceae UCG- 014 and Gastranaerophilales_uncultured bacterium were enriched in EC compared to VL cases; (ii) Pentatrichomonas was more abundant in VL cases than in EC, whereas the reverse occurred for Entamoeba. Across the cohort, high Escherichia-Shigella was associated with reduced bacterial diversity, while high Blastocystis was associated with high bacterial diversity and low Escherichia-Shigella. Individuals with high Blastocystis had low Bacteroidaceae and high Clostridiales vadin BB60 whereas the reverse held true for low Blastocystis. This scoping study provides useful baseline data upon which to develop a broader analysis of pathogenic enteric microflora and their influence on gut microbial health and NTDs generally.

Transketolase and vitamin B1 influence on ROS-dependent neutrophil extracellular traps (NETs) formation.

Neutrophil extracellular traps (NETs) are a recently identified, web-like, extracellular structure composed of decondensed nuclear DNA and associated antimicrobial granules. NETs are extruded into the extracellular environment via the reactive oxygen species (ROS)-dependent cell death pathway participating in inflammation and autoimmune diseases. Transketolase (TKT) is a thiamine pyrophosphate (vitamin B1)-dependent enzyme that links the pentose phosphate pathway with the glycolytic pathway by feeding excess sugar phosphates into the main carbohydrate metabolic pathways to generate biosynthetic reducing capacity in the form of NADPH as a substrate for ROS generation. In this work, TKT was selected as a lead candidate from 24 NET-associated proteins obtained by literature screening and knowledge gap assessment. Consequently, we determined whether TKT influenced NET formation in vitro. We firstly established that the release of ROS-dependent NETs was significantly decreased after purified human PMNs were pretreated with oxythiamine, a TKT inhibitor, and in a concentration dependent manner. As a cofactor for TKT reaction, we evaluated the release of NET formation either in vitamin B1 treatment or in combined use of oxythiamine and vitamin B1, and found that those treatments also exerted a significant suppressive effect on the amount of NET-DNA and ROS production. The regulation of TKT by oxythiamine and/or vitamin B1 may therefore be associated with response to the modulation of NET formation by preventing generation of excessive NETs in inflammatory diseases.

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).

Publisher Correction: Toxoplasma Modulates Signature Pathways of Human Epilepsy, Neurodegeneration & Cancer.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper.

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.

Epigenetic dysregulation of host gene expression in Toxoplasma infection with specific reference to dopamine and amyloid pathways.

Recent interest has focussed on the influence of infectious disease organisms on the host epigenome. Toxoplasma gondii infection acquired congenitally or in early life is associated with severe ocular and brain developmental anomalies, while persistent asymptomatic infection is a proposed risk factor for neurodegenerative and psychiatric disorders, including Parkinson's and Alzheimer's Diseases, and schizophrenia. Genome-wide analysis of the host methylome and transcriptome following T. gondii infection in a retinal cell line identified genes (132, 186 and 128 genes at 2, 6 and 24 h post-infection) concordant for methylation and expression, i.e. hypermethylated and decreased expression or hypomethylated and increased expression. Pathway analyses showed perturbation of two neurologically-associated pathways: dopamine-DARPP32 feedback in cAMP signalling (p-value = 8.3 × 10-5; adjusted p-value = 0.020); and amyloid processing (p-value = 1.0 × 10-3; adjusted p-value = 0.043). Amyloid Precursor Protein (APP) decreased in level following T. gondii infection. These results are of interest given the expression of APP early in nervous system development affecting neural migration and the role of amyloid processing in Alzheimer's disease, while dopamine has roles in the developing retina as well as in Parkinson's disease and schizophrenia. Our results provide a possible functional link between T. gondii infection and congenital/early life and adult neurological clinical signs.

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.

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.

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.