Knockdown by shRNA identifies SLC44A5 as a potential therapeutic target in hepatocellular carcinoma

Comparative RNA-Sequencing Analysis Reveals High Complexity and Heterogeneity of Transcriptomic and Immune Profiles in Hepatocellular Carcinoma Tumors of Viral (HBV, HCV) and Non-Viral Etiology

Background and Objectives: Hepatocellular carcinoma (HCC), the most common type of primary liver cancer, is the leading cause of cancer-related mortality. It arises and progresses against fibrotic or cirrhotic backgrounds mainly due to infection with hepatitis viruses B (HBV) or C (HCV) or non-viral causes that lead to chronic inflammation and genomic changes. A better understanding of molecular and immune mechanisms in HCC subtypes is needed. Materials and Methods: To identify transcriptional changes in primary HCC tumors with or without hepatitis viral etiology, we analyzed the transcriptomes of 24 patients by next-generation sequencing. Results: We identified common and unique differentially expressed genes for each etiological tumor group and analyzed the expression of SLC, ATP binding cassette, cytochrome 450, cancer testis, and heat shock protein genes. Metascape functional enrichment analysis showed mainly upregulated cell-cycle pathways in HBV and HCV and upregulated cell response to stress in non-viral infection. GeneWalk analysis identified regulator, hub, and moonlighting genes and highlighted CCNB1, ACTN2, BRCA1, IGF1, CDK1, AURKA, AURKB, and TOP2A in the HCV group and HSF1, HSPA1A, HSP90AA1, HSPB1, HSPA5, PTK2, and AURKB in the group without viral infection as hub genes. Immune infiltrate analysis showed that T cell, cytotoxic, and natural killer cell markers were significantly more highly expressed in HCV than in non-viral tumors. Genes associated with monocyte activation had the highest expression levels in HBV, while high expression of genes involved in primary adaptive immune response and complement receptor activity characterized tumors without viral infection. Conclusions: Our comprehensive study underlines the high degree of complexity of immune profiles in the analyzed groups, which adds to the heterogeneous HCC genomic landscape. The biomarkers identified in each HCC group might serve as therapeutic targets.

Gene expression profiling using targeted RNA-sequencing to elucidate the progression from histologically normal lung tissues to non-invasive lesions in invasive lung adenocarcinoma

Lung adenocarcinoma (LUAD) shows heterogeneous morphological features and the stepwise progression from adenocarcinoma in situ to minimally invasive adenocarcinoma to invasive LUAD. Although multiple genetic alterations have been linked to the progression, the differences between the gene expression profiles of non-invasive lesions (non-ILs) and adjacent histologically normal lung (aNL) tissues within invasive LUAD have not been investigated. Herein, we analyzed differentially expressed genes (DEGs) specific to early-stage carcinogenesis in LUAD. Invasive LUAD tissue samples containing both non-ILs and aNL tissues were obtained from seven patients with pathological stage I LUAD, and each component was subjected to microdissection. Gene expression profiles of each component were determined using targeted RNA-sequencing. In total, 2536 DEGs, including 863 upregulated and 1673 downregulated genes, were identified in non-ILs. In non-ILs, the expression of SLC44A5, a choline transporter-like protein-coding gene, was significantly upregulated, and that of TMEM100, a gene encoding a transmembrane protein, was significantly downregulated. Reportedly, SLC44A5 plays an important role in the development and progression of hepatocellular carcinoma, whereas TMEM100 functions as a tumor suppressor in non-small cell lung cancer. Gene set enrichment analysis showed that DEGs in non-ILs were negatively enriched in cell death and immune response. Immunohistochemical analysis revealed that increased SLC44A5 expression and decreased TMEM100 expression were maintained in ILs. A protein-protein interaction (PPI) network analysis identified several upregulated and downregulated hub genes with high degrees in non-ILs. In conclusion, several new DEGs and key PPI network hub genes were identified in non-ILs, contributing to understanding of early-stage carcinogenesis in LUAD.

Genomic and transcriptomic correlates of Richter transformation in chronic lymphocytic leukemia

The transformation of chronic lymphocytic leukemia (CLL) to high-grade B-cell lymphoma is known as Richter syndrome (RS), a rare event with dismal prognosis. In this study, we conducted whole-genome sequencing (WGS) of paired circulating CLL (PB-CLL) and RS biopsies (tissue-RS) from 17 patients recruited into a clinical trial (CHOP-O). We found that tissue-RS was enriched for mutations in poor-risk CLL drivers and genes in the DNA damage response (DDR) pathway. In addition, we identified genomic aberrations not previously implicated in RS, including the protein tyrosine phosphatase receptor (PTPRD) and tumor necrosis factor receptor-associated factor 3 (TRAF3). In the noncoding genome, we discovered activation-induced cytidine deaminase-related and unrelated kataegis in tissue-RS affecting regulatory regions of key immune-regulatory genes. These include BTG2, CXCR4, NFATC1, PAX5, NOTCH-1, SLC44A5, FCRL3, SELL, TNIP2, and TRIM13. Furthermore, differences between the global mutation signatures of pairs of PB-CLL and tissue-RS samples implicate DDR as the dominant mechanism driving transformation. Pathway-based clonal deconvolution analysis showed that genes in the MAPK and DDR pathways demonstrate high clonal-expansion probability. Direct comparison of nodal-CLL and tissue-RS pairs from an independent cohort confirmed differential expression of the same pathways by RNA expression profiling. Our integrated analysis of WGS and RNA expression data significantly extends previous targeted approaches, which were limited by the lack of germline samples, and it facilitates the identification of novel genomic correlates implicated in RS transformation, which could be targeted therapeutically. Our results inform the future selection of investigative agents for a UK clinical platform study. This trial was registered at www.clinicaltrials.gov as #NCT03899337.

Tea and coffee consumption in relation to DNA methylation in four European cohorts

Lifestyle factors, such as food choices and exposure to chemicals, can alter DNA methylation and lead to changes in gene activity. Two such exposures with pharmacologically active components are coffee and tea consumption. Both coffee and tea have been suggested to play an important role in modulating disease-risk in humans by suppressing tumour progression, decreasing inflammation and influencing estrogen metabolism. These mechanisms may be mediated by changes in DNA methylation. To investigate if DNA methylation in blood is associated with coffee and tea consumption, we performed a genome-wide DNA methylation study for coffee and tea consumption in four European cohorts (N = 3,096). DNA methylation was measured from whole blood at 421,695 CpG sites distributed throughout the genome and analysed in men and women both separately and together in each cohort. Meta-analyses of the results and additional regional-level analyses were performed. After adjusting for multiple testing, the meta-analysis revealed that two individual CpG-sites, mapping to DNAJC16 and TTC17, were differentially methylated in relation to tea consumption in women. No individual sites were associated with men or with the sex-combined analysis for tea or coffee. The regional analysis revealed that 28 regions were differentially methylated in relation to tea consumption in women. These regions contained genes known to interact with estradiol metabolism and cancer. No significant regions were found in the sex-combined and male-only analysis for either tea or coffee consumption.