High-throughput sequencing of small RNA transcriptomes reveals critical biological features targeted by microRNAs in cell models used for squamous cell cancer research

Abstract Background The implication of post-transcriptional regulation by microRNAs in molecular mechanisms underlying cancer disease is well documented. However, their interference at the cellular level is not fully explored. Functional in vitro studies are fundamental for the comprehension of their role; nevertheless results are highly dependable on the adopted cellular model. Next generation small RNA transcriptomic sequencing data of a tumor cell line and keratinocytes derived from primary culture was generated in order to characterize the microRNA content of these systems, thus helping in their understanding. Both constitute cell models for functional studies of microRNAs in head and neck squamous cell carcinoma (HNSCC), a smoking-related cancer. Known microRNAs were quantified and analyzed in the context of gene regulation. New microRNAs were investigated using similarity and structural search, ab initio classification, and prediction of the location of mature microRNAs within would-be precursor sequences. Results were compared with small RNA transcriptomic sequences from HNSCC samples in order to access the applicability of these cell models for cancer phenotype comprehension and for novel molecule discovery. Results Ten miRNAs represented over 70% of the mature molecules present in each of the cell types. The most expressed molecules were miR-21, miR-24 and miR-205, Accordingly; miR-21 and miR-205 have been previously shown to play a role in epithelial cell biology. Although miR-21 has been implicated in cancer development, and evaluated as a biomarker in HNSCC progression, no significant expression differences were seen between cell types. We demonstrate that differentially expressed mature miRNAs target cell differentiation and apoptosis related biological processes, indicating that they might represent, with acceptable accuracy, the genetic context from which they derive. Most miRNAs identified in the cancer cell line and in keratinocytes were present in tumor samples and cancer-free samples, respectively, with miR-21, miR-24 and miR-205 still among the most prevalent molecules at all instances. Thirteen miRNA-like structures, containing reads identified by the deep sequencing, were predicted from putative miRNA precursor sequences. Strong evidences suggest that one of them could be a new miRNA. This molecule was mostly expressed in the tumor cell line and HNSCC samples indicating a possible biological function in cancer. Conclusions Critical biological features of cells must be fully understood before they can be chosen as models for functional studies. Expression levels of miRNAs relate to cell type and tissue context. This study provides insights on miRNA content of two cell models used for cancer research. Pathways commonly deregulated in HNSCC might be targeted by most expressed and also by differentially expressed miRNAs. Results indicate that the use of cell models for cancer research demands careful assessment of underlying molecular characteristics for proper data interpretation. Additionally, one new miRNA-like molecule with a potential role in cancer was identified in the cell lines and clinical samples.

A NEW TUMOR SUPPRESSOR GENE CANDIDATE REGULATED BY THE NON-CODING RNA PCA3 IN HUMAN PROSTATE CANCER

Prostate cancer is the second leading cause of cancer-related death and the most common non-skin cancer in men in the USA. Considerable advancements in the practice of medicine have allowed a significant improvement in the diagnosis and treatment of this disease and, in recent years, both incidence and mortality rates have been slightly declining. However, it is still estimated that 1 man in 6 will be diagnosed with prostate cancer during his lifetime, and 1 man in 35 will die of the disease. In order to identify novel strategies and effective therapeutic approaches in the fight against prostate cancer, it is imperative to improve our understanding of its complex biology since many aspects of prostate cancer initiation and progression still remain elusive. The study of tumor biomarkers, due to their specific altered expression in tumor versus normal tissue, is a valid tool for elucidating key aspects of cancer biology, and may provide important insights into the molecular mechanisms underlining the tumorigenesis process of prostate cancer. PCA3, is considered the most specific prostate cancer biomarker, however its biological role, until now, remained unknown. PCA3 is a long non-coding RNA (ncRNA) expressed from chromosome 9q21 and its study led us to the discovery of a novel human gene, PC-TSGC, transcribed from the opposite strand and in an antisense orientation to PCA3. With the work presented in this thesis, we demonstrate that PCA3 exerts a negative regulatory role over PC-TSGC, and we propose PC-TSGC to be a new tumor suppressor gene that contrasts the transformation of prostate cells by inhibiting Rho-GTPases signaling pathways. Our findings provide a biological role for PCA3 in prostate cancer and suggest a new mechanism of tumor suppressor gene inactivation mediated by non-coding RNA. Also, the characterization of PCA3 and PC-TSGC led us to propose a new molecular pathway involving both genes in the transformation process of the prostate, thus providing a new piece of the jigsaw puzzle representing the complex biology of prostate cancer.

Spatial genomic heterogeneity within localized, multifocal prostate cancer

Herein we provide a detailed molecular analysis of the spatial heterogeneity of clinically localized, multifocal prostate cancer to delineate new oncogenes or tumor suppressors. We initially determined the copy number aberration (CNA) profiles of 74 patients with index tumors of Gleason score 7. Of these, 5 patients were subjected to whole-genome sequencing using DNA quantities achievable in diagnostic biopsies, with detailed spatial sampling of 23 distinct tumor regions to assess intraprostatic heterogeneity in focal genomics. Multifocal tumors are highly heterogeneous for single-nucleotide variants (SNVs), CNAs and genomic rearrangements. We identified and validated a new recurrent amplification of MYCL, which is associated with TP53 deletion and unique profiles of DNA damage and transcriptional dysregulation. Moreover, we demonstrate divergent tumor evolution in multifocal cancer and, in some cases, tumors of independent clonal origin. These data represent the first systematic relation of intraprostatic genomic heterogeneity to predicted clinical outcome and inform the development of novel biomarkers that reflect individual prognosis.

Genetic Predictors of Clinical Outcomes in Non-Small Cell Lung Cancer Patients

Lung cancer is the leading cause of cancer-related mortality in the US. Emerging evidence has shown that host genetic factors can interact with environmental exposures to influence patient susceptibility to the diseases as well as clinical outcomes, such as survival and recurrence. We aimed to identify genetic prognostic markers for non-small cell lung cancer (NSCLC), a major (85%) subtype of lung cancer, and also in other subgroups. With the fast evolution of genotyping technology, genetic association studies have went through candidate gene approach, to pathway-based approach, to the genome wide association study (GWAS). Even in the era of GWAS, pathway-based approach has its own advantages on studying cancer clinical outcomes: it is cost-effective, requiring a smaller sample size than GWAS easier to identify a validation population and explore gene-gene interactions. In the current study, we adopted pathway-based approach focusing on two critical pathways - miRNA and inflammation pathways. MicroRNAs (miRNA) post-transcriptionally regulate around 30% of human genes. Polymorphisms within miRNA processing pathways and binding sites may influence patients’ prognosis through altered gene regulation. Inflammation plays an important role in cancer initiation and progression, and also has shown to impact patients’ clinical outcomes. We first evaluated 240 single nucleotide polymorphisms (SNPs) in miRNA biogenesis genes and predicted binding sites in NSCLC patients to determine associations with clinical outcomes in early-stage (stage I and II) and late-stage (stage III and IV) lung cancer patients, respectively. First, in 535 early-stage patients, after correcting multiple comparisons, FZD4:rs713065 (hazard ratio [HR]:0.46, 95% confidence interval [CI]:0.32-0.65) showed a significant inverse association with survival in early stage surgery-only patients. SP1:rs17695156 (HR:2.22, 95% CI:1.44-3.41) and DROSHA:rs6886834 (HR:6.38, 95% CI:2.49-16.31) conferred increased risk of progression in the all patients and surgery-only populations, respectively. FAS:rs2234978 was significantly associated with improved survival in all patients (HR:0.59, 95% CI:0.44-0.77) and in the surgery plus chemotherapy populations (HR:0.19, 95% CI:0.07-0.46).. Functional genomics analysis demonstrated that this variant creates a miR-651 binding site resulting in altered miRNA regulation of FAS, providing biological plausibility for the observed association. We then analyzed these associations in 598 late-stage patients. After multiple comparison corrections, no SNPs remained significant in the late stage group, while the top SNP NAT1:rs15561 (HR=1.98, 96%CI=1.32-2.94) conferred a significantly increased risk of death in the chemotherapy subgroup. To test the hypothesis that genetic variants in the inflammation-related pathways may be associated with survival in NSCLC patients, we first conducted a three-stage study. In the discovery phase, we investigated a comprehensive panel of 11,930 inflammation-related SNPs in three independent lung cancer populations. A missense SNP (rs2071554) in HLA-DOB was significantly associated with poor survival in the discovery population (HR: 1.46, 95% CI: 1.02-2.09), internal validation population (HR: 1.51, 95% CI: 1.02-2.25), and external validation (HR: 1.52, 95% CI: 1.01-2.29) population. Rs2900420 in KLRK1 was significantly associated with a reduced risk for death in the discovery (HR: 0.76, 95% CI: 0.60-0.96) and internal validation (HR: 0.77, 95% CI: 0.61-0.99) populations, and the association reached borderline significance in the external validation population (HR: 0.80, 95% CI: 0.63-1.02). We also evaluated these inflammation-related SNPs in NSCLC patients in never smokers. Lung cancer in never smokers has been increasingly recognized as distinct disease from that in ever-smokers. A two-stage study was performed using a discovery population from MD Anderson (411 patients) and a validation population from Mayo Clinic (311 patients). Three SNPs (IL17RA:rs879576, BMP8A:rs698141, and STK:rs290229) that were significantly associated with survival were validated (pCD74:rs1056400 and CD38:rs10805347) were borderline significant (p=0.08) in the Mayo Clinic population. In the combined analysis, IL17RA:rs879576 resulted in a 40% reduction in the risk for death (p=4.1 × 10-5 [p=0.61, heterogeneity test]). We also validated a survival tree created in MD Anderson population in the Mayo Clinic population. In conclusion, our results provided strong evidence that genetic variations in specific pathways that examined (miRNA and inflammation pathways) influenced clinical outcomes in NSCLC patients, and with further functional studies, the novel loci have potential to be translated into clinical use.

Predictive gene lists for breast cancer prognosis: a topographic visualisation study

Background: The controversy surrounding the non-uniqueness of predictive gene lists (PGL) of small selected subsets of genes from very large potential candidates as available in DNA microarray experiments is now widely acknowledged 1. Many of these studies have focused on constructing discriminative semi-parametric models and as such are also subject to the issue of random correlations of sparse model selection in high dimensional spaces. In this work we outline a different approach based around an unsupervised patient-specific nonlinear topographic projection in predictive gene lists. Methods: We construct nonlinear topographic projection maps based on inter-patient gene-list relative dissimilarities. The Neuroscale, the Stochastic Neighbor Embedding(SNE) and the Locally Linear Embedding(LLE) techniques have been used to construct two-dimensional projective visualisation plots of 70 dimensional PGLs per patient, classifiers are also constructed to identify the prognosis indicator of each patient using the resulting projections from those visualisation techniques and investigate whether a-posteriori two prognosis groups are separable on the evidence of the gene lists. A literature-proposed predictive gene list for breast cancer is benchmarked against a separate gene list using the above methods. Generalisation ability is investigated by using the mapping capability of Neuroscale to visualise the follow-up study, but based on the projections derived from the original dataset. Results: The results indicate that small subsets of patient-specific PGLs have insufficient prognostic dissimilarity to permit a distinction between two prognosis patients. Uncertainty and diversity across multiple gene expressions prevents unambiguous or even confident patient grouping. Comparative projections across different PGLs provide similar results. Conclusion: The random correlation effect to an arbitrary outcome induced by small subset selection from very high dimensional interrelated gene expression profiles leads to an outcome with associated uncertainty. This continuum and uncertainty precludes any attempts at constructing discriminative classifiers. However a patient's gene expression profile could possibly be used in treatment planning, based on knowledge of other patients' responses. We conclude that many of the patients involved in such medical studies are intrinsically unclassifiable on the basis of provided PGL evidence. This additional category of 'unclassifiable' should be accommodated within medical decision support systems if serious errors and unnecessary adjuvant therapy are to be avoided.

Low frequency of cigarette smoking and the risk of head and neck cancer in the INHANCE consortium pooled analysis

Funding • The pooled data coordination team (PBoffetta, MH, YCAL) were supported by National Cancer Institute grant R03CA113157 and by National Institute of Dental and Craniofacial Research grant R03DE016611 • The Milan study (CLV) was supported by the Italian Association for Research on Cancer (Grant no. 10068). • The Aviano study (LDM) was supported by a grant from the Italian Association for Research on Cancer (AIRC), Italian League Against Cancer and Italian Ministry of Research • The Italy Multicenter study (DS) was supported by the Italian Association for Research on Cancer (AIRC), Italian League Against Cancer and Italian Ministry of Research. • The Study from Switzerland (FL) was supported by the Swiss League against Cancer and the Swiss Research against Cancer/Oncosuisse [KFS-700, OCS-1633]. • The central Europe study (PBoffetta, PBrenan, EF, JL, DM, PR, OS, NS-D) was supported by the World Cancer Research Fund and the European Commission INCOCOPERNICUS Program [Contract No. IC15- CT98-0332] • The New York multicentre study (JM) was supported by a grant from National Institute of Health [P01CA068384 K07CA104231]. • The study from the Fred Hutchison Cancer Research Center from Seattle (CC, SMS) was supported by a National Institute of Health grant [R01CA048996, R01DE012609]. • The Iowa study (ES) was supported by National Instituteof Health [NIDCR R01DE011979, NIDCR R01DE013110, FIRCA TW001500] and Veterans Affairs Merit Review Funds. • The North Carolina studies (AFO) were supported by National Institute of Health [R01CA061188], and in part by a grant from the National Institute of Environmental Health Sciences [P30ES010126]. • The Tampa study (PLazarus, JM) was supported by National Institute of Health grants [P01CA068384, K07CA104231, R01DE013158] • The Los Angeles study (Z-F Z, HM) was supported by grants from National Institute of Health [P50CA090388, R01DA011386, R03CA077954, T32CA009142, U01CA096134, R21ES011667] and the Alper Research Program for Environmental Genomics of the UCLA Jonsson Comprehensive Cancer Center. • The Houston study (EMS, GL) was supported by a grant from National Institute of Health [R01ES011740, R01CA100264]. • The Puerto Rico study (RBH, MPP) was supported by a grant from National Institutes of Health (NCI) US and NIDCR intramural programs. • The Latin America study (PBoffetta, PBrenan, MV, LF, MPC, AM, AWD, SK, VW-F) was supported by Fondo para la Investigacion Cientifica y Tecnologica (FONCYT) Argentina, IMIM (Barcelona), Fundaco de Amparo a‘ Pesquisa no Estado de Sao Paulo (FAPESP) [No 01/01768-2], and European Commission [IC18-CT97-0222] • The IARC multicentre study (SF, RH, XC) was supported by Fondo de Investigaciones Sanitarias (FIS) of the Spanish Government [FIS 97/ 0024, FIS 97/0662, BAE 01/5013], International Union Against Cancer (UICC), and Yamagiwa-Yoshida Memorial International Cancer Study Grant. • The Boston study (KKelsey, MMcC) was supported by a grant from National Institute of Health [R01CA078609, R01CA100679]. • The Rome study (SB, GC) was supported by AIRC (Italian Agency for Research on Cancer). • The US multicentre study (BW) was supported by The Intramural Program of the National Cancer Institute, National Institute of Health, United States. • The Sao Paolo study (V W-F) was supported by Fundacao de Ampara a Pesquisa no Estado de Sao Paulo (FAPESP No 10/51168-0) • The MSKCC study (SS, G-P Y) was supported by a grant from National Institute of Health [R01CA051845]. • The Seattle-Leo stud (FV) was supported by a grant from National Institute of Health [R01CA030022] • The western Europe Study (PBoffetta, IH, WA, PLagiou, DS, LS, FM, CH, KKjaerheim, DC, TMc, PT, AA, AZ) was supported by European Community (5th Frame work Programme) grant no QLK1-CT-2001- 00182. • The Germany Heidelberg study (HR) was supported by the grant No. 01GB9702/3 from the German Ministry of Education and Research.

Shooting for the Moon or Flying Too Near the Sun? Crossing the Value Rubicon in Precision Cancer Care

In his last two State of the Union addresses, President Barack Obama has focused on the need to deliver innovative solutions to improve human health, through the Precision Medicine Initiative in 2015 and the recently announced Cancer Moonshot in 2016. Precision cancer care has delivered clear patient benefit, but even for high-impact medicines such as imatinib mesylate (Glivec) in chronic myeloid leukaemia, the excitement at the success of this practice-changing clinical intervention has been somewhat tempered by the escalating price of this 'poster child' for precision cancer medicine (PCM). Recent studies on the costs of cancer drugs have revealed significant price differentials, which are a major causative factor behind disparities in the access to new generations of immunological and molecularly targeted agents. In this perspective, we will discuss the benefits of PCM to modern cancer control, but also emphasise how increasing costs are rendering the current approaches to integrating the paradigm of PCM unsustainable. Despite the ever increasing pressure on cancer and health care budgets, innovation will and must continue. Value-based frameworks offer one of the most rational approaches for policymakers committed to improving cancer outcomes through a public health approach.

Molecular and biochemical mechanisms involved in the toxicity of a marine cyanobacteria compound on cancer cell lines

In recent years marine biotechnology has revealed a crucial role in the future of bioindustry. Among the many marine resources, cyanobacteria have shown great potential in the production of bioactive compounds with diverse applicability. The pharmacological potential of these organisms has been one of the most explored areas in particular its antibacterial, antifungal and anticancer potential. This work was based on the assessment of potential anticancer compound E13010 F 5.4 isolated from marine cyanobacteria strain Synechocystis salina LEGE 06099. Thus the aim of this work was to explore molecular and biochemical mechanisms underlying the bioactivity detected in human cancer cells, specifically in lines RKO colon carcinoma and HT-29. The isolation of the compound was performed from biomass obtained by large-scale culture. To obtain the compound fractionation was carried and confirmation and isolation performed by Nuclear Magnetic Resonance (NMR), Thin Layer Chromatography (TLC) and High-Performance Liquid Chromatography (HPLC). Cell viability assays were performed based on reduction of 3- (4,5-dimetiltiaziol-2-yl) -2,5-diphenyltetrazolium bromide (MTT) to assess the cytotoxic potential of the compound. From the battery of cell lines RKO (colon carcinoma), HT-29 (colorectal adenocarcinoma), MG-63 (osteosarcoma) and T47D (breast carcinoma) the cell lines RKO and HT-29 were selected for elucidation of mechanisms of cytotoxicity. For the elucidation of the mechanisms involved in cytotoxicity the cell lines RKO and HT29 were exposed to the compound. A genomic approach based in the mRNA expression of genes involved in apoptosis and cell cycle by Real-Time PCR and a proteomic approach based on the separation of proteins by two-dimensional electrophoresis (2DGE) was performed. For mRNA expression were selected the genes RPL8, HPRT1, VDAC, SHMT2, CCNE, CCNB1, P21CIP, BCL-2 and BAD and for proteomics isoelectric focussing between 3 – 10 and molecular weight of 19 – 117 kDa separated by polyacrylamide gels (2DGE). The MTT results confirmed the reduction of the cell viability. The RT-PCR results for the expression of genes studied were not yet fully elucidative. For the cell line RKO there was a significant reduction in the expression of the gene P21CIP, and a tendency for reduction in the BAD gene expression and for increased expression of gene CCNB1, pointing to an effort for cell proliferation. In HT-29 cell line, there was a tendency for increase in the expression of P21CIP and BAD, which may explain the reduction in cell viability. The 2DGE results indicate proteomic patterns with differentially altered spots in the treated and control cells with both qualitative and quantitative differences, and differences in response between the RKO and HT-29 cell lines.

No evidence for association of ataxia-telangiectasia mutated gene T2119C and C3161G amino acid substitution variants with risk of breast cancer

Background There is evidence that certain mutations in the double-strand break repair pathway ataxia-telangiectasia mutated gene act in a dominant-negative manner to increase the risk of breast cancer. There are also some reports to suggest that the amino acid substitution variants T2119C Ser707Pro and C3161G Pro1054Arg may be associated with breast cancer risk. We investigate the breast cancer risk associated with these two nonconservative amino acid substitution variants using a large Australian population-based case–control study. Methods The polymorphisms were genotyped in more than 1300 cases and 600 controls using 5' exonuclease assays. Case–control analyses and genotype distributions were compared by logistic regression. Results The 2119C variant was rare, occurring at frequencies of 1.4 and 1.3% in cases and controls, respectively (P = 0.8). There was no difference in genotype distribution between cases and controls (P = 0.8), and the TC genotype was not associated with increased risk of breast cancer (adjusted odds ratio = 1.08, 95% confidence interval = 0.59–1.97, P = 0.8). Similarly, the 3161G variant was no more common in cases than in controls (2.9% versus 2.2%, P = 0.2), there was no difference in genotype distribution between cases and controls (P = 0.1), and the CG genotype was not associated with an increased risk of breast cancer (adjusted odds ratio = 1.30, 95% confidence interval = 0.85–1.98, P = 0.2). This lack of evidence for an association persisted within groups defined by the family history of breast cancer or by age. Conclusion The 2119C and 3161G amino acid substitution variants are not associated with moderate or high risks of breast cancer in Australian women.