Single Nucleotide Polymorphisms at a Distance from Aryl Hydrocarbon Receptor (AHR) Binding Sites Influence AHR Ligand-Dependent Gene Expression

Novel Oral Adjuvant to Enhance Cytotoxic Memory-Like NK Cell Responses in an HIV Vaccine Platform

Antibody-dependent cell-mediated cytotoxicity, mediated by natural killer (NK) cells and antibodies, emerged as a secondary correlate of protection in the RV144 HIV vaccine clinical trial, the only vaccine thus far demonstrating some efficacy in human. Therefore, leveraging NK cells with enhanced cytotoxic effector responses may bolster vaccine induced protection against HIV. Here, we investigated the effect of orally administering indole-3-carbinol (I3C), an aryl hydrocarbon receptor (AHR) agonist, as an adjuvant to an RV144-like vaccine platform in a mouse model. We demonstrate the expansion of KLRG1-expressing NK cells induced by the vaccine together with I3C. This NK cell subset exhibited enhanced vaccine antigen-specific cytotoxic memory-like features. Our study underscores the potential of incorporating I3C as an oral adjuvant to HIV vaccine platforms to enhance antigen-specific (memory-like) cytotoxicity of NK cells against HIV-infected cells. This approach may contribute to enhancing the protective efficacy of HIV preventive vaccines against HIV acquisition.

Comprehensive pan-cancer analysis of 33 human cancers reveals immunotherapeutic value of focal adhesion tyrosine kinase

The immune environment in tumors is the key factor affecting the survival and immunotherapeutic response of patients. This research aimed to explore the underlying association between focal adhesion tyrosine kinase (FAK/PTK2) and cancer immunotherapy in 33 human cancers. Gene expression data and clinical features of 33 cancers were retrieved from the Cancer Genome Atlas Database. The immunotherapy cohorts included GSE67501, GSE78220, and IMVIGOR210, which were derived from the comprehensive gene expression database or from previous studies. Clinical parameters including patient age, gender, survival rate, and tumor stage were analyzed to evaluate the prognostic value of FAK/PTK2. FAK/PTK2 activity was detected by single-sample gene set enrichment analysis and used to compare the difference between FAK/PTK2 transcriptome and protein expression levels. To better understand the role of FAK/PTK2 in cancer immunotherapy, we analyzed its correlations with tumor microenvironment and with immune processes/elements (e.g., immune cell infiltration, immunosuppressants, and stimulants) and major histocompatible complexes. Potential pathways associated with FAK/PTK2 signaling in cancers were also explored. Correlations between FAK/PTK2 and 2 immunotherapeutic biomarkers (tumor mutation load and microsatellite instability) were studied. Finally, the 3 independent immunotherapy cohorts were used to study the relationship between FAK/PTK2 and immunotherapeutic response. Although FAK/PTK2 is not closely associated with age (13/33), gender (5/33), or tumor stage (5/33) in any of the studied human cancers, it has potential prognostic value for predicting patient survival. Consistency between FAK/PTK2 activity and expression exists in some cancers (3/33). Generally, FAK/PTK2 is robustly correlated with immune cell infiltration, immune modulators, and immunotherapeutic markers. Moreover, high FAK/PTK2 expression is significantly related to immune-relevant pathways. However, FAK/PTK2 is not significantly correlated with the immunotherapeutic response. Research on the immunotherapeutic value of FAK/PTK2 in 33 human cancers provides evidence regarding the function of FAK/PTK2 and its role in clinical treatment. However, given the use of a bioinformatics approach, our results are preliminary and require further validation.

Cytochrome P450 Transcriptional Regulation by Testis-Specific Y-Encoded-Like Protein: Identification of Novel Upstream Transcription Factors

Cytochrome P450s (CYPs) display significant inter-individual variation in expression, much of which remains unexplained by known CYP single-nucleotide polymorphisms (SNPs). Testis-specific Y-encoded-like proteins (TSPYLs) are transcriptional regulators for several drug-metabolizing CYPs including CYP3A4 However, transcription factors (TFs) that might influence CYP expression through an effect on TSPYL expression are unknown. Therefore, we studied regulators of TSPYL expression in hepatic cell lines and their possible SNP-dependent variation. Specifically, we identified candidate TFs that might influence TSPYL expression using the ENCODE ChIPseq database. Subsequently, the expression of TSPYL1/2/4 as well as that of selected CYP targets for TSPYL regulation were assayed in hepatic cell lines before and after knockdown of TFs that might influence CYP expression through TSPYL-dependent mechanisms. Those results were confirmed by studies of TF binding to TSPYL1/2/4 gene promoter regions. In hepatic cell lines, knockdown of the REST and ZBTB7A TFs resulted in decreased TSPYL1 and TSPYL4 expression and increased CYP3A4 expression, changes reversed by TSPYL1/4 overexpression. Potential binding sites for REST and ZBTB7A on the promoters of TSPYL1 and TSPYL4 were confirmed by chromatin immunoprecipitation. Finally, common SNP variants in upstream binding sites on the TSPYL1/4 promoters were identified and luciferase reporter constructs confirmed SNP-dependent modulation of TSPYL1/4 gene transcription. In summary, we identified REST and ZBTB7A as regulators of the expression of TSPYL genes which themselves can contribute to regulation of CYP expression and-potentially-of drug metabolism. SNP-dependent modulation of TSPYL transcription may contribute to individual variation in both CYP expression and-downstream-drug response phenotypes. SIGNIFICANCE STATEMENT: Testis-specific Y-encoded-like proteins (TSPYLs) are transcriptional regulators of cytochrome P450 (CYP) gene expression. Here, we report that variation in TSPYL expression as a result of the effects of genetically regulated TSPYL transcription factors is an additional factor that could result in downstream variation in CYP expression and potentially, as a result, variation in drug biotransformation.

Glucocorticoids unmask silent non-coding genetic risk variants for common diseases

Understanding the function of non-coding genomic sequence variants represents a challenge for biomedicine. Many diseases are products of gene-by-environment interactions with complex mechanisms. This study addresses these themes by mechanistic characterization of non-coding variants that influence gene expression only after drug or hormone exposure. Using glucocorticoid signaling as a model system, we integrated genomic, transcriptomic, and epigenomic approaches to unravel mechanisms by which variant function could be revealed by hormones or drugs. Specifically, we identified cis-regulatory elements and 3D interactions underlying ligand-dependent associations between variants and gene expression. One-quarter of the glucocorticoid-modulated variants that we identified had already been associated with clinical phenotypes. However, their affected genes were 'unmasked' only after glucocorticoid exposure and often with function relevant to the disease phenotypes. These diseases involved glucocorticoids as risk factors or therapeutic agents and included autoimmunity, metabolic and mood disorders, osteoporosis and cancer. For example, we identified a novel breast cancer risk gene, MAST4, with expression that was repressed by glucocorticoids in cells carrying the risk genotype, repression that correlated with MAST4 expression in breast cancer and treatment outcomes. These observations provide a mechanistic framework for understanding non-coding genetic variant-chemical environment interactions and their role in disease risk and drug response.

Identification of Two Genetic Loci Associated with Leukopenia after Chemotherapy in Patients with Breast Cancer

PURPOSE

To identify molecular predictors of grade 3/4 neutropenic or leukopenic events (NLE) after chemotherapy using a genome-wide association study (GWAS).

EXPERIMENTAL DESIGN

A GWAS was performed on patients in the phase III chemotherapy study SUCCESS-A (n = 3,322). Genotyping was done using the Illumina HumanOmniExpress-12v1 array. Findings were functionally validated with cell culture models and the genotypes and gene expression of possible causative genes were correlated with clinical treatment response and prognostic outcomes.

RESULTS

One locus on chromosome 16 (rs4784750; NLRC5; P = 1.56E-8) and another locus on chromosome 13 (rs16972207; TNFSF13B; P = 3.42E-8) were identified at a genome-wide significance level. Functional validation revealed that expression of these two genes is altered by genotype-dependent and chemotherapy-dependent activity of two transcription factors. Genotypes also showed an association with disease-free survival in patients with an NLE.

CONCLUSIONS

Two loci in NLRC5 and TNFSF13B are associated with NLEs. The involvement of the MHC I regulator NLRC5 implies the possible involvement of immuno-oncological pathways.

Hexokinase 2 is a transcriptional target and a positive modulator of AHR signalling

The aryl hydrocarbon receptor (AHR) regulates the expression of numerous genes in response to activation by agonists including xenobiotics. Although it is well appreciated that environmental signals and cell intrinsic features may modulate this transcriptional response, how it is mechanistically achieved remains poorly understood. We show that hexokinase 2 (HK2) a metabolic enzyme fuelling cancer cell growth, is a transcriptional target of AHR as well as a modulator of its activity. Expression of HK2 is positively regulated by AHR upon exposure to agonists both in human cells and in mice lung tissues. Conversely, over-expression of HK2 regulates the abundance of many proteins involved in the regulation of AHR signalling and these changes are linked with altered AHR expression levels and transcriptional activity. HK2 expression also shows a negative correlation with AHR promoter methylation in tumours, and these tumours with high HK2 expression and low AHR methylation are associated with a worse overall survival in patients. In sum, our study provides novel insights into how AHR signalling is regulated which may help our understanding of the context-specific effects of this pathway and may have implications in cancer.

A Comprehensive Pan-Cancer Analysis of 33 Human Cancers Reveals the Immunotherapeutic Value of Aryl Hydrocarbon Receptor

Background

Previous studies have reported the potential of aryl hydrocarbon receptor (AhR) in cancer immunotherapy. However, the mechanisms underpinning its therapeutic value have yet to be comprehensively investigated. Thus, this research aimed to explore the underlying association between AhR and cancer immunotherapy in 33 human cancers.

Methods

The gene expression data and clinical characteristics of 33 cancers were retrieved from The Cancer Genome Atlas database. The immunotherapeutic cohorts included GSE67501 and GSE78220 as well as IMvigor210, which were obtained from the Gene Expression Omnibus database and included in a previously published study respectively. Clinical parameters, including patient age, gender, survival, and tumor stage were analyzed to assess the prognostic value of AhR. The activity of AhR was generated by single sample gene set enrichment analysis and used to evaluate the difference between the AhR transcriptome and protein expression level. To better understand the role of AhR in cancer immunotherapy, the correlation between AhR and tumor microenvironment, as well as its relation to immune processes/elements, such as immune cell infiltration, immune inhibitors and stimulators, and the major histocompatibility complex were analyzed. The relevant underlying pathways associated with AhR signaling in cancer were also explored. Furthermore, the correlation between AhR and two immunotherapeutic biomarkers (tumor mutational burden and microsatellite instability) was investigated. Finally, the relationship between AhR and immunotherapeutic response was explored using three independent immunotherapeutic cohorts.

Results

Although AhR was not closely associated with age (5/33), gender (3/33), or tumor stage (3/21) in any of the studied human cancers, it exhibited potential prognostic value for predicting patient survival. Consistency has been observed between AhR activity and expression in some cancers (7/33). Generally, AhR presented a robust correlation with immune cell infiltration, immune modulators, and immunotherapeutic markers. Moreover, high AhR expression was significantly related to immune-relevant pathways. However, no significant correlation was observed between AhR and the immunotherapeutic response.

Conclusions

This research investigated the immunotherapeutic value of AhR in 33 human cancers, providing evidence regarding the function of AhR and its role in clinical treatment. However, considering that a bioinformatics approach was adopted, the current results are preliminary and require further validation.

Assessing the consequences of environmental exposures on the expression of the human receptor and proteases involved in SARS-CoV-2 cell-entry

The role of environmental condition on the infection by the novel pathogenic SARS-CoV-2 virus remains uncertain. In here, exploiting a large panel of publicly available genome-wide data, we investigated whether the human receptor ACE2 and human proteases TMPRSS2, FURIN and CATHEPSINs (B, L and V), which are involved in SARS-CoV-2 cell entry, are transcriptionally regulated by environmental cues. We report that more than 50 chemicals modulate the expression of ACE2 or human proteases important for SARS-CoV-2 cell entry. We further demonstrate that transcription factor AhR, which is commonly activated by pollutants, binds to the promoter of TMPRSS2 and enhancers and/or promoters of Cathepsin B, L and V encoding genes. Our exploratory study documents an influence of environmental exposures on the expression of genes involved in SARS-CoV-2 cell entry. These results could be conceptually and medically relevant to our understanding of the COVID-19 disease, and should be further explored in laboratory and epidemiologic studies.

TCF7L2 lncRNA: a link between bipolar disorder and body mass index through glucocorticoid signaling

Bipolar disorder (BD) and obesity are highly comorbid. We previously performed a genome-wide association study (GWAS) for BD risk accounting for the effect of body mass index (BMI), which identified a genome-wide significant single-nucleotide polymorphism (SNP) in the gene encoding the transcription factor 7 like 2 (TCF7L2). However, the molecular function of TCF7L2 in the central nervous system (CNS) and its possible role in the BD and BMI interaction remained unclear. In the present study, we demonstrated by studying human induced pluripotent stem cell (hiPSC)-derived astrocytes, cells that highly express TCF7L2 in the CNS, that the BD-BMI GWAS risk SNP is associated with glucocorticoid-dependent repression of the expression of a previously uncharacterized TCF7L2 transcript variant. That transcript is a long non-coding RNA (lncRNA-TCF7L2) that is highly expressed in the CNS but not in peripheral tissues such as the liver and pancreas that are involved in metabolism. In astrocytes, knockdown of the lncRNA-TCF7L2 resulted in decreased expression of the parent gene, TCF7L2, as well as alterations in the expression of a series of genes involved in insulin signaling and diabetes. We also studied the function of TCF7L2 in hiPSC-derived astrocytes by integrating RNA sequencing data after TCF7L2 knockdown with TCF7L2 chromatin-immunoprecipitation sequencing (ChIP-seq) data. Those studies showed that TCF7L2 directly regulated a series of BD risk genes. In summary, these results support the existence of a CNS-based mechanism underlying BD-BMI genetic risk, a mechanism based on a glucocorticoid-dependent expression quantitative trait locus that regulates the expression of a novel TCF7L2 non-coding transcript.

Prognostic value of preoperative serum bilirubin levels in ovarian cancer

Bilirubin is a promising prognostic factor for non-liver disease-related deaths in various cancers. We investigated the association between preoperative serum bilirubin levels and oncological outcomes in patients with ovarian cancer. We retrospectively analyzed the clinical data of 282 patients with epithelial ovarian carcinoma (EOC), and grouped them according to optimal threshold values of total bilirubin (TBIL), direct bilirubin (DBIL), and indirect bilirubin (IBL) measured by receiver operating characteristic curve analysis. Univariate and multivariate Cox proportional hazards regression analyses were used to evaluate various parameters that might affect overall survival (OS) and progression-free survival (PFS) in patients with EOC. The optimal cutoff values for TBIL, DBIL, and IBIL levels were 9.65 µmol/L, 2.95 µmol/L, and 6.75 µmol/L, respectively. Increased TBIL, DBIL, and IBIL levels correlated with the serum carbohydrate antigen (CA)-125 levels, International Federation of Gynecology and Obstetrics stage, and pathological differentiation (all P<0.05). Univariate analysis revealed longer OS and PFS in patients with high TBIL (≥9.65 µmol/L) and IBIL (≥6.75 µmol/L) levels (P<0.05). Multivariate analysis showed that patients with high IBIL levels (≥6.75 µmol/L) had significantly longer OS and PFS than those with low IBIL levels (<6.75 µmol/L) [hazard ratio (HR) = 0.333, 95% confidence interval (CI): 0.123~0.904, P<0.05; HR = 1.814, 95% CI: 1.169~2.816, P<0.05]. Therefore, IBIL is a potential independent prognostic factor for OS and PFS in patients with EOC. The higher the IBL level, the better the prognosis of patients with EOC.