Integrative multi-omics analysis identifies a prognostic miRNA signature and a targetable miR-21-3p/TSC2/mTOR axis in metastatic pheochromocytoma/paraganglioma.

Rationale: Pheochromocytomas and paragangliomas (PPGLs) are rare neuroendocrine tumors that present variable outcomes. To date, no effective therapies or reliable prognostic markers are available for patients who develop metastatic PPGL (mPPGL). Our aim was to discover robust prognostic markers validated through in vitro models, and define specific therapeutic options according to tumor genomic features. Methods: We analyzed three PPGL miRNome datasets (n=443), validated candidate markers and assessed them in serum samples (n=36) to find a metastatic miRNA signature. An integrative study of miRNome, transcriptome and proteome was performed to find miRNA targets, which were further characterized in vitro. Results: A signature of six miRNAs (miR-21-3p, miR-183-5p, miR-182-5p, miR-96-5p, miR-551b-3p, and miR-202-5p) was associated with metastatic risk and time to progression. A higher expression of five of these miRNAs was also detected in PPGL patients' liquid biopsies compared with controls. The combined expression of miR-21-3p/miR-183-5p showed the best power to predict metastasis (AUC=0.804, P=4.67·10-18), and was found associated in vitro with pro-metastatic features, such as neuroendocrine-mesenchymal transition phenotype, and increased cell migration rate. A pan-cancer multi-omic integrative study correlated miR-21-3p levels with TSC2 expression, mTOR pathway activation, and a predictive signature for mTOR inhibitor-sensitivity in PPGLs and other cancers. Likewise, we demonstrated in vitro a TSC2 repression and an enhanced rapamycin sensitivity upon miR-21-3p expression. Conclusions: Our findings support the assessment of miR-21-3p/miR-183-5p, in tumors and liquid biopsies, as biomarkers for risk stratification to improve the PPGL patients' management. We propose miR-21-3p to select mPPGL patients who may benefit from mTOR inhibitors.

TWEAK increases CD74 expression and sensitizes to DDT proinflammatory actions in tubular cells.

CD74 is a multifunctional protein and a receptor for Macrophage Migration Inhibitory Factor (MIF) and MIF-2 / D-dopachrome tautomerase (DDT) cytokines, upregulated in diabetic kidney disease. However, the drivers of CD74 expression and DDT function in kidney cells are poorly characterized. TWEAK is a proinflammatory cytokine that promotes kidney injury. We have now identified CD74 gene expression as upregulated in the kidneys in response to systemic TWEAK administration in mice, and have characterized the in vivo CD74 expression and the functional consequences in cultured cells. TWEAK administration to mice resulted in a progressive time-dependent (up to 24h) upregulation of kidney CD74 mRNA (RT-PCR) and protein (Western blot). Furthermore, the CD74 ligands MIF and DDT were also upregulated at the protein level 24h after TWEAK administration. Immunohistochemistry localized the increased CD74, MIF and DDT expression to tubular cells. In cultured tubular cells, TWEAK increased CD74 mRNA and protein expression dose-dependently, with a temporal pattern similar to in vivo. TWEAK-induced CD74 localized to the cell membrane, where it can function as a cytokine receptor. For the first time, we explored the actions of DDT in tubular cells and found that DDT amplified the increase in MCP-1 and RANTES expression in response to TWEAK. By contrast, DDT did not significantly modify TWEAK-induced Klotho downregulation. In conclusion, TWEAK upregulates CD74 and its ligands MIF and DDT in renal tubular cells. This may have functional consequences for kidney injury since DDT amplified the inflammatory response to TWEAK.

Saa3 is a key mediator of the protumorigenic properties of cancer-associated fibroblasts in pancreatic tumors.

Pancreatic ductal adenocarcinoma (PDAC) is characterized by the presence of abundant desmoplastic stroma primarily composed of cancer-associated fibroblasts (CAFs). It is generally accepted that CAFs stimulate tumor progression and might be implicated in drug resistance and immunosuppression. Here, we have compared the transcriptional profile of PDGFRα+ CAFs isolated from genetically engineered mouse PDAC tumors with that of normal pancreatic fibroblasts to identify genes potentially implicated in their protumorigenic properties. We report that the most differentially expressed gene, Saa3, a member of the serum amyloid A (SAA) apolipoprotein family, is a key mediator of the protumorigenic activity of PDGFRα+ CAFs. Whereas Saa3-competent CAFs stimulate the growth of tumor cells in an orthotopic model, Saa3-null CAFs inhibit tumor growth. Saa3 also plays a role in the cross talk between CAFs and tumor cells. Ablation of Saa3 in pancreatic tumor cells makes them insensitive to the inhibitory effect of Saa3-null CAFs. As a consequence, germline ablation of Saa3 does not prevent PDAC development in mice. The protumorigenic activity of Saa3 in CAFs is mediated by Mpp6, a member of the palmitoylated membrane protein subfamily of the peripheral membrane-associated guanylate kinases (MAGUK). Finally, we interrogated whether these observations could be translated to a human scenario. Indeed, SAA1, the ortholog of murine Saa3, is overexpressed in human CAFs. Moreover, high levels of SAA1 in the stromal component correlate with worse survival. These findings support the concept that selective inhibition of SAA1 in CAFs may provide potential therapeutic benefit to PDAC patients.

Bicycle: a bioinformatics pipeline to analyze bisulfite sequencing data.

Summary: High-throughput sequencing of bisulfite-converted DNA is a technique used to measure DNA methylation levels. Although a considerable number of computational pipelines have been developed to analyze such data, none of them tackles all the peculiarities of the analysis together, revealing limitations that can force the user to manually perform additional steps needed for a complete processing of the data. This article presents bicycle, an integrated, flexible analysis pipeline for bisulfite sequencing data. Bicycle analyzes whole genome bisulfite sequencing data, targeted bisulfite sequencing data and hydroxymethylation data. To show how bicycle overtakes other available pipelines, we compared them on a defined number of features that are summarized in a table. We also tested bicycle with both simulated and real datasets, to show its level of performance, and compared it to different state-of-the-art methylation analysis pipelines. Availability and implementation: Bicycle is publicly available under GNU LGPL v3.0 license at http://www.sing-group.org/bicycle. Users can also download a customized Ubuntu LiveCD including bicycle and other bisulfite sequencing data pipelines compared here. In addition, a docker image with bicycle and its dependencies, which allows a straightforward use of bicycle in any platform (e.g. Linux, OS X or Windows), is also available. Contact: ograna@cnio.es or dgpena@uvigo.es. Supplementary information: Supplementary data are available at Bioinformatics online.

Systems analysis identifies melanoma-enriched pro-oncogenic networks controlled by the RNA binding protein CELF1.

Melanomas are well-known for their altered mRNA expression profiles. Yet, the specific contribution of mRNA binding proteins (mRBPs) to melanoma development remains unclear. Here we identify a cluster of melanoma-enriched genes under the control of CUGBP Elav-like family member 1 (CELF1). CELF1 was discovered with a distinct prognostic value in melanoma after mining the genomic landscape of the 692 known mRBPs across different cancer types. Genome-wide transcriptomic, proteomic, and RNA-immunoprecipitation studies, together with loss-of-function analyses in cell lines, and histopathological evaluation in clinical biopsies, revealed an intricate repertoire of CELF1-RNA interactors with minimal overlap with other malignancies. This systems approach uncovered the oncogene DEK as an unexpected target and downstream effector of CELF1. Importantly, CELF1 and DEK were found to represent early-induced melanoma genes and adverse indicators of overall patient survival. These results underscore novel roles of CELF1 in melanoma, illustrating tumor type-restricted functions of RBPs in cancer.

Inactivation of Capicua in adult mice causes T-cell lymphoblastic lymphoma.

CIC (also known as Capicua) is a transcriptional repressor negatively regulated by RAS/MAPK signaling. Whereas the functions of Cic have been well characterized in Drosophila, little is known about its role in mammals. CIC is inactivated in a variety of human tumors and has been implicated recently in the promotion of lung metastases. Here, we describe a mouse model in which we inactivated Cic by selectively disabling its DNA-binding activity, a mutation that causes derepression of its target genes. Germline Cic inactivation causes perinatal lethality due to lung differentiation defects. However, its systemic inactivation in adult mice induces T-cell acute lymphoblastic lymphoma (T-ALL), a tumor type known to carry CIC mutations, albeit with low incidence. Cic inactivation in mice induces T-ALL by a mechanism involving derepression of its well-known target, Etv4 Importantly, human T-ALL also relies on ETV4 expression for maintaining its oncogenic phenotype. Moreover, Cic inactivation renders T-ALL insensitive to MEK inhibitors in both mouse and human cell lines. Finally, we show that Ras-induced mouse T-ALL as well as human T-ALL carrying mutations in the RAS/MAPK pathway display a genetic signature indicative of Cic inactivation. These observations illustrate that CIC inactivation plays a key role in this human malignancy.

Whole-body imaging of lymphovascular niches identifies pre-metastatic roles of midkine.

Cutaneous melanoma is a type of cancer with an inherent potential for lymph node colonization, which is generally preceded by neolymphangiogenesis. However, sentinel lymph node removal does not necessarily extend the overall survival of patients with melanoma. Moreover, lymphatic vessels collapse and become dysfunctional as melanomas progress. Therefore, it is unclear whether (and how) lymphangiogenesis contributes to visceral metastasis. Soluble and vesicle-associated proteins secreted by tumours and/or their stroma have been proposed to condition pre-metastatic sites in patients with melanoma. Still, the identities and prognostic value of lymphangiogenic mediators remain unclear. Moreover, our understanding of lymphangiogenesis (in melanomas and other tumour types) is limited by the paucity of mouse models for live imaging of distal pre-metastatic niches. Injectable lymphatic tracers have been developed, but their limited diffusion precludes whole-body imaging at visceral sites. Vascular endothelial growth factor receptor 3 (VEGFR3) is an attractive 'lymphoreporter' because its expression is strongly downregulated in normal adult lymphatic endothelial cells, but is activated in pathological situations such as inflammation and cancer. Here, we exploit this inducibility of VEGFR3 to engineer mouse melanoma models for whole-body imaging of metastasis generated by human cells, clinical biopsies or endogenously deregulated oncogenic pathways. This strategy revealed early induction of distal pre-metastatic niches uncoupled from lymphangiogenesis at primary lesions. Analyses of the melanoma secretome and validation in clinical specimens showed that the heparin-binding factor midkine is a systemic inducer of neo-lymphangiogenesis that defines patient prognosis. This role of midkine was linked to a paracrine activation of the mTOR pathway in lymphatic endothelial cells. These data support the use of VEGFR3 reporter mice as a 'MetAlert' discovery platform for drivers and inhibitors of metastasis.

Programmed mitophagy is essential for the glycolytic switch during cell differentiation.

Retinal ganglion cells (RGCs) are the sole projecting neurons of the retina and their axons form the optic nerve. Here, we show that embryogenesis-associated mouse RGC differentiation depends on mitophagy, the programmed autophagic clearance of mitochondria. The elimination of mitochondria during RGC differentiation was coupled to a metabolic shift with increased lactate production and elevated expression of glycolytic enzymes at the mRNA level. Pharmacological and genetic inhibition of either mitophagy or glycolysis consistently inhibited RGC differentiation. Local hypoxia triggered expression of the mitophagy regulator BCL2/adenovirus E1B 19-kDa-interacting protein 3-like (BNIP3L, best known as NIX) at peak RGC differentiation. Retinas from NIX-deficient mice displayed increased mitochondrial mass, reduced expression of glycolytic enzymes and decreased neuronal differentiation. Similarly, we provide evidence that NIX-dependent mitophagy contributes to mitochondrial elimination during macrophage polarization towards the proinflammatory and more glycolytic M1 phenotype, but not to M2 macrophage differentiation, which primarily relies on oxidative phosphorylation. In summary, developmentally controlled mitophagy promotes a metabolic switch towards glycolysis, which in turn contributes to cellular differentiation in several distinct developmental contexts.

Genomic and Molecular Screenings Identify Different Mechanisms for Acquired Resistance to MET Inhibitors in Lung Cancer Cells.

The development of resistance to tyrosine kinase inhibitors (TKI) limits the long-term efficacy of cancer treatments involving them. We aimed to understand the mechanisms that underlie acquired resistance (AR) to MET inhibitors in lung cancer. EBC1 cells, which have MET amplification and are sensitive to TKIs against MET, were used to generate multiple clones with AR to a MET-TKI. Whole-exome sequencing, RNA sequencing, and global DNA methylation analysis were used to scrutinize the genetic and molecular characteristics of the resistant cells. AR to the MET-TKI involved changes common to all resistant cells, that is, phenotypic modifications, specific changes in gene expression, and reactivation of AKT, ERK, and mTOR. The gene expression, global DNA methylation, and mutational profiles distinguished at least two groups of resistant cells. In one of these, the cells have acquired sensitivity to erlotinib, concomitantly with mutations of the KIRREL, HDAC11, HIATL1, and MAPK1IP1L genes, among others. In the other group, some cells have acquired inactivation of neurofibromatosis type 2 (NF2) concomitantly with strong overexpression of NRG1 and a mutational profile that includes changes in LMLN and TOMM34 Multiple independent and simultaneous strategies lead to AR to the MET-TKIs in lung cancer cells. The acquired sensitivity to erlotinib supports the known crosstalk between MET and the HER family of receptors. For the first time, we show inactivation of NF2 during acquisition of resistance to MET-TKI that may explain the refractoriness to erlotinib in these cells. Mol Cancer Ther; 16(7); 1366-76. ©2017 AACR.

Hepatocellular Carcinomas Originate Predominantly from Hepatocytes and Benign Lesions from Hepatic Progenitor Cells.

Hepatocellular carcinoma (HCC) is an aggressive primary liver cancer. However, its origin remains a debated question. Using human data and various hepatocarcinogenesis mouse models, we show that, in early stages, transformed hepatocytes, independent of their proliferation status, activate hepatic progenitor cell (HPC) expansion. Genetic lineage tracing of HPCs and hepatocytes reveals that, in all models, HCC originates from hepatocytes. However, whereas in various models tumors do not emanate from HPCs, tracking of progenitors in a model mimicking human hepatocarcinogenesis indicates that HPCs can generate benign lesions (regenerative nodules and adenomas) and aggressive HCCs. Mechanistically, galectin-3 and α-ketoglutarate paracrine signals emanating from oncogene-expressing hepatocytes instruct HPCs toward HCCs. α-Ketoglutarate preserves an HPC undifferentiated state, and galectin-3 maintains HPC stemness, expansion, and aggressiveness. Pharmacological or genetic blockage of galectin-3 reduces HCC, and its expression in human HCC correlates with poor survival. Our findings may have clinical implications for liver regeneration and HCC therapy.

Liver carcinogenesis by FOS-dependent inflammation and cholesterol dysregulation.

Human hepatocellular carcinomas (HCCs), which arise on a background of chronic liver damage and inflammation, express c-Fos, a component of the AP-1 transcription factor. Using mouse models, we show that hepatocyte-specific deletion of c-Fos protects against diethylnitrosamine (DEN)-induced HCCs, whereas liver-specific c-Fos expression leads to reversible premalignant hepatocyte transformation and enhanced DEN-carcinogenesis. c-Fos-expressing livers display necrotic foci, immune cell infiltration, and altered hepatocyte morphology. Furthermore, increased proliferation, dedifferentiation, activation of the DNA damage response, and gene signatures of aggressive HCCs are observed. Mechanistically, c-Fos decreases expression and activity of the nuclear receptor LXRα, leading to increased hepatic cholesterol and accumulation of toxic oxysterols and bile acids. The phenotypic consequences of c-Fos expression are partially ameliorated by the anti-inflammatory drug sulindac and largely prevented by statin treatment. An inverse correlation between c-FOS and the LXRα pathway was also observed in human HCC cell lines and datasets. These findings provide a novel link between chronic inflammation and metabolic pathways important in liver cancer.

Genomic Profiling of Patient-Derived Xenografts for Lung Cancer Identifies B2M Inactivation Impairing Immunorecognition.

Purpose: We aimed to maximize the performance of detecting genetic alterations in lung cancer using high-throughput sequencing for patient-derived xenografts (PDXs).Experimental Design: We undertook an integrated RNA and whole-exome sequencing of 14 PDXs. We focused on the genetic and functional analysis of ß2-microglobulin (B2M), a component of the HLA class-I complex.Results: We identified alterations in genes involved in various functions, such as B2M involved in immunosurveillance. We extended the mutational analysis of B2M to about 230 lung cancers. Five percent of the lung cancers carried somatic mutations, most of which impaired the correct formation of the HLA-I complex. We also report that genes such as CALR, PDIA3, and TAP1, which are involved in the maturation of the HLA-I complex, are altered in lung cancer. By gene expression microarrays, we observed that restitution of B2M in lung cancer cells upregulated targets of IFNα/IFNγ. Furthermore, one third of the lung cancers lacked the HLA-I complex, which was associated with lower cytotoxic CD8+ lymphocyte infiltration. The levels of B2M and HLA-I proteins correlated with those of PD-L1. Finally, a deficiency in HLA-I complex and CD8+ infiltration tended to correlate with reduced survival of patients with lung cancer treated with anti-PD-1/anti-PD-L1.Conclusions: Here, we report recurrent inactivation of B2M in lung cancer. These observations, coupled with the mutations found at CALR, PDIA3, and TAP1, and the downregulation of the HLA-I complex, indicate that an abnormal immunosurveillance axis contributes to lung cancer development. Finally, our observations suggest that an impaired HLA-I complex affects the response to anti-PD-1/anti-PD-L1 therapies. Clin Cancer Res; 23(12); 3203-13. ©2016 AACR.

Lineage-specific roles of the cytoplasmic polyadenylation factor CPEB4 in the regulation of melanoma drivers.

Nuclear 3'-end-polyadenylation is essential for the transport, stability and translation of virtually all eukaryotic mRNAs. Poly(A) tail extension can also occur in the cytoplasm, but the transcripts involved are incompletely understood, particularly in cancer. Here we identify a lineage-specific requirement of the cytoplasmic polyadenylation binding protein 4 (CPEB4) in malignant melanoma. CPEB4 is upregulated early in melanoma progression, as defined by computational and histological analyses. Melanoma cells are distinct from other tumour cell types in their dependency on CPEB4, not only to prevent mitotic aberrations, but to progress through G1/S cell cycle checkpoints. RNA immunoprecipitation, sequencing of bound transcripts and poly(A) length tests link the melanoma-specific functions of CPEB4 to signalling hubs specifically enriched in this disease. Essential in these CPEB4-controlled networks are the melanoma drivers MITF and RAB7A, a feature validated in clinical biopsies. These results provide new mechanistic links between cytoplasmic polyadenylation and lineage specification in melanoma.

Deep sequencing reveals microRNAs predictive of antiangiogenic drug response.

The majority of metastatic renal cell carcinoma (RCC) patients are treated with tyrosine kinase inhibitors (TKI) in first-line treatment; however, a fraction are refractory to these antiangiogenic drugs. MicroRNAs (miRNAs) are regulatory molecules proven to be accurate biomarkers in cancer. Here, we identified miRNAs predictive of progressive disease under TKI treatment through deep sequencing of 74 metastatic clear cell RCC cases uniformly treated with these drugs. Twenty-nine miRNAs were differentially expressed in the tumors of patients who progressed under TKI therapy (P values from 6 × 10-9 to 3 × 10-3). Among 6 miRNAs selected for validation in an independent series, the most relevant associations corresponded to miR-1307-3p, miR-155-5p, and miR-221-3p (P = 4.6 × 10-3, 6.5 × 10-3, and 3.4 × 10-2, respectively). Furthermore, a 2 miRNA-based classifier discriminated individuals with progressive disease upon TKI treatment (AUC = 0.75, 95% CI, 0.64-0.85; P = 1.3 × 10-4) with better predictive value than clinicopathological risk factors commonly used. We also identified miRNAs significantly associated with progression-free survival and overall survival (P = 6.8 × 10-8 and 7.8 × 10-7 for top hits, respectively), and 7 overlapped with early progressive disease. In conclusion, this is the first miRNome comprehensive study, to our knowledge, that demonstrates a predictive value of miRNAs for TKI response and provides a new set of relevant markers that can help rationalize metastatic RCC treatment.

NSD2 contributes to oncogenic RAS-driven transcription in lung cancer cells through long-range epigenetic activation.

The histone methyltransferase NSD2/WHSC1/MMSET is overexpressed in a number of solid tumors but its contribution to the biology of these tumors is not well understood. Here, we describe that NSD2 contributes to the proliferation of a subset of lung cancer cell lines by supporting oncogenic RAS transcriptional responses. NSD2 knock down combined with MEK or BRD4 inhibitors causes co-operative inhibitory responses on cell growth. However, while MEK and BRD4 inhibitors converge in the downregulation of genes associated with cancer-acquired super-enhancers, NSD2 inhibition affects the expression of clusters of genes embedded in megabase-scale regions marked with H3K36me2 and that contribute to the RAS transcription program. Thus, combinatorial therapies using MEK or BRD4 inhibitors together with NSD2 inhibition are likely to be needed to ensure a more comprehensive inhibition of oncogenic RAS-driven transcription programs in lung cancers with NSD2 overexpression.

Telomeric RNAs are essential to maintain telomeres.

Telomeres are transcribed generating long non-coding RNAs known as TERRA. Deciphering the role of TERRA has been one of the unsolved issues of telomere biology in the past decade. This has been, in part, due to lack of knowledge on the TERRA loci, thus preventing functional genetic studies. Here, we describe that long non-coding RNAs with TERRA features are transcribed from the human 20q and Xp subtelomeres. Deletion of the 20q locus by using the CRISPR-Cas9 technology causes a dramatic decrease in TERRA levels, while deletion of the Xp locus does not result in decreased TERRA levels. Strikingly, 20q-TERRA ablation leads to dramatic loss of telomere sequences and the induction of a massive DNA damage response. These findings identify chromosome 20q as a main TERRA locus in human cells and represent the first demonstration in any organism of the essential role of TERRA in the maintenance of telomeres.

Chronic skin inflammation leads to bone loss by IL-17-mediated inhibition of Wnt signaling in osteoblasts.

Inflammation has important roles in tissue regeneration, autoimmunity, and cancer. Different inflammatory stimuli can lead to bone loss by mechanisms that are not well understood. We show that skin inflammation induces bone loss in mice and humans. In psoriasis, one of the prototypic IL-17A-mediated inflammatory human skin diseases, low bone formation and bone loss correlated with increased serum IL-17A levels. Similarly, in two mouse models with chronic IL-17A-mediated skin inflammation,K14-IL17A(ind)andJunB(Δep), strong inhibition of bone formation was observed, different from classical inflammatory bone loss where osteoclast activation leads to bone degradation. We show that under inflammatory conditions, skin-resident cells such as keratinocytes, γδ T cells, and innate lymphoid cells were able to express IL-17A, which acted systemically to inhibit osteoblast and osteocyte function by a mechanism involving Wnt signaling. IL-17A led to decreased Wnt signaling in vitro, and importantly, pharmacological blockade of IL-17A rescued Wnt target gene expression and bone formation in vivo. These data provide a mechanism where IL-17A affects bone formation by regulating Wnt signaling in osteoblasts and osteocytes. This study suggests that using IL-17A blocking agents in psoriasis could be beneficial against bone loss in these patients.

MYC/PGC-1α Balance Determines the Metabolic Phenotype and Plasticity of Pancreatic Cancer Stem Cells.

The anti-diabetic drug metformin targets pancreatic cancer stem cells (CSCs), but not their differentiated progenies (non-CSCs), which may be related to distinct metabolic phenotypes. Here we conclusively demonstrate that while non-CSCs were highly glycolytic, CSCs were dependent on oxidative metabolism (OXPHOS) with very limited metabolic plasticity. Thus, mitochondrial inhibition, e.g., by metformin, translated into energy crisis and apoptosis. However, resistant CSC clones eventually emerged during treatment with metformin due to their intermediate glycolytic/respiratory phenotype. Mechanistically, suppression of MYC and subsequent increase of PGC-1α were identified as key determinants for the OXPHOS dependency of CSCs, which was abolished in resistant CSC clones. Intriguingly, no resistance was observed for the mitochondrial ROS inducer menadione and resistance could also be prevented/reversed for metformin by genetic/pharmacological inhibition of MYC. Thus, the specific metabolic features of pancreatic CSCs are amendable to therapeutic intervention and could provide the basis for developing more effective therapies to combat this lethal cancer.

The pluripotency factor NANOG promotes the formation of squamous cell carcinomas.

NANOG is a key pluripotency factor in embryonic stem cells that is frequently expressed in squamous cell carcinomas (SCCs). However, a direct link between NANOG and SCCs remains to be established. Here, we show that inducible overexpression of NANOG in mouse skin epithelia favours the malignant conversion of skin papillomas induced by chemical carcinogenesis, leading to increased SCC formation. Gene expression analyses in pre-malignant skin indicate that NANOG induces genes associated to epithelial-mesenchymal transition (EMT). Some of these genes are directly activated by NANOG, including EMT-associated genes Zeb1, Zeb2, Twist1, Prrx1 and miR-21. Finally, endogenous NANOG binds to the promoters of theses genes in human SCC cells and, moreover, NANOG induces EMT features in primary keratinocytes. These results provide in vivo evidence for the oncogenic role of NANOG in squamous cell carcinomas.

The contribution of cohesin-SA1 to gene expression and chromatin architecture in two murine tissues.

Cohesin, which in somatic vertebrate cells consists of SMC1, SMC3, RAD21 and either SA1 or SA2, mediates higher-order chromatin organization. To determine how cohesin contributes to the establishment of tissue-specific transcriptional programs, we compared genome-wide cohesin distribution, gene expression and chromatin architecture in cerebral cortex and pancreas from adult mice. More than one third of cohesin binding sites differ between the two tissues and these show reduced overlap with CCCTC-binding factor (CTCF) and are enriched at the regulatory regions of tissue-specific genes. Cohesin/CTCF sites at active enhancers and promoters contain, at least, cohesin-SA1. Analyses of chromatin contacts at the Protocadherin (Pcdh) and Regenerating islet-derived (Reg) gene clusters, mostly expressed in brain and pancreas, respectively, revealed remarkable differences that correlate with the presence of cohesin. We could not detect significant changes in the chromatin contacts at the Pcdh locus when comparing brains from wild-type and SA1 null embryos. In contrast, reduced dosage of SA1 altered the architecture of the Reg locus and decreased the expression of Reg genes in the pancreas of SA1 heterozygous mice. Given the role of Reg proteins in inflammation, such reduction may contribute to the increased incidence of pancreatic cancer observed in these animals.