Adattamento morfologico e variazione dell'espressione genica dopo pneumonectomia: osservazioni su modello sperimentale animale

Introduzione La pneumonectomia su modello animale potrebbe essere un’utile piattaforma di studio per approfondire i meccanismi della risposta compensatoria al danno polmonare. Scopo dello studio è determinare la presenza di variazioni morfologiche e di espressione del trascrittoma dopo pneumonectomia. Materiali e metodi Undici suini sono stati sottoposti a pneumonectomia sinistra. Sono stati eseguiti prelievi sito-specifici intraoperatori su polmone sinistro e successivamente confrontati con prelievi sito-specifici su polmone destro dopo eutanasia a 60 giorni. I prelievi degli animali con decorso regolare sono stati sottoposti a RNA-sequencing e successiva analisi computazionale per valutare il peso funzionale del singolo gene o di clusters di geni. Risultati Un animale è stato escluso per insorgenza di ernia diaframmatica. In 7/10 è stata riscontrata apertura della pleura mediastinica con parziale erniazione del polmone controlaterale e shift mediastinico. L’istologia ha mostrato dilatazione degli spazi aerei, rottura dei setti interalveolari, lieve infiammazione, assenza di fibrosi, stiramento radiale dei bronchi e riduzione del letto capillare. L’analisi di bulk RNA-sequencing ha identificato 553 geni espressi in modo differenziale (DEG)(P<0,001) tra pre e post-pneumonectomia. I primi 10 DEG up-regolati: Edn1, Areg, Havcr2, Gadd45g, Depp1, Cldn4, Atf3, Myc, Gadd45b, Socs3; i primi 10 geni down-regolati: Obscn, Cdkn2b, ENSSSCG00000015738, Prrt2, Amer1, Flrt3, Efnb2, Tox3, Znf793, Znf365. Tra i DEG è stata riscontrata una predominanza di geni specifici dei macrofagi. L’analisi di gene ontology basata su DAVID ha mostrato un significativo arricchimento della "via di segnalazione apoptotica estrinseca"(FDR q=7,60x10 -3), della via di “Risposta all’insulina”(FDR q=7,60x10 -3) ed un arricchimento di geni “Regolatori negativi del segnale DDX58/IFIH1”(FDR q=7.50x10 -4). Conclusioni Il presente studio conferma la presenza di variazioni macroscopiche e microscopiche fenotipiche dopo pneumonectomia. L’RNA sequencing e lo studio di genomica traslazionale hanno mostrato l’esistenza di geni singoli e di network di geni disregolati dopo pneumonectomia, prevalentemente in determinate popolazioni cellulari.

Evaluation of the transcriptomic response of an Alzheimer's disease murine model induced at different ages: searching for predictors

Alzheimer’s disease (AD) is a chronic, progressive neurodegenerative disease, characterized by the impairment of mnesic and cognitive functions, that represents the most frequent type of dementia in older people worldwide. Aging is the most important risk factor for the sporadic form of the pathology and it is associated to the progressive impairment of the proteostasis network. The endoplasmic reticulum (ER), the main cellular actor involved in proteostasis, appears significantly compromised in AD due to the accumulation of β-amyloid (Aβ) protein and phosphorylated-tau protein. Increasing proteins misfolding activates a specific cellular response known as Unfolded Protein response (UPR) which orchestrates the recovery of ER function. The aim of the present study was to investigate the role of UPR and aging process in a murine model of AD induced by intracerebroventricular (i.c.v.) injection of Aβ1-42 oligomers at 3 or 18 months. The oligomers injection in aged animals caused the increased of memory impairment, oxidative stress, and the depletion of glutathione reserve. Furthermore, the RNA-sequencing analysis was performed and the bioinformatic analysis showed the enrichment of several pathways involved in neurodegeneration and protein regulations. The following analysis highlighted the significant dysregulation of the three branches of the UPR, the protein kinase RNA-like ER kinase (PERK), inositol-requiring protein 1α (IRE1α) and activating transcription factor 6 (ATF-6). In turn, ER stress affected the PI3K/Akt/Gsk3β and MAPK/ERK pathways, highlighting Mapkapk5 as a potential marker of the neurodegenerative process, which regulation could lead to the definition of new pharmacological and neuroprotective strategies to counteract AD.

Characterization of the molecular mechanisms mediated by the insulin-like growth factor-2 mRNA-binding proteins and Ephrin receptor A2 (EPHA2) signaling in the malignancy of CIC-DUX4 sarcomas

Ewing sarcoma (EWS) and CIC-DUX4 sarcoma (CDS) are pediatric fusion gene-driven tumors of mesenchymal origin characterized by an extremely stable genome and limited clinical solutions. Post-transcriptional regulatory mechanisms are crucial for understanding the development of this class of tumors. RNA binding proteins (RBPs) play a crucial role in the aggressiveness of these tumors. Numerous RBP families are dysregulated in cancer, including IGF2BPs. Among these, IGF2BP3 is a negative prognostic factor in EWS because it promotes cell growth, chemoresistence, and induces the metastatic process. Based on preliminary RNA sequencing data from clinical samples of EWS vs CDS patients, three major axes that are more expressed in CDS have been identified, two of which are dissected in this PhD work. The first involves the transcription factor HMGA2, IGF2BP2-3, and IGF2; the other involves the ephrin receptor system, particularly EphA2. EphA2 is involved in numerous cellular functions during embryonic stages, and its increased expression in adult tissues is often associated with pathological conditions. In tumors, its role is controversial because it can be associated with both pro- and anti-tumoral mechanisms. In EWS, it has been shown to play a role in promoting cell migration and neoangiogenesis. Our study has confirmed that the HMGA2/IGF2BPs/IGF2 axis contributes to CDS malignancy, and Akt hyperactivation has a strong impact on migration. Using loss/gain of function models for EphA2, we confirmed that it is a substrate of Akt, and Akt hyperactivation in CDS triggers ligand-independent activation of EphA2 through phosphorylation of S897. Moreover, the combination of Trabectedin and NVP/BEZ235 partially inhibits Akt/mTOR activation, resulting in reduced tumor growth in vivo. Inhibition of EphA2 through ALWII 41_27 significantly reduces migration in vitro. The project aim is the identification of target molecules in CDS that can distinguish it from EWS and thus develop new targeted therapeutic strategies.

Cancers of unknown primary site (CUPs) as a model of metastatic process.

Cancers of unknown primary site (CUPs) are a rare group of metastatic tumours, with a frequency of 3-5%, with an overall survival of 6-10 month. The identification of tumour primary site is usually reached by a combination of diagnostic investigations and immunohistochemical testing of the tumour tissue. In CUP patients, these investigations are inconclusive. Since international guidelines for treatment are based on primary site indication, CUP treatment requires a blind approach. As a consequence, CUPs are usually empiric treated with poorly effective. In this study, we applied a set of microRNAs using EvaGreen-based Droplet Digital PCR in a retrospective and prospective collection of formalin-fixed paraffin-embedded tissue samples. We assessed miRNA expression of 155 samples including primary tumours (N=94), metastases of known origin (N=10) and metastases of unknown origin (N=50). Then, we applied the shrunken centroids predictive algorithm to obtain the CUP’s site(s)-of-origin. The molecular test was successfully applied to all CUP samples and provided a site-of-origin identification for all samples, potentially within a one-week time frame from sample inclusion. In the second part of the study we derived two CUP cell lines, and corresponding patient-derived xenografts (PDXs). CUP cell lines and PDXs underwent histological, molecular, and genomic characterization confirming the features of the original tumour. Tissues-of-origin prediction was obtained from the tumour microRNA expression profile and confirmed by single cell RNA sequencing. Genomic testing analysis identified FGFR2 amplification in both models. Drug-screening assays were performed to test the activity of FGFR2-targeting drug and the combination treatment with the MEK inhibitor trametinib, which proved to be synergic and exceptionally active, both in vitro and in vivo. In conclusion, our study demonstrated that miRNA expression profiling could be employed as diagnostic test. Then we successfully derived two CUP models from patients, used for therapy tests, bringing personalized therapy closer to CUP patients.

Comparative analyses of classical phenotypic method and ribosomal RNA gene sequencing for identification of medically relevant Candida species

As the distribution of Candida species and their susceptibility to antifungal agents have changed, a new means of accurately and rapidly identifying these species is necessary for the successful early resolution of infection and the subsequent reduction of morbidity and mortality. The current work aimed to evaluate ribosomal RNA gene sequencing for the identification of medically relevant Candida species in comparison with a standard phenotypic method. Eighteen reference strains (RSs), 69 phenotypically identified isolates and 20 inconclusively identified isolates were examined. Internal transcribed spaces (ITSs) and D1/D2 of the 26S ribosomal RNA gene regions were used as targets for sequencing. Additionally, the sequences of the ITS regions were used to establish evolutionary relationships. The sequencing of the ITS regions was successful for 88% (94/107) of the RS and isolates, whereas 100% of the remaining 12% (13/107) of the samples were successfully analysed by sequencing the D1/D2 region. Similarly, genotypic analysis identified all of the RS and isolates, including the 20 isolates that were not phenotypically identified. Phenotypic analysis, however, misidentified 10% (7/69) of the isolates. Phylogenetic analysis allowed the confirmation of the relationships between evolutionarily close species. Currently, the use of genotypic methods is necessary for the correct identification of Candida species.

Comparative clinical sample preparation of DNA and RNA viral nucleic acids for a commercial deep sequencing system (Illumina MiSeq (R))

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

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

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

Metodologie per l'analisi statistica di dati di sequencing relativi a un esperimento di rna interference

La RNA interference è un processo attraverso il quale alcuni piccoli frammenti di RNA (19-25 nucleotidi) sono in grado di silenziare l'espressione genica. La sua scoperta, nel 1998, ha rivoluzionato le concezioni della biologia molecolare, minando le basi del cosiddetto Dogma Centrale. Si è visto che la RNAi riveste ruoli fondamentali in meccanismi di regolazione genica, nello spegnimento dell'espressione e funziona come meccanismo di difesa innata contro varie tipologie di virus. Proprio a causa di queste implicazioni richiama interesse non solo dal punto di vista scientifico, ma anche da quello medico, in quanto potrebbe essere impiegata per lo sviluppo di nuove cure. Nonostante la scoperta di tale azione desti la curiosità e l'interesse di molti, i vari processi coinvolti, soprattutto a livello molecolare, non sono ancora chiari. In questo lavoro si propongono i metodi di analisi di dati di un esperimento prodotto dall'Istituto di Biologia molecolare e cellulare di Strasburgo. Nell'esperimento in questione vengono studiate le funzioni che l'enzima Dicer-2 ha nel pathway - cioè la catena di reazioni biomolecolari - della RNA interference durante un'infezione virale nel moscerino della frutta Drosophila Melanogaster. Per comprendere in che modo Dicer-2 intervenga nel silenziamento bisogna capire in quali casi e quali parti di RNA vengono silenziate, a seconda del diverso tipo di mutazione dell'enzima stesso. Dunque è necessario sequenziare l'RNA nelle diverse condizioni sperimentali, ottenendo così i dati da analizzare. Parte dei metodi statistici che verranno proposti risultano poco convenzionali, come conseguenza della peculiarità e della difficoltà dei quesiti che l'esperimento mette in luce. Siccome le tematiche affrontate richiedono un approccio sempre più interdisciplinare, è aumentata considerevolmente la richiesta di esperti di altri settori scientifici come matematici, informatici, fisici, statistici e ingegneri. Questa collaborazione, grazie a una diversità di approccio ai problemi, può fornire nuovi strumenti di comprensione in ambiti che, fino a poco tempo fa, rientravano unicamente nella sfera di competenza dei biologi.

Next-generation sequencing of HIV-1 RNA genomes: determination of error rates and minimizing artificial recombination.

Next-generation sequencing (NGS) is a valuable tool for the detection and quantification of HIV-1 variants in vivo. However, these technologies require detailed characterization and control of artificially induced errors to be applicable for accurate haplotype reconstruction. To investigate the occurrence of substitutions, insertions, and deletions at the individual steps of RT-PCR and NGS, 454 pyrosequencing was performed on amplified and non-amplified HIV-1 genomes. Artificial recombination was explored by mixing five different HIV-1 clonal strains (5-virus-mix) and applying different RT-PCR conditions followed by 454 pyrosequencing. Error rates ranged from 0.04-0.66% and were similar in amplified and non-amplified samples. Discrepancies were observed between forward and reverse reads, indicating that most errors were introduced during the pyrosequencing step. Using the 5-virus-mix, non-optimized, standard RT-PCR conditions introduced artificial recombinants in a fraction of at least 30% of the reads that subsequently led to an underestimation of true haplotype frequencies. We minimized the fraction of recombinants down to 0.9-2.6% by optimized, artifact-reducing RT-PCR conditions. This approach enabled correct haplotype reconstruction and frequency estimations consistent with reference data obtained by single genome amplification. RT-PCR conditions are crucial for correct frequency estimation and analysis of haplotypes in heterogeneous virus populations. We developed an RT-PCR procedure to generate NGS data useful for reliable haplotype reconstruction and quantification.

DNA sequencing and genotyping by transcriptional synthesis of chain-terminated RNA ladders and MALDI-TOF mass spectrometry

Sets of RNA ladders can be synthesized by transcription of a bacteriophage-encoded RNA polymerase using 3′-deoxynucleotides as chain terminators. These ladders can be used for sequencing of DNA. Using a nicked form of phage SP6 RNA polymerase in this study substantially enhanced yields of transcriptional sequencing ladders. Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) of chain-terminated RNA ladders allowed DNA sequence determination of up to 56 nt. It is also demonstrated that A→G and C→T variations in heterozygous and homozygous samples can be unambiguously identified by the mass spectrometric analysis. As a step towards single-tube sequencing reactions, α-thiotriphosphate nucleotide analogs were used to overcome problems caused by chain terminator-independent, premature termination and by the small mass difference between natural pyrimidine nucleotides.

A score system for quality evaluation of RNA sequence tags: an improvement for gene expression profiling

Background: High-throughput molecular approaches for gene expression profiling, such as Serial Analysis of Gene Expression (SAGE), Massively Parallel Signature Sequencing (MPSS) or Sequencing-by-Synthesis (SBS) represent powerful techniques that provide global transcription profiles of different cell types through sequencing of short fragments of transcripts, denominated sequence tags. These techniques have improved our understanding about the relationships between these expression profiles and cellular phenotypes. Despite this, more reliable datasets are still necessary. In this work, we present a web-based tool named S3T: Score System for Sequence Tags, to index sequenced tags in accordance with their reliability. This is made through a series of evaluations based on a defined rule set. S3T allows the identification/selection of tags, considered more reliable for further gene expression analysis. Results: This methodology was applied to a public SAGE dataset. In order to compare data before and after filtering, a hierarchical clustering analysis was performed in samples from the same type of tissue, in distinct biological conditions, using these two datasets. Our results provide evidences suggesting that it is possible to find more congruous clusters after using S3T scoring system. Conclusion: These results substantiate the proposed application to generate more reliable data. This is a significant contribution for determination of global gene expression profiles. The library analysis with S3T is freely available at http://gdm.fmrp.usp.br/s3t/.S3T source code and datasets can also be downloaded from the aforementioned website.

Ultra-Deep Sequencing Reveals the microRNA Expression Pattern of the Human Stomach

Background: While microRNAs (miRNAs) play important roles in tissue differentiation and in maintaining basal physiology, little is known about the miRNA expression levels in stomach tissue. Alterations in the miRNA profile can lead to cell deregulation, which can induce neoplasia. Methodology/Principal Findings: A small RNA library of stomach tissue was sequenced using high-throughput SOLiD sequencing technology. We obtained 261,274 quality reads with perfect matches to the human miRnome, and 42% of known miRNAs were identified. Digital Gene Expression profiling (DGE) was performed based on read abundance and showed that fifteen miRNAs were highly expressed in gastric tissue. Subsequently, the expression of these miRNAs was validated in 10 healthy individuals by RT-PCR showed a significant correlation of 83.97% (P<0.05). Six miRNAs showed a low variable pattern of expression (miR-29b, miR-29c, miR-19b, miR-31, miR-148a, miR-451) and could be considered part of the expression pattern of the healthy gastric tissue. Conclusions/Significance: This study aimed to validate normal miRNA profiles of human gastric tissue to establish a reference profile for healthy individuals. Determining the regulatory processes acting in the stomach will be important in the fight against gastric cancer, which is the second-leading cause of cancer mortality worldwide.

Effects of bed materials on the performance of an anaerobic sequencing batch biofilm reactor treating domestic sewage

The objective of this study was to determine the best performance of an anaerobic sequencing batch biofilm reactor (AnSBBR) based on the use of four different bed materials as support for biomass immobilization. The bed materials utilized were Polyurethane foam (PU), vegetal carbon (VC), synthetic pumice (SP), and recycled low-density polyethylene (PE). The AnSBBR. with I total volume Of 7.2 L, was operated in 8-h batch cycles over 10 months, and fed with domestic sewage with an average influent chemical oxygen demand (COD) of 358 +/- 110 mg/L. The average effluent COD values were 121 +/- 31, 208 +/- 54, 233 +/- 52, and 227 +/- 51 mg/L. for PU, VC, SP, and PE, respectively. A modified first-order kinetic model was adjusted to temporal profiles of COD during a batch cycle, and the apparent kinetic constants were 0.52 +/- 0.05, 0.37 +/- 0.05, 0.80 +/- 0.04, and 0.30 +/- 0.021h(-1) for PU, VC, SP, and PE, respectively. Specific substrate utilization rates of 1.08, 0.11, and 0.86 mg COD/mg VS day were obtained for PU, VC, and PE, respectively. Although SP yielded the highest kinetic coefficient, PU was considered the best support, since SP presented loss of chemical constituents during the reactor`s operational phase. In addition, findings oil the microbial community were associated with the reactor`s performance data. Although PE did not show a satisfactory performance, an interesting microbial diversity was found oil its surface. Based oil the morphology and denaturing gradient gel electrophoresis (DGGE) results, PE showed the best capacity for promoting the attachment of methanogenic organisms, and is therefore a material that merits further analysis. PU was considered the Most suitable material showing the best performance in terms of efficiency of solids and COD removal. (C) 2007 Elsevier Ltd. All rights reserved.

Degradation of detergent (linear alkylbenzene sulfonate) in an anaerobic stirred sequencing-batch reactor containing granular biomass

This study aimed to determine the efficiency of an anaerobic stirred sequencing-batch reactor containing granular biomass for the degradation of linear alkylbenzene sulfonate (LAS), a surfactant present in household detergent. The bioreactor was monitored for LAS concentrations in the influent, effluent and sludge, pH, chemical oxygen demand, bicarbonate alkalinity, total solids, and volatile solids. The degradation of LAS was found to be higher in the absence of co-substrates (53%) than in their presence (24-37%). Using the polymerase chain reaction and denaturing gradient gel electrophoresis (PCR/DGGE), we identified populations of microorganisms from the Bacteria and Archaea domains. Among the bacteria, we identified uncultivated populations of Arcanobacterium spp. (94%) and Opitutus spp. (96%). Among the Archaea, we identified Methanospirillum spp. (90%), Methanosaeta spp. (98%), and Methanobacterium spp. (96%). The presence of methanogenic microorganisms shows that LAS did not inhibit anaerobic digestion. Sampling at the last stage of reactor operation recovered 61 clones belonging to the domain bacteria. These represented a variety of phyla: 34% shared significant homology with Bacteroidetes, 18% with Proteobacteria, 11% with Verrucomicrobia, 8% with Fibrobacteres, 2% with Acidobacteria, 3% with Chlorobi and Firmicutes, and 1% with Acidobacteres and Chloroflexi. A small fraction of the clones (13%) were not related to any phylum. Published by Elsevier Ltd.