Analysis of cancer metabolism with high-throughput technologies

<p>Abstract</p> <p>Background</p> <p>Recent advances in genomics and proteomics have allowed us to study the nuances of the Warburg effect – a long-standing puzzle in cancer energy metabolism – at an unprecedented level of detail. While modern next-generation sequencing technologies are extremely powerful, the lack of appropriate data analysis tools makes this study difficult. To meet this challenge, we developed a novel application for comparative analysis of gene expression and visualization of RNA-Seq data.</p> <p>Results</p> <p>We analyzed two biological samples (normal human brain tissue and human cancer cell lines) with high-energy, metabolic requirements. We calculated digital topology and the copy number of every expressed transcript. We observed subtle but remarkable qualitative and quantitative differences between the citric acid (TCA) cycle and glycolysis pathways. We found that in the first three steps of the TCA cycle, digital expression of aconitase 2 (<it>ACO2</it>) in the brain exceeded both citrate synthase (<it>CS</it>) and isocitrate dehydrogenase 2 (<it>IDH2</it>), while in cancer cells this trend was quite the opposite. In the glycolysis pathway, all genes showed higher expression levels in cancer cell lines; and most notably, digital gene expression of glyceraldehyde-3-phosphate dehydrogenase (<it>GAPDH</it>) and enolase (<it>ENO</it>) were considerably increased when compared to the brain sample.</p> <p>Conclusions</p> <p>The variations we observed should affect the rates and quantities of ATP production. We expect that the developed tool will provide insights into the subtleties related to the causality between the Warburg effect and neoplastic transformation. Even though we focused on well-known and extensively studied metabolic pathways, the data analysis and visualization pipeline that we developed is particularly valuable as it is global and pathway-independent.</p

Gene expression analyses in breast cancer epidemiology: the Norwegian Women and Cancer postgenome cohort study

Introduction The introduction of high-throughput technologies, also called -omics technologies, into epidemiology has raised the need for high-quality observational studies to reduce several sources of error and bias. Methods The Norwegian Women and Cancer (NOWAC) postgenome cohort study consists of approximately 50,000 women born between 1943 and 1957 who gave blood samples between 2003 and 2006 and filled out a two-page questionnaire. Blood was collected in such a way that RNA is preserved and can be used for gene expression analyses. The women are part of the NOWAC study consisting of 172,471 women 30 to 70 years of age at recruitment from 1991 to 2006 who answered one to three questionnaires on diet, medication use, and lifestyle. In collaboration with the Norwegian Breast Cancer Group, every NOWAC participant born between 1943 and 1957 who is admitted to a collaborating hospital for a diagnostic biopsy or for surgery of breast cancer will be asked to donate a tumor biopsy and two blood samples. In parallel, at least three controls are approached for each breast cancer case in order to obtain blood samples from at least two controls per case. The controls are drawn at random from NOWAC matched by time of follow-up and age. In addition, 400 normal breast tissues as well as blood samples will be collected among healthy women participating at the Norwegian Mammography Screening program at the Breast Imaging Center at the University Hospital of North-Norway, Tromsø. Results The NOWAC postgenome cohort offers a unique opportunity (a) to study blood-derived gene expression profiles as a diagnostic test for breast cancer in a nested case-control design with adjustment for confounding factors related to different exposures, (b) to improve the reliability and accuracy of this approach by adjusting for an individual's genotype (for example, variants in genes coding for hormone and drug-metabolizing and detoxifying enzymes), (c) to study gene expression profiles from peripheral blood as surrogate tissue to biomonitor defined exposure (for example, hormone) and its association with disease risk (that is, breast cancer), and (d) to study gene variants (single nucleotide polymorphisms and copy number variations) and environmental exposure (endogenous and exogenous hormones) and their influence on the incidence of different molecular subtypes of breast cancer. Conclusion The NOWAC postgenome cohort combining a valid epidemiological approach with richness of biological samples should make an important contribution to the study of the etiology and system biology of breast cancer

Revealing the missing expressed genes beyond the human reference genome by RNA-Seq

<p>Abstract</p> <p>Background</p> <p>The complete and accurate human reference genome is important for functional genomics researches. Therefore, the incomplete reference genome and individual specific sequences have significant effects on various studies.</p> <p>Results</p> <p>we used two RNA-Seq datasets from human brain tissues and 10 mixed cell lines to investigate the completeness of human reference genome. First, we demonstrated that in previously identified ~5 Mb Asian and ~5 Mb African novel sequences that are absent from the human reference genome of NCBI build 36, ~211 kb and ~201 kb of them could be transcribed, respectively. Our results suggest that many of those transcribed regions are not specific to Asian and African, but also present in Caucasian. Then, we found that the expressions of 104 RefSeq genes that are unalignable to NCBI build 37 in brain and cell lines are higher than 0.1 RPKM. 55 of them are conserved across human, chimpanzee and macaque, suggesting that there are still a significant number of functional human genes absent from the human reference genome. Moreover, we identified hundreds of novel transcript contigs that cannot be aligned to NCBI build 37, RefSeq genes and EST sequences. Some of those novel transcript contigs are also conserved among human, chimpanzee and macaque. By positioning those contigs onto the human genome, we identified several large deletions in the reference genome. Several conserved novel transcript contigs were further validated by RT-PCR.</p> <p>Conclusion</p> <p>Our findings demonstrate that a significant number of genes are still absent from the incomplete human reference genome, highlighting the importance of further refining the human reference genome and curating those missing genes. Our study also shows the importance of <it>de novo </it>transcriptome assembly. The comparative approach between reference genome and other related human genomes based on the transcriptome provides an alternative way to refine the human reference genome.</p

Relative impact of key sources of systematic noise in Affymetrix and Illumina gene-expression microarray experiments

<p>Abstract</p> <p>Background</p> <p>Systematic processing noise, which includes batch effects, is very common in microarray experiments but is often ignored despite its potential to confound or compromise experimental results. Compromised results are most likely when re-analysing or integrating datasets from public repositories due to the different conditions under which each dataset is generated. To better understand the relative noise-contributions of various factors in experimental-design, we assessed several Illumina and Affymetrix datasets for technical variation between replicate hybridisations of Universal Human Reference (UHRR) and individual or pooled breast-tumour RNA.</p> <p>Results</p> <p>A varying degree of systematic noise was observed in each of the datasets, however in all cases the relative amount of variation between standard control RNA replicates was found to be greatest at earlier points in the sample-preparation workflow. For example, 40.6% of the total variation in reported expressions were attributed to replicate extractions, compared to 13.9% due to amplification/labelling and 10.8% between replicate hybridisations. Deliberate probe-wise batch-correction methods were effective in reducing the magnitude of this variation, although the level of improvement was dependent on the sources of noise included in the model. Systematic noise introduced at the chip, run, and experiment levels of a combined Illumina dataset were found to be highly dependant upon the experimental design. Both UHRR and pools of RNA, which were derived from the samples of interest, modelled technical variation well although the pools were significantly better correlated (4% average improvement) and better emulated the effects of systematic noise, over all probes, than the UHRRs. The effect of this noise was not uniform over all probes, with low GC-content probes found to be more vulnerable to batch variation than probes with a higher GC-content.</p> <p>Conclusions</p> <p>The magnitude of systematic processing noise in a microarray experiment is variable across probes and experiments, however it is generally the case that procedures earlier in the sample-preparation workflow are liable to introduce the most noise. Careful experimental design is important to protect against noise, detailed meta-data should always be provided, and diagnostic procedures should be routinely performed prior to downstream analyses for the detection of bias in microarray studies.</p

A phase II study of sequential neoadjuvant gemcitabine plus doxorubicin followed by gemcitabine plus cisplatin in patients with operable breast cancer: prediction of response using molecular profiling

This study examined the pathological complete response (pCR) rate and safety of sequential gemcitabine-based combinations in breast cancer. We also examined gene expression profiles from tumour biopsies to identify biomarkers predictive of response. Indian women with large or locally advanced breast cancer received 4 cycles of gemcitabine 1200 mg m−2 plus doxorubicin 60 mg m−2 (Gem+Dox), then 4 cycles of gemcitabine 1000 mg m−2 plus cisplatin 70 mg m−2 (Gem+Cis), and surgery. Three alternate dosing sequences were used during cycle 1 to examine dynamic changes in molecular profiles. Of 65 women treated, 13 (24.5% of 53 patients with surgery) had a pCR and 22 (33.8%) had a complete clinical response. Patients administered Gem d1, 8 and Dox d2 in cycle 1 (20 of 65) reported more toxicities, with G3/4 neutropenic infection/febrile neutropenia (7 of 20) as the most common cycle-1 event. Four drug-related deaths occurred. In 46 of 65 patients, 10-fold cross validated supervised analyses identified gene expression patterns that predicted with ⩾73% accuracy (1) clinical complete response after eight cycles, (2) overall clinical complete response, and (3) pCR. This regimen shows strong activity. Patients receiving Gem d1, 8 and Dox d2 experienced unacceptable toxicity, whereas patients on other sequences had manageable safety profiles. Gene expression patterns may predict benefit from gemcitabine-containing neoadjuvant therapy

A verified genomic reference sample for assessing performance of cancer panels detecting small variants of low allele frequency

BackgroundOncopanel genomic testing, which identifies important somatic variants, is increasingly common in medical practice and especially in clinical trials. Currently, there is a paucity of reliable genomic reference samples having a suitably large number of pre-identified variants for properly assessing oncopanel assay analytical quality and performance. The FDA-led Sequencing and Quality Control Phase 2 (SEQC2) consortium analyze ten diverse cancer cell lines individually and their pool, termed Sample A, to develop a reference sample with suitably large numbers of coding positions with known (variant) positives and negatives for properly evaluating oncopanel analytical performance.ResultsIn reference Sample A, we identify more than 40,000 variants down to 1% allele frequency with more than 25,000 variants having less than 20% allele frequency with 1653 variants in COSMIC-related genes. This is 5-100x more than existing commercially available samples. We also identify an unprecedented number of negative positions in coding regions, allowing statistical rigor in assessing limit-of-detection, sensitivity, and precision. Over 300 loci are randomly selected and independently verified via droplet digital PCR with 100% concordance. Agilent normal reference Sample B can be admixed with Sample A to create new samples with a similar number of known variants at much lower allele frequency than what exists in Sample A natively, including known variants having allele frequency of 0.02%, a range suitable for assessing liquid biopsy panels.ConclusionThese new reference samples and their admixtures provide superior capability for performing oncopanel quality control, analytical accuracy, and validation for small to large oncopanels and liquid biopsy assays.Peer reviewe

Comparative Analysis of Human Protein-Coding and Noncoding RNAs between Brain and 10 Mixed Cell Lines by RNA-Seq

In their expression process, different genes can generate diverse functional products, including various protein-coding or noncoding RNAs. Here, we investigated the protein-coding capacities and the expression levels of their isoforms for human known genes, the conservation and disease association of long noncoding RNAs (ncRNAs) with two transcriptome sequencing datasets from human brain tissues and 10 mixed cell lines. Comparative analysis revealed that about two-thirds of the genes expressed between brain and cell lines are the same, but less than one-third of their isoforms are identical. Besides those genes specially expressed in brain and cell lines, about 66% of genes expressed in common encoded different isoforms. Moreover, most genes dominantly expressed one isoform and some genes only generated protein-coding (or noncoding) RNAs in one sample but not in another. We found 282 human genes could encode both protein-coding and noncoding RNAs through alternative splicing in the two samples. We also identified more than 1,000 long ncRNAs, and most of those long ncRNAs contain conserved elements across either 46 vertebrates or 33 placental mammals or 10 primates. Further analysis showed that some long ncRNAs differentially expressed in human breast cancer or lung cancer, several of those differentially expressed long ncRNAs were validated by RT-PCR. In addition, those validated differentially expressed long ncRNAs were found significantly correlated with certain breast cancer or lung cancer related genes, indicating the important biological relevance between long ncRNAs and human cancers. Our findings reveal that the differences of gene expression profile between samples mainly result from the expressed gene isoforms, and highlight the importance of studying genes at the isoform level for completely illustrating the intricate transcriptome

Transcriptional activation of the Axl and PDGFR-α by c-Met through a ras- and Src-independent mechanism in human bladder cancer

<p>Abstract</p> <p>Background</p> <p>A cross-talk between different receptor tyrosine kinases (RTKs) plays an important role in the pathogenesis of human cancers.</p> <p>Methods</p> <p>Both NIH-Met5 and T24-Met3 cell lines harboring an inducible human c-Met gene were established. C-Met-related RTKs were screened by RTK microarray analysis. The cross-talk of RTKs was demonstrated by Western blotting and confirmed by small interfering RNA (siRNA) silencing, followed by elucidation of the underlying mechanism. The impact of this cross-talk on biological function was demonstrated by Trans-well migration assay. Finally, the potential clinical importance was examined in a cohort of 65 cases of locally advanced and metastatic bladder cancer patients.</p> <p>Results</p> <p>A positive association of Axl or platelet-derived growth factor receptor-alpha (PDGFR-α) with c-Met expression was demonstrated at translational level, and confirmed by specific siRNA knock-down. The transactivation of c-Met on Axl or PDGFR-α <it>in vitro </it>was through a <it>ras</it>- and Src-independent activation of mitogen-activated protein kinase/extracellular signal-regulated kinase (MEK/ERK) pathway. In human bladder cancer, co-expression of these RTKs was associated with poor patient survival (<it>p </it>< 0.05), and overexpression of c-Met/Axl/PDGFR-α or c-Met alone showed the most significant correlation with poor survival (<it>p </it>< 0.01).</p> <p>Conclusions</p> <p>In addition to c-Met, the cross-talk with Axl and/or PDGFR-α also contributes to the progression of human bladder cancer. Evaluation of Axl and PDGFR-α expression status may identify a subset of c-Met-positive bladder cancer patients who may require co-targeting therapy.</p

Tobacco use induces anti-apoptotic, proliferative patterns of gene expression in circulating leukocytes of Caucasian males

Abstract Background Strong epidemiologic evidence correlates tobacco use with a variety of serious adverse health effects, but the biological mechanisms that produce these effects remain elusive. Results We analyzed gene transcription data to identify expression spectra related to tobacco use in circulating leukocytes of 67 Caucasian male subjects. Levels of cotinine, a nicotine metabolite, were used as a surrogate marker for tobacco exposure. Significance Analysis of Microarray and Gene Set Analysis identified 109 genes in 16 gene sets whose transcription levels were differentially regulated by nicotine exposure. We subsequently analyzed this gene set by hyperclustering, a technique that allows the data to be clustered by both expression ratio and gene annotation (e.g. Gene Ontologies). Conclusion Our results demonstrate that tobacco use affects transcription of groups of genes that are involved in proliferation and apoptosis in circulating leukocytes. These transcriptional effects include a repertoire of transcriptional changes likely to increase the incidence of neoplasia through an altered expression of genes associated with transcription and signaling, interferon responses and repression of apoptotic pathways

Cross-oncopanel study reveals high sensitivity and accuracy with overall analytical performance depending on genomic regions

BackgroundTargeted sequencing using oncopanels requires comprehensive assessments of accuracy and detection sensitivity to ensure analytical validity. By employing reference materials characterized by the U.S. Food and Drug Administration-led SEquence Quality Control project phase2 (SEQC2) effort, we perform a cross-platform multi-lab evaluation of eight Pan-Cancer panels to assess best practices for oncopanel sequencing.ResultsAll panels demonstrate high sensitivity across targeted high-confidence coding regions and variant types for the variants previously verified to have variant allele frequency (VAF) in the 5-20% range. Sensitivity is reduced by utilizing VAF thresholds due to inherent variability in VAF measurements. Enforcing a VAF threshold for reporting has a positive impact on reducing false positive calls. Importantly, the false positive rate is found to be significantly higher outside the high-confidence coding regions, resulting in lower reproducibility. Thus, region restriction and VAF thresholds lead to low relative technical variability in estimating promising biomarkers and tumor mutational burden.ConclusionThis comprehensive study provides actionable guidelines for oncopanel sequencing and clear evidence that supports a simplified approach to assess the analytical performance of oncopanels. It will facilitate the rapid implementation, validation, and quality control of oncopanels in clinical use.Peer reviewe