Integrated genomics point to immune vulnerabilities in pleural mesothelioma

Pleural mesothelioma is an aggressive malignancy with limited effective therapies. In order to identify therapeutic targets, we integrated SNP genotyping, sequencing and transcriptomics from tumours and low-passage patient-derived cells. Previously unrecognised deletions of SUFU locus (10q24.32), observed in 21% of 118 tumours, resulted in disordered expression of transcripts from Hedgehog pathways and the T-cell synapse including VISTA. Co-deletion of Interferon Type I genes and CDKN2A was present in half of tumours and was a predictor of poor survival. We also found previously unrecognised deletions in RB1 in 26% of cases and show sub-micromolar responses to downstream PLK1, CHEK1 and Aurora Kinase inhibitors in primary mesothelioma cells. Defects in Hippo pathways that included RASSF7 amplification and NF2 or LATS1/2 mutations were present in 50% of tumours and were accompanied by micromolar responses to the YAP1 inhibitor Verteporfin. Our results suggest new therapeutic avenues in mesothelioma and indicate targets and biomarkers for immunotherapy.

Integrated genomics and comprehensive validation reveal drivers of genomic evolution in esophageal adenocarcinoma

Esophageal adenocarcinoma (EAC) is associated with a marked genomic instability, which underlies disease progression and development of resistance to treatment. In this study, we used an integrated genomics approach to identify a genomic instability signature. Here we show that elevated expression of this signature correlates with poor survival in EAC as well as three other cancers. Knockout and overexpression screens establish the relevance of these genes to genomic instability. Indepth evaluation of three genes (TTK, TPX2 and RAD54B) confirms their role in genomic instability and tumor growth. Mutational signatures identified by whole genome sequencing and functional studies demonstrate that DNA damage and homologous recombination are common mechanisms of genomic instability induced by these genes. Our data suggest that the inhibitors of TTK and possibly other genes identified in this study have potential to inhibit/reduce growth and spontaneous as well as chemotherapy-induced genomic instability in EAC and possibly other cancers.

Functional consequences of a rare missense BARD1 c.403G>A germline mutation identified in a triple-negative breast cancer patient

We identified a rare missense germline mutation in BARD1 (c.403G>A or p.Asp135Asn) as pathogenic using integrated genomics and transcriptomics profiling of germline and tumor samples from an early-onset triple-negative breast cancer patient who later was administrated with a PARP inhibitor for 2 months. We demonstrated in cell and mouse models that, compared to the wild-type, (1) c.403G>A mutant cell lines were more sensitive to irradiation, a DNA damage agent, and a PARP inhibitor; (2) c.403G>A mutation inhibited interaction between BARD1 and RAD51 (but not BRCA1); and (3) c.403G>A mutant mice were hypersensitive to ionizing radiation. Our study shed lights on the clinical interpretation of rare germline mutations of BARD1.

Thyroid Hormone Dependent Transcriptional Programming by TRβ Requires SWI/SNF Chromatin Remodelers

Transcriptional regulation in response to thyroid hormone (T3) is a dynamic and cell-type specific process that maintains cellular homeostasis and identity in all tissues. However, our understanding of the mechanisms of thyroid hormone receptor (TR) actions at the molecular level are actively being refined. We used an integrated genomics approach to profile and characterize the cistrome of TRβ, map changes in chromatin accessibility, and capture the transcriptomic changes in response to T3 in normal thyroid cells. There are significant shifts in TRβ genomic occupancy in response to T3, which are associated with differential chromatin accessibility, and differential recruitment of SWI/SNF chromatin remodelers. We further demonstrate selective recruitment of BAF and PBAF SWI/SNF complexes to TRβ binding sites, revealing novel differential functions in regulating chromatin accessibility and gene expression. Our findings highlight three distinct modes of TRβ interaction with chromatin and coordination of coregulator activity.

High-Altitude Adaptation: Mechanistic Insights from Integrated Genomics and Physiology

Population genomic analyses of high-altitude humans and other vertebrates have identified numerous candidate genes for hypoxia adaptation, and the physiological pathways implicated by such analyses suggest testable hypotheses about underlying mechanisms. Studies of highland natives that integrate genomic data with experimental measures of physiological performance capacities and subordinate traits are revealing associations between genotypes (e.g., hypoxia-inducible factor gene variants) and hypoxia-responsive phenotypes. The subsequent search for causal mechanisms is complicated by the fact that observed genotypic associations with hypoxia-induced phenotypes may reflect second-order consequences of selection-mediated changes in other (unmeasured) traits that are coupled with the focal trait via feedback regulation. Manipulative experiments to decipher circuits of feedback control and patterns of phenotypic integration can help identify causal relationships that underlie observed genotype–phenotype associations. Such experiments are critical for correct inferences about phenotypic targets of selection and mechanisms of adaptation.

Molecular correlates of immune cytolytic subgroups in colorectal cancer by integrated genomics analysis

Although immune checkpoint inhibition (ICI) has shown promising results in metastatic dMMR/MSI-H colorectal cancer (CRC), the majority of pMMR/MSS patients do not respond to such therapies. To systematically evaluate the determinants of immune response in CRC, we explored whether patients with diverse levels of immune cytolytic activity (CYT) have different patterns of chromothripsis and kataegis. Analysis of CRC genomic data from the TCGA, indicated an excess of chromothriptic clusters among CYT-low colon adenocarcinomas, affecting known cancer drivers (APC, KRAS, BRAF, TP53 and FBXW7), immune checkpoints (CD274, PDCD1LG2, IDO1/2 and LAG3) and immune-related genes (ENTPD1, PRF1, NKG7, FAS, GZMA/B/H/K and CD73). CYT-high tumors were characterized by hypermutation, enrichment in APOBEC-associated mutations and kataegis events, as well as APOBEC activation. We also assessed differences in the most prevalent mutational signatures (SBS15, SBS20, SBS54 and DBS2) across cytolytic subgroups. Regarding the composition of immune cells in the tumor milieu, we found enrichment of M1 macrophages, CD8+ T cells and Tregs, as well as higher CD8+ T-cells/Tregs ratio among CYT-high tumors. CYT-high patients had higher immunophenoscores, which is predictive of their responsiveness if they were to be treated with anti-PD-1 alone or in combination with anti-CTLA-4 drugs. These results could have implications for patient responsiveness to immune checkpoint inhibitors.

Understanding the evolution of Mycobacterium tuberculosis lineages using an integrated genomics and metabolomics approach

Despite being the number one cause of death from an infectious disease, little is known of the 7 phylogenetic lineages of Mycobacterium tuberculosis (Mtb). These lineages are thought to have adapted differently to their human hosts, as they are geographically localised. As a result, they show variation at the phenotype level, such as virulence and the ability to develop antibiotic resistance, and at a genomic level, such as in single nucleotide polymorphisms (SNPs). We have linked the differences in SNPs between lineages to differences in metabolites (i.e. what is ultimately produced by a cell). Through multi-omic integration of these datasets we have discovered lineage-specific metabolomic changes, potentially as a result of genomic adaptation. The differences between lineages will provide insight into new biological pathways to target and manipulate in future research.

Integrated genomics and comprehensive validation reveal novel drivers of genomic evolution in esophageal adenocarcinoma

Identification of genes driving genomic evolution can provide novel targets for cancer treatment and prevention. Here we show identification of a genomic instability gene signature, using an integrated genomics approach. Elevated expression of this signature correlated with poor survival in esophageal adenocarcinoma (EAC) as well as three other human cancers. Knockout and overexpression screens confirmed the relevance of this signature to genomic instability. Indepth evaluation of TTK (a kinase), TPX2 (spindle assembly factor) and RAD54B (recombination protein) further confirmed their role in genomic instability and tumor growth. Mutational signatures identified by whole genome sequencing and functional studies demonstrated that DNA damage and homologous recombination were common mechanisms of genomic instability induced by these genes. Consistently, a TTK inhibitor impaired EAC cell growth in vivo, and increased chemotherapy-induced cytotoxicity while inhibiting genomic instability in surviving cells. Thus inhibitors of TTK and other genes identified in this study have potential to inhibit/delay genomic evolution and tumor growth. Such inhibitors also have potential to increase chemotherapy-induced cytotoxicity while reducing its harmful genomic impact in EAC and possibly other cancers.