Epigenetic patient stratification via contrastive machine learning refines hallmark biomarkers in minoritized children with asthma.

Identifying and refining clinically significant patient stratification is a critical step toward realizing the promise of precision medicine in asthma. Several peripheral blood hallmarks, including total peripheral blood eosinophil count (BEC) and immunoglobulin E (IgE) levels, are routinely used in asthma clinical practice for endotype classification and predicting response to state-of-the-art targeted biologic drugs. However, these biomarkers appear ineffective in predicting treatment outcomes in some patients, and they differ in distribution between racially and ethnically diverse populations, potentially compromising medical care and hindering health equity due to biases in drug eligibility. Here, we propose constructing an unbiased patient stratification score based on DNA methylation (DNAm) and utilizing it to refine the efficacy of hallmark biomarkers for predicting drug response. We developed Phenotype Aware Component Analysis (PACA), a novel contrastive machine-learning method for learning combinations of DNAm sites reflecting biomedically meaningful patient stratifications. Leveraging whole-blood DNAm from Latino (discovery; n=1,016) and African American (replication; n=756) pediatric asthma case-control cohorts, we applied PACA to refine the prediction of bronchodilator response (BDR) to the short-acting ß2-agonist albuterol, the most used drug to treat acute bronchospasm worldwide. While BEC and IgE correlate with BDR in the general patient population, our PACA-derived DNAm score renders these biomarkers predictive of drug response only in patients with high DNAm scores. BEC correlates with BDR in patients with upper-quartile DNAm scores (OR 1.12; 95% CI [1.04, 1.22]; P=7.9 e-4) but not in patients with lower-quartile scores (OR 1.05; 95% CI [0.95, 1.17]; P=0.21); and IgE correlates with BDR in above-median (OR for response 1.42; 95% CI [1.24, 1.63]; P=3.9e-7) but not in below-median patients (OR 1.05; 95% CI [0.92, 1.2]; P=0.57). These results hold within the commonly recognized type 2 (T2)-high asthma endotype but not in T2-low patients, suggesting that our DNAm score primarily represents an unknown variation of T2 asthma. Among T2-high patients with high DNAm scores, elevated BEC or IgE also corresponds to baseline clinical presentation that is known to benefit more from biologic treatment, including higher exacerbation scores, higher allergen sensitization, lower BMI, more recent oral corticosteroids prescription, and lower lung function. Our findings suggest that BEC and IgE, the traditional asthma biomarkers of T2-high asthma, are poor biomarkers for millions worldwide. Revisiting existing drug eligibility criteria relying on these biomarkers in asthma medical care may enhance precision and equity in treatment.

Germline prediction of immune checkpoint inhibitor discontinuation for immune-related adverse events.

Introduction: Immune checkpoint inhibitors (ICIs) can yield remarkable clinical responses in subsets of patients with solid tumors but can also often lead to immune-related adverse events (irAEs). Predictive features of clinically severe irAEs leading to cessation of ICIs have yet to be established. Using data from 1,327 patients with lung cancer treated with ICIs between 2009 and 2022 at four academic medical centers, we evaluated the association of a germline polygenic risk score for autoimmune disease and discontinuation of ICIs due to irAEs. Methods: Using Cox proportional hazards model, we assessed the association between a polygenic risk score for autoimmune disease (PRSAD) and cessation of ICI therapy due to irAEs. All models were adjusted for age at diagnosis, sex, lung cancer histology, type of therapy, recruiting center, and the first 5 principal components. To further understand the differential effects of type of therapy and disease stage on the association between PRSAD and cessation of ICI due to irAEs, we conducted stratified logistic regression analysis by type of ICI therapy and disease stage. Results: We found an association between PRSAD and ICI cessation due to irAEs (HR per SD = 1.18, 95% CI = 1.02 - 1.37, P = 0.03). This association was particularly strong in patients who had ICI cessation due to irAEs within three months of therapy initiation (HR per SD = 1.38, 95% CI = 1.08 - 1.78, P = 0.01). Individuals in the top 20th percentile of PRSAD had 7.2% ICI discontinuation for irAEs by three months, compared to 3.9% discontinuation by three months among patients in the bottom 80th percentile (log-rank P = 0.02). In addition, among patients who received combination PD-1/PD-L1 and CTLA-4 inhibitor therapy, PRSAD had an OR per SD of 1.86 (95% CI = 1.08 - 3.51, P = 0.04). Conclusions: We demonstrate an association between a polygenic risk score for autoimmune disease and early ICI discontinuation for irAEs, particularly among patients treated with combination ICI therapy. Our results suggest that germline genetics may be used as an adjunctive tool for risk stratification around ICI clinical decision-making in solid tumor oncology.

Association of ESR1 Germline Variants with TP53 Somatic Variants in Breast Tumors in a Genome-wide Study.

In breast tumors, somatic mutation frequencies in TP53 and PIK3CA vary by tumor subtype and ancestry. Emerging data suggest tumor mutation status is associated with germline variants and genetic ancestry. We aimed to identify germline variants that are associated with somatic TP53 or PIK3CA mutation status in breast tumors. A genome-wide association study was conducted in 2,850 women of European ancestry with breast cancer using TP53 and PIK3CA mutation status (positive or negative) as well as specific functional categories [e.g., TP53 gain-of-function (GOF) and loss-of-function, PIK3CA activating] as phenotypes. Germline variants showing evidence of association were selected for validation analyses and tested in multiple independent datasets. Discovery association analyses found five variants associated with TP53 mutation status with P values <1 × 10-6 and 33 variants with P values <1 × 10-5. Forty-four variants were associated with PIK3CA mutation status with P values <1 × 10-5. In validation analyses, only variants at the ESR1 locus were associated with TP53 mutation status after multiple comparisons corrections. Combined analyses in European and Malaysian populations found ESR1 locus variants rs9383938 and rs9479090 associated with the presence of TP53 mutations overall (P values 2 × 10-11 and 4.6 × 10-10, respectively). rs9383938 also showed association with TP53 GOF mutations (P value 6.1 × 10-7). rs9479090 showed suggestive evidence (P value 0.02) for association with TP53 mutation status in African ancestry populations. No other variants were significantly associated with TP53 or PIK3CA mutation status. Larger studies are needed to confirm these findings and determine if additional variants contribute to ancestry-specific differences in mutation frequency. SIGNIFICANCE: Emerging data show ancestry-specific differences in TP53 and PIK3CA mutation frequency in breast tumors suggesting that germline variants may influence somatic mutational processes. This study identified variants near ESR1 associated with TP53 mutation status and identified additional loci with suggestive association which may provide biological insight into observed differences.

The Puerto Rican Infant Metagenomic and Epidemiologic Study of Respiratory Outcomes (PRIMERO): Design and Baseline Characteristics for a Birth Cohort Study of Early-life Viral Respiratory Illnesses and Airway Dysfunction in Puerto Rican Children.

Epidemiologic studies demonstrate an association between early-life respiratory illnesses (RIs) and the development of childhood asthma. However, it remains uncertain whether these children are predisposed to both conditions or if early-life RIs induce alterations in airway function, immune responses, or other human biology that contribute to the development of asthma. Puerto Rican children experience a disproportionate burden of early-life RIs and asthma, making them an important population for investigating this complex interplay. PRIMERO, the Puerto Rican Infant Metagenomics and Epidemiologic Study of Respiratory Outcomes , recruited pregnant women and their newborns to investigate how the airways develop in early life among infants exposed to different viral RIs, and will thus provide a critical understanding of childhood asthma development. As the first asthma birth cohort in Puerto Rico, PRIMERO will prospectively follow 2,100 term healthy infants. Collected samples include post-term maternal peripheral blood, infant cord blood, the child's peripheral blood at the year two visit, and the child's nasal airway epithelium, collected using minimally invasive nasal swabs, at birth, during RIs over the first two years of life, and at annual healthy visits until age five. Herein, we describe the study's design, population, recruitment strategy, study visits and procedures, and primary outcomes.

Polygenic risk score for ulcerative colitis predicts immune checkpoint inhibitor-mediated colitis.

Immune checkpoint inhibitor-mediated colitis (IMC) is a common adverse event of treatment with immune checkpoint inhibitors (ICI). We hypothesize that genetic susceptibility to Crohn's disease (CD) and ulcerative colitis (UC) predisposes to IMC. In this study, we first develop a polygenic risk scores for CD (PRSCD) and UC (PRSUC) in cancer-free individuals and then test these PRSs on IMC in a cohort of 1316 patients with ICI-treated non-small cell lung cancer and perform a replication in 873 ICI-treated pan-cancer patients. In a meta-analysis, the PRSUC predicts all-grade IMC (ORmeta=1.35 per standard deviation [SD], 95% CI = 1.12-1.64, P = 2×10-03) and severe IMC (ORmeta=1.49 per SD, 95% CI = 1.18-1.88, P = 9×10-04). PRSCD is not associated with IMC. Furthermore, PRSUC predicts severe IMC among patients treated with combination ICIs (ORmeta=2.20 per SD, 95% CI = 1.07-4.53, P = 0.03). Overall, PRSUC can identify patients receiving ICI at risk of developing IMC and may be useful to monitor patients and improve patient outcomes.

A noncoding regulatory variant in IKZF1 increases acute lymphoblastic leukemia risk in Hispanic/Latino children.

Hispanic/Latino children have the highest risk of acute lymphoblastic leukemia (ALL) in the US compared to other racial/ethnic groups, yet the basis of this remains incompletely understood. Through genetic fine-mapping analyses, we identified a new independent childhood ALL risk signal near IKZF1 in self-reported Hispanic/Latino individuals, but not in non-Hispanic White individuals, with an effect size of ∼1.44 (95% confidence interval = 1.33-1.55) and a risk allele frequency of ∼18% in Hispanic/Latino populations and <0.5% in European populations. This risk allele was positively associated with Indigenous American ancestry, showed evidence of selection in human history, and was associated with reduced IKZF1 expression. We identified a putative causal variant in a downstream enhancer that is most active in pro-B cells and interacts with the IKZF1 promoter. This variant disrupts IKZF1 autoregulation at this enhancer and results in reduced enhancer activity in B cell progenitors. Our study reveals a genetic basis for the increased ALL risk in Hispanic/Latino children.

A comprehensive analysis of clinical and polygenic germline influences on somatic mutational burden.

Tumor mutational burden (TMB), the total number of somatic mutations in the tumor, and copy number burden (CNB), the corresponding measure of aneuploidy, are established fundamental somatic features and emerging biomarkers for immunotherapy. However, the genetic and non-genetic influences on TMB/CNB and, critically, the manner by which they influence patient outcomes remain poorly understood. Here, we present a large germline-somatic study of TMB/CNB with >23,000 individuals across 17 cancer types, of which 12,000 also have extensive clinical, treatment, and overall survival (OS) measurements available. We report dozens of clinical associations with TMB/CNB, observing older age and male sex to have a strong effect on TMB and weaker impact on CNB. We additionally identified significant germline influences on TMB/CNB, including fine-scale European ancestry and germline polygenic risk scores (PRSs) for smoking, tanning, white blood cell counts, and educational attainment. We quantify the causal effect of exposures on somatic mutational processes using Mendelian randomization. Many of the identified features associated with TMB/CNB were additionally associated with OS for individuals treated at a single tertiary cancer center. For individuals receiving immunotherapy, we observed a complex relationship between PRSs for educational attainment, self-reported college attainment, TMB, and survival, suggesting that the influence of this biomarker may be substantially modified by socioeconomic status. While the accumulation of somatic alterations is a stochastic process, our work demonstrates that it can be shaped by host characteristics including germline genetics.

High-throughput PRIME-editing screens identify functional DNA variants in the human genome.

Despite tremendous progress in detecting DNA variants associated with human disease, interpreting their functional impact in a high-throughput and single-base resolution manner remains challenging. Here, we develop a pooled prime-editing screen method, PRIME, that can be applied to characterize thousands of coding and non-coding variants in a single experiment with high reproducibility. To showcase its applications, we first identified essential nucleotides for a 716 bp MYC enhancer via PRIME-mediated single-base resolution analysis. Next, we applied PRIME to functionally characterize 1,304 genome-wide association study (GWAS)-identified non-coding variants associated with breast cancer and 3,699 variants from ClinVar. We discovered that 103 non-coding variants and 156 variants of uncertain significance are functional via affecting cell fitness. Collectively, we demonstrate that PRIME is capable of characterizing genetic variants at single-base resolution and scale, advancing accurate genome annotation for disease risk prediction, diagnosis, and therapeutic target identification.

Association of ESR1 germline variants with TP53 somatic variants in breast tumors in a genome-wide study.

Background: In breast tumors, somatic mutation frequencies in TP53 and PIK3CA vary by tumor subtype and ancestry. HER2 positive and triple negative breast cancers (TNBC) have a higher frequency of TP53 somatic mutations than other subtypes. PIK3CA mutations are more frequently observed in hormone receptor positive tumors. Emerging data suggest tumor mutation status is associated with germline variants and genetic ancestry. We aimed to identify germline variants that are associated with somatic TP53 or PIK3CA mutation status in breast tumors. Methods: A genome-wide association study was conducted using breast cancer mutation status of TP53 and PIK3CA and functional mutation categories including TP53 gain of function (GOF) and loss of function mutations and PIK3CA activating/hotspot mutations. The discovery analysis consisted of 2850 European ancestry women from three datasets. Germline variants showing evidence of association with somatic mutations were selected for validation analyses based on predicted function, allele frequency, and proximity to known cancer genes or risk loci. Candidate variants were assessed for association with mutation status in a multi-ancestry validation study, a Malaysian study, and a study of African American/Black women with TNBC. Results: The discovery Germline x Mutation (GxM) association study found five variants associated with one or more TP53 phenotypes with P values <1×10-6, 33 variants associated with one or more TP53 phenotypes with P values <1×10-5, and 44 variants associated with one or more PIK3CA phenotypes with P values <1×10-5. In the multi-ancestry and Malaysian validation studies, germline ESR1 locus variant, rs9383938, was associated with the presence of TP53 mutations overall (P values 6.8×10-5 and 9.8×10-8, respectively) and TP53 GOF mutations (P value 8.4×10-6). Multiple variants showed suggestive evidence of association with PIK3CA mutation status in the validation studies, but none were significant after correction for multiple comparisons. Conclusions: We found evidence that germline variants were associated with TP53 and PIK3CA mutation status in breast cancers. Variants near the estrogen receptor alpha gene, ESR1, were significantly associated with overall TP53 mutations and GOF mutations. Larger multi-ancestry studies are needed to confirm these findings and determine if these variants contribute to ancestry-specific differences in mutation frequency.

Novel insights into the whole-blood DNA methylome of asthma in ethnically diverse children and youth.

BACKGROUND: The epigenetic mechanisms of asthma remain largely understudied in African Americans and Hispanics/Latinos, two populations disproportionately affected by asthma. We aimed to identify markers, regions and processes with differential patterns of DNA methylation (DNAm) in whole blood by asthma status in ethnically diverse children and youth, and to assess their functional consequences. METHODS: DNAm levels were profiled with the Infinium MethylationEPIC or HumanMethylation450 BeadChip arrays among 1226 African Americans or Hispanics/Latinos and assessed for differential methylation per asthma status at the CpG and region (differentially methylated region (DMR)) level. Novel associations were validated in blood and/or nasal epithelium from ethnically diverse children and youth. The functional and biological implications of the markers identified were investigated by combining epigenomics with transcriptomics from study participants. RESULTS: 128 CpGs and 196 DMRs were differentially methylated after multiple testing corrections, including 92.3% and 92.8% novel associations, respectively. 41 CpGs were replicated in other Hispanics/Latinos, prioritising cg17647904 (NCOR2) and cg16412914 (AXIN1) as asthma DNAm markers. Significant DNAm markers were enriched in previous associations for asthma, fractional exhaled nitric oxide, bacterial infections, immune regulation or eosinophilia. Functional annotation highlighted epigenetically regulated gene networks involved in corticosteroid response, host defence and immune regulation. Several implicated genes are targets for approved or experimental drugs, including TNNC1 and NDUFA12. Many differentially methylated loci previously associated with asthma were validated in our study. CONCLUSIONS: We report novel whole-blood DNAm markers for asthma underlying key processes of the disease pathophysiology and confirm the transferability of previous asthma DNAm associations to ethnically diverse populations.

Epigenomic response to albuterol treatment in asthma-relevant airway epithelial cells.

BACKGROUND: Albuterol is the first-line asthma medication used in diverse populations. Although DNA methylation (DNAm) is an epigenetic mechanism involved in asthma and bronchodilator drug response (BDR), no study has assessed whether albuterol could induce changes in the airway epithelial methylome. We aimed to characterize albuterol-induced DNAm changes in airway epithelial cells, and assess potential functional consequences and the influence of genetic variation and asthma-related clinical variables. RESULTS: We followed a discovery and validation study design to characterize albuterol-induced DNAm changes in paired airway epithelial cultures stimulated in vitro with albuterol. In the discovery phase, an epigenome-wide association study using paired nasal epithelial cultures from Puerto Rican children (n = 97) identified 22 CpGs genome-wide associated with repeated-use albuterol treatment (p < 9 × 10-8). Albuterol predominantly induced a hypomethylation effect on CpGs captured by the EPIC array across the genome (probability of hypomethylation: 76%, p value = 3.3 × 10-5). DNAm changes on the CpGs cg23032799 (CREB3L1), cg00483640 (MYLK4-LINC01600), and cg05673431 (KSR1) were validated in nasal epithelia from 10 independent donors (false discovery rate [FDR] < 0.05). The effect on the CpG cg23032799 (CREB3L1) was cross-tissue validated in bronchial epithelial cells at nominal level (p = 0.030). DNAm changes in these three CpGs were shown to be influenced by three independent genetic variants (FDR < 0.05). In silico analyses showed these polymorphisms regulated gene expression of nearby genes in lungs and/or fibroblasts including KSR1 and LINC01600 (6.30 × 10-14 ≤ p ≤ 6.60 × 10-5). Additionally, hypomethylation at the CpGs cg10290200 (FLNC) and cg05673431 (KSR1) was associated with increased gene expression of the genes where they are located (FDR < 0.05). Furthermore, while the epigenetic effect of albuterol was independent of the asthma status, severity, and use of medication, BDR was nominally associated with the effect on the CpG cg23032799 (CREB3L1) (p = 0.004). Gene-set enrichment analyses revealed that epigenomic modifications of albuterol could participate in asthma-relevant processes (e.g., IL-2, TNF-α, and NF-κB signaling pathways). Finally, nine differentially methylated regions were associated with albuterol treatment, including CREB3L1, MYLK4, and KSR1 (adjusted p value < 0.05). CONCLUSIONS: This study revealed evidence of epigenetic modifications induced by albuterol in the mucociliary airway epithelium. The epigenomic response induced by albuterol might have potential clinical implications by affecting biological pathways relevant to asthma.

Promotores' perspectives on the virtual adaptation of a hereditary breast cancer education program.

Breast cancer is the most common cancer in women in the United States (U.S.) and the leading cause of cancer related death among U.S. Hispanics/Latinas (H/Ls). H/Ls have lower rates of screening and longer time to follow up after an abnormal mammogram. We developed a comprehensive community health educator (promotores)-led education and risk identification program for Spanish-speaking H/Ls in California to increase mammography screening, genetic testing, and the understanding of the impact of family history on cancer risk. Due to COVID-19, we adapted the program to a virtual platform. The experience of transforming the program to a virtual platform provided unique opportunities for collaboration between researchers, community partners, and participants. Promotores are major partners in community based participatory research and in the provision of health care services, but their voices are often excluded from scientific reports. This commentary is an effort to provide a platform for promotores to share their experiences and for the readers to understand their approach in bridging the gap between health care services and communities.

High throughput PRIME editing screens identify functional DNA variants in the human genome.

Despite tremendous progress in detecting DNA variants associated with human disease, interpreting their functional impact in a high-throughput and base-pair resolution manner remains challenging. Here, we develop a novel pooled prime editing screen method, PRIME, which can be applied to characterize thousands of coding and non-coding variants in a single experiment with high reproducibility. To showcase its applications, we first identified essential nucleotides for a 716 bp MYC enhancer via PRIME-mediated saturation mutagenesis. Next, we applied PRIME to functionally characterize 1,304 non-coding variants associated with breast cancer and 3,699 variants from ClinVar. We discovered that 103 non-coding variants and 156 variants of uncertain significance are functional via affecting cell fitness. Collectively, we demonstrate PRIME capable of characterizing genetic variants at base-pair resolution and scale, advancing accurate genome annotation for disease risk prediction, diagnosis, and therapeutic target identification.

Polygenic risk score for ulcerative colitis predicts immune checkpoint inhibitor-mediated colitis.

Immune checkpoint inhibitors (ICIs) are a remarkable advancement in cancer therapeutics; however, a substantial proportion of patients develop severe immune-related adverse events (irAEs). Understanding and predicting irAEs is a key to advancing precision immuno-oncology. Immune checkpoint inhibitor-mediated colitis (IMC) is a significant complication from ICI and can have life-threatening consequences. Based on clinical presentation, IMC mimics inflammatory bowel disease, however the link is poorly understood. We hypothesized that genetic susceptibility to Crohn's disease (CD) and ulcerative colitis (UC) may predispose to IMC. We developed and validated polygenic risk scores for CD (PRSCD) and UC (PRSUC) in cancer-free individuals and assessed the role of each of these PRSs on IMC in a cohort of 1,316 patients with non-small cell lung cancer who received ICIs. Prevalence of all-grade IMC in our cohort was 4% (55 cases), and for severe IMC, 2.5% (32 cases). The PRSUC predicted the development of all-grade IMC (HR=1.34 per standard deviation [SD], 95% CI=1.02-1.76, P=0.04) and severe IMC (HR=1.62 per SD, 95% CI=1.12-2.35, P=0.01). PRSCD was not associated with IMC or severe IMC. The association between PRSUC and IMC (all-grade and severe) was consistent in an independent pan-cancer cohort of patients treated with ICIs. Furthermore, PRSUC predicted severe IMC among patients treated with combination ICIs (OR = 2.20 per SD, 95% CI = 1.07-4.53, P=0.03). This is the first study to demonstrate the potential clinical utility of a PRS for ulcerative colitis in identifying patients receiving ICI at high risk of developing IMC, where risk reduction and close monitoring strategies could help improve overall patient outcomes.

Gene expression in African Americans, Puerto Ricans and Mexican Americans reveals ancestry-specific patterns of genetic architecture.

We explored ancestry-related differences in the genetic architecture of whole-blood gene expression using whole-genome and RNA sequencing data from 2,733 African Americans, Puerto Ricans and Mexican Americans. We found that heritability of gene expression significantly increased with greater proportions of African genetic ancestry and decreased with higher proportions of Indigenous American ancestry, reflecting the relationship between heterozygosity and genetic variance. Among heritable protein-coding genes, the prevalence of ancestry-specific expression quantitative trait loci (anc-eQTLs) was 30% in African ancestry and 8% for Indigenous American ancestry segments. Most anc-eQTLs (89%) were driven by population differences in allele frequency. Transcriptome-wide association analyses of multi-ancestry summary statistics for 28 traits identified 79% more gene-trait associations using transcriptome prediction models trained in our admixed population than models trained using data from the Genotype-Tissue Expression project. Our study highlights the importance of measuring gene expression across large and ancestrally diverse populations for enabling new discoveries and reducing disparities.

Development and testing of a polygenic risk score for breast cancer aggressiveness.

Aggressive breast cancers portend a poor prognosis, but current polygenic risk scores (PRSs) for breast cancer do not reliably predict aggressive cancers. Aggressiveness can be effectively recapitulated using tumor gene expression profiling. Thus, we sought to develop a PRS for the risk of recurrence score weighted on proliferation (ROR-P), an established prognostic signature. Using 2363 breast cancers with tumor gene expression data and single nucleotide polymorphism (SNP) genotypes, we examined the associations between ROR-P and known breast cancer susceptibility SNPs using linear regression models. We constructed PRSs based on varying p-value thresholds and selected the optimal PRS based on model r2 in 5-fold cross-validation. We then used Cox proportional hazards regression to test the ROR-P PRS's association with breast cancer-specific survival in two independent cohorts totaling 10,196 breast cancers and 785 events. In meta-analysis of these cohorts, higher ROR-P PRS was associated with worse survival, HR per SD = 1.13 (95% CI 1.06-1.21, p = 4.0 × 10-4). The ROR-P PRS had a similar magnitude of effect on survival as a comparator PRS for estrogen receptor (ER)-negative versus positive cancer risk (PRSER-/ER+). Furthermore, its effect was minimally attenuated when adjusted for PRSER-/ER+, suggesting that the ROR-P PRS provides additional prognostic information beyond ER status. In summary, we used integrated analysis of germline SNP and tumor gene expression data to construct a PRS associated with aggressive tumor biology and worse survival. These findings could potentially enhance risk stratification for breast cancer screening and prevention.

An integrated tumor, immune and microbiome atlas of colon cancer.

The lack of multi-omics cancer datasets with extensive follow-up information hinders the identification of accurate biomarkers of clinical outcome. In this cohort study, we performed comprehensive genomic analyses on fresh-frozen samples from 348 patients affected by primary colon cancer, encompassing RNA, whole-exome, deep T cell receptor and 16S bacterial rRNA gene sequencing on tumor and matched healthy colon tissue, complemented with tumor whole-genome sequencing for further microbiome characterization. A type 1 helper T cell, cytotoxic, gene expression signature, called Immunologic Constant of Rejection, captured the presence of clonally expanded, tumor-enriched T cell clones and outperformed conventional prognostic molecular biomarkers, such as the consensus molecular subtype and the microsatellite instability classifications. Quantification of genetic immunoediting, defined as a lower number of neoantigens than expected, further refined its prognostic value. We identified a microbiome signature, driven by Ruminococcus bromii, associated with a favorable outcome. By combining microbiome signature and Immunologic Constant of Rejection, we developed and validated a composite score (mICRoScore), which identifies a group of patients with excellent survival probability. The publicly available multi-omics dataset provides a resource for better understanding colon cancer biology that could facilitate the discovery of personalized therapeutic approaches.

Profiling the Somatic Mutational Landscape of Breast Tumors from Hispanic/Latina Women Reveals Conserved and Unique Characteristics.

Somatic mutational profiling is increasingly being used to identify potential targets for breast cancer. However, limited tumor-sequencing data from Hispanic/Latinas (H/L) are available to guide treatment. To address this gap, we performed whole-exome sequencing (WES) and RNA sequencing on 146 tumors and WES of matched germline DNA from 140 H/L women in California. Tumor intrinsic subtype, somatic mutations, copy-number alterations, and expression profiles of the tumors were characterized and compared with data from tumors of non-Hispanic White (White) women in The Cancer Genome Atlas (TCGA). Eight genes were significantly mutated in the H/L tumors including PIK3CA, TP53, GATA3, MAP3K1, CDH1, CBFB, PTEN, and RUNX1; the prevalence of mutations in these genes was similar to that observed in White women in TCGA. Four previously reported Catalogue of Somatic Mutations in Cancer (COSMIC) mutation signatures (1, 2, 3, 13) were found in the H/L dataset, along with signature 16 that has not been previously reported in other breast cancer datasets. Recurrent amplifications were observed in breast cancer drivers including MYC, FGFR1, CCND1, and ERBB2, as well as a recurrent amplification in 17q11.2 associated with high KIAA0100 gene expression that has been implicated in breast cancer aggressiveness. In conclusion, this study identified a higher prevalence of COSMIC signature 16 and a recurrent copy-number amplification affecting expression of KIAA0100 in breast tumors from H/L compared with White women. These results highlight the necessity of studying underrepresented populations. SIGNIFICANCE: Comprehensive characterization of genomic and transcriptomic alterations in breast tumors from Hispanic/Latina patients reveals distinct genetic alterations and signatures, demonstrating the importance of inclusive studies to ensure equitable care for patients. See related commentary by Schmit et al., p. 2443.

Spectrum of Clinical Presentations, Imaging Findings, and HLA Types in Immune Checkpoint Inhibitor-Induced Hypophysitis.

Context: Hypophysitis is a known immune-related adverse event (irAE) of immune checkpoint inhibitors (CPIs), commonly associated with CTLA-4 inhibitors and less often with PD-1/PD-L1 inhibitors. Objective: We aimed to determine clinical, imaging, and HLA characteristics of CPI-induced hypophysitis (CPI-hypophysitis). Methods: We examined the clinical and biochemical characteristics, magnetic resonance imaging (MRI) of the pituitary, and association with HLA type in patients with CPI-hypophysitis. Results: Forty-nine patients were identified. Mean age was 61.3 years, 61.2% were men, 81.6% were Caucasian, 38.8% had melanoma, and 44.5% received PD-1/PD-L1 inhibitor monotherapy while the remainder received CTLA-4 inhibitor monotherapy or CTLA-4/PD-1 inhibitor combination therapy. A comparison of CTLA-4 inhibitor exposure vs PD-1/PD-L1 inhibitor monotherapy revealed faster time to CPI-hypophysitis (median 84 vs 185 days, P < .01) and abnormal pituitary appearance on MRI (odds ratio 7.00, P = .03). We observed effect modification by sex in the association between CPI type and time to CPI-hypophysitis. In particular, anti-CTLA-4 exposed men had a shorter time to onset than women. MRI changes of the pituitary were most common at the time of hypophysitis diagnosis (55.6% enlarged, 37.0% normal, 7.4% empty or partially empty) but persisted in follow-up (23.8% enlarged, 57.1% normal, 19.1% empty or partially empty). HLA typing was done on 55 subjects; HLA type DQ0602 was over-represented in CPI-hypophysitis relative to the Caucasian American population (39.4% vs 21.5%, P = 0.01) and CPI population. Conclusion: The association of CPI-hypophysitis with HLA DQ0602 suggests a genetic risk for its development. The clinical phenotype of hypophysitis appears heterogenous, with differences in timing of onset, changes in thyroid function tests, MRI changes, and possibly sex related to CPI type. These factors may play an important role in our mechanistic understanding of CPI-hypophysitis.