Search

BackgroundRandom Forests (RF) are a widely used modelling tool, enabling feature-selection via a variable importance measure. For this, a threshold is required that separates label-associated features from false positives. In the absence of a good understanding of the characteristics of the variable importance measures, current approaches attempt to select features by training multiple random forest to generate statistical power via a permutation null, employ recursive feature elimination or a combination of both. However, for high-dimensional datasets, such as genome data with millions of variables, this is computationally infeasible.MethodHere we present RFlocalfdr, a statistical approach to a threshold that identifies which features are significantly associated with the prediction label and reduces false positives. It builds on the empirical Bayes argument of Efron (2005) and models the variable importance as mixture of two distributions – null and non-null “genes.”ResultWe demonstrate on synthetic data that RFlocalfdr has an equivalent performance to computationally more intensive approaches. But unlike these methods, RFlocalfdr successfully thresholds a 6 Million-feature dataset of 10,000 samples in real-time, being up to 100 times faster than other tools. RFlocalfdr can be applied to the output of any RF implementation returning variable importance and counts and we demonstrate its application to ranger and VariantSpark.ConclusionRFlocalfdr allows for robust feature selection by placing a confidence value on the predicted importance score. It does so without repeated fitting of the RF or the use of additional shadow variables and is thus usable for data sets with very large numbers of variables.

BackgroundAmyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease with phenotypic and genetic heterogeneity. Approximately 10% of cases are familial, while remaining cases are classified as sporadic. To date, >30 genes and several hundred genetic variants have been implicated in ALS.MethodsSeven hundred and fifty-seven sporadic ALS cases were recruited from Australian neurology clinics. Detailed clinical data and whole genome sequencing (WGS) data were available from 567 and 616 cases, respectively, of which 426 cases had both datasets available. As part of a comprehensive genetic analysis, 853 genetic variants previously reported as ALS-linked mutations or disease-associated alleles were interrogated in sporadic ALS WGS data. Statistical analyses were performed to identify correlation between clinical variables, and between phenotype and the number of ALS-implicated variants carried by an individual. Relatedness between individuals carrying identical variants was assessed using identity-by-descent analysis.ResultsForty-three ALS-implicated variants from 18 genes, including C9orf72, ATXN2, TARDBP, SOD1, SQSTM1 and SETX, were identified in Australian sporadic ALS cases. One-third of cases carried at least one variant and 6.82% carried two or more variants, implicating a potential oligogenic or polygenic basis of ALS. Relatedness was detected between two sporadic ALS cases carrying a SOD1 p.I114T mutation, and among three cases carrying a SQSTM1 p.K238E mutation. Oligogenic/polygenic sporadic ALS cases showed earlier age of onset than those with no reported variant.ConclusionWe confirm phenotypic associations among ALS cases, and highlight the contribution of genetic variation to all forms of ALS.

ObjectiveSince the first report of CHCHD10 gene mutations in amyotrophiclateral sclerosis (ALS)/frontotemporaldementia (FTD) patients, genetic variation in CHCHD10 has been inconsistently linked to disease. A pathological assessment of the CHCHD10 protein in patient neuronal tissue also remains to be reported. We sought to characterise the genetic and pathological contribution of CHCHD10 to ALS/FTD in Australia.MethodsWhole-exome and whole-genome sequencing data from 81 familial and 635 sporadic ALS, and 108 sporadic FTD cases, were assessed for genetic variation in CHCHD10. CHCHD10 protein expression was characterised by immunohistochemistry, immunofluorescence and western blotting in control, ALS and/or FTD postmortem tissues and further in a transgenic mouse model of TAR DNA-binding protein 43 (TDP-43) pathology.ResultsNo causal, novel or disease-associated variants in CHCHD10 were identified in Australian ALS and/or FTD patients. In human brain and spinal cord tissues, CHCHD10 was specifically expressed in neurons. A significant decrease in CHCHD10 protein level was observed in ALS patient spinal cord and FTD patient frontal cortex. In a TDP-43 mouse model with a regulatable nuclear localisation signal (rNLS TDP-43 mouse), CHCHD10 protein levels were unaltered at disease onset and early in disease, but were significantly decreased in cortex in mid-stage disease.ConclusionsGenetic variation in CHCHD10 is not a common cause of ALS/FTD in Australia. However, we showed that in humans, CHCHD10 may play a neuron-specific role and a loss of CHCHD10 function may be linked to ALS and/or FTD. Our data from the rNLS TDP-43 transgenic mice suggest that a decrease in CHCHD10 levels is a late event in aberrant TDP-43-induced ALS/FTD pathogenesis.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease caused by the progressive degeneration of motor neurons. Recently, genetic variants in GLT8D1 and ARPP21 were associated with ALS in a cohort of European descent. A synergistic relationship was proposed between ALS associated variants in GLT8D1 and ARPP21. We aimed to determine the prevalence of genetic variation in GLT8D1 and ARPP21 in an Australian cohort of familial (n = 81) and sporadic ALS (n = 618) cases using whole-exome and whole-genome sequencing data. No novel mutations were identified in either gene, nor was there significant enrichment of protein-altering sequence variation among ALS cases. GLT8D1 and ARPP21 mutations are not a common cause of ALS in Australian familial and sporadic cohorts.

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disorder characterised by the loss of upper and lower motor neurons resulting in paralysis and eventual death. Approximately 10% of ALS cases have a family history of disease, while the remainder present as apparently sporadic cases. Heritability studies suggest a significant genetic component to sporadic ALS, and although most sporadic cases have an unknown genetic aetiology, some familial ALS mutations have also been found in sporadic cases. This suggests that some sporadic cases may be unrecognised familial cases with reduced disease penetrance in their ancestors. A powerful strategy to uncover a familial link is identity-by-descent (IBD) analysis, which detects genomic regions that have been inherited from a common ancestor. IBD analysis was performed on 83 Australian familial ALS cases from 25 families and three sporadic ALS cases, each of whom carried one of three SOD1 mutations (p.I114T, p.V149G and p.E101G). We defined five unique 350-SNP haplotypes that carry these mutations in our cohort, indicative of five founder events. This included two founder haplotypes that carry SOD1 p.I114T; linking familial and sporadic cases. We found that SOD1 p.E101G arose independently in each family that carries this mutation and linked two families that carry SOD1 p.V149G. The age of disease onset varied between cases that carried each SOD1 p.I114T haplotype. Linking families with identical ALS mutations allows for larger sample sizes and increased statistical power to identify putative phenotypic modifiers.

Long-term in vitro expansion of bone marrow stromal (skeletal) stem cells (also known as human mesenchymal stem cells [hMSC]) is associated with replicative senescence and impaired functions. We have previously reported that telomerization of hMSC through hTERT overexpression led to bypassing a replicative senescence phenotype and improved in vitro and in vivo functions. However, the molecular consequence of telomerization is poorly characterized. Thus, we compared the molecular phenotype of a well-studied telomerized hMSC (hMSC-TERT) cell line with primary hMSC. At a cellular level, both cell populations exhibited strong concordance for the known hMSC CD markers, similar responses to osteoblast (OB) differentiation induction, and formed heterotopic bone in vivo. Overall gene expression was highly correlated between both cell types with an average Pearson's correlation coefficient (R) between the gene expression of all primary hMSC and all hMSC-TERT samples of 0.95 (range 0.93-0.96). Quantitative analysis of gene expression of CD markers, OB cell markers, and transcription factors (TF) showed a high degree of similarity between the two cell populations (72%, 77%, and 81%, respectively). The hMSC-TERT population was enriched mainly for genes associated with cell cycle and cell cycle signaling when compared with primary hMSC. Other enrichment was observed for genes involved in cell adhesion and skeletal system development and immune response pathways. Interestingly, hMSC-TERT shared a telomerization signature with upregulation of cancer/testis antigens, MAGE, and PAGE genes. Our data demonstrate that the enhanced biological characteristics of hMSC after telomerization are mainly due to enhanced expression of cell proliferation genes, whereas gene expression responses to differentiation are maintained. © 2018 The Authors. Published by Wiley Periodicals, Inc. on behalf of the American Society for Bone and Mineral Research.

BackgroundMany traits and diseases are thought to be driven by >1 gene (polygenic). Polygenic risk scores (PRS) hence expand on genome-wide association studies by taking multiple genes into account when risk models are built. However, PRS only considers the additive effect of individual genes but not epistatic interactions or the combination of individual and interacting drivers. While evidence of epistatic interactions ais found in small datasets, large datasets have not been processed yet owing to the high computational complexity of the search for epistatic interactions.FindingsWe have developed VariantSpark, a distributed machine learning framework able to perform association analysis for complex phenotypes that are polygenic and potentially involve a large number of epistatic interactions. Efficient multi-layer parallelization allows VariantSpark to scale to the whole genome of population-scale datasets with 100,000,000 genomic variants and 100,000 samples.ConclusionsCompared with traditional monogenic genome-wide association studies, VariantSpark better identifies genomic variants associated with complex phenotypes. VariantSpark is 3.6 times faster than ReForeSt and the only method able to scale to ultra-high-dimensional genomic data in a manageable time.

SummaryTRIBES is a user-friendly pipeline for relatedness detection in genomic data. TRIBES is the first tool which is both accurate up to 7th degree relatives (e.g. third cousins) and combines essential data processing steps into a single user-friendly pipeline. Furthermore, using a proof-of-principle cohort comprising amyotrophic lateral sclerosis cases with known relationship structures and a known causal mutation in SOD1, we demonstrated that TRIBES can successfully uncover disease susceptibility loci. TRIBES has multiple applications in addition to disease gene mapping, including sample quality control in genome wide association studies and avoiding consanguineous unions in family planning.AvailabilityTRIBES is freely available on GitHub: https://github.com/aehrc/TRIBES/Contactnatalie.twine@csiro.auSupplementary informationXXXX

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disorder characterised by the loss of upper and lower motor neurons resulting in paralysis and eventual death. Approximately 10% of ALS cases have a family history of disease, while the remaining cases present as apparently sporadic. Heritability studies suggest a significant genetic component to sporadic ALS, and although most sporadic cases have an unknown genetic etiology, some familial ALS mutations have also been found in sporadic cases. This suggests that some sporadic cases may be unrecognised familial cases with reduced disease penetrance. Identifying a familial basis of disease in apparently sporadic ALS cases has significant genetic counselling implications for immediate relatives. A powerful strategy to uncover a familial link is identity-by-descent (IBD) analysis which detects genomic regions that have been inherited from a common ancestor. We performed IBD analysis on 90 Australian familial ALS cases from 25 families and three sporadic ALS cases, each of whom carried one of three SOD1 mutations (p.I114T, p.V149G and p.E101G). We identified five unique haplotypes that carry these mutations in our cohort, indicative of five founder events. This included two different haplotypes that carry SOD1 p.I114T, where one haplotype was present in one sporadic case and 20 families, while the second haplotype was found in the remaining two sporadic cases and one family, thus linking these familial and sporadic cases. Furthermore, we linked two families that carry SOD1 p.V149G and found that SOD1 p.E101G arose independently in each family that carries this mutation.

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. Electronic supplementary material The online version of this article (10.1186/s12964-019-0368-9) contains supplementary material, which is available to authorized users.

This article provides an overview of Artificial Intelligence and Machine Learning methodologies and surveys representative applications in the bioinformatics field, such as for genomics, genetics, proteomics, systems biology and metagenomics. It provides examples of how to choose the appropriate Machine Learning algorithm and discusses Machine Learning in the era of “Big Data”.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease characterised by the loss of upper and lower motor neurons. ALS exhibits high phenotypic variability including age and site of onset, and disease duration. To uncover epigenetic and transcriptomic factors that may modify an ALS phenotype, we used a cohort of Australian monozygotic twins (n = 3 pairs) and triplets (n = 1 set) that are discordant for ALS and represent sporadic ALS and the two most common types of familial ALS, linked to C9orf72 and SOD1. Illumina Infinium HumanMethylation450K BeadChip, EpiTYPER and RNA-Seq analyses in these ALS-discordant twins/triplets and control twins (n = 2 pairs), implicated genes with consistent longitudinal differential DNA methylation and/or gene expression. Two identified genes, RAD9B and C8orf46, showed significant differential methylation in an extended cohort of >1000 ALS cases and controls. Combined longitudinal methylation-transcription analysis within a single twin set implicated CCNF, DPP6, RAMP3, and CCS, which have been previously associated with ALS. Longitudinal transcriptome data showed an 8-fold enrichment of immune function genes and under-representation of transcription and protein modification genes in ALS. Examination of these changes in a large Australian sporadic ALS cohort suggest a broader role in ALS. Furthermore, we observe that increased methylation age is a signature of ALS in older patients.

© 2018 Elsevier Inc., To identify novel genes associated with ALS, we undertook two lines of investigation. We carried out a genome-wide association study comparing 20,806 ALS cases and 59,804 controls. Independently, we performed a rare variant burden analysis comparing 1,138 index familial ALS cases and 19,494 controls. Through both approaches, we identified kinesin family member 5A (KIF5A) as a novel gene associated with ALS. Interestingly, mutations predominantly in the N-terminal motor domain of KIF5A are causative for two neurodegenerative diseases: hereditary spastic paraplegia (SPG10) and Charcot-Marie-Tooth type 2 (CMT2). In contrast, ALS-associated mutations are primarily located at the C-terminal cargo-binding tail domain and patients harboring loss-of-function mutations displayed an extended survival relative to typical ALS cases. Taken together, these results broaden the phenotype spectrum resulting from mutations in KIF5A and strengthen the role of cytoskeletal defects in the pathogenesis of ALS.

Inhibition of Poly(ADP-ribose) Polymerase1 (PARP1) impairs DNA damage repair, and early generation PARP1/2 inhibitors (olaparib, niraparib, etc.) have demonstrated clinical proof of concept for cancer treatment. Here, we describe the development of the novel PARP inhibitor E7449, a potent PARP1/2 inhibitor that also inhibits PARP5a/5b, otherwise known as tankyrase1 and 2 (TNKS1 and 2), important regulators of canonical Wnt/β-catenin signaling. E7449 inhibits PARP enzymatic activity and additionally traps PARP1 onto damaged DNA; a mechanism previously shown to augment cytotoxicity. Cells deficient in DNA repair pathways beyond homologous recombination were sensitive to E7449 treatment. Chemotherapy was potentiated by E7449 and single agent had significant antitumor activity in BRCA-deficient xenografts. Additionally, E7449 inhibited Wnt/β-catenin signaling in colon cancer cell lines, likely through TNKS inhibition. Consistent with this possibility, E7449 stabilized axin and TNKS proteins resulting in β-catenin de-stabilization and significantly altered expression of Wnt target genes. Notably, hair growth mediated by Wnt signaling was inhibited by E7449. A pharmacodynamic effect of E7449 on Wnt target genes was observed in tumors, although E7449 lacked single agent antitumor activity in vivo, a finding typical for selective TNKS inhibitors. E7449 antitumor activity was increased through combination with MEK inhibition. Particularly noteworthy was the lack of toxicity, most significantly the lack of intestinal toxicity reported for other TNKS inhibitors. E7449 represents a novel dual PARP1/2 and TNKS1/2 inhibitor which has the advantage of targeting Wnt/β-catenin signaling addicted tumors. E7449 is currently in early clinical development.

Long‐term in vitro expansion of bone marrow stromal (skeletal) stem cells (also known as human mesenchymal stem cells [hMSC]) is associated with replicative senescence and impaired functions. We have previously reported that telomerization of hMSC through hTERT overexpression led to bypassing a replicative senescence phenotype and improved in vitro and in vivo functions. However, the molecular consequence of telomerization is poorly characterized. Thus, we compared the molecular phenotype of a well‐studied telomerized hMSC (hMSC‐TERT) cell line with primary hMSC. At a cellular level, both cell populations exhibited strong concordance for the known hMSC CD markers, similar responses to osteoblast (OB) differentiation induction, and formed heterotopic bone in vivo. Overall gene expression was highly correlated between both cell types with an average Pearson's correlation coefficient (R2) between the gene expression of all primary hMSC and all hMSC‐TERT samples of 0.95 (range 0.93–0.96). Quantitative analysis of gene expression of CD markers, OB cell markers, and transcription factors (TF) showed a high degree of similarity between the two cell populations (72%, 77%, and 81%, respectively). The hMSC‐TERT population was enriched mainly for genes associated with cell cycle and cell cycle signaling when compared with primary hMSC. Other enrichment was observed for genes involved in cell adhesion and skeletal system development and immune response pathways. Interestingly, hMSC‐TERT shared a telomerization signature with upregulation of cancer/testis antigens, MAGE, and PAGE genes. Our data demonstrate that the enhanced biological characteristics of hMSC after telomerization are mainly due to enhanced expression of cell proliferation genes, whereas gene expression responses to differentiation are maintained. © 2018 The Authors. JBMR Plus Published by Wiley Periodicals, Inc. on behalf of the American Society for Bone and Mineral Research

Reprogramming of immunosuppressive tumor microenvironment (TME) by targeting alternatively activated tumor associated macrophages (M2TAM), myeloid-derived suppressor cells (MDSC), and regulatory T cells (Tregs), represents a promising strategy for developing novel cancer immunotherapy. Prostaglandin E2 (PGE2), an arachidonic acid pathway metabolite and mediator of chronic inflammation, has emerged as a powerful immunosuppressor in the TME through engagement with one or more of its 4 receptors (EP1-EP4). We have developed E7046, an orally bioavailable EP4-specific antagonist and show here that E7046 has specific and potent inhibitory activity on PGE2-mediated pro-tumor myeloid cell differentiation and activation. E7046 treatment reduced the growth or even rejected established tumors in vivo in a manner dependent on both myeloid and CD8+ T cells. Furthermore, co-administration of E7046 and E7777, an IL-2-diphtheria toxin fusion protein that preferentially kills Tregs, synergistically disrupted the myeloid and Treg immunosuppressive networks, resulting in effective and durable anti-tumor immune responses in mouse tumor models. In the TME, E7046 and E7777 markedly increased ratios of CD8+granzymeB+ cytotoxic T cells (CTLs)/live Tregs and of M1-like/M2TAM, and converted a chronic inflammation phenotype into acute inflammation, shown by substantial induction of STAT1/IRF-1 and IFNγ-controlled genes. Notably, E7046 also showed synergistic anti-tumor activity when combined with anti-CTLA-4 antibodies, which have been reported to diminish intratumoral Tregs. Our studies thus reveal a specific myeloid cell differentiation-modifying activity by EP4 blockade and a novel combination of E7046 and E7777 as a means to synergistically mitigate both myeloid and Treg-derived immunosuppression for cancer treatment in preclinical models.

Diatoms are of enormous ecological importance as they account for as much as 20% of global primary production, yet they are still understudied from a genomic perspective. The benthic diatom Ceratoneis closterium is well-characterized from an ecotoxicological perspective including its use in ecotoxicological risk assessments and investigating the mode-of-action of metal toxicity. However, this organism has little sequence information available. In this study, 454 pyrosequencing of the stressor-responsive transcriptome was undertaken. These transcripts could be used to characterize general physiological processes such as photosynthesis and respiration, as well as to enable a description of the ecotoxicogenomic responses of this organism. After a 96 h exposure to the concentration of toxicant that inhibited growth rate by 10% (IC 10 ) for the following common coastal contaminants: ammonia, copper, crude oil and simazine (a photosystem II inhibiting herbicide), diatom cells were harvested for RNA extraction and their transcriptomes characterized via 454 pyrosequencing. This resulted in 1.25 million reads, which were assembled into 4768 contigs, when contigs encoding rRNA were removed. More than 80% of the remaining contigs had an ortholog in the BLASTx protein databases. These contigs represented 1660 unique transcripts. The role of these transcripts in stress response, as well as photosynthesis and respiration is discussed. Overall, this study greatly enhances the genomic information available for this important taxonomic group.

Next generation sequencing using Roche's 454 pyrosequencing platform can be used to generate genomic information for non-model organisms, although there are bioinformatic challenges associated with these studies. These challenges are compounded by a lack of a standardized protocol to either assemble data or to evaluate the quality of a de novo transcriptome. This study presents an assembly of the control and toxicant responsive transcriptome of Melita plumulosa, an Australian amphipod commonly used in ecotoxicological studies. RNA was harvested from control amphipods, juvenile amphipods, and from amphipods exposed to either metal or diesel contaminated sediments. This RNA was used as the basis for a 454 based transcriptome sequencing effort. Sequencing generated 1.3 million reads from control, juvenile, metal-exposed and diesel-exposed amphipods. Different read filtering and assembly protocols were evaluated to generate an assembly that (i) had an optimal number of contigs; (ii) had long contigs; (iii) contained a suitable representation of conserved genes; and (iv) had long ortholog alignment lengths relative to the length of each contig. A final assembly, generated using fixed-length trimming based on the sequence quality scores, followed by assembly using the MIRA algorithm, produced the best results. The 26,625 contigs generated via this approach were annotated using Blast2GO, and the differential expression between treatments and control was determined by mapping with BWA followed by DESeq. Although the mapping generated low coverage, many differentially expressed contigs, including some with known developmental or toxicological function, were identified. This study demonstrated that 454 pyrosequencing is an effective means of generating reference transcriptome information for organisms, such as the amphipod M. plumulosa, that have no genomic information available in databases or in closely related sequenced species. It also demonstrated how optimization of read filtering protocols and assembly approaches changes the utility of results obtained from next generation sequencing studies, and establishes criteria to determine the quality of a de novo assembly in species lacking a reference genome. This new transcriptomic knowledge provides the genomic foundation for the creation of microarray and qPCR assays, serving as a reference transcriptome in future RNAseq studies, and allowing both the biology and ecotoxicology of this organism to be better understood. This approach will allow genomics-based methodology to be applied to a wider range of environmentally relevant species.

Sydney rock oysters (Saccostrea glomerata) were exposed to environmental stressors at contaminated field sites or in a controlled laboratory setting. RNA seq transcriptome data were generated for the gill and digestive gland using Roche's 454 pyrosequencing technology. 28,685 contigs were de novo assembled which encoded 11,671 different protein products. The data will act as a reference for future studies in ecology, immunology and environmental toxicology.

Krüppel-like factors 3 and 8 (KLF3 and KLF8) are highly related transcriptional regulators that bind to similar sequences of DNA. We have previously shown that in erythroid cells there is a regulatory hierarchy within the KLF family, whereby KLF1 drives the expression of both the Klf3 and Klf8 genes and KLF3 in turn represses Klf8 expression. While the erythroid roles of KLF1 and KLF3 have been explored, the contribution of KLF8 to this regulatory network has been unknown. To investigate this, we have generated a mouse model with disrupted KLF8 expression. While these mice are viable, albeit with a reduced life span, mice lacking both KLF3 and KLF8 die at around embryonic day 14.5 (E14.5), indicative of a genetic interaction between these two factors. In the fetal liver, Klf3 Klf8 double mutant embryos exhibit greater dysregulation of gene expression than either of the two single mutants. In particular, we observe derepression of embryonic, but not adult, globin expression. Taken together, these results suggest that KLF3 and KLF8 have overlapping roles in vivo and participate in the silencing of embryonic globin expression during development.

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.

Meeting abstracts Immunotherapies targeting the immune checkpoint receptors have shown great promise for a subset of cancer patients. However, robust and safe combination therapies are still needed to increase the benefit of cancer immunotherapy and bring it to broader patient populations. We have

Human proteome analysis now requires an understanding of protein isoforms. We recently published the PG Nexus pipeline, which facilitates high confidence validation of exons and splice junctions by integrating genomics and proteomics data. Here we comprehensively explore how RNA-seq transcriptomics data, and proteomic analysis of the same sample, can identify protein isoforms. RNA-seq data from human mesenchymal (hMSC) stem cells were analyzed with our new TranscriptCoder tool to generate a database of protein isoform sequences. MS/MS data from matching hMSC samples were then matched against the TranscriptCoder-derived database, along with Ensembl and the neXtProt database. Querying the TranscriptCoder-derived or Ensembl database could unambiguously identify ∼450 protein isoforms, with isoform-specific proteotypic peptides, including candidate hMSC-specific isoforms for the genes DPYSL2 and FXR1. Where isoform-specific peptides did not exist, groups of nonisoform-specific proteotypic peptides could specifically identify many isoforms. In both the above cases, isoforms will be detectable with targeted MS/MS assays. Unfortunately, our analysis also revealed that some isoforms will be difficult to identify unambiguously as they do not have peptides that are sufficiently distinguishing. We covisualize mRNA isoforms and peptides in a genome browser to illustrate the above situations. Mass spectrometry data is available via ProteomeXchange (PXD001449).

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., Highlights • hESC-derived MSC-like cells were compared to immortalised BM-MSC. • Comparison was performed using microarrays on non-induced and OB induced cells. • Analysis demonstrated close hierarchical relationships and molecular phenotypes. • 90.7% of genes were similarly expressed in non-induced cells. • 73% of OB induced genes for both cell lines correlated with GO ontology analysis.

Purpose: Given their accessibility, surrogate tissues, such as peripheral blood mononuclear cells (PBMC), may provide potential predictive biomarkers in clinical pharmacogenomic studies. In leukemias and lymphomas, the prognostic value of peripheral blast expression profiles is clear; however, it is unclear whether circulating mononuclear cells of patients with solid tumors might yield profiles with similar prognostic associations. Experimental Design: In this study, we evaluated the association of expression profiles in PBMCs with clinical outcomes in patients with advanced renal cell cancer. Transcriptional patterns in PBMCs of 45 renal cell cancer patients were compared with clinical outcome data at the conclusion of a phase II study of the mTOR kinase inhibitor CCI-779 to determine whether pretreatment transcriptional patterns in PBMCs were correlated with eventual patient outcomes. Results: Unsupervised hierarchical clustering of the PBMC profiles using all expressed genes identified clusters of patients with significant differences in survival. Cox proportional hazards modeling showed that the expression levels of many PBMC transcripts were predictors for the patient outcomes of time to progression and overall survival (time to death). Supervised class prediction approaches identified multivariate expression patterns in PBMCs capable of assigning favorable outcomes of time to death and time to progression in a test set of renal cancer patients, with overall performance accuracies of 72% and 85%, respectively. Conclusions: The present study provides the first example of gene expression profiling in peripheral blood, a clinically accessible surrogate tissue, for identifying patterns of gene expression associated with higher likelihoods of positive outcome in patients with a solid tumor.

Objective Our objective was to estimate the pharmacokinetic parameters of CCI-779 and its metabolite, sirolimus, and evaluate associations of exposure parameters with safety and clinical activity. Exposure parameters were also correlated with pharmacogenomic responses in peripheral blood mononuclear cells (PBMCs). Methods In this randomized, double-blind, multicenter trial, once-weekly intravenous doses of 25, 75, or 250 mg CCI-779 were administered to patients with advanced renal cancer. Whole blood for CCI-779 and sirolimus concentrations was drawn. Population pharmacokinetic analyses yielded Bayesian-predicted exposure metrics that were correlated with severity and duration of adverse events and survival. PBMC samples taken before and after treatment were examined for pharmacogenomic responses. Ribonucleic acid samples were converted to labeled probes and hybridized to oligonucleotide arrays containing more than 12,600 human sequences. Results The final population pharmacokinetic models of CCI-779 and sirolimus included 235 and 305 observations, respectively, from 50 patients. For CCI-779, dose, single versus multiple dose, and body surface area were significant pharmacokinetic covariates. For sirolimus, dose and hematocrit were significant covariates. Age, sex, or race did not influence drug disposition. CCI-779 area under the curve correlated with adverse event severity for thrombocytopenia (P = .007), pruritus (P = .011), and hyperlipemia (P = .040). Exposure (CCI-779 cumulative area under the curve) correlated with a specific subset of gene transcripts in PBMCs following 16 weeks after therapy (P < .001, Spearman correlation). Conclusions Concentrations of CCI-779 and sirolimus were adequately described with a population model incorporating factors for dose, attenuated exposure of multiple doses, body surface area, and hematocrit. Correlations with adverse event severity and duration profiles were provided to aid in the detection of treatment-emergent effects. Pharmacogenomic profiling of PBMCs identified altered ribonucleic acid transcript expression levels that correlate with exposure. These transcripts represent potential biomarkers of CCI-779 exposure in peripheral blood. Clinical Pharmacology & Therapeutics (2005) 77, 76–89; doi: 10.1016/j.clpt.2004.08.025

Objective Levels of interleukin-6 (IL-6) have been shown to correlate with the fever and disease activity of systemic juvenile idiopathic arthritis (JIA). In a previous case–control study, a significant association between the IL-6 −174 nucleotide variant and systemic JIA was noted, and HeLa cell transfection assays show functional differences in levels of transcription of the IL-6 −174 alleles. The present study was undertaken to confirm the previous findings and to assess possible association with variations of the AnTn tract in the promoter. Methods We studied a cohort of JIA families from 3 countries, using transmission disequilibrium testing. Genotyping of the −174 nucleotide variant was done by restriction fragment length polymorphism, heteroduplex analysis, or allelic discrimination. The AnTn tract at −392 to −373 was typed using DNA sequencing. Statistical analysis was performed using the programs Transmit and EHplus. Results There was a significant excess transmission of the −174G allele in the systemic JIA families (P = 0.041). The excess transmission was only to systemic JIA patients with age at onset >5 years (P = 0.007). No significant association with the other subtypes was found. No AnTn alleles or −174/AnTn haplotypes were significantly associated with systemic JIA. Conclusion This study confirms that the IL-6 –174 nucleotide variant is significantly associated with systemic JIA. The significant excess transmission to patients with age at onset >5 years but not to those with age at onset ≤5 years suggests that there may be genetic heterogeneity between the 2 groups.

The phenotype of osteoblastic (OB) cells in culture is currently defined using a limited number of markers of low sensitivity and specificity. For the clinical use of human skeletal (stromal, mesenchymal) stem cells (hMSC) in therapy, there is also a need to identify a set of gene markers that predict in vivo bone forming capacity. Thus, we used RNA sequencing to examine changes in expression for a set of skeletally-related genes across 8 time points between 0 and 12days of ex vivo OB differentiation of hMSC. We identified 123 genes showing significant temporal expression change. Hierarchical clustering and Pearson's correlation generated 4 groups of genes: early stage differentiation genes (peak expression: 0-24h, n=28) which were enriched for extracellular matrix organisation, e.g. COL1A1, LOX, and SERPINH1; middle stage differentiating genes (peak expression days: 3 and 6, n=20) which were enriched for extracellular matrix/skeletal system development e.g. BMP4, CYP24A1, and TGFBR2; and late stage differentiation genes (peak expression days: 9 and 12, n=27) which were enriched for bone development/osteoblast differentiation, e.g. BMP2 and IGF2. In addition, we identified 13 genes with bimodal temporal expression (2 peaks of expression: days 0 and 12) including VEGFA, PDGFA and FGF2. We examined the specificity of the 123 genes' expression in skeletal tissues and thus propose a set of ex vivo differentiation-stage-specific markers (n=21). In an independent analysis, we identified a subset of genes (n=20, e.g. ELN, COL11A1, BMP4) to predict the bone forming capacity of hMSC and another set (n=20, e.g. IGF2, TGFB2, SMAD3) associated with the ex vivo phenotype of hMSC obtained from osteoporotic patients.

The discovery that circulating nucleic acid-containing complexes in the serum of autoimmune lupus patients can stimulate B cells and plasmacytoid dendritic cells via Toll-like receptors 7 and 9 suggested that agents that block these receptors might be useful therapeutics. We identified two compounds, AT791 {3-[4-(6-(3-(dimethylamino)propoxy)benzo[d]oxazol-2-yl)phenoxy]-N,N-dimethylpropan-1-amine} and E6446 {6-[3-(pyrrolidin-1-yl)propoxy)-2-(4-(3-(pyrrolidin-1-yl)propoxy)phenyl]benzo[d]oxazole}, that inhibit Toll-like receptor (TLR)7 and 9 signaling in a variety of human and mouse cell types and inhibit DNA-TLR9 interaction in vitro. When administered to mice, these compounds suppress responses to challenge doses of cytidine-phosphate-guanidine (CpG)-containing DNA, which stimulates TLR9. When given chronically in spontaneous mouse lupus models, E6446 slowed development of circulating antinuclear antibodies and had a modest effect on anti-double-stranded DNA titers but showed no observable impact on proteinuria or mortality. We discovered that the ability of AT791 and E6446 to inhibit TLR7 and 9 signaling depends on two properties: weak interaction with nucleic acids and high accumulation in the intracellular acidic compartments where TLR7 and 9 reside. Binding of the compounds to DNA prevents DNA-TLR9 interaction in vitro and modulates signaling in vivo. Our data also confirm an earlier report that this same mechanism may explain inhibition of TLR7 and 9 signaling by hydroxychloroquine (Plaquenil; Sanofi-Aventis, Bridgewater, NJ), a drug commonly prescribed to treat lupus. Thus, very different structural classes of molecules can inhibit endosomal TLRs by essentially identical mechanisms of action, suggesting a general mechanism for targeting this group of TLRs.

Introduction: In a small neoadjuvant study in patients with triple negative breast cancer (TNBC) the combination of eribulin plus carboplatin was effective, with a pathologic complete response rate of 43% following 4 cycles of treatment. Significant numbers of sporadic TNBC tumors are deficient in DNA repair capacity and share clinical and pathological features with hereditary BRCA1 mutant disease. PARP inhibitors have demonstrated synthetic lethality in cancer cells with defective DNA repair and have therapeutic potential for TNBC. In this study we describe the combination of PARP inhibitor E7449 with eribulin +/- carboplatin in preclinical models of TNBC. Methods: E7449, an orally available PARP inhibitor, was administered in combination with eribulin +/- carboplatin to 4 s.c. xenograft models of TNBC: MDA-MB-436 (BRCA1 mutant, PTEN deficient), MDA-MB-468 (BRCA wild type, PTEN deficient), HCC1806 and MDA-MB-231 (BRCA and PTEN wild type). Results and Discussion: Addition of E7449 to eribulin significantly delayed tumor progression in PTEN deficient MDA-MB-468 xenografts. In the BRCA1 mutant and PTEN deficient MDA-MB-436 xenograft model, combination of E7449 with eribulin enhanced antitumor activity versus eribulin alone. Similar potentiation was observed for carboplatin upon combination with E7449. Treatment of MDA-MB-436 xenografts with the triple combination of E7449 + eribulin + carboplatin was more efficacious than any double combination and was well tolerated at the doses examined. In contrast, no significant combination activity was observed for E7449 plus eribulin in the BRCA and PTEN wild type xenografts HCC1806 and MDA-MB-231, and similarly no potentiation of carboplatin was observed in an MDA-MB-231 xenograft. Notably, combination activity was observed in the BRCA1 mutant (MDA-MB-436) and PTEN deficient (MDA-MB-436 and MDA-MB-468) xenografts and not in the BRCA and PTEN wild-type models (HCC1806 and MDA-MB-231). Data from ongoing studies to evaluate the combination activity of E7449 + eribulin in patient-derived xenograft (PDx) models of TNBC will be presented at the meeting. Potential biomarkers of sensitivity to the combination are under investigation in both cell line xenograft and PDx models and will be described. Conclusion: The addition of E7449 to eribulin +/- carboplatin increased antitumor activity in a subset of TNBC models. Biomarker studies aimed at a better understanding of the underlying cause of sensitivity are underway. The preclinical data support assessment of E7449 + eribulin + carboplatin combination therapy in the current phase I/II clinical trial. Citation Format: Sharon McGonigle, Jiayi Wu, Donna Kolber-Simonds, Natalie C Twine, Jue-lon Shie, Noel Taylor, Sergei Agoulnik, Zoltan Dezso, Shannon McGrath, Mark Matijevic, Shanqin Xu, Galina Kuznetsov, Mary Woodall-Jappe, Kenichi Nomoto. Combination of the PARP inhibitor E7449 with eribulin +/- carboplatin in preclinical models of triple negative breast cancer [abstract]. In: Proceedings of the Thirty-Seventh Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2014 Dec 9-13; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2015;75(9 Suppl):Abstract nr P5-06-03.

Direct links between proteomic and genomic/transcriptomic data are not frequently made, partly because of lack of appropriate bioinformatics tools. To help address this, we have developed the PG Nexus pipeline. The PG Nexus allows users to covisualize peptides in the context of genomes or genomic contigs, along with RNA-seq reads. This is done in the Integrated Genome Viewer (IGV). A Results Analyzer reports the precise base position where LC-MS/MS-derived peptides cover genes or gene isoforms, on the chromosomes or contigs where this occurs. In prokaryotes, the PG Nexus pipeline facilitates the validation of genes, where annotation or gene prediction is available, or the discovery of genes using a "virtual protein"-based unbiased approach. We illustrate this with a comprehensive proteogenomics analysis of two strains of Campylobacter concisus . For higher eukaryotes, the PG Nexus facilitates gene validation and supports the identification of mRNA splice junction boundaries and splice variants that are protein-coding. This is illustrated with an analysis of splice junctions covered by human phosphopeptides, and other examples of relevance to the Chromosome-Centric Human Proteome Project. The PG Nexus is open-source and available from https://github.com/IntersectAustralia/ap11_Samifier . It has been integrated into Galaxy and made available in the Galaxy tool shed.

The hippocampus and cerebellum represent anatomically and functionally distinct parts of the human brain. The RNA-Seq technique makes it possible to investigate the human transcriptome with unprecedented resolution, allowing identification of differential mRNA splicing and promoter usage on a genome-wide scale. We undertook whole-mRNA sequencing of samples from the human hippocampus and cerebellum. A bioinformatic analysis revealed distinct expression patterns of genes related to the molecular physiology of neurons and glial cells. Upregulated genes in hippocampal tissue included serpin peptidase inhibitor, clade A (SERPINA3), lymphocyte antigen 6 complex, locus H (LY6H) and transthyretin (TTR). In cerebellum, the cerebellin 3 precursor (CLBN3) and Zic family member 4 (ZIC4) genes were significantly upregulated. These changes were validated in independent donor samples by qRT-PCR. The hippocampus and the cerebellum showed striking differences in splicing patterns and promoter usage. A notable example of this was the gene for NGFI-A binding protein 2 (NAB2), which displayed tissue-specific isoforms which may affect its function as a transcriptional repressor.

The CACCC-box binding protein erythroid Kruppel-like factor (EKLF/KLF1) is a master regulator that directs the expression of many important erythroid genes. We have previously shown that EKLF drives transcription of the gene for a second KLF, basic Kruppel-like factor, or KLF3. We have now tested the in vivo role of KLF3 in erythroid cells by examining Klf3 knockout mice. KLF3-deficient adults exhibit a mild compensated anemia, including enlarged spleens, increased red pulp, and a higher percentage of erythroid progenitors, together with elevated reticulocytes and abnormal erythrocytes in the peripheral blood. Impaired erythroid maturation is also observed in the fetal liver. We have found that KLF3 levels rise as erythroid cells mature to become TER119(+). Consistent with this, microarray analysis of both TER119(-) and TER119(+) erythroid populations revealed that KLF3 is most critical at the later stages of erythroid maturation and is indeed primarily a transcriptional repressor. Notably, many of the genes repressed by KLF3 are also known to be activated by EKLF. However, the majority of these are not currently recognized as erythroid-cell-specific genes. These results reveal the molecular and physiological function of KLF3, defining it as a feedback repressor that counters the activity of EKLF at selected target genes to achieve normal erythropoiesis.

We investigated functional epigenetic changes that occur in primary human T lymphocytes during entry into the cell cycle and mapped these at the single-nucleosome level by ChIP-chip on tiling arrays for chromosomes 1 and 6. We show that nucleosome loss and flanking active histone marks define active transcriptional start sites (TSSs). Moreover, these signatures are already set at many inducible genes in quiescent cells prior to cell stimulation. In contrast, there is a dearth of the inactive histone mark H3K9me3 at the TSS, and under-representation of H3K9me2 and H3K9me3 defines the body of active genes. At the DNA level, cytosine methylation (meC) is enriched for nucleosomes that remain at the TSS, whereas in general there is a dearth of meC at TSSs. Furthermore, a drop in meC also marks 3′ transcription termination, and a peak of meC occurs at stop codons. This mimics the 3′ nucleosomal distribution in yeast, which we show does not occur in human T cells.

We report here that des-methyl, des-amino pateamine A (DMDA-PatA), a structurally simplified analogue of the marine natural product pateamine A, has potent antiproliferative activity against a wide variety of human cancer cell lines while showing relatively low cytotoxicity against nonproliferating, quiescent human fibroblasts. DMDA-PatA retains almost full in vitro potency in P-glycoprotein-overexpressing MES-SA/Dx5-Rx1 human uterine sarcoma cells that are significantly resistant to paclitaxel, suggesting that DMDA-PatA is not a substrate for P-glycoprotein-mediated drug efflux. Treatment of proliferating cells with DMDA-PatA leads to rapid shutdown of DNA synthesis in the S phase of the cell cycle. Cell-free studies show that DMDA-PatA directly inhibits DNA polymerases α and γ in vitro albeit at concentrations considerably higher than those that inhibit cell proliferation. DMDA-PatA shows potent anticancer activity in several human cancer xenograft models in nude mice, including significant regressions observed in the LOX and MDA-MB-435 melanoma models. DMDA-PatA thus represents a promising natural product-based anticancer agent that warrants further investigation.[Mol Cancer Ther 2009;8(5):1250–60]

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.

Low-risk myelodysplastic syndrome (MDS) with normal cytogenetics accounts for approximately 50% of MDS patients. There are no pathognomonic markers in these cases and the diagnosis rests on cytomorphologic abnormalities in bone marrow and/or peripheral blood. Affymetrix high-resolution single-nucleotide polymorphism (SNP) genotyping microarrays allow detection of cytogenetically cryptic genomic aberrations. We have studied 119 low-risk MDS patients (refractory anemia [RA] = 22; refractory cytopenia with multilineage dysplasia [RCMD] = 51; refractory anemia with ringed sideroblasts [RARS] = 12; refractory cytopenia with multilineage dysplasia with ringed sideroblasts [RCMD-RS] = 12; 5q− syndrome = 16; refractory anemia with excess blasts [RAEB] = 6) using SNP microarrays to seek chromosomal markers undetected by conventional cytogenetics. Loss of heterozygosity (LOH) detected by 50K arrays was verified using 250K and 500K arrays. We demonstrate the presence of uniparental disomy (UPD) in 46%, deletions in 10%, and amplifications in 8% of cases. Copy number (CN) changes were acquired, whereas UPDs were also detected in constitutional DNA. UPD on 4q was identified in 25% of RARS, 12% of RCMD with normal cytogenetics, 17% of RAEB, and 6% of 5q− syndrome cases. Univariate analysis showed deletions (P = .04) and International Prognostic Scoring System (IPSS; P < .001) scores correlated with overall survival; however, on multivariate analysis only IPSS scores retained prognostic significance (P < .001). We show, for the first time, that SNP microarray analysis in low-risk MDS patients reveals hitherto unrecognized UPD and CN changes that may allow stratification of these patients for early therapeutic interventions.