Search

Abstract Alzheimer’s disease (AD) is a complex genetic disease, and variants identified through genome-wide association studies (GWAS) explain only part of its heritability. Epistasis has been proposed as a major contributor to this ‘missing heritability’, however, many current methods are limited to only modelling additive effects. We use VariantSpark, a machine learning approach to GWAS, and BitEpi, a tool for epistasis detection, to identify AD associated variants and interactions across two independent cohorts, ADNI and UK Biobank. By incorporating significant epistatic interactions, we captured 10.41% more phenotypic variance than logistic regression (LR). We validate the well-established AD loci, APOE, and identify two novel genome-wide significant AD associated loci in both cohorts, SH3BP4 and SASH1, which are also in significant epistatic interactions with APOE. We show that the SH3BP4 SNP has a modulating effect on the known pathogenic APOE SNP, demonstrating a possible protective mechanism against AD. SASH1 is involved in a triplet interaction with pathogenic APOE SNP and ACOT11, where the SASH1 SNP lowered the pathogenic interaction effect between ACOT11 and APOE. Finally, we demonstrate that VariantSpark detects disease associations with 80% fewer controls than LR, unlocking discoveries in well annotated but smaller cohorts.

Random forests (RFs) are a widely used modelling tool capable of feature selection via a variable importance measure (VIM), however, a threshold is needed to control for false positives. In the absence of a good understanding of the characteristics of VIMs, many current approaches attempt to select features associated to the response by training multiple RFs to generate statistical power via a permutation null, by employing recursive feature elimination, or through a combination of both. However, for high-dimensional datasets these approaches become computationally infeasible. In this paper, we present RFlocalfdr, a statistical approach, built on the empirical Bayes argument of Efron, for thresholding mean decrease in impurity (MDI) importances. It identifies features significantly associated with the response while controlling the false positive rate. Using synthetic data and real-world data in health, we demonstrate that RFlocalfdr has equivalent accuracy to currently published approaches, while being orders of magnitude faster. We show that RFlocalfdr can successfully threshold a dataset of 106 datapoints, establishing its usability for large-scale datasets, like genomics. Furthermore, RFlocalfdr is compatible with any RF implementation that returns a VIM and counts, making it a versatile feature selection tool that reduces false discoveries.

Hyperactivation of the toll-like receptors (TLRs) have been implicated as risk factors for more severe forms of disease in COVID-19 and metabolic (dysfunction) associated fatty liver disease (MAFLD). Here the authors report that MBOAT7 is reduced in macrophages of patients with MAFLD and COVID-19, and acts as a negative regulator of TLR signalling.

New SARS-CoV-2 variants emerge as part of the virus’ adaptation to the human host. The Health Organizations are monitoring newly emerging variants with suspected impact on disease or vaccination efficacy as Variants Being Monitored (VBM), like Delta and Omicron. Genetic changes (SNVs) compared to the Wuhan variant characterize VBMs with current emphasis on the spike protein and lineage markers. However, monitoring VBMs in such a way might miss SNVs with functional effect on disease.Here we introduce a lineage-agnostic genome-wide approach to identify SNVs associated with disease. We curated a case-control dataset of 10,520 samples and identified 117 SNVs significantly associated with adverse patient outcome. While 40% (47) SNV are already monitored and 36% (43) are in the spike protein, we also identified 70 new SNVs that are associated with disease outcome. 31 of these are disease-worsening and predominantly located in the 3′-5′ exonuclease (NSP14) with structural modelling revealing a concise cluster in the Zn binding domain that has known host-immune modulating function. Furthermore, we generate clade-independent VBM groupings by identifying interacting SNVs (epistasis). We find 37 sets of higher-order epistatic interactions joining 5 genomic regions (nsp3, nsp14, Spike S1, ORF3a, N). Structural modelling of these regions provides insights into potential mechanistic pathways of increased virulence as well as orthogonal methods of validation.Clade-independent monitoring of functionally interacting (epistasis, co-evolution) SNVs detected emerging VBM a week before they were flagged by Health Organizations and in conjunction with structural modelling provides faster, mechanistic insight into emerging strains to guide public health interventions.

Clinical genomics is a system of multiple stakeholders and institutions. Yet, studies focusing on the comparative perspectives of these stakeholders are limited. This study engages four groups of professionals (clinical geneticists, genetic counselors, laboratory professionals, and researchers) working in clinical genomics to investigate their perceptions of the benefits and risks of using genomics in Australian healthcare. The study is underpinned by a risk governance approach. For data collection, qualitative semi-structured interviews were used. Our results show that all professionals unanimously identified that the benefit of clinical genomics lies in improving health outcomes for patients. However, the risks associated with delivering this benefit differed by professional category. We found that the further the profession was from the patient (e.g. researcher) the narrower the perceived risks were amongst the individuals interviewed. However, “privacy” as a perceived risk was ranked highly by all professions indicating a shared desire for responsible data governance practices.

Abstract Complex genetic diseases may be modulated by a large number of epistatic interactions affecting a polygenic phenotype. Identifying these interactions is difficult due to computational complexity, especially in the case of higher-order interactions where more than two genomic variants are involved. In this paper, we present BitEpi, a fast and accurate method to test all possible combinations of up to four bi-allelic variants (i.e. Single Nucleotide Variant or SNV for short). BitEpi introduces a novel bitwise algorithm that is 1.7 and 56 times faster for 3-SNV and 4-SNV search, than established software. The novel entropy statistic used in BitEpi is 44% more accurate to identify interactive SNVs, incorporating a p-value-based significance testing. We demonstrate BitEpi on real world data of 4900 samples and 87,000 SNPs. We also present EpiExplorer to visualize the potentially large number of individual and interacting SNVs in an interactive Cytoscape graph. EpiExplorer uses various visual elements to facilitate the discovery of true biological events in a complex polygenic environment.

The essential amino acid tryptophan (TRP) is the initiating metabolite of the kynurenine pathway (KP), which can be upregulated by inflammatory conditions in cells. Neuroinflammation-triggered activation of the KP and excessive production of the KP metabolite quinolinic acid are common features of multiple neurodegenerative diseases, including amyotrophic lateral sclerosis (ALS). In addition to its role in the KP, genes involved in TRP metabolism, including its incorporation into proteins, and synthesis of the neurotransmitter serotonin, have also been genetically and functionally linked to these diseases. ALS is a late onset neurodegenerative disease that is classified as familial or sporadic, depending on the presence or absence of a family history of the disease. Heritability estimates support a genetic basis for all ALS, including the sporadic form of the disease. However, the genetic basis of sporadic ALS (SALS) is complex, with the presence of multiple gene variants acting to increase disease susceptibility and is further complicated by interaction with potential environmental factors. We aimed to determine the genetic contribution of 18 genes involved in TRP metabolism, including protein synthesis, serotonin synthesis and the KP, by interrogating whole-genome sequencing data from 614 Australian sporadic ALS cases. Five genes in the KP (AFMID, CCBL1, GOT2, KYNU, HAAO) were found to have either novel protein-altering variants, and/or a burden of rare protein-altering variants in SALS cases compared to controls. Four genes involved in TRP metabolism for protein synthesis (WARS) and serotonin synthesis (TPH1, TPH2, MAOA) were also found to carry novel variants and/or gene burden. These variants may represent ALS risk factors that act to alter the KP and lead to neuroinflammation. These findings provide further evidence for the role of TRP metabolism, the KP and neuroinflammation in ALS disease pathobiology.

Human bone marrow-derived stromal (skeletal) stem cells (BM-hMSC) are being employed in an increasing number of clinical trials for tissue regeneration. A limiting factor for their clinical use is the inability to obtain sufficient cell numbers. Human embryonic stem cells (hESC) can provide an unlimited source of clinical grade cells for therapy. We have generated MSC-like cells from hESC (called here hESC-stromal) that exhibit surface markers and differentiate to osteoblasts and adipocytes, similar to BM-hMSC. In the present study, we used microarray analysis to compare the molecular phenotype of hESC-stromal and immortalised BM-hMSC cells (hMSC-TERT). Of the 7379 genes expressed above baseline, only 9.3% of genes were differentially expressed between undifferentiated hESC-stromal and BM-hMSC. Following ex vivo osteoblast induction, 665 and 695 genes exhibited ≥2-fold change (FC) in hESC-stromal and BM-hMSC, respectively with 172 genes common to both cell types. Functional annotation of significantly changing genes revealed similarities in gene ontology between the two cell types. Interestingly, genes in categories of cell adhesion/motility and epithelial–mesenchymal transition (EMT) were highly enriched in hESC-stromal whereas genes associated with cell cycle processes were enriched in hMSC-TERT. This data suggests that while hESC-stromal cells exhibit a similar molecular phenotype to hMSC-TERT, differences exist that can be explained by ontological differences between these two cell types. hESC-stromal cells can thus be considered as a possible alternative candidate cells for hMSC, to be employed in regenerative medicine protocols.

Background Aberrant or impaired repair of double-strand DNA breaks is a common feature of de novo acute myeloid leukemia and myelodysplastic syndromes. Since poly (ADP-ribose) polymerase (PARP) inhibitors have been recently shown to selectively target cells with defects in double-strand DNA repair, the aim of this study was to explore the possibility of exploiting defects in DNA repair in leukemic cells using PARP inhibitors.Design and Methods Leukemic cell lines were exposed to various PARP inhibitors alone and in combination with non-cytotoxic concentrations of DNA methyltransferase inhibitor, 5’ aza-2’-deoxycytidine and/or the histone deacetylase inhibitor, MS275, to test for potentiation of apoptosis with these agents.Results PARP inhibitors, KU-0058948 and PJ34, induced cell cycle arrest and apoptosis of primary myeloid leukemic cells and myeloid leukemic cell lines in vitro. Immunofluorescence analysis also revealed that PARP inhibitor sensitivity in these leukemic cells was due to a defect in homologous recombination DNA repair. Addition of 5’ aza-2’-deoxycytidine failed to increase the cytotoxicity of PARP inhibitors. In contrast, MS275 potentiated the cytotoxic effect of KU-0058948 and PJ34 in all PARP inhibitor-sensitive leukemic cells. Immunofluorescence analysis supported the idea that histone deacetylase inhibitors potentiate cytotoxicity by inhibiting DNA repair processes.Conclusions On the basis of the data presented here, we suggest that PARP inhibitors can potentially exploit defects in double-strand DNA break repair in leukemic cells, paving the way for testing the therapeutic potential of these agents in myelodysplastic syndromes and acute myeloid leukemia.

Abstract Background Chlorotoxin (Cltx) isolated from scorpion venom is an established tumor targeting and antiangiogenic peptide. Radiolabeled Cltx therapeutic (131I-TM601) yielded promising results in human glioma clinical studies, and the imaging agent tozuleristide, is under investigation in CNS cancer studies. Several binding targets have previously been proposed for Cltx but none effectively explain its pleiotropic effects; its true target remains ambiguous and is the focus of this study. Methods A peptide-drug conjugate (ER-472) composed of Cltx linked to cryptophycin as warhead was developed as a tool to probe the molecular target and mechanism of action of Cltx, using multiple xenograft models. Results Neuropilin-1 (NRP1), an endocytic receptor on tumor and endothelial cells, was identified as a novel Cltx target, and NRP1 binding by Cltx increased drug uptake into tumor. Metabolism of Cltx to peptide bearing free C-terminal arginine, a prerequisite for NRP1 binding, took place in the tumor microenvironment, while native scorpion Cltx with amidated C-terminal arginine did not bind NRP1, and instead acts as a cryptic peptide. Antitumor activity of ER-472 in xenografts correlated to tumor NRP1 expression. Potency was significantly reduced by treatment with NRP1 blocking antibodies or knockout in tumor cells, confirming a role for NRP1-binding in ER-472 activity. Higher cryptophycin metabolite levels were measured in NRP1-expressing tumors, evidence of NRP1-mediated enhanced drug uptake and presumably responsible for the superior antitumor efficacy. Conclusions NRP1 was identified as a novel Cltx target which enhances tumor drug uptake. This finding should facilitate tumor selection for chlorotoxin-based therapeutics and diagnostics.

Folate receptor alpha (FOLR1/FRA) is reported to be overexpressed in epithelial ovarian cancers (EOC), especially the serous histotype. Further, while dysregulation of the folate-dependent 1-carbon cycle has been implicated in tumorogenesis, little is known relative to the potential mechanism of action of FOLR1 expression in these processes. We therefore investigated the expression of FOLR1, other folate receptors, and genes within the 1-carbon cycle in samples of EOC, normal ovary and fallopian tube on a custom TaqMan Low Density Array. Also included on this array were known markers of EOC such as MSLN, MUC16 and HE4. While few differences were observed in the expression profiles of genes in the 1-carbon cycle, genes previously considered to be overexpressed in EOC (e.g., FOLR1, MSLN, MUC16 and HE4) showed significantly increased expression when comparing EOC to normal ovary. However, when the comparator was changed to normal fallopian tube, these differences were abolished, supporting the hypothesis that EOC derives from fallopian fimbriae and, further, that markers previously considered to be upregulated or overexpressed in EOC are most likely not of ovarian origin, but fallopian in derivation. Our findings therefore support the hypothesis that the cell of origin of EOC is tubal epithelium.

Recent studies strongly indicate that aberrations in the control of gene expression might contribute to the initiation and progression of Alzheimer's disease (AD). In particular, alternative splicing has been suggested to play a role in spontaneous cases of AD. Previous transcriptome profiling of AD models and patient samples using microarrays delivered conflicting results. This study provides, for the first time, transcriptomic analysis for distinct regions of the AD brain using RNA-Seq next-generation sequencing technology. Illumina RNA-Seq analysis was used to survey transcriptome profiles from total brain, frontal and temporal lobe of healthy and AD post-mortem tissue. We quantified gene expression levels, splicing isoforms and alternative transcript start sites. Gene Ontology term enrichment analysis revealed an overrepresentation of genes associated with a neuron's cytological structure and synapse function in AD brain samples. Analysis of the temporal lobe with the Cufflinks tool revealed that transcriptional isoforms of the apolipoprotein E gene, APOE-001, -002 and -005, are under the control of different promoters in normal and AD brain tissue. We also observed differing expression levels of APOE-001 and -002 splice variants in the AD temporal lobe. Our results indicate that alternative splicing and promoter usage of the APOE gene in AD brain tissue might reflect the progression of neurodegeneration.