Experimental and Computational Approaches to Measure Telomere Length: Recent Advances and Future Directions.

PURPOSE OF REVIEW: The length of telomeres, protective structures at the chromosome ends, is a well-established biomarker for pathological conditions including multisystemic syndromes called telomere biology disorders. Approaches to measure telomere length (TL) differ on whether they estimate average, distribution, or chromosome-specific TL, and each presents their own advantages and limitations. RECENT FINDINGS: The development of long-read sequencing and publication of the telomere-to-telomere human genome reference has allowed for scalable and high-resolution TL estimation in pre-existing sequencing datasets but is still impractical as a dedicated TL test. As sequencing costs continue to fall and strategies for selectively enriching telomere regions prior to sequencing improve, these approaches may become a promising alternative to classic methods. Measurement methods rely on probe hybridization, qPCR or more recently, computational methods using sequencing data. Refinements of existing techniques and new approaches have been recently developed but a test that is accurate, simple, and scalable is still lacking.

CRISPR/Cas9-mediated excision of ALS/FTD-causing hexanucleotide repeat expansion in C9ORF72 rescues major disease mechanisms in vivo and in vitro.

A GGGGCC24+ hexanucleotide repeat expansion (HRE) in the C9ORF72 gene is the most common genetic cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD), fatal neurodegenerative diseases with no cure or approved treatments that substantially slow disease progression or extend survival. Mechanistic underpinnings of neuronal death include C9ORF72 haploinsufficiency, sequestration of RNA-binding proteins in the nucleus, and production of dipeptide repeat proteins. Here, we used an adeno-associated viral vector system to deliver CRISPR/Cas9 gene-editing machineries to effectuate the removal of the HRE from the C9ORF72 genomic locus. We demonstrate successful excision of the HRE in primary cortical neurons and brains of three mouse models containing the expansion (500-600 repeats) as well as in patient-derived iPSC motor neurons and brain organoids (450 repeats). This resulted in a reduction of RNA foci, poly-dipeptides and haploinsufficiency, major hallmarks of C9-ALS/FTD, making this a promising therapeutic approach to these diseases.

Telogator: a method for reporting chromosome-specific telomere lengths from long reads.

MOTIVATION: Telomeres are the repetitive sequences found at the ends of eukaryotic chromosomes and are often thought of as a 'biological clock,' with their average length shortening during division in most cells. In addition to their association with senescence, abnormal telomere lengths are well known to be associated with multiple cancers, short telomere syndromes and as risk factors for a broad range of diseases. While a majority of methods for measuring telomere length will report average lengths across all chromosomes, it is known that aberrations in specific chromosome arms are biomarkers for certain diseases. Due to their repetitive nature, characterizing telomeres at this resolution is prohibitive for short read sequencing approaches, and is challenging still even with longer reads. RESULTS: We present Telogator: a method for reporting chromosome-specific telomere length from long read sequencing data. We demonstrate Telogator's sensitivity in detecting chromosome-specific telomere length in simulated data across a range of read lengths and error rates. Telogator is then applied to 10 germline samples, yielding a high correlation with short read methods in reporting average telomere length. In addition, we investigate common subtelomere rearrangements and identify the minimum read length required to anchor telomere/subtelomere boundaries in samples with these haplotypes. AVAILABILITY AND IMPLEMENTATION: Telogator is written in Python3 and is available at github.com/zstephens/telogator. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Combining copy number, methylation markers, and mutations as a panel for endometrial cancer detection via intravaginal tampon collection.

OBJECTIVE: We aimed to assess whether endometrial cancer (EC) can be detected in shed DNA collected with vaginal tampon by analyzing copy number, methylation markers, and mutations. METHODS: Tampons were collected prior to hysterectomy from 38 EC patients and 28 women with benign indications. Extracted tampon DNA underwent the following: 1) low-coverage whole genome sequencing (LC-WGS) to assess copy number, 2) pyrosequencing to measure percent promotor methylation of HOXA9, RASSF1, and CDH13 and 3) next generation sequencing (NGS) to identify mutations in 19 genes associated with EC identified through The Cancer Genome Atlas. Sensitivity and specificity for each test and test combinations were calculated. RESULTS: Methylation analysis yielded the highest specificities but lowest sensitivities (37-40% sensitivity; 100% specificity for HOXA9, RASSF1 and HTR1B) while mutation analysis had improved sensitivity (50% sensitivity; 83% specificity). Only one "false positive" result for copy number variants was identified among women with benign surgical indications, which was based on detection of copy number changes, and associated with a leiomyosarcoma that was only recognized at hysterectomy. Considering any of the 3 biomarker classes as a positive, resulted in a sensitivity of 92% and specificity of 86%. Mutation analysis did not add sensitivity to the combination of analysis of copy number and methylation. CONCLUSIONS: This study demonstrates a proof-of-principle for non-invasive yet precise detection of endometrial cancer. We propose that with improved biomarker testing, it may be possible to develop a clinically useful test for detecting EC.

Prognostic subclass of intrahepatic cholangiocarcinoma by integrative molecular-clinical analysis and potential targeted approach.

BACKGROUND: Although molecular characterization of iCCA has been studied recently, integrative analysis of molecular and clinical characterization has not been fully established. If molecular features of iCCA can be predicted based on clinical findings, we can approach to distinguish targeted treatment. We analyzed RNA sequencing data annotated with clinicopathologic data to clarify molecular-specific clinical features and to evaluate potential therapies for molecular subtypes. METHODS: We performed next-generation RNA sequencing of 30 surgically resected iCCA from Korean patients and the clinicopathologic features were analyzed. The RNA sequences from 32 iCCA resected from US patients were used for validation. RESULTS: Patients were grouped into two subclasses on the basis of unsupervised clustering, which showed a difference in 5-year survival rates (48.5% vs 14.2%, p = 0.007) and similar survival outcome in the US samples. In subclass B (poor prognosis), both data sets were similar in higher carcinoembryonic antigen and cancer antigen 19-9 levels, underlying cholangitis, and bile duct-type pathology; in subclass A (better prognosis), there was more frequent viral hepatitis and cholangiolar-type pathology. On pathway analysis, subclass A had enriched liver-related signatures. Subclass B had enriched inflammation-related and TP53 pathways, with more frequent KRAS mutations. CCA cell lines with similar gene expression patterns of subclass A were sensitive to gemcitabine. CONCLUSIONS: Two molecular subtypes of iCCA with distinct clinicopathological differences were identified. Knowledge of clinical and pathologic characteristics can predict molecular subtypes, and knowledge of different subtype signaling pathways may lead to more rational, targeted approaches to treatment.

Ccne1 Overexpression Causes Chromosome Instability in Liver Cells and Liver Tumor Development in Mice.

BACKGROUND & AIMS: The CCNE1 locus, which encodes cyclin E1, is amplified in many types of cancer cells and is activated in hepatocellular carcinomas (HCCs) from patients infected with hepatitis B virus or adeno-associated virus type 2, due to integration of the virus nearby. We investigated cell-cycle and oncogenic effects of cyclin E1 overexpression in tissues of mice. METHODS: We generated mice with doxycycline-inducible expression of Ccne1 (Ccne1T mice) and activated overexpression of cyclin E1 from age 3 weeks onward. At 14 months of age, livers were collected from mice that overexpress cyclin E1 and nontransgenic mice (controls) and analyzed for tumor burden and by histology. Mouse embryonic fibroblasts (MEFs) and hepatocytes from Ccne1T and control mice were analyzed to determine the extent to which cyclin E1 overexpression perturbs S-phase entry, DNA replication, and numbers and structures of chromosomes. Tissues from 4-month-old Ccne1T and control mice (at that age were free of tumors) were analyzed for chromosome alterations, to investigate the mechanisms by which cyclin E1 predisposes hepatocytes to transformation. RESULTS: Ccne1T mice developed more hepatocellular adenomas and HCCs than control mice. Tumors developed only in livers of Ccne1T mice, despite high levels of cyclin E1 in other tissues. Ccne1T MEFs had defects that promoted chromosome missegregation and aneuploidy, including incomplete replication of DNA, centrosome amplification, and formation of nonperpendicular mitotic spindles. Whereas Ccne1T mice accumulated near-diploid aneuploid cells in multiple tissues and organs, polyploidization was observed only in hepatocytes, with losses and gains of whole chromosomes, DNA damage, and oxidative stress. CONCLUSIONS: Livers, but not other tissues of mice with inducible overexpression of cyclin E1, develop tumors. More hepatocytes from the cyclin E1-overexpressing mice were polyploid than from control mice, and had losses or gains of whole chromosomes, DNA damage, and oxidative stress; all of these have been observed in human HCC cells. The increased risk of HCC in patients with hepatitis B virus or adeno-associated virus type 2 infection might involve activation of cyclin E1 and its effects on chromosomes and genomes of liver cells.

A data science approach for the classification of low-grade and high-grade ovarian serous carcinomas.

BACKGROUND: Copy Number Alternations (CNAs) is defined as somatic gain or loss of DNA regions. The profiles of CNAs may provide a fingerprint specific to a tumor type or tumor grade. Low-coverage sequencing for reporting CNAs has recently gained interest since successfully translated into clinical applications. Ovarian serous carcinomas can be classified into two largely mutually exclusive grades, low grade and high grade, based on their histologic features. The grade classification based on the genomics may provide valuable clue on how to best manage these patients in clinic. Based on the study of ovarian serous carcinomas, we explore the methodology of combining CNAs reporting from low-coverage sequencing with machine learning techniques to stratify tumor biospecimens of different grades. RESULTS: We have developed a data-driven methodology for tumor classification using the profiles of CNAs reported by low-coverage sequencing. The proposed method called Bag-of-Segments is used to summarize fixed-length CNA features predictive of tumor grades. These features are further processed by machine learning techniques to obtain classification models. High accuracy is obtained for classifying ovarian serous carcinoma into high and low grades based on leave-one-out cross-validation experiments. The models that are weakly influenced by the sequence coverage and the purity of the sample can also be built, which would be of higher relevance for clinical applications. The patterns captured by Bag-of-Segments features correlate with current clinical knowledge: low grade ovarian tumors being related to aneuploidy events associated to mitotic errors while high grade ovarian tumors are induced by DNA repair gene malfunction. CONCLUSIONS: The proposed data-driven method obtains high accuracy with various parametrizations for the ovarian serous carcinoma study, indicating that it has good generalization potential towards other CNA classification problems. This method could be applied to the more difficult task of classifying ovarian serous carcinomas with ambiguous histology or in those with low grade tumor co-existing with high grade tumor. The closer genomic relationship of these tumor samples to low or high grade may provide important clinical value.

Indel sensitive and comprehensive variant/mutation detection from RNA sequencing data for precision medicine.

BACKGROUND: RNA-seq is the most commonly used sequencing application. Not only does it measure gene expression but it is also an excellent media to detect important structural variants such as single nucleotide variants (SNVs), insertion/deletion (Indels) or fusion transcripts. However, detection of these variants is challenging and complex from RNA-seq. Here we describe a sensitive and accurate analytical pipeline which detects various mutations at once for translational precision medicine. METHODS: The pipeline incorporates most sensitive aligners for Indels in RNA-Seq, the best practice for data preprocessing and variant calling, and STAR-fusion is for chimeric transcripts. Variants/mutations are annotated, and key genes can be extracted for further investigation and clinical actions. Three datasets were used to evaluate the performance of the pipeline for SNVs, indels and fusion transcripts. RESULTS: For the well-defined variants from NA12878 by GIAB project, about 95% and 80% of sensitivities were obtained for SNVs and indels, respectively, in matching RNA-seq. Comparison with other variant specific tools showed good performance of the pipeline. For the lung cancer dataset with 41 known and oncogenic mutations, 39 were detected by the pipeline with STAR aligner and all by the GSNAP aligner. An actionable EML4 and ALK fusion was also detected in one of the tumors, which also demonstrated outlier ALK expression. For 9 fusions spiked-into RNA-seq libraries with different concentrations, the pipeline was able to detect all in unfiltered results although some at very low concentrations may be missed when filtering was applied. CONCLUSIONS: The new RNA-seq workflow is an accurate and comprehensive mutation profiler from RNA-seq. Key or actionable mutations are reliably detected from RNA-seq, which makes it a practical alternative source for personalized medicine.

HGT-ID: an efficient and sensitive workflow to detect human-viral insertion sites using next-generation sequencing data.

BACKGROUND: Transfer of genetic material from microbes or viruses into the host genome is known as horizontal gene transfer (HGT). The integration of viruses into the human genome is associated with multiple cancers, and these can now be detected using next-generation sequencing methods such as whole genome sequencing and RNA-sequencing. RESULTS: We designed a novel computational workflow, HGT-ID, to identify the integration of viruses into the human genome using the sequencing data. The HGT-ID workflow primarily follows a four-step procedure: i) pre-processing of unaligned reads, ii) virus detection using subtraction approach, iii) identification of virus integration site using discordant and soft-clipped reads and iv) HGT candidates prioritization through a scoring function. Annotation and visualization of the events, as well as primer design for experimental validation, are also provided in the final report. We evaluated the tool performance with the well-understood cervical cancer samples. The HGT-ID workflow accurately detected known human papillomavirus (HPV) integration sites with high sensitivity and specificity compared to previous HGT methods. We applied HGT-ID to The Cancer Genome Atlas (TCGA) whole-genome sequencing data (WGS) from liver tumor-normal pairs. Multiple hepatitis B virus (HBV) integration sites were identified in TCGA liver samples and confirmed by HGT-ID using the RNA-Seq data from the matched liver pairs. This shows the applicability of the method in both the data types and cross-validation of the HGT events in liver samples. We also processed 220 breast tumor WGS data through the workflow; however, there were no HGT events detected in those samples. CONCLUSIONS: HGT-ID is a novel computational workflow to detect the integration of viruses in the human genome using the sequencing data. It is fast and accurate with functions such as prioritization, annotation, visualization and primer design for future validation of HGTs. The HGT-ID workflow is released under the MIT License and available at http://kalarikrlab.org/Software/HGT-ID.html .

Correction to: the impact of pharmacokinetic gene profiles across human cancers.

Correction to: BMC Cancer (2018) 18:577 DOI https://doi.org/10.1186/s12885-018-4345-2.

The impact of pharmacokinetic gene profiles across human cancers.

BACKGROUND: The right drug to the right patient at the right time is one of the ideals of Individualized Medicine (IM) and remains one of the most compelling promises of the post-genomic age. The addition of genomic information is expected to increase the precision of an individual patient's treatment, resulting in improved outcomes. While pilot studies have been encouraging, key aspects of interpreting tumor genomics information, such as somatic activation of drug transport or metabolism, have not been systematically evaluated. METHODS: In this work, we developed a simple rule-based approach to classify the therapies administered to each patient from The Cancer Genome Atlas PanCancer dataset (n = 2858) as effective or ineffective. Our Therapy Efficacy model used each patient's drug target and pharmacokinetic (PK) gene expression profile; the specific genes considered for each patient depended on the therapies they received. Patients who received predictably ineffective therapies were considered at high-risk of cancer-related mortality and those who did not receive ineffective therapies were considered at low-risk. The utility of our Therapy Efficacy model was assessed using per-cancer and pan-cancer differential survival. RESULTS: Our simple rule-based Therapy Efficacy model classified 143 (5%) patients as high-risk. High-risk patients had age ranges comparable to low-risk patients of the same cancer type and tended to be later stage and higher grade (odds ratios of 1.6 and 1.4, respectively). A significant pan-cancer association was identified between predictions of our Therapy Efficacy model and poorer overall survival (hazard ratio, HR = 1.47, p = 6.3 × 10- 3). Individually, drug export (HR = 1.49, p = 4.70 × 10- 3) and drug metabolism (HR = 1.73, p = 9.30 × 10- 5) genes demonstrated significant survival associations. Survival associations for target gene expression are mechanism-dependent. Similar results were observed for event-free survival. CONCLUSIONS: While the resolution of clinical information within the dataset is not ideal, and modeling the relative contribution of each gene to the activity of each therapy remains a challenge, our approach demonstrates that somatic PK alterations should be integrated into the interpretation of somatic transcriptomic profiles as they likely have a significant impact on the survival of specific patients. We believe that this approach will aid the prospective design of personalized therapeutic strategies.

Prevalent Homozygous Deletions of Type I Interferon and Defensin Genes in Human Cancers Associate with Immunotherapy Resistance.

Purpose: Homozygous deletions play important roles in carcinogenesis. The genome-wide screening for homozygously deleted genes in many different cancer types with a large number of patient specimens representing the tumor heterogeneity has not been done.Experimental Design: We performed integrative analyses of the copy-number profiles of 10,759 patients across 31 cancer types from The Cancer Genome Atlas project.Results: We found that the type-I interferon, α-, and ß-defensin genes were homozygously deleted in 19 cancer types with high frequencies (7%-31%, median = 12%; interquartile range = 10%-16.5%). Patients with homozygous deletion of interferons exhibited significantly shortened overall or disease-free survival time in a number of cancer types, whereas patients with homozygous deletion of defensins did not significantly associate with worse overall or disease-free survival. Gene expression analyses suggested that homozygous deletion of interferon and defensin genes could activate genes involved in oncogenic and cell-cycle pathways but repress other genes involved in immune response pathways, suggesting their roles in promoting tumorigenesis and helping cancer cells evade immune surveillance. Further analysis of the whole exomes of 109 patients with melanoma demonstrated that the homozygous deletion of interferon (P = 0.0029, OR = 11.8) and defensin (P = 0.06, OR = 2.79) genes are significantly associated with resistance to anti-CTLA4 immunotherapy.Conclusions: Our analysis reveals that the homozygous deletion of interferon and defensin genes is prevalent in human cancers, and importantly this feature can be used as a novel prognostic biomarker for immunotherapy resistance. Clin Cancer Res; 24(14); 3299-308. ©2018 AACR.

Early genetic aberrations in patients with sporadic colorectal cancer.

Chromosome instability (CIN) is widely observed in both sporadic and hereditary colorectal cancer (CRC). Defects in APC and WNT signaling are primarily associated with CIN in hereditary CRC, but the genetic causes for CIN in sporadic CRC remain elusive. Using high-density SNP array and exome data from The Cancer Genome Atlas (TCGA), we characterized loss of heterozygosity (LOH) and copy number variation (CNV) in the peripheral blood, normal colon, and corresponding tumor tissue in 15 CRC patients with proficient mismatch repair (MMR) and 24 CRC patients with deficient MMR. We found a high frequency of 18q LOH in tumors and arm-specific enrichment of genetic aberrations on 18q in the normal colon (primarily copy neutral LOH) and blood (primarily copy gain). These aberrations were specific to the sporadic, pMMR CRC. Though in tumor samples genetic aberrations were observed for genes commonly mutated in hereditary CRC (eg, APC, CTNNB1, SMAD4, BRAF), none of them showed LOH or CNV in the normal colon or blood. DCC located on 18q21.1 topped the list of genes with genetic aberrations in the tumor. In an independent cohort of 13 patients subjected to Whole Genome Sequencing (WGS), we found LOH and CNV on 18q in adenomatous polyp and tumor tissues. Our data suggests that patients with sporadic CRC may have genetic aberrations preferentially enriched on 18q in their blood, normal colon epithelium, and non-malignant polyp lesions that may prove useful as a clinical marker for sporadic CRC detection and risk assessment.

Comprehensive Genomic Profiling of a Rare Thyroid Follicular Dendritic Cell Sarcoma.

We previously reported an extremely rare case of follicular dendritic cell sarcoma (FDCS) presented as a thyroid mass. Given the rarity of this disease, there are no personalized and molecularly targeted treatment options due to the lack of knowledge in the genomic makeup of the tumor. A 44-year-old white woman was diagnosed with an extranodal FDCS in thyroid. The patient underwent a total thyroidectomy, central compartment dissection, parathyroid re-implantation, and adjuvant radiation therapy. Tumor DNA sequencing of 236 genes by FoundationOne panel found truncating mutations in PTEN and missense mutations in RET and TP53. However, patient-matched germline DNA was not sequenced which is critical for identification of true somatic mutations. Furthermore, the FoundationOne panel doesn't measure genomic rearrangements which have been shown to be abundant in sarcomas and are associated with sarcoma tumorigenesis and progression. In the current study, we carried out comprehensive genomic sequencing of the tumor, adjacent normal tissues, and patient-matched blood, in an effort to understand the genomic makeup of this rare extranodal FDCS and to identify potential therapeutic targets. Eighty-one somatic point mutations were identified in tumor but not in adjacent normal tissues or blood. A clonal truncating mutation in the CLTCL1 gene, which stabilizes the mitotic spindle, was likely a driver mutation of tumorigenesis and could explain the extensive copy number aberrations (CNAs) and genomic rearrangements in the tumor including a chr15/chr17 local chromothripsis resulted in 6 expressed fusion genes. The fusion gene HDGFRP3→SHC4 led to a 200-fold increase in the expression of oncogene SHC4 which is a potential target of the commercial drug Dasatinib. Missense mutations in ATM and splice-site mutation in VEGFR1 were also detected in addition to the TP53 missense mutation reported by FoundationOne.

Experience with precision genomics and tumor board, indicates frequent target identification, but barriers to delivery.

BACKGROUND: The ability to analyze the genomics of malignancies has opened up new possibilities for off-label targeted therapy in cancers that are refractory to standard therapy. At Mayo Clinic these efforts are organized through the Center for Individualized Medicine (CIM). RESULTS: Prior to GTB, datasets were analyzed and integrated by a team of bioinformaticians and cancer biologists. Therapeutically actionable mutations were identified in 65% (92/141) of the patients tested with 32% (29/92) receiving genomically targeted therapy with FDA approved drugs or in an independent clinical trial with 45% (13/29) responding. Standard of care (SOC) options were continued by 15% (14/92) of patients tested before exhausting SOC options, with 71% (10/14) responding to treatment. Over 35% (34/92) of patients with actionable targets were not treated with 65% (22/34) choosing comfort measures or passing away. MATERIALS AND METHODS: Patients (N = 165) were referred to the CIM Clinic between October 2012 and December 2015. All patients received clinical genomic panel testing with selected subsets receiving array comparative genomic hybridization and clinical whole exome sequencing to complement and validate panel findings. A genomic tumor board (GTB) reviewed results and, when possible, developed treatment recommendations. CONCLUSIONS: Treatment decisions driven by tumor genomic analysis can lead to significant clinical benefit in a minority of patients. The success of genomically driven therapy depends both on access to drugs and robustness of bioinformatics analysis. While novel clinical trial designs are increasing the utility of genomic testing, robust data sharing of outcomes is needed to optimize clinical benefit for all patients.

Indel detection from RNA-seq data: tool evaluation and strategies for accurate detection of actionable mutations.

Driver somatic mutations are a hallmark of a tumor that can be used for diagnosis and targeted therapy. Mutations are primarily detected from tumor DNA. As dynamic molecules of gene activities, transcriptome profiling by RNA sequence (RNA-seq) is becoming increasingly popular, which not only measures gene expression but also structural variations such as mutations and fusion transcripts. Although single-nucleotide variants (SNVs) can be easily identified from RNA-seq, intermediate long insertions/deletions (indels > 2 bases and less than sequence reads) cause significant challenges and are ignored by most RNA-seq analysis tools. This study evaluates commonly used RNA-seq analysis programs along with variant and somatic mutation callers in a series of data sets with simulated and known indels. The aim is to develop strategies for accurate indel detection. Our results show that the RNA-seq alignment is the most important step for indel identification and the evaluated programs have a wide range of sensitivity to map sequence reads with indels, from not at all to decently sensitive. The sensitivity is impacted by sequence read lengths. Most variant calling programs rely on hard evidence indels marked in the alignment and the programs with realignment may use soft-clipped reads for indel inferencing. Based on the observations, we have provided practical recommendations for indel detection when different RNA-seq aligners are used and demonstrated the best option with highly reliable results. With careful customization of bioinformatics algorithms, RNA-seq can be reliably used for both SNV and indel mutation detection that can be used for clinical decision-making.

Molecular Modeling and Functional Analysis of Exome Sequencing-Derived Variants of Unknown Significance Identify a Novel, Constitutively Active FGFR2 Mutant in Cholangiocarcinoma.

PURPOSE: Genomic testing has increased the quantity of information available to oncologists. Unfortunately, many identified sequence alterations are variants of unknown significance (VUSs), which thus limit the clinician's ability to use these findings to inform treatment. We applied a combination of in silico prediction and molecular modeling tools and laboratory techniques to rapidly define actionable VUSs. MATERIALS AND METHODS: Exome sequencing was conducted on 308 tumors from various origins. Most single nucleotide alterations within gene coding regions were VUSs. These VUSs were filtered to identify a subset of therapeutically targetable genes that were predicted with in silico tools to be altered in function by their variant sequence. A subset of receptor tyrosine kinase VUSs was characterized by laboratory comparison of each VUS versus its wild-type counterpart in terms of expression and signaling activity. RESULTS: The study identified 4,327 point mutations of which 3,833 were VUSs. Filtering for mutations in genes that were therapeutically targetable and predicted to affect protein function reduced these to 522VUSs of interest, including a large number of kinases. Ten receptortyrosine kinase VUSs were selected to explore in the laboratory. Of these, seven were found to be functionally altered. Three VUSs (FGFR2 F276C, FGFR4 R78H, and KDR G539R) showed increased basal or ligand-stimulated ERK phosphorylation compared with their wild-type counterparts, which suggests that they support transformation. Treatment of a patient who carried FGFR2 F276C with an FGFR inhibitor resulted in significant and sustained tumor response with clinical benefit. CONCLUSION: The findings demonstrate the feasibility of rapid identification of the biologic relevance of somatic mutations, which thus advances clinicians' ability to make informed treatment decisions.

Measure transcript integrity using RNA-seq data.

BACKGROUND: Stored biological samples with pathology information and medical records are invaluable resources for translational medical research. However, RNAs extracted from the archived clinical tissues are often substantially degraded. RNA degradation distorts the RNA-seq read coverage in a gene-specific manner, and has profound influences on whole-genome gene expression profiling. RESULT: We developed the transcript integrity number (TIN) to measure RNA degradation. When applied to 3 independent RNA-seq datasets, we demonstrated TIN is a reliable and sensitive measure of the RNA degradation at both transcript and sample level. Through comparing 10 prostate cancer clinical samples with lower RNA integrity to 10 samples with higher RNA quality, we demonstrated that calibrating gene expression counts with TIN scores could effectively neutralize RNA degradation effects by reducing false positives and recovering biologically meaningful pathways. When further evaluating the performance of TIN correction using spike-in transcripts in RNA-seq data generated from the Sequencing Quality Control consortium, we found TIN adjustment had better control of false positives and false negatives (sensitivity = 0.89, specificity = 0.91, accuracy = 0.90), as compared to gene expression analysis results without TIN correction (sensitivity = 0.98, specificity = 0.50, accuracy = 0.86). CONCLUSION: TIN is a reliable measurement of RNA integrity and a valuable approach used to neutralize in vitro RNA degradation effect and improve differential gene expression analysis.

PANDA: pathway and annotation explorer for visualizing and interpreting gene-centric data.

Objective. Bringing together genomics, transcriptomics, proteomics, and other -omics technologies is an important step towards developing highly personalized medicine. However, instrumentation has advances far beyond expectations and now we are able to generate data faster than it can be interpreted. Materials and Methods. We have developed PANDA (Pathway AND Annotation) Explorer, a visualization tool that integrates gene-level annotation in the context of biological pathways to help interpret complex data from disparate sources. PANDA is a web-based application that displays data in the context of well-studied pathways like KEGG, BioCarta, and PharmGKB. PANDA represents data/annotations as icons in the graph while maintaining the other data elements (i.e., other columns for the table of annotations). Custom pathways from underrepresented diseases can be imported when existing data sources are inadequate. PANDA also allows sharing annotations among collaborators. Results. In our first use case, we show how easy it is to view supplemental data from a manuscript in the context of a user's own data. Another use-case is provided describing how PANDA was leveraged to design a treatment strategy from the somatic variants found in the tumor of a patient with metastatic sarcomatoid renal cell carcinoma. Conclusion. PANDA facilitates the interpretation of gene-centric annotations by visually integrating this information with context of biological pathways. The application can be downloaded or used directly from our website: http://bioinformaticstools.mayo.edu/research/panda-viewer/.