AGBT-2018-Andy-Pang.Pdf

Total Page:16

File Type:pdf, Size:1020Kb

AGBT-2018-Andy-Pang.Pdf Efficient Structural Variation Detection and Annotation Using Bionano Genome Mapping A. W. C. Pang1, J. Lee1, A. Hastie1, E. Lam1, E. Chan2, V. Hayes2, H. Cao1, M. Borodkin1 1) Bionano Genomics, San Diego, California, USA; 2) Garvan Institute of Medical Research, Darlinghurst, Australia Abstract Background Structural variation (SV) detection is fundamental to understanding cancer workflow can compare against a control sample, and examines whether the Generating high-quality finished genomes replete with accurate genomes. While karyotyping and conventional molecular detection approaches cancer mutations are present in low fraction among the control’s molecules. identification of structural variation and high completion (minimal are robust, they can be manually intensive, biased towards targeted loci, and Using this pipeline, whose runtime is only a few hours, we can efficiently focus gaps) remains challenging using short read sequencing cannot determine the copy number of long repeats. on several dozen significant candidates for further analysis. technologies alone. Bionano mapping provides direct visualization of long DNA molecules in their native state, bypassing the statistical inference needed to align Bionano Genomics’ Saphyr System offers a sensitive method for detecting large We ran multiple solid and hematologic cancer samples. First, we generated a paired-end reads with an uncertain insert size distribution. These long labeled SVs. DNA molecules larger than 100 kbp are extracted, labelled at specific highly contiguous genome map assembly on a solid tumor sample. The ultra molecules are de novo assembled into physical maps spanning the whole genome. motifs, and linearized through NanoChannels arrays for subsequent long maps – with lengths encompassing entire chromosomal arms – were able The resulting order and orientation of sequence elements in the map can be used visualization. Molecule images are digitized and the genome de novo to capture large inversion events, and balanced and unbalanced translocations. for anchoring NGS contigs and structural variation detection. assembled, creating megabase long Bionano maps. While our calling algorithm In a separate study, we discovered known rearrangements such as the can sensitively detect all types of SVs in these long maps, we also use a variant t(9;22)(q34.12;q11.23) in chronic myeloid leukemia (CML). Moreover, we annotation workflow to specifically uncover rare and sample-specific mutations. uncovered a large megabases long deletion on chr13 by using our copy number To determine variant frequency in a tumor normal experimental design, it detection tool. In conclusion, Bionano maps can discover functional SVs and compares the cancer sample’s calls against over 600,000 SVs from >160 improve our understanding of the mechanisms of diseases. humans with no observable diseases. To identify somatic mutations, the Methods 1 2 3 4 5 6 Label DNA at specific sequence Saphyr Chip linearizes DNA in Saphyr automates imaging of single Molecules and labels detected in Bionano Access software Extraction of long DNA molecules motifs NanoChannel arrays molecules in NanoChannel arrays images by instrument software assembles optical maps Free DNA Solution DNA in a Microchannel DNA in a Nanochannel Blood Cell Tissue Microbes Gaussian Coil Partially Elongated Linearized Position (kbp) (1) Long molecules of DNA are labeled with Bionano reagents by (2) incorporation of fluorophores at a specific sequence motif throughout the genome. (3) The labeled genomic DNA is then linearized in the Saphyr Chip using NanoChannel arrays (4) Single molecules are imaged by Saphyr and then digitized. (5) Molecules are uniquely identifiable by distinct distribution of sequence motif labels (6) and then assembled by pairwise alignment into de novo genome maps. Study 1: Direct labeling experiment of a solid tumor sample Study 2: SV results on samples with hematological malignancy Cancer SV DLE1 labeled sample Solid tumor 135.5 Mbp chr10 Data collected (molecules >150 kb) 367 Gbp Step 1 Tumor • Compare with > 160 human control samples SV from genome mapping maps 38.8 Mbp 91.7 Mbp • Examine the coverage and chimeric quality scores Assembly size (haplotype-aware) 5.87 Gbp • Overlap with gene annotation 78.1 Mbp 133,653,967 Genome map N50 71.0 Mbp chr18 hg19 chr9 ABL1 Step 2 Tumor • Compare the cancer sample calls with a match control sample calls SV called against hg19 maps 15.6 Mbp 59.5 Mbp Map 23,631,928 Insertions 3,894 Two examples of ultra long genome maps spanning entire arms of chromosomes 10 and 18. The maps hg19 chr22 BCR are only separated by centromeres. Step 3 • Align the match control sample molecules to the cancer assembly Deletions 1,799 • Re-align the cancer molecules to the cancer assembly • Determine the number of molecules supporting the cancer variant call Inversion breakpoints 95 er. chromosomal er. fragments RBM6 d chr22 chr9 Genes Annotated cancer Interchromosomal translocations 33 SV A Philadelphia translocation t(9;22) was detected in leukemia cancer cells. The A flowchart illustrating the variant annotation pipeline for chr3 Intrachromosomal translocations 21 map was aligned to the public reference assembly hg19, and the resulting alignments cancer-normal studies. show a conjoined junction between chr9 and chr22, creating a fusion gene called BCR- PKD1P1 XYLT1 NPIP Tumor ABL1. Genes map chr17:26.6 M chr16 chrX 50.6 M chr17 fusion BspQI map Copy number chr3 Tumor map Tumor map Molecule pileup Molecule chrX 1 3 5 7 9 11 13 15 17 19 21 X pileup A whole genome depth of coverage shoring a large chr13 deletion. The green dots represent raw copy number count, while the blue dots represent smoothed copy (Top) A 596.9 kbp genome map captures a t(3;X) translocation. The translocation breakpoint interrupts number count after the removal of local noise the gene RBM6 on chr3. RBM6 codes for a RNA-binding protein 6, and mutations and copy number mutations have been found in many types of cancer. A 2.3 Mbp inversion on chr16 captures by a 32.2 Mbp de novo assembled tumor genome map. Note (Bottom) An alternate map is able to capture the reciprocal translocation event. that the breakpoints contain some additional sequences that are novel to the reference genome hg19. The An intra-chromosomal fusion event detected in chr17. The two genes presence of overlapping alignments indicates homology (in opposite direction) between the two inversion located at the junctions are TAOK1 and HOXB3. The red arrows indicate the breakpoints. variant junction on the genome map. Conclusions Reference 1) Cao, H., et al., Rapid detection of structural variation in a human genome using NanoChannel- based genome mapping technology. Gigascience (2014); 3(1):34 We demonstrate that the Saphyr system can be used to accurately detect genetic mutation hallmarks in samples with solid tumor and 2) Hastie, A.R., et al. Rapid Genome Mapping in NanoChannel Arrays for Highly Complete and hematologic malignancies. These includes large rearrangements ranging from interchromosomal translocations, within chromosome Accurate De Novo Sequence Assembly of the Complex Aegilops tauschii Genome. PLoS ONE fusions, to copy number alterations. Researchers can perform mapping experiments to uncover somatic variants by comparing with (2013); 8(2): e55864. our control sample database and a matching non-tumor sample. Our results shown here indicate that the Saphyr system can capture a 3) Lam, E.T., et al. Genome mapping on NanoChannel arrays for structural variation analysis and broad spectrum of variation with functional importance, and can provide easy solutions for cancer studies. sequence assembly. Nature Biotechnology (2012); 10: 2303 4) Das, S. K., et al. Single molecule linear analysis of DNA in nano-channel labeled with sequence specific fluorescent probes. Nucleic Acids Research (2010); 38: 8 5) Xiao, M et. al. Rapid DNA mapping by fluorescent single molecule detection. Nucleic Acids ©2018 Bionano Genomics. All ©2018 All rights reserved. Bionano Genomics. Research (2007); 35:e16..
Recommended publications
  • Whole-Genome Microarray Detects Deletions and Loss of Heterozygosity of Chromosome 3 Occurring Exclusively in Metastasizing Uveal Melanoma
    Anatomy and Pathology Whole-Genome Microarray Detects Deletions and Loss of Heterozygosity of Chromosome 3 Occurring Exclusively in Metastasizing Uveal Melanoma Sarah L. Lake,1 Sarah E. Coupland,1 Azzam F. G. Taktak,2 and Bertil E. Damato3 PURPOSE. To detect deletions and loss of heterozygosity of disease is fatal in 92% of patients within 2 years of diagnosis. chromosome 3 in a rare subset of fatal, disomy 3 uveal mela- Clinical and histopathologic risk factors for UM metastasis noma (UM), undetectable by fluorescence in situ hybridization include large basal tumor diameter (LBD), ciliary body involve- (FISH). ment, epithelioid cytomorphology, extracellular matrix peri- ϩ ETHODS odic acid-Schiff-positive (PAS ) loops, and high mitotic M . Multiplex ligation-dependent probe amplification 3,4 5 (MLPA) with the P027 UM assay was performed on formalin- count. Prescher et al. showed that a nonrandom genetic fixed, paraffin-embedded (FFPE) whole tumor sections from 19 change, monosomy 3, correlates strongly with metastatic death, and the correlation has since been confirmed by several disomy 3 metastasizing UMs. Whole-genome microarray analy- 3,6–10 ses using a single-nucleotide polymorphism microarray (aSNP) groups. Consequently, fluorescence in situ hybridization were performed on frozen tissue samples from four fatal dis- (FISH) detection of chromosome 3 using a centromeric probe omy 3 metastasizing UMs and three disomy 3 tumors with Ͼ5 became routine practice for UM prognostication; however, 5% years’ metastasis-free survival. to 20% of disomy 3 UM patients unexpectedly develop metas- tases.11 Attempts have therefore been made to identify the RESULTS. Two metastasizing UMs that had been classified as minimal region(s) of deletion on chromosome 3.12–15 Despite disomy 3 by FISH analysis of a small tumor sample were found these studies, little progress has been made in defining the key on MLPA analysis to show monosomy 3.
    [Show full text]
  • Methods Mouse Strains and Housing Male Wild-Type Mice
    Methods Mouse strains and housing Male wild-type mice (C57BL/6J background) and B6.Cg-Tg (APPSwFlLon, PSEN1*M146L*L286V) 6799Vas/Mmjax (5xFAD, JAX 008730) were purchased from the Jackson Laboratory. μMT-/- mice which are deficient in B cells were purchased from Shanghai Model Organisms Center, Inc. Both B6.Cg-Tg and μMT-/- mice were bred on a C57BL/6J background. 5xFAD mice and μMT-/- mice were crossed to generate μMT-/-/5xFAD mice. The mice were used at different ages that are indicated throughout the manuscript. Mice of all strains were specific pathogen free environment with controlled temperature and humidity, on 12 h light: dark cycles (lights on at 7:00), and fed with regular rodent’s chow and sterilized tap water ad libitum. All experiments were approved by the Institutional Animal Care and Use Committee of the Nanjing Medical University. Human samples Frontal cortex tissues were obtained from 4 cases within 4-6 h of death, via informed donation for the Medical Education and Research of Nanjing Medical University, with corresponding written consents prepared by the donors and their families. Cases that were died from brain associated diseases were excluded from this study. The utilization of human tissues was approved by the Ethics Committee of Nanjing Medical University. All obtained samples were fixed in a 4% formalin solution and kept in paraffin blocks until further sectioning. Animal surgery DcLN ligation: The procedure of surgical ligation of the lymphatics afferent to the dcLNs was according to published literature (3, 30). In brief, mice were anaesthetized by i.p. injection with ketamine and xylazine in saline, and fixed on a stereotaxic apparatus in a supine position.
    [Show full text]
  • Supplemental Figure and Table Legends
    Supplemental figure and table legends Supplementary Figure 1: KIAA1841 is well conserved among vertebrates. NCBI HomoloGene pairwise alignment scores of human KIAA1841 sequence compared to other vertebrate orthologs. Supplementary Figure 2: µ-germline transcripts (GLT) and AID mRNA expression are not affected by overexpression of KIAA1841. Splenic B cells were isolated from wild-type mice, and transduced with retroviral vector control (pMIG) or a vector expressing KIAA1841. Levels of µ-GLT and AID mRNA were determined at 72h post-infection by RT-qPCR, and normalized to -actin mRNA and the pMIG control. The mean of three independent experiments +/- SD is shown. NS, p = not significant, p 0.05, two-tailed paired student’s t-test. Supplementary Figure 3: Overexpression of untagged and Xpress-tagged KIAA1841 does not affect cell proliferation. Splenic B cells were isolated from wild-type mice, stimulated with LPS+IL4, and transduced with retroviral vector control (pMIG) or a vector expressing KIAA1841 or Xpress (Xp)-tagged KIAA1841. Cells are labeled with seminaphthorhodafluor (SNARF) cell tracking dye and SNARF intensity was measured at 0h, 24h, and 48h after retroviral infection. Histograms of transduced cells (GFP+) for pMIG control, KIAA1841 and Xp-KIAA1841 were superimposed at each time point. Three independent retroviral infection experiments are shown. Supplementary Figure 4: Sequence alignment of the putative SANT domain of KIAA1841 with the SANT domain of SWI3. Alignment was performed using ClustalOmega; *, conserved residue, :, strongly similar residues, ., weakly similar residues. Numbers indicate amino acid residues in each sequence. Helix 3, which has been reported to be important for the chromatin remodeling function of SWI3 (Boyer et.
    [Show full text]
  • Transcriptional Recapitulation and Subversion Of
    Open Access Research2007KaiseretVolume al. 8, Issue 7, Article R131 Transcriptional recapitulation and subversion of embryonic colon comment development by mouse colon tumor models and human colon cancer Sergio Kaiser¤*, Young-Kyu Park¤†, Jeffrey L Franklin†, Richard B Halberg‡, Ming Yu§, Walter J Jessen*, Johannes Freudenberg*, Xiaodi Chen‡, Kevin Haigis¶, Anil G Jegga*, Sue Kong*, Bhuvaneswari Sakthivel*, Huan Xu*, Timothy Reichling¥, Mohammad Azhar#, Gregory P Boivin**, reviews Reade B Roberts§, Anika C Bissahoyo§, Fausto Gonzales††, Greg C Bloom††, Steven Eschrich††, Scott L Carter‡‡, Jeremy E Aronow*, John Kleimeyer*, Michael Kleimeyer*, Vivek Ramaswamy*, Stephen H Settle†, Braden Boone†, Shawn Levy†, Jonathan M Graff§§, Thomas Doetschman#, Joanna Groden¥, William F Dove‡, David W Threadgill§, Timothy J Yeatman††, reports Robert J Coffey Jr† and Bruce J Aronow* Addresses: *Biomedical Informatics, Cincinnati Children's Hospital Medical Center, Cincinnati, OH 45229, USA. †Departments of Medicine, and Cell and Developmental Biology, Vanderbilt University and Department of Veterans Affairs Medical Center, Nashville, TN 37232, USA. ‡McArdle Laboratory for Cancer Research, University of Wisconsin, Madison, WI 53706, USA. §Department of Genetics and Lineberger Cancer Center, University of North Carolina, Chapel Hill, NC 27599, USA. ¶Molecular Pathology Unit and Center for Cancer Research, Massachusetts deposited research General Hospital, Charlestown, MA 02129, USA. ¥Division of Human Cancer Genetics, The Ohio State University College of Medicine, Columbus, Ohio 43210-2207, USA. #Institute for Collaborative BioResearch, University of Arizona, Tucson, AZ 85721-0036, USA. **University of Cincinnati, Department of Pathology and Laboratory Medicine, Cincinnati, OH 45267, USA. ††H Lee Moffitt Cancer Center and Research Institute, Tampa, FL 33612, USA. ‡‡Children's Hospital Informatics Program at the Harvard-MIT Division of Health Sciences and Technology (CHIP@HST), Harvard Medical School, Boston, Massachusetts 02115, USA.
    [Show full text]
  • Use of IBM Watson to Identify Additional RNA-Binding Proteins
    Acta Neuropathol DOI 10.1007/s00401-017-1785-8 ORIGINAL PAPER Artifcial intelligence in neurodegenerative disease research: use of IBM Watson to identify additional RNA‑binding proteins altered in amyotrophic lateral sclerosis Nadine Bakkar1 · Tina Kovalik1 · Ileana Lorenzini1 · Scott Spangler2 · Alix Lacoste3 · Kyle Sponaugle1 · Philip Ferrante1 · Elenee Argentinis3 · Rita Sattler1 · Robert Bowser1 Received: 29 September 2017 / Revised: 4 November 2017 / Accepted: 4 November 2017 © The Author(s) 2017. This article is an open access publication Abstract Amyotrophic lateral sclerosis (ALS) is a dev- to controls. Overall, we successfully used IBM Watson to astating neurodegenerative disease with no efective treat- help identify additional RBPs altered in ALS, highlighting ments. Numerous RNA-binding proteins (RBPs) have been the use of artifcial intelligence tools to accelerate scientifc shown to be altered in ALS, with mutations in 11 RBPs discovery in ALS and possibly other complex neurological causing familial forms of the disease, and 6 more RBPs disorders. showing abnormal expression/distribution in ALS albeit without any known mutations. RBP dysregulation is widely Keywords Amyotrophic lateral sclerosis · RNA-binding accepted as a contributing factor in ALS pathobiology. There protein · Artifcial intelligence · Protein aggregation · are at least 1542 RBPs in the human genome; therefore, Motor neuron other unidentifed RBPs may also be linked to the patho- genesis of ALS. We used IBM ­Watson® to sieve through all RBPs in the genome and identify new RBPs linked to ALS Introduction (ALS-RBPs). IBM Watson extracted features from pub- lished literature to create semantic similarities and identify Amyotrophic lateral sclerosis (ALS) is characterized by new connections between entities of interest.
    [Show full text]
  • Mouse Rbm6 Knockout Project (CRISPR/Cas9)
    http://www.alphaknockout.com/ Mouse Rbm6 Knockout Project (CRISPR/Cas9) Objective: To create a Rbm6 knockout Mouse model (C57BL/6N) by CRISPR/Cas-mediated genome engineering. Strategy summary: The Rbm6 gene (NCBI Reference Sequence: NM_029169 ; Ensembl: ENSMUSG00000032582 ) is located on Mouse chromosome 9. 21 exons are identified, with the ATG start codon in exon 3 and the TAA stop codon in exon 21 (Transcript: ENSMUST00000035201). Exon 3~17 will be selected as target site. Cas9 and gRNA will be co-injected into fertilized eggs for KO Mouse production. The pups will be genotyped by PCR followed by sequencing analysis. Note: Exon 3 starts from about 0.1% of the coding region. Exon 3~17 covers 85.6% of the coding region. The size of effective KO region: ~71065 bp. The KO region does not have any other known gene. Page 1 of 9 http://www.alphaknockout.com/ Overview of the Targeting Strategy Wildtype allele 5' gRNA region 15 3' 13 gRNA region 10 12 17 1 3 4 5 6 7 8 9 11 14 16 21 Legends Exon of mouse Rbm6 Knockout region Page 2 of 9 http://www.alphaknockout.com/ Overview of the Dot Plot (up) Window size: 15 bp Forward Reverse Complement Sequence 12 Note: The 2000 bp section of Exon 3 is aligned with itself to determine if there are tandem repeats. No significant tandem repeat is found in the dot plot matrix. So this region is suitable for PCR screening or sequencing analysis. Overview of the Dot Plot (down) Window size: 15 bp Forward Reverse Complement Sequence 12 Note: The 2000 bp section of Exon 17 is aligned with itself to determine if there are tandem repeats.
    [Show full text]
  • Interplay of RNA-Binding Proteins and Micrornas in Neurodegenerative Diseases
    International Journal of Molecular Sciences Review Interplay of RNA-Binding Proteins and microRNAs in Neurodegenerative Diseases Chisato Kinoshita 1,* , Noriko Kubota 1,2 and Koji Aoyama 1,* 1 Department of Pharmacology, Teikyo University School of Medicine, 2-11-1 Kaga, Itabashi, Tokyo 173-8605, Japan; [email protected] 2 Teikyo University Support Center for Women Physicians and Researchers, 2-11-1 Kaga, Itabashi, Tokyo 173-8605, Japan * Correspondence: [email protected] (C.K.); [email protected] (K.A.); Tel.: +81-3-3964-3794 (C.K.); +81-3-3964-3793 (K.A.) Abstract: The number of patients with neurodegenerative diseases (NDs) is increasing, along with the growing number of older adults. This escalation threatens to create a medical and social crisis. NDs include a large spectrum of heterogeneous and multifactorial pathologies, such as amyotrophic lateral sclerosis, frontotemporal dementia, Alzheimer’s disease, Parkinson’s disease, Huntington’s disease and multiple system atrophy, and the formation of inclusion bodies resulting from protein misfolding and aggregation is a hallmark of these disorders. The proteinaceous components of the pathological inclusions include several RNA-binding proteins (RBPs), which play important roles in splicing, stability, transcription and translation. In addition, RBPs were shown to play a critical role in regulating miRNA biogenesis and metabolism. The dysfunction of both RBPs and miRNAs is Citation: Kinoshita, C.; Kubota, N.; often observed in several NDs. Thus, the data about the interplay among RBPs and miRNAs and Aoyama, K. Interplay of RNA-Binding Proteins and their cooperation in brain functions would be important to know for better understanding NDs and microRNAs in Neurodegenerative the development of effective therapeutics.
    [Show full text]
  • Supporting Information
    Supporting Information Friedman et al. 10.1073/pnas.0812446106 SI Results and Discussion intronic miR genes in these protein-coding genes. Because in General Phenotype of Dicer-PCKO Mice. Dicer-PCKO mice had many many cases the exact borders of the protein-coding genes are defects in additional to inner ear defects. Many of them died unknown, we searched for miR genes up to 10 kb from the around birth, and although they were born at a similar size to hosting-gene ends. Out of the 488 mouse miR genes included in their littermate heterozygote siblings, after a few weeks the miRBase release 12.0, 192 mouse miR genes were found as surviving mutants were smaller than their heterozygote siblings located inside (distance 0) or in the vicinity of the protein-coding (see Fig. 1A) and exhibited typical defects, which enabled their genes that are expressed in the P2 cochlear and vestibular SE identification even before genotyping, including typical alopecia (Table S2). Some coding genes include huge clusters of miRNAs (in particular on the nape of the neck), partially closed eyelids (e.g., Sfmbt2). Other genes listed in Table S2 as coding genes are [supporting information (SI) Fig. S1 A and C], eye defects, and actually predicted, as their transcript was detected in cells, but weakness of the rear legs that were twisted backwards (data not the predicted encoded protein has not been identified yet, and shown). However, while all of the mutant mice tested exhibited some of them may be noncoding RNAs. Only a single protein- similar deafness and stereocilia malformation in inner ear HCs, coding gene that is differentially expressed in the cochlear and other defects were variable in their severity.
    [Show full text]
  • Skeletal Muscle Function, Morphology, and Biochemistry in Ts65dn Mice: a Model of Down Syndrome
    Syracuse University SURFACE Exercise Science - Dissertations School of Education 12-2011 Skeletal Muscle Function, Morphology, and Biochemistry in Ts65Dn Mice: A Model of Down Syndrome Patrick Michael Cowley Syracuse University Follow this and additional works at: https://surface.syr.edu/ppe_etd Part of the Kinesiology Commons Recommended Citation Cowley, Patrick Michael, "Skeletal Muscle Function, Morphology, and Biochemistry in Ts65Dn Mice: A Model of Down Syndrome" (2011). Exercise Science - Dissertations. 5. https://surface.syr.edu/ppe_etd/5 This Dissertation is brought to you for free and open access by the School of Education at SURFACE. It has been accepted for inclusion in Exercise Science - Dissertations by an authorized administrator of SURFACE. For more information, please contact [email protected]. ABSTRACT A common clinical observation of persons with Down syndrome at all developmental stages is hypotonia and generalized muscle weakness. The cause of muscle weakness in Down syndrome is not known and there is an immediate need to establish an acceptable animal model to explore the muscle dysfunction that is widely reported in the human population. Using a combination of functional, histological, and biochemical analyses this dissertation provides the initial characterization of skeletal muscle from the Ts65Dn mouse, a model of Down syndrome. The experiments revealed that Ts65Dn muscle over-expresses SOD1 protein but this did not lead to oxidative stress. Ts65Dn soleus muscles displayed normal force generation in the unfatigued state, but exhibited muscle weakness following fatiguing contractions. We show that a reduction in cytochrome c oxidase expression may contribute to the impaired muscle performance in Ts65Dn soleus. These findings support the use of the Ts65Dn mouse model of Down syndrome to delineate mechanisms of muscle dysfunction in the human condition.
    [Show full text]
  • Novel Gene Signatures for Stage Classification of the Squamous Cell
    www.nature.com/scientificreports OPEN Novel gene signatures for stage classifcation of the squamous cell carcinoma of the lung Angel Juarez‑Flores, Gabriel S. Zamudio & Marco V. José* The squamous cell carcinoma of the lung (SCLC) is one of the most common types of lung cancer. As GLOBOCAN reported in 2018, lung cancer was the frst cause of death and new cases by cancer worldwide. Typically, diagnosis is made in the later stages of the disease with few treatment options available. The goal of this work was to fnd some key components underlying each stage of the disease, to help in the classifcation of tumor samples, and to increase the available options for experimental assays and molecular targets that could be used in treatment development. We employed two approaches. The frst was based in the classic method of diferential gene expression analysis, network analysis, and a novel concept known as network gatekeepers. The second approach was using machine learning algorithms. From our combined approach, we identifed two sets of genes that could function as a signature to identify each stage of the cancer pathology. We also arrived at a network of 55 nodes, which according to their biological functions, they can be regarded as drivers in this cancer. Although biological experiments are necessary for their validation, we proposed that all these genes could be used for cancer development treatments. As GLOBOCAN reported in 2018, lung cancer was the frst cause of deaths and new cases by cancer worldwide 1. Squamous cell carcinoma of the lung (SCC) is one type of lung cancer which comprises approximately 30% of all lung cancer cases.
    [Show full text]
  • Supplemental Table S-2. Genes That Were Differentially Regulated by 2-Fold Or Greater Between Male and Female Tendons
    Supplemental Table S-2. Genes that were differentially regulated by 2-fold or greater between male and female tendons. Log Fold Affymetrix Change Male vs Probe ID Entrez ID Symbol Description Female 17543785 213742 Xist inactive X specific transcripts -5.47 17426022 100041687 Mup8 major urinary protein 8 -3.73 17426062 100048885 Gm2083 major urinary protein LOC100048885 -3.53 17426043 100048885 Gm2083 major urinary protein LOC100048885 -3.50 17297591 NA NA NA -3.50 17426011 100048885 Gm2083 major urinary protein LOC100048885 -3.46 17425990 100048885 Gm2083 major urinary protein LOC100048885 -3.46 17425915 100048885 Gm2083 major urinary protein LOC100048885 -3.46 17425941 100048885 Gm2083 major urinary protein LOC100048885 -3.40 17517360 66402 Sln sarcolipin -3.18 17426032 100048885 Gm2083 major urinary protein LOC100048885 -3.12 17346069 240120 Zfp119b zinc finger protein 119b -3.07 17426000 100048885 Gm2083 major urinary protein LOC100048885 -3.03 17425900 100048885 Gm2083 major urinary protein LOC100048885 -2.91 17425888 100048885 Gm2083 major urinary protein LOC100048885 -2.85 17315002 70152 Mettl7a1 methyltransferase like 7A1 -2.62 17248348 NA NA NA -2.41 17524209 NA NA NA -2.39 17443890 24135 Zfp68 zinc finger protein 68 -2.38 17238834 NA NA NA -2.33 17343251 108121 U2af1 U2 small nuclear ribonucleoprotein auxiliary factor (U2AF) 1 -2.32 17406421 99889 Arfip1 ADP-ribosylation factor interacting protein 1 -2.30 17497612 93747 Echs1 enoyl Coenzyme A hydratase, short chain, 1, mitochondrial -2.22 17238892 NA NA NA -2.22 17238934 64009
    [Show full text]
  • UTR Shortening Has a Limited Effect on Protein Abundance in Proliferating T Cells
    ARTICLE Received 17 Apr 2014 | Accepted 3 Oct 2014 | Published 21 Nov 2014 DOI: 10.1038/ncomms6465 Global 30 UTR shortening has a limited effect on protein abundance in proliferating T cells Andreas R. Gruber1,*, Georges Martin1,*, Philipp Mu¨ller2,3, Alexander Schmidt4, Andreas J. Gruber1, Rafal Gumienny1, Nitish Mittal1, Rajesh Jayachandran3, Jean Pieters3, Walter Keller1, Erik van Nimwegen1 & Mihaela Zavolan1 Alternative polyadenylation is a cellular mechanism that generates mRNA isoforms differing in their 30 untranslated regions (30 UTRs). Changes in polyadenylation site usage have been described upon induction of proliferation in resting cells, but the underlying mechanism and functional significance of this phenomenon remain largely unknown. To understand the functional consequences of shortened 30 UTR isoforms in a physiological setting, we used 30 end sequencing and quantitative mass spectrometry to determine polyadenylation site usage, mRNA and protein levels in murine and human naive and activated T cells. Although 30 UTR shortening in proliferating cells is conserved between human and mouse, orthologous genes do not exhibit similar expression of alternative 30 UTR isoforms. We generally find that 30 UTR shortening is not accompanied by a corresponding change in mRNA and protein levels. This suggests that although 30 UTR shortening may lead to changes in the RNA-binding protein interactome, it has limited effects on protein output. 1 Computational and Systems Biology, Biozentrum, University of Basel and Swiss Institute of Bioinformatics, Basel CH-4056, Switzerland. 2 Department of Biomedicine, University of Basel and University Hospital Basel, Hebelstrasse 20, Basel CH-4031, Switzerland. 3 Infection Biology, Biozentrum, University of Basel, Basel CH-4056, Switzerland.
    [Show full text]