Surgical Tumor Resection Deregulates Hallmarks of Cancer in Resected Tissue and the Surrounding Microenvironment.

Surgery exposes tumor tissue to severe hypoxia and mechanical stress leading to rapid gene expression changes in the tumor and its microenvironment, which remain poorly characterized. We biopsied tumor and adjacent normal tissues from patients with breast (n = 81) and head/neck squamous cancers (HNSC; n = 10) at the beginning (A), during (B), and end of surgery (C). Tumor/normal RNA from 46/81 patients with breast cancer was subjected to mRNA-Seq using Illumina short-read technology, and from nine patients with HNSC to whole-transcriptome microarray with Illumina BeadArray. Pathways and genes involved in 7 of 10 known cancer hallmarks, namely, tumor-promoting inflammation (TNF-A, NFK-B, IL18 pathways), activation of invasion and migration (various extracellular matrix-related pathways, cell migration), sustained proliferative signaling (K-Ras Signaling), evasion of growth suppressors (P53 signaling, regulation of cell death), deregulating cellular energetics (response to lipid, secreted factors, and adipogenesis), inducing angiogenesis (hypoxia signaling, myogenesis), and avoiding immune destruction (CTLA4 and PDL1) were significantly deregulated during surgical resection (time points A vs. B vs. C). These findings were validated using NanoString assays in independent pre/intra/post-operative breast cancer samples from 48 patients. In a comparison of gene expression data from biopsy (analogous to time point A) with surgical resection samples (analogous to time point C) from The Cancer Genome Atlas study, the top deregulated genes were the same as identified in our analysis, in five of the seven studied cancer types. This study suggests that surgical extirpation deregulates the hallmarks of cancer in primary tumors and adjacent normal tissue across different cancers. IMPLICATIONS: Surgery deregulates hallmarks of cancer in human tissue.

Building flexible and robust analysis frameworks for molecular subtyping of cancers.

Molecular subtyping is essential to infer tumor aggressiveness and predict prognosis. In practice, tumor profiling requires in-depth knowledge of bioinformatics tools involved in the processing and analysis of the generated data. Additionally, data incompatibility (e.g., microarray versus RNA sequencing data) and technical and uncharacterized biological variance between training and test data can pose challenges in classifying individual samples. In this article, we provide a roadmap for implementing bioinformatics frameworks for molecular profiling of human cancers in a clinical diagnostic setting. We describe a framework for integrating several methods for quality control, normalization, batch correction, classification and reporting, and develop a use case of the framework in breast cancer.

The mutational landscape of early- and typical-onset oral tongue squamous cell carcinoma.

BACKGROUND: The incidence of oral tongue squamous cell carcinoma (OTSCC) is increasing among younger birth cohorts. The etiology of early-onset OTSCC (diagnosed before the age of 50 years) and cancer driver genes remain largely unknown. METHODS: The Sequencing Consortium of Oral Tongue Cancer was established through the pooling of somatic mutation data of oral tongue cancer specimens (n = 227 [107 early-onset cases]) from 7 studies and The Cancer Genome Atlas. Somatic mutations at microsatellite loci and Catalog of Somatic Mutations in Cancer mutation signatures were identified. Cancer driver genes were identified with the MutSigCV and WITER algorithms. Mutation comparisons between early- and typical-onset OTSCC were evaluated via linear regression with adjustments for patient-related factors. RESULTS: Two novel driver genes (ATXN1 and CDC42EP1) and 5 previously reported driver genes (TP53, CDKN2A, CASP8, NOTCH1, and FAT1) were identified. Six recurrent mutations were identified, with 4 occurring in TP53. Early-onset OTSCC had significantly fewer nonsilent mutations even after adjustments for tobacco use. No associations of microsatellite locus mutations and mutation signatures with the age of OTSCC onset were observed. CONCLUSIONS: This international, multicenter consortium is the largest study to characterize the somatic mutational landscape of OTSCC and the first to suggest differences by age of onset. This study validates multiple previously identified OTSCC driver genes and proposes 2 novel cancer driver genes. In analyses by age, early-onset OTSCC had a significantly smaller somatic mutational burden that was not explained by differences in tobacco use. LAY SUMMARY: This study identifies 7 specific areas in the human genetic code that could be responsible for promoting the development of tongue cancer. Tongue cancer in young patients (under the age of 50 years) has fewer overall changes to the genetic code in comparison with tongue cancer in older patients, but the authors do not think that this is due to differences in smoking rates between the 2 groups. The cause of increasing cases of tongue cancer in young patients remains unclear.

Functional genomics screen with pooled shRNA library and gene expression profiling with extracts of Azadirachta indica identify potential pathways for therapeutic targets in head and neck squamous cell carcinoma.

Tumor suppression by the extracts of Azadirachta indica (neem) works via anti-proliferation, cell cycle arrest, and apoptosis, demonstrated previously using cancer cell lines and live animal models. However, very little is known about the molecular targets and pathways that neem extracts and their associated compounds act through. Here, we address this using a genome-wide functional pooled shRNA screen on head and neck squamous cell carcinoma cell lines treated with crude neem leaf extracts, known for their anti-tumorigenic activity. We analyzed differences in global clonal sizes of the shRNA-infected cells cultured under no treatment and treatment with neem leaf extract conditions, assayed using next-generation sequencing. We found 225 genes affected the cancer cell growth in the shRNA-infected cells treated with neem extract. Pathway enrichment analyses of whole-genome gene expression data from cells temporally treated with neem extract revealed important roles played by the TGF-ß pathway and HSF-1-related gene network. Our results indicate that neem extract affects various important molecular signaling pathways in head and neck cancer cells, some of which may be therapeutic targets for this devastating tumor.

Detection of High-Risk Human Papillomavirus in Oral Cavity Squamous Cell Carcinoma Using Multiple Analytes and Their Role in Patient Survival.

PURPOSE: Accurate detection of human papillomavirus (HPV) in oral cavity squamous cell carcinoma (OSCC) is essential to understanding the role of HPV in disease prognosis and management of patients. We used different analytes and methods to understand the true prevalence of HPV in a cohort of patients with OSCC with different molecular backgrounds, and we correlated HPV data with patient survival. METHODS: We integrated data from multiple analytes (HPV DNA, HPV RNA, and p16), assays (immunohistochemistry, polymerase chain reaction [PCR], quantitative PCR [qPCR], and digital PCR), and molecular changes (somatic mutations and DNA methylation) from 153 patients with OSCC to correlate p16 expression, HPV DNA, and HPV RNA with HPV incidence and patient survival. RESULTS: High prevalence (33% to 58%) of HPV16/18 DNA did not correlate with the presence of transcriptionally active viral genomes (15%) in tumors. Eighteen percent of the tumors were p16 positive and only 6% were both HPV DNA and HPV RNA positive. Most tumors with relatively high copy number HPV DNA and/or HPV RNA, but not with HPV DNA alone (irrespective of copy number), were wild-type for TP53 and CASP8 genes. In our study, p16 protein, HPV DNA, and HPV RNA, either alone or in combination, did not correlate with patient survival. Nine HPV-associated genes stratified the virus-positive from the virus-negative tumor group with high confidence ( P < .008) when HPV DNA copy number and/or HPV RNA were considered to define HPV positivity, and not HPV DNA alone, irrespective of copy number ( P < .2). CONCLUSION: In OSCC, the presence of both HPV RNA and p16 is rare. HPV DNA alone is not an accurate measure of HPV positivity and therefore may not be informative. HPV DNA, HPV RNA, and p16 do not correlate with patients' outcome.

A minimal set of internal control genes for gene expression studies in head and neck squamous cell carcinoma.

Selection of the right reference gene(s) is crucial in the analysis and interpretation of gene expression data. The aim of the present study was to discover and validate a minimal set of internal control genes in head and neck tumor studies. We analyzed data from multiple sources (in house whole-genome gene expression microarrays, previously published quantitative real-time PCR (qPCR) data and RNA-seq data from TCGA) to come up with a list of 18 genes (discovery set) that had the lowest variance, a high level of expression across tumors, and their matched normal samples. The genes in the discovery set were ranked using four different algorithms (BestKeeper, geNorm, NormFinder, and comparative delta Ct) and a web-based comparative tool, RefFinder, for their stability and variance in expression across tissues. Finally, we validated their expression using qPCR in an additional set of tumor:matched normal samples that resulted in five genes (RPL30, RPL27, PSMC5, MTCH1, and OAZ1), out of which RPL30 and RPL27 were most stable and were abundantly expressed across the tissues. Our data suggest that RPL30 or RPL27 in combination with either PSMC5 or MTCH1 or OAZ1 can be used as a minimal set of control genes in head and neck tumor gene expression studies.

CAFE MOCHA: An Integrated Platform for Discovering Clinically Relevant Molecular Changes in Cancer-An Example of Distant Metastasis- and Recurrence-Linked Classifiers in Head and Neck Squamous Cell Carcinoma.

PURPOSE: With large amounts of multidimensional molecular data on cancers generated and deposited into public repositories such as The Cancer Genome Atlas and International Cancer Genome Consortium, a cancer type agnostic and integrative platform will help to identify signatures with clinical relevance. We devised such a platform and showcase it by identifying a molecular signature for patients with metastatic and recurrent (MR) head and neck squamous cell carcinoma (HNSCC). METHODS: We devised a statistical framework accompanied by a graphical user interface-driven application, Clinical Association of Functionally Established MOlecular CHAnges ( CAFE MOCHA; https://github.com/binaypanda/CAFEMOCHA), to discover molecular signatures linked to a specific clinical attribute in a cancer type. The platform integrates mutations and indels, gene expression, DNA methylation, and copy number variations to discover a classifier first and then to predict an incoming tumor for the same by pulling defined class variables into a single framework that incorporates a coordinate geometry-based algorithm called complete specificity margin-based clustering, which ensures maximum specificity. CAFE MOCHA classifies an incoming tumor sample using either its matched normal or a built-in database of normal tissues. The application is packed and deployed using the install4j multiplatform installer. We tested CAFE MOCHA in HNSCC tumors (n = 513) followed by validation in tumors from an independent cohort (n = 18) for discovering a signature linked to distant MR. RESULTS: CAFE MOCHA identified an integrated signature, MR44, associated with distant MR HNSCC, with 80% sensitivity and 100% specificity in the discovery stage and 100% sensitivity and 100% specificity in the validation stage. CONCLUSION: CAFE MOCHA is a cancer type and clinical attribute agnostic statistical framework to discover integrated molecular signatures.

Comparative analyses of putative toxin gene homologs from an Old World viper, Daboia russelii.

Availability of snake genome sequences has opened up exciting areas of research on comparative genomics and gene diversity. One of the challenges in studying snake genomes is the acquisition of biological material from live animals, especially from the venomous ones, making the process cumbersome and time-consuming. Here, we report comparative sequence analyses of putative toxin gene homologs from Russell's viper (Daboia russelii) using whole-genome sequencing data obtained from shed skin. When compared with the major venom proteins in Russell's viper studied previously, we found 45-100% sequence similarity between the venom proteins and their putative homologs in the skin. Additionally, comparative analyses of 20 putative toxin gene family homologs provided evidence of unique sequence motifs in nerve growth factor (NGF), platelet derived growth factor (PDGF), Kunitz/Bovine pancreatic trypsin inhibitor (Kunitz BPTI), cysteine-rich secretory proteins, antigen 5, andpathogenesis-related1 proteins (CAP) and cysteine-rich secretory protein (CRISP). In those derived proteins, we identified V11 and T35 in the NGF domain; F23 and A29 in the PDGF domain; N69, K2 and A5 in the CAP domain; and Q17 in the CRISP domain to be responsible for differences in the largest pockets across the protein domain structures in crotalines, viperines and elapids from the in silico structure-based analysis. Similarly, residues F10, Y11 and E20 appear to play an important role in the protein structures across the kunitz protein domain of viperids and elapids. Our study highlights the usefulness of shed skin in obtaining good quality high-molecular weight DNA for comparative genomic studies, and provides evidence towards the unique features and evolution of putative venom gene homologs in vipers.

Sample preparation method considerations for integrated transcriptomic and proteomic analysis of tumors.

Sample processing protocols that enable compatible recovery of differentially expressed transcripts and proteins are necessary for integration of the multiomics data applied in the analysis of tumors. In this pilot study, we compared two different isolation methods for extracting RNA and protein from laryngopharyngeal tumor tissues and the corresponding adjacent normal sections. In Method 1, RNA and protein were isolated from a single tissue section sequentially and in Method 2, the extraction was carried out using two different sections and two independent and parallel protocols for RNA and protein. RNA and protein from both methods were subjected to RNA-seq and iTRAQ-based LC-MS/MS analysis, respectively. Analysis of data revealed that a higher number of differentially expressed transcripts and proteins were concordant in their regulation trends in Method 1 as compared to Method 2. Cross-method comparison of concordant entities revealed that RNA and protein extraction from the same tissue section (Method 1) recovered more concordant entities that are missed in the other extraction method (Method 2) indicating heterogeneity in distribution of these entities in different tissue sections. Method 1 could thus be the method of choice for integrated analysis of transcriptome and proteome data.

Comparative proteomics of the superior and inferior spikelets at the early grain filling stage in rice cultivars contrast for panicle compactness and ethylene evolution.

The breeding programmes in rice aimed at increasing the number of spikelets per panicle have been accompanied by poor grain filling in the inferior spikelets of large panicle rice, leading to yield disadvantage. The present study attempted to understand the reason for differential grain filling in the inferior and superior spikelets by comparative proteomics considering a compact-panicle rice cultivar Mahalaxmi and a lax-panicle rice cultivar Upahar, which show poor and good grain filling, respectively. An initial study of two rice cultivars for panicle compactness and grain filling revealed an inverse correlation between the two parameters. It was further observed that the panicle compactness in Mahalaxmi was associated with a higher evolution of ethylene by the spikelets, both superior and inferior, compared with the lax-panicle Upahar. The proteomic studies revealed that the superior and inferior spikelets of Mahalaxmi differentially expressed 21 proteins that were also expressed in Upahar. However, in Upahar, only two of these proteins were differentially expressed between the superior and inferior spikelets, indicating that the metabolic activities of the spikelets occupying the superior and inferior positions on the panicle were very different in Mahalaxmi compared with those in Upahar. Among the proteins that were downregulated in the inferior spikelets compared with the superior ones in Mahalaxmi were importin-α, elongation factor 1-ß and cell division control protein 48, which are essential for cell cycle progression and cell division. Low expression of these proteins might inhibit endosperm cell division in the inferior spikelets, limiting their sink capacity and leading to poor grain filling compared to that in the superior spikelets. The poor grain filling in Mahalaxmi may also be a result of the high evolution of ethylene in the inferior spikelets, which is reflected from the observation that these spikelets showed significantly higher expression of S-adenosylmethionine synthase and the gene encoding the enzyme than the superior spikelets in this cultivar, but not in Upahar; S-adenosynlmethionine synthase catalyses the synthesis of S-adenosylmethionine, the precursor of ethylene biosynthesis.

A Minimal DNA Methylation Signature in Oral Tongue Squamous Cell Carcinoma Links Altered Methylation with Tumor Attributes.

UNLABELLED: Oral tongue squamous cell carcinomas (OTSCC) are a homogenous group of aggressive tumors in the head and neck region that spread early to lymph nodes and have a higher incidence of regional failure. In addition, there is a rising incidence of oral tongue cancer in younger populations. Studies on functional DNA methylation changes linked with altered gene expression are critical for understanding the mechanisms underlying tumor development and metastasis. Such studies also provide important insight into biomarkers linked with viral infection, tumor metastasis, and patient survival in OTSCC. Therefore, we performed genome-wide methylation analysis of tumors (N = 52) and correlated altered methylation with differential gene expression. The minimal tumor-specific DNA 5-methylcytosine signature identified genes near 16 different differentially methylated regions, which were validated using genomic data from The Cancer Genome Atlas cohort. In our cohort, hypermethylation of MIR10B was significantly associated with the differential expression of its target genes NR4A3 and BCL2L11 (P = 0.0125 and P = 0.014, respectively), which was inversely correlated with disease-free survival (P = 9E-15 and P = 2E-15, respectively) in patients. Finally, differential methylation in FUT3, TRIM5, TSPAN7, MAP3K8, RPS6KA2, SLC9A9, and NPAS3 genes was found to be predictive of certain clinical and epidemiologic parameters. IMPLICATIONS: This study reveals a functional minimal methylation profile in oral tongue tumors with associated risk habits, clinical, and epidemiologic outcomes. In addition, NR4A3 downregulation and correlation with patient survival suggests a potential target for therapeutic intervention in oral tongue tumors. Data from the current study are deposited in the NCBI Geo database (accession number GSE75540). Mol Cancer Res; 14(9); 805-19. ©2016 AACR.

Identification and Characterization of Differentially Expressed Genes in Inferior and Superior Spikelets of Rice Cultivars with Contrasting Panicle-Compactness and Grain-Filling Properties.

Breeding programs for increasing spikelet number in rice have resulted in compactness of the panicle, accompanied by poor grain filling in inferior spikelets. Although the inefficient utilization of assimilate has been indicated as responsible for this poor grain filling, the underlying cause remains elusive. The current study utilized the suppression subtractive hybridization technique to identify 57 and 79 genes that overexpressed in the superior and inferior spikelets (with respect to each other), respectively, of the compact-panicle rice cultivar Mahalaxmi. Functional categorization of these differentially expressed genes revealed a marked metabolic difference between the spikelets according to their spatial location on the panicle. The expression of genes encoding seed storage proteins was dominant in inferior spikelets, whereas genes encoding regulatory proteins, such as serine-threonine kinase, zinc finger protein and E3 ligase, were highly expressed in superior spikelets. The expression patterns of these genes in the inferior and superior spikelets of Mahalaxmi were similar to those observed in another compact-panicle cultivar, OR-1918, but differed from those obtained in two lax-panicle cultivars, Upahar and Lalat. The results first suggest that the regulatory proteins abundantly expressed in the superior spikelets of compact-panicle cultivars and in both the superior and inferior spikelets of lax-panicle cultivars but poorly expressed in the inferior spikelets of compact-panicle cultivars promote grain filling. Second, the high expression of seed-storage proteins observed in the inferior spikelets of compact-panicle cultivars appears to inhibit the grain filling process. Third, the low expression of enzymes of the Krebs cycle in inferior spikelets compared with superior spikelets of compact-panicle cultivars is bound to lead to poor ATP generation in the former and consequently limit starch biosynthesis, an ATP-consuming process, resulting in poor grain filling.

Gene and miRNA expression changes in squamous cell carcinoma of larynx and hypopharynx.

Laryngo-pharyngeal squamous cell carcinomas are one of the most common head and neck cancers. Despite the presence of a large body of information, molecular biomarkers are not currently used in the diagnosis, treatment and management of patients for this group of cancer. Here, we have profiled expression of genes and microRNAs of larynx and hypopharynx tumors using high-throughput sequencing experiments. We found that matrix metalloproteinases along with SCEL, CRNN, KRT4, SPINK5, and TGM3 among others have significantly altered expression in these tumors. Alongside gene expression, the microRNAs hsa-miR-139, hsa-miR-203 and the hsa-miR-424/503 cluster have aberrant expression in these cancers. Using target genes for these microRNAs, we found the involvement of pathways linked to cell cycle, p53 signaling, and viral carcinogenesis significant (P-values 10(-13), 10(-9) and 10(-7) respectively). Finally, using an ensemble machine-learning tool, we discovered a unique 8-gene signature for this group of cancers that differentiates the group from the other tumor subsites of head and neck region. We investigated the role of promoter methylation in one of these genes, WIF1, and found no correlation between DNA methylation and down-regulation of WIF1. We validated our findings of gene expression, 8-gene signature and promoter methylation using q-PCR, data from TCGA and q-MSP respectively. Data presented in this manuscript has been submitted to the NCBI Geo database with the accession number GSE67994.

Spikelet-specific variation in ethylene production and constitutive expression of ethylene receptors and signal transducers during grain filling of compact- and lax-panicle rice (Oryza sativa) cultivars.

Grain yields in modern super rice cultivars do not always meet the expectations because many spikelets are located on secondary branches in closely packed homogeneous distribution in these plants, and they do not fill properly. The factors limiting grain filling of such spikelets, especially in the lower panicle branches, are elusive. Two long-duration rice cultivars differing in panicle density, Mahalaxmi (compact) and Upahar (lax), were cultivated in an open field plot. Grain filling, ethylene production and constitutive expression of ethylene receptors and ethylene signal transducers in apical and basal spikelets of the panicle were compared during the early post-anthesis stage, which is the most critical period for grain development. In another experiment, a similar assessment was made for the medium-duration cultivars compact-panicle OR-1918 and lax-panicle Lalat. Grain weight of the apical spikelets was always higher than that of the basal spikelets. This gradient of grain weight was wide in the compact-panicle cultivars and narrow in the lax-panicle cultivars. Compared to apical spikelets, the basal spikelets produced more ethylene at anthesis and retained the capacity for post-anthesis expression of ethylene receptors and ethylene signal transducers longer. High ethylene production enhanced the expression of the RSR1 gene, but reduced expression of the GBSS1 gene. Ethylene inhibited the partitioning of assimilates of developing grains resulting in low starch biosynthesis and high accumulation of soluble carbohydrates. It is concluded that an increase in grain/spikelet density in rice panicles reduces apical dominance to the detriment of grain filling by production of ethylene and/or enhanced perception of the ethylene signal. Ethylene could be a second messenger for apical dominance in grain filling. The manipulation of the ethylene signal would possibly improve rice grain yield.

Integrated analysis of oral tongue squamous cell carcinoma identifies key variants and pathways linked to risk habits, HPV, clinical parameters and tumor recurrence.

Oral tongue squamous cell carcinomas (OTSCC) are a homogeneous group of tumors characterized by aggressive behavior, early spread to lymph nodes and a higher rate of regional failure. Additionally, the incidence of OTSCC among younger population (<50yrs) is on the rise; many of whom lack the typical associated risk factors of alcohol and/or tobacco exposure. We present data on single nucleotide variations (SNVs), indels, regions with loss of heterozygosity (LOH), and copy number variations (CNVs) from fifty-paired oral tongue primary tumors and link the significant somatic variants with clinical parameters, epidemiological factors including human papilloma virus (HPV) infection and tumor recurrence. Apart from the frequent somatic variants harbored in TP53, CASP8, RASA1, NOTCH and CDKN2A genes, significant amplifications and/or deletions were detected in chromosomes 6-9, and 11 in the tumors. Variants in CASP8 and CDKN2A were mutually exclusive. CDKN2A, PIK3CA, RASA1 and DMD variants were exclusively linked to smoking, chewing, HPV infection and tumor stage. We also performed a whole-genome gene expression study that identified matrix metalloproteases to be highly expressed in tumors and linked pathways involving arachidonic acid and NF-k-B to habits and distant metastasis, respectively. Functional knockdown studies in cell lines demonstrated the role of CASP8 in a HPV-negative OTSCC cell line. Finally, we identified a 38-gene minimal signature that predicts tumor recurrence using an ensemble machine-learning method. Taken together, this study links molecular signatures to various clinical and epidemiological factors in a homogeneous tumor population with a relatively high HPV prevalence.

Role of the N6-methyladenosine RNA mark in gene regulation and its implications on development and disease.

Epigenetics is a field that encompasses chemical modifications of DNA and histone proteins, both of which alter gene expression without changing the underlying nucleotide sequence. DNA methylation and modifications of histone tails have been studied in detail and are now known to be global gene regulatory mechanisms. An analogous post-transcriptional modification is chemical modification of specific nucleotides in RNA. Study of RNA modifications is a nascent field as yet, and the significance of these marks in controlling cell growth and differentiation is just beginning to be appreciated. The addition of a methyl group to adenosine (N-methyl-6-adenosine) or m6A is the most abundant modification in mammalian mRNAs. Though identified four decades ago, interest in this particular modification was set off by the discovery that the obesity gene FTO was an RNA demethylase. Since then, many studies have investigated m6A modification in different species. In this review, we summarize the current literature and hypotheses about the presence and function of this ubiquitous RNA modification in mammals, viruses, yeast and plants in terms of the consensus sequence and the methyltransferase/demethylation machinery identified thus far. We discuss its potential role in regulating molecular and physiological processes in each of these organisms, especially its role in RNA splicing, RNA degradation and development. We also enlist the methodologies developed so far, both locus-specific and transcriptome-wide, to study this modification. Lastly, we discuss whether m6A alterations have consequences on modulating disease aetiology, and speculate about its potential role in cancer.

COPS: a sensitive and accurate tool for detecting somatic Copy Number Alterations using short-read sequence data from paired samples.

Copy Number Alterations (CNAs) such as deletions and duplications; compose a larger percentage of genetic variations than single nucleotide polymorphisms or other structural variations in cancer genomes that undergo major chromosomal re-arrangements. It is, therefore, imperative to identify cancer-specific somatic copy number alterations (SCNAs), with respect to matched normal tissue, in order to understand their association with the disease. We have devised an accurate, sensitive, and easy-to-use tool, COPS, COpy number using Paired Samples, for detecting SCNAs. We rigorously tested the performance of COPS using short sequence simulated reads at various sizes and coverage of SCNAs, read depths, read lengths and also with real tumor:normal paired samples. We found COPS to perform better in comparison to other known SCNA detection tools for all evaluated parameters, namely, sensitivity (detection of true positives), specificity (detection of false positives) and size accuracy. COPS performed well for sequencing reads of all lengths when used with most upstream read alignment tools. Additionally, by incorporating a downstream boundary segmentation detection tool, the accuracy of SCNA boundaries was further improved. Here, we report an accurate, sensitive and easy to use tool in detecting cancer-specific SCNAs using short-read sequence data. In addition to cancer, COPS can be used for any disease as long as sequence reads from both disease and normal samples from the same individual are available. An added boundary segmentation detection module makes COPS detected SCNA boundaries more specific for the samples studied. COPS is available at ftp://115.119.160.213 with username "cops" and password "cops".

Customisation of the exome data analysis pipeline using a combinatorial approach.

The advent of next generation sequencing (NGS) technologies have revolutionised the way biologists produce, analyse and interpret data. Although NGS platforms provide a cost-effective way to discover genome-wide variants from a single experiment, variants discovered by NGS need follow up validation due to the high error rates associated with various sequencing chemistries. Recently, whole exome sequencing has been proposed as an affordable option compared to whole genome runs but it still requires follow up validation of all the novel exomic variants. Customarily, a consensus approach is used to overcome the systematic errors inherent to the sequencing technology, alignment and post alignment variant detection algorithms. However, the aforementioned approach warrants the use of multiple sequencing chemistry, multiple alignment tools, multiple variant callers which may not be viable in terms of time and money for individual investigators with limited informatics know-how. Biologists often lack the requisite training to deal with the huge amount of data produced by NGS runs and face difficulty in choosing from the list of freely available analytical tools for NGS data analysis. Hence, there is a need to customise the NGS data analysis pipeline to preferentially retain true variants by minimising the incidence of false positives and make the choice of right analytical tools easier. To this end, we have sampled different freely available tools used at the alignment and post alignment stage suggesting the use of the most suitable combination determined by a simple framework of pre-existing metrics to create significant datasets.