Pseudoalignment tools as an efficient alternative to detect repeated transposable elements in scRNAseq data.

MOTIVATION: Transposable elements (TE) have played a major role in configuring the structures of mammalian genomes through evolution. In normal conditions, the expression of these elements is repressed by different epigenetic regulation mechanisms such as DNA methylation, histone modification and regulation by small RNAs. TE re-activation is associated with stemness potential acquisition, regulation of innate immunity and disease, such as cancer. However, the vast majority of current knowledge in the field is based on bulk expression studies, and very little is known on cell-type- or state-specific expression of TE-derived transcripts. Therefore, cost-efficient single-cell-resolution TE expression analytical approaches are needed. RESULTS: We have implemented an analytical approach based on pseudoalignment to consensus sequences to incorporate TE expression information to scRNAseq data. AVAILABILITY AND IMPLEMENTATION: All the data and code implemented are available as Supplementary data and in: https://github.com/jmzvillarreal/kallisto_TE_scRNAseq. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Secondary resistance to anti-EGFR therapy by transcriptional reprogramming in patient-derived colorectal cancer models.

BACKGROUND: The development of secondary resistance (SR) in metastatic colorectal cancer (mCRC) treated with anti-epidermal growth factor receptor (anti-EGFR) antibodies is not fully understood at the molecular level. Here we tested in vivo selection of anti-EGFR SR tumors in CRC patient-derived xenograft (PDX) models as a strategy for a molecular dissection of SR mechanisms. METHODS: We analyzed 21 KRAS, NRAS, BRAF, and PI3K wildtype CRC patient-derived xenograft (PDX) models for their anti-EGFR sensitivity. Furthermore, 31 anti-EGFR SR tumors were generated via chronic in vivo treatment with cetuximab. A multi-omics approach was employed to address molecular primary and secondary resistance mechanisms. Gene set enrichment analyses were used to uncover SR pathways. Targeted therapy of SR PDX models was applied to validate selected SR pathways. RESULTS: In vivo anti-EGFR SR could be established with high efficiency. Chronic anti-EGFR treatment of CRC PDX tumors induced parallel evolution of multiple resistant lesions with independent molecular SR mechanisms. Mutations in driver genes explained SR development in a subgroup of CRC PDX models, only. Transcriptional reprogramming inducing anti-EGFR SR was discovered as a common mechanism in CRC PDX models frequently leading to RAS signaling pathway activation. We identified cAMP and STAT3 signaling activation, as well as paracrine and autocrine signaling via growth factors as novel anti-EGFR secondary resistance mechanisms. Secondary resistant xenograft tumors could successfully be treated by addressing identified transcriptional changes by tailored targeted therapies. CONCLUSIONS: Our study demonstrates that SR PDX tumors provide a unique platform to study molecular SR mechanisms and allow testing of multiple treatments for efficient targeting of SR mechanisms, not possible in the patient. Importantly, it suggests that the development of anti-EGFR tolerant cells via transcriptional reprogramming as a cause of anti-EGFR SR in CRC is likely more prevalent than previously anticipated. It emphasizes the need for analyses of SR tumor tissues at a multi-omics level for a comprehensive molecular understanding of anti-EGFR SR in CRC.

Pancreatic duct ligation reduces premalignant pancreatic lesions in a Kras model of pancreatic adenocarcinoma in mice.

Pancreatic duct ligation (PDL) in the murine model has been described as an exocrine pancreatic atrophy-inducing procedure. However, its influence has scarcely been described on premalignant lesions. This study describes the histological changes of premalignant lesions and the gene expression in a well-defined model of pancreatic ductal adenocarcinoma by PDL. Selective ligation of the splenic lobe of the pancreas was performed in Ptf1a-Cre(+/ki); K-ras LSLG12Vgeo(+/ki) mice (PDL-Kras mice). Three experimental groups were evaluated: PDL group, controls and shams. The presence and number of premalignant lesions (PanIN 1-3 and Atypical Flat Lesions-AFL) in proximal (PP) and distal (DP) pancreas were studied for each group over time. Microarray analysis was performed to find differentially expressed genes (DEG) between PP and PD. Clinical human specimens after pancreaticoduodenectomy with ductal occlusion were also evaluated. PDL-Kras mice showed an intense pattern of atrophy in DP which was shrunk to a minimal portion of tissue. Mice in control and sham groups had a 7 and 10-time increase respectively of risk of high-grade PanIN 2 and 3 and AFL in their DP than PDL-Kras mice. Furthermore, PDL-Kras mice had significantly less PanIN 1 and 2 and AFL lesions in DP compared to PP. We identified 38 DEGs comparing PP and PD. Among them, several mapped to protein secretion and digestion while others such as Nupr1 have been previously associated with PanIN and PDAC. PDL in Ptf1a-Cre(+/ki); K-ras LSLG12Vgeo(+/ki) mice induces a decrease in the presence of premalignant lesions in the ligated DP. This could be a potential line of research of interest in some cancerous risk patients.