Pan-cancer analysis of pre-diagnostic blood metabolite concentrations in the European Prospective Investigation into Cancer and Nutrition.

BACKGROUND: Epidemiological studies of associations between metabolites and cancer risk have typically focused on specific cancer types separately. Here, we designed a multivariate pan-cancer analysis to identify metabolites potentially associated with multiple cancer types, while also allowing the investigation of cancer type-specific associations. METHODS: We analysed targeted metabolomics data available for 5828 matched case-control pairs from cancer-specific case-control studies on breast, colorectal, endometrial, gallbladder, kidney, localized and advanced prostate cancer, and hepatocellular carcinoma nested within the European Prospective Investigation into Cancer and Nutrition (EPIC) cohort. From pre-diagnostic blood levels of an initial set of 117 metabolites, 33 cluster representatives of strongly correlated metabolites and 17 single metabolites were derived by hierarchical clustering. The mutually adjusted associations of the resulting 50 metabolites with cancer risk were examined in penalized conditional logistic regression models adjusted for body mass index, using the data-shared lasso penalty. RESULTS: Out of the 50 studied metabolites, (i) six were inversely associated with the risk of most cancer types: glutamine, butyrylcarnitine, lysophosphatidylcholine a C18:2, and three clusters of phosphatidylcholines (PCs); (ii) three were positively associated with most cancer types: proline, decanoylcarnitine, and one cluster of PCs; and (iii) 10 were specifically associated with particular cancer types, including histidine that was inversely associated with colorectal cancer risk and one cluster of sphingomyelins that was inversely associated with risk of hepatocellular carcinoma and positively with endometrial cancer risk. CONCLUSIONS: These results could provide novel insights for the identification of pathways for cancer development, in particular those shared across different cancer types.

A New Pipeline for the Normalization and Pooling of Metabolomics Data.

Pooling metabolomics data across studies is often desirable to increase the statistical power of the analysis. However, this can raise methodological challenges as several preanalytical and analytical factors could introduce differences in measured concentrations and variability between datasets. Specifically, different studies may use variable sample types (e.g., serum versus plasma) collected, treated, and stored according to different protocols, and assayed in different laboratories using different instruments. To address these issues, a new pipeline was developed to normalize and pool metabolomics data through a set of sequential steps: (i) exclusions of the least informative observations and metabolites and removal of outliers; imputation of missing data; (ii) identification of the main sources of variability through principal component partial R-square (PC-PR2) analysis; (iii) application of linear mixed models to remove unwanted variability, including samples' originating study and batch, and preserve biological variations while accounting for potential differences in the residual variances across studies. This pipeline was applied to targeted metabolomics data acquired using Biocrates AbsoluteIDQ kits in eight case-control studies nested within the European Prospective Investigation into Cancer and Nutrition (EPIC) cohort. Comprehensive examination of metabolomics measurements indicated that the pipeline improved the comparability of data across the studies. Our pipeline can be adapted to normalize other molecular data, including biomarkers as well as proteomics data, and could be used for pooling molecular datasets, for example in international consortia, to limit biases introduced by inter-study variability. This versatility of the pipeline makes our work of potential interest to molecular epidemiologists.

Stabilization of G-quadruplex DNA structures in Schizosaccharomyces pombe causes single-strand DNA lesions and impedes DNA replication.

G-quadruplex (G4) structures are stable non-canonical DNA structures that are implicated in the regulation of many cellular pathways. We show here that the G4-stabilizing compound PhenDC3 causes growth defects in Schizosaccharomyces pombe cells, especially during S-phase in synchronized cultures. By visualizing individual DNA molecules, we observed shorter DNA fragments of newly replicated DNA in the PhenDC3-treated cells, suggesting that PhenDC3 impedes replication fork progression. Furthermore, a novel single DNA molecule damage assay revealed increased single-strand DNA lesions in the PhenDC3-treated cells. Moreover, chromatin immunoprecipitation showed enrichment of the leading-strand DNA polymerase at sites of predicted G4 structures, suggesting that these structures impede DNA replication. We tested a subset of these sites and showed that they form G4 structures, that they stall DNA synthesis in vitro and that they can be resolved by the breast cancer-associated Pif1 family helicases. Our results thus suggest that G4 structures occur in S. pombe and that stabilized/unresolved G4 structures are obstacles for the replication machinery. The increased levels of DNA damage might further highlight the association of the human Pif1 helicase with familial breast cancer and the onset of other human diseases connected to unresolved G4 structures.

A whole-genome sequenced control population in northern Sweden reveals subregional genetic differences.

The number of national reference populations that are whole-genome sequenced are rapidly increasing. Partly driving this development is the fact that genetic disease studies benefit from knowing the genetic variation typical for the geographical area of interest. A whole-genome sequenced Swedish national reference population (n = 1000) has been recently published but with few samples from northern Sweden. In the present study we have whole-genome sequenced a control population (n = 300) (ACpop) from Västerbotten County, a sparsely populated region in northern Sweden previously shown to be genetically different from southern Sweden. The aggregated variant frequencies within ACpop are publicly available (DOI 10.17044/NBIS/G000005) to function as a basic resource in clinical genetics and for genetic studies. Our analysis of ACpop, representing approximately 0.11% of the population in Västerbotten, indicates the presence of a genetic substructure within the county. Furthermore, a demographic analysis showed that the population from which samples were drawn was to a large extent geographically stationary, a finding that was corroborated in the genetic analysis down to the level of municipalities. Including ACpop in the reference population when imputing unknown variants in a Västerbotten cohort resulted in a strong increase in the number of high-confidence imputed variants (up to 81% for variants with minor allele frequency < 5%). ACpop was initially designed for cancer disease studies, but the genetic structure within the cohort will be of general interest for all genetic disease studies in northern Sweden.

Specific functions for Mediator complex subunits from different modules in the transcriptional response of Arabidopsis thaliana to abiotic stress.

Adverse environmental conditions are detrimental to plant growth and development. Acclimation to abiotic stress conditions involves activation of signaling pathways which often results in changes in gene expression via networks of transcription factors (TFs). Mediator is a highly conserved co-regulator complex and an essential component of the transcriptional machinery in eukaryotes. Some Mediator subunits have been implicated in stress-responsive signaling pathways; however, much remains unknown regarding the role of plant Mediator in abiotic stress responses. Here, we use RNA-seq to analyze the transcriptional response of Arabidopsis thaliana to heat, cold and salt stress conditions. We identify a set of common abiotic stress regulons and describe the sequential and combinatorial nature of TFs involved in their transcriptional regulation. Furthermore, we identify stress-specific roles for the Mediator subunits MED9, MED16, MED18 and CDK8, and putative TFs connecting them to different stress signaling pathways. Our data also indicate different modes of action for subunits or modules of Mediator at the same gene loci, including a co-repressor function for MED16 prior to stress. These results illuminate a poorly understood but important player in the transcriptional response of plants to abiotic stress and identify target genes and mechanisms as a prelude to further biochemical characterization.

A geographically matched control population efficiently limits the number of candidate disease-causing variants in an unbiased whole-genome analysis.

Whole-genome sequencing is a promising approach for human autosomal dominant disease studies. However, the vast number of genetic variants observed by this method constitutes a challenge when trying to identify the causal variants. This is often handled by restricting disease studies to the most damaging variants, e.g. those found in coding regions, and overlooking the remaining genetic variation. Such a biased approach explains in part why the genetic causes of many families with dominantly inherited diseases, in spite of being included in whole-genome sequencing studies, are left unsolved today. Here we explore the use of a geographically matched control population to minimize the number of candidate disease-causing variants without excluding variants based on assumptions on genomic position or functional predictions. To exemplify the benefit of the geographically matched control population we apply a typical disease variant filtering strategy in a family with an autosomal dominant form of colorectal cancer. With the use of the geographically matched control population we end up with 26 candidate variants genome wide. This is in contrast to the tens of thousands of candidates left when only making use of available public variant datasets. The effect of the local control population is dual, it (1) reduces the total number of candidate variants shared between affected individuals, and more importantly (2) increases the rate by which the number of candidate variants are reduced as additional affected family members are included in the filtering strategy. We demonstrate that the application of a geographically matched control population effectively limits the number of candidate disease-causing variants and may provide the means by which variants suitable for functional studies are identified genome wide.

Heterozygous colon cancer-associated mutations of SAMHD1 have functional significance.

Even small variations in dNTP concentrations decrease DNA replication fidelity, and this observation prompted us to analyze genomic cancer data for mutations in enzymes involved in dNTP metabolism. We found that sterile alpha motif and histidine-aspartate domain-containing protein 1 (SAMHD1), a deoxyribonucleoside triphosphate triphosphohydrolase that decreases dNTP pools, is frequently mutated in colon cancers, that these mutations negatively affect SAMHD1 activity, and that several SAMHD1 mutations are found in tumors with defective mismatch repair. We show that minor changes in dNTP pools in combination with inactivated mismatch repair dramatically increase mutation rates. Determination of dNTP pools in mouse embryos revealed that inactivation of one SAMHD1 allele is sufficient to elevate dNTP pools. These observations suggest that heterozygous cancer-associated SAMHD1 mutations increase mutation rates in cancer cells.

Expression of CXCL10 is associated with response to radiotherapy and overall survival in squamous cell carcinoma of the tongue.

Five-year survival for patients with oral cancer has been disappointingly stable during the last decades, creating a demand for new biomarkers and treatment targets. Lately, much focus has been set on immunomodulation as a possible treatment or an adjuvant increasing sensitivity to conventional treatments. The objective of this study was to evaluate the prognostic importance of response to radiotherapy in tongue carcinoma patients as well as the expression of the CXC-chemokines in correlation to radiation response in the same group of tumours. Thirty-eight patients with tongue carcinoma that had received radiotherapy followed by surgery were included. The prognostic impact of pathological response to radiotherapy, N-status, T-stage, age and gender was evaluated using Cox's regression models, Kaplan-Meier survival curves and chi-square test. The expression of 23 CXC-chemokine ligands and their receptors were evaluated in all patients using microarray and qPCR and correlated with response to treatment using logistic regression. Pathological response to radiotherapy was independently associated to overall survival with a 2-year survival probability of 81% for patients showing a complete pathological response, while patients with a non-complete response only had a probability of 42% to survive for 2 years (p = 0.016). The expression of one CXC-chemokine, CXCL10, was significantly associated with response to radiotherapy and the group of patients with the highest CXCL10 expression responded, especially poorly (p = 0.01). CXCL10 is a potential marker for response to radiotherapy and overall survival in patients with squamous cell carcinoma of the tongue.

Transcriptional profiling of formalin fixed paraffin embedded tissue: pitfalls and recommendations for identifying biologically relevant changes.

Expression profiling techniques have been used to study the biology of many types of cancer but have been limited to some extent by the requirement for collection of fresh tissue. In contrast, formalin fixed paraffin embedded (FFPE) samples are widely available and represent a vast resource of potential material. The techniques used to handle the degraded and modified RNA from these samples are relatively new and all the pitfalls and limitations of this material for whole genome expression profiling are not yet clarified. Here, we analyzed 70 FFPE tongue carcinoma samples and 17 controls using the whole genome DASL array covering nearly 21000 genes. We identified that sample age is related to quality of extracted RNA and that sample quality influences apparent expression levels in a non-random manner related to gene probe sequence, leading to spurious results. However, by removing sub-standard samples and analysing only those 28 cancers and 15 controls that had similar quality we were able to generate a list of 934 genes significantly altered in tongue cancer compared to control samples of tongue. This list contained previously identified changes and was enriched for genes involved in many cancer-related processes such as tissue remodelling, inflammation, differentiation and apoptosis. Four novel genes of potential importance in tongue cancer development and maintenance, SH3BGL2, SLC2A6, SLC16A3 and CXCL10, were independently confirmed, validating our data. Hence, gene expression profiling can be performed usefully on archival material if appropriate quality assurance steps are taken to ensure sample consistency and we present some recommendations for the use of FFPE material based on our findings.

Tubulin α-6 chain is a stably expressed reference gene in archival samples of normal oral tissue and oral squamous cell carcinoma.

One of the most critical factors in gene expression studies using quantitative real-time PCR is the choice of reference gene. Many of the commonly used reference genes have been shown to vary during a number of biological processes as well as between tissues. It is therefore important to always verify the stability of the gene of choice for all new tissues and experimental conditions. Here, we used two publicly available computer software packages (GeNorm and NormFinder) to investigate the stability of eight potential reference genes in formalin-fixed paraffin-embedded (FFPE) samples from normal oral tissue of different origin as well as from oral squamous cell carcinomas. Both programs found the tubulin α-6 chain (TUBA6) and ribosomal protein S13 (RPS13) to have the most stable expression between malignant and non-malignant tissue. NormFinder also found TUBA6 to be the most stable gene when samples were grouped according to tissue origin. FFPE samples constitute a large research resource, which considerably increases the number of samples available for analysis, leading to more reliable conclusions. Verification of a proper reference gene in oral FFPE tissue is therefore of great importance for future studies.

Gene expression profiling of archival tongue squamous cell carcinomas provides sub-classification based on DNA repair genes.

A subgroup of patients with squamous cell carcinoma of the head and neck (SCCHN) comprise young persons under the age of 40, who have not been heavily exposed to the classical risk factors, smoking and alcohol. The number of SCCHN in young adults, particularly tongue tumours, is increasing in several parts of the world. Here we employed a novel gene expression array methodology specifically developed for analysis of degraded RNA and investigated the expression of 502 cancer-related genes in archival paraffin-embedded SCCHN of the tongue from young (< or =40) and elderly patients (> or =50). Genes detected as de-regulated in tumours compared to non-malignant controls were in concordance with results from earlier studies of fresh frozen material. No genes were detected as significantly differentially expressed between young and old patients suggesting that the overall pathobiology of SCCHN is similar in young and old. Unsupervised clustering divided tumours into three groups, irrespective of age, where several differentially expressed DNA repair genes were a prominent separation factor. High levels of DNA repair genes associated with impaired therapeutic response to radiation, suggesting that DNA repair genes play a role in clinical outcome after radiotherapy.

Enhancer requirement for histone methylation linked with gene activation.

Enhancers cause a high level of transcription and activation of chromatin structure at target genes. Hyperacetylation of histones H3 and H4, a mark of active chromatin, is established broadly across target loci by enhancers that function over long distances. In the present study, we studied the role of an enhancer in methylation of various lysine residues on H3 by comparing a model gene locus having an active enhancer with one in which the enhancer has been inactivated within the context of minichromosomes. The intact enhancer affected histone methylation at K4, K9 and K36 in distinct ways depending on the methylation level and the location in the locus. All three lysine residues were highly tri-methylated in the coding region of the gene linked to the active enhancer but not the inactive enhancer. However di-methylation of K9 and K36 was not affected by the enhancer. The enhancer region itself was marked by mono-methylation at K4 and K9, distinguishing it from the methyl marks in the gene coding region. These results indicate that an enhancer has roles in establishing active histone methylation patterns linked with gene transcription rather than removing methylation linked with gene inactivation.