Clinical Significance of Apurinic/Apyrimidinic Endodeoxyribonuclease 1 and MicroRNA Axis in Hepatocellular Carcinoma

Background and Aims

Identification of prognostic factors for hepatocellular carcinoma (HCC) opens new perspectives for therapy. Circulating and cellular onco-miRNAs are noncoding RNAs which can control the expression of genes involved in oncogenesis through post-transcriptional mechanisms. These microRNAs (miRNAs) are considered novel prognostic and predictive factors in HCC. The apurinic/apyrimidinic endodeoxyribonuclease 1 (APE1) contributes to the quality control and processing of specific onco-miRNAs and is a negative prognostic factor in several tumors. The present work aims to: a) define APE1 prognostic value in HCC; b) identify miRNAs regulated by APE1 and their relative target genes and c) study their prognostic value.

Methods

We used The Cancer Genome Atlas (commonly known as TCGA) data analysis to evaluate the expression of APE1 in HCC. To identify differentially-expressed miRNAs (DEmiRNAs) upon APE1 depletion through specific small interfering RNA, we used NGS and nanostring approaches in the JHH-6 HCC tumor cell line. Bioinformatics analyses were performed to identify signaling pathways involving APE1-regulated miRNAs. Microarray analysis was performed to identify miRNAs correlating with serum APE1 expression.

Results

APE1 is considerably overexpressed in HCC tissues compared to normal liver, according to the TCGA-liver HCC (known as LIHC) dataset. Enrichment analyses showed that APE1-regulated miRNAs are implicated in signaling and metabolic pathways linked to cell proliferation, transformation, and angiogenesis, identifying Cyclin Dependent Kinase 6 and Lysosomal Associated Membrane Protein 2 as targets. miR-33a-5p, miR-769, and miR-877 are related to lower overall survival in HCC patients. Through array profiling, we identified eight circulating DE-miRNAs associated with APE1 overexpression. A training phase identified positive association between sAPE1 and miR-3180-3p and miR-769.

Conclusions

APE1 regulates specific miRNAs having prognostic value in HCC.

Shallow Whole-Genome Sequencing of Aedes japonicus and Aedes koreicus from Italy and an Updated Picture of Their Evolution Based on Mitogenomics and Barcoding

Aedes japonicus and Aedes koreicus are two invasive mosquitoes native to East Asia that are quickly establishing in temperate regions of Europe. Both species are vectors of arboviruses, but we currently lack a clear understanding of their evolution. Here, we present new short-read, shallow genome sequencing of A. japonicus and A. koreicus individuals from northern Italy, which we used for downstream phylogenetic and barcode analyses. We explored associated microbial DNA and found high occurrences of Delftia bacteria in both samples, but neither Asaia nor Wolbachia. We then assembled complete mitogenomes and used these data to infer divergence times estimating the split of A. japonicus from A. koreicus in the Oligocene, which was more recent than that previously reported using mitochondrial markers. We recover a younger age for most other nodes within Aedini and other Culicidae. COI barcoding and phylogenetic analyses indicate that A. japonicus yaeyamensis, A. japonicus amamiensis, and the two A. koreicus sampled from Europe should be considered as separate species within a monophyletic species complex. Our studies further clarify the evolution of A. japonicus and A. koreicus, and indicate the need to obtain whole-genome data from putative species in order to disentangle their complex patterns of evolution.

Alpha-1 Adrenergic Antagonists Sensitize Neuroblastoma to Therapeutic Differentiation

Neuroblastoma (NB) is an aggressive childhood tumor, with high-risk cases having a 5-year overall survival probability of approximately 50%. The multimodal therapeutic approach for NB includes treatment with the retinoid isotretinoin (13-cis retinoic acid; 13cRA), which is used in the post-consolidation phase as an antiproliferation and prodifferentiation agent to minimize residual disease and prevent relapse. Through small-molecule screening, we identified isorhamnetin (ISR) as a synergistic compound with 13cRA in inhibiting up to 80% of NB cell viability. The synergistic effect was accompanied by a marked increase in the expression of the adrenergic receptor α1B (ADRA1B) gene. Genetic knockout of ADRA1B or its specific blockade using α1/α1B adrenergic antagonists led to selective sensitization of MYCN-amplified NB cells to cell viability reduction and neural differentiation induced by 13cRA, thus mimicking ISR activity. Administration of doxazosin, a safe α1-antagonist used in pediatric patients, in combination with 13cRA in NB xenografted mice exerted marked control of tumor growth, whereas each drug alone was ineffective. Overall, this study identified the α1B adrenergic receptor as a pharmacologic target in NB, supporting the evaluation of adding α1-antagonists to the post-consolidation therapy of NB to more efficiently control residual disease.

SIGNIFICANCE

Targeting α-adrenergic receptors synergizes with isotretinoin to suppress growth and to promote differentiation of neuroblastoma, revealing a combinatorial approach for more effective management of the disease and prevention of relapse.

Next Generation Sequencing Analysis of MODY-X Patients: A Case Report Series

BACKGROUND

Classic criteria for a maturity-onset diabetes of the young (MODY) diagnosis are often unable to identify all subjects, and traditional Sanger sequencing, using a candidate gene approach, leads to a high prevalence of missed genetic diagnosis, classified as MODY-X. Next generation sequencing (NGS) panels provide a highly sensitive method even for rare forms.

METHODS

We investigated 28 pediatric subjects suspected for MODY-X, utilizing a 15-gene NGS panel for monogenic diabetes (MD).

RESULTS

NGS detected variants of uncertain significance (VUS), likely pathogenic or pathogenic for rarer subtypes of MODY, in six patients. We found variants in the wolframin gene (WFS1), traditionally not considered in MD genetic screening panels, in three patients; KCNJ11 gene mutation, typically responsible for neonatal diabetes and rarely causing isolated diabetes in adolescents; INS gene mutation; a variant in the HNF1B gene in a young male with diabetes on sulfonylurea treatment.

CONCLUSION

In our cohort, the availability of an NGS panel for MD was determined for the correct identification of MD subtypes in six patients with MODY-X. Our study underlines how a precise diagnosis utilizing NGS may have an impact on the management of different forms of MODY and, thus, lead to a tailored treatment and enable genetic counselling of other family members.

Human dyskerin binds to cytoplasmic H/ACA-box-containing transcripts affecting nuclear hormone receptor dependence

Background

Dyskerin is a nuclear protein involved in H/ACA box snoRNA-guided uridine modification of RNA. In humans, its defective function is associated with cancer development and induces specific post-transcriptional alterations of gene expression. In this study, we seek to unbiasedly identify mRNAs regulated by dyskerin in human breast cancer-derived cells.

Results

We find that dyskerin depletion affects the expression and the association with polysomes of selected mRNA isoforms characterized by the retention of H/ACA box snoRNA-containing introns. These snoRNA retaining transcripts (snoRTs) are bound by dyskerin in the cytoplasm in the form of shorter 3′ snoRT fragments. We then characterize the whole cytoplasmic dyskerin RNA interactome and find both H/ACA box snoRTs and protein-coding transcripts which may be targeted by the snoRTs’ guide properties. Since a fraction of these protein-coding transcripts is involved in the nuclear hormone receptor binding, we test to see if this specific activity is affected by dyskerin. Obtained results indicate that dyskerin dysregulation may alter the dependence on nuclear hormone receptor ligands in breast cancer cells. These results are paralleled by consistent observations on the outcome of primary breast cancer patients stratified according to their tumor hormonal status. Accordingly, experiments in nude mice show that the reduction of dyskerin levels in estrogen-dependent cells favors xenograft development in the absence of estrogen supplementation.

Conclusions

Our work suggests a cytoplasmic function for dyskerin which could affect mRNA post-transcriptional networks relevant for nuclear hormone receptor functions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02746-3.

APE1 controls DICER1 expression in NSCLC through miR-33a and miR-130b

Increasing evidence suggests different, not completely understood roles of microRNA biogenesis in the development and progression of lung cancer. The overexpression of the DNA repair protein apurinic/apyrimidinic endodeoxyribonuclease 1 (APE1) is an important cause of poor chemotherapeutic response in lung cancer and its involvement in onco-miRNAs biogenesis has been recently described. Whether APE1 regulates miRNAs acting as prognostic biomarkers of lung cancer has not been investigated, yet. In this study, we analyzed miRNAs differential expression upon APE1 depletion in the A549 lung cancer cell line using high-throughput methods. We defined a signature of 13 miRNAs that strongly correlate with APE1 expression in human lung cancer: miR-1246, miR-4488, miR-24, miR-183, miR-660, miR-130b, miR-543, miR-200c, miR-376c, miR-218, miR-146a, miR-92b and miR-33a. Functional enrichment analysis of this signature revealed its biological relevance in cancer cell proliferation and survival. We validated DICER1 as a direct functional target of the APE1-regulated miRNA-33a-5p and miR-130b-3p. Importantly, IHC analyses of different human tumors confirmed a negative correlation existing between APE1 and Dicer1 protein levels. DICER1 downregulation represents a prognostic marker of cancer development but the mechanisms at the basis of this phenomenon are still completely unknown. Our findings, suggesting that APE1 modulates DICER1 expression via miR-33a and miR-130b, reveal new mechanistic insights on DICER1 regulation, which are of relevance in lung cancer chemoresistance and cancer invasiveness.

A Community Study of SARS-CoV-2 Detection by RT-PCR in Saliva: A Reliable and Effective Method

Efficient, wide-scale testing for SARS-CoV-2 is crucial for monitoring the incidence of the infection in the community. The gold standard for COVID-19 diagnosis is the molecular analysis of epithelial secretions from the upper respiratory system captured by nasopharyngeal (NP) or oropharyngeal swabs. Given the ease of collection, saliva has been proposed as a possible substitute to support testing at the population level. Here, we used a novel saliva collection device designed to favour the safe and correct acquisition of the sample, as well as the processivity of the downstream molecular analysis. We tested 1003 nasopharyngeal swabs and paired saliva samples self-collected by individuals recruited at a public drive-through testing facility. An overall moderate concordance (68%) between the two tests was found, with evidence that neither system can diagnose the infection in 100% of the cases. While the two methods performed equally well in symptomatic individuals, their discordance was mainly restricted to samples from convalescent subjects. The saliva test was at least as effective as NP swabs in asymptomatic individuals recruited for contact tracing. Our study describes a testing strategy of self-collected saliva samples, which is reliable for wide-scale COVID-19 screening in the community and is particularly effective for contact tracing.

The altered transcriptome of pediatric myelodysplastic syndrome revealed by RNA sequencing

Pediatric myelodysplastic syndrome (PMDS) is a very rare and still poorly characterized disorder. In this work, we identified novel potential targets of PMDS by determining genes with aberrant expression, which can be correlated with PMDS pathogenesis. We identified 291 differentially expressed genes (DEGs) in PMDS patients, comprising genes involved in the regulation of apoptosis and the cell cycle, ribosome biogenesis, inflammation and adaptive immunity. Ten selected DEGs were then validated, confirming the sequencing data. These DEGs will potentially represent new molecular biomarkers and therapeutic targets for PMDS.

ddSeeker: a tool for processing Bio-Rad ddSEQ single cell RNA-seq data

Background

New single-cell isolation technologies are facilitating studies on the transcriptomics of individual cells. Bio-Rad ddSEQ is a droplet-based microfluidic system that, when coupled with downstream Illumina library preparation and sequencing, enables the monitoring of thousands of genes per cell. Sequenced reads show unique features that do not permit the use of freely available tools to perform single cell demultiplexing.

Results

We present ddSeeker, a tool to perform initial processing and quality metrics of reads generated through Bio-Rad ddSEQ/Illumina experiments. Its application to the Illumina test dataset demonstrates that ddSeeker performs better than Illumina BaseSpace software, enabling a higher recovery of valid reads. We also show its utility in the analysis of an in-house dataset including two read sets characterized by low and high sequencing quality. ddSeeker and its source code are available at https://github.com/cgplab/ddSeeker.

Conclusions

ddSeeker is a freely available tool to perform initial processing and quality metrics of reads generated through Bio-Rad ddSEQ/Illumina single cell transcriptomic experiments.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5249-x) contains supplementary material, which is available to authorized users.

Mother-to-Infant Microbial Transmission from Different Body Sites Shapes the Developing Infant Gut Microbiome

The acquisition and development of the infant microbiome are key to establishing a healthy host-microbiome symbiosis. The maternal microbial reservoir is thought to play a crucial role in this process. However, the source and transmission routes of the infant pioneering microbes are poorly understood. To address this, we longitudinally sampled the microbiome of 25 mother-infant pairs across multiple body sites from birth up to 4 months postpartum. Strain-level metagenomic profiling showed a rapid influx of microbes at birth followed by strong selection during the first few days of life. Maternal skin and vaginal strains colonize only transiently, and the infant continues to acquire microbes from distinct maternal sources after birth. Maternal gut strains proved more persistent in the infant gut and ecologically better adapted than those acquired from other sources. Together, these data describe the mother-to-infant microbiome transmission routes that are integral in the development of the infant microbiome.

The new phylogeny of the genus Mycobacterium: The old and the news

BACKGROUND

Phylogenetic studies of bacteria have been based so far either on a single gene (usually the 16S rRNA) or on concatenated housekeeping genes. For what concerns the genus Mycobacterium these approaches support the separation of rapidly and slowly growing species and the clustering of most species in well-defined phylogenetic groups. The advent of high-throughput shotgun sequencing leads us to revise conventional taxonomy of mycobacteria on the light of genomic data. For this purpose we investigated 88 newly sequenced species in addition to 60 retrieved from GenBank and used the Average Nucleotide Identity pairwise scores to reconstruct phylogenetic relationships within this genus.

RESULTS

Our analysis confirmed the separation of slow and rapid growers and the intermediate position occupied by the M. terrae complex. Among the rapid growers, the species of the M. chelonae-abscessus complex belonged to the most ancestral cluster. Other major clades of rapid growers included the species related to M. fortuitum and M. smegmatis and a large grouping containing mostly environmental species rarely isolated from humans. The members of the M. terrae complex appeared as the most ancestral slow growers. Among slow growers two deep branches led to the clusters of species related to M. celatum and M. xenopi and to a large group harboring most of the species more frequently responsible of disease in humans, including the major pathogenic mycobacteria (M. tuberculosis, M. leprae, M. ulcerans). The species previously grouped in the M. simiae complex were allocated in a number of sub-clades; of them, only the one including the species M. simiae identified the real members of this complex. The other clades included also species previously not considered related to M. simiae. The ANI analysis, in most cases supported by Genome to Genome Distance and by Genomic Signature-Delta Difference, showed that a number of species with standing in literature were indeed synonymous.

CONCLUSIONS

Genomic data revealed to be much more informative in comparison with phenotype. We believe that the genomic revolution enabled by high-throughput shotgun sequencing should now be considered in order to revise the conservative approaches still informing taxonomic sciences.

Genomic characterization of Nontuberculous Mycobacteria

Mycobacterium tuberculosis and Mycobacterium leprae have remained, for many years, the primary species of the genus Mycobacterium of clinical and microbiological interest. The other members of the genus, referred to as nontuberculous mycobacteria (NTM), have long been underinvestigated. In the last decades, however, the number of reports linking various NTM species with human diseases has steadily increased and treatment difficulties have emerged. Despite the availability of whole genome sequencing technologies, limited effort has been devoted to the genetic characterization of NTM species. As a consequence, the taxonomic and phylogenetic structure of the genus remains unsettled and genomic information is lacking to support the identification of these organisms in a clinical setting. In this work, we widen the knowledge of NTMs by reconstructing and analyzing the genomes of 41 previously uncharacterized NTM species. We provide the first comprehensive characterization of the genomic diversity of NTMs and open new venues for the clinical identification of opportunistic pathogens from this genus.

Abstract 2883: Impact of novel CDKN2A/p16INK4a 5’UTR variants predisposing to melanoma on p16 translational regulation

CDKN2A/p16INK4a is an important tumor-suppressor gene whose dysregulation is associated with melanoma. We have recently demonstrated that the p16INK4a expression can be modulated by IRES-dependent mRNA translation and this regulation might have an important role during tumorigenesis. We also identified YBX1 as RNA binding protein targeting and stimulating p16 mRNA translation efficiency. The identified IRES function was particularly active in hypoxia and inhibition of mTOR. This post-transcriptional regulatory mechanism would potentially be a target for mutational inactivation during melanomagenesis. Indeed, we previously showed that p16INK4a 5’UTR variants found in melanoma patients with a family predisposition can have a negative effect on p16 translation. Single Nucleotide Variants (SNVs) have been investigated during routine testing of CDKN2A/p16INK4a in Italian, English and French melanoma patients followed by GenoMEL, the International Melanoma Genetics Consortium and a total of 17 germline variants were identified in a cohort of nearly 6000 patients. These p16INK4a 5’UTR 17 variants were studied with multiple approaches, that included mono- and bi-cistronic reporter assays, western blot of endogenous protein, and quantification of allelic representation after polysomal profiling to investigate their impact on p16INK4a mRNA translation regulation. We devised a classification score based on the concordance between functional assays and the extent of dysfunction displayed by each variant compared to the wild type. Variants are classified as neutral (score 0) when no difference was observed in at least 3 assays. This applied to: c.-14C>T, c.-20A>G, c.-25C>T+c.-180G>A, c.-30G>A, c.-40C>T, c.-45G>A, c.-59C>G, c.-87T>A, c.-252A>T. Variants were considered as potential mutations when a defect was measured in either one or two assays (score 1-2; c.-42T>A and c.-67G>C variants). Finally, we classified variants as causal mutations when three or more assays showed impairment (score >3). This applied to: c.-27del23, c.-56G>T, c.-93-91delAGG, as well as to c.-21C>T and c.-34G>T, which were already considered as a causal mutations. We have also determined the structure of the wild type p16INK4a 5’UTR by Selective Hydroxyl Acylation or SHAPE assay. The variant c.-42T>A encompassing the predicted YBX1 binding site was also studied and shown to induce a local change in conformation. The structure of all p16INK4a 5’UTR variants examined so far is being investigated as well as their impact on the interaction of RNA binding proteins such as YBX1, SRSF1 and RBM4. Our data indicate that the sequencing of the entire p16INK4a 5’UTR should be included in routine screening of melanoma families as nearly half of SNVs tested so far displayed a negative impact on the p16 mRNA translation efficiency.

Citation Format: Alessandra Bisio, Elisa Latorre, Virginia Andreotti, Brigitte Bressac-de Paillerets, Mark Harland, Odile Cabaret, Julia Newton-Bishop, Lorenza Pastorino, William Bruno, Roberto Bertorelli, Veronica De Sanctis, Chiara Menin, Gilberto Fronza, Paola Queirolo, Giovanna Bianchi Scarrà, Robert C. Spitale, Alessandro Provenzani, Alberto Inga, Paola Ghiorzo. Impact of novel CDKN2A/p16INK4a 5’UTR variants predisposing to melanoma on p16 translational regulation. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 2883.

The CDKN2A/p16(INK) (4a) 5'UTR sequence and translational regulation: impact of novel variants predisposing to melanoma

Many variants of uncertain functional significance in cancer susceptibility genes lie in regulatory regions, and clarifying their association with disease risk poses significant challenges. We studied 17 germline variants (nine of which were novel) in the CDKN2A 5'UTR with independent approaches, which included mono and bicistronic reporter assays, Western blot of endogenous protein, and allelic representation after polysomal profiling to investigate their impact on CDKN2A mRNA translation regulation. Two of the novel variants (c.-27del23, c.-93-91delAGG) were classified as causal mutations (score ≥3), along with the c.-21C>T, c.-34G>T, and c.-56G>T, which had already been studied by a subset of assays. The novel c.-42T>A as well as the previously described c.-67G>C were classified as potential mutations (score 1 or 2). The remaining variants (c.-14C>T, c.-20A>G, c.-25C>T+c.-180G>A, c.-30G>A, c.-40C>T, c.-45G>A, c.-59C>G, c.-87T>A, c.-252A>T) were classified as neutral (score 0). In conclusion, we found evidence that nearly half of the variants found in this region had a negative impact on CDKN2A mRNA translation, supporting the hypothesis that 5'UTR can act as a cellular Internal Ribosome Entry Site (IRES) to modulate p16(INK) (4a) translation.

Characterization of 17 strains belonging to the Mycobacterium simiae complex and description of Mycobacterium paraense sp. nov

Fourteen mycobacterial strains isolated from pulmonary samples of independent patients in the state of Pará (Brazil), and three strains isolated in Italy, were characterized using a polyphasic approach. Thorough genetic investigation, including whole-genome sequencing, demonstrated that the strains belong to the M. simiae complex, being most closely related to Mycobacterium interjectum. For 14 of the strains, evidence emerged supporting their inclusion in a previously unreported species of the genus Mycobacterium, for which the name Mycobacterium paraense sp. nov. is proposed (type strain, IEC26(T) = DSM 46749(T) = CCUG 66121(T)). The novel species is characterized by slow growth, unpigmented or pale yellow scotochromogenic colonies, and a HPLC mycolic acid profile different from other known mycobacteria. In different genetic regions, high sequence microheterogeneity was detected.

OncomiR detection in circulating body fluids: a PDMS microdevice perspective

There is an increasing interest in circulating microRNAs (miRNAs) as potential minimally invasive diagnostic biomarkers in oncology. Considerable efforts are being made in the development of lab-on-a-chip devices for biomedical applications to purify and detect miRNAs from biological fluids. Here, we report the development of an innovative polydimethylsiloxane (PDMS)-based parallel device whose internal surface can opportunely be functionalized with positively charged 3-aminopropyltriethoxysilane (APTES) alone or mixed with two different neutral poly(ethylene glycol) silanes (PEG-s). The differently functionalized internal surfaces of the PDMS chip were characterized with s-SDTB (sulfosuccinimidyl-4-o-(4,4-dimethoxytrityl) butyrate) and the portion of the surface able to adsorb a synthetic fluorescently labeled miRNA was determined. Interestingly, the adsorbed miRNA (both synthetic and cell supernatant-derived) was found mainly on the bottom surface of the chip and could be reverse transcribed into cDNA directly on the same PDMS chip used for its purification, saving hours with respect to the use of standard purification kits. We identified 0.1% APTES/0.9% PEG-silane as the most efficient PDMS functionalization to capture both synthetic and extracellular miRNA. Moreover, the amount of captured miRNA was increased by treating the cell supernatant with a commercially available lysis buffer for RNA extraction. We assessed that the available miRNA binding sites on the functionalized surface were efficiently saturated with only one incubation, shortening the time and greatly simplifying the protocol for miRNA purification from biological samples. Finally, the extracellular miRNA purification efficiency of the PDMS functionalized multichip determined via real-time quantitative polymerase chain reaction (RT-qPCR) was confirmed by droplet digital PCR (ddPCR) quantification. This work shows an innovative, rapid and easy to use microdevice for the purification and reverse transcription of circulating miRNAs, approaching the realization of diagnostic and prognostic oncomiR-based assays.

Enhanced microbial diversity in the saliva microbiome induced by short-term probiotic intake revealed by 16S rRNA sequencing on the IonTorrent PGM platform

Microbial communities populating several human body habitats are important determinants of human health. Cultivation-free community-wide approaches like bacterial 16S rRNA sequencing recently revolutionized the study of such human-associated microbial diversity, and the continuously decreasing cost/throughput ratio of current sequencing platforms is further enhancing the availability and effectiveness of microbiome research. The IonTorrent PGM platform is among the latest available commercial high-throughput sequencing tools, but it is just starting to be used for 16S rRNA surveys with only episodic assessments of its performance. We present here the first IonTorrent profiling of the human saliva microbiome collected from 12 healthy individuals. In this cohort, a subset of the volunteers was asked to assume a probiotic product, in order to investigate its impact on the composition and the structure of the saliva microbiome. Analysis of the generated dataset suggests the suitability of the IonTorrent platform for 16S rRNA surveys, even though some platform-specific choices are required to optimize the consistency of the obtained bacterial profiles. Interestingly, we found a marked and statistically significant increase of the overall bacterial diversity in the saliva of individuals who received the probiotic product compared to the control group, suggesting a short-term effect of probiotic product administration on oral microbiome composition.

Identification of new p53 target microRNAs by bioinformatics and functional analysis

BACKGROUND

The tumor suppressor p53 is a sequence-specific transcription factor that regulates an extensive network of coding genes, long non-coding RNAs and microRNAs, that establish intricate gene regulatory circuits influencing many cellular responses beyond the prototypical control of cell cycle, apoptosis and DNA repair.

METHODS

Using bioinformatic approaches, we identified an additional group of candidate microRNAs (miRs) under direct p53 transcriptional control. To validate p53 family-mediated responsiveness of the newly predicted target miRs we first evaluated the potential for wild type p53, p63β and p73β to transactivate from p53 response elements (REs) mapped in the miR promoters, using an established yeast-based assay.

RESULTS

The REs found in miR-10b, -23b, -106a, -151a, -191, -198, -202, -221, -320, -1204, -1206 promoters were responsive to p53 and 8 of them were also responsive to p63β or p73β. The potential for germline p53 mutations to drive transactivation at selected miR-associated REs was also examined. Chromatin Immuno-Precipitation (ChIP) assays conducted in doxorubicin-treated MCF7 cells and HCT116 p53+/+ revealed moderate induction of p53 occupancy at the miR-202, -1204, -1206, -10b RE-containing sites, while weak occupancy was observed for the miR-23b-associated RE only in MCF7 cells. RT-qPCR analyses cells showed modest doxorubicin- and/or Nutlin-dependent induction of the levels of mature miR-10b, -23b, -151a in HCT116 p53+/+ and MCF7 cells. The long noncoding RNA PVT1 comprising miR-1204 and -1206 was weakly induced only in HCT116 p53+/+ cells, but the mature miRs were not detected. miR-202 expression was not influenced by p53-activating stimuli in our cell systems.

CONCLUSIONS

Our study reveals additional miRs, particularly miR-10b and miR-151a, that could be directly regulated by the p53-family of transcription factors and contribute to the tuning of p53-induced responses.

Interaction between p53 and estradiol pathways in transcriptional responses to chemotherapeutics

Estrogen receptors (ERs) and p53 can interact via cis-elements to regulate the angiogenesis-related VEGFR-1 (FLT1) gene, as we reported previously. Here, we address cooperation between these transcription factors on a global scale. Human breast adenocarcinoma MCF7 cells were exposed to single or combinatorial treatments with the chemotherapeutic agent doxorubicin and the ER ligand 17β-estradiol (E2). Whole-genome transcriptome changes were measured by expression microarrays. Nearly 200 differentially expressed genes were identified that showed limited responsiveness to either doxorubicin treatment or ER ligand alone but were upregulated in a greater than additive manner following combined treatment. Based on exposure to 5-fuorouracil and nutlin-3a, the combined responses were treatment-specific. Among 16 genes chosen for validation using quantitative real-time PCR, seven (INPP5D, TLR5, KRT15, EPHA2, GDNF, NOTCH1, SOX9) were confirmed to be novel direct targets of p53, based on responses in MCF7 cells silenced for p53 or cooperative targets of p53 and ER. Promoter pattern searches and chromatin IP experiments for the INPP5D, TLR5, KRT15 genes supported direct, cis-mediated p53 and/or ER regulation through canonical and noncanonical p53 and ER response elements. Collectively, we establish that combinatorial activation of p53 and ER can induce novel gene expression programs that have implications for cell-cell communications, adhesion, cell differentiation, development and inflammatory responses as well as cancer treatments.

Abstract 2286: p53-miR-dependent post-transcriptional circuits: mechanisms, targets and inter-individual variation

Using bioinformatics approaches, we identified a group of candidate microRNAs (miRs) for direct p53 transcriptional control. To validate p53 family-mediated control of the newly predicted target miRs we evaluated the potential for wild type p53, p63α and p73α to transactivate from p53 response elements (REs) mapped at the miR promoters. Results obtained with an established yeast-based assay indicated that the REs found at miR-10b, -23b, -106a, -151, -191, -198, -202, -221, -320, -1204, -1026 genes were responsive to p53 of which 8 were also responsive to p63α or p73α. ChIP assays in doxorubicin-treated HCT116 cells -p53+/+ vs p53-/– revealed moderate p53 occupancy at the miR-202, -1206 sites, and weak occupancy at miR-10b, -191. RT-qPCR analyses in HCT116 and MCF7 cells showed modest doxorubicin- and/or p53-dependent regulation for miR-23b, -151, -191. We then focused on miR-191 and its potential impact on MDM4, an important modulator of p53. The A>C SNP (rs4245739) in the proximal portion of the MDM4 3′UTR occurs in a predicted binding site for miR-191, with the A allele resulting in a mismatch at the predicted miR-191 seed. This SNP is in strong association with the tag-SNP (rs2369244) within the MDM4 haplotype that was reported to impact on cancer risk. Using MCF7 and HCT116 cells that are heterozygous for rs4245739, we examined at the endogenous MDM4 gene level the impact of varying the expression of miR-191. A slight allelic imbalance was observed by allele-specific qPCR when miR-191 was over-expressed, and the co-transfection of miR-34a, for which seed sequences are predicted in the distal portion of the MDM4-3′UTR, had an additive effect, leading to ∼20% reduction of the C-allele. Western blot analyses indicated that steady-state MDM4 protein levels were reduced by miR-191 and/or miR-34a overexpression. Our study reveals additional miRs that could be directly regulated by the p53-family of transcription factors and could contribute to the tuning of p53-induced responses with the possibility of inter-individual variations due to functional SNPs.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 2286. doi:1538-7445.AM2012-2286

Gene expression pathways induced by axotomy and decentralization of rat superior cervical ganglion neurons

To identify genes potentially involved in remodelling synaptic connections, we induced the temporary detachment of pre- and post-synaptic elements by axotomy or denervation of rat superior cervical ganglion neurons. cDNA microarray analysis followed by stringent selection criteria allowed the identification of a panel of genes whose expression was modulated by axotomy at various time points after injury. Among these genes, 11 were validated by real-time reverse transcriptase-polymerase chain reaction on independently prepared samples after superior cervical ganglion neuron axotomy (1, 3 and 6 days) and compared with the effect of decentralization (8 h, 1 and 3 days). These genes code for extracellular matrix/space [apolipoprotein D (apoD), decorin, collagen alpha1 type I, collagen alpha1 type III] and intermediate filament (vimentin) proteins, for modulators of neurite outgrowth (thrombin receptor, plasminogen activator inhibitor-1, bone morphogenetic protein 4, annexin II and S-100-related protein, clone 42C) and for a nerve cell transcription factor (brain finger protein). Eight of these 11 genes showed significant and persistent modulations after both types of injury. Finally, protein levels of apoD were shown to increase in superior cervical ganglion after axotomy. Our results identify hitherto unrecorded genes responsive to axotomy and decentralization of superior cervical ganglion neurons, and probably involved in synapse formation, remodelling and elimination.

Distinct DNA elements contribute to Rap1p affinity for its binding sites

The essential Saccharomyces cerevisiae regulatory protein Rap1 contains two tandem Myb-like DNA binding sub-domains that interact with two defined DNA "hemisites", separated by a trinucleotide linker sequence. We have mapped the thermodynamically defined DNA-binding site of Rap1 by a primer extension method coupled with electrophoretic separation of bound and unbound DNAs. Relative to published consensus sequences, we detect binding interactions that extend 3 bp beyond the 5'-end of the putative DNA-binding site. This new site of interaction is located where the DNA minor groove faces the protein, and may account for the major DNA bending induced by Rap1p that previous studies have mapped to a site immediately upstream of the consensus binding site. In addition, we show that a minimal DNA-binding site made of one single consensus hemisite, preceded or followed by a spacer trinucleotide that interacts with the unstructured protein linker between the two Rap1p DNA binding domains, is able to bind the protein, although at lower affinity. These findings may explain the observed in vivo binding properties of Rap1p at many promoters that lack canonical binding sites.

In vivo topography of Rap1p-DNA complex at Saccharomyces cerevisiae TEF2 UAS(RPG) during transcriptional regulation

We have analyzed in detail the structure of RAP1-UAS(RPG) complexes in Saccharomyces cerevisiae cells using multi-hit KMnO(4), UV and micrococcal nuclease high-resolution footprinting. Three copies of the Rap1 protein are bound to the promoter simultaneously in exponentially growing cells, as shown by KMnO(4) multi-hit footprinting analysis, causing extended and diagnostic changes in the DNA structure of the region containing the UAS(RPG). Amino acid starvation does not cause loss of Rap1p from the complex; however, in vivo UV-footprinting reveals the occurrence of structural modifications of the complex. Moreover, low-resolution micrococcal nuclease digestion shows that the chromatin of the entire region is devoid of positioned nucleosomes but is susceptible to changes in accessibility to the nuclease upon amino acid starvation. The implications of these results for the mechanism of Rap1p action are discussed.