The genomic and evolutionary landscapes of anaplastic thyroid carcinoma.

Anaplastic thyroid carcinoma is arguably the most lethal human malignancy. It often co-occurs with differentiated thyroid cancers, yet the molecular origins of its aggressivity are unknown. We sequenced tumor DNA from 329 regions of thyroid cancer, including 213 from patients with primary anaplastic thyroid carcinomas. We also whole genome sequenced 9 patients using multi-region sequencing of both differentiated and anaplastic thyroid cancer components. Using these data, we demonstrate thatanaplastic thyroid carcinomas have a higher burden of mutations than other thyroid cancers, with distinct mutational signatures and molecular subtypes. Further, different cancer driver genes are mutated in anaplastic and differentiated thyroid carcinomas, even those arising in a single patient. Finally, we unambiguously demonstrate that anaplastic thyroid carcinomas share a genomic origin with co-occurring differentiated carcinomas and emerge from a common malignant field through acquisition of characteristic clonal driver mutations.

Early Cancer Detection in Li-Fraumeni Syndrome with Cell-Free DNA.

People with Li-Fraumeni syndrome (LFS) harbor a germline pathogenic variant in the TP53 tumor suppressor gene, face a near 100% lifetime risk of cancer, and routinely undergo intensive surveillance protocols. Liquid biopsy has become an attractive tool for a range of clinical applications, including early cancer detection. Here, we provide a proof-of-principle for a multimodal liquid biopsy assay that integrates a targeted gene panel, shallow whole-genome, and cell-free methylated DNA immunoprecipitation sequencing for the early detection of cancer in a longitudinal cohort of 89 LFS patients. Multimodal analysis increased our detection rate in patients with an active cancer diagnosis over uni-modal analysis and was able to detect cancer-associated signal(s) in carriers prior to diagnosis with conventional screening (positive predictive value = 67.6%, negative predictive value = 96.5%). Although adoption of liquid biopsy into current surveillance will require further clinical validation, this study provides a framework for individuals with LFS. SIGNIFICANCE: By utilizing an integrated cell-free DNA approach, liquid biopsy shows earlier detection of cancer in patients with LFS compared with current clinical surveillance methods such as imaging. Liquid biopsy provides improved accessibility and sensitivity, complementing current clinical surveillance methods to provide better care for these patients. See related commentary by Latham et al., p. 23. This article is featured in Selected Articles from This Issue, p. 5.

Loss of LRP1B expression drives acquired chemo and radio-resistance in HPV-positive head and neck cancer.

OBJECTIVES: Although human papillomavirus positive (HPV+) oropharyngeal squamous cell carcinoma (OPSCC) patients typically experience excellent survival, 15-20 % of patients recur after treatment with chemotherapy and radiation. Therefore, there is a need for biomarkers of treatment failure to guide treatment intensity. MATERIALS AND METHODS: Whole genome sequencing was carried out on HPV+OPSCC patients who were primarily treated with concurrent chemotherapy (cisplatin) and radiation. We then explored whether the loss of LRP1Bwas sufficient to drive an aggressive phenotype, and promote a resistance to cisplatin and radiation therapy both in vitro using HPV+ cell lines (93VU147T, UMSCC47, UWO37 and UWO23) and in vivo. RESULTS: Through integrative genomic analysis of three HPV+OPSCC tumour datasets, we identified that deletion of LRP1B was enriched in samples that recurred following chemo-radiation. Knockdown using siRNA in four HPV+ cell lines (UWO23, UWO37, UMSCC47 and 93VU147T) resulted in increased proliferation of all cases. CRISPR/Cas9 deletion of LRP1B in the same cell line panel demonstrated increased proliferation, clonogenic growth and migration, as well as resistance to both cisplatin and radiation in LRP1B deleted cells compared to their respective non-targeting control cells. Cell line derived xenograft studies indicated that the LRP1B knockout tumours were more resistant to cisplatin and radiation therapy compared to their controls invivo. CONCLUSION: Taken together, our work implicates LRP1B deletion as a potential biomarker for identifying treatment resistant HPV+ OPSCC cases.

Germline whole genome sequencing in adults with multiple primary tumors.

Multiple primary tumors (MPTs) are a harbinger of hereditary cancer syndromes. Affected individuals often fit genetic testing criteria for a number of hereditary cancer genes and undergo multigene panel testing. Other genomic testing options, such as whole exome (WES) and whole genome sequencing (WGS) are available, but the utility of these genomic approaches as a second-tier test for those with uninformative multigene panel testing has not been explored. Here, we report our germline sequencing results from WGS in 9 patients with MPTs who had non-informative multigene panel testing. Following germline WGS, sequence (agnostic or 735 selected genes) and copy number variant (CNV) analysis was performed according to the American College of Medical Genetics (ACMG) standards and guidelines for interpreting sequence variants and reporting CNVs. In this cohort, WGS, as a second-tier test, did not identify additional pathogenic or likely pathogenic variants in cancer predisposition genes. Although we identified a CHEK2 likely pathogenic variant and a MUTYH pathogenic variant, both were previously identified in the multigene panels and were not explanatory for the presented type of tumors. CNV analysis also failed to identify any pathogenic or likely pathogenic variants in cancer predisposition genes. In summary, after multigene panel testing, WGS did not reveal any additional pathogenic variants in patients with MPTs. Our study, based on a small cohort of patients with MPT, suggests that germline gene panel testing may be sufficient to investigate these cases. Future studies with larger sample sizes may further elucidate the additional utility of WGS in MPTs.

Identifying Mechanisms of Resistance by Circulating Tumor DNA in EVOLVE, a Phase II Trial of Cediranib Plus Olaparib for Ovarian Cancer at Time of PARP Inhibitor Progression.

PURPOSE: To evaluate the use of blood cell-free DNA (cfDNA) to identify emerging mechanisms of resistance to PARP inhibitors (PARPi) in high-grade serous ovarian cancer (HGSOC). EXPERIMENTAL DESIGN: We used targeted sequencing (TS) to analyze 78 longitudinal cfDNA samples collected from 30 patients with HGSOC enrolled in a phase II clinical trial evaluating cediranib (VEGF inhibitor) plus olaparib (PARPi) after progression on PARPi alone. cfDNA was collected at baseline, before treatment cycle 2, and at end of treatment. These were compared with whole-exome sequencing (WES) of baseline tumor tissues. RESULTS: At baseline (time of initial PARPi progression), cfDNA tumor fractions were 0.2% to 67% (median, 3.25%), and patients with high ctDNA levels (>15%) had a higher tumor burden (sum of target lesions; P = 0.043). Across all timepoints, cfDNA detected 74.4% of mutations known from prior tumor WES, including three of five expected BRCA1/2 reversion mutations. In addition, cfDNA identified 10 novel mutations not detected by WES, including seven TP53 mutations annotated as pathogenic by ClinVar. cfDNA fragmentation analysis attributed five of these novel TP53 mutations to clonal hematopoiesis of indeterminate potential (CHIP). At baseline, samples with significant differences in mutant fragment size distribution had shorter time to progression (P = 0.001). CONCLUSIONS: Longitudinal testing of cfDNA by TS provides a noninvasive tool for detection of tumor-derived mutations and mechanisms of PARPi resistance that may aid in directing patients to appropriate therapeutic strategies. With cfDNA fragmentation analyses, CHIP was identified in several patients and warrants further investigation.

Integrated, Longitudinal Analysis of Cell-free DNA in Uveal Melanoma.

Uveal melanomas are rare tumors arising from melanocytes that reside in the eye. Despite surgical or radiation treatment, approximately 50% of patients with uveal melanoma will progress to metastatic disease, most often to the liver. Cell-free DNA (cfDNA) sequencing is a promising technology due to the minimally invasive sample collection and ability to infer multiple aspects of tumor response. We analyzed 46 serial cfDNA samples from 11 patients with uveal melanoma over a 1-year period following enucleation or brachytherapy (n = ∼4/patient) using targeted panel, shallow whole genome, and cell-free methylated DNA immunoprecipitation sequencing. We found detection of relapse was highly variable using independent analyses (P = 0.06-0.46), whereas a logistic regression model integrating all cfDNA profiles significantly improved relapse detection (P = 0.02), with greatest power derived from fragmentomic profiles. This work provides support for the use of integrated analyses to improve the sensitivity of circulating tumor DNA detection using multi-modal cfDNA sequencing. Significance: Here, we demonstrate integrated, longitudinal cfDNA sequencing using multi-omic approaches is more effective than unimodal analysis. This approach supports the use of frequent blood testing using comprehensive genomic, fragmentomic, and epigenomic techniques.

Molecular Genetic Characteristics of FANCI, a Proposed New Ovarian Cancer Predisposing Gene.

FANCI was recently identified as a new candidate ovarian cancer (OC)-predisposing gene from the genetic analysis of carriers of FANCI c.1813C>T; p.L605F in OC families. Here, we aimed to investigate the molecular genetic characteristics of FANCI, as they have not been described in the context of cancer. We first investigated the germline genetic landscape of two sisters with OC from the discovery FANCI c.1813C>T; p.L605F family (F1528) to re-affirm the plausibility of this candidate. As we did not find other conclusive candidates, we then performed a candidate gene approach to identify other candidate variants in genes involved in the FANCI protein interactome in OC families negative for pathogenic variants in BRCA1, BRCA2, BRIP1, RAD51C, RAD51D, and FANCI, which identified four candidate variants. We then investigated FANCI in high-grade serous ovarian carcinoma (HGSC) from FANCI c.1813C>T carriers and found evidence of loss of the wild-type allele in tumour DNA from some of these cases. The somatic genetic landscape of OC tumours from FANCI c.1813C>T carriers was investigated for mutations in selected genes, copy number alterations, and mutational signatures, which determined that the profiles of tumours from carriers were characteristic of features exhibited by HGSC cases. As other OC-predisposing genes such as BRCA1 and BRCA2 are known to increase the risk of other cancers including breast cancer, we investigated the carrier frequency of germline FANCI c.1813C>T in various cancer types and found overall more carriers among cancer cases compared to cancer-free controls (p = 0.007). In these different tumour types, we also identified a spectrum of somatic variants in FANCI that were not restricted to any specific region within the gene. Collectively, these findings expand on the characteristics described for OC cases carrying FANCI c.1813C>T; p.L605F and suggest the possible involvement of FANCI in other cancer types at the germline and/or somatic level.

Serological abnormalities that predict progression to systemic autoimmune rheumatic diseases in antinuclear antibody-positive individuals.

OBJECTIVE: We investigated the autoantibody (autoAb) profiles in ANA+ individuals lacking systemic autoimmune rheumatic disease (SARD) and early SARD patients to determine the key differences between these groups and identify factors that are associated with an increased risk of symptomatic progression within the next 2 years in ANA+ individuals. METHODS: Using custom antigen (Ag) microarrays, 144 IgM and IgG autoAbs were surveyed in 84 asymptomatic and 123 symptomatic (48 UCTD and 75 SARD patients) ANA+ individuals. AutoAbs were compared in ANA+ individuals lacking a SARD diagnosis with ≥2 years follow-up (n = 52), including all those who demonstrated progression (n = 14) during this period, with changes over time assessed in a representative subset. RESULTS: We show that ANA+ individuals have autoAb to many self-Ags that are not being captured by current screening techniques and very high levels of these autoAbs are predominantly restricted to early SARD patients, with SLE patients displaying reactivity to many more autoAgs than the other groups. In general, the symptoms that developed in progressors mirrored those seen in SARD patients with similar patterns of autoAbs. Only anti-Ro52 Abs were found to predict progression (positive predictive value 46%, negative predictive value 89%). Surprisingly, over 2 years of follow-up the levels of autoAbs remained remarkably stable regardless of whether individuals progressed or not. CONCLUSION: Our findings strongly argue that development of assays with an expanded set of auto-Ags and enhanced dynamic range would improve the diagnostic and prognostic ability of autoAb testing.

Whole-genome profiling of nasopharyngeal carcinoma reveals viral-host co-operation in inflammatory NF-κB activation and immune escape.

Interplay between EBV infection and acquired genetic alterations during nasopharyngeal carcinoma (NPC) development remains vague. Here we report a comprehensive genomic analysis of 70 NPCs, combining whole-genome sequencing (WGS) of microdissected tumor cells with EBV oncogene expression to reveal multiple aspects of cellular-viral co-operation in tumorigenesis. Genomic aberrations along with EBV-encoded LMP1 expression underpin constitutive NF-κB activation in 90% of NPCs. A similar spectrum of somatic aberrations and viral gene expression undermine innate immunity in 79% of cases and adaptive immunity in 47% of cases; mechanisms by which NPC may evade immune surveillance despite its pro-inflammatory phenotype. Additionally, genomic changes impairing TGFBR2 promote oncogenesis and stabilize EBV infection in tumor cells. Fine-mapping of CDKN2A/CDKN2B deletion breakpoints reveals homozygous MTAP deletions in 32-34% of NPCs that confer marked sensitivity to MAT2A inhibition. Our work concludes that NPC is a homogeneously NF-κB-driven and immune-protected, yet potentially druggable, cancer.

Flavopiridol causes cell cycle inhibition and demonstrates anti-cancer activity in anaplastic thyroid cancer models.

Anaplastic thyroid cancer (ATC) is a rare, but nearly uniformly fatal disease that is typically resistant to chemotherapy and radiation. Alternative strategies to target this cancer at a molecular level are necessary in order to improve dismal outcomes for ATC patients. We examined the effects of flavopiridol, a CDK inhibitor, in a panel of ATC cell lines. When cell lines were treated over a ten-point concentration range, CAL62, KMH2 and BHT-101 cell lines had a sub micromolar half-maximal inhibitory concentration, while no effect was seen in the non-cancerous cell line IMR-90. Flavopiridol treatment resulted in decreased levels of the cell cycle proteins CDK9 and MCL1, and induced cell cycle arrest. Flavopiridol also decreased the in vitro ability of ATC cells to form colonies and impeded migration using a transwell migration assay. In vivo, flavopiridol decreased tumor weight and tumor volume over time in a patient-derived xenograft model of ATC. Given the observed in vitro and in vivo activity, flavopiridol warrants further investigation for treatment of ATC.

Sex differences in oncogenic mutational processes.

Sex differences have been observed in multiple facets of cancer epidemiology, treatment and biology, and in most cancers outside the sex organs. Efforts to link these clinical differences to specific molecular features have focused on somatic mutations within the coding regions of the genome. Here we report a pan-cancer analysis of sex differences in whole genomes of 1983 tumours of 28 subtypes as part of the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium. We both confirm the results of exome studies, and also uncover previously undescribed sex differences. These include sex-biases in coding and non-coding cancer drivers, mutation prevalence and strikingly, in mutational signatures related to underlying mutational processes. These results underline the pervasiveness of molecular sex differences and strengthen the call for increased consideration of sex in molecular cancer research.

Transgenerational epigenetic and transcriptomic effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin exposure in rat.

In rats, direct exposure to TCDD causes myriad toxicities. Exposed rats experience hepatotoxicity, wasting syndrome and immune suppression, amongst others. "Inherited exposure", as occurs in the F3 generation of directly exposed F0 animals, has also been shown to cause toxicity: both male and female F3 rats demonstrate an increased incidence of adult onset disease, females also display reproductive abnormalities and increased incidence of ovarian diseases while males show increased incidence of kidney disease and an altered sperm epigenome. Here, we explore the hepatic transcriptomic profile of male and female F3 Sprague-Dawley rats bred through the paternal germ line from F0 dams exposed to a single dose of TCDD (0, 30, 100, 300 or 1000 ng/kg body weight) by oral gavage. We hypothesize that RNA transcripts with altered abundance in livers of unexposed F3 progeny of treated F0 Sprague-Dawley rats may result from epigenetic modifications to the genome. We further survey patterns of differential methylation within male F3 rat testis. Female F3 rats demonstrated more TCDD-mediated hepatic transcriptomic changes than males, with differences primarily in the lowest dose group. In testis from male F3 rats, multiple olfactory receptors displayed patterns of differential methylation. Hypermethylation of Egfr and Mc5r among testes from TCDD lineage rats was observed, but without corresponding changes in hepatic mRNA abundance. Further studies examining these differences in other tissue types are warranted.

Sex disparities in head & neck cancer driver genes: An analysis of the TCGA dataset.

OBJECTIVES: Survival in head and neck squamous cell carcinoma (HNSCC) has been associated with patient sex, typically with males experiencing poorer outcomes. It is unclear if this disparity is based in divergent tumor biology. We analyzed the TCGA HNSCC cohort to uncover disparities in the somatic single nucleotide variation (SNV), copy number alteration (CNA) and mRNA abundance profiles between males and females. Critically, we stratified our results by tumor HPV status to control for this significant confounder. METHODS: SNV, CNA and mRNA abundance differences between males and females were compared separately for the HPV-positive (n = 67) and negative (n = 431) TCGA HNSCC cohorts. Overall survival outcomes were compared in males and females in both HPV-positive and HPV-negative subsets of patients. RESULTS: Females were found to have poorer overall survival than males (p = 0.048), largely due to higher rates of HPV-positive disease among men. SNV analysis revealed that in HPV-positive disease, there were no differences by sex after accounting for the false discovery rate (FDR). In HPV-negative tumors, BRWD3 mutations occurred more frequently in the tumors of female patients compared to males after adjusting for the FDR (p = 0.02). Further, HPV-negative BRWD3 mutant tumors were found to have significantly worse 5-year overall survival compared to wildtype on multivariate analysis (p = 0.02). There were 88 heterozygous deletions and 14 amplifications that were differentially altered between male and female HPV-negative tumors and associated with expression changes. Pathway analysis of these genes revealed that tumors from males were enriched in five pathways including chemokine and phosphophatidylinositol signaling. CONCLUSIONS: Reanalysis of the TCGA HNSCC dataset stratified by sex revealed that males in this cohort had a significant survival advantage, due to a higher proportion of HPV-positive disease. Mutations in BRWD3 were more frequent in HPV-negative tumors of females and were associated with poorer overall survival. BRWD3 may represent a novel biomarker of patient outcomes, but will require additional validation.

Genomic and human papillomavirus profiling of an oral cancer cohort identifies TP53 as a predictor of overall survival.

BACKGROUND: The genomic landscape of head and neck cancer has been reported through The Cancer Genome Atlas project. We attempt to determine if high-risk human papillomavirus (HPV) or frequently mutated genes are correlated with survival in an oral cancer cohort. METHODS: Patient demographic data along with data from final pathology was collected. Tumor DNA was analyzed using a custom Illumina targeted sequencing panel. Five high-risk HPV types were tested by qPCR. Statistical analyses were used to identify associations between patient outcome and mutational status. RESULTS: High-risk HPV types were identified in 7% of cases; HPV status was not associated with survival. Mutations were identified in TP53, TERT promoter, & PIK3CA. Mutations in TP53 were significantly associated with poorer overall survival on multi-variate analysis (p = 0.03). CONCLUSIONS: Mutations in TP53 were associated with poor patient survival. Expanding our sample size may identify further predictors of outcome to direct customized cancer care.

Rethinking Lupus Nephritis Classification on a Molecular Level.

The International Society of Nephrology/Renal Pathology Society (ISN/RPS) lupus nephritis (LN) classification is under reconsideration, given challenges with inter-rater reliability and resultant inconsistent relationship with treatment response. Integration of molecular classifiers into histologic evaluation can improve diagnostic precision and identify therapeutic targets. This study described the relationship between histological and molecular phenotypes and clinical responses in LN. Renal compartmental mRNA abundance was measured in 54 biopsy specimens from LN patients and correlated to ISN/RPS classification and individual histologic lesions. A subset of transcripts was also evaluated in sequential biopsies of a separate longitudinal cohort of 36 patients with paired samples obtained at the time of flare and at follow up. Unsupervised clustering based on mRNA abundance did not demonstrate a relationship with the (ISN/RPS) classification, nor did univariate statistical analysis. Exploratory analyses suggested a correlation with individual histologic lesions. Glomerular FN1 (fibronectin), SPP1 (secreted phosphoprotein 1), and LGALS3 (galectin 3) abundance correlated with disease activity and changed following treatment. Exploratory analyses suggested relationships between specific transcripts and individual histologic lesions, with the important representation of interferon-regulated genes. Our findings suggested that the current LN classification could be refined by the inclusion of molecular descriptors. Combining molecular and pathologic kidney biopsy phenotypes may hold promise to better classify disease and identify actionable treatment targets and merits further exploration in larger cohorts.

Comparative toxicoproteogenomics of mouse and rat liver identifies TCDD-resistance genes.

The aryl hydrocarbon receptor (AHR) mediates many toxic effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). However, the AHR alone does not explain the widely different outcomes among organisms. To identify the other factors involved, we evaluated three transgenic mouse lines, each expressing a different rat AHR isoform (rWT, DEL, and INS) providing widely different resistance to TCDD toxicity, as well as C57BL/6 and DBA/2 mice which exhibit a ~ tenfold divergence in TCDD sensitivity (exposures of 5-1000 µg/kg TCDD). We supplement these with whole-genome sequencing, together with transcriptomic and proteomic analyses of the corresponding rat models, Long-Evans (L-E) and Han/Wistar (H/W) rats (having a ~ 1000-fold difference in their TCDD sensitivities; 100 µg/kg TCDD), to identify genes associated with TCDD-response phenotypes. Overall, we identified up to 50% of genes with altered mRNA abundance following TCDD exposure are associated with a single AHR isoform (33.8%, 11.7%, 5.2% and 0.3% of 3076 genes altered unique to rWT, DEL, C57BL/6 and INS respectively following 1000 µg/kg TCDD). Hepatic Pxdc1 was significantly repressed in all three TCDD-sensitive animal models (C57BL/6 and rWT mice, and L-E rat) after TCDD exposure. Three genes, including Cxxc5, Sugp1 and Hgfac, demonstrated different AHRE-1 (full) motif occurrences within their promoter regions between rat strains, as well as different patterns of mRNA abundance. Several hepatic proteins showed parallel up- or downward alterations with their RNAs, with three genes (SNRK, IGTP and IMPA2) showing consistent, strain-dependent changes. These data show the value of integrating genomic, transcriptomic and proteomic evidence across multi-species models in toxicologic studies.

2,3,7,8-Tetrachlorodibenzo-p-dioxin modifies alternative splicing in mouse liver.

Alternative splicing is a co-transcriptional mechanism that generates protein diversity by including or excluding exons in different combinations, thereby expanding the diversity of protein isoforms of a single gene. Abnormalities in this process can result in deleterious effects to human health, and several xenobiotics are known to interfere with splicing regulation through multiple mechanisms. These changes could lead to human diseases such as cancer, neurological disorders, autoimmune diseases, and developmental disorders. 2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is an environmental contaminant generated as a byproduct of various industrial activities. Exposure to this dioxin has been linked to a wide range of pathologies through the alteration of multiple cellular processes. However, the effects of TCDD exposure on alternative splicing have not yet been studied. Here, we investigated whether a single po. dose of 5 µg/kg or 500 µg/kg TCDD influence hepatic alternative splicing in adult male C57BL/6Kou mouse. We identified several genes whose alternative splicing of precursor messenger RNAs was modified following TCDD exposure. In particular, we demonstrated that alternative splicing of Cyp1a1, Ahrr, and Actn1 was significantly altered after TCDD treatment. These findings show that the exposure to TCDD has an impact on alternative-splicing, and suggest a new avenue for understanding TCDD-mediated toxicity and pathogenesis.

A readout of metabolic efficiency in arylamine N-acetyltransferase-deficient mice reveals minor energy metabolism changes.

Recent studies have revealed a possible link between the activities of polymorphic arylamine N-acetyltransferases (NATs) and energy metabolism. We used a Nat1/Nat2 double knockout (KO) mouse model to demonstrate that ablation of the two Nat genes is associated with modest, intermittent alterations in respiratory exchange rate. Pyruvate tolerance tests show that double KO mice have attenuated hepatic gluconeogenesis when maintained on a high-fat/high-sucrose diet. Absence of the two Nat genes also leads to an increase in the hepatic concentration of coenzyme A in mice fed a high-fat/high-sucrose diet. Our results suggest a modest involvement of NAT in energy metabolism in mice, which is consistent with the absence of major phenotypic deregulation of energy metabolism in slow human acetylators.

Transcriptomic Impact of IMA-08401, a Novel AHR Agonist Resembling Laquinimod, on Rat Liver.

IMA-08401 (C2) is a novel aryl hydrocarbon receptor (AHR) agonist and selective AHR modulator (SAHRM) that is structurally similar to laquinimod (LAQ). Both compounds are converted to the AHR-active metabolite DELAQ (IMA-06201) in vivo. SAHRMs have been proposed as therapeutic options for various autoimmune disorders. Clinical trials on LAQ have not reported any significant toxic outcomes and C2 has shown low toxicity in rats; however, their functional resemblance to the highly toxic AHR agonist 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) raises questions. Here, we characterize the hepatic transcriptomic changes induced by acute (single-dose) and subacute exposure (repeated dosing for 5 days followed by a 5-day recovery period) to C2 in Sprague-Dawley rats. Exposure to C2 leads to activation of the AHR, as shown by altered transcription of Cyp1a1. We identify a heightened response early after exposure that drops off by day 10. Acute exposure to C2 leads to changes to transcription of genes involved in antiviral and antibacterial responses, which highlights the immunomodulator effects of this AHR agonist. Subacute exposure causes an oxidative stress response in the liver, the consequences of which require further study on target tissues such as the CNS and immune system, both of which may be compromised in this patient population.