Oncogenic Kras drives invasion and maintains metastases in colorectal cancer.

Human colorectal cancer (CRC) is a major cause of cancer mortality and frequently harbors activating mutations in the KRAS gene. To understand the role of oncogenic KRAS in CRC, we engineered a mouse model of metastatic CRC that harbors an inducible oncogenic Kras allele (Krasmut ) and conditional null alleles of Apc and Trp53 (iKAP). The iKAP model recapitulates tumor progression from adenoma through metastases. Whole-exome sequencing revealed that the Krasmut allele was heterogenous in primary tumors yet homogenous in metastases, a pattern consistent with activated Krasmut signaling being a driver of progression to metastasis. System-level and functional analyses revealed the TGF-ß pathway as a key mediator of Krasmut -driven invasiveness. Genetic extinction of Krasmut resulted in specific elimination of the Krasmut subpopulation in primary and metastatic tumors, leading to apoptotic elimination of advanced invasive and metastatic disease. This faithful CRC model provides genetic evidence that Krasmut drives CRC invasion and maintenance of metastases.

Metagenomes of rectal swabs in larger, advanced stage cervical cancers have enhanced mucus degrading functionalities and distinct taxonomic structure.

BACKGROUND: Gut microbiome community composition differs between cervical cancer (CC) patients and healthy controls, and increased gut diversity is associated with improved outcomes after treatment. We proposed that functions of specific microbial species adjoining the mucus layer may directly impact the biology of CC. METHOD: Metagenomes of rectal swabs in 41 CC patients were examined by whole-genome shotgun sequencing to link taxonomic structures, molecular functions, and metabolic pathway to patient's clinical characteristics. RESULTS: Significant association of molecular functions encoded by the metagenomes was found with initial tumor size and stage. Profiling of the molecular function abundances and their distributions identified 2 microbial communities co-existing in each metagenome but having distinct metabolism and taxonomic structures. Community A (Clostridia and Proteobacteria predominant) was characterized by high activity of pathways involved in stress response, mucus glycan degradation and utilization of degradation byproducts. This community was prevalent in patients with larger, advanced stage tumors. Conversely, community B (Bacteroidia predominant) was characterized by fast growth, active oxidative phosphorylation, and production of vitamins. This community was prevalent in patients with smaller, early-stage tumors. CONCLUSIONS: In this study, enrichment of mucus degrading microbial communities in rectal metagenomes of CC patients was associated with larger, more advanced stage tumors.

Genomic profiling of multifocal intrahepatic cholangiocarcinoma reveals intraindividual concordance of genetic alterations.

In multifocal intrahepatic cholangiocarcinoma (IHC), intrahepatic metastases (IM) represent a contraindication to surgical resection, whereas satellite nodules (SN) do not. However, no consensus criteria exist to distinguish IM from SN. The purpose of this study was to determine genetic alterations and clonal relationships in surgically resected multifocal IHC. Next-generation sequencing of 34 spatially separated IHC tumors was performed using a targeted panel of 201 cancer-associated genes. Proposed definitions in the literature were applied of SN located in the same liver segment and ≤2 cm from the primary tumor; and IM located in a different liver segment and/or >2 cm from the primary tumor. Somatic point mutations concordant across tumors from individual patients included BAP1, SMARCA4 and IDH1. Small insertions and deletions (indels) present at the same genome positions among all tumors from individuals included indels in DNA repair genes, CHEK1, ERCC5, ATR and MSH6. Copy number alterations were also similar between all tumors in each patient. In this cohort of multifocal IHC, genomic profiles were concordant across all tumors in each patient, suggesting a common progenitor cell origin, regardless of the location of tumors in the liver. The decision to perform surgery should not be based upon a perceived distinction between IM and SN.

Novel algorithmic approach predicts tumor mutation load and correlates with immunotherapy clinical outcomes using a defined gene mutation set.

BACKGROUND: While clinical outcomes following immunotherapy have shown an association with tumor mutation load using whole exome sequencing (WES), its clinical applicability is currently limited by cost and bioinformatics requirements. METHODS: We developed a method to accurately derive the predicted total mutation load (PTML) within individual tumors from a small set of genes that can be used in clinical next generation sequencing (NGS) panels. PTML was derived from the actual total mutation load (ATML) of 575 distinct melanoma and lung cancer samples and validated using independent melanoma (n = 312) and lung cancer (n = 217) cohorts. The correlation of PTML status with clinical outcome, following distinct immunotherapies, was assessed using the Kaplan-Meier method. RESULTS: PTML (derived from 170 genes) was highly correlated with ATML in cutaneous melanoma and lung adenocarcinoma validation cohorts (R2 = 0.73 and R2 = 0.82, respectively). PTML was strongly associated with clinical outcome to ipilimumab (anti-CTLA-4, three cohorts) and adoptive T-cell therapy (1 cohort) clinical outcome in melanoma. Clinical benefit from pembrolizumab (anti-PD-1) in lung cancer was also shown to significantly correlate with PTML status (log rank P value < 0.05 in all cohorts). CONCLUSIONS: The approach of using small NGS gene panels, already applied to guide employment of targeted therapies, may have utility in the personalized use of immunotherapy in cancer.

Insights from 20 years of bacterial genome sequencing.

Since the first two complete bacterial genome sequences were published in 1995, the science of bacteria has dramatically changed. Using third-generation DNA sequencing, it is possible to completely sequence a bacterial genome in a few hours and identify some types of methylation sites along the genome as well. Sequencing of bacterial genome sequences is now a standard procedure, and the information from tens of thousands of bacterial genomes has had a major impact on our views of the bacterial world. In this review, we explore a series of questions to highlight some insights that comparative genomics has produced. To date, there are genome sequences available from 50 different bacterial phyla and 11 different archaeal phyla. However, the distribution is quite skewed towards a few phyla that contain model organisms. But the breadth is continuing to improve, with projects dedicated to filling in less characterized taxonomic groups. The clustered regularly interspaced short palindromic repeats (CRISPR)-Cas system provides bacteria with immunity against viruses, which outnumber bacteria by tenfold. How fast can we go? Second-generation sequencing has produced a large number of draft genomes (close to 90 % of bacterial genomes in GenBank are currently not complete); third-generation sequencing can potentially produce a finished genome in a few hours, and at the same time provide methlylation sites along the entire chromosome. The diversity of bacterial communities is extensive as is evident from the genome sequences available from 50 different bacterial phyla and 11 different archaeal phyla. Genome sequencing can help in classifying an organism, and in the case where multiple genomes of the same species are available, it is possible to calculate the pan- and core genomes; comparison of more than 2000 Escherichia coli genomes finds an E. coli core genome of about 3100 gene families and a total of about 89,000 different gene families. Why do we care about bacterial genome sequencing? There are many practical applications, such as genome-scale metabolic modeling, biosurveillance, bioforensics, and infectious disease epidemiology. In the near future, high-throughput sequencing of patient metagenomic samples could revolutionize medicine in terms of speed and accuracy of finding pathogens and knowing how to treat them.

Metabolic environments and genomic features associated with pathogenic and mutualistic interactions between bacteria and plants.

Genomic characteristics discriminating parasitic and mutualistic relationship of bacterial symbionts with plants are poorly understood. This study comparatively analyzed the genomes of 54 mutualists and pathogens to discover genomic markers associated with the different phenotypes. Using metabolic network models, we predict external environments associated with free-living and symbiotic lifestyles and quantify dependences of symbionts on the host in terms of the consumed metabolites. We show that specific differences between the phenotypes are pronounced at the levels of metabolic enzymes, especially carbohydrate active, and protein functions. Overall, biosynthetic functions are enriched and more diverse in plant mutualists whereas processes and functions involved in degradation and host invasion are enriched and more diverse in pathogens. A distinctive characteristic of plant pathogens is a putative novel secretion system with a circadian rhythm regulator. A specific marker of plant mutualists is the co-residence of genes encoding nitrogenase and ribulose bisphosphate carboxylase/oxygenase (RuBisCO). We predict that RuBisCO is likely used in a putative metabolic pathway to supplement carbon obtained heterotrophically with low-cost assimilation of carbon from CO2. We validate results of the comparative analysis by predicting correct phenotype, pathogenic or mutualistic, for 20 symbionts in an independent set of 30 pathogens, mutualists, and commensals.

Analyzing large biological datasets with association networks.

Due to advances in high-throughput biotechnologies biological information is being collected in databases at an amazing rate, requiring novel computational approaches that process collected data into new knowledge in a timely manner. In this study, we propose a computational framework for discovering modular structure, relationships and regularities in complex data. The framework utilizes a semantic-preserving vocabulary to convert records of biological annotations of an object, such as an organism, gene, chemical or sequence, into networks (Anets) of the associated annotations. An association between a pair of annotations in an Anet is determined by the similarity of their co-occurrence pattern with all other annotations in the data. This feature captures associations between annotations that do not necessarily co-occur with each other and facilitates discovery of the most significant relationships in the collected data through clustering and visualization of the Anet. To demonstrate this approach, we applied the framework to the analysis of metadata from the Genomes OnLine Database and produced a biological map of sequenced prokaryotic organisms with three major clusters of metadata that represent pathogens, environmental isolates and plant symbionts.

Mutant alcohol dehydrogenase leads to improved ethanol tolerance in Clostridium thermocellum.

Clostridium thermocellum is a thermophilic, obligately anaerobic, gram-positive bacterium that is a candidate microorganism for converting cellulosic biomass into ethanol through consolidated bioprocessing. Ethanol intolerance is an important metric in terms of process economics, and tolerance has often been described as a complex and likely multigenic trait for which complex gene interactions come into play. Here, we resequence the genome of an ethanol-tolerant mutant, show that the tolerant phenotype is primarily due to a mutated bifunctional acetaldehyde-CoA/alcohol dehydrogenase gene (adhE), hypothesize based on structural analysis that cofactor specificity may be affected, and confirm this hypothesis using enzyme assays. Biochemical assays confirm a complete loss of NADH-dependent activity with concomitant acquisition of NADPH-dependent activity, which likely affects electron flow in the mutant. The simplicity of the genetic basis for the ethanol-tolerant phenotype observed here informs rational engineering of mutant microbial strains for cellulosic ethanol production.

Intratumoral microbiome of adenoid cystic carcinomas and comparison with other head and neck cancers.

Adenoid cystic carcinoma (ACC) is a rare, usually slow-growing yet aggressive head and neck malignancy. Despite its clinical significance, our understanding of the cellular evolution and microenvironment in ACC remains limited. We investigated the intratumoral microbiomes of 50 ACC tumor tissues and 33 adjacent normal tissues using 16S rRNA gene sequencing. This allowed us to characterize the bacterial communities within the ACC and explore potential associations between the bacterial community structure, patient clinical characteristics, and tumor molecular features obtained through RNA sequencing. The bacterial composition in the ACC was significantly different from that in adjacent normal salivary tissue, and the ACC exhibited diverse levels of species richness. We identified two main microbial subtypes within the ACC: oral-like and gut-like. Oral-like microbiomes, characterized by increased diversity and abundance of Neisseria, Leptotrichia, Actinomyces, Streptococcus, Rothia, and Veillonella (commonly found in healthy oral cavities), were associated with a less aggressive ACC-II molecular subtype and improved patient outcomes. Notably, we identified the same oral genera in oral cancer and head and neck squamous cell carcinomas. In both cancers, they were part of shared oral communities associated with a more diverse microbiome, less aggressive tumor phenotype, and better survival that reveal the genera as potential pancancer biomarkers for favorable microbiomes in ACC and other head and neck cancers. Conversely, gut-like intratumoral microbiomes, which feature low diversity and colonization by gut mucus layer-degrading species, such as Bacteroides, Akkermansia, Blautia, Bifidobacterium, and Enterococcus, were associated with poorer outcomes. Elevated levels of Bacteroides thetaiotaomicron were independently associated with significantly worse survival and positively correlated with tumor cell biosynthesis of glycan-based cell membrane components.

Exploratory analysis of the cervix tumoral HPV antigen-specific T-cell repertoire during chemoradiation and after brachytherapy.

BACKGROUND: Chemoradiation (CRT) may modulate the immune milieu as an in-situ vaccine. Rapid dose delivery of brachytherapy has unclear impact on T-cell repertoires. HPV-associated cancers express viral oncoproteins E6/E7, which enable tracking antigen/tumor-specific immunity during CRT. METHODS: Thirteen cervical cancer patients on a multi-institutional prospective protocol from 1/2020-1/2023 underwent standard-of-care CRT with pulsed-dose-rate brachytherapy boost (2 fractions). Cervix swabs at various timepoints underwent multiplex DNA deep sequencing of the TCR-ß/CDR3 region with immunoSEQ. Separately, HPV-responsive T-cell clones were also expanded ex vivo. Statistical analysis was via Mann-Whitney-U. RESULTS: TCR productive clonality, templates, frequency, or rearrangements increased post-brachytherapy in 8 patients. Seven patients had E6/E7-responsive evolution over CRT with increased productive templates (ranges: 1.2-50.2 fold-increase from baseline), frequency (1.2-1.7), rearrangements (1.2-40.2), and clonality (1.2-15.4). Five patients had HPV-responsive clonal expansion post-brachytherapy, without changes in HPV non-responsive clones. Epitope mapping revealed VDJ rearrangements targeting cervical cancer-associated antigens in 5 patients. The only two patients with disease recurrence lacked response in all metrics. A lack of global TCR remodeling correlated with worse recurrence-free survival, p = 0.04. CONCLUSION: CRT and brachytherapy alters the cervical cancer microenvironment to facilitate the expansion of specific T-cell populations, which may contribute to treatment efficacy.

Synthetic PET from CT improves diagnosis and prognosis for lung cancer: Proof of concept.

[18F]Fluorodeoxyglucose positron emission tomography (FDG-PET) and computed tomography (CT) are indispensable components in modern medicine. Although PET can provide additional diagnostic value, it is costly and not universally accessible, particularly in low-income countries. To bridge this gap, we have developed a conditional generative adversarial network pipeline that can produce FDG-PET from diagnostic CT scans based on multi-center multi-modal lung cancer datasets (n = 1,478). Synthetic PET images are validated across imaging, biological, and clinical aspects. Radiologists confirm comparable imaging quality and tumor contrast between synthetic and actual PET scans. Radiogenomics analysis further proves that the dysregulated cancer hallmark pathways of synthetic PET are consistent with actual PET. We also demonstrate the clinical values of synthetic PET in improving lung cancer diagnosis, staging, risk prediction, and prognosis. Taken together, this proof-of-concept study testifies to the feasibility of applying deep learning to obtain high-fidelity PET translated from CT.

Enhancing NSCLC recurrence prediction with PET/CT habitat imaging, ctDNA, and integrative radiogenomics-blood insights.

While we recognize the prognostic importance of clinicopathological measures and circulating tumor DNA (ctDNA), the independent contribution of quantitative image markers to prognosis in non-small cell lung cancer (NSCLC) remains underexplored. In our multi-institutional study of 394 NSCLC patients, we utilize pre-treatment computed tomography (CT) and 18F-fluorodeoxyglucose positron emission tomography (FDG-PET) to establish a habitat imaging framework for assessing regional heterogeneity within individual tumors. This framework identifies three PET/CT subtypes, which maintain prognostic value after adjusting for clinicopathologic risk factors including tumor volume. Additionally, these subtypes complement ctDNA in predicting disease recurrence. Radiogenomics analysis unveil the molecular underpinnings of these imaging subtypes, highlighting downregulation in interferon alpha and gamma pathways in the high-risk subtype. In summary, our study demonstrates that these habitat imaging subtypes effectively stratify NSCLC patients based on their risk levels for disease recurrence after initial curative surgery or radiotherapy, providing valuable insights for personalized treatment approaches.

BESC knowledgebase public portal.

UNLABELLED: The BioEnergy Science Center (BESC) is undertaking large experimental campaigns to understand the biosynthesis and biodegradation of biomass and to develop biofuel solutions. BESC is generating large volumes of diverse data, including genome sequences, omics data and assay results. The purpose of the BESC Knowledgebase is to serve as a centralized repository for experimentally generated data and to provide an integrated, interactive and user-friendly analysis framework. The Portal makes available tools for visualization, integration and analysis of data either produced by BESC or obtained from external resources. AVAILABILITY: http://besckb.ornl.gov.

Serial Genotyping of the Human Papillomavirus in Cervical Cancer: An Insight Into Virome Dynamics During Chemoradiation Therapy.

PURPOSE: Human papillomavirus (HPV) is the primary driver of cervical cancer. Although studies in other malignancies correlated peripheral blood DNA clearance with favorable outcomes, research on the prognostic value of HPV clearance in gynecologic cancers using intratumoral HPV is scarce. We aimed to quantify the intratumoral HPV virome in patients undergoing chemoradiation therapy (CRT) and associate this with clinical characteristics and outcomes. METHODS AND MATERIALS: This prospective study enrolled 79 patients with stage IB-IVB cervical cancer undergoing definitive CRT. Cervical tumor swabs collected at baseline and week 5 (end of intensity modulated radiation therapy) were sent for shotgun metagenome sequencing and processed via VirMAP, a viral genome sequencing and identification tool for all known HPV types. The data were categorized into HPV groups (16, 18, high risk [HR], and low risk [LR]). We used independent t tests and Wilcoxon signed-rank to compare continuous variables and χ2 and Fisher exact tests to compare categorical variables. Kaplan-Meier survival modeling was performed with log-rank testing. HPV genotyping was verified using quantitative polymerase chain reaction to validate VirMAP results using receiver operating characteristic curve and Cohen's kappa. RESULTS: At baseline, 42%, 12%, 25%, and 16% of patients were positive for HPV 16, HPV 18, HPV HR, and HPV LR, respectively, and 8% were HPV negative. HPV type was associated with insurance status and CRT response. Patients with HPV 16+ and other HPV HR+ tumors were significantly more likely to have a complete response to CRT versus patients with HPV 18 and HPV LR/HPV-negative tumors. Overall HPV viral loads predominantly decreased throughout CRT, except for HPV LR viral load. CONCLUSIONS: Rarer, less well-studied HPV types in cervical tumors are clinically significant. HPV 18 and HPV LR/negative tumors are associated with poor CRT response. This feasibility study provides a framework for a larger study of intratumoral HPV profiling to predict outcomes in patients with cervical cancer.

Immune environment and antigen specificity of the T cell receptor repertoire of malignant ascites in ovarian cancer.

We evaluated the association of disease outcome with T cell immune-related characteristics and T cell receptor (TCR) repertoire in malignant ascites from patients with high-grade epithelial ovarian cancer. Ascitic fluid samples were collected from 47 high-grade epithelial ovarian cancer patients and analyzed using flow cytometry and TCR sequencing to characterize the complementarity determining region 3 TCR ß-chain. TCR functions were analyzed using the McPAS-TCR and VDJ databases. TCR clustering was implemented using Grouping of Lymphocyte Interactions by Paratope Hotspots software. Patients with poor prognosis had ascites characterized by an increased ratio of CD8+ T cells to regulatory T cells, which correlated with an increased productive frequency of the top 100 clones and decreased productive entropy. TCRs enriched in patients with an excellent or good prognosis were more likely to recognize cancer antigens and contained more TCR reads predicted to recognize epithelial ovarian cancer antigens. In addition, a TCR motif that is predicted to bind the TP53 neoantigen was identified, and this motif was enriched in patients with an excellent or good prognosis. Ascitic fluid in high-grade epithelial ovarian cancer patients with an excellent or good prognosis is enriched with TCRs that may recognize ovarian cancer-specific neoantigens, including mutated TP53 and TEAD1. These results suggest that an effective antigen-specific immune response in ascites is vital for a good outcome in high-grade epithelial ovarian cancer.

HPV-related anal cancer is associated with changes in the anorectal microbiome during cancer development.

Background: Squamous cell carcinoma of the anus (SCCA) is a rare gastrointestinal cancer. Factors associated with progression of HPV infection to anal dysplasia and cancer are unclear and screening guidelines and approaches for anal dysplasia are less clear than for cervical dysplasia. One potential contributing factor is the anorectal microbiome. In this study, we aimed to identify differences in anal microbiome composition in the settings of HPV infection, anal dysplasia, and anal cancer in this rare disease. Methods: Patients were enrolled in two prospective studies. Patients with anal dysplasia were part of a cross-sectional cohort that enrolled women with high-grade lower genital tract dysplasia. Anorectal tumor swabs were prospectively collected from patients with biopsy-confirmed locally advanced SCCA prior to receiving standard-of-care chemoradiotherapy (CRT). Patients with high-grade lower genital tract dysplasia without anal dysplasia were considered high-risk (HR Normal). 16S V4 rRNA Microbiome sequencing was performed for anal swabs. Alpha and Beta Diversity and composition were compared for HR Normal, anal dysplasia, and anal cancer. Results: 60 patients with high-grade lower genital tract dysplasia were initially enrolled. Seven patients had concurrent anal dysplasia and 44 patients were considered HR Normal. Anorectal swabs from 21 patients with localized SCCA were included, sequenced, and analyzed in the study. Analysis of weighted and unweighted UniFrac distances demonstrated significant differences in microbial community composition between anal cancer and HR normal (p=0.018). LEfSe identified that all three groups exhibited differential enrichment of specific taxa. Peptoniphilus (p=0.028), Fusobacteria (p=0.0295), Porphyromonas (p=0.034), and Prevotella (p=0.029) were enriched in anal cancer specimens when compared to HR normal. Conclusion: Although alpha diversity was similar between HR Normal, dysplasia and cancer patients, composition differed significantly between the three groups. Increased anorectal Peptoniphilus, Fusobacteria, Porphyromonas, and Prevotella abundance were associated with anal cancer. These organisms have been reported in various gastrointestinal cancers with roles in facilitating the proinflammatory microenvironment and neoplasia progression. Future work should investigate a potential role of microbiome analysis in screening for anal dysplasia and investigation into potential mechanisms of how these microbial imbalances influence the immune system and anal carcinogenesis.

Tumor-resident Lactobacillus iners confer chemoradiation resistance through lactate-induced metabolic rewiring.

Tumor microbiota can produce active metabolites that affect cancer and immune cell signaling, metabolism, and proliferation. Here, we explore tumor and gut microbiome features that affect chemoradiation response in patients with cervical cancer using a combined approach of deep microbiome sequencing, targeted bacterial culture, and in vitro assays. We identify that an obligate L-lactate-producing lactic acid bacterium found in tumors, Lactobacillus iners, is associated with decreased survival in patients, induces chemotherapy and radiation resistance in cervical cancer cells, and leads to metabolic rewiring, or alterations in multiple metabolic pathways, in tumors. Genomically similar L-lactate-producing lactic acid bacteria commensal to other body sites are also significantly associated with survival in colorectal, lung, head and neck, and skin cancers. Our findings demonstrate that lactic acid bacteria in the tumor microenvironment can alter tumor metabolism and lactate signaling pathways, causing therapeutic resistance. Lactic acid bacteria could be promising therapeutic targets across cancer types.

Feasibility of a novel non-invasive swab technique for serial whole-exome sequencing of cervical tumors during chemoradiation therapy.

BACKGROUND: Clinically relevant genetic predictors of radiation response for cervical cancer are understudied due to the morbidity of repeat invasive biopsies required to obtain genetic material. Thus, we aimed to demonstrate the feasibility of a novel noninvasive cervical swab technique to (1) collect tumor DNA with adequate throughput to (2) perform whole-exome sequencing (WES) at serial time points over the course of chemoradiation therapy (CRT). METHODS: Cervical cancer tumor samples from patients undergoing chemoradiation were collected at baseline, at week 1, week 3, and at the completion of CRT (week 5) using a noninvasive swab-based biopsy technique. Swab samples were analyzed with whole-exome sequencing (WES) with mutation calling using a custom pipeline optimized for shallow whole-exome sequencing with low tumor purity (TP). Tumor mutation changes over the course of treatment were profiled. RESULTS: 216 samples were collected and successfully sequenced for 70 patients (94% of total number of tumor samples collected). A total of 33 patients had a complete set of samples at all four time points. The mean mapping rate was 98% for all samples, and the mean target coverage was 180. Estimated TP was greater than 5% for all samples. Overall mutation frequency decreased during CRT but mapping rate and mean target coverage remained at >98% and >180 reads at week 5. CONCLUSION: This study demonstrates the feasibility and application of a noninvasive swab-based technique for WES analysis which may be applied to investigate dynamic tumor mutational changes during treatment to identify novel genes which confer radiation resistance.

Fusobacterium is enriched in oral cancer and promotes induction of programmed death-ligand 1 (PD-L1).

Recently, increased number of studies have demonstrated a relationship between the oral microbiome and development of head and neck cancer, however, there are few studies to investigate the role of oral bacteria in the context of the tumor microenvironment in a single head and neck subsite. Here, paired tumor and adjacent normal tissues from thirty-seven oral tongue squamous cell carcinoma (SCC) patients were subjected to 16S rRNA gene sequencing and whole exome sequencing (WES), in addition to RNA sequencing for tumor samples. We observed that Fusobacterium was significantly enriched in oral tongue cancer and that Rothia and Streptococcus were enriched in adjacent normal tissues. A decrease in alpha diversity was found in tumor when compared to adjacent normal tissues. While increased Fusobacterium in tumor samples was not associated with changes in immune cell infiltration, it was associated with increased PD-L1 mRNA expression. Therefore, we examined the effects of Fusobacterium on PD-L1 expression in head and neck SCC cell lines. We demonstrated that infection with Fusobacterium species can increase both PD-L1 mRNA and surface PD-L1 protein expression on head and neck cancer cell lines. The correlation between Fusobacterium and PD-L1 expression in oral tongue SCC, in conjunction with the ability of the bacterium to induce PD-L1 expression in vitro suggests a potential role for Fusobacterium on modulation of the tumor immune microenvironment in head and neck cancer.

Genomic Correlates of Outcome in Tumor-Infiltrating Lymphocyte Therapy for Metastatic Melanoma.

PURPOSE: Adoptive cell therapy (ACT) of tumor-infiltrating lymphocytes (TIL) historically yields a 40%-50% response rate in metastatic melanoma. However, the determinants of outcome are largely unknown. EXPERIMENTAL DESIGN: We investigated tumor-based genomic correlates of overall survival (OS), progression-free survival (PFS), and response to therapy by interrogating tumor samples initially collected to generate TIL infusion products. RESULTS: Whole-exome sequencing (WES) data from 64 samples indicated a positive correlation between neoantigen load and OS, but not PFS or response to therapy. RNA sequencing analysis of 34 samples showed that expression of PDE1C, RTKN2, and NGFR was enriched in responders who had improved PFS and OS. In contrast, the expression of ELFN1 was enriched in patients with unfavorable response, poor PFS and OS, whereas enhanced methylation of ELFN1 was observed in patients with favorable outcomes. Expression of ELFN1, NGFR, and PDE1C was mainly found in cancer-associated fibroblasts and endothelial cells in tumor tissues across different cancer types in publicly available single-cell RNA sequencing datasets, suggesting a role for elements of the tumor microenvironment in defining the outcome of TIL therapy. CONCLUSIONS: Our findings suggest that transcriptional features of melanomas correlate with outcomes after TIL therapy and may provide candidates to guide patient selection.