Somatic inactivating PTPRJ mutations and dysregulated pathways identified in canine malignant melanoma by integrated comparative genomic analysis.

Canine malignant melanoma, a significant cause of mortality in domestic dogs, is a powerful comparative model for human melanoma, but little is known about its genetic etiology. We mapped the genomic landscape of canine melanoma through multi-platform analysis of 37 tumors (31 mucosal, 3 acral, 2 cutaneous, and 1 uveal) and 17 matching constitutional samples including long- and short-insert whole genome sequencing, RNA sequencing, array comparative genomic hybridization, single nucleotide polymorphism array, and targeted Sanger sequencing analyses. We identified novel predominantly truncating mutations in the putative tumor suppressor gene PTPRJ in 19% of cases. No BRAF mutations were detected, but activating RAS mutations (24% of cases) occurred in conserved hotspots in all cutaneous and acral and 13% of mucosal subtypes. MDM2 amplifications (24%) and TP53 mutations (19%) were mutually exclusive. Additional low-frequency recurrent alterations were observed amidst low point mutation rates, an absence of ultraviolet light mutational signatures, and an abundance of copy number and structural alterations. Mutations that modulate cell proliferation and cell cycle control were common and highlight therapeutic axes such as MEK and MDM2 inhibition. This mutational landscape resembles that seen in BRAF wild-type and sun-shielded human melanoma subtypes. Overall, these data inform biological comparisons between canine and human melanoma while suggesting actionable targets in both species.

Long insert whole genome sequencing for copy number variant and translocation detection.

As next-generation sequencing continues to have an expanding presence in the clinic, the identification of the most cost-effective and robust strategy for identifying copy number changes and translocations in tumor genomes is needed. We hypothesized that performing shallow whole genome sequencing (WGS) of 900-1000-bp inserts (long insert WGS, LI-WGS) improves our ability to detect these events, compared with shallow WGS of 300-400-bp inserts. A priori analyses show that LI-WGS requires less sequencing compared with short insert WGS to achieve a target physical coverage, and that LI-WGS requires less sequence coverage to detect a heterozygous event with a power of 0.99. We thus developed an LI-WGS library preparation protocol based off of Illumina's WGS library preparation protocol and illustrate the feasibility of performing LI-WGS. We additionally applied LI-WGS to three separate tumor/normal DNA pairs collected from patients diagnosed with different cancers to demonstrate our application of LI-WGS on actual patient samples for identification of somatic copy number alterations and translocations. With the evolution of sequencing technologies and bioinformatics analyses, we show that modifications to current approaches may improve our ability to interrogate cancer genomes.

Germline mutations in HOXB13 and prostate-cancer risk.

BACKGROUND: Family history is a significant risk factor for prostate cancer, although the molecular basis for this association is poorly understood. Linkage studies have implicated chromosome 17q21-22 as a possible location of a prostate-cancer susceptibility gene. METHODS: We screened more than 200 genes in the 17q21-22 region by sequencing germline DNA from 94 unrelated patients with prostate cancer from families selected for linkage to the candidate region. We tested family members, additional case subjects, and control subjects to characterize the frequency of the identified mutations. RESULTS: Probands from four families were discovered to have a rare but recurrent mutation (G84E) in HOXB13 (rs138213197), a homeobox transcription factor gene that is important in prostate development. All 18 men with prostate cancer and available DNA in these four families carried the mutation. The carrier rate of the G84E mutation was increased by a factor of approximately 20 in 5083 unrelated subjects of European descent who had prostate cancer, with the mutation found in 72 subjects (1.4%), as compared with 1 in 1401 control subjects (0.1%) (P=8.5x10(-7)). The mutation was significantly more common in men with early-onset, familial prostate cancer (3.1%) than in those with late-onset, nonfamilial prostate cancer (0.6%) (P=2.0x10(-6)). CONCLUSIONS: The novel HOXB13 G84E variant is associated with a significantly increased risk of hereditary prostate cancer. Although the variant accounts for a small fraction of all prostate cancers, this finding has implications for prostate-cancer risk assessment and may provide new mechanistic insights into this common cancer. (Funded by the National Institutes of Health and others.).

Identification of extracellular miRNA in human cerebrospinal fluid by next-generation sequencing.

There has been a growing interest in using next-generation sequencing (NGS) to profile extracellular small RNAs from the blood and cerebrospinal fluid (CSF) of patients with neurological diseases, CNS tumors, or traumatic brain injury for biomarker discovery. Small sample volumes and samples with low RNA abundance create challenges for downstream small RNA sequencing assays. Plasma, serum, and CSF contain low amounts of total RNA, of which small RNAs make up a fraction. The purpose of this study was to maximize RNA isolation from RNA-limited samples and apply these methods to profile the miRNA in human CSF by small RNA deep sequencing. We systematically tested RNA isolation efficiency using ten commercially available kits and compared their performance on human plasma samples. We used RiboGreen to quantify total RNA yield and custom TaqMan assays to determine the efficiency of small RNA isolation for each of the kits. We significantly increased the recovery of small RNA by repeating the aqueous extraction during the phenol-chloroform purification in the top performing kits. We subsequently used the methods with the highest small RNA yield to purify RNA from CSF and serum samples from the same individual. We then prepared small RNA sequencing libraries using Illumina's TruSeq sample preparation kit and sequenced the samples on the HiSeq 2000. Not surprisingly, we found that the miRNA expression profile of CSF is substantially different from that of serum. To our knowledge, this is the first time that the small RNA fraction from CSF has been profiled using next-generation sequencing.

Open-access synthetic spike-in mRNA-seq data for cancer gene fusions.

BACKGROUND: Oncogenic fusion genes underlie the mechanism of several common cancers. Next-generation sequencing based RNA-seq analyses have revealed an increasing number of recurrent fusions in a variety of cancers. However, absence of a publicly available gene-fusion focused RNA-seq data impedes comparative assessment and collaborative development of novel gene fusions detection algorithms. We have generated nine synthetic poly-adenylated RNA transcripts that correspond to previously reported oncogenic gene fusions. These synthetic RNAs were spiked at known molarity over a wide range into total RNA prior to construction of next-generation sequencing mRNA libraries to generate RNA-seq data. RESULTS: Leveraging a priori knowledge about replicates and molarity of each synthetic fusion transcript, we demonstrate utility of this dataset to compare multiple gene fusion algorithms' detection ability. In general, more fusions are detected at higher molarity, indicating that our constructs performed as expected. However, systematic detection differences are observed based on molarity or algorithm-specific characteristics. Fusion-sequence specific detection differences indicate that for applications where specific sequences are being investigated, additional constructs may be added to provide quantitative data that is specific for the sequence of interest. CONCLUSIONS: To our knowledge, this is the first publicly available synthetic RNA-seq data that specifically leverages known cancer gene-fusions. The proposed method of designing multiple gene-fusion constructs over a wide range of molarity allows granular performance analyses of multiple fusion-detection algorithms. The community can leverage and augment this publicly available data to further collaborative development of analytical tools and performance assessment frameworks for gene fusions from next-generation sequencing data.

Whole-genome sequencing of multiple myeloma from diagnosis to plasma cell leukemia reveals genomic initiating events, evolution, and clonal tides.

The longitudinal evolution of a myeloma genome from diagnosis to plasma cell leukemia has not previously been reported. We used whole-genome sequencing (WGS) on 4 purified tumor samples and patient germline DNA drawn over a 5-year period in a t(4;14) multiple myeloma patient. Tumor samples were acquired at diagnosis, first relapse, second relapse, and end-stage secondary plasma cell leukemia (sPCL). In addition to the t(4;14), all tumor time points also shared 10 common single-nucleotide variants (SNVs) on WGS comprising shared initiating events. Interestingly, we observed genomic sequence variants that waxed and waned with time in progressive tumors, suggesting the presence of multiple independent, yet related, clones at diagnosis that rose and fell in dominance. Five newly acquired SNVs, including truncating mutations of RB1 and ZKSCAN3, were observed only in the final sPCL sample suggesting leukemic transformation events. This longitudinal WGS characterization of the natural history of a high-risk myeloma patient demonstrated tumor heterogeneity at diagnosis with shifting dominance of tumor clones over time and has also identified potential mutations contributing to myelomagenesis as well as transformation from myeloma to overt extramedullary disease such as sPCL.

Effect of selection of QTc formula on eligibility of cancer patients for phase I clinical trials.

A retrospective analysis of 130 patients was conducted in a Phase I oncology clinic to assess the effect of QTc formula selection on clinical trial eligibility. QTc values were calculated from screening electrocardiograms using 7 formulae (Bazett, Fridericia, Framingham, Hodges, Mayeda, Van de Water and Wohlfart). QTc values > 470 ms for females and > 450 ms for males were used to define prolongation. Concomitant medication potential for QTc prolongation was determined using a public database (AzCert). Ineligibility rates ranged from 3.1 % to 17.7 % (Framingham: 3.1 %, Van de Water: 3.1 %, Hodges: 3.1 %, Wohlfart: 3.1 %, Fridericia: 3.9 %, Bazett: 10.8 % and Mayeda: 17.7 %). A consistent ineligibility rate was achieved by using formulae-specific thresholds. Fifty one percent of patients were taking concomitant medications with QTc prolongation potential. The proportion of concomitant medications with the potential to prolong QTc was 11.57 % (96 of 830). Uniform criteria and guidelines for selection of QTc formulae need to be developed. Formulae-specific QTc thresholds also need to be specified.

Next-generation sequencing of Coccidioides immitis isolated during cluster investigation.

Next-generation sequencing enables use of whole-genome sequence typing (WGST) as a viable and discriminatory tool for genotyping and molecular epidemiologic analysis. We used WGST to confirm the linkage of a cluster of Coccidioides immitis isolates from 3 patients who received organ transplants from a single donor who later had positive test results for coccidioidomycosis. Isolates from the 3 patients were nearly genetically identical (a total of 3 single-nucleotide polymorphisms identified among them), thereby demonstrating direct descent of the 3 isolates from an original isolate. We used WGST to demonstrate the genotypic relatedness of C. immitis isolates that were also epidemiologically linked. Thus, WGST offers unique benefits to public health for investigation of clusters considered to be linked to a single source.

GRM7 variants confer susceptibility to age-related hearing impairment.

Age-related hearing impairment (ARHI), or presbycusis, is the most prevalent sensory impairment in the elderly. ARHI is a complex disease caused by an interaction between environmental and genetic factors. Here we describe the results of the first whole genome association study for ARHI. The study was performed using 846 cases and 846 controls selected from 3434 individuals collected by eight centers in six European countries. DNA pools for cases and controls were allelotyped on the Affymetrix 500K GeneChip for each center separately. The 252 top-ranked single nucleotide polymorphisms (SNPs) identified in a non-Finnish European sample group (1332 samples) and the 177 top-ranked SNPs from a Finnish sample group (360 samples) were confirmed using individual genotyping. Subsequently, the 23 most interesting SNPs were individually genotyped in an independent European replication group (138 samples). This resulted in the identification of a highly significant and replicated SNP located in GRM7, the gene encoding metabotropic glutamate receptor type 7. Also in the Finnish sample group, two GRM7 SNPs were significant, albeit in a different region of the gene. As the Finnish are genetically distinct from the rest of the European population, this may be due to allelic heterogeneity. We performed histochemical studies in human and mouse and showed that mGluR7 is expressed in hair cells and in spiral ganglion cells of the inner ear. Together these data indicate that common alleles of GRM7 contribute to an individual's risk of developing ARHI, possibly through a mechanism of altered susceptibility to glutamate excitotoxicity.

G-SQZ: compact encoding of genomic sequence and quality data.

SUMMARY: Large volumes of data generated by high-throughput sequencing instruments present non-trivial challenges in data storage, content access and transfer. We present G-SQZ, a Huffman coding-based sequencing-reads-specific representation scheme that compresses data without altering the relative order. G-SQZ has achieved from 65% to 81% compression on benchmark datasets, and it allows selective access without scanning and decoding from start. This article focuses on describing the underlying encoding scheme and its software implementation, and a more theoretical problem of optimal compression is out of scope. The immediate practical benefits include reduced infrastructure and informatics costs in managing and analyzing large sequencing data. AVAILABILITY: http://public.tgen.org/sqz. Academic/non-profit: SOURCE: available at no cost under a non-open-source license by requesting from the web-site; Binary: available for direct download at no cost. For-Profit: Submit request for for-profit license from the web-site.

Resolving individuals contributing trace amounts of DNA to highly complex mixtures using high-density SNP genotyping microarrays.

We use high-density single nucleotide polymorphism (SNP) genotyping microarrays to demonstrate the ability to accurately and robustly determine whether individuals are in a complex genomic DNA mixture. We first develop a theoretical framework for detecting an individual's presence within a mixture, then show, through simulations, the limits associated with our method, and finally demonstrate experimentally the identification of the presence of genomic DNA of specific individuals within a series of highly complex genomic mixtures, including mixtures where an individual contributes less than 0.1% of the total genomic DNA. These findings shift the perceived utility of SNPs for identifying individual trace contributors within a forensics mixture, and suggest future research efforts into assessing the viability of previously sub-optimal DNA sources due to sample contamination. These findings also suggest that composite statistics across cohorts, such as allele frequency or genotype counts, do not mask identity within genome-wide association studies. The implications of these findings are discussed.

Multimarker analysis and imputation of multiple platform pooling-based genome-wide association studies.

For many genome-wide association (GWA) studies individually genotyping one million or more SNPs provides a marginal increase in coverage at a substantial cost. Much of the information gained is redundant due to the correlation structure inherent in the human genome. Pooling-based GWA studies could benefit significantly by utilizing this redundancy to reduce noise, improve the accuracy of the observations and increase genomic coverage. We introduce a measure of correlation between individual genotyping and pooling, under the same framework that r(2) provides a measure of linkage disequilibrium (LD) between pairs of SNPs. We then report a new non-haplotype multimarker multi-loci method that leverages the correlation structure between SNPs in the human genome to increase the efficacy of pooling-based GWA studies. We first give a theoretical framework and derivation of our multimarker method. Next, we evaluate simulations using this multimarker approach in comparison to single marker analysis. Finally, we experimentally evaluate our method using different pools of HapMap individuals on the Illumina 450S Duo, Illumina 550K and Affymetrix 5.0 platforms for a combined total of 1 333 631 SNPs. Our results show that use of multimarker analysis reduces noise specific to pooling-based studies, allows for efficient integration of multiple microarray platforms and provides more accurate measures of significance than single marker analysis. Additionally, this approach can be extended to allow for imputing the association significance for SNPs not directly observed using neighboring SNPs in LD. This multimarker method can now be used to cost-effectively complete pooling-based GWA studies with multiple platforms across over one million SNPs and to impute neighboring SNPs weighted for the loss of information due to pooling.

Oligonucleotide fingerprint identification for microarray-based pathogen diagnostic assays.

MOTIVATION: Advances in DNA microarray technology and computational methods have unlocked new opportunities to identify 'DNA fingerprints', i.e. oligonucleotide sequences that uniquely identify a specific genome. We present an integrated approach for the computational identification of DNA fingerprints for design of microarray-based pathogen diagnostic assays. We provide a quantifiable definition of a DNA fingerprint stated both from a computational as well as an experimental point of view, and the analytical proof that all in silico fingerprints satisfying the stated definition are found using our approach. RESULTS: The presented computational approach is implemented in an integrated high-performance computing (HPC) software tool for oligonucleotide fingerprint identification termed TOFI. We employed TOFI to identify in silico DNA fingerprints for several bacteria and plasmid sequences, which were then experimentally evaluated as potential probes for microarray-based diagnostic assays. Results and analysis of approximately 150 in silico DNA fingerprints for Yersinia pestis and 250 fingerprints for Francisella tularensis are presented. AVAILABILITY: The implemented algorithm is available upon request.

SNiPer-HD: improved genotype calling accuracy by an expectation-maximization algorithm for high-density SNP arrays.

MOTIVATION: The technology to genotype single nucleotide polymorphisms (SNPs) at extremely high densities provides for hypothesis-free genome-wide scans for common polymorphisms associated with complex disease. However, we find that some errors introduced by commonly employed genotyping algorithms may lead to inflation of false associations between markers and phenotype. RESULTS: We have developed a novel SNP genotype calling program, SNiPer-High Density (SNiPer-HD), for highly accurate genotype calling across hundreds of thousands of SNPs. The program employs an expectation-maximization (EM) algorithm with parameters based on a training sample set. The algorithm choice allows for highly accurate genotyping for most SNPs. Also, we introduce a quality control metric for each assayed SNP, such that poor-behaving SNPs can be filtered using a metric correlating to genotype class separation in the calling algorithm. SNiPer-HD is superior to the standard dynamic modeling algorithm and is complementary and non-redundant to other algorithms, such as BRLMM. Implementing multiple algorithms together may provide highly accurate genotyping calls, without inflation of false positives due to systematically miss-called SNPs. A reliable and accurate set of SNP genotypes for increasingly dense panels will eliminate some false association signals and false negative signals, allowing for rapid identification of disease susceptibility loci for complex traits. AVAILABILITY: SNiPer-HD is available at TGen's website: http://www.tgen.org/neurogenomics/data.

A somatic reference standard for cancer genome sequencing.

Large-scale multiplexed identification of somatic alterations in cancer has become feasible with next generation sequencing (NGS). However, calibration of NGS somatic analysis tools has been hampered by a lack of tumor/normal reference standards. We thus performed paired PCR-free whole genome sequencing of a matched metastatic melanoma cell line (COLO829) and normal across three lineages and across separate institutions, with independent library preparations, sequencing, and analysis. We generated mean mapped coverages of 99X for COLO829 and 103X for the paired normal across three institutions. Results were combined with previously generated data allowing for comparison to a fourth lineage on earlier NGS technology. Aggregate variant detection led to the identification of consensus variants, including key events that represent hallmark mutation types including amplified BRAF V600E, a CDK2NA small deletion, a 12 kb PTEN deletion, and a dinucleotide TERT promoter substitution. Overall, common events include >35,000 point mutations, 446 small insertion/deletions, and >6,000 genes affected by copy number changes. We present this reference to the community as an initial standard for enabling quantitative evaluation of somatic mutation pipelines across institutions.

Genome-wide analysis of hepatic lipid content in extreme obesity.

AIMS: Individuals with type 2 diabetes have an increased risk of developing non-alcoholic fatty liver disease (NAFLD), and NAFLD patients are also at greater risk for developing type 2 diabetes. Although the relationship between type 2 diabetes and NAFLD is highly interconnected, the pathogenic mechanisms linking the two diseases are poorly understood. The goal of this study was to identify genetic determinants of hepatic lipid accumulation through association analysis using histological phenotypes in obese individuals. METHODS: Using the Illumina HumanOmniExpress BeadChip assay, we genotyped 2,300 individuals on whom liver biopsy data were available. RESULTS: We analyzed total bilirubin levels, which are linked to fatty liver in severe obesity, and observed the strongest evidence for association with rs4148325 in UGT1A (P < 5.0 × 10(-93)), replicating previous findings. We assessed hepatic fat level and found strong evidence for association with rs4823173, rs2896019, and rs2281135, all located in PNPLA3 and rs10401969 in SUGP1. Analysis of liver transcript levels of 20 genes residing at the SUGP1/NCAN locus identified a 1.6-fold change in the expression of the LPAR2 gene in fatty liver. We also observed suggestive evidence for association between low-grade fat accumulation and rs10859525 and rs1294908, located upstream from SOCS2 and RAMP3, respectively. SOCS2 was differentially expressed between fatty and normal liver. CONCLUSIONS: These results replicate findings for several hepatic phenotypes in the setting of extreme obesity and implicate new loci that may play a role in the pathophysiology of hepatic lipid accumulation.

Personalized treatment of Sézary syndrome by targeting a novel CTLA4:CD28 fusion.

Matching molecularly targeted therapies with cancer subtype-specific gene mutations is revolutionizing oncology care. However, for rare cancers this approach is problematic due to the often poor understanding of the disease's natural history and phenotypic heterogeneity, making treatment of these cancers a particularly unmet medical need in clinical oncology. Advanced Sézary syndrome (SS), an aggressive, exceedingly rare variant of cutaneous T-cell lymphoma (CTCL) is a prototypical example of a rare cancer. Through whole genome and RNA sequencing (RNA-seq) of a SS patient's tumor we discovered a highly expressed gene fusion between CTLA4 (cytotoxic T lymphocyte antigen 4) and CD28 (cluster of differentiation 28), predicting a novel stimulatory molecule on the surface of tumor T cells. Treatment with the CTLA4 inhibitor ipilimumab resulted in a rapid clinical response. Our findings suggest a novel driver mechanism for SS, and cancer in general, and exemplify an emerging model of cancer treatment using exploratory genomic analysis to identify a personally targeted treatment option when conventional therapies are exhausted.

Toward precision medicine in glioblastoma: the promise and the challenges.

Integrated sequencing strategies have provided a broader understanding of the genomic landscape and molecular classifications of multiple cancer types and have identified various therapeutic opportunities across cancer subsets. Despite pivotal advances in the characterization of genomic alterations in glioblastoma, targeted agents have shown minimal efficacy in clinical trials to date, and patient survival remains poor. In this review, we highlight potential reasons why targeting single alterations has yielded limited clinical efficacy in glioblastoma, focusing on issues of tumor heterogeneity and pharmacokinetic failure. We outline strategies to address these challenges in applying precision medicine to glioblastoma and the rationale for applying targeted combination therapy approaches that match genomic alterations with compounds accessible to the central nervous system.

Detection of an atypical teratoid rhabdoid brain tumor gene deletion in circulating blood using next-generation sequencing.

Circulating biomarkers such as somatic chromosome mutations are novel diagnostic tools to detect cancer noninvasively. We describe focal deletions found in a patient with atypical teratoid rhabdoid tumor, a highly aggressive early childhood pediatric tumor. First, we used magnetic resonance imaging (MRI) and histopathology to study the tumor anatomy. Next, we used whole genome sequencing (Next Gen Sequencing) and Bioinformatics interrogation to discover the presence of 3 focal deletions in tumor tissue and 2 of these 3 focal deletions in patient's blood also. About 20% of the blood DNA sequencing reads matched the tumor DNA reads at the SMARCB1 gene locus. Circulating, tumor-specific DNA aberrations are a promising biomarker for atypical teratoid rhabdoid tumor patients. The high percentage of tumor DNA detected in blood indicates that either circulating brain tumor cells lyse in the blood or that contents of brain tumor cells traverse a possibly compromised blood-brain barrier in this patient.

Profiles of extracellular miRNA in cerebrospinal fluid and serum from patients with Alzheimer's and Parkinson's diseases correlate with disease status and features of pathology.

The discovery and reliable detection of markers for neurodegenerative diseases have been complicated by the inaccessibility of the diseased tissue--such as the inability to biopsy or test tissue from the central nervous system directly. RNAs originating from hard to access tissues, such as neurons within the brain and spinal cord, have the potential to get to the periphery where they can be detected non-invasively. The formation and extracellular release of microvesicles and RNA binding proteins have been found to carry RNA from cells of the central nervous system to the periphery and protect the RNA from degradation. Extracellular miRNAs detectable in peripheral circulation can provide information about cellular changes associated with human health and disease. In order to associate miRNA signals present in cell-free peripheral biofluids with neurodegenerative disease status of patients with Alzheimer's and Parkinson's diseases, we assessed the miRNA content in cerebrospinal fluid and serum from postmortem subjects with full neuropathology evaluations. We profiled the miRNA content from 69 patients with Alzheimer's disease, 67 with Parkinson's disease and 78 neurologically normal controls using next generation small RNA sequencing (NGS). We report the average abundance of each detected miRNA in cerebrospinal fluid and in serum and describe 13 novel miRNAs that were identified. We correlated changes in miRNA expression with aspects of disease severity such as Braak stage, dementia status, plaque and tangle densities, and the presence and severity of Lewy body pathology. Many of the differentially expressed miRNAs detected in peripheral cell-free cerebrospinal fluid and serum were previously reported in the literature to be deregulated in brain tissue from patients with neurodegenerative disease. These data indicate that extracellular miRNAs detectable in the cerebrospinal fluid and serum are reflective of cell-based changes in pathology and can be used to assess disease progression and therapeutic efficacy.