Genomic data in the All of Us Research Program

Comprehensively mapping the genetic basis of human disease across diverse individuals is a long-standing goal for the field of human genetics1-4. The All of Us Research Program is a longitudinal cohort study aiming to enrol a diverse group of at least one million individuals across the USA to accelerate biomedical research and improve human health5,6. Here we describe the programme's genomics data release of 245,388 clinical-grade genome sequences. This resource is unique in its diversity as 77% of participants are from communities that are historically under-represented in biomedical research and 46% are individuals from under-represented racial and ethnic minorities. All of Us identified more than 1 billion genetic variants, including more than 275 million previously unreported genetic variants, more than 3.9 million of which had coding consequences. Leveraging linkage between genomic data and the longitudinal electronic health record, we evaluated 3,724 genetic variants associated with 117 diseases and found high replication rates across both participants of European ancestry and participants of African ancestry. Summary-level data are publicly available, and individual-level data can be accessed by researchers through the All of Us Researcher Workbench using a unique data passport model with a median time from initial researcher registration to data access of 29 hours. We anticipate that this diverse dataset will advance the promise of genomic medicine for all.

The frequency of pathogenic variation in the All of Us cohort reveals ancestry-driven disparities

Disparities in data underlying clinical genomic interpretation is an acknowledged problem, but there is a paucity of data demonstrating it. The All of Us Research Program is collecting data including whole-genome sequences, health records, and surveys for at least a million participants with diverse ancestry and access to healthcare, representing one of the largest biomedical research repositories of its kind. Here, we examine pathogenic and likely pathogenic variants that were identified in the All of Us cohort. The European ancestry subgroup showed the highest overall rate of pathogenic variation, with 2.26% of participants having a pathogenic variant. Other ancestry groups had lower rates of pathogenic variation, including 1.62% for the African ancestry group and 1.32% in the Latino/Admixed American ancestry group. Pathogenic variants were most frequently observed in genes related to Breast/Ovarian Cancer or Hypercholesterolemia. Variant frequencies in many genes were consistent with the data from the public gnomAD database, with some notable exceptions resolved using gnomAD subsets. Differences in pathogenic variant frequency observed between ancestral groups generally indicate biases of ascertainment of knowledge about those variants, but some deviations may be indicative of differences in disease prevalence. This work will allow targeted precision medicine efforts at revealed disparities.

The All of Us Data and Research Center: Creating a Secure, Scalable, and Sustainable Ecosystem for Biomedical Research

The All of Us Research Program's Data and Research Center (DRC) was established to help acquire, curate, and provide access to one of the world's largest and most diverse datasets for precision medicine research. Already, over 500,000 participants are enrolled in All of Us, 80% of whom are underrepresented in biomedical research, and data are being analyzed by a community of over 2,300 researchers. The DRC created this thriving data ecosystem by collaborating with engaged participants, innovative program partners, and empowered researchers. In this review, we first describe how the DRC is organized to meet the needs of this broad group of stakeholders. We then outline guiding principles, common challenges, and innovative approaches used to build the All of Us data ecosystem. Finally, we share lessons learned to help others navigate important decisions and trade-offs in building a modern biomedical data platform.

A validated lineage-derived somatic truth data set enables benchmarking in cancer genome analysis

Existing cancer benchmark data sets for human sequencing data use germline variants, synthetic methods, or expensive validations, none of which are satisfactory for providing a large collection of true somatic variation across a whole genome. Here we propose a data set, Lineage derived Somatic Truth (LinST), of short somatic mutations in the HT115 colon cancer cell-line, that are validated using a known cell lineage that includes thousands of mutations and a high confidence region covering 2.7 gigabases per sample.

Megan Shand et al. present Lineage derived Somatic Truth (LinST), a validated data set of somatic mutations from a colon cancer cell line with a known lineage tree structure. They show that LinST can be used to benchmark true-positive and false-positive rates in somatic variant-calling pipelines applied to cancer genomic data.

Interplay of somatic alterations and immune infiltration modulates response to PD-1 blockade in advanced clear cell renal cell carcinoma

PD-1 blockade has transformed the management of advanced clear cell renal cell carcinoma (ccRCC), but the drivers and resistors of the PD-1 response remain incompletely elucidated. Here, we analyzed 592 tumors from patients with advanced ccRCC enrolled in prospective clinical trials of treatment with PD-1 blockade by whole-exome and RNA sequencing, integrated with immunofluorescence analysis, to uncover the immunogenomic determinants of the therapeutic response. Although conventional genomic markers (such as tumor mutation burden and neoantigen load) and the degree of CD8⁺ T cell infiltration were not associated with clinical response, we discovered numerous chromosomal alterations associated with response or resistance to PD-1 blockade. These advanced ccRCC tumors were highly CD8⁺ T cell infiltrated, with only 27% having a non-infiltrated phenotype. Our analysis revealed that infiltrated tumors are depleted of favorable PBRM1 mutations and enriched for unfavorable chromosomal losses of 9p21.3, as compared with non-infiltrated tumors, demonstrating how the potential interplay of immunophenotypes with somatic alterations impacts therapeutic efficacy.

Abstract 5108: Somatic small variant and copy number alteration calling with the Genome Analysis Toolkit

Somatic small mutations, SNVs or Indels, and copy number alterations are the two categories of mutations with the largest impact on cancer tumors. The Broad Institute has released somatic variant calling workflows for small mutations (M2) and copy number alterations (ModelSegments) based on the Genome Analysis Toolkit (GATK). The suite of workflows can call variants in capture or whole-genome sequencing data and will include functional annotations (Funcotator), such as protein change (for small variants) and impacted gene (for all variants). Common artifacts in sequencing data, such as those arising from oxidative DNA damage, FFPE/deamination, or mapping errors, are corrected automatically. Evaluation of the workflows is standardized and repeatable, which allows tracking of performance across versions, both detection performance (e.g. sensitivity, precision), as well as runtime performance (e.g. CPU and RAM usage). A matched normal is not required for a given tumor sample, since the workflows can leverage pre-processed panels of normals (PoNs). The workflows are freely available, are portable (i.e. can be run on local, on-prem, or cloud compute), are optimized for cost reduction, and can be tuned to optimally leverage available compute.The measured sensitivity of M2 was at least 0.93 for small somatic nucleotide variants (SNVs) and 0.83 for small insertions/deletions (Indels) on DREAM1, DREAM2, and DREAM3 challenges, and on a titrated mixture of germline samples (>=100x depth, AF = 0.2). The measured precision of M2 ranged from 0.91 to 0.98 on DREAM1, DREAM2, and DREAM3 for both SNVs and Indels. The false positive rate (FPR) of M2 was between 0.03 and 0.21 FP/Mb for SNVs, and between 0.0 and 0.1 FP/Mb for indels, on twelve paired, replicate normal-normal samples. The cost of the M2 workflow is about USD$1.15 for a pair of 35x WGS matched tumor-normal samples, using Google Cloud Compute, and required about 32 hours of CPU time on a single core with 3GB RAM.

The measured sensitivity of ModelSegments was at least 0.91 for deletions and amplifications across three cohorts of TCGA whole-exome samples (Stomach adenocarcinoma N=39, Thyroid carcinoma N=50, and Lung adenocarcinoma N=60). The measured specificity for the same set of cohorts was at least 0.96 for both deletions and amplifications. All results reported here were using the corresponding SNP Array results as a truth set.

GATK MS cost was approximately USD$0.65 on a 30x WGS pair using Google Cloud Compute and required about 6 hours of CPU time with a single core. The RAM usage was varied automatically in the workflow to minimize cost, but was in the range of 2-13GB.

Citation Format: Lee Lichtenstein, Jonn Smith, David Benjamin, Aaron Chevalier, Kristian Cibulskis, Samuel K. Lee, Eric Banks. Somatic small variant and copy number alteration calling with the Genome Analysis Toolkit [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 5108.

Abstract 2287: Precise common and rare germline CNV calling with GATK

We propose, implement, and evaluate a novel method (GATK gCNV) for accurate discovery of rare and common copy-number variations (CNVs) from read-depth data obtained from whole genome sequencing (WGS), whole exome sequencing (WES), or custom gene panels. GATK gCNV utilizes a sophisticated Bayesian model to learn bias factors arising from sequencing and library preparation. This model accounts for ploidy of sex chromosomes and autosomal aneuploidies, treats GC bias probabilistically, and automatically determines the necessary level of model complexity in a data-driven manner. Unlike most existing read-depth methods, GATK gCNV maintains a high level of sensitivity in common CNV regions, due to a hierarchical hidden Markov model used for accurate genotyping of multi-allelic loci. Furthermore, GATK gCNV performs bias modeling and CNV discovery simultaneously and self-consistently, resulting in significantly improved sensitivity and precision. Our implementation utilizes the PyMC3/Theano framework for performing automatic differentiation variational inference (ADVI). In addition, GATK gCNV automatically scatters large tasks across multiple machines using the Cromwell/WDL framework, enabling the scalable processing of large cohorts.

We use GATK gCNV to compute copy-number transmission and de novo rates in a cohort of WES trios and observe consistency with observed population metrics. Furthermore, we benchmark GATK gCNV for sensitivity, precision, and reproducibility of both rare and common CNV calls. Using a cohort of WES blood normal samples from The Cancer Genome Atlas (TCGA), we show that GATK gCNV calls are in remarkable concordance with Genome STRiP calls on matched WGS data and gCNV substantially outperforms XHMM and CODEX. We also validate GATK gCNV calls on WGS data against calls made using PacBio long reads.

Citation Format: Mehrtash Babadi, Samuel K. Lee, Andrey Smirnov, Lee Lichtenstein, Laura D. Gauthier, Daniel P. Howrigan, Timothy Poterba. Precise common and rare germline CNV calling with GATK [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 2287.

Abstract 3581: GATK ACNV: allelic copy-number variation discovery from SNPs and coverage data

The presence of somatic copy-number alterations in tumor genomes can be used to predict both patient sensitivity to treatments as well as outcomes. The inclusion of allelic data improves statistical power to detect copy-number events and allows for discovery of copy-neutral events. We present GATK ACNV, an allelic copy-number variation method built on the Genome Analysis Toolkit. ACNV is a tool for detecting somatic copy-number activity from whole exome and whole genome sequencing data by segmenting the genome into regions of constant copy number and estimating copy ratio and minor-allele fraction in those regions.

ACNV uses a novel probabilistic model to account for reference bias (optionally using a panel of normals), which improves the estimation of minor-allele fraction. We combine this with the coverage model from GATK CNV by segmenting with a unified hidden Markov model, improving the statistical power to detect copy-number variation.

We validate ACNV using a purity series of the cell line HCC1143 and cancer samples from The Cancer Genome Atlas. Our results show that ACNV is able to discover regions of somatic copy-number activity accurately and with high resolution in both whole exome and whole genome sequencing data.

Citation Format: Aaron Chevalier, Lee Lichtenstein, Andrey Smirnov, Samuel K. Lee, Mehrtash Babidi, David I. Benjamin, Valentin Ruano-Rubio. GATK ACNV: allelic copy-number variation discovery from SNPs and coverage data [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 3581. doi:10.1158/1538-7445.AM2017-3581

Comprehensive, Integrative Genomic Analysis of Diffuse Lower-Grade Gliomas

BACKGROUND

Diffuse low-grade and intermediate-grade gliomas (which together make up the lower-grade gliomas, World Health Organization grades II and III) have highly variable clinical behavior that is not adequately predicted on the basis of histologic class. Some are indolent; others quickly progress to glioblastoma. The uncertainty is compounded by interobserver variability in histologic diagnosis. Mutations in IDH, TP53, and ATRX and codeletion of chromosome arms 1p and 19q (1p/19q codeletion) have been implicated as clinically relevant markers of lower-grade gliomas.

METHODS

We performed genomewide analyses of 293 lower-grade gliomas from adults, incorporating exome sequence, DNA copy number, DNA methylation, messenger RNA expression, microRNA expression, and targeted protein expression. These data were integrated and tested for correlation with clinical outcomes.

RESULTS

Unsupervised clustering of mutations and data from RNA, DNA-copy-number, and DNA-methylation platforms uncovered concordant classification of three robust, nonoverlapping, prognostically significant subtypes of lower-grade glioma that were captured more accurately by IDH, 1p/19q, and TP53 status than by histologic class. Patients who had lower-grade gliomas with an IDH mutation and 1p/19q codeletion had the most favorable clinical outcomes. Their gliomas harbored mutations in CIC, FUBP1, NOTCH1, and the TERT promoter. Nearly all lower-grade gliomas with IDH mutations and no 1p/19q codeletion had mutations in TP53 (94%) and ATRX inactivation (86%). The large majority of lower-grade gliomas without an IDH mutation had genomic aberrations and clinical behavior strikingly similar to those found in primary glioblastoma.

CONCLUSIONS

The integration of genomewide data from multiple platforms delineated three molecular classes of lower-grade gliomas that were more concordant with IDH, 1p/19q, and TP53 status than with histologic class. Lower-grade gliomas with an IDH mutation either had 1p/19q codeletion or carried a TP53 mutation. Most lower-grade gliomas without an IDH mutation were molecularly and clinically similar to glioblastoma. (Funded by the National Institutes of Health.).

Oncotator: cancer variant annotation tool

Oncotator is a tool for annotating genomic point mutations and short nucleotide insertions/deletions (indels) with variant- and gene-centric information relevant to cancer researchers. This information is drawn from 14 different publicly available resources that have been pooled and indexed, and we provide an extensible framework to add additional data sources. Annotations linked to variants range from basic information, such as gene names and functional classification (e.g. missense), to cancer-specific data from resources such as the Catalogue of Somatic Mutations in Cancer (COSMIC), the Cancer Gene Census, and The Cancer Genome Atlas (TCGA). For local use, Oncotator is freely available as a python module hosted on Github (https://github.com/broadinstitute/oncotator). Furthermore, Oncotator is also available as a web service and web application at http://www.broadinstitute.org/oncotator/.

Abstract 4692: Relationships between somatic genomic alterations, tumor stage and progression-free survival in cervical cancer

Background:

Cervical cancer is a major public health problem worldwide. We have recently identified novel significantly recurrent somatic mutations in HLA-B, ERBB2 and MAPK1 in cervical squamous cell carcinomas. However, the significance of somatic mutations and copy number alterations for clinical phenotype in cervical cancer is not well understood. This study seeks to identify relationships between somatic genomic alterations, epidemiological exposures, tumor stage (FIGO) and progression-free survival in cervical carcinomas.

Methods:

Cervical tumors were surgically resected or biopsied from 100 Norwegian women with tumor stages I - IV. Patients were subsequently followed for 0-109 months (mean=25 months). Whole exome sequencing (WES) was performed on DNA extracted from tumors and corresponding normal blood. Somatic single nucleotide variants and small insertion/deletions were identified by the MuTect and Indelocator algorithms respectively. The ABSOLUTE algorithm was used to classify mutations as clonal or subclonal. Somatic copy number (CN) data were derived from WES data using the CapSeg algorithm, and significantly recurrent CN alterations were identified by GISTIC2.0 analysis (q<0.25). HPV typing was done by the multiplex fluorescent f-HPV DNA and MassARRAY assays. The log-rank test was used to compare survival curves.

Results:

Non-localized tumors (FIGO stages ≥ II) were associated with focal amplification of the FGFR2 gene on chromosomal cytoband 10q26 (GISTIC q = 0.18531). Indeed, 6 of 8 (75%) tumors with FGFR2 amplification were non-localized, in contrast to 16 of 92 (17%) tumors without FGFR2 amplification (p = 0.001). In addition, patients with somatic ERBB2 mutations and/or amplifications (p = 0.04), somatic TP53 mutations and/or deletions (p = 0.04), or infection with multiple HPV types (p = 0.02) had poorer prognosis (progression-free survival). We also observed a trend for higher frequency of subclonal driver events in patients with poorer survival (p = 0.07).

Conclusion:

We have identified novel relationships between somatic genomic alterations, tumor stage and patient prognosis in cervical cancer. Our data suggest a potential for exploring FGFR2 inhibition in non-localized cervical carcinomas with FGFR2 alterations in a clinical trial context.

Citation Format: Akinyemi I. Ojesina, Bradley Murray, Line Bjorge, Kathrine Woie, Camilla Krakstad, Lee Lichtenstein, Chandra Sekhar Pedamallu, Amaro Taylor-Weiner, Samuel S. Freeman, Andrew D. Cherniack, Michael S. Lawrence, Kristian Cibulskis, Scott L. Carter, Heather Walline, Thomas E. Carey, Olav K. Vintermyr, Bjorn Bertelsen, Christopher P. Crum, Gad Getz, Matthew Meyerson, Helga B. Salvesen. Relationships between somatic genomic alterations, tumor stage and progression-free survival in cervical cancer. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 4692. doi:10.1158/1538-7445.AM2014-4692

Landscape of genomic alterations in cervical carcinomas

Cervical cancer is responsible for 10-15% of cancer-related deaths in women worldwide. The aetiological role of infection with high-risk human papilloma viruses (HPVs) in cervical carcinomas is well established. Previous studies have also implicated somatic mutations in PIK3CA, PTEN, TP53, STK11 and KRAS as well as several copy-number alterations in the pathogenesis of cervical carcinomas. Here we report whole-exome sequencing analysis of 115 cervical carcinoma-normal paired samples, transcriptome sequencing of 79 cases and whole-genome sequencing of 14 tumour-normal pairs. Previously unknown somatic mutations in 79 primary squamous cell carcinomas include recurrent E322K substitutions in the MAPK1 gene (8%), inactivating mutations in the HLA-B gene (9%), and mutations in EP300 (16%), FBXW7 (15%), NFE2L2 (4%), TP53 (5%) and ERBB2 (6%). We also observe somatic ELF3 (13%) and CBFB (8%) mutations in 24 adenocarcinomas. Squamous cell carcinomas have higher frequencies of somatic nucleotide substitutions occurring at cytosines preceded by thymines (Tp*C sites) than adenocarcinomas. Gene expression levels at HPV integration sites were statistically significantly higher in tumours with HPV integration compared with expression of the same genes in tumours without viral integration at the same site. These data demonstrate several recurrent genomic alterations in cervical carcinomas that suggest new strategies to combat this disease.

Mutational heterogeneity in cancer and the search for new cancer-associated genes

Major international projects are underway that are aimed at creating a comprehensive catalogue of all the genes responsible for the initiation and progression of cancer. These studies involve the sequencing of matched tumour-normal samples followed by mathematical analysis to identify those genes in which mutations occur more frequently than expected by random chance. Here we describe a fundamental problem with cancer genome studies: as the sample size increases, the list of putatively significant genes produced by current analytical methods burgeons into the hundreds. The list includes many implausible genes (such as those encoding olfactory receptors and the muscle protein titin), suggesting extensive false-positive findings that overshadow true driver events. We show that this problem stems largely from mutational heterogeneity and provide a novel analytical methodology, MutSigCV, for resolving the problem. We apply MutSigCV to exome sequences from 3,083 tumour-normal pairs and discover extraordinary variation in mutation frequency and spectrum within cancer types, which sheds light on mutational processes and disease aetiology, and in mutation frequency across the genome, which is strongly correlated with DNA replication timing and also with transcriptional activity. By incorporating mutational heterogeneity into the analyses, MutSigCV is able to eliminate most of the apparent artefactual findings and enable the identification of genes truly associated with cancer.

Abstract 5135: Analysis of formalin-fixed paraffin-embedded (FFPE) samples

The availability of large numbers of formalin-fixed, paraffin-embedded (FFPE) tumor DNA samples, together with accurate and specific clinical annotations, can be a great resource for in-depth analysis of correlations of somatic and germline DNA alterations with clinical outcome. However, the effects of this tissue preservation technique on next-generation DNA sequencing technologies and downstream analyses are still being investigated. Furthermore, usually the archives of such FFPE samples do not include matched normal tissue, which complicates the process of identifying somatic alterations. We are in the process of analyzing lung adenocarcinoma, colorectal and prostate patient datasets, which contain both an FFPE and a fresh-frozen tumor together with a matched blood normal sample. The availability of both a frozen and an FFPE tumor sample enables us to compare the detected somatic point mutations and indels, taking into account the fact that these pairs are produced from different aliquots of DNA. This analysis includes calculating the power to detect variants in the FFPE sample given that they have been observed in the fresh frozen tumor, as well as the validation rate of clonal SNVs between the two samples. We are also developing a method for classification of SNVs as germline or somatic without a paired normal sample by taking advantage of the fact that most tumor samples contain a substantial fraction of normal cells. Because stromal contamination has a different effect on the allelic fraction of somatic vs. germline SNVs, we can apply algorithms to estimate tumor purity and absolute somatic copy numbers in order to distinguish the two types of events. We present preliminary results from these analyses.

Citation Format: Petar Stojanov, Scott L. Carter, Ivan Imaz Rosshandler, Andrey Sivachenko, Alberto Salido Guadarrama, Karla Vazquez, Sandra Romero Cordoba, Kristian Cibulskis, Carrie Sougnez, Douglas Voet, Gordon Saksena, Lee Lichtenstein, Lihua Zou, Scott Frazer, Chip Stewart, Rameen Beroukhim, Matthew Meyerson, Michael S. Lawrence, Gad Getz. Analysis of formalin-fixed paraffin-embedded (FFPE) samples. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 5135. doi:10.1158/1538-7445.AM2013-5135

Abstract 4604: Landscape of human and viral genomic alterations in cervical carcinomas

Background: Cervical cancer is a major public health problem worldwide. The etiological role of human papilloma virus (HPV) infections in cervical cancer is well established. However, HPV infection is insufficient to account for the development of cervical cancers because only 2 % of women infected with HPV eventually develop invasive carcinomas. We have therefore initiated a large scale sequencing effort to provide comprehensive data on the global landscape of genomic aberrations and HPV variants that contribute to cervical cancer. Methods: We have carried out comprehensive whole exome sequencing analyses on 120 tumor-normal paired samples from Mexico and Norway. We also carried out transcriptome and whole genome sequencing analyses on a subset of the patients (81 and 14 individuals respectively). Results: The aggregate mutation rate across the dataset was 3.8 per megabase (Mb), with the most common mutations being C to T/G in the Tp*C dinucleotide context, at a rate of 15 mutations per Mb. In all, 18,037 mutations were found across the entire dataset, including 11,536 missense, 984 nonsense, 4691 silent, 231 splice site, 32 translation start site mutations, as well 443 deletions and 142 insertions. MutSig analyses to identify genes that were mutated at statistically significant frequencies across our dataset revealed 11 genes to be recurrently mutated with a false discovery rate of q<0.1 after correction for multiple hypothesis testing (and RNASeq-based evidence of robust gene expression). The most significantly mutated genes encode for members of the PIK3CA/PTEN and RAS/RAF/MAPK signaling pathways, as well as the major histocompatibility complex (MHC). We have also uncovered novel patterns of HPV transcript abundance and sites of recurrent HPV integration in cell cycle related genes. In addition, our whole genome sequencing data suggests that HPV-negative p53-mutant tumors harbor high frequencies of genomic rearrangements. Conclusion: The comprehensive catalogue of genomic alterations provided by this project reveals potential novel therapeutic targets in cervical carcinomas. Our data also sets the stage for improving diagnostic and preventive strategies, especially in resource-limited settings with the highest incidence of cervical cancer.

Citation Format: Akinyemi I. Ojesina, Lee Lichtenstein, Lauren Ambrogio, Kristian Cibulskis, Samuel Freeman, Chandra Sekhar Pedamallu, Bjørn Bertelsen, Ivan Imaz, Karla Vazquez, Alberto Salido Guadarrama, Victor Treviño, Sandra Romero-Cordoba, Fujiko Duke, Bethany Kaplan, Iram Rodriguez, Magali Espinosa Castilla, Katherine Woie, Line Bjorge, Elisabeth Wik, Mari K. Halle, Erling Høivik, Camilla Krakstad, Gabriela Gómez Macías, María de Lourdes Garza Rodríguez, Jorge Vazquez, Carlos Rodea, Adrian Cravioto, Maria L. Cortes, Heidi Greulich, Christopher P. Crum, Lars Akslen, Hugo Barrera Saldaña, Jorge Melendez-Zajgla, Gad Getz, Helga B. Salvesen, Matthew L. Meyerson. Landscape of human and viral genomic alterations in cervical carcinomas. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 4604. doi:10.1158/1538-7445.AM2013-4604

Discovery and characterization of artifactual mutations in deep coverage targeted capture sequencing data due to oxidative DNA damage during sample preparation

As researchers begin probing deep coverage sequencing data for increasingly rare mutations and subclonal events, the fidelity of next generation sequencing (NGS) laboratory methods will become increasingly critical. Although error rates for sequencing and polymerase chain reaction (PCR) are well documented, the effects that DNA extraction and other library preparation steps could have on downstream sequence integrity have not been thoroughly evaluated. Here, we describe the discovery of novel C > A/G > T transversion artifacts found at low allelic fractions in targeted capture data. Characteristics such as sequencer read orientation and presence in both tumor and normal samples strongly indicated a non-biological mechanism. We identified the source as oxidation of DNA during acoustic shearing in samples containing reactive contaminants from the extraction process. We show generation of 8-oxoguanine (8-oxoG) lesions during DNA shearing, present analysis tools to detect oxidation in sequencing data and suggest methods to reduce DNA oxidation through the introduction of antioxidants. Further, informatics methods are presented to confidently filter these artifacts from sequencing data sets. Though only seen in a low percentage of reads in affected samples, such artifacts could have profoundly deleterious effects on the ability to confidently call rare mutations, and eliminating other possible sources of artifacts should become a priority for the research community.

Abstract 5056: A glimpse into the somatic mutation landscape of melanoma through exome sequencing of 121 tumor-normal pairs

Melanoma is an aggressive skin cancer of melanocytic origin characterized by high metastatic potential and mutation rate. Affording a survey of the wide breadth of genomic lesions found in melanoma, we present here an analysis of the somatic mutations discovered in the sequenced exomes of 121 melanoma tumor-normal pairs. We identify frequent genomic alterations both in genes previously implicated in melanoma (BRAF, NRAS, TP53, CDKN2A, PTEN) as well as in several genes whose role in melanoma tumorigenesis has not yet been established and thus are of particular interest. To do so we implement a novel method to increase the identification of genes that are significantly recurrently mutated in melanoma in the setting of its exceptionally high mutation rate. A preponderance of C>T transitions (∼85%) in the observed mutational profile reflects a history of DNA damage due to UV radiation, though the majority of somatic mutations in known melanoma genes are not C>T events. Our study broadens understanding of the genomic lesions involved in melanoma tumorigenesis, and we expect our analysis approach to inform future genomic studies of cancer lineages with similarly high mutation rates.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 5056. doi:1538-7445.AM2012-5056

Growth and structure of crystalline silica sheet on Ru(0001)

Thin SiO₂ films were grown on a Ru(0001) single crystal and studied by photoelectron spectroscopy, infrared spectroscopy and scanning probe microscopy. The experimental results in combination with density functional theory calculations provide compelling evidence for the formation of crystalline, double-layer sheet silica weakly bound to a metal substrate.

Transcription of the human uncoupling protein 3 gene is governed by a complex interplay between the promoter and intronic sequences

AIMS/HYPOTHESIS

Uncoupling protein (UCP) 3 is an inner mitochondrial membrane transporter mainly produced in skeletal muscle in humans. UCP3 plays a role in fatty acid metabolism and energy homeostasis and modulates insulin sensitivity. In humans, UCP3 content is higher in fast-twitch glycolytic muscle than in slow-twitch oxidative muscle and is dysregulated in type 2 diabetes. Here, we studied the molecular mechanisms determining human UCP3 levels in skeletal muscle and their regulation by fasting in transgenic mice.

METHODS

We produced a series of transgenic lines with constructs bearing different putative regulatory regions of the human UCP3 gene, including promoter and intron sequences. UCP3 mRNA and reporter gene expression and activity were measured in different skeletal muscles and tissues.

RESULTS

The profile of expression and the response to fasting and thyroid hormone of human UCP3 mRNA in transgenic mice with 16 kb of the human UCP3 gene were similar to that of the endogenous human gene. Various parts of the UCP3 promoter did not confer expression in transgenic lines. Inclusion of intron 1 resulted in an expression profile in skeletal muscle that was identical to that of human UCP3 mRNA. Further dissection of intron 1 revealed that distinct regions were involved in skeletal muscle expression, distribution among fibre types and response to fasting.

CONCLUSIONS/INTERPRETATION

The control of human UCP3 transcription in skeletal muscle is not solely conferred by the promoter, but depends on several cis-acting elements in intron 1, suggesting a complex interplay between the promoter and intronic sequences.

A Feedback-Controlled Interface for Treadmill Locomotion in Virtual Environments

Virtual environments (VEs) allow safe, repeatable, and controlled evaluations of obstacle avoidance and navigation performance of people with visual impairments using visual aids. Proper simulation of mobility in a VE requires an interface, which allows subjects to set their walking pace. Using conventional treadmills, the subject can change their walking speed by pushing the tread with their feet, while leveraging handrails or ropes (self-propelled mode). We developed a feedback-controlled locomotion interface that allows the VE workstation to control the speed of the treadmill, based on the position of the user. The position and speed information is also used to implement automated safety measures, so that the treadmill can be halted in case of erratic behavior. We compared the feedback-controlled mode to the self-propelled mode by using speed-matching tasks (follow a moving object or match the speed of an independently moving scene) to measure the efficacy of each mode in maintaining constant subject position, subject control of the treadmill, and subject pulse rates. Additionally, we measured the perception of speed in the VE on each mode. The feedback-controlled mode required less physical exertion than self-propelled. The average position of subjects on the feedback-controlled treadmill was always within a centimeter of the desired position. There was a smaller standard deviation in subject position when using the self-propelled mode than when using the feedback-controlled mode, but the difference averaged less than six centimeters across all subjects walking at a constant speed. Although all subjects underestimated the speed of an independently moving scene at higher speeds, their estimates were more accurate when using the feedback-controlled treadmill than the self-propelled.

Relationship between histamine and physiological changes during the early response to nasal antigen provocation

To investigate the temporal relationships of mediator release and physiological changes during the early response to allergen, we challenged allergic individuals intranasally with antigen and followed their responses. This was done by using small filter paper disks to challenge one nostril and collect secretions from both the challenged and the contralateral nostril, thus enabling us to evaluate the nasonasal reflex. There was a significant increase in sneezing after allergen challenge that peaked within 2 min and returned to baseline. The weights of nasal secretions as well as nasal symptoms increased immediately and remained significantly elevated for 20 min in both nostrils. Nasal airway resistance increased slowly, reaching its peak at approximately 6 min after challenge on the ipsilateral side, but it did not change on the contralateral side. Histamine levels peaked 30 s after removal of the allergen disk on the side of challenge, whereas albumin levels peaked after those of histamine. Lactoferrin paralleled the increase in secretion weights and occurred in both nostrils. Increasing doses of antigen produced dose-dependent increases in all parameters, whereas control challenges produced no response. These studies describe a human model for the evaluation of the allergic response that is capable of simultaneously measuring mediator release and the physiological response, including the nasonasal reflex. This model should prove useful in studying the mechanism of allergic rhinitis in humans.

In vivo release of inflammatory mediators by hyperosmolar solutions

Hyperosmolar environments induce histamine release from mast cells and basophils in vitro. To assess whether the same stimulus induces mediator release in vivo, 15 healthy human volunteers underwent nasal challenges with instilled solutions of differing osmolalities: lactated Ringer's solution (257 +/- 3 mOsm/kg), isosmolar mannitol (277 +/- 6 mOsm/kg), and hyperosmolar mannitol (869 +/- 8 mOsm/kg). The effect of these challenges on the volume, osmolality, and inflammatory mediator content of subsequent 5-ml isosmolar lavages was determined. The volumes of lavages returned after hyperosmolar challenges were significantly greater than those after isosmolar challenges (5.5 +/- 0.2 ml versus 4.2 +/- 0.1 ml; p less than 0.01) and these lavage solutions had higher osmolalities. Even when corrected for increased volumes, the lavages after hyperosmolar challenges contained significantly higher quantities of inflammatory mediators such as histamine (29.0 versus 10.1 ng; p less than 0.01), TAME-esterase activity (32.7 versus 11.1 cpm x 10(-3); p less than 0.01), and immunoreactive leukotrienes (9.9 versus 3.4 ng; p less than 0.01). The changes in mediators were dose dependent in that incremental increase in challenge osmolality were associated with incremental increases in histamine release. Therefore, when exposed to hyperosmolar stimuli in vivo, the nasal respiratory airway releases inflammatory mediators and fluid rapidly shifts into the airway lumen. It has been suggested that the mediator release observed on breathing cold and dry air is due to increased osmolality of airway secretions; the present data confirm that osmotic variations at the airway surface can provide an adequate stimulus for cell activation.

The effect of a topical tricyclic antihistamine on the response of the nasal mucosa to challenge with cold, dry air and histamine

We have previously demonstrated that azatadine, a tricyclic antihistamine, known also to inhibit mediator release from mast cells and basophils in vitro, inhibits the early release of histamine and other mediators after nasal challenge with antigen. In this article, we studied the effect of azatadine on preventing the release of histamine after nasal challenge with cold, dry air (CDA) and its effect on antagonizing nasal challenge with histamine. With histamine challenge, azatadine inhibited symptoms (sneezing, nasal congestion, and rhinorrhea) and the increase in the level of albumin in nasal secretions (p less than 0.01 all). With challenge with CDA, the drug had no effect on either symptoms or histamine and N-tosyl-L-arginine methyl ester-esterase release. Although the patterns of mediator release after CDA and after the early reaction to antigen are similar, the pharmacologic control differs, suggesting different mechanisms of induction of histamine release from mast cells.

In vitro studies of antigen-induced bronchospasm: effect of antihistamine and SRS-A antagonist on response of sensitized guinea pig and human airways to antigen

Exposure of sensitized guinea pig tracheal rings or human bronchial strips to specific antigen in vitro resulted in a rapidly developing, prolonged contraction that was resistant to washing. Treatment of the tissue with diphenhydramine, a histamine H1 antagonist, before antigen delayed the onset and decreased the amplitude of the initial phase of the contraction but did not reduce the duration. Diphenhydramine treatment after development of the contraction did not relax the airway tissue. Antigen-induced histamine release from guinea pig trachea and from human bronchus was complete within the initial 15% of the duration of the contraction. Treatment of sensitized airway tissue with FPL 55712, a SRS-A antagonist, before antigen selectively inhibited the prolonged phase of the response. FPL 55712 administration after the development of antigen-induced contraction resulted in relaxation. These data suggest that both histamine and SRS-A are involved in the response of sensitized guinea pig and human airway tissue to antigen, with histamine mediating the early phase of the contraction and SRS-A primarily mediating the protracted phase.

Ewing's Sarcoma of Bone

Images

Necrotizing Arterial Lesions Resembling Those of Periarteritis Nodosa and Focal Visceral Necrosis Following Administration of Sulfathiazole: Report of a Case

Images

Chondrosarcoma of Bone

Images

Benign Chondroblastoma of Bone: A Reinterpretation of the So-Called Calcifying or Chondromatous Giant Cell Tumor

Images

Non-osteogenic fibroma of bone

Images

Eosinophilic granuloma of bone: With report of a case

Images

Pathological changes following therapeutic hyperthermia: Report of a Case

Images

The Virus of Lymphogranuloma Inguinale

Images