Divergence of mutational signatures in association with breast cancer subtype.

Abnormal molecular processes occurring throughout the genome leave distinct somatic mutational patterns termed mutational signatures. Exploring the associations between mutational signatures and clinicopathological features can unravel potential mechanisms driving tumorigenic processes. We analyzed whole genome sequencing (WGS) data of tumor and peripheral blood samples from 37 primary breast cancer (BC) patients receiving neoadjuvant chemotherapy. Comprehensive clinico-pathologic features were correlated with genomic profiles and mutational signatures. Somatic mutational landscapes were highly concordant with known BC data sets. Remarkably, we observed a divergence of dominant mutational signatures in association with BC subtype. Signature 5 was overrepresented in hormone receptor positive (HR+) patients, whereas triple-negative tumors mostly lacked Signature 5, but expectedly overrepresented Signature 3. We validated these findings in a large WGS data set of BC, demonstrating dominance of Signature 5 in HR+ patients, mostly in luminal A subtype. We further investigated the association between Signature 5 and gene expression signatures, and identified potential networks, likely related to estrogen regulation. Our results suggest that the yet elusive Signature 5 represents an alternative mechanism for mutation accumulation in HR+ BC, independent of the homologous recombination repair machinery related to Signature 3. This study provides theoretical basis for further elucidating the processes promoting hormonal breast carcinogenesis.

MiR-199a-3p Induces Mesenchymal to Epithelial Transition of Keratinocytes by Targeting RAP2B.

Cutaneous squamous cell carcinoma (CSCC) is an epidermal skin cancer that evolves from normal epidermis along several pre-malignant stages. Previously we found specific miRNAs alterations in each step along these stages. miR-199a-3p expression decreases at the transition to later stages. A crucial step for epithelial carcinoma cells to acquire invasive capacity is the disruption of cell-cell contacts and the gain of mesenchymal motile phenotype, a process known as epithelial-to-mesenchymal transition (EMT). This study aims to study the role of decreased expression of miR-199a-3p in keratinocytes' EMT towards carcinogenesis. First, we measured miR-199a-3p in different stages of epidermal carcinogenesis. Then, we applied Photoactivatable Ribonucleoside-Enhanced Crosslinking and Immunoprecipitation (PAR-CLIP) assay to search for possible biochemical targets of miR-199a-3p and verified that Ras-associated protein B2 (RAP2B) is a bona-fide target of miR-199a-3p. Next, we analyzed RAP2B expression, in CSCC biopsies. Last, we evaluated possible mechanisms leading to decreased miR-199a-3p expression. miR-199a-3p induces a mesenchymal to epithelial transition (MET) in CSSC cells. Many of the under-expressed genes in CSCC overexpressing miR-199a-3p, are possible targets of miR-199a-3p and play roles in EMT. RAP2B is a biochemical target of miR-199a-3p. Overexpression of miR-199a-3p in CSCC results in decreased phosphorylated focal adhesion kinase (FAK). In addition, inhibiting FAK phosphorylation inhibits EMT marker genes' expression. In addition, we proved that DNA methylation is part of the mechanism by which miR-199a-3p expression is inhibited. However, it is not by the methylation of miR-199a putative promoter. These findings suggest that miR-199a-3p inhibits the EMT process by targeting RAP2B. Inhibitors of RAP2B or FAK may be effective therapeutic agents for CSCC.

NF-κB-miR-155 axis activation mediates ovulation-induced oncogenic effects in fallopian tube epithelium.

The fallopian tube secretory epithelial cells (FTSECs) are the cell-of-origin of most high-grade serous ovarian carcinomas (HGSOC). FTSECs are repeatedly exposed to inflammation induced by follicular fluid (FF) that is released with every ovulation cycle throughout a woman's reproductive years. Uninterrupted ovulation cycles are an established risk factor for HGSOC. Stimuli present in the FF induce an inflammatory environment which may cause DNA damage eventually leading to serous tumorigenesis. With the aim of elucidating possible mechanistic pathways, we established an 'ex vivo persistent ovulation model' mimicking the repeated exposure of human benign fallopian tube epithelium (FTE) to FF. We performed mass spectrometry analysis of the secretome of the ex vivo cultures as well as confirmatory targeted expressional and functional analyses. We demonstrated activation of the NF-κB pathway and upregulation of miR-155 following short-term exposure of FTE to human FF. Increased expression of miR-155 was also detected in primary HGSOC tumors compared with benign primary human FTE and corresponded with changes in the expression of miR-155 target genes. The phenotype of miR-155 overexpression in FTSEC cell line is of increased migratory and altered adhesion capacities. Overall, activation of the NF-κB-miR-155 axis in FTE may represent a possible link between ovulation-induced inflammation, DNA damage, and transcriptional changes that may eventually lead to serious carcinogenesis.

Two phases of mitogenic signaling unveil roles for p53 and EGR1 in elimination of inconsistent growth signals.

Normal cells require continuous exposure to growth factors in order to cross a restriction point and commit to cell-cycle progression. This can be replaced by two short, appropriately spaced pulses of growth factors, where the first pulse primes a process, which is completed by the second pulse, and enables restriction point crossing. Through integration of comprehensive proteomic and transcriptomic analyses of each pulse, we identified three processes that regulate restriction point crossing: (1) The first pulse induces essential metabolic enzymes and activates p53-dependent restraining processes. (2) The second pulse eliminates, via the PI3K/AKT pathway, the suppressive action of p53, as well as (3) sets an ERK-EGR1 threshold mechanism, which digitizes graded external signals into an all-or-none decision obligatory for S phase entry. Together, our findings uncover two gating mechanisms, which ensure that cells ignore fortuitous growth factors and undergo proliferation only in response to consistent mitogenic signals.

Topology of the human and mouse m6A RNA methylomes revealed by m6A-seq.

An extensive repertoire of modifications is known to underlie the versatile coding, structural and catalytic functions of RNA, but it remains largely uncharted territory. Although biochemical studies indicate that N(6)-methyladenosine (m(6)A) is the most prevalent internal modification in messenger RNA, an in-depth study of its distribution and functions has been impeded by a lack of robust analytical methods. Here we present the human and mouse m(6)A modification landscape in a transcriptome-wide manner, using a novel approach, m(6)A-seq, based on antibody-mediated capture and massively parallel sequencing. We identify over 12,000 m(6)A sites characterized by a typical consensus in the transcripts of more than 7,000 human genes. Sites preferentially appear in two distinct landmarks--around stop codons and within long internal exons--and are highly conserved between human and mouse. Although most sites are well preserved across normal and cancerous tissues and in response to various stimuli, a subset of stimulus-dependent, dynamically modulated sites is identified. Silencing the m(6)A methyltransferase significantly affects gene expression and alternative splicing patterns, resulting in modulation of the p53 (also known as TP53) signalling pathway and apoptosis. Our findings therefore suggest that RNA decoration by m(6)A has a fundamental role in regulation of gene expression.

A-to-I RNA editing occurs at over a hundred million genomic sites, located in a majority of human genes.

RNA molecules transmit the information encoded in the genome and generally reflect its content. Adenosine-to-inosine (A-to-I) RNA editing by ADAR proteins converts a genomically encoded adenosine into inosine. It is known that most RNA editing in human takes place in the primate-specific Alu sequences, but the extent of this phenomenon and its effect on transcriptome diversity are not yet clear. Here, we analyzed large-scale RNA-seq data and detected ∼1.6 million editing sites. As detection sensitivity increases with sequencing coverage, we performed ultradeep sequencing of selected Alu sequences and showed that the scope of editing is much larger than anticipated. We found that virtually all adenosines within Alu repeats that form double-stranded RNA undergo A-to-I editing, although most sites exhibit editing at only low levels (<1%). Moreover, using high coverage sequencing, we observed editing of transcripts resulting from residual antisense expression, doubling the number of edited sites in the human genome. Based on bioinformatic analyses and deep targeted sequencing, we estimate that there are over 100 million human Alu RNA editing sites, located in the majority of human genes. These findings set the stage for exploring how this primate-specific massive diversification of the transcriptome is utilized.

ADAR1 deletion induces NFκB and interferon signaling dependent liver inflammation and fibrosis.

Adenosine deaminase acting on RNA (ADAR) 1 binds and edits double-stranded (ds) RNA secondary structures found mainly within untranslated regions of many transcripts. In the current research, our aim was to study the role of ADAR1 in liver homeostasis. As previous studies show a conserved immunoregulatory function for ADAR1 in mammalians, we focused on its role in preventing chronic hepatic inflammation and the associated activation of hepatic stellate cells to produce extracellular matrix and promote fibrosis. We show that hepatocytes specific ADAR1 knock out (KO) mice display massive liver damage with multifocal inflammation and fibrogenesis. The bioinformatics analysis of the microarray gene-expression datasets of ADAR1 KO livers reveled a type-I interferons signature and an enrichment for immune response genes compared to control littermate livers. Furthermore, we found that in vitro silencing of ADAR1 expression in HepG2 cells leads to enhanced transcription of NFκB target genes, foremost of the pro-inflammatory cytokines IL6 and IL8. We also discovered immune cell-independent paracrine signaling among ADAR1-depleted HepG2 cells and hepatic stellate cells, leading to the activation of the latter cell type to adopt a profibrogenic phenotype. This paracrine communication dependent mainly on the production and secretion of the cytokine IL6 induced by ADAR1 silencing in hepatocytes. Thus, our findings shed a new light on the vital regulatory role of ADAR1 in hepatic immune homeostasis, chiefly its inhibitory function on the crosstalk between the NFκB and type-I interferons signaling cascades, restraining the development of liver inflammation and fibrosis.

A role for APETALA1/fruitfull transcription factors in tomato leaf development.

Flexible maturation rates underlie part of the diversity of leaf shape, and tomato (Solanum lycopersicum) leaves are compound due to prolonged organogenic activity of the leaf margin. The CINCINNATA-teosinte branched1, cycloidea, PCF (CIN-TCP) transcription factor lanceolate (LA) restricts this organogenic activity and promotes maturation. Here, we show that tomato APETALA1/fruitfull (AP1/FUL) MADS box genes are involved in tomato leaf development and are repressed by LA. AP1/FUL expression is correlated negatively with LA activity and positively with the organogenic activity of the leaf margin. LA binds to the promoters of the AP1/FUL genes MBP20 and TM4. Overexpression of MBP20 suppressed the simple-leaf phenotype resulting from upregulation of LA activity or from downregulation of class I knotted like homeobox (KNOXI) activity. Overexpression of a dominant-negative form of MBP20 led to leaf simplification and partly suppressed the increased leaf complexity of plants with reduced LA activity or increased KNOXI activity. Tomato plants overexpressing miR319, a negative regulator of several CIN-TCP genes including LA, flower with fewer leaves via an SFT-dependent pathway, suggesting that miR319-sensitive CIN-TCPs delay flowering in tomato. These results identify a role for AP1/FUL genes in vegetative development and show that leaf and plant maturation are regulated via partially independent mechanisms.

A module of negative feedback regulators defines growth factor signaling.

Signaling pathways invoke interplays between forward signaling and feedback to drive robust cellular response. In this study, we address the dynamics of growth factor signaling through profiling of protein phosphorylation and gene expression, demonstrating the presence of a kinetically defined cluster of delayed early genes that function to attenuate the early events of growth factor signaling. Using epidermal growth factor receptor signaling as the major model system and concentrating on regulation of transcription and mRNA stability, we demonstrate that a number of genes within the delayed early gene cluster function as feedback regulators of immediate early genes. Consistent with their role in negative regulation of cell signaling, genes within this cluster are downregulated in diverse tumor types, in correlation with clinical outcome. More generally, our study proposes a mechanistic description of the cellular response to growth factors by defining architectural motifs that underlie the function of signaling networks.

MiR-377 targets E2F3 and alters the NF-kB signaling pathway through MAP3K7 in malignant melanoma.

BACKGROUND: The incidence of cutaneous malignant melanoma continues to rise, and once the disease metastasizes it is almost inevitably fatal. We recently reported that a large miRNAs cluster on human chromosome 14q32, implicated in many types of cancers, is significantly down-regulated in melanoma. miR-377, one of the miRNAs located within this cluster, was studied here. METHODS: qRT-pCR was used to quantify miR-377 levels in melanoma cell lines and samples. Melanoma cell lines ectopically expressing miR-377 were generated by stable transfection, mRNA expression was assessed using mRNA arrays and protein expression was assessed by Western blot analysis. Potential targets of miR-377 were identified through luciferase reporter assays. Cellular proliferation, migration and soft-agar colony formation were monitored in control and miR-377-expressing cells using cell biology techniques. RESULTS: miR-377 is expressed in normal melanocytes but not in melanoma cell lines or samples. Its ectopic stable expression in melanoma cell lines decreased their proliferative and migratory capacity and their colony-forming capability. mRNA arrays of melanoma cells over-expressing miR-377 pointed to several down-regulated mRNAs that have putative binding sites for miR-377 in their 3'UTR, of which both E2F3 and MAP3K7 were found to be direct targets of miR-377. E2F3, a potent transcriptional inducer of cell-cycle progression, was found to be elevated in melanoma cell lines, but decreased following ectopic expression of miR-377. Ectopic miR-377 also led to a decrease in the activity of a reporter plasmid containing three E2F DNA-binding sites linked to a luciferase cDNA sequence, demonstrating that miR-377 down-regulates E2F3-induced transcription. MAP3K7 (known as TAK1), a serine/threonine kinase along the MAPK signaling pathway, was over-expressed in melanoma but decreased following ectopic expression of miR-377. MAP3K7 is involved in the activation of NF-κB. MiR-377 over-expression led to decreased activity of a reporter plasmid containing two NF-κB DNA-binding sites and to decreased output along the NF-kB signaling pathway. CONCLUSION: Our results suggest that miR-377 is an important negative regulator of E2F and MAP3K7/NF-kB signaling pathway in melanoma cells; it is tempting to speculate that its silencing in melanoma promotes the tumorigenic and metastatic potential of the cells through activation of these pathways.

Global regulation of alternative splicing by adenosine deaminase acting on RNA (ADAR).

Alternative mRNA splicing is a major mechanism for gene regulation and transcriptome diversity. Despite the extent of the phenomenon, the regulation and specificity of the splicing machinery are only partially understood. Adenosine-to-inosine (A-to-I) RNA editing of pre-mRNA by ADAR enzymes has been linked to splicing regulation in several cases. Here we used bioinformatics approaches, RNA-seq and exon-specific microarray of ADAR knockdown cells to globally examine how ADAR and its A-to-I RNA editing activity influence alternative mRNA splicing. Although A-to-I RNA editing only rarely targets canonical splicing acceptor, donor, and branch sites, it was found to affect splicing regulatory elements (SREs) within exons. Cassette exons were found to be significantly enriched with A-to-I RNA editing sites compared with constitutive exons. RNA-seq and exon-specific microarray revealed that ADAR knockdown in hepatocarcinoma and myelogenous leukemia cell lines leads to global changes in gene expression, with hundreds of genes changing their splicing patterns in both cell lines. This global change in splicing pattern cannot be explained by putative editing sites alone. Genes showing significant changes in their splicing pattern are frequently involved in RNA processing and splicing activity. Analysis of recently published RNA-seq data from glioblastoma cell lines showed similar results. Our global analysis reveals that ADAR plays a major role in splicing regulation. Although direct editing of the splicing motifs does occur, we suggest it is not likely to be the primary mechanism for ADAR-mediated regulation of alternative splicing. Rather, this regulation is achieved by modulating trans-acting factors involved in the splicing machinery.

Genome-scale expression and transcription factor binding profiles reveal therapeutic targets in transgenic ERG myeloid leukemia.

The ETS transcription factor ERG plays a central role in definitive hematopoiesis, and its overexpression in acute myeloid leukemia (AML) is associated with a stem cell signature and poor prognosis. Yet how ERG causes leukemia is unclear. Here we show that pan-hematopoietic ERG expression induces an early progenitor myeloid leukemia in transgenic mice. Integrated genome-scale analysis of gene expression and ERG binding profiles revealed that ERG activates a transcriptional program similar to human AML stem/progenitor cells and to human AML with high ERG expression. This transcriptional program was associated with activation of RAS that was required for leukemia cells growth in vitro and in vivo. We further show that ERG induces expression of the Pim1 kinase oncogene through a novel hematopoietic enhancer validated in transgenic mice and human CD34(+) normal and leukemic cells. Pim1 inhibition disrupts growth and induces apoptosis of ERG-expressing leukemic cells. The importance of the ERG/PIM1 axis is further underscored by the poorer prognosis of AML highly expressing ERG and PIM1. Thus, integrative genomic analysis demonstrates that ERG causes myeloid progenitor leukemia characterized by an induction of leukemia stem cell transcriptional programs. Pim1 and the RAS pathway are potential therapeutic targets of these high-risk leukemias.

Systematic identification of rhythmic genes reveals camk1gb as a new element in the circadian clockwork.

A wide variety of biochemical, physiological, and molecular processes are known to have daily rhythms driven by an endogenous circadian clock. While extensive research has greatly improved our understanding of the molecular mechanisms that constitute the circadian clock, the links between this clock and dependent processes have remained elusive. To address this gap in our knowledge, we have used RNA sequencing (RNA-seq) and DNA microarrays to systematically identify clock-controlled genes in the zebrafish pineal gland. In addition to a comprehensive view of the expression pattern of known clock components within this master clock tissue, this approach has revealed novel potential elements of the circadian timing system. We have implicated one rhythmically expressed gene, camk1gb, in connecting the clock with downstream physiology of the pineal gland. Remarkably, knockdown of camk1gb disrupts locomotor activity in the whole larva, even though it is predominantly expressed within the pineal gland. Therefore, it appears that camk1gb plays a role in linking the pineal master clock with the periphery.

Pharyngeal mesoderm regulatory network controls cardiac and head muscle morphogenesis.

The search for developmental mechanisms driving vertebrate organogenesis has paved the way toward a deeper understanding of birth defects. During embryogenesis, parts of the heart and craniofacial muscles arise from pharyngeal mesoderm (PM) progenitors. Here, we reveal a hierarchical regulatory network of a set of transcription factors expressed in the PM that initiates heart and craniofacial organogenesis. Genetic perturbation of this network in mice resulted in heart and craniofacial muscle defects, revealing robust cross-regulation between its members. We identified Lhx2 as a previously undescribed player during cardiac and pharyngeal muscle development. Lhx2 and Tcf21 genetically interact with Tbx1, the major determinant in the etiology of DiGeorge/velo-cardio-facial/22q11.2 deletion syndrome. Furthermore, knockout of these genes in the mouse recapitulates specific cardiac features of this syndrome. We suggest that PM-derived cardiogenesis and myogenesis are network properties rather than properties specific to individual PM members. These findings shed new light on the developmental underpinnings of congenital defects.

Blastocyst implantation failure relates to impaired translational machinery gene expression.

Oocyte quality is a well-established determinant of embryonic fate. However, the molecular participants and biological markers that affect and may predict adequate embryonic development are largely elusive. Our aim was to identify the components of the oocyte molecular machinery that part take in the production of a healthy embryo. For this purpose, we used an animal model, generated by us previously, the oocytes of which do not express Cx43 (Cx43(del/del)). In these mice, oogenesis appears normal, fertilisation does occur, early embryonic development is successful but implantation fails. We used magnetic resonance imaging analysis combined with histological examination to characterise the embryonic developmental incompetence. Reciprocal embryo transfer confirmed that the blastocyst evolved from the Cx43(del/del) oocyte is responsible for the implantation disorder. In order to unveil the genes, the impaired expression of which brings about the development of defective embryos, we carried out a genomic screening of both the oocytes and the resulting blastocysts. This microarray analysis revealed a low expression of Egr1, Rpl21 and Eif4a1 in Cx43(del/del) oocytes and downregulation of Rpl15 and Eif4g2 in the resulting blastocysts. We propose that global deficiencies in genes related to the expression of ribosomal proteins and translation initiation factors in apparently normal oocytes bring about accumulation of defects, which significantly compromise their developmental capacity. The blastocysts resulting from such oocytes, which grow within a confined space until implantation, may be unable to generate enough biological mass to allow their expansion. This information could be implicated to diagnosis and treatment of infertility, particularly to IVF.

A reciprocal tensin-3-cten switch mediates EGF-driven mammary cell migration.

Cell migration driven by the epidermal growth factor receptor (EGFR) propels morphogenesis and involves reorganization of the actin cytoskeleton. Although de novo transcription precedes migration, transcript identity remains largely unknown. Through their actin-binding domains, tensins link the cytoskeleton to integrin-based adhesion sites. Here we report that EGF downregulates tensin-3 expression, and concomitantly upregulates cten, a tensin family member that lacks the actin-binding domain. Knockdown of cten or tensin-3, respectively, impairs or enhances mammary cell migration. Furthermore, cten displaces tensin-3 from the cytoplasmic tail of integrin beta1, thereby instigating actin fibre disassembly. In invasive breast cancer, cten expression correlates not only with high EGFR and HER2, but also with metastasis to lymph nodes. Moreover, treatment of inflammatory breast cancer patients with an EGFR/HER2 dual-specificity kinase inhibitor significantly downregulated cten expression. In conclusion, a transcriptional tensin-3-cten switch may contribute to the metastasis of mammary cancer.

Genome-wide activation of latent donor splice sites in stress and disease.

Sequences that conform to the 5' splice site (5'SS) consensus are highly abundant in mammalian introns. Most of these sequences are preceded by at least one in-frame stop codon; thus, their use for splicing would result in pre-maturely terminated aberrant mRNAs. In normally grown cells, such intronic 5'SSs appear not to be selected for splicing. However, under heat shock conditions aberrant splicing involving such latent 5'SSs occurred in a number of specific gene transcripts. Using a splicing-sensitive microarray, we show here that stress-induced (e.g. heat shock) activation of latent splicing is widespread across the human transcriptome, thus highlighting the possibility that latent splicing may underlie certain diseases. Consistent with this notion, our analyses of data from the Gene Expression Omnibus (GEO) revealed widespread activation of latent splicing in cells grown under hypoxia and in certain cancers such as breast cancer and gliomas. These changes were found in thousands of transcripts representing a wide variety of functional groups; among them are genes involved in cell proliferation and differentiation. The GEO analysis also revealed a set of gene transcripts in oligodendroglioma, in which the level of activation of latent splicing increased with the severity of the disease.

Genotoxicity Associated with Retroviral CAR Transduction of ATM-Deficient T Cells.

Somatic variants in DNA damage response genes such as ATM are widespread in hematologic malignancies. ATM protein is essential for double-strand DNA break repair. Germline ATM deficiencies underlie ataxia-telangiectasia (A-T), a disease manifested by radiosensitivity, immunodeficiency, and predisposition to lymphoid malignancies. Patients with A-T diagnosed with malignancies have poor tolerance to chemotherapy or radiation. In this study, we investigated chimeric antigen receptor (CAR) T cells using primary T cells from patients with A-T (ATM-/-), heterozygote donors (ATM+/-), and healthy donors. ATM-/- T cells proliferate and can be successfully transduced with CARs, though functional impairment of ATM-/- CAR T-cells was observed. Retroviral transduction of the CAR in ATM-/- T cells resulted in high rates of chromosomal lesions at CAR insertion sites, as confirmed by next-generation long-read sequencing. This work suggests that ATM is essential to preserve genome integrity of CAR T-cells during retroviral manufacturing, and its lack poses a risk of chromosomal translocations and potential leukemogenicity. Significance: CAR T-cells are clinically approved genetically modified cells, but the control of genome integrity remains largely uncharacterized. This study demonstrates that ATM deficiency marginally impairs CAR T-cell function and results in high rates of chromosomal aberrations after retroviral transduction, which may be of concern in patients with DNA repair deficiencies.

Detection of BCR-ABL1 mutations in chronic myeloid leukaemia by massive parallel sequencing.

The ability to sequence nucleic acids at an unprecedented pace and decreased costs using massive parallel sequencing (MPS) strongly affects biomedical research. Here we applied MPS for the detection of rare, clinically relevant mutations in a chronic myeloid leukaemia (CML) patient. Tyrosine kinase inhibitors revolutionized CML therapy but in some patients the disease progresses due to resistance-conferring mutations. MPS was applied herein to monitor such mutations in BCR-ABL1 transcripts at different time points. The large volume of sequencing data increases sensitivity compared to direct sequencing and allows detection of marginally represented and previously uncharacterized mutations. We detected changes in the frequency of mutated clones including the emergence and disappearance of the resistance-associated ABL1 T315I mutation. We also observed correlation in appearance of adjacent mutations, and exploited this observation to demonstrate the existence of mutated clones at the time of diagnosis. A tool is provided for detection of low frequency single nucleotide variants/mutations from deep coverage MPS data, applicable to clinical translation of advanced sequencing technologies.