School of Mathematics and Statistics - Research Publications

Permanent URI for this collection

Search Results

Now showing 1 - 10 of 44
  • Item
    Thumbnail Image
    G protein-linked signaling pathways in bipolar and major depressive disorders.
    Tomita, H ; Ziegler, ME ; Kim, HB ; Evans, SJ ; Choudary, PV ; Li, JZ ; Meng, F ; Dai, M ; Myers, RM ; Neal, CR ; Speed, TP ; Barchas, JD ; Schatzberg, AF ; Watson, SJ ; Akil, H ; Jones, EG ; Bunney, WE ; Vawter, MP (Frontiers Media SA, 2013)
    The G-protein linked signaling system (GPLS) comprises a large number of G-proteins, G protein-coupled receptors (GPCRs), GPCR ligands, and downstream effector molecules. G-proteins interact with both GPCRs and downstream effectors such as cyclic adenosine monophosphate (cAMP), phosphatidylinositols, and ion channels. The GPLS is implicated in the pathophysiology and pharmacology of both major depressive disorder (MDD) and bipolar disorder (BPD). This study evaluated whether GPLS is altered at the transcript level. The gene expression in the dorsolateral prefrontal (DLPFC) and anterior cingulate (ACC) were compared from MDD, BPD, and control subjects using Affymetrix Gene Chips and real time quantitative PCR. High quality brain tissue was used in the study to control for confounding effects of agonal events, tissue pH, RNA integrity, gender, and age. GPLS signaling transcripts were altered especially in the ACC of BPD and MDD subjects. Transcript levels of molecules which repress cAMP activity were increased in BPD and decreased in MDD. Two orphan GPCRs, GPRC5B and GPR37, showed significantly decreased expression levels in MDD, and significantly increased expression levels in BPD. Our results suggest opposite changes in BPD and MDD in the GPLS, "activated" cAMP signaling activity in BPD and "blunted" cAMP signaling activity in MDD. GPRC5B and GPR37 both appear to have behavioral effects, and are also candidate genes for neurodegenerative disorders. In the context of the opposite changes observed in BPD and MDD, these GPCRs warrant further study of their brain effects.
  • Item
    No Preview Available
    Evolving gene/transcript definitions significantly alter the interpretation of GeneChip data
    Dai, MH ; Wang, PL ; Boyd, AD ; Kostov, G ; Athey, B ; Jones, EG ; Bunney, WE ; Myers, RM ; Speed, TP ; Akil, H ; Watson, SJ ; Meng, F (OXFORD UNIV PRESS, 2005)
    Genome-wide expression profiling is a powerful tool for implicating novel gene ensembles in cellular mechanisms of health and disease. The most popular platform for genome-wide expression profiling is the Affymetrix GeneChip. However, its selection of probes relied on earlier genome and transcriptome annotation which is significantly different from current knowledge. The resultant informatics problems have a profound impact on analysis and interpretation the data. Here, we address these critical issues and offer a solution. We identified several classes of problems at the individual probe level in the existing annotation, under the assumption that current genome and transcriptome databases are more accurate than those used for GeneChip design. We then reorganized probes on more than a dozen popular GeneChips into gene-, transcript- and exon-specific probe sets in light of up-to-date genome, cDNA/EST clustering and single nucleotide polymorphism information. Comparing analysis results between the original and the redefined probe sets reveals approximately 30-50% discrepancy in the genes previously identified as differentially expressed, regardless of analysis method. Our results demonstrate that the original Affymetrix probe set definitions are inaccurate, and many conclusions derived from past GeneChip analyses may be significantly flawed. It will be beneficial to re-analyze existing GeneChip data with updated probe set definitions.
  • Item
    No Preview Available
    A Quartet of PIF bHLH Factors Provides a Transcriptionally Centered Signaling Hub That Regulates Seedling Morphogenesis through Differential Expression-Patterning of Shared Target Genes in Arabidopsis
    Zhang, Y ; Mayba, O ; Pfeiffer, A ; Shi, H ; Tepperman, JM ; Speed, TP ; Quail, PH ; Copenhaver, GP (PUBLIC LIBRARY SCIENCE, 2013-01)
    Dark-grown seedlings exhibit skotomorphogenic development. Genetic and molecular evidence indicates that a quartet of Arabidopsis Phytochrome (phy)-Interacting bHLH Factors (PIF1, 3, 4, and 5) are critically necessary to maintaining this developmental state and that light activation of phy induces a switch to photomorphogenic development by inducing rapid degradation of the PIFs. Here, using integrated ChIP-seq and RNA-seq analyses, we have identified genes that are direct targets of PIF3 transcriptional regulation, exerted by sequence-specific binding to G-box (CACGTG) or PBE-box (CACATG) motifs in the target promoters genome-wide. In addition, expression analysis of selected genes in this set, in all triple pif-mutant combinations, provides evidence that the PIF quartet members collaborate to generate an expression pattern that is the product of a mosaic of differential transcriptional responsiveness of individual genes to the different PIFs and of differential regulatory activity of individual PIFs toward the different genes. Together with prior evidence that all four PIFs can bind to G-boxes, the data suggest that this collective activity may be exerted via shared occupancy of binding sites in target promoters.
  • Item
    No Preview Available
    International network of cancer genome projects
    Hudson, TJ ; Anderson, W ; Aretz, A ; Barker, AD ; Bell, C ; Bernabe, RR ; Bhan, MK ; Calvo, F ; Eerola, I ; Gerhard, DS ; Guttmacher, A ; Guyer, M ; Hemsley, FM ; Jennings, JL ; Kerr, D ; Klatt, P ; Kolar, P ; Kusuda, J ; Lane, DP ; Laplace, F ; Lu, Y ; Nettekoven, G ; Ozenberger, B ; Peterson, J ; Rao, TS ; Remacle, J ; Schafer, AJ ; Shibata, T ; Stratton, MR ; Vockley, JG ; Watanabe, K ; Yang, H ; Yuen, MMF ; Knoppers, M ; Bobrow, M ; Cambon-Thomsen, A ; Dressler, LG ; Dyke, SOM ; Joly, Y ; Kato, K ; Kennedy, KL ; Nicolas, P ; Parker, MJ ; Rial-Sebbag, E ; Romeo-Casabona, CM ; Shaw, KM ; Wallace, S ; Wiesner, GL ; Zeps, N ; Lichter, P ; Biankin, AV ; Chabannon, C ; Chin, L ; Clement, B ; de Alava, E ; Degos, F ; Ferguson, ML ; Geary, P ; Hayes, DN ; Johns, AL ; Nakagawa, H ; Penny, R ; Piris, MA ; Sarin, R ; Scarpa, A ; van de Vijver, M ; Futreal, PA ; Aburatani, H ; Bayes, M ; Bowtell, DDL ; Campbell, PJ ; Estivill, X ; Grimmond, SM ; Gut, I ; Hirst, M ; Lopez-Otin, C ; Majumder, P ; Marra, M ; Ning, Z ; Puente, XS ; Ruan, Y ; Stunnenberg, HG ; Swerdlow, H ; Velculescu, VE ; Wilson, RK ; Xue, HH ; Yang, L ; Spellman, PT ; Bader, GD ; Boutros, PC ; Flicek, P ; Getz, G ; Guigo, R ; Guo, G ; Haussler, D ; Heath, S ; Hubbard, TJ ; Jiang, T ; Jones, SM ; Li, Q ; Lopez-Bigas, N ; Luo, R ; Pearson, JV ; Quesada, V ; Raphael, BJ ; Sander, C ; Speed, TP ; Stuart, JM ; Teague, JW ; Totoki, Y ; Tsunoda, T ; Valencia, A ; Wheeler, DA ; Wu, H ; Zhao, S ; Zhou, G ; Stein, LD ; Lathrop, M ; Ouellette, BFF ; Thomas, G ; Yoshida, T ; Axton, M ; Gunter, C ; McPherson, JD ; Miller, LJ ; Kasprzyk, A ; Zhang, J ; Haider, SA ; Wang, J ; Yung, CK ; Cross, A ; Liang, Y ; Gnaneshan, S ; Guberman, J ; Hsu, J ; Chalmers, DRC ; Hasel, KW ; Kaan, TSH ; Knoppers, BM ; Lowrance, WW ; Masui, T ; Rodriguez, LL ; Vergely, C ; Cloonan, N ; Defazio, A ; Eshleman, JR ; Etemadmoghadam, D ; Gardiner, BA ; Kench, JG ; Sutherland, RL ; Tempero, MA ; Waddell, NJ ; Wilson, PJ ; Gallinger, S ; Tsao, M-S ; Shaw, PA ; Petersen, GM ; Mukhopadhyay, D ; DePinho, RA ; Thayer, S ; Muthuswamy, L ; Shazand, K ; Beck, T ; Sam, M ; Timms, L ; Ballin, V ; Ji, J ; Zhang, X ; Chen, F ; Hu, X ; Yang, Q ; Tian, G ; Zhang, L ; Xing, X ; Li, X ; Zhu, Z ; Yu, Y ; Yu, J ; Tost, J ; Brennan, P ; Holcatova, I ; Zaridze, D ; Brazma, A ; Egevad, L ; Prokhortchouk, E ; Banks, RE ; Uhlen, M ; Viksna, J ; Ponten, F ; Skryabin, K ; Birney, E ; Borg, A ; Borresen-Dale, A-L ; Caldas, C ; Foekens, JA ; Martin, S ; Reis-Filho, JS ; Richardson, AL ; Sotiriou, C ; van't Veer, L ; Birnbaum, D ; Blanche, H ; Boucher, P ; Boyault, S ; Masson-Jacquemier, JD ; Pauporte, I ; Pivot, X ; Vincent-Salomon, A ; Tabone, E ; Theillet, C ; Treilleux, I ; Bioulac-Sage, P ; Decaens, T ; Franco, D ; Gut, M ; Samuel, D ; Zucman-Rossi, J ; Eils, R ; Brors, B ; Korbel, JO ; Korshunov, A ; Landgraf, P ; Lehrach, H ; Pfister, S ; Radlwimmer, B ; Reifenberger, G ; Taylor, MD ; von Kalle, C ; Majumder, PP ; Pederzoli, P ; Lawlor, RT ; Delledonne, M ; Bardelli, A ; Gress, T ; Klimstra, D ; Zamboni, G ; Nakamura, Y ; Miyano, S ; Fujimoto, A ; Campo, E ; de Sanjose, S ; Montserrat, E ; Gonzalez-Diaz, M ; Jares, P ; Himmelbaue, H ; Bea, S ; Aparicio, S ; Easton, DF ; Collins, FS ; Compton, CC ; Lander, ES ; Burke, W ; Green, AR ; Hamilton, SR ; Kallioniemi, OP ; Ley, TJ ; Liu, ET ; Wainwright, BJ (NATURE PORTFOLIO, 2010-04-15)
    The International Cancer Genome Consortium (ICGC) was launched to coordinate large-scale cancer genome studies in tumours from 50 different cancer types and/or subtypes that are of clinical and societal importance across the globe. Systematic studies of more than 25,000 cancer genomes at the genomic, epigenomic and transcriptomic levels will reveal the repertoire of oncogenic mutations, uncover traces of the mutagenic influences, define clinically relevant subtypes for prognosis and therapeutic management, and enable the development of new cancer therapies.
  • Item
    Thumbnail Image
    Investigating and Correcting Plasma DNA Sequencing Coverage Bias to Enhance Aneuploidy Discovery
    Chandrananda, D ; Thorne, NP ; Ganesamoorthy, D ; Bruno, DL ; Benjamini, Y ; Speed, TP ; Slater, HR ; Bahlo, M ; Zhou, F (PUBLIC LIBRARY SCIENCE, 2014-01-29)
    Pregnant women carry a mixture of cell-free DNA fragments from self and fetus (non-self) in their circulation. In recent years multiple independent studies have demonstrated the ability to detect fetal trisomies such as trisomy 21, the cause of Down syndrome, by Next-Generation Sequencing of maternal plasma. The current clinical tests based on this approach show very high sensitivity and specificity, although as yet they have not become the standard diagnostic test. Here we describe improvements to the analysis of the sequencing data by reducing GC bias and better handling of the genomic repeats. We show substantial improvements in the sensitivity of the standard trisomy 21 statistical tests, which we measure by artificially reducing read coverage. We also explore the bias stemming from the natural cleavage of plasma DNA by examining DNA motifs and position specific base distributions. We propose a model to correct this fragmentation bias and observe that incorporating this bias does not lead to any further improvements in the detection of fetal trisomy. The improved bias corrections that we demonstrate in this work can be readily adopted into existing fetal trisomy detection protocols and should also lead to improvements in sub-chromosomal copy number variation detection.
  • Item
    Thumbnail Image
    Lineage-specific expansion of proteins exported to erythrocytes in malaria parasites
    Sargeant, TJ ; Marti, M ; Caler, E ; Carlton, JM ; Simpson, K ; Speed, TP ; Cowman, AF (BMC, 2006)
    BACKGROUND: The apicomplexan parasite Plasmodium falciparum causes the most severe form of malaria in humans. After invasion into erythrocytes, asexual parasite stages drastically alter their host cell and export remodeling and virulence proteins. Previously, we have reported identification and functional analysis of a short motif necessary for export of proteins out of the parasite and into the red blood cell. RESULTS: We have developed software for the prediction of exported proteins in the genus Plasmodium, and identified exported proteins conserved between malaria parasites infecting rodents and the two major causes of human malaria, P. falciparum and P. vivax. This conserved 'exportome' is confined to a few subtelomeric chromosomal regions in P. falciparum and the synteny of these and surrounding regions is conserved in P. vivax. We have identified a novel gene family PHIST (for Plasmodium helical interspersed subtelomeric family) that shares a unique domain with 72 paralogs in P. falciparum and 39 in P. vivax; however, there is only one member in each of the three species studied from the P. berghei lineage. CONCLUSION: These data suggest radiation of genes encoding remodeling and virulence factors from a small number of loci in a common Plasmodium ancestor, and imply a closer phylogenetic relationship between the P. vivax and P. falciparum lineages than previously believed. The presence of a conserved 'exportome' in the genus Plasmodium has important implications for our understanding of both common mechanisms and species-specific differences in host-parasite interactions, and may be crucial in developing novel antimalarial drugs to this infectious disease.
  • Item
    Thumbnail Image
    Integrative analysis of RUNX1 downstream pathways and target genes
    Michaud, J ; Simpson, KM ; Escher, R ; Buchet-Poyau, K ; Beissbarth, T ; Carmichael, C ; Ritchie, ME ; Schuetz, F ; Cannon, P ; Liu, M ; Shen, X ; Ito, Y ; Raskind, WH ; Horwitz, MS ; Osato, M ; Turner, DR ; Speed, TP ; Kavallaris, M ; Smyth, GK ; Scott, HS (BMC, 2008-07-31)
    BACKGROUND: The RUNX1 transcription factor gene is frequently mutated in sporadic myeloid and lymphoid leukemia through translocation, point mutation or amplification. It is also responsible for a familial platelet disorder with predisposition to acute myeloid leukemia (FPD-AML). The disruption of the largely unknown biological pathways controlled by RUNX1 is likely to be responsible for the development of leukemia. We have used multiple microarray platforms and bioinformatic techniques to help identify these biological pathways to aid in the understanding of why RUNX1 mutations lead to leukemia. RESULTS: Here we report genes regulated either directly or indirectly by RUNX1 based on the study of gene expression profiles generated from 3 different human and mouse platforms. The platforms used were global gene expression profiling of: 1) cell lines with RUNX1 mutations from FPD-AML patients, 2) over-expression of RUNX1 and CBFbeta, and 3) Runx1 knockout mouse embryos using either cDNA or Affymetrix microarrays. We observe that our datasets (lists of differentially expressed genes) significantly correlate with published microarray data from sporadic AML patients with mutations in either RUNX1 or its cofactor, CBFbeta. A number of biological processes were identified among the differentially expressed genes and functional assays suggest that heterozygous RUNX1 point mutations in patients with FPD-AML impair cell proliferation, microtubule dynamics and possibly genetic stability. In addition, analysis of the regulatory regions of the differentially expressed genes has for the first time systematically identified numerous potential novel RUNX1 target genes. CONCLUSION: This work is the first large-scale study attempting to identify the genetic networks regulated by RUNX1, a master regulator in the development of the hematopoietic system and leukemia. The biological pathways and target genes controlled by RUNX1 will have considerable importance in disease progression in both familial and sporadic leukemia as well as therapeutic implications.
  • Item
    Thumbnail Image
    Copy Number Variation in Patients with Disorders of Sex Development Due to 46,XY Gonadal Dysgenesis
    White, S ; Ohnesorg, T ; Notini, A ; Roeszler, K ; Hewitt, J ; Daggag, H ; Smith, C ; Turbitt, E ; Gustin, S ; van den Bergen, J ; Miles, D ; Western, P ; Arboleda, V ; Schumacher, V ; Gordon, L ; Bell, K ; Bengtsson, H ; Speed, T ; Hutson, J ; Warne, G ; Harley, V ; Koopman, P ; Vilain, E ; Sinclair, A ; Orban, L (PUBLIC LIBRARY SCIENCE, 2011-03-07)
    Disorders of sex development (DSD), ranging in severity from mild genital abnormalities to complete sex reversal, represent a major concern for patients and their families. DSD are often due to disruption of the genetic programs that regulate gonad development. Although some genes have been identified in these developmental pathways, the causative mutations have not been identified in more than 50% 46,XY DSD cases. We used the Affymetrix Genome-Wide Human SNP Array 6.0 to analyse copy number variation in 23 individuals with unexplained 46,XY DSD due to gonadal dysgenesis (GD). Here we describe three discrete changes in copy number that are the likely cause of the GD. Firstly, we identified a large duplication on the X chromosome that included DAX1 (NR0B1). Secondly, we identified a rearrangement that appears to affect a novel gonad-specific regulatory region in a known testis gene, SOX9. Surprisingly this patient lacked any signs of campomelic dysplasia, suggesting that the deletion affected expression of SOX9 only in the gonad. Functional analysis of potential SRY binding sites within this deleted region identified five putative enhancers, suggesting that sequences additional to the known SRY-binding TES enhancer influence human testis-specific SOX9 expression. Thirdly, we identified a small deletion immediately downstream of GATA4, supporting a role for GATA4 in gonad development in humans. These CNV analyses give new insights into the pathways involved in human gonad development and dysfunction, and suggest that rearrangements of non-coding sequences disturbing gene regulation may account for significant proportion of DSD cases.
  • Item
    Thumbnail Image
    Differential splicing using whole-transcript microarrays
    Robinson, MD ; Speed, TP (BMC, 2009-05-22)
    BACKGROUND: The latest generation of Affymetrix microarrays are designed to interrogate expression over the entire length of every locus, thus giving the opportunity to study alternative splicing genome-wide. The Exon 1.0 ST (sense target) platform, with versions for Human, Mouse and Rat, is designed primarily to probe every known or predicted exon. The smaller Gene 1.0 ST array is designed as an expression microarray but still interrogates expression with probes along the full length of each well-characterized transcript. We explore the possibility of using the Gene 1.0 ST platform to identify differential splicing events. RESULTS: We propose a strategy to score differential splicing by using the auxiliary information from fitting the statistical model, RMA (robust multichip analysis). RMA partitions the probe-level data into probe effects and expression levels, operating robustly so that if a small number of probes behave differently than the rest, they are downweighted in the fitting step. We argue that adjacent poorly fitting probes for a given sample can be evidence of differential splicing and have designed a statistic to search for this behaviour. Using a public tissue panel dataset, we show many examples of tissue-specific alternative splicing. Furthermore, we show that evidence for putative alternative splicing has a strong correspondence between the Gene 1.0 ST and Exon 1.0 ST platforms. CONCLUSION: We propose a new approach, FIRMAGene, to search for differentially spliced genes using the Gene 1.0 ST platform. Such an analysis complements the search for differential expression. We validate the method by illustrating several known examples and we note some of the challenges in interpreting the probe-level data.Software implementing our methods is freely available as an R package.
  • Item
    Thumbnail Image
    A High Force of Plasmodium vivax Blood-Stage Infection Drives the Rapid Acquisition of Immunity in Papua New Guinean Children
    Koepfli, C ; Colborn, KL ; Kiniboro, B ; Lin, E ; Speed, TP ; Siba, PM ; Felger, I ; Mueller, I ; McCarthy, JS (PUBLIC LIBRARY SCIENCE, 2013-09)
    BACKGROUND: When both parasite species are co-endemic, Plasmodium vivax incidence peaks in younger children compared to P. falciparum. To identify differences in the number of blood stage infections of these species and its potential link to acquisition of immunity, we have estimated the molecular force of blood-stage infection of P. vivax ((mol)FOB, i.e. the number of genetically distinct blood-stage infections over time), and compared it to previously reported values for P. falciparum. METHODS: P. vivax (mol)FOB was estimated by high resolution genotyping parasites in samples collected over 16 months in a cohort of 264 Papua New Guinean children living in an area highly endemic for P. falciparum and P. vivax. In this cohort, P. vivax episodes decreased three-fold over the age range of 1-4.5 years. RESULTS: On average, children acquired 14.0 new P. vivax blood-stage clones/child/year-at-risk. While the incidence of clinical P. vivax illness was strongly associated with mol FOB (incidence rate ratio (IRR) = 1.99, 95% confidence interval (CI95) [1.80, 2.19]), (mol)FOB did not change with age. The incidence of P. vivax showed a faster decrease with age in children with high (IRR = 0.49, CI95 [0.38, 0.64] p<0.001) compared to those with low exposure (IRR = 0.63, CI95[0.43, 0.93] p = 0.02). CONCLUSION: P. vivax (mol)FOB is considerably higher than P. falciparum (mol)FOB (5.5 clones/child/year-at-risk). The high number of P. vivax clones that infect children in early childhood contribute to the rapid acquisition of immunity against clinical P. vivax malaria.