Search

Enhancing the efficacy of immunotherapy in brain metastases (BrM) requires an improved understanding of the immune composition of BrM and how this is affected by radiation and dexamethasone. Our two-arm pilot study (NCT04895592) allocated 26 patients with BrM to either low (Arm A) or high (Arm B) dose peri-operative dexamethasone followed by pre-operative stereotactic radiosurgery (pSRS) and resection (n= 13 per arm). The primary endpoint, a safety analysis at 4 months, was met. The secondary clinical endpoints of overall survival, distant brain failure, leptomeningeal disease and local recurrence at 12-months were 66%, 37.3%, 6%, and 0% respectively and were not significantly different between arms (p= 0.7739, p= 0.3884, p= 0.3469). Immunological data from two large retrospective BrM datasets and confirmed by correlates from both arms of this pSRS prospective trial revealed that BrM CD8 T cells were composed of predominantly PD1+ TCF1+ stem-like and PD1+ TCF1-TIM3+ effector-like cells. Clustering of TCF1+ CD8 T cells with antigen presenting cells in immune niches was prognostic for local control, even without pSRS. Following pSRS, CD8 T cell and immune niche density were transiently reduced compared to untreated BrM, followed by a rebound 6+ days post pSRS with an increased frequency of TCF1- effector-like cells. In sum, pSRS is safe and therapeutically beneficial, and these data provide a framework for how pSRS may be leveraged to maximize intracranial CD8 T cell responses.

SummaryHarmonizing variant indexing and allele assignments across datasets is crucial for data integrity in cross-dataset studies such as multi-cohort genome-wide association studies, meta-analyses, and the development, validation, and application of polygenic risk scores. Ensuring this indexing and allele consistency is a laborious, time-consuming, and error-prone process requiring a certain degree of computational proficiency. Here, we introduce GRIEVOUS, a command-line tool for cross-dataset variant homogenization. By means of an internal database and a custom indexing methodology, GRIEVOUS identifies, formats, and aligns all biallelic single nucleotide polymorphisms (SNPs) across all summary statistic and genotype files of interest. Upon completion of dataset harmonization, GRIEVOUS can also be used to extract the maximal set of biallelic SNPs common to all datasets.Availability and implementationGRIEVOUS and all supporting documentation and tutorials can be found at https://github.com/jvtalwar/GRIEVOUS. It is freely and publicly available under the MIT license and can be installed via pip.

BackgroundThe US government considers veterans to have been exposed to Agent Orange if they served in Vietnam while the carcinogen was in use, and these veterans are often deemed at high risk of prostate cancer (PCa). Here, we assess whether presumed Agent Orange exposure is independently associated with increased risk of any metastatic or fatal PCa in a diverse Veteran cohort still alive in the modern era (at least 2011), when accounting for race/ethnicity, family history, and genetic risk.Patients and methodsParticipants in the Million Veteran Program (MVP; enrollment began in 2011) who were on active duty during the Vietnam War era (August 1964-April 1975) were included (n = 301,470). Agent Orange exposure was determined using the US government definition. Genetic risk was assessed via a validated polygenic hazard score. Associations with age at diagnosis of any PCa, metastatic PCa, and death from PCa were assessed via Cox proportional hazards models.Results and interpretationOn univariable analysis, exposure to Agent Orange was not associated with increased PCa (hazard ratio [HR]: 1.02, 95% confidence interval [CI]: 1.00-1.04, p = 0.06), metastatic PCa (HR: 0.98, 95% CI: 0.91-1.05, p = 0.55), or fatal PCa (HR: 0.94, 95% CI: 0.79-1.09, p = 0.41). When accounting for race/ethnicity and family history, Agent Orange exposure was independently associated with slightly increased risk of PCa (HR: 1.06, 95% CI: 1.04-1.09,

BackgroundImpaired DNA damage response (DDR) can affect immune checkpoint inhibitors (ICI) efficacy and lead to heightened immune activation. We assessed the impact of pathogenic or likely pathogenic (P/LP) germline DDR mutations on ICI response and toxicity.Materials and methodsA retrospective analysis of 131 cancer patients with germline DNA testing and ICI treatment was performed.ResultsNinety-two patients were DDR-negative (DDR-), and 39 had ≥1 DDR mutation (DDR+). DDR+ patients showed higher objective response rates (ORRs) compared to DDR- in univariate and multivariable analyses, adjusting for age and metastatic disease (62% vs. 23%, unadjusted OR = 5.41; 95% CI, 2.41-12.14; adjusted OR 5.94; 95% CI, 2.35-15.06). Similar results were seen in mismatch repair (MMR), DDR pathways with intact MMR (DDR+MMRi), and homologous recombination (HR) subgroups versus DDR- (adjusted OR MMR = 24.52; 95% CI 2.72-221.38, DDR+MMRi = 4.26; 95% CI, 1.57-11.59, HR = 4.74; 95% CI, 1.49-15.11). DDR+ patients also had higher ORRs with concurrent chemotherapy (82% vs. 39% DDR-, p=0.03) or concurrent tyrosine kinase inhibitors (50% vs. 5% DDR-, p=0.03). No significant differences in immune-related adverse events were observed between DDR+ and DDR- cohorts.ConclusionP/LP germline DDR mutations may enhance ICI response without significant additional toxicity.

Circular extrachromosomal DNA (ecDNA) in patient tumors is an important driver of oncogenic gene expression, evolution of drug resistance and poor patient outcomes. Applying computational methods for the detection and reconstruction of ecDNA across a retrospective cohort of 481 medulloblastoma tumors from 465 patients, we identify circular ecDNA in 82 patients (18%). Patients with ecDNA-positive medulloblastoma were more than twice as likely to relapse and three times as likely to die within 5 years of diagnosis. A subset of tumors harbored multiple ecDNA lineages, each containing distinct amplified oncogenes. Multimodal sequencing, imaging and CRISPR inhibition experiments in medulloblastoma models reveal intratumoral heterogeneity of ecDNA copy number per cell and frequent putative 'enhancer rewiring' events on ecDNA. This study reveals the frequency and diversity of ecDNA in medulloblastoma, stratified into molecular subgroups, and suggests copy number heterogeneity and enhancer rewiring as oncogenic features of ecDNA.

Immune checkpoint blockade (ICB) targeting PD-1 and CTLA-4 has revolutionized cancer treatment. However, many cancers do not respond to ICB, prompting the search for additional strategies to achieve durable responses. G-protein-coupled receptors (GPCRs) are the most intensively studied drug targets but are underexplored in immuno-oncology. Here, we cross-integrated large singe-cell RNA-sequencing datasets from CD8+ T cells covering 19 distinct cancer types and identified an enrichment of Gαs-coupled GPCRs on exhausted CD8+ T cells. These include EP2, EP4, A2AR, β1AR and β2AR, all of which promote T cell dysfunction. We also developed transgenic mice expressing a chemogenetic CD8-restricted Gαs-DREADD to activate CD8-restricted Gαs signaling and show that a Gαs-PKA signaling axis promotes CD8+ T cell dysfunction and immunotherapy failure. These data indicate that Gαs-GPCRs are druggable immune checkpoints that might be targeted to enhance the response to ICB immunotherapies.

Autoimmunity and cancer represent two different aspects of immune dysfunction. Autoimmunity is characterized by breakdowns in immune self-tolerance, while impaired immune surveillance can allow for tumorigenesis. The class I major histocompatibility complex (MHC-I), which displays derivatives of the cellular peptidome for immune surveillance by CD8+ T cells, serves as a common genetic link between these conditions. As melanoma-specific CD8+ T cells have been shown to target melanocyte-specific peptide antigens more often than melanoma-specific antigens, we investigated whether vitiligo- and psoriasis-predisposing MHC-I alleles conferred a melanoma-protective effect. In individuals with cutaneous melanoma from both The Cancer Genome Atlas (n = 451) and an independent validation set (n = 586), MHC-I autoimmune-allele carrier status was significantly associated with a later age of melanoma diagnosis. Furthermore, MHC-I autoimmune-allele carriers were significantly associated with decreased risk of developing melanoma in the Million Veteran Program (OR = 0.962, p = 0.024). Existing melanoma polygenic risk scores (PRSs) did not predict autoimmune-allele carrier status, suggesting these alleles provide orthogonal risk-relevant information. Mechanisms of autoimmune protection were neither associated with improved melanoma-driver mutation association nor improved gene-level conserved antigen presentation relative to common alleles. However, autoimmune alleles showed higher affinity relative to common alleles for particular windows of melanocyte-conserved antigens and loss of heterozygosity of autoimmune alleles caused the greatest reduction in presentation for several conserved antigens across individuals with loss of HLA alleles. Overall, this study presents evidence that MHC-I autoimmune-risk alleles modulate melanoma risk unaccounted for by current PRSs.

BackgroundGenetic scores may provide an objective measure of prostate cancer risk and thus inform screening decisions. We evaluated whether a polygenic hazard score based on 290 genetic variants (PHS290) is associated with prostate cancer risk in a diverse population, including Black men, who have higher average risk of prostate cancer death but are often treated as a homogeneously high-risk group.MethodsThis was a retrospective analysis of the Million Veteran Program, a national, population-based cohort study of US military veterans conducted 2011-2021. Cox proportional hazards analyses tested for association of genetic and other risk factors (including self-reported race and ethnicity and family history) with age at death from prostate cancer, age at diagnosis of metastatic (nodal or distant) prostate cancer, and age at diagnosis of any prostate cancer.ResultsA total of 590 750 male participants were included. Median age at last follow-up was 69 years. PHS290 was associated with fatal prostate cancer in the full cohort and for each racial and ethnic group (P

Ependymoma is a tumor of the brain or spinal cord. The two most common and aggressive molecular groups of ependymoma are the supratentorial ZFTA-fusion associated and the posterior fossa ependymoma group A. In both groups, tumors occur mainly in young children and frequently recur after treatment. Although molecular mechanisms underlying these diseases have recently been uncovered, they remain difficult to target and innovative therapeutic approaches are urgently needed. Here, we use genome-wide chromosome conformation capture (Hi-C), complemented with CTCF and H3K27ac ChIP-seq, as well as gene expression and DNA methylation analysis in primary and relapsed ependymoma tumors, to identify chromosomal conformations and regulatory mechanisms associated with aberrant gene expression. In particular, we observe the formation of new topologically associating domains ('neo-TADs') caused by structural variants, group-specific 3D chromatin loops, and the replacement of CTCF insulators by DNA hyper-methylation. Through inhibition experiments, we validate that genes implicated by these 3D genome conformations are essential for the survival of patient-derived ependymoma models in a group-specific manner. Thus, this study extends our ability to reveal tumor-dependency genes by 3D genome conformations even in tumors that lack targetable genetic alterations.

With the continued promise of immunotherapy for treating cancer, understanding how host genetics contributes to the tumor immune microenvironment (TIME) is essential to tailoring cancer screening and treatment strategies. Here, we study 1084 eQTLs affecting the TIME found through analysis of The Cancer Genome Atlas and literature curation. These TIME eQTLs are enriched in areas of active transcription, and associate with gene expression in specific immune cell subsets, such as macrophages and dendritic cells. Polygenic score models built with TIME eQTLs reproducibly stratify cancer risk, survival and immune checkpoint blockade (ICB) response across independent cohorts. To assess whether an eQTL-informed approach could reveal potential cancer immunotherapy targets, we inhibit CTSS, a gene implicated by cancer risk and ICB response-associated polygenic models; CTSS inhibition results in slowed tumor growth and extended survival in vivo. These results validate the potential of integrating germline variation and TIME characteristics for uncovering potential targets for immunotherapy.

Considerable research has shown that testosterone regulates many physiological systems, modulates clinical disorders, and contributes to health outcome. However, studies on the interaction of testosterone levels with depression and the antidepressant effect of testosterone replacement therapy in hypogonadal men with depression have been inconclusive. Current findings indicate that low circulating levels of total testosterone meeting stringent clinical criteria for hypogonadism and testosterone deficiency induced by androgen deprivation therapy are associated with increased risk for depression and current depressive symptoms. The benefits of testosterone replacement therapy in men with major depressive disorder and low testosterone levels in the clinically defined hypogonadal range remain uncertain and require further investigation. Important considerations going forward are that major depressive disorder is a heterogeneous phenotype with depressed individuals differing in inherited polygenic determinants, onset and clinical course, symptom complexes, and comorbidities that contribute to potential multifactorial differences in pathophysiology. Furthermore, polygenic mechanisms are likely to be critical to the biological heterogeneity that influences testosterone-depression interactions. A genetically informed precision medicine approach using genes regulating testosterone levels and androgen receptor sensitivity will likely be essential in gaining critical insight into the role of testosterone in depression.

BackgroundProstate cancer (PrCa) is one of the most genetically driven solid cancers with heritability estimates as high as 57%. Men of African ancestry are at an increased risk of PrCa; however, current polygenic risk score (PRS) models are based on European ancestry groups and may not be broadly applicable. The objective of this study was to construct an African ancestry-specific PrCa PRS (PRState) and evaluate its performance.MethodsAfrican ancestry group of 4,533 individuals in ELLIPSE consortium was used for discovery of African ancestry-specific PrCa SNPs. PRState was constructed as weighted sum of genotypes and effect sizes from genome-wide association study (GWAS) of PrCa in African ancestry group. Performance was evaluated using ROC-AUC analysis.ResultsWe identified African ancestry-specific PrCa risk loci on chromosomes 3, 8, and 11 and constructed a polygenic risk score (PRS) from 10 African ancestry-specific PrCa risk SNPs, achieving an AUC of 0.61 [0.60-0.63] and 0.65 [0.64-0.67], when combined with age and family history. Performance dropped significantly when using ancestry-mismatched PRS models but remained comparable when using trans-ancestry models. Importantly, we validated the PRState score in the Million Veteran Program (MVP), demonstrating improved prediction of PrCa and metastatic PrCa in individuals of African ancestry.ConclusionsAfrican ancestry-specific PRState improves PrCa prediction in African ancestry groups in ELLIPSE consortium and MVP. This study underscores the need for inclusion of individuals of African ancestry in gene variant discovery to optimize PRSs and identifies African ancestry-specific variants for use in future studies.

PurposeThe contribution of common genetic variants to pre-cancer progression is understudied due to long follow-up time, rarity of poor outcomes and lack of available germline DNA collection. Alternatively, DNA from diagnostic archival tissue is available, but its somatic nature, limited quantity and suboptimal quality would require an accurate cost-effective genome-wide germline genotyping methodology.Experimental designBlood and tissue DNA from 10 individuals were used to benchmark the accuracy of Single Nucleotide Polymorphisms (SNP) genotypes, Polygenic Risk Scores (PRS) or HLA haplotypes using low-coverage whole-genome sequencing (lc-WGS) and genotype imputation. Tissue-derived PRS were further evaluated for 36 breast cancer patients (11.7 years median follow-up time) diagnosed with DCIS and used to model the risk of Breast Cancer Subsequent Events (BCSE).ResultsTissue-derived germline DNA profiling resulted in accurate genotypes at common SNPs (blood correlation r2 > 0.94) and across 22 disease-related polygenic risk scores (PRS, mean correlation r = 0.93). Imputed Class I and II HLA haplotypes were 96.7% and 82.5% concordant with clinical-grade blood HLA haplotypes, respectively. In DCIS patients, tissue-derived PRS was significantly associated with BCSE (HR = 2, 95% CI 1.2-3.8). The top and bottom decile patients had an estimated 28% and 5% chance of BCSE at 10 years, respectively.ConclusionsArchival tissue DNA germline profiling using lc-WGS and imputation, represents a cost and resource-effective alternative in the retrospective design of long-term disease genetic studies. Initial results in breast cancer suggest that common risk variants contribute to pre-cancer progression.

Prostate cancer (PrCa) is one of the most genetically driven solid cancers with heritability estimates as high as 57%. African American men are at an increased risk of PrCa; however, current risk prediction models are based on European ancestry groups and may not be broadly applicable. In this study, we define an African ancestry group of 4,533 individuals to develop an African ancestry-specific PrCa polygenic risk score (PRState). We identified risk loci on chromosomes 3, 8, and 11 in the African ancestry group GWAS and constructed a polygenic risk score (PRS) from 10 African ancestry-specific PrCa risk SNPs, achieving an AUC of 0.61 [0.60-0.63] and 0.65 [0.64-0.67], when combined with age and family history. Performance dropped significantly when using ancestry-mismatched PRS models but remained comparable when using trans-ancestry models. Importantly, we validated the PRState score in the Million Veteran Program, demonstrating improved prediction of PrCa and metastatic PrCa in African American individuals. This study underscores the need for inclusion of individuals of African ancestry in gene variant discovery to optimize PRS.

BackgroundDNA sequencing is increasingly incorporated into the routine care of cancer patients, many of whom also carry inherited, moderate/high-penetrance variants associated with other diseases. Yet, the prevalence and consequence of such variants remain unclear.MethodsWe analyzed the germline genomes of 10,389 adult cancer cases in the TCGA cohort, identifying pathogenic/likely pathogenic variants in autosomal-dominant genes, autosomal-recessive genes, and 59 medically actionable genes curated by the American College of Molecular Genetics (i.e., the ACMG 59 genes). We also analyzed variant- and gene-level expression consequences in carriers.ResultsThe affected genes exhibited varying pan-ancestry and population-specific patterns, and overall, the European population showed the highest frequency of pathogenic/likely pathogenic variants. We further identified genes showing expression consequence supporting variant functionality, including altered gene expression, allelic specific expression, and mis-splicing determined by a massively parallel splicing assay.ConclusionsOur results demonstrate that expression-altering variants are found in a substantial fraction of cases and illustrate the yield of genomic risk assessments for a wide range of diseases across diverse populations.

BackgroundDNA sequencing is increasingly incorporated into the routine care of cancer patients, many of whom also carry inherited, moderate/high-penetrance variants associated with other diseases. Yet, the prevalence and consequence of such variants remain unclear.MethodsWe analyzed the germline genomes of 10,389 adult cancer cases in the TCGA cohort, identifying pathogenic/likely pathogenic variants in autosomal-dominant genes, autosomal-recessive genes, and 59 medically actionable genes curated by the American College of Molecular Genetics (i.e., the ACMG 59 genes). We also analyzed variant- and gene-level expression consequences in carriers.ResultsThe affected genes exhibited varying pan-ancestry and population-specific patterns, and overall, the European population showed the highest frequency of pathogenic/likely pathogenic variants. We further identified genes showing expression consequence supporting variant functionality, including altered gene expression, allelic specific expression, and mis-splicing determined by a massively parallel splicing assay.ConclusionsOur results demonstrate that expression-altering variants are found in a substantial fraction of cases and illustrate the yield of genomic risk assessments for a wide range of diseases across diverse populations.

Insights into oncogenesis derived from cancer susceptibility loci (SNP) hold the potential to facilitate better cancer management and treatment through precision oncology. However, therapeutic insights have thus far been limited by our current lack of understanding regarding both interactions of these loci with somatic cancer driver mutations and their influence on tumorigenesis. For example, although both germline and somatic genetic variation to the p53 tumor suppressor pathway are known to promote tumorigenesis, little is known about the extent to which such variants cooperate to alter pathway activity. Here we hypothesize that cancer risk-associated germline variants interact with somatic TP53 mutational status to modify cancer risk, progression, and response to therapy. Focusing on a cancer risk SNP (rs78378222) with a well-documented ability to directly influence p53 activity as well as integration of germline datasets relating to cancer susceptibility with tumor data capturing somatically-acquired genetic variation provided supportive evidence for this hypothesis. Integration of germline and somatic genetic data enabled identification of a novel entry point for therapeutic manipulation of p53 activities. A cluster of cancer risk SNPs resulted in increased expression of prosurvival p53 target gene KITLG and attenuation of p53-mediated responses to genotoxic therapies, which were reversed by pharmacologic inhibition of the prosurvival c-KIT signal. Together, our results offer evidence of how cancer susceptibility SNPs can interact with cancer driver genes to affect cancer progression and identify novel combinatorial therapies. SIGNIFICANCE: These results offer evidence of how cancer susceptibility SNPs can interact with cancer driver genes to affect cancer progression and present novel therapeutic targets.

BackgroundGenome-wide association studies (GWASs) have enriched the fields of genomics and drug development. Adrenocortical carcinoma (ACC) is a rare cancer with a bimodal age distribution and inadequate treatment options. Paediatric ACC is frequently associated with TP53 mutations, with particularly high incidence in Southern Brazil due to the TP53 p.R337H (R337H) germline mutation. The heterogeneous risk among carriers suggests other genetic modifiers could exist.MethodsWe analysed clinical, genotype and gene expression data derived from paediatric ACC, R337H carriers, and adult ACC patients. We restricted our analyses to single nucleotide polymorphisms (SNPs) previously identified in GWASs to associate with disease or human traits.ResultsA SNP, rs971074, in the alcohol dehydrogenase 7 gene significantly and reproducibly associated with allelic differences in ACC age-of-onset in both cohorts. Patients homozygous for the minor allele were diagnosed up to 16 years earlier. This SNP resides in a gene involved in the retinoic acid (RA) pathway and patients with differing levels of RA pathway gene expression in their tumours associate with differential ACC progression.ConclusionsThese results identify a novel genetic component to ACC development that resides in the retinoic acid pathway, thereby informing strategies to develop management, preventive and therapeutic treatments for ACC.

BackgroundGenome-wide association studies (GWASs) have enriched the fields of genomics and drug development. Adrenocortical carcinoma (ACC) is a rare cancer with a bimodal age distribution and inadequate treatment options. Paediatric ACC is frequently associated with TP53 mutations, with particularly high incidence in Southern Brazil due to the TP53 p.R337H (R337H) germline mutation. The heterogeneous risk among carriers suggests other genetic modifiers could exist.MethodsWe analysed clinical, genotype and gene expression data derived from paediatric ACC, R337H carriers, and adult ACC patients. We restricted our analyses to single nucleotide polymorphisms (SNPs) previously identified in GWASs to associate with disease or human traits.ResultsA SNP, rs971074, in the alcohol dehydrogenase 7 gene significantly and reproducibly associated with allelic differences in ACC age-of-onset in both cohorts. Patients homozygous for the minor allele were diagnosed up to 16 years earlier. This SNP resides in a gene involved in the retinoic acid (RA) pathway and patients with differing levels of RA pathway gene expression in their tumours associate with differential ACC progression.ConclusionsThese results identify a novel genetic component to ACC development that resides in the retinoic acid pathway, thereby informing strategies to develop management, preventive and therapeutic treatments for ACC.

T cells engineered to express a tumor-specific αβ T cell receptor (TCR) mediate anti-tumor immunity. However, mispairing of the therapeutic αβ chains with endogenous αβ chains reduces therapeutic TCR surface expression and generates self-reactive TCRs. We report a general strategy to prevent TCR mispairing: swapping constant domains between the α and β chains of a therapeutic TCR. When paired, domain-swapped (ds)TCRs assemble with CD3, express on the cell surface, and mediate antigen-specific T cell responses. By contrast, dsTCR chains mispaired with endogenous chains cannot properly assemble with CD3 or signal, preventing autoimmunity. We validate this approach in cell-based assays and in a mouse model of TCR gene transfer-induced graft-versus-host disease. We also validate a related approach whereby replacement of αβ TCR domains with corresponding γδ TCR domains yields a functional TCR that does not mispair. This work enables the design of safer TCR gene therapies for cancer immunotherapy.

T cells engineered to express a tumor-specific αβ T cell receptor (TCR) mediate anti-tumor immunity. However, mispairing of the therapeutic αβ chains with endogenous αβ chains reduces therapeutic TCR surface expression and generates self-reactive TCRs. We report a general strategy to prevent TCR mispairing: swapping constant domains between the α and β chains of a therapeutic TCR. When paired, domain-swapped (ds)TCRs assemble with CD3, express on the cell surface, and mediate antigen-specific T cell responses. By contrast, dsTCR chains mispaired with endogenous chains cannot properly assemble with CD3 or signal, preventing autoimmunity. We validate this approach in cell-based assays and in a mouse model of TCR gene transfer-induced graft-versus-host disease. We also validate a related approach whereby replacement of αβ TCR domains with corresponding γδ TCR domains yields a functional TCR that does not mispair. This work enables the design of safer TCR gene therapies for cancer immunotherapy.