RESUMEN
Genome-wide enrichment of methylated DNA followed by sequencing (MeDIP-seq) offers a reasonable compromise between experimental costs and genomic coverage. However, the computational analysis of these experiments is complex, and quantification of the enrichment signals in terms of absolute levels of methylation requires specific transformation. In this work, we present QSEA, Quantitative Sequence Enrichment Analysis, a comprehensive workflow for the modelling and subsequent quantification of MeDIP-seq data. As the central part of the workflow we have developed a Bayesian statistical model that transforms the enrichment read counts to absolute levels of methylation and, thus, enhances interpretability and facilitates comparison with other methylation assays. We suggest several calibration strategies for the critical parameters of the model, either using additional data or fairly general assumptions. By comparing the results with bisulfite sequencing (BS) validation data, we show the improvement of QSEA over existing methods. Additionally, we generated a clinically relevant benchmark data set consisting of methylation enrichment experiments (MeDIP-seq), BS-based validation experiments (Methyl-seq) as well as gene expression experiments (RNA-seq) derived from non-small cell lung cancer patients, and show that the workflow retrieves well-known lung tumour methylation markers that are causative for gene expression changes, demonstrating the applicability of QSEA for clinical studies. QSEA is implemented in R and available from the Bioconductor repository 3.4 (www.bioconductor.org/packages/qsea).
Asunto(s)
Metilación de ADN , Genómica/métodos , Análisis de Secuencia de ADN/métodos , Animales , Teorema de Bayes , Regulación de la Expresión Génica , Humanos , Neoplasias Pulmonares/genética , Ratones , Regiones Promotoras Genéticas , Sulfitos , Flujo de TrabajoRESUMEN
BACKGROUND: Non-small cell lung cancer (NSCLC) is the most common cause of cancer-related deaths worldwide and is primarily treated with radiation, surgery, and platinum-based drugs like cisplatin and carboplatin. The major challenge in the treatment of NSCLC patients is intrinsic or acquired resistance to chemotherapy. Molecular markers predicting the outcome of the patients are urgently needed. METHODS: Here, we employed patient-derived xenografts (PDXs) to detect predictive methylation biomarkers for platin-based therapies. We used MeDIP-Seq to generate genome-wide DNA methylation profiles of 22 PDXs, their parental primary NSCLC, and their corresponding normal tissues and complemented the data with gene expression analyses of the same tissues. Candidate biomarkers were validated with quantitative methylation-specific PCRs (qMSP) in an independent cohort. RESULTS: Comprehensive analyses revealed that differential methylation patterns are highly similar, enriched in PDXs and lung tumor-specific when comparing differences in methylation between PDXs versus primary NSCLC. We identified a set of 40 candidate regions with methylation correlated to carboplatin response and corresponding inverse gene expression pattern even before therapy. This analysis led to the identification of a promoter CpG island methylation of LDL receptor-related protein 12 (LRP12) associated with increased resistance to carboplatin. Validation in an independent patient cohort (n = 35) confirmed that LRP12 methylation status is predictive for therapeutic response of NSCLC patients to platin therapy with a sensitivity of 80% and a specificity of 84% (p < 0.01). Similarly, we find a shorter survival time for patients with LRP12 hypermethylation in the TCGA data set for NSCLC (lung adenocarcinoma). CONCLUSIONS: Using an epigenome-wide sequencing approach, we find differential methylation patterns from primary lung cancer and PDX-derived cancers to be very similar, albeit with a lower degree of differential methylation in primary tumors. We identify LRP12 DNA methylation as a powerful predictive marker for carboplatin resistance. These findings outline a platform for the identification of epigenetic therapy resistance biomarkers based on PDX NSCLC models.
Asunto(s)
Biomarcadores de Tumor/genética , Carboplatino/uso terapéutico , Carcinoma de Pulmón de Células no Pequeñas/tratamiento farmacológico , Carcinoma de Pulmón de Células no Pequeñas/genética , Metilación de ADN/genética , Epigenómica , Proteína 1 Relacionada con Receptor de Lipoproteína de Baja Densidad/genética , Ensayos Antitumor por Modelo de Xenoinjerto , Animales , Biomarcadores de Tumor/metabolismo , Carboplatino/farmacología , Supervivencia sin Enfermedad , Resistencia a Antineoplásicos/genética , Genes Supresores de Tumor , Genoma Humano , Humanos , Proteína 1 Relacionada con Receptor de Lipoproteína de Baja Densidad/metabolismo , Neoplasias Pulmonares/genética , Ratones Desnudos , Regiones Promotoras Genéticas , Resultado del TratamientoRESUMEN
Scalable and genetically stable recombinant adeno-associated virus (rAAV) production systems combined with facile adaptability for an extended repertoire of AAV serotypes are required to keep pace with the rapidly increasing clinical demand. For scalable high-titer production of the full range of rAAV serotypes 1-12, we developed OneBac, consisting of stable insect Sf9 cell lines harboring silent copies of AAV1-12 rep and cap genes induced upon infection with a single baculovirus that also carries the rAAV genome. rAAV burst sizes reach up to 5 × 10(5) benzonase-resistant, highly infectious genomic particles per cell, exceeding typical yields of current rAAV production systems. In contrast to recombinant rep/cap baculovirus strains currently employed for large-scale rAAV production, the Sf9rep/cap cell lines are genetically stable, leading to undiminished rAAV burst sizes over serial passages. Thus, OneBac combines full AAV serotype options with the capacity for stable scale-up production, the current bottleneck for the transition of AAV from gene therapy trials to routine clinical treatment.