Search

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., Genome-wide analysis identifies common noncoding putative regulatory, as well as common and rare potentially deleterious variants directly influencing protein coding, in genes and pathways related to cilium structure and function that may contribute to severe otitis media in Aboriginal Australians.

Background Over 400 million people worldwide are living with a rare disease. Next Generation Sequencing (NGS) identifies potential disease causative genetic variants. However, many are identified as variants of uncertain significance (VUS) and require functional laboratory validation to determine pathogenicity, and this creates major diagnostic delays. Methods In this study we test a rapid genetic variant assessment pipeline using CRISPR homology directed repair to introduce single nucleotide variants into inducible pluripotent stem cells (iPSCs), followed by neuronal disease modelling, and functional genomics on amplicon and RNA sequencing, to determine cellular changes to support patient diagnosis and identify disease mechanism. Results As proof-of-principle, we investigated an EHMT1 (Euchromatin histone methyltransferase 1; EHMT1 c.3430C > T; p.Gln1144*) genetic variant pathogenic for Kleefstra syndrome and determined changes in gene expression during neuronal progenitor cell differentiation. This pipeline rapidly identified Kleefstra syndrome in genetic variant cells compared to healthy cells, and revealed novel findings potentially implicating the key transcription factors REST and SP1 in disease pathogenesis. Conclusion The study pipeline is a rapid, robust method for genetic variant assessment that will support rare diseases patient diagnosis. The results also provide valuable information on genome wide perturbations key to disease mechanism that can be targeted for drug treatments.

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., Measurement(s)Aboriginal Australian • DNA • sequence feature annotationTechnology Type(s)Whole Exome Sequencing • DNA sequencing • sequence annotationFactor Type(s)ancestry • sex • ageSample Characteristic - OrganismHomo sapiensSample Characteristic - LocationNorthern Territory Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.12040638

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.

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.

Background Rheumatic heart disease (RHD) after group A streptococcus (GAS) infections is heritable and prevalent in Indigenous populations. Molecular mimicry between human and GAS proteins triggers proinflammatory cardiac valve-reactive T cells. Methods Genome-wide genetic analysis was undertaken in 1263 Aboriginal Australians (398 RHD cases; 865 controls). Single-nucleotide polymorphisms were genotyped using Illumina HumanCoreExome BeadChips. Direct typing and imputation was used to fine-map the human leukocyte antigen (HLA) region. Epitope binding affinities were mapped for human cross-reactive GAS proteins, including M5 and M6. Results The strongest genetic association was intronic to HLA-DQA1 (rs9272622; P = 1.86 × 10-7). Conditional analyses showed rs9272622 and/or DQA1*AA16 account for the HLA signal. HLA-DQA1*0101_DQB1*0503 (odds ratio [OR], 1.44; 95% confidence interval [CI], 1.09-1.90; P = 9.56 × 10-3) and HLA-DQA1*0103_DQB1*0601 (OR, 1.27; 95% CI, 1.07-1.52; P = 7.15 × 10-3) were risk haplotypes; HLA_DQA1*0301-DQB1*0402 (OR 0.30, 95%CI 0.14-0.65, P = 2.36 × 10-3) was protective. Human myosin cross-reactive N-terminal and B repeat epitopes of GAS M5/M6 bind with higher affinity to DQA1/DQB1 alpha/beta dimers for the 2-risk haplotypes than the protective haplotype. Conclusions Variation at HLA_DQA1-DQB1 is the major genetic risk factor for RHD in Aboriginal Australians studied here. Cross-reactive epitopes bind with higher affinity to alpha/beta dimers formed by risk haplotypes, supporting molecular mimicry as the key mechanism of RHD pathogenesis.

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, Author summary Visceral leishmaniasis (VL), also known as kala-azar, is a potentially fatal disease caused by intracellular parasites of the Leishmania donovani complex. VL is a serious public health problem in rural India, causing high morbidity and mortality, as well as major costs to local and national health budgets. Amphotericin B provides improved therapy for VL with single dose liposomal-encapsulated Ambisome, now affordable through WHO-negotiated price reductions, providing the best cure rates. The VL elimination program aims to reduce the incidence rate in the Indian subcontinent to

Conventional dendritic cells (cDC) resident in the lymphoid organs of mice have been classically divided into CD8+ and CD8neg subsets. It is well-established that CD8+ dendritic cells (DCs) and their migratory counterparts in the periphery comprise the cross-presenting cDC1 subset. In contrast, CD8neg DCs are grouped together in the heterogeneous cDC2 subset. CD8neg DCs are relatively poor cross-presenters and drive more prominent CD4+ T cell responses against exogenous antigens. The discovery of the X-C motif chemokine receptor 1 (XCR1) as a specific marker of cross-presenting DCs, has led to the identification of a divergent subset of CD8+ DCs that lacks the ability to cross-present. Here, we report that these poorly characterized CD8+XCR1neg DCs have a gene expression profile that is consistent with both plasmacytoid DCs (pDCs) and cDC2. Our data demonstrate that CD8+XCR1neg DCs possess a unique pattern of endocytic receptors and a restricted toll-like receptor (TLR) profile that is particularly enriched for TLR5, giving them a unique position within the DC immunosurveillance network.

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.

Upon invasion of host cells, the ubiquitous pathogen Toxoplasma gondii manipulates several host processes, including re-organization of host organelles, to create a replicative niche. Host mitochondrial association to T. gondii parasitophorous vacuoles is rapid and has roles in modulating host immune responses. Here gene expression profiling of T. gondii infected cells reveals enrichment of genes involved in oxidative phosphorylation (OXPHOS) and mitochondrial dysfunction 6 h post-infection. We identified 11 hub genes (HIF-1α, CASP8, FN1, POU5F1, CD44, ISG15, HNRNPA1, MDM2, RPL35, VHL, and NUPR1) and 10 predicted upstream regulators, including 4 endogenous regulators RICTOR, KDM5A, RB1, and D-glucose. We characterized a number of mitochondrial parameters in T. gondii infected human foreskin fibroblast cells over a 36 h time-course. In addition to the usual rapid recruitment and apparent enlargement of mitochondria around the parasitophorous vacuole we observed fragmented host mitochondria in infected cells, not linked to cellular apoptosis, from 24 h post-infection. An increase in mitochondrial superoxide levels in T. gondii infected cells was observed that required active parasite invasion and peaked at 30 h post-infection. Measurement of OXPHOS proteins showed decreased expression of Complex IV in infected cells at 24 h post-infection, followed by decreased expression of Complexes I and II at 36 h post-infection. No change occurred in Complex V. No difference in host mitochondrial membrane potential between infected and mock-infected cells was observed at any time. Our results show perturbation of host mitochondrial function following T. gondii infection that likely impacts on pathogenesis of disease.

BackgroundRheumatic heart disease (RHD) following Group A Streptococcus (GAS) infections is heritable and prevalent in Indigenous populations. Molecular mimicry between human and GAS proteins triggers pro-inflammatory cardiac valve-reactive T-cells.MethodsGenome-wide genetic analysis was undertaken in 1263 Aboriginal Australians (398 RHD cases; 865 controls). Single nucleotide polymorphisms (SNPs) were genotyped using Illumina HumanCoreExome BeadChips. Direct typing and imputation was used to fine-map the human leukocyte antigen (HLA) region. Epitope binding affinities were mapped for human cross-reactive GAS proteins, including M5 and M6.ResultsThe strongest genetic association was intronic to HLA-DQA1 (rs9272622; P=1.86x10−7). Conditional analyses showed rs9272622 and/or DQA1*AA16 account for the HLA signal. HLA-DQA1*0101_DQB1*0503 (OR 1.44, 95%CI 1.09-1.90, P=9.56x10−3) and HLA-DQA1*0103_DQB1*0601 (OR 1.27, 95%CI 1.07-1.52, P=7.15x10−3) were risk haplotypes; HLA_DQA1*0301-DQB1*0402 (OR 0.30, 95%CI 0.14-0.65, P=2.36x10−3) was protective. Human myosin cross-reactive N-terminal and B repeat epitopes of GAS M5/M6 bind with higher affinity to DQA1/DQB1 alpha/beta dimers for the two risk haplotypes than the protective haplotype.ConclusionsVariation at HLA_DQA1-DQB1 is the major genetic risk factor for RHD in Aboriginal Australians studied here. Cross-reactive epitopes bind with higher affinity to alpha/beta dimers formed by risk haplotypes, supporting molecular mimicry as the key mechanism of RHD pathogenesis.

Visceral leishmaniasis (VL) in Sudan caused by Leishmania donovani is fatal in susceptible individuals if untreated. Treatment with sodium stibogluconate (SSG) leads to post-kala-azar dermal leishmaniasis (PKDL) in 58% of patients. Here, Affymetrix microarrays were used to identify genes differentially expressed in lymph nodes (N=9 paired samples) pre- and post-treatment with SSG. Using the Bioconductor package limma, 438 genes from 28 869 post-quality-control probe sets were differentially expressed (Pnominal ≤.02) post- vs pretreatment. Canonical pathway analysis using Ingenuity Pathway Analysis™ identified "role of nuclear factor of activated T-cell in regulation of immune response" (Pnominal =1.35×10-5 ; PBH-adjusted =4.79×10-3 ), "B-cell development" (Pnominal =2.04×10-4 ; PBH-adjusted =.024), "Fcγ receptor-mediated phagocytosis in macrophages and monocytes" (Pnominal =2.04×10-4 ; PBH-adjusted =.024) and "OX40 signalling" (Pnominal =2.82×10-4 ; PBH-adjusted =.025) as pathways differentially regulated post- vs pretreatment. Major network hub genes included TP53, FN1, MYC, BCL2, JUN, SYK, RUNX2, MMP1 and ACTA2. Top endogenous upstream regulators included IL-7 (P=2.28×10-6 ), TNF (P=4.26×10-6 ), Amyloid Precursor Protein (P=4.23×10-5 ) and SPI1/PI.1 (P=1.17×10-7 ). Top predicted chemical drug regulators included the flavonoid genistein (P=4.56×10-7 ) and the quinoline alkaloid camptothecin (P=5.14×10-5 ). These results contribute to our understanding of immunopathology associated with VL and response to SSG treatment. Further replication could identify novel therapeutic strategies that improve on SSG treatment and reduce the likelihood of progression to PKDL.

Genetic analyses, including genome-wide association studies and whole exome sequencing (WES), provide powerful tools for the analysis of complex and rare genetic diseases. To date there are no 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 72 Aboriginal individuals to a depth of 20X coverage in ∼80% of the sequenced nucleotides. We determined 320,976 single nucleotide variants (SNVs) and 47,313 insertions/deletions using the Genome Analysis Toolkit. We had previously genotyped a subset of the Aboriginal individuals (70/72) using the Illumina Omni2.5 BeadChip platform and found ~99% concordance at overlapping sites, which suggests high quality genotyping. Finally, we compared our SNVs to six publicly available variant databases, such as dbSNP and the Exome Sequencing Project, and 70,115 of our SNVs did not overlap any of the single nucleotide polymorphic sites in all the databases. Our data set provides a useful reference point for genomic studies on Aboriginal Australians.

Viruses, bacteria, and parasites have developed strategies to invade and establish long-term infections in their hosts. They have been shown to manipulate their hosts' processes such as the prevention of apoptosis and suppression of the immune response in order to successfully colonize their hosts. Recently, there is increasing interest in the study of epigenetics during an infection, particularly the changes induced in the hosts' DNA methylome, histone marks, and microRNA profiles. In this chapter, we summarize the literature published about viral-, bacterial-, and parasite-induced changes in the host epigenome. We also discuss some pathogens and the potential of the changes they induce in the host epigenome as epigenetic biomarkers to predict clinical outcomes of infection, such as pathogen-induced cancers, or to differentiate between active or latent infection. This might provide the basis to develop more efficient treatments or complement the diagnostic assays currently used clinically.