Faculty Bibliography

Lung cancer remains the leading cause of cancer death worldwide1. With new treatment modalities, there has been a shift in focus to how we can predict who may respond to targeted treatments. Current data suggest that the human microbiome can affect lung cancer treatment through its effects on the systemic immune tone. Our group has shown that the lower airway microbiota of lung cancer patients is characterized by enrichment with oral commensals2 which triggers transcriptomic signatures (PI3K, MAPK) previously described in NSCLC 2,3. The impact of local lung dysbiosis on lung cancer progression and survival is unknown. Patients with suspicious nodules on imaging who underwent bronchoscopy were recruited. High-throughput sequencing of bacterial 16S rRNA-encoding gene amplicons was performed. Clustering was based on Dirichlet-Multinomial mixtures (DMM) modeling. RNAseq was performed on bronchial epithelial cells obtained by airway brushing. We focused our analysis on 83 NSCLC samples. Overall alpha-diversity showed that advanced stage (IIIb-VI) lower airway samples were more similar to buccal samples than local stage (I-IIIa), p<0.0001. In addition, worse 6-month and 1-year survival was associated with more similar alpha-diversity between lower airway and buccal samples (Figure 1A-D). Utilizing DMM two clusters were identified, Supraglottic-Predominant-Taxa (SPT) and Background-Predominant-Taxa (BPT). There was a significant increase in percentage of SPT in advance stage compared to local stage (p<0.008) Kaplan-Meir survival analysis shows worse survival in those with NSCLC who were clustered into the SPT group compared to BPT (p=0.0003, Figure 1E). With RNAseq, differentially expressed genes between advanced stage vs. local stage and 6-month vs. 1-year survival were not as pronounced as SPT vs. BPT (Figure 1F) suggesting that globally, transcriptomic changes between different stage and NSCLC survival were difficult to detect as compared to when airway microbiome were differentiated. In lung cancer, dysbiosis within the lower airway microenvironment, possibly by micro-aspiration, is associated with a worse 6-month and 1-year survival. This change is also associated with transcriptome changes in the local environment

Focal radiation therapy enhances systemic responses to anti-CTLA-4 antibodies in preclinical studies and in some patients with melanoma^1-3, but its efficacy in inducing systemic responses (abscopal responses) against tumors unresponsive to CTLA-4 blockade remained uncertain. Radiation therapy promotes the activation of anti-tumor T cells, an effect dependent on type I interferon induction in the irradiated tumor^4-6. The latter is essential for achieving abscopal responses in murine cancers⁶. The mechanisms underlying abscopal responses in patients treated with radiation therapy and CTLA-4 blockade remain unclear. Here we report that radiation therapy and CTLA-4 blockade induced systemic anti-tumor T cells in chemo-refractory metastatic non-small-cell lung cancer (NSCLC), where anti-CTLA-4 antibodies had failed to demonstrate significant efficacy alone or in combination with chemotherapy^7,8. Objective responses were observed in 18% of enrolled patients, and 31% had disease control. Increased serum interferon-β after radiation and early dynamic changes of blood T cell clones were the strongest response predictors, confirming preclinical mechanistic data. Functional analysis in one responding patient showed the rapid in vivo expansion of CD8 T cells recognizing a neoantigen encoded in a gene upregulated by radiation, supporting the hypothesis that one explanation for the abscopal response is radiation-induced exposure of immunogenic mutations to the immune system.

BACKGROUND:In lung cancer, upregulation of the PI3K pathway is an early event that contributes to cell proliferation, survival, and tissue invasion. Upregulation of this pathway was recently described as associated with enrichment of the lower airways with bacteria identified as oral commensals. We hypothesize that host-microbe interactions in the lower airways of subjects with lung cancer affect known cancer pathways. METHODS:Airway brushes were collected prospectively from subjects with lung nodules at time of diagnostic bronchoscopy, including 39 subjects with final lung cancer diagnoses and 36 subjects with non-cancer diagnosis. Additionally, samples from 10 healthy control subjects were included. 16S rRNA gene amplicon sequencing and paired transcriptome sequencing (RNAseq) were performed on all airway samples. In addition, an in vitro model with airway epithelial cells exposed to bacteria/bacterial products was performed. RESULTS:The composition of the lower airway transcriptome in the cancer patients was significantly different from the controls, which included upregulation of ERK and PI3K signaling pathways. The lower airways of lung cancer patients were enriched for oral taxa (Streptococcus and Veillonella), which was associated with upregulation of the ERK and PI3K signaling pathways. In vitro exposure of airway epithelial cells to Veillonella, Prevotella, and Streptococcus led to upregulation of these same signaling pathways. CONCLUSIONS:The data presented here shows that several transcriptomic signatures previously identified as relevant to lung cancer pathogenesis are associated with enrichment of the lower airway microbiota with oral commensals.

Venous thromboembolism (VTE), comprised of deep vein thrombosis and pulmonary embolism, is a common health problem both in the United States and worldwide, with significant associated morbidity and mortality. Despite multiple known genetic and situational risk factors, an estimated 30% of all events remain classified as idiopathic, demonstrating a significant knowledge gap in the pathophysiology VTE. While platelets are well established as an essential contributor to thrombus formation and there has been recent interest in the role of neutrophil extracellular traps, specific cell types and pathways involved in the pathogenesis of VTE remain uncertain. In this study, our primary aims were to define a unique transcriptional signature for VTE and to identify the types of cells and specific pathways involved in development of VTE. Whole blood was collected in PAX gene tubes and RNA sequencing for coding mRNA was performed in an unbiased manner in 201 patients with prevalent VTE as well as 43 healthy controls. We used a bioinformatics approach to develop a unique signature for VTE by identifying differentially expressed genes, developing cell-type modules, and ascertaining pathways driving differentially expressed transcripts. We performed additional analyses on subgroups of patients with idiopathic VTE, patients with incident VTE, and VTE patients matched to healthy controls by age and sex. We went on to use machine learning methods to learn models that best differentiate VTE patients from healthy controls and validated it on a left out test set within our VTE population. Genes specific to neutrophils, erythrocytes, and platelets, in that order, were most significantly upregulated in patients with VTE compared to healthy controls. Genes related to T-cells were downregulated. Pathway analysis revealed upregulated neutrophil activation and degranulation, erythrocyte differentiation and homeostasis, and platelet degranulation. A gene signature of 217 transcripts was outstanding at differentiating patients with VTE versus healthy controls (AUC 0.94). Following adjustment for age, sex, and race/ethnicity our genetic signature remained significantly robust at differentiating patients with VTE versus controls (AUC 0.83). Our expression signature remained stable across patients with idiopathic VTE (AUC 0.93), and in patients who went on to develop future VTE events (AUC 0.95). In summary, we have demonstrated a whole blood transcriptional signature for prevalent and incident VTE. Genes related to neutrophils, erythrocytes, and platelets are upregulated in patients with VTE and genes related to T-cells were downregulated. These findings suggest an active role of cell types once thought to be passively entrapped within thrombus and provide new areas of study to establish the pathophysiology of VTE

Unveiling the transcriptome of human blastocysts can provide a wealth of important information regarding early embryonic ontology. Comparing the mRNA production of embryos with normal and abnormal karyotypes allows for a deeper understanding of the protein pathways leading to viability and aberrant fetal development. In addition, identifying transcripts specific for normal or abnormal chromosome copy number could aid in the search for secreted substances that could be used to non-invasively identify embryos best suited for IVF embryo transfer. Using RNA-seq, we characterized the transcriptome of 71 normally developing human blastocysts that were karyotypically normal vs. trisomic or monosomic. Every monosomy and trisomy of the autosomal and sex chromosomes were evaluated, mostly in duplicate. We first mapped the transcriptome of three normal embryos and found that a common core of more than 3,000 genes is expressed in all embryos. These genes represent pathways related to actively dividing cells, such as ribosome biogenesis and function, spliceosome, oxidative phosphorylation, cell cycle and metabolic pathways. We then compared transcriptome profiles of aneuploid embryos to those of normal embryos. We observed that non-viable embryos had a large number of dysregulated genes, some showing a hundred-fold difference in expression. On the contrary, sex chromosome abnormalities, XO and XXX displayed transcriptomes more closely mimicking those embryos with 23 normal chromosome pairs. Intriguingly, we identified a set of commonly deregulated genes in the majority of both trisomies and monosomies. This is the first paper demonstrating a comprehensive transcriptome delineation of karyotypic abnormalities found in the human pre-implantation embryo. We believe that this information will contribute to the development of new pre-implantation genetic screening methods as well as a better understanding of the underlying developmental abnormalities of abnormal embryos, fetuses and children.

Purpose: Acrylamide, a probable human carcinogen (Group 2A), is ubiquitously present in the human environment. The compound is found in heated starchy foods, coffee, as well as cigarette smoke. Humans are also exposed to acrylamide occupationally. The carcinogenicity of acrylamide has been attributed to the effects of glycidamide, its reactive and mutagenic epoxide metabolite shown to be an inducer of rodent tumors at various anatomical sites. Covalent adducts with haemoglobin were reported in humans exposed to acrylamide, and DNA adducts in experimental systems. Methods: In order to characterize the pre-mutagenic DNA lesions and global mutation spectra induced by acrylamide and glycidamide, we combined DNA-adduct and wholeexome sequencing analyses in an established exposureclonal immortalization system based on mouse embryonic fibroblasts. Results: Sequencing and computational analysis of cell clones immortalized after exposure revealed a unique mutational signature of glycidamide, characterized by predominant T:A>A:T transversions, T:A>C:G transitions and C:G>A:T transversions in specific trinucleotide contexts, and with significant transcription strand bias. Computational interrogation of the human pan-cancer genome sequencing data (PCAWG data set) using experimentally derived glycidamide signature revealed the presence of the glycidamide signature in lung adenocarcinomas, lung squamous cell carcinomas and liver cancer, manifesting as a sub-component of the COSMIC tobacco smoking-related signature 4, as well as possibly associated with dietary intake of acrylamide. Conclusions: Our study demonstrates how the controlled experimental characterization of specific genetic damage associated with glycidamide exposure facilitates identifying corresponding patterns in cancer genome data, thereby underscoring how mutational signature laboratory experimentation contributes to the elucidation of cancer causation

Pineoblastoma is a rare and highly aggressive brain cancer of childhood, histologically belonging to the spectrum of primitive neuroectodermal tumors. Patients with germline mutations in DICER1, a ribonuclease involved in microRNA processing, have increased risk of pineoblastoma, but genetic drivers of sporadic pineoblastoma remain unknown. Here, we analyzed pediatric and adult pineoblastoma samples (n = 23) using a combination of genome-wide DNA methylation profiling and whole-exome sequencing or whole-genome sequencing. Pediatric and adult pineoblastomas showed distinct methylation profiles, the latter clustering with lower-grade pineal tumors and normal pineal gland. Recurrent variants were found in genes involved in PKA- and NF-κB signaling, as well as in chromatin remodeling genes. We identified recurrent homozygous deletions of DROSHA, acting upstream of DICER1 in microRNA processing, and a novel microduplication involving chromosomal region 1q21 containing PDE4DIP (myomegalin), comprising the ancient DUF1220 protein domain. Expresion of PDE4DIP and DUF1220 proteins was present exclusively in pineoblastoma with PDE4DIP gain.

Introduction The analysis of somatic mutational signatures in a single tumour can reveal clues on the natural history of this tumour, including past exposures to carcinogens. Characterising the mutational signatures of potential carcinogens used in industry may thus help identify cancers that are due to occupational exposure. Here, we used experimental models to define the mutational signatures of dichloromethane (DCM) and 1,2- dichloropropane (DCP), two solvents suspected to cause cholangiocarcinomas in occupational settings. Material and methods Experimentally induced tumours and matching normal tissues from mice exposed to DCP (by gavage), DCM (by gavage or inhalation) or DCP +DCM (by gavage), and spontaneous tumours from vehicle/sham-exposed mice were analysed by whole-exome sequencing. Somatic mutation calling was performed to define exome-wide mutation patterns induced by DCP and DCM in these assays. Mutation data obtained in biliary tract cancers of workers in the printing industry who have been exposed to DCP and DCM as single agents or as mixtures, as well as with public somatic mutation data from biliary tract cancers were mined for the presence of the experimentally defined DCP and DCM mutational signatures. Results and discussions Liver tumours from DCP and DCM exposed mice had distinct somatic mutations patterns compared to spontaneous liver tumours from vehicle/sham-exposed mice. While mutations in DCP-exposed mice were dominated by C:G>T:A transitions, the most frequent types of mutations in DCM-exposed mice were T:A>A:T and T:A>C:G substitutions. An average of 10.3 somatic mutations was observed per Mb in tumours from DCM-exposed mice, approximately 3- fold higher than in DCP-exposed or sham-treated mice. The mutation patterns found in DCP-exposed mice, but not DCMexposed mice, presented some similarities with the mutational signature observed in cholangiocarcinomas of Japanese patients with occupational DCP/DCM exposure history. Conclusion These results show that the analysis of tumour genomes from mouse carcinogenicity assays can support the characterisation of mutational signatures of carcinogenic compounds relevant to human exposures

BACKGROUND: Pineoblastoma is a rare and highly aggressive brain cancer of childhood, histologically belonging to the spectrum of primitive neuroectodermal tumors. Patients with germline mutations in DICER1, a ribonuclease involved in microRNA processing, have increased risk of pineoblastoma, but genetic drivers of sporadic pineoblastoma remain unknown. METHODS: We analyzed pediatric and adult pineoblastoma samples (n=23) using integrated genomic studies, including genome-wide DNA methylation profiling, whole-exome or whole-genome sequencing, and whole-transcriptome analysis. RESULTS: Pediatric and adult pineoblastomas showed distinct methylation profiles, the latter clustering with lower grade pineal tumors and normal pineal gland. Recurrent somatic mutations were found in genes involved in PKA-and NF-kappaB signaling, as well as in chromatin remodeling genes. We identified recurrent homozygous deletions of DROSHA, acting upstream of DICER1 in microRNA processing, and a novel microduplication involving chromosomal region 1q21 containing PDE4DIP (myomegalin), comprising the ancient DUF1220 protein domain. Expression of PDE4DIP and DUF1220 proteins was present exclusively in pineoblastoma with PDE4DIP gain. Whole-transcriptome analysis showed that homozygous loss of DROSHA led to distinct changes in RNA expression profile. Disruption of the DROSHA locus in human neural stem cells using the CRISPR/Cas9 system, led to decrease of the DROSHA protein, and massive loss of miRNAs. CONCLUSION: We identified recurrent homozygous deletions of DROSHA in pineoblastoma, suggesting that different mechanisms disrupting miRNA processing are involved in the pathogenesis of familial versus sporadic pineoblastoma. Furthermore, a novel microduplication of PDE4DIP leading to upregulation of DUF1220 protein suggests DUF1220 as a novel oncogenic driver in pineoblastoma