RITAN: rapid integration of term annotation and network resources

Ehmt2 inactivation in pancreatic epithelial cells shapes the transcriptional landscape and inflammation response of the whole pancreas

Introduction: The Euchromatic Histone Methyl Transferase Protein 2 (EHMT2), also known as G9a, deposits transcriptionally repressive chromatin marks that play pivotal roles in the maturation and homeostasis of multiple organs. Recently, we have shown that Ehmt2 inactivation in the mouse pancreas alters growth and immune gene expression networks, antagonizing Kras-mediated pancreatic cancer initiation and promotion. Here, we elucidate the essential role of Ehmt2 in maintaining a transcriptional landscape that protects organs from inflammation. Methods: Comparative RNA-seq studies between normal postnatal and young adult pancreatic tissue from Ehmt2 conditional knockout animals (Ehmt2 fl/fl ) targeted to the exocrine pancreatic epithelial cells (Pdx1-Cre and P48 Cre/+ ), reveal alterations in gene expression networks in the whole organ related to injury-inflammation-repair, suggesting an increased predisposition to damage. Thus, we induced an inflammation repair response in the Ehmt2 fl/fl pancreas and used a data science-based approach to integrate RNA-seq-derived pathways and networks, deconvolution digital cytology, and spatial transcriptomics. We also analyzed the tissue response to damage at the morphological, biochemical, and molecular pathology levels. Results and discussion: The Ehmt2 fl/fl pancreas displays an enhanced injury-inflammation-repair response, offering insights into fundamental molecular and cellular mechanisms involved in this process. More importantly, these data show that conditional Ehmt2 inactivation in exocrine cells reprograms the local environment to recruit mesenchymal and immunological cells needed to mount an increased inflammatory response. Mechanistically, this response is an enhanced injury-inflammation-repair reaction with a small contribution of specific Ehmt2-regulated transcripts. Thus, this new knowledge extends the mechanisms underlying the role of the Ehmt2-mediated pathway in suppressing pancreatic cancer initiation and modulating inflammatory pancreatic diseases.

Transcriptional Profiling Underscores the Role of Preprocurement Allograft Metabolism and Innate Immune Status on Outcomes in Human Liver Transplantation

Objective

The adverse effects of ischemia-reperfusion injury (IRI) remain a principal barrier to a successful outcome after lifesaving orthotopic liver transplantation (OLT). Gene expression during different phases of IRI is dynamic and modified by individual exposures, making it attractive for identifying potential therapeutic targets for improving the number of suitable organs for transplantation and patient outcomes. However, data remain limited on the functional landscape of gene expression during liver graft IRI, spanning procurement to reperfusion and recovery. Therefore, we sought to characterize transcriptomic profiles of IRI during multiple phases in human OLT.

Methods

We conducted clinical data analyses, histologic evaluation, and RNA sequencing of 17 consecutive human primary OLT. We performed liver allograft biopsies at 4 time points: baseline (B, before donor cross-clamp), at the end of cold ischemia (CI), during early reperfusion (ER, after revascularization), and during late reperfusion (LR). Data were generated and then recipients grouped by post-OLT outcomes categories: immediate allograft function (IAF; n = 11) versus early allograft dysfunction (EAD; n = 6) groups.

Results

We observed that CI (vs B) modified a transcriptomic landscape enriched for a metabolic and immune process. Expression levels of hallmark inflammatory response genes were higher transitioning from CI to ER and decreased from ER to LR. IAF group predominantly showed higher bile and fatty acid metabolism activity during LR compared with EAD group, while EAD group maintained more immunomodulatory activities. Throughout all time points, EAD specimens exhibited decreased metabolic activity in both bile and fatty acid pathways.

Conclusions

We report transcriptomic profiles of human liver allograft IRI from prepreservation in the donor to posttransplantation in the recipient. Immunomodulatory and metabolic landscapes across ER and LR phases were different between IAF and EAD allografts. Our study also highlights marker genes for these biological processes that we plan to explore as novel therapeutic targets or surrogate markers for severe allograft injury in clinical OLT.

Changes in Daily Apparent Diffusion Coefficient on Fully Quantitative Magnetic Resonance Imaging Correlate With Established Genomic Pathways of Radiation Sensitivity and Reveal Novel Biologic Associations

PURPOSE

Changes in quantitative magnetic resonance imaging (qMRI) are frequently observed during chemotherapy or radiation therapy (RT). It is hypothesized that qMRI features are reflective of underlying tissue responses. It's unknown what underlying genomic characteristics underly qMRI changes. We hypothesized that qMRI changes may correlate with DNA damage response (DDR) capacity within human tumors. Therefore, we designed the current study to correlate qMRI changes from daily RT treatment with underlying tumor transcriptomic profiles.

METHODS AND MATERIALS

Study participants were prospectively enrolled (National Clinical Trial 03500081). RNA expression levels for 757 genes from pretreatment biopsies were obtained using a custom panel that included signatures of radiation sensitivity and DDR. Daily qMRI data were obtained from a 1.5 Tesla MR linear accelerator. Using these images, d-slow, d-star, perfusion, and apparent diffusion coefficient-mean values in tumors were plotted per-fraction, over time, and associated with genomic pathways.

RESULTS

A total of 1022 qMRIs were obtained from 39 patients and both genomic data and qMRI data from 27 total patients. For 20 of those patients, we also generated normal tissue transcriptomic data. Radio sensitivity index values most closely associated with tissue of origin. Multiple genomic pathways including DNA repair, peroxisome, late estrogen receptor responses, KRAS signaling, and UV response were significantly associated with qMRI feature changes (P < .001).

CONCLUSIONS

Genomic pathway associations across metabolic, RT sensitivity, and DDR pathways indicate common tumor biology that may correlate with qMRI changes during a course of treatment. Such data provide hypothesis-generating novel mechanistic insight into the biologic meaning of qMRI changes during treatment and enable optimal selection of imaging biomarkers for biologically MR-guided RT.

EHMT2 Inactivation in Pancreatic Epithelial Cells Shapes the Transcriptional Landscape and Inflammation Response of the Whole Pancreas

The Euchromatic Histone Methyl Transferase Protein 2 (EHMT2), also known as G9a, deposits transcriptionally repressive chromatin marks that play pivotal roles in the maturation and homeostasis of multiple organs. Recently, we have shown that EHMT2 inactivation alters growth and immune gene expression networks, antagonizing KRAS-mediated pancreatic cancer initiation and promotion. Here, we elucidate the essential role of EHMT2 in maintaining a transcriptional landscape that protects organs from inflammation. Comparative RNA-seq studies between normal postnatal and young adult pancreatic tissue from EHMT2 conditional knockout animals ( EHMT2 fl/fl ) targeted to the exocrine pancreatic epithelial cells ( Pdx1-Cre and P48 Cre/+ ), reveal alterations in gene expression networks in the whole organ related to injury-inflammation-repair, suggesting an increased predisposition to damage. Thus, we induced an inflammation repair response in the EHMT2 fl/fl pancreas and used a data science-based approach to integrate RNA-seq-derived pathways and networks, deconvolution digital cytology, and spatial transcriptomics. We also analyzed the tissue response to damage at the morphological, biochemical, and molecular pathology levels. The EHMT2 fl/fl pancreas displays an enhanced injury-inflammation-repair response, offering insights into fundamental molecular and cellular mechanisms involved in this process. More importantly, these data show that conditional EHMT2 inactivation in exocrine cells reprograms the local environment to recruit mesenchymal and immunological cells needed to mount an increased inflammatory response. Mechanistically, this response is an enhanced injury-inflammation-repair reaction with a small contribution of specific EHMT2-regulated transcripts. Thus, this new knowledge extends the mechanisms underlying the role of the EHMT2-mediated pathway in suppressing pancreatic cancer initiation and modulating inflammatory pancreatic diseases.

Hypomethylation and Overexpression of Th17-Associated Genes is a Hallmark of Intestinal CD4+ Lymphocytes in Crohn's Disease

BACKGROUND

The development of Crohn's disease [CD] involves immune cell signalling pathways regulated by epigenetic modifications. Aberrant DNA methylation has been identified in peripheral blood and bulk intestinal tissue from CD patients. However, the DNA methylome of disease-associated intestinal CD4+ lymphocytes has not been evaluated.

MATERIALS AND METHODS

Genome-wide DNA methylation sequencing was performed from terminal ileum CD4+ cells from 21 CD patients and 12 age- and sex-matched controls. Data were analysed for differentially methylated CpGs [DMCs] and methylated regions [DMRs]. Integration was performed with RNA-sequencing data to evaluate the functional impact of DNA methylation changes on gene expression. DMRs were overlapped with regions of differentially open chromatin [by ATAC-seq] and CCCTC-binding factor [CTCF] binding sites [by ChIP-seq] between peripherally derived Th17 and Treg cells.

RESULTS

CD4+ cells in CD patients had significantly increased DNA methylation compared to those from the controls. A total of 119 051 DMCs and 8113 DMRs were detected. While hypermethylated genes were mostly related to cell metabolism and homeostasis, hypomethylated genes were significantly enriched within the Th17 signalling pathway. The differentially enriched ATAC regions in Th17 cells [compared to Tregs] were hypomethylated in CD patients, suggesting heightened Th17 activity. There was significant overlap between hypomethylated DNA regions and CTCF-associated binding sites.

CONCLUSIONS

The methylome of CD patients shows an overall dominant hypermethylation yet hypomethylation is more concentrated in proinflammatory pathways, including Th17 differentiation. Hypomethylation of Th17-related genes associated with areas of open chromatin and CTCF binding sites constitutes a hallmark of CD-associated intestinal CD4+ cells.

Comparative Analysis of Proteomes and Phosphoproteomes in Patients with Prostate Cancer Using Different Surgical Conditions

Purpose

To establish the standard of procedure in preparing benign and cancerous prostate tissues and evaluate the quality of proteomics and phosphoproteomics during transurethral resection of the prostate (TUR-P) with different surgical conditions.

Materials and Methods

TUR-P tissue samples from three patients, two diagnosed with prostate cancer and one with benign prostatic hyperplasia, were each analyzed under three different conditions, based on differences in energy values, tissue locations, and surgical techniques. Global- and phosphorylated proteomic profiles of prostate tissues were analyzed by liquid chromatography-tandem mass spectrometry.

Results

A total of 6,019 global proteins and 4,280 phosphorylated peptides were identified in the nine tissues. The quantitative distributions of proteins and phosphorylation in tissues from the same patient were not affected by changes in the surgical conditions, but indirect relative comparisons differed among patients. Phosphorylation levels, especially of proteins involved in the androgen receptor pathway, important in prostate cancer, were preserved in each patient.

Conclusions

Proteomic profiles of prostate tissue collected by TUR-P were not significantly affected by energy levels, tissue location, or surgical technique. In addition, since protein denaturation of samples through TUR-P is rarely confirmed in this study, we think that it will be an important guide for tissue samples in castration resistant prostate cancer patients, where it is difficult to obtain tissue. This result is the first report about proteomic and phosphoproteomic results with TUR-P samples in prostate cancer and will be theoretical basis in protein analysis research with prostate cancer tissues.

Comprehensive methylome sequencing reveals prognostic epigenetic biomarkers for prostate cancer mortality

BACKGROUND

Prostate cancer is a clinically heterogeneous disease with a subset of patients rapidly progressing to lethal-metastatic prostate cancer. Current clinicopathological measures are imperfect predictors of disease progression. Epigenetic changes are amongst the earliest molecular changes in tumourigenesis. To find new prognostic biomarkers to enable earlier intervention and improved outcomes, we performed methylome sequencing of DNA from patients with localised prostate cancer and long-term clinical follow-up.

METHODS

We used whole-genome bisulphite sequencing (WGBS) to comprehensively map and compare DNA methylation of radical prostatectomy tissue between patients with lethal disease (n = 7) and non-lethal (n = 8) disease (median follow-up 19.5 years). Validation of differentially methylated regions (DMRs) was performed in an independent cohort (n = 185, median follow-up 15 years) using targeted multiplex bisulphite sequencing of candidate regions. Survival was assessed via univariable and multivariable analyses including clinicopathological measures (log-rank and Cox regression models).

RESULTS

WGBS data analysis identified cancer-specific methylation patterns including CpG island hypermethylation, and hypomethylation of repetitive elements, with increasing disease risk. We identified 1420 DMRs associated with prostate cancer-specific mortality (PCSM), which showed enrichment for gene sets downregulated in prostate cancer and de novo methylated in cancer. Through comparison with public prostate cancer datasets, we refined the DMRs to develop an 18-gene prognostic panel. Applying this panel to an independent cohort, we found significant associations between PCSM and hypermethylation at EPHB3, PARP6, TBX1, MARCH6 and a regulatory element within CACNA2D4. Strikingly in a multivariable model, inclusion of CACNA2D4 methylation was a better predictor of PCSM versus grade alone (Harrell's C-index: 0.779 vs. 0.684).

CONCLUSIONS

Our study provides detailed methylome maps of non-lethal and lethal prostate cancer and identifies novel genic regions that distinguish these patient groups. Inclusion of our DNA methylation biomarkers with existing clinicopathological measures improves prognostic models of prostate cancer mortality, and holds promise for clinical application.

KrasG12D induces changes in chromatin territories that differentially impact early nuclear reprogramming in pancreatic cells

BACKGROUND

Pancreatic ductal adenocarcinoma initiation is most frequently caused by Kras mutations.

RESULTS

Here, we apply biological, biochemical, and network biology methods to validate GEMM-derived cell models using inducible KrasG12D expression. We describe the time-dependent, chromatin remodeling program that impacts function during early oncogenic signaling. We find that the KrasG12D-induced transcriptional response is dominated by downregulated expression concordant with layers of epigenetic events. More open chromatin characterizes the ATAC-seq profile associated with a smaller group of upregulated genes and epigenetic marks. RRBS demonstrates that promoter hypermethylation does not account for the silencing of the extensive gene promoter network. Moreover, ChIP-Seq reveals that heterochromatin reorganization plays little role in this early transcriptional program. Notably, both gene activation and silencing primarily depend on the marking of genes with a combination of H3K27ac, H3K4me3, and H3K36me3. Indeed, integrated modeling of all these datasets shows that KrasG12D regulates its transcriptional program primarily through unique super-enhancers and enhancers, and marking specific gene promoters and bodies. We also report chromatin remodeling across genomic areas that, although not contributing directly to cis-gene transcription, are likely important for KrasG12D functions.

CONCLUSIONS

In summary, we report a comprehensive, time-dependent, and coordinated early epigenomic program for KrasG12D in pancreatic cells, which is mechanistically relevant to understanding chromatin remodeling events underlying transcriptional outcomes needed for the function of this oncogene.

Inactivation of the Euchromatic Histone-Lysine N-Methyltransferase 2 Pathway in Pancreatic Epithelial Cells Antagonizes Cancer Initiation and Pancreatitis-Associated Promotion by Altering Growth and Immune Gene Expression Networks

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive, painful disease with a 5-year survival rate of only 9%. Recent evidence indicates that distinct epigenomic landscapes underlie PDAC progression, identifying the H3K9me pathway as important to its pathobiology. Here, we delineate the role of Euchromatic Histone-lysine N-Methyltransferase 2 (EHMT2), the enzyme that generates H3K9me, as a downstream effector of oncogenic KRAS during PDAC initiation and pancreatitis-associated promotion. EHMT2 inactivation in pancreatic cells reduces H3K9me2 and antagonizes Kras G12D -mediated acinar-to-ductal metaplasia (ADM) and Pancreatic Intraepithelial Neoplasia (PanIN) formation in both the Pdx1-Cre and P48 Cre/+ Kras G12D mouse models. Ex vivo acinar explants also show impaired EGFR-KRAS-MAPK pathway-mediated ADM upon EHMT2 deletion. Notably, Kras G12D increases EHMT2 protein levels and EHMT2-EHMT1-WIZ complex formation. Transcriptome analysis reveals that EHMT2 inactivation upregulates a cell cycle inhibitory gene expression network that converges on the Cdkn1a/p21-Chek2 pathway. Congruently, pancreas tissue from Kras G12D animals with EHMT2 inactivation have increased P21 protein levels and enhanced senescence. Furthermore, loss of EHMT2 reduces inflammatory cell infiltration typically induced during Kras G12D -mediated initiation. The inhibitory effect on Kras G12D -induced growth is maintained in the pancreatitis-accelerated model, while simultaneously modifying immunoregulatory gene networks that also contribute to carcinogenesis. This study outlines the existence of a novel KRAS-EHMT2 pathway that is critical for mediating the growth-promoting and immunoregulatory effects of this oncogene in vivo, extending human observations to support a pathophysiological role for the H3K9me pathway in PDAC.

A Proteomics-Based Analysis of Blood Biomarkers for the Diagnosis of COPD Acute Exacerbation

Purpose

The identification of blood biomarkers to diagnose acute exacerbation of chronic obstructive pulmonary disease (AECOPD) will have clinical utility. Here, we used a proteomics-based approach to identify biomarkers capable of identifying AECOPD.

Patients and Methods

This prospective, single-center pilot study enrolled 12 patients who came to Asan Medical Center (South Korea) via the outpatient clinic or emergency department with symptoms of AECOPD and were follow-up in the outpatient clinic during convalescence between 2015 and 2017. Paired blood samples collected from each patient during the treatment naïve AECOPD and convalescence stages were analyzed. A sequential window acquisition of all theoretical fragmentation spectra-mass spectrometry (SWATH-MS)-based proteome analysis was performed and a subset of the data were verified by ELISA.

Results

The SWATH-MS analysis identified 226 plasma proteins across all samples examined. The median coefficient of variation for triplicate technical replicates of each sample was 1.13 ± 1.38%, indicating high precision of the technique. Fold-change and paired t-test analyses revealed that 14 proteins were present at higher levels in the AECOPD samples than in the convalescence samples. A gene ontology analysis revealed that these proteins are involved in the acute-phase response. A total of 15 proteins were present at higher levels during the recovery (convalescence) stage than during the acute exacerbation phase, and gene ontology analysis revealed that these proteins are related to lipid metabolism and transport. Verification of the SWATH-MS data was performed using ELISAs for three proteins that were up-regulated in AECOPD, namely, LBP, ORM2, and SERPINA3. Among them, SERPINA3 (p = 0.005) was up-regulated significantly in AECOPD compared with the convalescence state.

Conclusion

Potential plasma biomarkers of AECOPD were discovered using the SWATH-MS proteomics method, and functional molecular associations were investigated. SERPINA3 could be a promising diagnostic biomarker for the early identification and tracking of AECOPD.

PyGNA: a unified framework for geneset network analysis

BACKGROUND

Gene and protein interaction experiments provide unique opportunities to study the molecular wiring of a cell. Integrating high-throughput functional genomics data with this information can help identifying networks associated with complex diseases and phenotypes.

RESULTS

Here we introduce an integrated statistical framework to test network properties of single and multiple genesets under different interaction models. We implemented this framework as an open-source software, called Python Geneset Network Analysis (PyGNA). Our software is designed for easy integration into existing analysis pipelines and to generate high quality figures and reports. We also developed PyGNA to take advantage of multi-core systems to generate calibrated null distributions on large datasets. We then present the results of extensive benchmarking of the tests implemented in PyGNA and a use case inspired by RNA sequencing data analysis, showing how PyGNA can be easily integrated to study biological networks. PyGNA is available at http://github.com/stracquadaniolab/pygna and can be easily installed using the PyPi or Anaconda package managers, and Docker.

CONCLUSIONS

We present a tool for network-aware geneset analysis. PyGNA can either be readily used and easily integrated into existing high-performance data analysis pipelines or as a Python package to implement new tests and analyses. With the increasing availability of population-scale omic data, PyGNA provides a viable approach for large scale geneset network analysis.

Inorganic arsenic promotes luminal to basal transition and metastasis of breast cancer

Inorganic arsenic (iAs/As₂ O₃ ^2- ) is an environmental toxicant found in watersheds around the world including in densely populated areas. iAs is a class I carcinogen known to target the skin, lungs, bladder, and digestive organs, but its role as a primary breast carcinogen remains controversial. Here, we examined a different possibility: that exposure to iAs promotes the transition of well-differentiated epithelial breast cancer cells characterized by estrogen and progesterone receptor expression (ER+/PR+), to more basal phenotypes characterized by active proliferation, and propensity to metastasis in vivo. Our results indicate two clear phenotypic responses to low-level iAs that depend on the duration of the exposure. Short-term pulses of iAs activate ER signaling, consistent with its reported pseudo-estrogen activity, but longer-term, chronic treatments for over 6 months suppresses both ER and PR expression and signaling. In fact, washout of these chronically exposed cells for up to 1 month failed to fully reverse the transcriptional and phenotypic effects of prolonged treatments, indicating durable changes in cellular physiologic identity. RNA-seq studies found that chronic iAs drives the transition toward more basal phenotypes characterized by impaired hormone receptor signaling despite the conservation of estrogen receptor expression. Because treatments for breast cancer patients are largely designed based on the detection of hormone receptor expression, our results suggest greater scrutiny of ER+ cancers in patients exposed to iAs, because these tumors may spawn more aggressive phenotypes than unexposed ER+ tumors, in particular, basal subtypes that tend to develop therapy resistance and metastasis.

Genome-wide resolution peripheral blood methylome profiling reveals signatures for cholestatic liver disease

Aim:

To profile DNA methylation changes of primary biliary cholangitis (PBC) and primary sclerosing cholangitis (PSC).

Materials & methods:

Patients with: PBC, PSC with inflammatory bowel disease (IBD), PSC without IBD, and age-, sex-matched controls were profiled for methylomes of peripheral blood by reduced representation bisulfite sequencing. Differentially methylated CpG (DMC) and differentially methylated region (DMR) were detected and compared.

Results:

We identified consistently altered DMCs and DMRs across diseases with involvement in key pathways. Many similarities noted between two subtypes of PSC, interestingly few existed between PBC and PSC. DMRs were highly enriched with transcription factor binding. Top DMC changes were validated in liver tissue of an independent cohort.

Conclusion:

Methylome profiling provides insights to PBC and PSC.

Transcriptome Profiling of Acquired Gefitinib Resistant Lung Cancer Cells Reveals Dramatically Changed Transcription Programs and New Treatment Targets

Background: Targeted therapy for lung cancer with epidermal growth factor receptor (EGFR) mutations with tyrosine kinase inhibitors (TKIs) represents one of the major breakthroughs in lung cancer management. However, gradually developed resistance to these drugs prevents sustained clinical benefits and calls for resistant mechanism research and identification of new therapeutic targets. Acquired T790M mutation accounts for the majority of resistance cases, yet transcriptome changes in these cells are less characterized, and it is not known if new treatment targets exist by available drugs.

Methods: Transcriptome profiling was performed for lung cancer cell line PC9 and its resistant line PC9GR after long-term exposure to gefitinib through RNA sequencing. Differentially expressed genes and changed pathways were identified along with existing drugs targeting these upregulated genes. Using 144 lung cancer cell lines with both gene expression and drug response data from the cancer cell line encyclopedia (CCLE) and Cancer Therapeutics Response Portal (CTRP), we screened 549 drugs whose response was correlated with these upregulated genes in PC9GR cells, and top drugs were evaluated for their response in both PC9 and PC9GR cells.

Results: In addition to the acquired T790M mutation, the resistant PC9GR cells had very different transcription programs from the sensitive PC9 cells. Multiple pathways were changed with the top ones including TNFA signaling, androgen/estrogen response, P53 pathway, MTORC1 signaling, hypoxia, and epithelial mesenchymal transition. Thirty-two upregulated genes had available drugs that can potentially be effective in treating the resistant cells. From the response profiles of CCLE, we found 17 drugs whose responses were associated with at least four of these upregulated genes. Among the four drugs evaluated (dasatinib, KPT-185, trametinib, and pluripotin), all except trametinib demonstrated strong inhibitory effects on the resistant PC9GR cells, among which KPT185 was the most potent. KPT-185 suppressed growth, caused apoptosis, and inhibited migration of the PC9GR cells at similar (or better) rates as the sensitive PC9 cells in a dose-dependent manner.

Conclusions: Acquired TKI-resistant lung cancer cells (PC9GR) have dramatically changed transcription and pathway regulation, which expose new treatment targets. Existing drugs may be repurposed to treat those patients with developed resistance to TKIs.

RNA Sequencing Reveals Cancer-Associated Changes in Laryngeal Cells Exposed to Non-Acid Pepsin

OBJECTIVE

Laryngopharyngeal reflux (LPR) is a common affliction that contributes to laryngeal inflammation, symptoms that impact quality of life, and life-threatening illnesses such as cancer. Effective treatment strategies for LPR are lacking. Pepsin is a proinflammatory and carcinogenic element of refluxate. Investigation of molecular pathways involved in pepsin-mediated damage may lead to identification of novel biomarkers and therapeutic targets for LPR. In this study, RNA sequencing was used to examine changes in human laryngeal epithelial cells following brief pepsin insult. Cells were immortalized to generate a model to aid future study of laryngeal injury and therapeutics.

STUDY DESIGN

In vitro translational.

METHODS

Laryngeal epithelial cells were cultured from a patient without signs or symptoms of LPR or laryngeal cancer. Cells were treated with 0.1 mg/ml pepsin for 1 hour or normal growth media (control) prior to RNA sequencing. Cells were immortalized via HPV E6/7 and characterized by microscopy, immunohistochemistry, G-banding, and soft agar assay.

RESULTS

Three hundred ninety-seven genes exhibited differences in expression with pepsin treatment (P < .05). Pathway analysis revealed association with cancer and related signaling processes including dysregulation of cancer-associated molecules, Metastasis-Associated Lung Adenocarcinoma Transcript 1 and KRT82, and the long-noncoding RNA, lipoprotein receptor-related protein 1 (LRP1)-AS, which regulates the putative pepsin receptor LRP1.

CONCLUSIONS

A single, brief exposure to pepsin activated cancer-associated signaling pathways in laryngeal cells in vitro, revealing novel mechanisms by which chronic reflux may contribute to carcinogenesis. The cell line developed herein represents a novel tool in which to investigate pepsin-dysregulated pathways identified by RNA sequencing and disparities of tumor proneness of laryngeal subsites.

LEVEL OF EVIDENCE

N/A Laryngoscope, 131:121-129, 2021.