Exploring combined effects of DNA methylation and copy number on gene expression with a two-stage approach.

DNA methylation and copy number may be associated with each other to some extent, in positive or negative ways. Whether differential methylation and copy number variation have combined effects on gene expression is largely unknown. We use a multivariate linear model to formulate the relationship among the three genomic measurements, gene expression, copy number, and methylation levels. We propose a method that combines a distance covariance measure and the group LASSO to analyze multiple types of genomic data collectively for the purpose of insightfully revealing how gene expression is potentially affected by both copy number variation and differential methylation levels in cellular process. Our approach is of two stages, the first is a variable screening process in which a variable selection method is utilized by employing the concept of joint distance covariance (JdCov) of random vectors, and the second is to implement a penalized regression approach, a group LASSO, on the screened data of much lower dimension. The two-stage approach is tested in extensive simulation studies and shown to be effective. The approach is then applied to the TCGA Melanoma data, which consists of gene expression, methylation and copy number measurements of more than 300 patients and relationship of some genes with methylation and copy number measurements were revealed for the involved subjects.

Subpopulation commensalism promotes Rac1-dependent invasion of single cells via laminin-332.

Phenotypic heterogeneity poses a significant hurdle for cancer treatment but is under-characterized in the context of tumor invasion. Amidst the range of phenotypic heterogeneity across solid tumor types, collectively invading cells and single cells have been extensively characterized as independent modes of invasion, but their intercellular interactions have rarely been explored. Here, we isolated collectively invading cells and single cells from the heterogeneous 4T1 cell line and observed extensive transcriptional and epigenetic diversity across these subpopulations. By integrating these datasets, we identified laminin-332 as a protein complex exclusively secreted by collectively invading cells. Live-cell imaging revealed that laminin-332 derived from collectively invading cells increased the velocity and directionality of single cells. Despite collectively invading and single cells having similar expression of the integrin α6ß4 dimer, single cells demonstrated higher Rac1 activation upon laminin-332 binding to integrin α6ß4. This mechanism suggests a novel commensal relationship between collectively invading and single cells, wherein collectively invading cells promote the invasive potential of single cells through a laminin-332/Rac1 axis.

Associations of differentially expressed genes with psychoneurological symptoms in patients with head and neck cancer: A longitudinal study.

OBJECTIVE: Patients with head and neck cancer (HNC) experience psychoneurological symptoms (PNS, i.e., depression, fatigue, sleep disturbance, pain, and cognitive dysfunction) during intensity-modulated radiotherapy (IMRT) that negatively impact their functional status, quality of life, and overall survival. The underlying mechanisms for PNS are still not fully understood. This study aimed to examine differentially expressed genes and pathways related to PNS for patients undergoing IMRT (i.e., before, end of, 6 months, and 12 months after IMRT). METHODS: Participants included 142 patients with HNC (mean age 58.9 ± 10.3 years, 72.5% male, 83.1% White). Total RNA extracted from blood leukocytes were used for genome-wide gene expression assays. Linear mixed effects model was used to examine the association between PNS and gene expression across time. Gene Ontology (GO) enrichment analysis was employed to identify pathways related to PNS. RESULTS: A total of 1352 genes (162 upregulated, 1190 downregulated) were significantly associated with PNS across time (false discovery rate (FDR) < 0.05). Among these genes, 112 GO terms were identified (FDR < 0.05). The top 20 GO terms among the significant upregulated genes were related to immune and inflammatory responses, while the top 20 GO terms among the significant downregulated genes were associated with telomere maintenance. CONCLUSION: This study is the first to identify genes and pathways linked to immune and inflammatory responses and telomere maintenance that are associated with PNS in patients with HNC receiving IMRT. Inflammation and aging markers may be candidate biomarkers for PNS. Understanding biological markers may produce targets for novel interventions.

Circulating short chain fatty acids and fatigue in patients with head and neck cancer: A longitudinal prospective study.

Fatigue among patients with head and neck cancer (HNC) has been associated with higher inflammation. Short-chain fatty acids (SCFAs) have been shown to have anti-inflammatory and immunoregulatory effects. Therefore, this study aimed to examine the association between SCFAs and fatigue among patients with HNC undergoing treatment with radiotherapy with or without concurrent chemotherapy. Plasma SCFAs and the Multidimensional Fatigue Inventory-20 were collected prior to and one month after the completion of treatment in 59 HNC patients. The genome-wide gene expression profile was obtained from blood leukocytes prior to treatment. Lower butyrate concentrations were significantly associated with higher fatigue (p = 0.013) independent of time of assessment, controlling for covariates. A similar relationship was observed for iso/valerate (p = 0.025). Comparison of gene expression in individuals with the top and bottom 33% of butyrate or iso/valerate concentrations prior to radiotherapy revealed 1,088 and 881 significantly differentially expressed genes, respectively (raw p < 0.05). The top 10 Gene Ontology terms from the enrichment analyses revealed the involvement of pathways related to cytokines and lipid and fatty acid biosynthesis. These findings suggest that SCFAs may regulate inflammatory and immunometabolic responses and, thereby, reduce inflammatory-related symptoms, such as fatigue.

Live-cell invasive phenotyping uncovers the ALK2/BMP6 iron homeostasis pathway as a therapeutic vulnerability in LKB1-mutant lung cancer.

The acquisition of invasive properties is a prerequisite for tumor progression and metastasis. Molecular subtypes of KRAS-driven lung cancer exhibit distinct modes of invasion that likely contribute to unique growth properties and therapeutic susceptibilities. Despite this, pre-clinical discovery strategies designed to exploit invasive phenotypes are lacking. To address this, we designed an experimental system to screen for targetable signaling pathways linked to active early invasion phenotypes in the two most prominent molecular subtypes, TP53 and LKB1, of KRAS-driven lung adenocarcinoma (LUAD). By combining live-cell imaging of human bronchial epithelial cells in a 3D invasion matrix with RNA transcriptome profiling, we identified the LKB1-specific upregulation of bone morphogenetic protein 6 (BMP6). Examination of early-stage lung cancer patients confirmed upregulation of BMP6 in LKB1-mutant lung tumors. At the molecular level, we find that the canonical iron regulatory hormone Hepcidin is induced via BMP6 signaling upon LKB1 loss, where intact LKB1 kinase activity is necessary to maintain signaling homeostasis. Furthermore, pre-clinical studies in a novel Kras/Lkb1-mutant syngeneic mouse model show that potent growth suppression was achieved by inhibiting the ALK2/BMP6 signaling axis with single agents that are currently in clinical trials. We show that alterations in the iron homeostasis pathway are accompanied by simultaneous upregulation of ferroptosis protection proteins. Thus, LKB1 is sufficient to regulate both the 'gas' and 'breaks' to finely tune iron-regulated tumor progression.

Loss of the endocytic tumor suppressor HD-PTP phenocopies LKB1 and promotes RAS-driven oncogenesis.

Oncogenic RAS mutations drive aggressive cancers that are difficult to treat in the clinic, and while direct inhibition of the most common KRAS variant in lung adenocarcinoma (G12C) is undergoing clinical evaluation, a wide spectrum of oncogenic RAS variants together make up a large percentage of untargetable lung and GI cancers. Here we report that loss-of-function alterations (mutations and deep deletions) in the gene that encodes HD-PTP (PTPN23) occur in up to 14% of lung cancers in the ORIEN Avatar lung cancer cohort, associate with adenosquamous histology, and occur alongside an altered spectrum of KRAS alleles. Furthermore, we show that in publicly available early-stage NSCLC studies loss of HD-PTP is mutually exclusive with loss of LKB1, which suggests they restrict a common oncogenic pathway in early lung tumorigenesis. In support of this, knockdown of HD-PTP in RAS-transformed lung cancer cells is sufficient to promote FAK-dependent invasion. Lastly, knockdown of the Drosophila homolog of HD-PTP (dHD-PTP/Myopic) synergizes to promote RAS-dependent neoplastic progression. Our findings highlight a novel tumor suppressor that can restrict RAS-driven lung cancer oncogenesis and identify a targetable pathway for personalized therapeutic approaches for adenosquamous lung cancer.