Multiple Gene Sets Research Articles

Gene expression variation across genetically identical individuals predicts reproductive traits.

In recent decades, genome-wide association studies (GWAS) have been the major approach to understand the biological basis of individual differences in traits and diseases. However, GWAS approaches have limited predictive power to explain individual differences, particularly for complex traits and diseases in which environmental factors play a substantial role in their etiology. Indeed, individual differences persist even in genetically identical individuals, although fully separating genetic and environmental causation is difficult in most organisms. To understand the basis of individual differences in the absence of genetic differences, we measured two quantitative reproductive traits in 180 genetically identical young adult Caenorhabditis elegans roundworms in a shared environment and performed single-individual transcriptomics on each worm. We identified hundreds of genes for which expression variation was strongly associated with reproductive traits, some of which depended on individuals' historical environments and some of which was random. Multiple small sets of genes together were highly predictive of reproductive traits, explaining on average over half and over a quarter of variation in the two traits. We manipulated mRNA levels of predictive genes to identify a set of causal genes, demonstrating the utility of this approach for both prediction and understanding underlying biology. Finally, we found that the chromatin environment of predictive genes was enriched for H3K27 trimethylation, suggesting that gene expression variation may be driven in part by chromatin structure. Together, this work shows that individual, non-genetic differences in gene expression are both highly predictive and causal in shaping reproductive traits.

Isolation, Sphalerite Bioleaching, and Whole Genome Sequencing of Acidithiobacillus ferriphilus QBS3 from Zinc-Rich Sulfide Mine Drainage.

The genus Acidithiobacillus has been widely used in bioleaching, and novel strains in this genus, such as A. ferriphilus, have also been confirmed to possess bioleaching capabilities. In this study, an Acidithiobacillus ferriphilus strain, QBS3, was isolated from zinc-rich sulfide mine drainage using the gradient dilution method. QBS3 is a Gram-negative, 1.3 µm rod-shaped bacterium with small red colonies. It showed a high iron oxidation efficiency of 0.361 g/(L·h) and a sulfur oxidation efficiency of 0.206 g/(L·d). QBS3 has sphalerite bioleaching ability; using QBS3 for pure sphalerite bioleaching, 18.8% of zinc was extracted in 14 days at 1% pulp density. Whole genome sequencing was performed on QBS3. Functional prediction showed that 9.13% of the genes were involved in replication, recombination, and repair. Bioleaching-related genes were analyzed, including iron and sulfur oxidation genes, and carbon and nitrogen fixation genes. For iron oxidation, the Cyc2→RusA pathway and Iro→RusB pathway were found in QBS3. In terms of sulfur oxidation, QBS3 has an incomplete SOX system and lacks the SDO gene, but Rho and Trx may complement the SOX system, enabling QBS3 to oxidize sulfur. QBS3 has multiple sets of carbon fixation genes, and nitrogen fixation genes were also identified. A hypothetical sphalerite bioleaching model is proposed; this study provides a theoretical basis for the zinc sulfide ore bioleaching industry.

Abstract 1350: ESR1 expression as a predictive biomarker for immune checkpoint inhibitor efficacy in estrogen receptor-positive, HER2-negative breast cancer

Abstract Introduction: Immune checkpoint inhibitors (ICIs) hold promise for preventing late recurrence, a major cause of mortality in estrogen receptor-positive/HER2-negative (ER+/HER2-) subtype, which constitutes the majority of breast cancer (BC) cases. The response to ICIs is worse in this subtype compared to triple-negative breast cancer (TNBC). Recent clinical trials, CheckMate 7FL and KEYNOTE-756, have uniformly demonstrated that ER+/HER2- BC with less Estrogen Receptor (ER)-α staining in cancer cells shows a better response to ICIs. Considering the relatively precise quantification achievable through transcriptomic analysis, we hypothesized that ER-α gene (ESR1) expression level is a negative predictive biomarker of ER+/HER2- BC response to ICIs. Methods: Total of 4138 BC patients from ten cohorts (GSE271080, I-SPY2 [GSE173839], TCGA, METABRIC, SCAN-B [GSE96058], GSE246613, GSE255107/GSE248288, GSE161529, GSE167036) were analyzed. High expression of ESR1 was defined as the top two-thirds of expression levels for each gene. Results: In the ER+/HER2- subtype, high ESR1 expression was associated with significantly less intratumoral-Tumor Infiltrating Lymphocytes in both TCGA and METABRIC cohorts (p<0.001). High ESR1 expression was also associated with less infiltrations of CD8+ T cells, dendritic cells (DCs), cancer-associated fibroblasts (CAFs), a lower CD8/Treg ratio, and a lower cytolytic activity consistently in TCGA, METABRIC, SCAN-B cohorts (all p<0.01). Further, the low ESR1 group significantly enriched multiple immune-related gene sets, including IL2 signaling, TNF-α signaling, and interferon-γ response. In the I-SPY2 trial, high ESR1 expression correlated with lower pathological complete response rates after neoadjuvant immunotherapy (n=50, p<0.001). The predictive accuracy of ESR1 expression (AUC=0.88) exceeded that of PDL1 (AUC=0.60) in HR+/HER2- BC. Single-cell analysis (GSE246613) revealed negative associations between ICI response and ESR1 expression in cancer/normal epithelial cells, endothelial cells, and plasmablasts, as well as activation of estrogen signaling in these cell types. Notably, CAFs showed a positive association between ESR1 expression, estrogen signaling activation, and ICI response. Conclusion: ESR1 expression is a promising predictive biomarker for ICI response in ER+/HER2- BC. Additionally, CAF infiltration and estrogen signaling activation in CAFs may enhance ICI response, highlighting their potential as therapeutic targets. Citation Format: Jun Arima, Kohei Chida, Rongrong Wu, Kohei Taniguchi, Mitsuhiko Iwamoto, Takashi Ishikawa, Kenichi Hakamada, Sang-Woong Lee, Kazuaki Takabe. ESR1 expression as a predictive biomarker for immune checkpoint inhibitor efficacy in estrogen receptor-positive, HER2-negative breast cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2025; Part 1 (Regular Abstracts); 2025 Apr 25-30; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2025;85(8_Suppl_1):Abstract nr 1350.

Hybrid lizards with introgressed MtDNA show increased resistance to DNA damage from reactive oxygen species

Hybrids can experience genetic incompatibilities that lead to cellular dysfunction and infertility. Selection against hybrids can reinforce species boundaries, while compensatory mechanisms that support hybrid survival can lead to gene flow across species boundaries. Previous analyses of mitochondrial function of the hybrid lizards between Urosaurus graciosus and Urosaurus ornatus, which have introgressed mtDNA, showed evidence of mitochondrial dysfunction (e.g. higher oxygen consumption during ATP production) likely leading to higher reactive oxygen species (ROS) production, one of the main causes of damage to DNA. Hybrid populations of these lizards have persisted so we predicted they may have mechanisms for protecting their DNA from damage. To test this hypothesis, we exposed liver cells of these lizards to hydrogen peroxide, a source of ROS, and compared DNA damage using single-cell gel electrophoresis. Cells from hybrid individuals showed significantly less DNA damage. We then examined potential mechanisms to explain this result by comparing melanin content of livers and expression levels of antioxidant and DNA repair genes using gene set enrichment analyses. Melanin content of livers was highly variable among U. graciosus individuals and therefore not highest in hybrids as expected but was lowest in the U. ornatus which had the highest levels of DNA damage. The overall gene expression profiles of hybrids and U. graciosus were very similar while that of the mtDNA donor species, U. ornatus, was highly divergent from both. Three DNA repair and antioxidant related gene sets showed significant levels of enrichment in hybrids (Reactome gap filling DNA repair synthesis and ligation in global genome nucleotide excision repair, Glutathione metabolism, and GOMF Antioxidant activity). Antioxidant genes contributing to high enrichment scores in multiple gene sets in hybrids warrant further study. These results show how species can compensate for mtDNA introgression and highlight the complexity of genetic interactions encountered when closely related species exchange genetic material in zones of secondary contact.

Drinking motives and alcohol sensitivity mediate multidimensional genetic influences on alcohol use behaviors.

Genetic influences account for a substantial proportion of individual differences in alcohol use behaviors (AUBs). However, multiple distinct sets of genes are linked to different AUBs via uncertain causal links. Here, we explore whether intermediate neurobiological traits mediate the relationship between polygenic scores (PGSs) and multiple AUBs, with the aim to better understand processes captured by different genetic profiles. We derived four alcohol-related PGSs in participants from Spit for Science, a longitudinal study of college students in the United States (n = 4549). Using linear regression, we tested the relationship between PGSs and 22 potential mediators, including personality, alcohol expectancies, drinking motives, and alcohol sensitivity. Nominally significant effects were carried forward to a multiple mediation model to estimate direct and indirect effects on four measured AUBs (frequency, quantity, alcohol use disorder symptoms [AUDsx], and maximum drinks in 24 h). In univariable regression, PGSs indexing genetic effects on drinks per week (DPW) and problematic alcohol use (PAU) predicted higher levels of impulsivity and drinking motives as well as lower alcohol sensitivity. BeerPref PGSs (indexing a variable pattern of alcohol problems and preference for beer) predicted higher negative urgency and lower alcohol sensitivity. Mediational models indicated direct and indirect effects of DPW PGSs on multiple AUBs via social/enhancement drinking motives and alcohol sensitivity, indirect effects of PAU PGSs on AUDsx, and indirect effects of BeerPref PGS on drinking frequency and AUDsx via the joint effect of mediators including alcohol sensitivity. These findings provide evidence that the genetic influences on AUBs are associated with and partially mediated by intermediate neurobiological and cognitive factors, which may be more amenable to intervention. Greater focus on drinking motives and alcohol sensitivity is warranted in genetic research, as well as attention to the heterogeneous pathways linking genes to alcohol use outcomes.

5-hydroxymethylcytosine sequencing of plasma cell-free DNA identifies epigenomic features in prostate cancer patients receiving androgen deprivation therapies

BackgroundWe evaluated whether 5hmC signatures in cell-free DNA (cfDNA) are associated with treatment failure to androgen-deprivation therapies (ADT) among men with hormone-naive prostate cancer.MethodsWe collected a total of 139 serial plasma samples from 55 prostate cancer patients receiving ADT at 3 time points including baseline (before initiating ADT, n = 55); 3 months (after initiating ADT, n = 55); and disease progression (n = 15) within 24 months or 24 months if no progression was detected (n = 14). We used selective chemical labeling sequencing to quantify 5hmC abundance across the genome and Kaplan–Meier analysis to assess survival association.ResultsHere we show a significant 5hmC difference in 1642 of 23433 genes between patients with and without progression (false discovery rate [FDR] < 0.1) in baseline plasma samples. Patients with progression demonstrate significant 5hmC enrichments in multiple hallmark gene sets, with androgen responses as the top enriched gene-set (FDR = 1.19E−13). We further show a significant association between high activity scores in these gene sets and poor progression-free survival (P < 0.05), even after adjusting for circulating tumor DNA fraction and prostate-specific antigen values. Additionally, our longitudinal analysis shows that the high activity score is significantly reduced after 3 months of initiating ADT (P = 0.0004) but returns to higher levels when the disease progresses (P = 0.0317).Conclusions5hmC-based activity scores from gene-sets involved in AR responses show great potential in assessing treatment resistance, monitoring disease progression, and identifying patients who would benefit from upfront treatment intensification. However, further studies are needed to validate these findings.

Analysis of fatty acid metabolism in prostate cancer and discovery of a new programmed cell death-associated luminal cell subpopulation

IntroductionThis study focuses on the role of fatty acid metabolism in prostate cancer, particularly in oncogenic luminal cells associated with programmed cell death under the influence of metabolic reprogramming.Materials and MethodsProstate cancer was analyzed using single-cell transcriptomics and spatial transcriptomics data. Fatty acid metabolism levels in the tumor microenvironment were quantified by multiple gene set scoring methods, and data were processed using NMF and deconvolution methods to identify different cell populations and their interactions in the tumor microenvironment.ResultsLuminal cells have significantly increased activity in fatty acid metabolism, which is associated with the aggressiveness and metastatic capability of tumors. Luminal cell subpopulations have been found to play a key role in the development of prostate cancer, especially their close association with programmed cell death.ConclusionThis study deepens the understanding of the role of fatty acid metabolism in prostate cancer, identifies fatty acid metabolism-related luminal cell subtypes, and proposes new therapeutic targets, providing new insights into prostate cancer treatment.

GeneSetCart: assembling, augmenting, combining, visualizing, and analyzing gene sets.

Converting multiomics datasets into gene sets facilitates data integration that leads to knowledge discovery. Although there are tools developed to analyze gene sets, only a few offer the management of gene sets from multiple sources. GeneSetCart is an interactive web-based platform that enables investigators to gather gene sets from various sources; augment these sets with gene-gene coexpression correlations and protein-protein interactions; perform set operations on these sets such as union, consensus, and intersection; and visualize and analyze these gene sets, all in one place. GeneSetCart supports the upload of single or multiple gene sets, as well as fetching gene sets by searching PubMed for genes comentioned with terms in publications. Venn diagrams, heatmaps, Uniform Manifold Approximation and Projection (UMAP) plots, SuperVenn diagrams, and UpSet plots can visualize the gene sets in a GeneSetCart session to summarize the similarity and overlap among the sets. Users of GeneSetCart can also perform enrichment analysis on their assembled gene sets with external tools. All gene sets in a session can be saved to a user account for reanalysis and sharing with collaborators. GeneSetCart has a gene set library crossing feature that enables analysis of gene sets created from several National Institutes of Health Common Fund programs. For the top overlapping sets from pairs of programs, a large language model (LLM) is prompted to propose possible reasons for the high overlap. Using this feature, two use cases are presented. In addition, users of GeneSetCart can produce publication-ready reports from their uploaded sets. Text in these reports is also supplemented with an LLM. Overall, GeneSetCart is a useful resource enabling biologists without programming expertise to facilitate data integration for hypothesis generation.

Low tristetraprolin expression activates phenotypic plasticity and primes transition to lethal prostate cancer in mice.

Phenotypic plasticity is a hallmark of cancer and is increasingly realized as a mechanism of resistance to androgen receptor-targeted (AR-targeted) therapy. Now that many prostate cancer (PCa) patients are treated upfront with AR-targeted agents, it is critical to identify actionable mechanisms that drive phenotypic plasticity, to prevent the emergence of resistance. We showed that loss of tristetraprolin (TTP; gene ZFP36) increased NF-κB activation, and was associated with higher rates of aggressive disease and early recurrence in primary PCa. We also examined the clinical and biological impact of ZFP36 loss with co-loss of PTEN, a known driver of PCa. Analysis of multiple independent primary PCa cohorts demonstrated that PTEN and ZFP36 co-loss was associated with increased recurrence risk. Engineering prostate-specific Zfp36 deletion in vivo induced prostatic intraepithelial neoplasia, and, with Pten codeletion, resulted in rapid progression to castration-resistant adenocarcinoma. Zfp36 loss altered the cell state driven by Pten loss, as demonstrated by enrichment of epithelial-mesenchymal transition (EMT), inflammation, TNF-α/NF-κB, and IL-6-JAK/STAT3 gene sets. Additionally, our work revealed that ZFP36 loss also induced enrichment of multiple gene sets involved in mononuclear cell migration, chemotaxis, and proliferation. Use of the NF-κB inhibitor dimethylaminoparthenolide (DMAPT) induced marked therapeutic responses in tumors with PTEN and ZFP36 co-loss and reversed castration resistance.

5-hydroxymethylcytosine sequencing in plasma cell-free DNA identifies unique epigenomic features in prostate cancer patients resistant to androgen deprivation therapies.

Currently, no biomarkers are available to identify resistance to androgen-deprivation therapies (ADT) in men with hormone-naive prostate cancer. Since 5-hydroxymethylcytosines (5hmC) in gene body are associated with gene activation, in this study, we evaluated whether 5hmC signatures in cell-free DNA (cfDNA) predicts early resistance to ADT. We collected a total of 139 serial plasma samples from 55 prostate cancer patients receiving ADT at three time points including baseline (prior to initiating ADT, N=55), 3-month (after initiating ADT, N=55), and disease progression (N=15) within 24 months or 24-month if no progression was detected (N=14). To quantify 5hmC abundance across the genome, we used selective chemical labeling sequencing and mapped sequence reads to individual genes. Differential methylation analysis in baseline samples identified significant 5hmC difference in 1,642 of 23,433 genes between patients with and without progression (false discovery rate, FDR<0.1). Patients with disease progression showed significant 5hmC enrichments in multiple hallmark gene sets with androgen responses as top enriched gene set (FDR=1.19E-13). Interestingly, this enrichment was driven by a subgroup of patients featuring a significant 5hmC hypermethylation in the gene sets involving AR , FOXA1 and GRHL2 . To quantify overall activities of these gene sets, we developed a gene set activity scoring algorithm and observed significant association of high activity scores with poor progression-free survival (P<0.05). Longitudinal analysis showed that the high activity scores were significantly reduced after 3-months of initiating ADT (P<0.0001) but returned to higher levels when the disease was progressed (P<0.05). This study demonstrates that 5hmC-based activity scores from gene sets involved in AR , FOXA1 and GRHL2 may be used as biomarkers to determine early treatment resistance, monitor disease progression, and potentially identify patients who would benefit from upfront treatment intensification.

Drinking motives and alcohol sensitivity mediate multi-dimensional genetic influences on alcohol use behaviors

Background:Genetic influences account for a substantial proportion of individual differences in alcohol use behaviors (AUBs). However, multiple distinct sets of genes are linked to different AUBs, which may explain their dramatic variability in risk factors and manifestations. In this study, we explore whether intermediate neurobiological traits and alcohol-related cognitions mediate the relationship between polygenic scores (PGS) and multiple AUBs, with the aim to better understand processes captured by different genetic profiles.Methods:Using results from prior genome-wide association studies, we derived PGS for 6 AUBs in participants from Spit for Science, a longitudinal study of college students in the U.S. (n=4,549). Self-report measures included personality traits, alcohol expectancies, drinking motivations, and alcohol sensitivity measures as well as drinking frequency, drinking quantity, alcohol use disorder (AUD) symptoms, and maximum drinks in 24 hours. Using linear regression and multiple mediation models, we investigated the direct and indirect effects of PGS on AUBs.Results:In univariable regression results, PGSs indexing broad AUB dimensions such as drinks per week (DPW) and AUD predicted higher levels of sensation-seeking and multiple drinking motives, while BeerPref PGSs (indexing a variable pattern of alcohol problems associated with a preference for beer) predicted higher negative urgency and lower alcohol sensitivity. Mediational models indicated strong direct and indirect effects of DPW PGSs on multiple AUBs via social/enhancement drinking motives and alcohol sensitivity, indirect effects of AUD PGSs on AUD symptoms via coping motives, and indirect effects of BeerPref PGS on all AUBs via the joint effect of mediators including alcohol sensitivity.Conclusions:These findings provide initial evidence that the genetic influences on different AUBs are associated with and partially mediated by intermediate neurobiological and cognitive factors, which may be more amenable to intervention. Greater focus on drinking motives and alcohol sensitivity is warranted in genetic research, as well as attention to the heterogeneous pathways linking genes to alcohol use outcomes.

High-fidelity CRISPR/Cas12a dual-crRNA screening reveals novel synergistic interactions in hepatocellular carcinoma.

: CRISPR/Cas12a-based combinational screening has shown remarkable potential for identifying genetic interactions. Here, we describe an innovative method for combinational genetic screening with rapid construction of a dual-CRISPR RNA (crRNA) library using gene splicing through overlap extension PCR (SOE PCR) and the adoption of CeCas12a, which we previously identified with strict PAM recognition and low off-targeting to guarantee fidelity and efficiency. The custom-pooled SOE crRNA array (SOCA) library for double-knockout screening could be conveniently constructed in the laboratory for widespread use, and the CeCas12a-mediated high-fidelity screen displayed good performance even under a negative selection screen. By designing a SOCA dual-crRNA library that covered most of the kinase and metabolism-associated gene targets of FDA-approved drugs implicated in hepatocellular carcinoma (HCC) tumourigenesis, novel cross-talk between the two gene sets was negatively selected to inhibit HCC cell growth in vitro and in vivo and was validated using virtual double-knockdown screening based on TCGA databases. Thus, this rapid, efficient and high-fidelity double-knockout screening system is promising for systemically identifying potential genetic interactions between multiple gene sets or combinations of FDA- approved drugs for clinical translational medicine in the future.

Using de novo transcriptomes to decipher the relationships in cutthroat trout subspecies (Oncorhynchus clarkii)

For almost 200 years, the taxonomy of cutthroat trout (Oncorhynchus clarkii), a salmonid native to Western North America, has been in flux as ichthyologists and fisheries biologists have tried to describe the diversity within these fishes. Starting in the 1950s, Robert Behnke reexamined the cutthroat trout and identified 14 subspecies based on morphological traits, Pleistocene events, and modern geographic ranges. His designations became instrumental in recognizing and preserving the remaining diversity of cutthroat trout. Over time, molecular techniques (i.e. karyotypes, allozymes, mitochondrial DNA, SNPs, and microsatellite arrays) have largely reinforced Behnke's phylogenies, but have also revealed that some relationships are consistently weakly supported. To further resolve these relationships, we generated de novo transcriptomes for nine cutthroat subspecies, as well as a Bear River Bonneville form and two Colorado River lineages (blue and green). We present phylogenies of these subspecies generated from multiple sets of orthologous genes extracted from our transcriptomes. We confirm many of the relationships identified in previous morphological and molecular studies, as well as discuss the importance of significant differences apparent in our phylogenies from these studies within a geological perspective. Specific findings include three distinct clades: (1) Bear River Bonneville form and Yellowstone cutthroat trout; (2) Bonneville cutthroat trout (n = 2); and (3) Greenback and Rio Grande cutthroat trout. We also identify potential gene transfer between Bonneville cutthroat trout and a population of Colorado River green lineage cutthroat trout. Using these findings, it appears that additional groups warrant species‐level consideration if other recent species elevations are retained.

Abstract PO2-09-05: Building genomic and transcriptomic data among African American and Black patients with triple-negative inflammatory breast cancer

Abstract Background: Inflammatory breast cancer (IBC) is an aggressive breast cancer associated with poor response and early metastases. It represents a health disparity, with higher incidence and worse outcomes in Black compared to White patients. Further, IBC is enriched for the more aggressive triple-negative receptor subtype, TN-IBC, and this subtype is more prevalent among Black IBC patients. It is critical to identify genetic and molecular features that may guide better treatments for these aggressive cancers across races, and thus imperative to identify and increase racial representation among bioinformatics datasets. We previously performed comprehensive genomic and transcriptomic analyses for TN-IBC patient samples and here examine the differences among self-identified African American or Black (AA/B) and other races among IBC patients. Methods: As previously described, we collected matched blood and baseline tumor samples before treatment from 19 patients with primary TN-IBC. We performed whole exome and RNA sequencing (RNA-Seq) on these samples and compared mutations and differentially expressed genes (DEGs). Results were compared between 3 AA/B patient tissues vs the remaining 16 (Other) cases. Results: Mutations shared by all AA/B patients included amplification of the five most commonly mutated hallmark genes in the complete cohort, ARNT, BCL9, DDR2, FCGR2B and LMNA (all 100% in AA/B patients versus all 50% in the Other cohort). All AA/B patients had a Notch1 mutation (two deletions and one missense mutation, 100% AA/B vs 6% Other (one deletion)). Comparing differentially expressed genes and gene set enrichment analyses of the Broad molecular signatures database (MSigDB) demonstrates the majority of DEGs are downregulated in these AA/B patients. Six of the top 20 downregulated gene sets from the MSigDB C2 group relate to breast cancer subtype and among this TNIBC cohort AA/B tumor expression was downregulated for luminal and normal subtype related genes in multiple gene sets and enriched in the Lien metaplastic related genes. In the MSigDB C5 sets, over half of the top 20 enriched gene ontology sets among AA/B patients related to chromosomes, chromatin, spindle assembly or histone binding, while 16 of the top 20 downregulated enriched sets related to immune function. Multiple KEGG pathways enriched in AA/B cases related to RNA processing and expression and DNA replication and repair. Lastly, hallmark gene sets enriched in AA/B patients were cell cycle related, and targets of Myc, Hedgehog, B-catenin, TGFb, and TNFa via NFKb, while top significantly downregulated sets related to metabolism and immune response. Conclusion: While findings from a small sample size can only be considered hypothesis generating, some patterns emerge from the data among this rare cohort of AA/B TN IBC patients. These tumors appear to be devoid of luminal signals, characterized by cell cycle and chromatin remodeling signals, enriched for signaling typical of stem cell self-renewal, and perhaps more immune-deprived than other TN-IBC. Additional work to build racially diverse IBC transcriptomic data are needed to develop strategies to understand and improve outcomes for these patients. Citation Format: Wencai Ma, Li Zhao, Gayathri Devi, Savitri Krishnamurthy, Larry Coffer, Angela Alexander, Megumi Kai, Xiaoping Wang, Huong Le-Petross, Miral Patel, Bisrat Debeb, Chandra Bartholomeusz, Jianhua Zhang, Xingzhi Song, Andrew Futreal, Naoto Ueno, Rachel Layman, Anthony Lucci, Jing Wang, Wendy Woodward. Building genomic and transcriptomic data among African American and Black patients with triple-negative inflammatory breast cancer [abstract]. In: Proceedings of the 2023 San Antonio Breast Cancer Symposium; 2023 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2024;84(9 Suppl):Abstract nr PO2-09-05.

Abstract PO1-24-12: Investigating the immunological function of alpha-2-glycoprotein 1, zinc-binding in regulating tumor response in the breast cancer microenvironment

Abstract Purpose: Tumor immunology has attracted considerable attention as an innovative therapeutic strategy for various types of cancer. Elucidating the unique immunoregulatory mechanisms in the breast cancer microenvironment will assist in the development of novel treatment strategies. Recently, we and other researchers have found that androgen receptor (AR) expression is associated with the immunosuppressive phenotype in the breast cancer microenvironment, suggesting some immunoregulatory function in breast cancer; however, the mechanism remains unclear. Here, we focused on AR-dependent secreted protein, alpha-2-glycoprotein 1, zinc-binding (ZAG) encoded by the AZGP1 gene in breast cancer, which is structurally similar to HLA class I and is implicated in immune regulation. In this study, we investigated the immunological function of AZGP1/ZAG in the breast cancer microenvironment. Methods: We performed a gene set enrichment analysis (GSEA) to screen the biological processes associated with AZGP1 expression using a gene expression profile dataset of the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC). Subsequently, we analyzed the correlation between AZGP1 expression and the immune cell composition in breast cancer tissues estimated with the CIBERSORTx using METABRIC and another dataset of The Sweden Cancerome Analysis Network-Breast. In our previous study of 45 breast cancer tissue samples, we evaluated the infiltration of 11 types of immune cells using flow cytometry (FCM). ZAG expression was further evaluated by immunohistochemistry, and the relationship between ZAG expression and the immune cell composition in breast cancer tissue was analyzed. Based on the results of this analysis (shown in the results section) we hypothesized that ZAG has some effect on the macrophage (Mφ). The action of ZAG in M1/M2 polarization models constructed using primary culture of human peripheral blood mononuclear cells (PBMC)-derived Mφ was assessed based on the expression of M1/M2 differentiation markers (CD86, CD80/CD163, MRC1) and HLA class I/II expression evaluated by FCM (n=15 each). Results: GSEA demonstrated that AZGP1 expression was negatively correlated with multiple gene sets representing immunological processes, including inflammatory response, allograft rejection, interferon gamma response, IL6/JAK/STAT3 signaling, complement, and IL2/STAT5 signaling. Analysis by CIBERSORTx showed that AZGP1 expression was negatively correlated with the absolute score (the absolute abundance of total immune cell infiltration), Mφ M1, NK cells activated, CD4+ T memory activated, and CD8+ T (r&amp;lt;-3 and p&amp;lt; 0.05). Analyses of our in-house dataset using breast cancer tissue showed that ZAG expression was associated with decreased infiltration of monocytes/macrophages, non-classical monocytes, and myeloid-derived suppressor cells in breast cancer tissues assessed by FCM. In the in vitro analyses, ZAG decreased the expression of CD80, CD163, MRC1, and HLA classes I and II in the M1 polarization model and the expression of CD163 and MRC1 in the M2 polarization model. Conclusion: AZGP1/ZAG was associated with an immunosuppressive phenotype and reduced infiltration of specific immune cell subsets, particularly Mφ, into breast cancer tissues. In the in vitro analysis, ZAG demonstrated some regulatory effects on the phenotypic change of Mφ. Our findings strongly suggest ZAG is a novel mediator of AR-dependent immunomodulation in the breast cancer microenvironment, laying the foundation for future studies to elucidate the immunological role of ZAG in breast cancer. Citation Format: Toru Hanamura, Kozue Yokoyama, Shigehisa Kitano, Hiroshi Kagamu, Makiko Yamashita, Mayako Terao, Takuho Okamura, Nobue Kumaki, Katsuto Hozumi, Takayuki Iwamoto, Chikako Honda, Sasagu Kurozumi, Jennifer Richer, Naoki Niikura. Investigating the immunological function of alpha-2-glycoprotein 1, zinc-binding in regulating tumor response in the breast cancer microenvironment [abstract]. In: Proceedings of the 2023 San Antonio Breast Cancer Symposium; 2023 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2024;84(9 Suppl):Abstract nr PO1-24-12.

Enhanced epithelial-mesenchymal transition signatures are linked with adverse tumor microenvironment, angiogenesis and worse survival in gastric cancer.

Epithelial-mesenchymal transition (EMT) is a crucial mechanism that facilitates cancer cell metastasis. Despite its importance, the clinical significance of EMT in gastric cancer (GC) patients has yet to be clearly demonstrated. For gauging the extent of EMT in GC, we employed gene set variation analysis to score 807 patient samples from two large cohorts: TCGA and GSE84437. In both cohorts, EMT high GC showed a significant association with worse overall survival (hazard ratio (HR) = 1.74, p = 0.011 and HR = 2.01, p < 0.001, respectively). This association was stronger when considering the EMT signature score compared to the individual expressions of EMT-related genes (CDH1, CDH2, VIM, and FN1). While the EMT signature level did not differ among various cancers, high EMT signature specifically correlated with survival in GC alone. Mucinous and diffuse histological types exhibited higher EMT levels compared to others (p < 0.001), and the EMT signature level was correlated with tumor depth and AJCC stage (all p < 0.001). Interestingly, the EMT score was an independent factor for overall and disease-specific survival (multivariate; p = 0.006 and 0.032, respectively). EMT high GC displayed a lower fraction of Th1 cells and a higher fraction of dendritic cells, M1 macrophages and several stromal cells. EMT high GC exhibited an inverse correlation with cell proliferation-related gene sets. While they significantly enriched multiple pro-cancerous gene sets, such as TGF-β signaling, hypoxia, and angiogenesis. The presence of EMT signature in a bulk tumor was linked to TGF-β signaling, hypoxia, and angiogenesis, and was also associated with poorer survival outcomes in GC patients.

Abstract 7044: Transcriptome and TCR repertoire in human telomerase reverse transcriptase peptide vaccine treated metastatic prostate cancer patients

Abstract Human telomerase reverse transcriptase (hTERT) vaccine has found to be effective for most of the participated de novo metastatic prostate cancer patients in clinical studies. Nevertheless, the molecular mechanisms behind the immune response induced by the hTERT vaccine are not fully elucidated. Among the cohort of 21 patients treated with androgen deprivation therapy, radiotherapy and hTERT vaccine, we used RNA-sequencing and TCR-sequencing (20 and 14 pre-treatment and 20 and 21 post-treatment samples, respectively) to inspect the changes in transcriptome and TCR-repertoire after the hTERT vaccine treatment and differences between patients that had early response and patients that had late or no response. We found multiple hallmarks being shifted after the treatment and being different between early and late responders. Also, a change in activity of multiple gene set signatures was detected after the treatment and activity of part of the gene set signatures predicted survival. Deconvolution of cell type composition from RNA-sequencing data revealed changes in immune cell abundances after the treatment. TCR-sequencing revealed changes in composition and diversity of TCR repertoires after the treatment and these differed between early and late responders. The results of this study will help to predict patients’ responses to the cancer vaccine and to further understand the mechanisms behind successful cancer vaccines. Citation Format: Reetta Nätkin, Julia Elisabeth Simensen, Nikolai Engedal, Espen Basmo Ellingsen, Marita Westhrin, Ingunn Westgaard, Eivind Hovig, Wolfgang Lilleby, Matti Nykter, Alfonso Urbanucci. Transcriptome and TCR repertoire in human telomerase reverse transcriptase peptide vaccine treated metastatic prostate cancer patients [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 7044.

Genome-wide transcriptional response to silver stress in extremely halophilic archaeon Haloferax alexandrinus DSM 27206 T

BackgroundThe extremely halophilic archaeon Haloferax (Hfx.) alexandrinus DSM 27206 T was previously documented for the ability to biosynthesize silver nanoparticles while mechanisms underlying its silver tolerance were overlooked. In the current study, we aimed to assess the transcriptional response of this haloarchaeon to varying concentrations of silver, seeking a comprehensive understanding of the molecular determinants underpinning its heavy metal tolerance.ResultsThe growth curves confirmed the capacity of Hfx. alexandrinus to surmount silver stress, while the SEM–EDS analysis illustrated the presence of silver nanoparticles in cultures exposed to 0.5 mM silver nitrate. The RNA-Seq based transcriptomic analysis of Hfx. alexandrinus cells exposed to 0.1, 0.25, and 0.5 mM silver nitrate revealed the differential expression of multiple sets of genes potentially employed in heavy-metal stress response, genes mostly related to metal transporters, basic metabolism, oxidative stress response and cellular motility. The RT-qPCR analysis of selected transcripts was conducted to verify and validate the generated RNA-Seq data.ConclusionsOur results indicated that copA, encoding the copper ATPase, is essential for the survival of Hfx. alexandrinus cells in silver-containing saline media. The silver-exposed cultures underwent several metabolic adjustments that enabled the activation of enzymes involved in the oxidative stress response and impairment of the cellular movement capacity. To our knowledge, this study represents the first comprehensive analysis of gene expression in halophillic archaea facing increased levels of heavy metals.

Sex-dimorphic expression of extracellular matrix genes in mouse bone marrow neutrophils.

The mammalian innate immune system is sex-dimorphic. Neutrophils are the most abundant leukocyte in humans and represent innate immunity's first line of defense. We previously found that primary mouse bone marrow neutrophils show widespread sex-dimorphism throughout life, including at the transcriptional level. Extracellular matrix [ECM]-related terms were observed among the top sex-dimorphic genes. Since the ECM is emerging as an important regulator of innate immune responses, we sought to further investigate the transcriptomic profile of primary mouse bone marrow neutrophils at both the bulk and single-cell level to understand how biological sex may influence ECM component expression in neutrophils throughout life. Here, using curated gene lists of ECM components and unbiased weighted gene co-expression network analysis [WGCNA], we find that multiple ECM-related gene sets show widespread female-bias in expression in primary mouse neutrophils. Since many immune-related diseases (e.g., rheumatoid arthritis) are more prevalent in females, our work may provide insights into the pathogenesis of sex-dimorphic inflammatory diseases.

Novel genome-wide associations for effort valuation and psychopathology in children and adults.

The Research Domain Criteria (RDoC) initiative was established by the US National Institute of Mental Health as a multilevel, disorder-agnostic framework for analysis of human psychopathology through designated domains and constructs, including the "Positive Valence Systems" domain focused on reward-related behavior. This study investigates the reward valuation subconstruct of "effort" and its association with genetic markers, functional neurobiological pathways, and polygenic risk scores for psychopathology in 1215 children aged 6-12 and their parents (n = 1044). All participants completed the effort expenditure for rewards task (EEfRT), which assesses "effort" according to two quantitative measures: hard-task choice and reward sensitivity. Genetic association analyses were undertaken in MAGMA, utilizing EEfRT outcome variables as genome-wide association studies phenotypes to compute SNP and gene-level associations. Genome-wide association analyses found two distinct genetic loci that were significantly associated with measures of reward sensitivity and a separate genetic locus associated with hard task choice. Gene-set enrichment analysis yielded significant associations between "effort" and multiple gene sets involved in reward processing-related pathways, including dopamine receptor signaling, limbic system and forebrain development, and biological response to cocaine. These results serve to establish "effort" as a relevant construct for understanding reward-related behavior at the genetic level and support the RDoC framework for assessing disorder-agnostic psychopathology.