Ultrasensitive plasma-based monitoring of tumor burden using machine-learning-guided signal enrichment.

In solid tumor oncology, circulating tumor DNA (ctDNA) is poised to transform care through accurate assessment of minimal residual disease (MRD) and therapeutic response monitoring. To overcome the sparsity of ctDNA fragments in low tumor fraction (TF) settings and increase MRD sensitivity, we previously leveraged genome-wide mutational integration through plasma whole-genome sequencing (WGS). Here we now introduce MRD-EDGE, a machine-learning-guided WGS ctDNA single-nucleotide variant (SNV) and copy-number variant (CNV) detection platform designed to increase signal enrichment. MRD-EDGESNV uses deep learning and a ctDNA-specific feature space to increase SNV signal-to-noise enrichment in WGS by ~300× compared to previous WGS error suppression. MRD-EDGECNV also reduces the degree of aneuploidy needed for ultrasensitive CNV detection through WGS from 1 Gb to 200 Mb, vastly expanding its applicability within solid tumors. We harness the improved performance to identify MRD following surgery in multiple cancer types, track changes in TF in response to neoadjuvant immunotherapy in lung cancer and demonstrate ctDNA shedding in precancerous colorectal adenomas. Finally, the radical signal-to-noise enrichment in MRD-EDGESNV enables plasma-only (non-tumor-informed) disease monitoring in advanced melanoma and lung cancer, yielding clinically informative TF monitoring for patients on immune-checkpoint inhibition.

Motive and Opportunity: MYC rearrangements in high-grade B-cell lymphoma with MYC and BCL2 rearrangements-an LLMPP study.

Rearrangements that place the oncogenes MYC, BCL2, or BCL6 adjacent to superenhancers are common in mature B-cell lymphomas. Lymphomas with diffuse large B-cell lymphoma (DLBCL) or high-grade morphology with both MYC and BCL2 rearrangements are classified as high-grade B-cell lymphoma with MYC and BCL2 rearrangements ("double hit": HGBCL-DH-BCL2) and are associated with aggressive disease and poor outcomes. Although it is established that MYC rearrangements involving immunoglobulin (IG) loci are associated with inferior outcomes relative to those involving other non-IG superenhancers, the frequency of, and mechanisms driving, IG vs non-IG MYC rearrangements have not been elucidated. Here we used custom targeted capture and/or whole genome sequencing to characterize oncogene rearrangements across 883 mature B-cell lymphomas including Burkitt lymphoma, follicular lymphoma, DLBCL, and HGBCL-DH-BCL2 tumors. We demonstrate that, while BCL2 rearrangement topology is consistent across entities, HGBCL-DH-BCL2 have distinct MYC rearrangement architecture relative to tumors with single MYC rearrangements or with both MYC and BCL6 rearrangements (HGBCL-DH-BCL6), including both a higher frequency of non-IG rearrangements and different architecture of MYC::IGH rearrangements. The distinct MYC rearrangement patterns in HGBCL-DH-BCL2 occur on the background of high levels of somatic hypermutation across MYC partner loci in HGBCL-DH-BCL2, creating more opportunity to form these rearrangements. Furthermore, because one IGH allele is already disrupted by the existing BCL2 rearrangement, the MYC rearrangement architecture in HGBCL-DH-BCL2 likely reflects selective pressure to preserve both BCL2 and B cell receptor expression. These data provide new mechanistic explanations for the distinct patterns of MYC rearrangements observed across different lymphoma entities.

Whole genome and transcriptome integrated analyses guide clinical care of pediatric poor prognosis cancers.

The role for routine whole genome and transcriptome analysis (WGTA) for poor prognosis pediatric cancers remains undetermined. Here, we characterize somatic mutations, structural rearrangements, copy number variants, gene expression, immuno-profiles and germline cancer predisposition variants in children and adolescents with relapsed, refractory or poor prognosis malignancies who underwent somatic WGTA and matched germline sequencing. Seventy-nine participants with a median age at enrollment of 8.8 y (range 6 months to 21.2 y) are included. Germline pathogenic/likely pathogenic variants are identified in 12% of participants, of which 60% were not known prior. Therapeutically actionable variants are identified by targeted gene report and whole genome in 32% and 62% of participants, respectively, and increase to 96% after integrating transcriptome analyses. Thirty-two molecularly informed therapies are pursued in 28 participants with 54% achieving a clinical benefit rate; objective response or stable disease ≥6 months. Integrated WGTA identifies therapeutically actionable variants in almost all tumors and are directly translatable to clinical care of children with poor prognosis cancers.

Exploration of Germline Correlates and Risk of Immune-Related Adverse Events in Advanced Cancer Patients Treated with Immune Checkpoint Inhibitors.

Immune checkpoint inhibitors (ICIs) are increasingly used in the treatment of many tumor types, and durable responses can be observed in select populations. However, patients may exhibit significant immune-related adverse events (irAEs) that may lead to morbidity. There is limited information on whether the presence of specific germline mutations may highlight those at elevated risk of irAEs. We evaluated 117 patients with metastatic solid tumors or hematologic malignancies who underwent genomic analysis through the ongoing Personalized OncoGenomics (POG) program at BC Cancer and received an ICI during their treatment history. Charts were reviewed for irAEs. Whole genome sequencing of a fresh biopsy and matched normal specimens (blood) was performed at the time of POG enrollment. Notably, we found that MHC class I alleles in the HLA-B27 family, which have been previously associated with autoimmune conditions, were associated with grade 3 hepatitis and pneumonitis (q = 0.007) in patients treated with combination PD-1/PD-L1 and CTLA-4 inhibitors, and PD-1 inhibitors in combination with IDO-1 inhibitors. These data highlight that some patients may have a genetic predisposition to developing irAEs.

SETD1A variant-associated psychosis: A systematic review of the clinical literature and description of two new cases.

OBJECTIVE: SETD1A encodes a histone methyltransferase involved in various cell cycle regulatory processes. Loss-of-function SETD1A variants have been associated with numerous neurodevelopmental phenotypes, including intellectual disability and schizophrenia. While the association between rare coding variants in SETD1A and schizophrenia has achieved genome-wide significance by rare variant burden testing, only a few studies have described the psychiatric phenomenology of such individuals in detail. This systematic review and case report aims to characterize the neurodevelopmental and psychiatric phenotypes of SETD1A variant-associated schizophrenia. METHODS: A PubMed search was completed in July 2022 and updated in May 2023. Only studies that reported individuals with a SETD1A variant as well as a primary psychotic disorder were ultimately included. Additionally, another two previously unpublished cases of SETD1A variant-associated psychosis from our own sequencing cohort are described. RESULTS: The search yielded 32 articles. While 15 articles met inclusion criteria, only five provided case descriptions. In total, phenotypic information was available for 11 individuals, in addition to our own two unpublished cases. Our findings suggest that although individuals with SETD1A variant-associated schizophrenia may share a number of common features, phenotypic variability nonetheless exists. Moreover, although such individuals may exhibit numerous other neurodevelopmental features suggestive of the syndrome, their psychiatric presentations appear to be similar to those of general schizophrenia populations. CONCLUSIONS: Loss-of-function SETD1A variants may underlie the development of psychosis in a small percentage of individuals with schizophrenia. Identifying such individuals may become increasingly important, given the potential for advances in precision medicine treatment approaches.

Genomic structures and regulation patterns at HPV integration sites in cervical cancer.

Human papillomavirus (HPV) integration has been implicated in transforming HPV infection into cancer, but its genomic consequences have been difficult to study using short-read technologies. To resolve the dysregulation associated with HPV integration, we performed long-read sequencing on 63 cervical cancer genomes. We identified six categories of integration events based on HPV-human genomic structures. Of all HPV integrants, defined as two HPV-human breakpoints bridged by an HPV sequence, 24% contained variable copies of HPV between the breakpoints, a phenomenon we termed heterologous integration. Analysis of DNA methylation within and in proximity to the HPV genome at individual integration events revealed relationships between methylation status of the integrant and its orientation and structure. Dysregulation of the human epigenome and neighboring gene expression in cis with the HPV-integrated allele was observed over megabase-ranges of the genome. By elucidating the structural, epigenetic, and allele-specific impacts of HPV integration, we provide insight into the role of integrated HPV in cervical cancer.

Assembly and annotation of the black spruce genome provide insights on spruce phylogeny and evolution of stress response.

Black spruce (Picea mariana [Mill.] B.S.P.) is a dominant conifer species in the North American boreal forest that plays important ecological and economic roles. Here, we present the first genome assembly of P. mariana with a reconstructed genome size of 18.3 Gbp and NG50 scaffold length of 36.0 kbp. A total of 66,332 protein-coding sequences were predicted in silico and annotated based on sequence homology. We analyzed the evolutionary relationships between P. mariana and 5 other spruces for which complete nuclear and organelle genome sequences were available. The phylogenetic tree estimated from mitochondrial genome sequences agrees with biogeography; specifically, P. mariana was strongly supported as a sister lineage to P. glauca and 3 other taxa found in western North America, followed by the European Picea abies. We obtained mixed topologies with weaker statistical support in phylogenetic trees estimated from nuclear and chloroplast genome sequences, indicative of ancient reticulate evolution affecting these 2 genomes. Clustering of protein-coding sequences from the 6 Picea taxa and 2 Pinus species resulted in 34,776 orthogroups, 560 of which appeared to be specific to P. mariana. Analysis of these specific orthogroups and dN/dS analysis of positive selection signatures for 497 single-copy orthogroups identified gene functions mostly related to plant development and stress response. The P. mariana genome assembly and annotation provides a valuable resource for forest genetics research and applications in this broadly distributed species, especially in relation to climate adaptation.

A high-throughput pipeline for DNA/RNA/small RNA purification from tissue samples for sequencing.

High-throughput total nucleic acid (TNA) purification methods based on solid-phase reversible immobilization (SPRI) beads produce TNA suitable for both genomic and transcriptomic applications. Even so, small RNA species, including miRNA, bind weakly to SPRI beads under standard TNA purification conditions, necessitating a separate workflow using column-based methods that are difficult to automate. Here, an SPRI-based high-throughput TNA purification protocol that recovers DNA, RNA and small RNA, called GSC-modified RLT+ Aline bead-based protocol (GRAB-ALL), which incorporates modifications to enhance small RNA recovery is presented. GRAB-ALL was benchmarked against existing nucleic acid purification workflows and GRAB-ALL efficiently purifies TNA, including small RNA, for next-generation sequencing applications in a plate-based format suitable for automated high-throughput sample preparation.

The genome sequence of the Loggerhead sea turtle, Caretta caretta Linnaeus 1758.

We present a genome assembly of Caretta caretta (the Loggerhead sea turtle; Chordata, Testudines, Cheloniidae), generated from genomic data from two unrelated females. The genome sequence is 2.13 gigabases in size. The assembly has a busco completion score of 96.1% and N50 of 130.95 Mb. The majority of the assembly is scaffolded into 28 chromosomal representations with a remaining 2% of the assembly being excluded from these.

Genetic subdivisions of follicular lymphoma defined by distinct coding and noncoding mutation patterns.

Follicular lymphoma (FL) accounts for ∼20% of all new lymphoma cases. Increases in cytological grade are a feature of the clinical progression of this malignancy, and eventual histologic transformation (HT) to the aggressive diffuse large B-cell lymphoma (DLBCL) occurs in up to 15% of patients. Clinical or genetic features to predict the risk and timing of HT have not been described comprehensively. In this study, we analyzed whole-genome sequencing data from 423 patients to compare the protein coding and noncoding mutation landscapes of untransformed FL, transformed FL, and de novo DLBCL. This revealed 2 genetically distinct subgroups of FL, which we have named DLBCL-like (dFL) and constrained FL (cFL). Each subgroup has distinguishing mutational patterns, aberrant somatic hypermutation rates, and biological and clinical characteristics. We implemented a machine learning-derived classification approach to stratify patients with FL into cFL and dFL subgroups based on their genomic features. Using separate validation cohorts, we demonstrate that cFL status, whether assigned with this full classifier or a single-gene approximation, is associated with a reduced rate of HT. This implies distinct biological features of cFL that constrain its evolution, and we highlight the potential for this classification to predict HT from genetic features present at diagnosis.

Defining the heterogeneity of unbalanced structural variation underlying breast cancer susceptibility by nanopore genome sequencing.

Germline structural variants (SVs) are challenging to resolve by conventional genetic testing assays. Long-read sequencing has improved the global characterization of SVs, but its sensitivity at cancer susceptibility loci has not been reported. Nanopore long-read genome sequencing was performed for nineteen individuals with pathogenic copy number alterations in BRCA1, BRCA2, CHEK2 and PALB2 identified by prior clinical testing. Fourteen variants, which spanned single exons to whole genes and included a tandem duplication, were accurately represented. Defining the precise breakpoints of SVs in BRCA1 and CHEK2 revealed unforeseen allelic heterogeneity and informed the mechanisms underlying the formation of recurrent deletions. Integrating read-based and statistical phasing further helped define extended haplotypes associated with founder alleles. Long-read sequencing is a sensitive method for characterizing private, recurrent and founder SVs underlying breast cancer susceptibility. Our findings demonstrate the potential for nanopore sequencing as a powerful genetic testing assay in the hereditary cancer setting.

Molecular determinants of clinical outcomes in a real-world diffuse large B-cell lymphoma population.

Molecular heterogeneity of diffuse large B-cell lymphoma (DLBCL) underlies the variable outcomes achieved with immunochemotherapy. However, outcomes of gene expression profiling (GEP)-defined molecular subgroups in a real-world DLBCL population remain unknown. Here we examined the prevalence and outcomes of molecular subgroups in an unselected population of 1149 patients with de novo DLBCL in British Columbia, Canada. Evaluable biopsies were profiled by fluorescence in situ hybridization (FISH), immunohistochemistry, and digital GEP to assign cell-of-origin and the so-called "double-hit signature" (DHITsig)-a signature originally described as being characteristic for high-grade B-cell lymphoma with MYC and BCL2 rearrangements (HGBCL-DH-BCL2). DHITsig was expressed in 21% of 431 germinal center B-cell-like (GCB)-DLBCL and all 55 Burkitt lymphomas examined. Reflecting this latter finding, DHITsig has been renamed the "dark zone signature" (DZsig). DZsigpos-DLBCL, non-DZsigpos GCB-DLBCL and activated B-cell-like (ABC)-DLBCL were associated with a 2 year overall survival of 57%, 89%, and 71%, respectively. 62% of DZsigpos tumors were negative for HGBCL-DH-BCL2 by FISH, but were associated with outcomes similar to HGBCL-DH-BCL2. A small group of HGBCL-DH-BCL2 that lacked DZsig expression had different molecular features compared with DZsig-expressing HGBCL-DH-BCL2 and were associated with favorable outcomes comparable to DLBCL, not otherwise specified. DZsigpos and ABC-DLBCL had a shorter diagnosis-to-treatment interval (DTI) than GCB-DLBCL, with this metric being associated with outcome. In conclusion, DZsig expression extends beyond HGBCL-DH-BCL2 and captures a poor-prognosis DLBCL subgroup with short DTI, including patients unidentifiable by routine FISH testing, that should be considered for treatment intensification or novel therapies in prospective trials.

Genetic subgroups inform on pathobiology in adult and pediatric Burkitt lymphoma.

Burkitt lymphoma (BL) accounts for most pediatric non-Hodgkin lymphomas, being less common but significantly more lethal when diagnosed in adults. Much of the knowledge of the genetics of BL thus far has originated from the study of pediatric BL (pBL), leaving its relationship to adult BL (aBL) and other adult lymphomas not fully explored. We sought to more thoroughly identify the somatic changes that underlie lymphomagenesis in aBL and any molecular features that associate with clinical disparities within and between pBL and aBL. Through comprehensive whole-genome sequencing of 230 BL and 295 diffuse large B-cell lymphoma (DLBCL) tumors, we identified additional significantly mutated genes, including more genetic features that associate with tumor Epstein-Barr virus status, and unraveled new distinct subgroupings within BL and DLBCL with 3 predominantly comprising BLs: DGG-BL (DDX3X, GNA13, and GNAI2), IC-BL (ID3 and CCND3), and Q53-BL (quiet TP53). Each BL subgroup is characterized by combinations of common driver and noncoding mutations caused by aberrant somatic hypermutation. The largest subgroups of BL cases, IC-BL and DGG-BL, are further characterized by distinct biological and gene expression differences. IC-BL and DGG-BL and their prototypical genetic features (ID3 and TP53) had significant associations with patient outcomes that were different among aBL and pBL cohorts. These findings highlight shared pathogenesis between aBL and pBL, and establish genetic subtypes within BL that serve to delineate tumors with distinct molecular features, providing a new framework for epidemiologic, diagnostic, and therapeutic strategies.

Integrative analysis of KRAS wildtype metastatic pancreatic ductal adenocarcinoma reveals mutation and expression-based similarities to cholangiocarcinoma.

Oncogenic KRAS mutations are absent in approximately 10% of patients with metastatic pancreatic ductal adenocarcinoma (mPDAC) and may represent a subgroup of mPDAC with therapeutic options beyond standard-of-care cytotoxic chemotherapy. While distinct gene fusions have been implicated in KRAS wildtype mPDAC, information regarding other types of mutations remain limited, and gene expression patterns associated with KRAS wildtype mPDAC have not been reported. Here, we leverage sequencing data from the PanGen trial to perform comprehensive characterization of the molecular landscape of KRAS wildtype mPDAC and reveal increased frequency of chr1q amplification encompassing transcription factors PROX1 and NR5A2. By leveraging data from colorectal adenocarcinoma and cholangiocarcinoma samples, we highlight similarities between cholangiocarcinoma and KRAS wildtype mPDAC involving both mutation and expression-based signatures and validate these findings using an independent dataset. These data further establish KRAS wildtype mPDAC as a unique molecular entity, with therapeutic opportunities extending beyond gene fusion events.

Spruce giga-genomes: structurally similar yet distinctive with differentially expanding gene families and rapidly evolving genes.

Spruces (Picea spp.) are coniferous trees widespread in boreal and mountainous forests of the northern hemisphere, with large economic significance and enormous contributions to global carbon sequestration. Spruces harbor very large genomes with high repetitiveness, hampering their comparative analysis. Here, we present and compare the genomes of four different North American spruces: the genome assemblies for Engelmann spruce (Picea engelmannii) and Sitka spruce (Picea sitchensis) together with improved and more contiguous genome assemblies for white spruce (Picea glauca) and for a naturally occurring introgress of these three species known as interior spruce (P. engelmannii × glauca × sitchensis). The genomes were structurally similar, and a large part of scaffolds could be anchored to a genetic map. The composition of the interior spruce genome indicated asymmetric contributions from the three ancestral genomes. Phylogenetic analysis of the nuclear and organelle genomes revealed a topology indicative of ancient reticulation. Different patterns of expansion of gene families among genomes were observed and related with presumed diversifying ecological adaptations. We identified rapidly evolving genes that harbored high rates of non-synonymous polymorphisms relative to synonymous ones, indicative of positive selection and its hitchhiking effects. These gene sets were mostly distinct between the genomes of ecologically contrasted species, and signatures of convergent balancing selection were detected. Stress and stimulus response was identified as the most frequent function assigned to expanding gene families and rapidly evolving genes. These two aspects of genomic evolution were complementary in their contribution to divergent evolution of presumed adaptive nature. These more contiguous spruce giga-genome sequences should strengthen our understanding of conifer genome structure and evolution, as their comparison offers clues into the genetic basis of adaptation and ecology of conifers at the genomic level. They will also provide tools to better monitor natural genetic diversity and improve the management of conifer forests. The genomes of four closely related North American spruces indicate that their high similarity at the morphological level is paralleled by the high conservation of their physical genome structure. Yet, the evidence of divergent evolution is apparent in their rapidly evolving genomes, supported by differential expansion of key gene families and large sets of genes under positive selection, largely in relation to stimulus and environmental stress response.

The impact of whole genome and transcriptome analysis (WGTA) on predictive biomarker discovery and diagnostic accuracy of advanced malignancies.

In this study, we evaluate the impact of whole genome and transcriptome analysis (WGTA) on predictive molecular profiling and histologic diagnosis in a cohort of advanced malignancies. WGTA was used to generate reports including molecular alterations and site/tissue of origin prediction. Two reviewers analyzed genomic reports, clinical history, and tumor pathology. We used National Comprehensive Cancer Network (NCCN) consensus guidelines, Food and Drug Administration (FDA) approvals, and provincially reimbursed treatments to define genomic biomarkers associated with approved targeted therapeutic options (TTOs). Tumor tissue/site of origin was reassessed for most cases using genomic analysis, including a machine learning algorithm (Supervised Cancer Origin Prediction Using Expression [SCOPE]) trained on The Cancer Genome Atlas data. WGTA was performed on 652 cases, including a range of primary tumor types/tumor sites and 15 malignant tumors of uncertain histogenesis (MTUH). At the time WGTA was performed, alterations associated with an approved TTO were identified in 39 (6%) cases; 3 of these were not identified through routine pathology workup. In seven (1%) cases, the pathology workup either failed, was not performed, or gave a different result from the WGTA. Approved TTOs identified by WGTA increased to 103 (16%) when applying 2021 guidelines. The histopathologic diagnosis was reviewed in 389 cases and agreed with the diagnostic consensus after WGTA in 94% of non-MTUH cases (n = 374). The remainder included situations where the morphologic diagnosis was changed based on WGTA and clinical data (0.5%), or where the WGTA was non-contributory (5%). The 15 MTUH were all diagnosed as specific tumor types by WGTA. Tumor board reviews including WGTA agreed with almost all initial predictive molecular profile and histopathologic diagnoses. WGTA was a powerful tool to assign site/tissue of origin in MTUH. Current efforts focus on improving therapeutic predictive power and decreasing cost to enhance use of WGTA data as a routine clinical test.

Optimization of magnetic bead-based nucleic acid extraction for SARS-CoV-2 testing using readily available reagents.

The COVID-19 pandemic has highlighted the need for generic reagents and flexible systems in diagnostic testing. Magnetic bead-based nucleic acid extraction protocols using 96-well plates on open liquid handlers are readily amenable to meet this need. Here, one such approach is rigorously optimized to minimize cross-well contamination while maintaining sensitivity.

Proteotranscriptomic classification and characterization of pancreatic neuroendocrine neoplasms.

Pancreatic neuroendocrine neoplasms (PNENs) are biologically and clinically heterogeneous. Here, we use a multi-omics approach to uncover the molecular factors underlying this heterogeneity. Transcriptomic analysis of 84 PNEN specimens, drawn from two cohorts, is substantiated with proteomic profiling and identifies four subgroups: Proliferative, PDX1-high, Alpha cell-like and Stromal/Mesenchymal. The Proliferative subgroup, consisting of both well- and poorly differentiated specimens, is associated with inferior overall survival probability. The PDX1-high and Alpha cell-like subgroups partially resemble previously described subtypes, and we further uncover distinctive metabolism-related features in the Alpha cell-like subgroup. The Stromal/Mesenchymal subgroup exhibits molecular characteristics of YAP1/WWTR1(TAZ) activation suggestive of Hippo signaling pathway involvement in PNENs. Whole-exome sequencing reveals subgroup-enriched mutational differences, supported by activity inference analysis, and identifies hypermorphic proto-oncogene variants in 14.3% of sequenced PNENs. Our study reveals differences in cellular signaling axes that provide potential directions for PNEN patient stratification and treatment strategies.

Modulation of the Host Cell Transcriptome and Epigenome by Fusobacterium nucleatum.

Fusobacterium nucleatum is a ubiquitous opportunistic pathogen with an emerging role as an oncomicrobe in colorectal cancer and other cancer settings. F. nucleatum can adhere to and invade host cells in a manner that varies across F. nucleatum strains and host cell phenotypes. Here, we performed pairwise cocultures between three F. nucleatum strains and two immortalized primary host cell types (human colonic epithelial [HCE] cells and human carotid artery endothelial [HCAE] cells) followed by transcriptome sequencing (RNA-seq) and chromatin immunoprecipitation sequencing (ChIP-seq) to investigate transcriptional and epigenetic host cell responses. We observed that F. nucleatum-induced host cell transcriptional modulation involves strong upregulation of genes related to immune migration and inflammatory processes, such as TNF, CXCL8, CXCL1, and CCL20. Furthermore, we identified genes strongly upregulated in a cell line-specific manner. In HCE cells, overexpressed genes included UBD and DUOX2/DUOXA2, associated with p53 degradation-mediated proliferation and intestinal reactive oxygen species (ROS) production, respectively. In HCAE cells, overexpressed genes included EFNA1 and LIF, two genes commonly upregulated in colorectal cancer and associated with poor patient outcomes, and PTGS2 (COX2), a gene associated with the protective effect of aspirin in the colorectal cancer setting. Interestingly, we also observed downregulation of numerous histone modification genes upon F. nucleatum exposure. We used the ChIP-seq data to annotate chromatin states genome wide and found significant chromatin remodeling following F. nucleatum exposure in HCAE cells, with increased frequencies of active enhancer and low-signal/quiescent states. Thus, our results highlight increased inflammation and chemokine gene expression as conserved host cell responses to F. nucleatum exposure and extensive host cell epigenomic changes specific to host cell type. IMPORTANCE Fusobacterium nucleatum is a bacterium normally found in the healthy oral cavity but also has an emerging role in colorectal cancer and other cancer settings. The host-microbe interactions of F. nucleatum and its involvement in tumor initiation, progression, and treatment resistance are not fully understood. We explored host cell changes that occur in response to F. nucleatum. We identified key genes differentially expressed in response to various conditions of F. nucleatum exposure and determined that the conserved host cell response to F. nucleatum was dominated by increased inflammation and chemokine gene expression. Additionally, we found extensive host cell epigenomic changes as a novel aspect of host modulation associated with F. nucleatum exposure. These results extend our understanding of F. nucleatum as an emerging pathogen and highlight the importance of considering strain heterogeneity and host cell phenotypic variation when exploring pathogenic mechanisms of F. nucleatum.