Division of Computational Biology, Mayo Clinic, Rochester, Minnesota.
Citations 42
Discovery and Validation of Molecular Biomarkers for Differentiation of Nondysplastic Barrett’s Esophagus from High-grade Dysplasia and Esophageal Adenocarcinoma
Abstract Aberrant DNA methylation and copy-number alterations (CNA) drive Barrett’s esophagus progression to esophageal adenocarcinoma; however, their combined utility for early detection is unclear. We aimed to identify and validate methylated DNA markers (MDM) and CNAs to distinguish esophageal adenocarcinoma/high-grade dysplasia (HGD) from nondysplastic Barrett’s esophagus (NDBE). In this multiphase, multicenter study, we discovered and validated MDMs and quantified CNAs utilizing whole-genome methylation sequencing of esophageal brushings. DNA biomarkers identified from discovery were further validated in independent patients with paired esophageal brushing and swallowed capsule sponge samples. MDMs were filtered against a reduced representation bisulfite sequencing dataset obtained from independent tissue samples to advance only concordant candidates. CNA burden was quantified using ichorCNA-derived aneuploidy scores (AS). Two hundred MDMs discovered in HGD (N = 18) and esophageal adenocarcinoma (N = 18) versus NDBE brushing samples (N = 18) were tested in independent samples (N = 146). A 52-MDM panel achieved a cross-validated AUC of 0.88 [95% confidence interval (CI), 0.82–0.95]; the addition of AS improved discrimination of HGD/esophageal adenocarcinoma from NDBE to 0.91 (95% CI, 0.86–0.97) AUC. At 80% specificity, the combined model detected 93% of esophageal adenocarcinoma and 88% of HGD cases. In paired capsule sponge samples, a 58-MDM panel achieved a cross-validated AUC of 0.77 (95% CI, 0.66–0.88); a combined 58-MDM and AS model achieved AUC 0.80 (95% CI, 0.7–0.9). MDMs and AS discerned HGD/esophageal adenocarcinoma from normal esophagus/NDBE in endoscopic brushing and capsule sponge samples. This approach may improve Barrett’s esophagus surveillance. Prevention Relevance: This study demonstrates that combining epigenetic and genomic biomarkers across minimally invasive sampling methods can accurately distinguish HGD/esophageal adenocarcinoma from nondysplastic Barrett’s esophagus, offering promising, less invasive strategies to improve BE surveillance and enable endoscopic therapy for esophageal adenocarcinoma prevention and treatment.
Abnormal ERV Expression and Its Clinical Relevance in Colon Cancer
Background/Objectives: Human endogenous retroviruses (ERVs) are genomic sequences integrated into the human genome from ancestral exogenous retroviruses and are epigenetically silenced under normal conditions. Growing evidence has shown that they can be reactivated in human diseases such as cancers and autoimmune diseases. However, their clinical implications in colon cancer are yet to be explored. Methods: RNA-seq data were downloaded from RNA Atlas and TCGA for cell lines and tissue samples, respectively. After alignment, ERV expression was quantified against comprehensively compiled ERVs (3220). ERV expression profiles were compared between sequencing protocols, cancer and normal cells, and matched tumor and normal tissue pairs. Unsupervised clustering was used to identify ERV-defined tumor subtypes and their associations with clinical and other molecular features. ERV association with disease-specific survival (DSS) was performed using the Cox regression model. Results: PolyA and total RNA protocols were comparable in ERV expression detection. Cancer cells had significantly increased ERV expression and reactivation. Upregulated ERVs were significantly enriched in viral protein interactions with cytokine and cytokine receptors. ERV expression-defined tumor classes were significantly associated with tumor mutation burden and immuno-phenotypes such as antigen processing and presenting machinery and tumor immune infiltration score. Survival analysis identified 152 ERVs to be independently associated with DSS. Conclusions: ERV abnormal expression is common in colon cancer. The ERV-defined subtypes are associated with tumor immunity, and some ERVs are independently associated with patient outcomes.
Epigenetic disease markers in primary sclerosing cholangitis and primary biliary cholangitis—methylomics of cholestatic liver disease
Background: The epigenome, the set of modifications to DNA and associated molecules that control gene expression, cellular identity, and function, plays a major role in mediating cellular responses to outside factors. Thus, evaluation of the epigenetic state can provide insights into cellular adaptions occurring over the course of disease. Methods: We performed epigenome-wide association studies of primary sclerosing cholangitis (PSC) and primary biliary cholangitis (PBC) using the Illumina MethylationEPIC Bead Chip. Results: We found evidence of increased epigenetic age acceleration and differences in predicted immune cell composition in patients with PSC and PBC. Epigenetic profiles demonstrated differences in predicted protein levels including increased levels of tumor necrosis factor receptor superfamily member 1B in patients with cirrhotic compared to noncirrhotic PSC and PBC. Epigenome-wide association studies of PSC discovered strongly associated 5′-C-phosphate-G-3′ sites in genes including vacuole membrane protein 1 and SOCS3, and epigenome-wide association studies of PBC found strong 5′-C-phosphate-G-3′ associations in genes including NOD-like receptor family CARD domain containing 5, human leukocyte antigen-E, and PSMB8. Analyses identified disease-associated canonical pathways and upstream regulators involved with immune signaling and activation of macrophages and T-cells. A comparison of PSC and PBC data found relatively little overlap at the 5′-C-phosphate-G-3′ and gene levels with slightly more overlap at the level of pathways and upstream regulators. Conclusions: This study provides insights into methylation profiles of patients that support current concepts of disease mechanisms and provide novel data to inspire future research. Studies to corroborate our findings and expand into other -omics layers will be invaluable to further our understanding of these rare diseases with the goal to improve and individualize prognosis and treatment.
Performance Comparisons of Methylation and Structural Variants From Low-Input Whole-Genome Methylation Sequencing
Long non-coding RNAs (lncRNAs) are a large class of gene transcripts that do not code proteins; however, their functions are largely unknown and many new lncRNAs are yet to be discovered. Taking advantage of our previously developed, super-fast, novel lncRNA discovery pipeline, UClncR, and rich resources of GTEx RNA-seq data, we performed systematic novel lincRNA discovery for over 8000 samples across 30 tissue types. We conducted novel detection for each major tissue type first and then consolidated the novel discoveries from all tissue types. These novel lincRNs were profiled and analyzed along with known genes to identify tissue-specific genes in 30 major human tissue types. Thirteen sub-brain regions were also analyzed in a similar manner. Our analysis revealed thousands to tens of thousands of novel lincRNAs for each tissue type. These lincRNAs could define each tissue type’s identity and demonstrated their reliability and tissue-specific expression. Tissue-specific genes were identified for each major tissue type and sub-brain region. The tissue-specific genes clearly defined each respective tissue’s unique function and could be used to expand the interpretation of non-coding SNPs from genome-wide association (GWAS) studies.
Spectrum of EGFR aberrations and potential clinical implications: insights from integrative pan‐cancer analysis
Abstract Background Human epidermal growth factor receptor ( EGFR ) is an oncogenic gene and one of top targets of precision therapy in lung cancer with EGFR mutations. Although there are many reports for some individual cancers, comprehensive profiling of EGFR mutations, overexpression, amplification, DNA methylation, and their clinical associations across many different cancers simultaneously was not available. This study aimed to fill the gap and provide insights to the alteration spectrum of EGFR and its therapeutic and prognostic implications. Methods The Cancer Genome Atlas (TCGA) datasets for 32 cancer types involving 11,314 patients were analyzed for alterations (mutations and amplification/deletion), abnormal expression and DNA methylation in EGFR gene. Mutation frequency, genomic location distribution, functional impact, and clinical targeted therapy implication were compared among different cancer types, and their associations with patient survival were analyzed. Results EGFR alteration frequency, mutation sites across functional domains, amplification, overexpression, and DNA methylation patterns differed greatly among different cancer types. The overall mutation frequency in all cancers combined was relatively low. Targetable mutations, mainly in lung cancer, were primarily found in the Pkinase_Tyr domain. Glioblastoma multiforme had the highest rate of alterations, but it was dominated by gene amplification and most mutations were in the Furin‐like domain where targeted therapy was less effective. Low‐grade glioma often had gene amplification and increased EGFR expression which was associated with poor outcome. Colon and pancreatic adenocarcinoma had very few EGFR mutations; however, high EGFR expression was significantly associated with short patient survival. Squamous cell carcinoma regardless of their sites (the head and neck, lung, or esophagus) exhibited similar characteristics with an alteration frequency of about 5.0%, was dominated by gene amplification, and had increased EGFR expression generally associated with short patient survival. DNA methylation was highly associated with EGFR expression and patient outcomes in some cancers. Conclusions EGFR aberration type, frequency, distribution in functional domains, and expression vary from cancer to cancer. While mutations in the Pkinase_Tyr domain are more important for treatment selection, increased expression from amplification or deregulation affects more tumor types and leads to worse outcome, which calls for new treatment strategies for these EGFR‐ driven tumors.
Androgen Mediation—and Antiandrogens Mitigation—of the Epidermal Growth Factor Receptor (EGFR) Inhibitor–Induced Rash: Results From a Pilot Randomized Trial and Small Translational Case Series
Background: Although 50% to 90% of patients who receive epidermal growth factor receptor (EGFR) inhibitors develop a rash, options for rash prevention or palliation remain limited. This issue is particularly important from a palliative care standpoint because these agents are prescribed only to patients with incurable cancer. Here, we report (1) gene expression profiling of skin biopsies from patients with an EGFR inhibitor–induced rash and (2) a randomized, placebo-controlled feasibility trial with the antiandrogen, spironolactone. Both investigations were undertaken to begin to explore the hypothesis that androgens mediate EGFR inhibitor–induced rash and that antiandrogens palliate it. Methods/Results: First, 4 skin biopsies from patients with EGFR inhibitor–induced rash (3 men and 1 woman) were subject to gene expression microarray profiling. A public data set of normal skin gene expression (Gene Expression Omnibus, GSE22998) served as a reference. Sixty percent of commonly interrogated androgen receptor genes (207 of 308 between the 2 data sets) were differentially expressed ( P < .05) in the rash samples. Second, in a 17-patient double-blinded, placebo-controlled trial with topical spironolactone applied to the face, although the primary feasibility end point was not achieved, patients in the spironolactone arm received more doses of EGFR inhibitor, and anecdotal photographic evidence suggested salutatory effects of spironolactone on rash. Conclusions: Epidermal growth factor receptor inhibitor–induced rash appears to be androgen-mediated; antiandrogen therapy merits further study for rash prevention/palliation.
Molecular profiling reveals immunogenic cues in anaplastic large cell lymphomas with DUSP22 rearrangements
Key PointsDUSP22-rearranged ALCLs belong to a distinct subset of ALCLs lacking activated STAT3. DUSP22-rearranged ALCLs have a unique molecular signature characterized by DNA hypomethylation and an immunogenic phenotype.
Pathways Impacted by Genomic Alterations in Pulmonary Carcinoid Tumors