|J Pathol Inform 2022,
SCHOOL: Software for Clinical Health in Oncology for Omics Laboratories
Chelsea K Raulerson1, Erika C Villa2, Jeremy A Mathews2, Benjamin Wakeland2, Yan Xu3, Jeffrey Gagan3, Brandi L Cantarel2
1 Bioinformatics Core Facility, University of Texas Southwestern Medical Center, Dallas, Texas, USA; Department of Pathology, University of Texas Southwestern Medical Center, Dallas, Texas, USA
2 Bioinformatics Core Facility, University of Texas Southwestern Medical Center, Dallas, Texas, USA; Department of Bioinformatics, University of Texas Southwestern Medical Center, Dallas, Texas, USA
3 Department of Pathology, University of Texas Southwestern Medical Center, Dallas, Texas, USA
|Date of Submission||18-Mar-2021|
|Date of Decision||01-Sep-2021|
|Date of Acceptance||04-Oct-2021|
|Date of Web Publication||05-Jan-2022|
Dr. Brandi L Cantarel
Bioinformatics Core Facility, Department of Bioinformatics, UT Southwestern Medical Center, 5323 Harry Hines Boulevard, EB3.212, MC 9365, Dallas, Texas 75390-9365.
Source of Support: None, Conflict of Interest: None
| Abstract|| |
Bioinformatics analysis is a key element in the development of in-house next-generation sequencing assays for tumor genetic profiling that can include both tumor DNA and RNA with comparisons to matched-normal DNA in select cases. Bioinformatics analysis encompasses a computationally heavy component that requires a high-performance computing component and an assay-dependent quality assessment, aggregation, and data cleaning component. Although there are free, open-source solutions and fee-for-use commercial services for the computationally heavy component, these solutions and services can lack the options commonly utilized in increasingly complex genomic assays. Additionally, the cost to purchase commercial solutions or implement and maintain open-source solutions can be out of reach for many small clinical laboratories. Here, we present Software for Clinical Health in Oncology for Omics Laboratories (SCHOOL), a collection of genomics analysis workflows that (i) can be easily installed on any platform; (ii) run on the cloud with a user-friendly interface; and (iii) include the detection of single nucleotide variants, insertions/deletions, copy number variants (CNVs), and translocations from RNA and DNA sequencing. These workflows contain elements for customization based on target panel and assay design, including somatic mutational analysis with a matched-normal, microsatellite stability analysis, and CNV analysis with a single nucleotide polymorphism backbone. All of the features of SCHOOL have been designed to run on any computer system, where software dependencies have been containerized. SCHOOL has been built into apps with workflows that can be run on a cloud platform such as DNANexus using their point-and-click graphical interface, which could be automated for high-throughput laboratories.
Keywords: Bioinformatics, cancer, NGS
|How to cite this article:|
Raulerson CK, Villa EC, Mathews JA, Wakeland B, Xu Y, Gagan J, Cantarel BL. SCHOOL: Software for Clinical Health in Oncology for Omics Laboratories. J Pathol Inform 2022;13:1
|How to cite this URL:|
Raulerson CK, Villa EC, Mathews JA, Wakeland B, Xu Y, Gagan J, Cantarel BL. SCHOOL: Software for Clinical Health in Oncology for Omics Laboratories. J Pathol Inform [serial online] 2022 [cited 2022 Jan 28];13:1. Available from: https://www.jpathinformatics.org/text.asp?2022/13/1/1/334714
| Introduction|| |
The rapid expansion and decreasing cost of next-generation sequencing (NGS) technology present an opportunity to improve the diagnosis and treatment of cancer through identifying tumor-specific mutations and enabling physicians to adapt treatment plans that suit the unique molecular profile of each patient., To address the evolving list of clinically actionable and prognostic biomarkers in the treatment of cancer, academic clinical laboratories have developed sequencing assays with varying size gene panels (100–1600 genes), with consistent quality to detect relevant genetic variants.
Bioinformatics analysis of sequence data includes two phases: (i) primary analysis, which converts the raw sequencing reads into predicted genetic variants and read abundances, and (ii) secondary analysis, which is customized for each clinical assay to maximize sensitivity and specificity by identifying artifact and poor-quality variant predictions. For a quick turn-around-time, the primary analysis is computationally demanding and requires computational resources with high memory (>32 GB) and multiple processors using a local high-performance cluster or cloud-computing resources. Furthermore, the primary analysis requires multiple elements for complete somatic variant detection including single nucleotide variants (SNVs), insertions/deletions (indels), copy number variants (CNVs), and structural variants such as translocations, large deletions, internal tandem duplications, and differences in microsatellite length. Because each assay might include different elements for detecting these different variant types, the primary analysis should be customizable for a variety of assays. Secondary analysis is much less resource-intensive, can be run on a desktop computer, and should be tailored to the needs of the specific assay.
Both commercial and open-source solutions have been introduced to address primary analysis needs in cancer genomics. Commercially developed bioinformatics pipelines are proprietary, often have limited options for customization, and require licensing, which can increase computational costs. By contrast, open-source solutions are usually customizable and lack licensing costs., However, common best-practice tools for variant detection, such as BWA and GATK4, require computational programming expertise to run in a Linux command line environment. Some commonly used open-source tools for more complex variant detection lack thorough documentation, continued support for development, or the flexibility to process varied data types (tumor-only samples versus matched tumor/normal control). Additionally, many existing software tools are difficult to install and maintain, due to the sometimes difficult installation of software dependencies, which make them sensitive to updates and changes to default programs. Finally, these tools are not natively packaged as an end-to-end analysis pipeline, which starts with raw sequencing reads and results in predicted variants. A user-friendly interface is critical in a customizable, open-source bioinformatics pipeline that is easy to install and run without specialized computational training.
In order to address the need for an end-to-end customizable bioinformatics pipeline for the primary analysis of sequence data, we have developed a collection of analysis workflows for NGS data and the detection of genetic alterations in cancer called Software for Clinical Health in Oncology for Omics Laboratories (SCHOOL) [Figure 1]. SCHOOL: (i) can detect SNVs, indels, CNVs, and translocations from RNA and DNA sequencing, (ii) has tools for mutational profiling and omics integration, and (iii) is designed to be easy to run on local computing resources or the cloud, with all software packaged alongside dependencies, so that they can work on any system where singularity or docker programming packages are available. We have optimized each step to use a minimal amount of RAM and processors to reduce computation costs on the cloud. These workflows contain the steps necessary to complete primary NGS analysis, including variant detection and annotation. Furthermore, these workflows can execute on a local cluster using Nextflow, a command-line workflow manager, or on the cloud, https://platform.dnanexus.com/panx/projects/FvPKK200Y9g81KqkKjJ9X818/data/, using the DNANexus applets and workflows code.
|Figure 1: Overview of SCHOOL workflow from sequencing through reporting. In SCHOOL, data flow from the sequencer into the primary analysis pipeline, which includes quality control, alignment, and variant calling appropriate for the sample type. Then, in secondary analysis, the variants can be annotated for eventual clinical reports|
Click here to view
| Technical Background|| |
The primary analysis of sequence data for the detection of somatic variants in tumor samples requires five main steps including (i) alignment of the raw sequencing data to a reference genome, (ii) identification of SNVs and indels in DNA, (iii) identification of CNVs in DNA, (iv) identification of copy number and structural variants in RNA and DNA, and (v) annotation and the prediction of effect. For each step, there are many considerations for the bioinformatics workflow that can affect accuracy.
For alignments, the user should consider the genome reference, removal of duplicate reads, and the alignment program. There are currently two available versions of the reference human genome: GRCh37 (hg19), released in 2009, and GRCh38 (hg38), released in 2013. GRCh38 has been shown to produce more accurate alignments. Sequence duplicates are caused when polymerase chain reaction (PCR) errors are amplified. Duplicates can be marked or removed, so that they are ignored in downstream steps, Picard (http://broadinstitute.github.io/picard/) can be used to mark duplicates, and Samtools can be used to remove duplicates. When unique molecular barcodes are included in the sequencing adapter, added during sample preparation, FGBio (see URLs) can be used to create consensus sequences of duplicates. BWA-MEM is the most used tool for sequence alignment of sequencing reads to the human reference genome. The output of BWA must be converted to BAM, sorted, and indexed using Samtools in order to be used in variant detection tools.
For SNVs and indels, there are many different open-source tools that can be used for the detection of these variant types. Somatic variants can be detected in somatic or tumor-only mode using the following tools: Strelka2, Freebayes, MuTect2, Pindel, BCFtools call, LoFreq, VarScan, Platypus, GATK4, and Scapel. Some tools such as GATK and BCFtools are designed to identify germline variants, whereas other tools, such as LoFreq and MuTect2, are designed to identify low-frequency variants common in tumor samples. A comparison of nine somatic variant calling programs found that Mutect2, Strelka2, and Virmid were among the most accurate.
There are also many different methods for the detection of copy number and structural variation. Some methods use sequence depth of coverage, or the average number of reads overlapping each region of the genome, to determine copy number changes. Because biases in coverage differ in each region of the genome, this coverage is not uniform, but can be normalized with healthy control samples. Additionally, some methods can take into account the allele frequency of common polymorphisms, called b-allele frequency, to correct these biases in the depth of coverage. A comparison of four CNV detection tools found that CNVKit had high sensitivity but a lower specificity relative to other programs like Control-FREEC. Although there are also many methods for the detection of structural variants, the accuracy of detection of structural variants is low compared with SNV and indels for short-read data. The strategy that most laboratories employ is the validation of known clinically actionable structural variants such as the FLT3 internal tandem duplication or known gene fusion events.
Computation interoperability and graphical interfaces
A software container is a freestanding unit that comprises the software of interest and all of its dependencies. Containerization allows for better maintenance and stability of computer software because it supports: (i) the deployment of the same code on the same running environment on any computer system (on premises or cloud); (ii) the separation of software packages to their own environment to satisfy specific running conditions; and (iii) dependencies and easier implementation of new packages or tools. Two popular containerization tools are docker, which was designed to work in a cloud environment, and singularity, which was designed to work in a local high-performance computing environment. Fortunately, containers that are created using docker can be converted to run using singularity on a local cluster. Additionally, containers can be versioned and easily shared, making them ideal to distribute code for bioinformatics pipelines.
Whether it is running in a high-performance computing environment locally or on the cloud, a bioinformatics pipeline has several important elements including automated advancement from step to step and parallelization. Pipelines are designed to take the input of each step as it becomes available and run blocks of code to convert it into predetermined output files. Some steps are serialized, meaning that they are dependent on the successful completion of earlier steps, such as alignment being dependent on the completion of read trimming. Other processes, such as variant calling using different callers, can be run concurrently to save time in a process called parallelization. Pipelines can be controlled on local computing resources using languages such as Nextflow, WDL, and Snakemate. In the cloud, commercial cloud frameworks exist to allow biologists with limited computational expertise to run pipelines in a point-and-click environment, using platforms such as DNANexus, Seven Bridges, and Illumina Connected Analytics.
| Approach|| |
To implement our pipelines, we (i) tested several software packages for analysis accuracy, (ii) created docker containers for each self-contained step in our workflow, (iii) implemented input options for customizing each step, and (iv) designed an end–end workflow for primary analysis that could be run on an internal high-performance computing cluster or run on the cloud using a point-and-click graphical interface.
We tested 10 variant calling methods, including Strelka2, Freebayes, MuTect2, Pindel, BCFtools call, LoFreq, VarScan, Platypus, GATK4, and Scapel, using data generated from an engineered cell line that was designed with variants with low allele frequency and validated by quantitative PCR ([Table S1]). Consistent with previous studies, we found high sensitivity with MuTect2 and Strelka2. Of the 21 SNVs and small indels (<50 bp) present, Freebayes detected all 21 variants, MuTect2 detected 20 and LoFreq detected 16; each caller was within 10% of the expected variants allele frequency (VAF). Callers such as Strelka2, GATK2, Platypus, and VarScan, using default parameters, detected the variants with >20% VAF, as designed ([Table S1 [Additional file 1]]). Pindel detected the 300 bp internal tandem duplicate (ITD) in the FLT3 gene, along with two other indels at 35%–40% VAF.
Tools implemented in docker and customization
Our workflow allows users to choose for SNV and indel detection Strelka2, Freebayes, MuTect2, and Pindel or a combination of the four. For MuTect2, alignments are first recalibrated using GATK4 with BaseRecalibrator and ApplyBQSR. If matched tumor/normal pairs are sequenced, users can ensure that these samples originate from the same patient using BCFtools and NGS Checkmate.
We implemented several tools for the detection of copy number and structural variants. For CNVs, we implemented CNVKit. Because CNVKit works best with a panel of healthy normal control samples, we have also implemented a container and tool to generate this healthy control reference. To detect ITDs, users can use Pindel and ITDSeeker., When microsatellite-specific baits are included, microsatellite stability can be estimated with an MSI-Sensor pro. Gene fusions (translocations) can be detected using DNA- and RNA-specific tools. Star-Fusion is implemented for RNASeq data, and DELLY and SVABA are implemented for DNA, sequencing data.
Because expression can also be assessed with the RNASeq data, we have implemented the steps necessary for assessing gene expression. Reads can be aligned with HiSAT2, and expression values are determined using FeatureCount and StringTie., Variants in the RNA can be determined using Freebayes and BamReadCt (see URLs).
Variants can be annotated using gnomAD for the detection of common mutations and snpEff for gene effect. Other sources of annotations include the database of oncoKB hotspots, Encode repeat regions, the database of non-synonymous of functional predictions (dbNSFP), and variant databases dbSNP, Clinvar, and COSMIC.
The bioinformatics workflow contains three elements: (i) a software container created using Docker, which contains all software dependencies for each step, (ii) scripts written in bash that contain software commands necessary to complete each step of the workflow, and (iii) the workflow script and configuration that defines the inputs and outputs of each step, the compute requirements, the bash script parameters, and the container used for each step. The workflow script and configuration were implemented for execution on a local high-performance computing cluster, using the workflow management program Nextflow, and on the cloud, using the DNANexus Toolkit.
For users with bioinformatics and computing expertise, Nextflow can be configured to run on a variety of platforms locally and on the cloud. Nextflow readily submits jobs to commonly used compute cluster scheduling software such as SGE, PBS, and SLURM but also can be configured to submit jobs to cloud systems on Amazon Web Services (AWS) and Google Cloud. Users can create a Nextflow configuration file to customize the workflow for their hardware. Additionally, Nextflow allows for users to resume failed jobs and has extensive logging of each step, making troubleshooting easy to document. Finally, these Nextflow workflows can be configured to run on individual tumor samples or tumor/normal sample pairs or in a batch mode for processing the data for an entire sequencing run.
In order to make these workflows accessible to users with limited computational expertise, we transformed the workflows to run on AWS resources on the DNANexus platform. DNANexus has a point-and-click user interface for running data analysis. Each step of the workflow was transformed into a DNANexus App, and apps were combined into a DNANexus Workflow. Users can run these pipelines with their raw sequence files in FastQ, a DNA reference tar gzip file and a gene panel reference tar gzip file. To reduce the price to run each step, every DNANexus app was run on test files to determine the minimal resources necessary, largely through monitoring memory and processor consumption and increasing resources incrementally when tools reached memory or disk usage limits ([Table S2 [Additional file 2]]). We then set these resource requirements as the default settings for each app. Users can alter these settings to decrease computing time. The cost of analysis is highly dependent on the size of the data set and the machines chosen to do the analysis, where the user will want to balance cost and computational time.
| Conclusion|| |
We have developed SCHOOL, a set of bioinformatics analysis tools and pipelines for the analysis of NGS data in an academic clinical oncology laboratory ([Figure S1 [Additional file 3]]), which has been in use at the CAP/CLIA laboratory at UT Southwestern Medical Center for four years. Additionally, SCHOOL pipelines have been used in over 20 research studies ranging from basic science to case reports.,,,, SCHOOL includes tools and methods for primary analysis of sequencing data from raw reads to finished variant calls, accommodating germline and somatic DNA and tumor RNA. We further include tools for panel-specific customization, including extensions for: copy number analysis, microsatellite stability, integration between DNA and RNA data, and structural variant calling to detect gene fusions in DNA and ITDs.
The SCHOOL pipeline can be optimized for the panel used in the assay. For example, when a panel of normal samples is included in assay development, a panel customization pipeline will align each sample and generate panel reference samples and the input BED files for copy number analysis. The aggregate normal sample VCF file will also be created to use when running MuTect2 to remove artifacts and rare variant sequences. Lastly, the normal samples can be used as a microsatellite reference for predicting microsatellite stability in the absence of a matched-normal sample.
For groups running an RNA Sequencing assay, there are several opportunities for data integrations, including comparison of RNA and DNA breakpoints in gene fusion events, comparison of splice site alteration using RNA data, and independent confirmation of variants by concurrent expression in RNA. The presence of variants in both the tumor RNA and tumor DNA provides enhanced confidence that a variant is not an artifact of the assay. Additionally, the presence of variants in RNA could indicate that the aberrant gene mutation is expressed in the tumor tissue. This could provide additional support of the importance of the variant, particularly for a suspected gain of function variants in oncogenic drivers or potential splice site variants that result in exon exclusion or intron inclusion. However, it is still important to note that variants in regulatory regions or gene deletions could prevent RNA expression.
These workflows represent a user-friendly, inexpensive, and flexible way to implement NGS bioinformatics for mutation detection and annotation in CAP/CLIA laboratories. SCHOOL can be easily installed on a local computing cluster that uses a queueing system such as SLURM or SGE with Nextflow with minimal dependency on pre-installed packages or runs on the cloud for laboratories that lack local resources and expertise with a user-friendly point-and-click graphic interface. We estimate the costs for analysis using cloud resources is a fraction of the costs of data generation and for some labs will reduce the need for additional computational expertise and resources on site.
These pipelines implemented in SCHOOL perform computationally heavy analysis in variant detection, which we consider to be the primary analysis. Users should, in the course of their validation studies, determine the filtering parameters of these results to maximize sensitivity and specificity for their assay using a set of quality metrics for the variants based on variant type, including a number of alternate reads, percent of alternate reads, strand bias, and other quality scores. In this secondary analysis step, in addition to filtering variants and removing artifacts, the user can determine mutational profiling metrics like tumor mutational burden and distributions of SNV by codon change. These secondary analysis steps often need tuning based on the assay and are less computationally intensive, meaning they can be done locally on a PC or laptop computer.
All the code used in SCHOOL is available at the UTSW Clinical Laboratory github wiki site: https://medforomics.github.io/schoolwiki/ using the repos: school for pipelines, process_scripts used in each docker container, and dnanexus_applets for running on DNANexus.
Financial support and sponsorship
Conflicts of interest
There are no conflicts of interest.
| References|| |
Kamps R, Brandão RD, Bosch BJvd, Paulussen ADC, Xanthoulea S, Blok MJ, et al
. Next-generation sequencing in oncology: Genetic diagnosis, risk prediction and cancer classification. Int J Mol Sci 2017;18:308.
Surrey LF, Luo M, Chang F, Li MM The genomic era of clinical oncology: Integrated genomic analysis for precision cancer care. Cytogenet Genome Res 2016;150:162-75.
Chang L, BChir MB, Chang M, Chang HM, Chang F Microsatellite instability: A predictive biomarker for cancer immunotherapy. Appl Immunohistochem Mol Morphol 2017;26:e15-21.
Salipante SJ, Scroggins SM, Hampel HL, Turner EH, Pritchard CC Microsatellite instability detection by next generation sequencing. Clin Chem 2014;60:1192-9.
Stenzinger A, Allen JD, Maas J, Stewart MD, Merino DM, Wempe MM, et al
. Tumor mutational burden standardization initiatives: Recommendations for consistent tumor mutational burden assessment in clinical samples to guide immunotherapy treatment decisions. Genes Chromosomes Cancer 2019;58:578-88.
SoRelle JA, Wachsmann M, Cantarel BL Assembling and validating bioinformatic pipelines for next-generation sequencing clinical assays. Arch Pathol Lab Med 2020;144:1118-30.
Ellrott K, Bailey MH, Saksena G, Covington KR, Kandoth C, Stewart C, et al
; MC3 Working Group; Cancer Genome Atlas Research Network. Scalable open science approach for mutation calling of tumor exomes using multiple genomic pipelines. Cell Syst 2018;6:271-281.e7.
Christensen PA, Ni Y, Bao F, Hendrickson HL, Greenwood M, Thomas JS, et al
. Houston methodist variant viewer: An application to support clinical laboratory interpretation of next-generation sequencing data for cancer. J Pathol Inform 2017;8:44.
Di Tommaso P, Chatzou M, Floden EW, Barja PP, Palumbo E, Notredame C Nextflow enables reproducible computational workflows. Nat Biotechnol 2017;35:316-9.
Zook JM, Catoe D, McDaniel J, Vang L, Spies N, Sidow A, et al
. Extensive sequencing of seven human genomes to characterize benchmark reference materials. Sci Data 2016;3:160025.
Li H Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv 2013:1-3.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al
; 1000 Genome Project Data Processing Subgroup. The sequence alignment/map format and Samtools. Bioinformatics 2009;25:2078-9.
Pan B, Kusko R, Xiao W, Zheng Y, Liu Z, Xiao C, et al
. Correction to: Similarities and differences between variants called with human reference genome HG19 or HG38. BMC Bioinformatics 2019;20:252.
Benjamin D, Sato T, Cibulskis K, Getz G, Stewart C, Lichtenstein L Calling somatic SNVs and Indels with Mutect2. bioRxiv 2019:1-8. doi:10.1101/861054.
DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, et al
. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet 2011;43:491-8.
Kim S, Scheffler K, Halpern AL, Bekritsky MA, Noh E, Källberg M, et al
. Strelka2: Fast and accurate calling of germline and somatic variants. Nat Methods 2018;15:591-4.
Garrison E, Marth G Haplotype-based variant detection from short-read sequencing. arXiv 2012:1-9.
Ye K, Schulz MH, Long Q, Apweiler R, Ning Z Pindel: A pattern growth approach to detect break points of large deletions and medium sized insertions from paired-end short reads. Bioinformatics 2009;25:2865-71.
Li H A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics 2011;27:2987-93.
Wilm A, Aw PP, Bertrand D, Yeo GH, Ong SH, Wong CH, et al
. Lofreq: A sequence-quality aware, ultra-sensitive variant caller for uncovering cell-population heterogeneity from high-throughput sequencing datasets. Nucleic Acids Res 2012;40:11189-201.
Koboldt DC, Chen K, Wylie T, Larson DE, McLellan MD, Mardis ER, et al
. Varscan: Variant detection in massively parallel sequencing of individual and pooled samples. Bioinformatics 2009;25:2283-5.
Rimmer A, Phan H, Mathieson I, Iqbal Z, Twigg SRF, Wilkie AOM, et al
; WGS500 Consortium. Integrating mapping-, assembly- and haplotype-based approaches for calling variants in clinical sequencing applications. Nat Genet 2014;46:912-8.
Fang H, Bergmann EA, Arora K, Vacic V, Zody MC, Iossifov I, et al
. Indel variant analysis of short-read sequencing data with scalpel. Nat Protoc 2016;11:2529-48.
Cingolani P, Platts A, Wang le L, Coon M, Nguyen T, Wang L, et al
. A program for annotating and predicting the effects of single nucleotide polymorphisms, SNPeff: SNPs in the genome of drosophila melanogaster strain w1118; iso-2; iso-3. Fly (Austin) 2012;6:80-92.
Karczewski KJ, Francioli LC, Tiao G, Cummings BB, Alföldi J, Wang Q, et al
; Genome Aggregation Database Consortium. The mutational constraint spectrum quantified from variation in 141,456 humans. Nature 2020;581:434-43.
Tate JG, Bamford S, Jubb HC, Sondka Z, Beare DM, Bindal N, et al
. COSMIC: The catalogue of somatic mutations in cancer. Nucleic Acids Res 2019;47:D941-7.
Rausch T, Zichner T, Schlattl A, Stütz AM, Benes V, Korbel JO DELLY: Structural variant discovery by integrated paired-end and split-read analysis. Bioinformatics 2012;28:i333-9.
Wala JA, Bandopadhayay P, Greenwald NF, O’Rourke R, Sharpe T, Stewart C, et al
. Svaba: Genome-wide detection of structural variants and indels by local assembly. Genome Res 2018;28:581-91.
Talevich E, Shain AH, Botton T, Bastian BC Cnvkit: Genome-wide copy number detection and visualization from targeted DNA sequencing. PloS Comput Biol 2016;12:e1004873.
Au CH, Wa A, Ho DN, Chan TL, Ma ES Clinical evaluation of panel testing by next-generation sequencing (NGS) for gene mutations in myeloid neoplasms. Diagn Pathol 2016;11:11.
Lee S, Lee J, Chae S, Moon Y, Lee HY, Park B, et al
. Multi-dimensional histone methylations for coordinated regulation of gene expression under hypoxia. Nucleic Acids Res 2017;45:11643-57.
Jia P, Yang X, Guo L, Liu B, Lin J, Liang H, et al
. Msisensor-pro: Fast, accurate, and matched-normal-sample-free detection of microsatellite instability. Genom Proteom Bioinform 2020;18:65-71.
Kim D, Paggi JM, Park C, Bennett C, Salzberg SL Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol 2019;37:907-15.
Liao Y, Smyth GK, Shi W Featurecounts: An efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 2014;30:923-30.
Kovaka S, Zimin AV, Pertea GM, Razaghi R, Salzberg SL, Pertea M Transcriptome assembly from long-read RNA-seq alignments with stringtie2. Genome Biol 2019;20:278.
Feng Y-Y, Cotto KC, Ramu A, Skidmore ZL, Kunisaki J, Richters M, et al
. RegTools: Integrated analysis of genomic and transcriptomic data for discovery of splicing variants in cancer. bioRxiv 2018:1-21. doi:10.1101/436634.
Mose LE, Perou CM, Parker JS Improved indel detection in DNA and RNA via realignment with ABRA2. Bioinformatics 2019;35:2966-73.
Lagunas-Rangel FA, Chávez-Valencia V FLT3–ITD and its current role in acute myeloid leukaemia. Med Oncol 2017;34:114.
Zaritsky A, Jamieson AR, Welf ES, Nevarez A, Cillay J, Eskiocak U, et al
. Interpretable deep learning uncovers cellular properties in label-free live cell images that are predictive of highly metastatic melanoma. Cell Syst 2021;12:733-47.e6.
Zhang W, Williams TA, Bhagwath AS, Hiermann JS, Peacock CD, Watkins CN, et al
. GEAMP, a novel gastroesophageal junction carcinoma cell line derived from a malignant pleural effusion. Lab Investig 2020;100:16-26.
Bishop JA, Gagan J, Krane JF, Jo VY Low-grade apocrine intraductal carcinoma: Expanding the morphologic and molecular spectrum of an enigmatic salivary gland tumor. Head Neck Pathol 2020;14:869-75.
Rooper LM, Agaimy A, Dickson BC, Dueber JC, Eberhart CG, Gagan J, et al
. DEK-AFF2 carcinoma of the sinonasal region and skull base: Detailed clinicopathologic characterization of a distinctive entity. Am J Surg Pathol 2021;45:1682-93.
Argani P, Palsgrove DN, Anders RA, Smith SC, Saoud C, Kwon R, et al
. A novel NIPBL-NACC1 gene fusion is characteristic of the cholangioblastic variant of intrahepatic cholangiocarcinoma. Am J Surg Pathol 2021;45:1550-60.