Bystander CD8+ T cells are abundant and phenotypically distinct in human tumour infiltrates.

Various forms of immunotherapy, such as checkpoint blockade immunotherapy, are proving to be effective at restoring T cell-mediated immune responses that can lead to marked and sustained clinical responses, but only in some patients and cancer types1-4. Patients and tumours may respond unpredictably to immunotherapy partly owing to heterogeneity of the immune composition and phenotypic profiles of tumour-infiltrating lymphocytes (TILs) within individual tumours and between patients5,6. Although there is evidence that tumour-mutation-derived neoantigen-specific T cells play a role in tumour control2,4,7-10, in most cases the antigen specificities of phenotypically diverse tumour-infiltrating T cells are largely unknown. Here we show that human lung and colorectal cancer CD8+ TILs can not only be specific for tumour antigens (for example, neoantigens), but also recognize a wide range of epitopes unrelated to cancer (such as those from Epstein-Barr virus, human cytomegalovirus or influenza virus). We found that these bystander CD8+ TILs have diverse phenotypes that overlap with tumour-specific cells, but lack CD39 expression. In colorectal and lung tumours, the absence of CD39 in CD8+ TILs defines populations that lack hallmarks of chronic antigen stimulation at the tumour site, supporting their classification as bystanders. Expression of CD39 varied markedly between patients, with some patients having predominantly CD39- CD8+ TILs. Furthermore, frequencies of CD39 expression among CD8+ TILs correlated with several important clinical parameters, such as the mutation status of lung tumour epidermal growth factor receptors. Our results demonstrate that not all tumour-infiltrating T cells are specific for tumour antigens, and suggest that measuring CD39 expression could be a straightforward way to quantify or isolate bystander T cells.

GWAS meta-analysis of 16 790 patients with Barrett's oesophagus and oesophageal adenocarcinoma identifies 16 novel genetic risk loci and provides insights into disease aetiology beyond the single marker level.

OBJECTIVE: Oesophageal cancer (EC) is the sixth leading cause of cancer-related deaths. Oesophageal adenocarcinoma (EA), with Barrett's oesophagus (BE) as a precursor lesion, is the most prevalent EC subtype in the Western world. This study aims to contribute to better understand the genetic causes of BE/EA by leveraging genome wide association studies (GWAS), genetic correlation analyses and polygenic risk modelling. DESIGN: We combined data from previous GWAS with new cohorts, increasing the sample size to 16 790 BE/EA cases and 32 476 controls. We also carried out a transcriptome wide association study (TWAS) using expression data from disease-relevant tissues to identify BE/EA candidate genes. To investigate the relationship with reported BE/EA risk factors, a linkage disequilibrium score regression (LDSR) analysis was performed. BE/EA risk models were developed combining clinical/lifestyle risk factors with polygenic risk scores (PRS) derived from the GWAS meta-analysis. RESULTS: The GWAS meta-analysis identified 27 BE and/or EA risk loci, 11 of which were novel. The TWAS identified promising BE/EA candidate genes at seven GWAS loci and at five additional risk loci. The LDSR analysis led to the identification of novel genetic correlations and pointed to differences in BE and EA aetiology. Gastro-oesophageal reflux disease appeared to contribute stronger to the metaplastic BE transformation than to EA development. Finally, combining PRS with BE/EA risk factors improved the performance of the risk models. CONCLUSION: Our findings provide further insights into BE/EA aetiology and its relationship to risk factors. The results lay the foundation for future follow-up studies to identify underlying disease mechanisms and improving risk prediction.

An integrative model of pathway convergence in genetically heterogeneous blast crisis chronic myeloid leukemia.

Targeted therapies against the BCR-ABL1 kinase have revolutionized treatment of chronic phase (CP) chronic myeloid leukemia (CML). In contrast, management of blast crisis (BC) CML remains challenging because BC cells acquire complex molecular alterations that confer stemness features to progenitor populations and resistance to BCR-ABL1 tyrosine kinase inhibitors. Comprehensive models of BC transformation have proved elusive because of the rarity and genetic heterogeneity of BC, but are important for developing biomarkers predicting BC progression and effective therapies. To better understand BC, we performed an integrated multiomics analysis of 74 CP and BC samples using whole-genome and exome sequencing, transcriptome and methylome profiling, and chromatin immunoprecipitation followed by high-throughput sequencing. Employing pathway-based analysis, we found the BC genome was significantly enriched for mutations affecting components of the polycomb repressive complex (PRC) pathway. While transcriptomically, BC progenitors were enriched and depleted for PRC1- and PRC2-related gene sets respectively. By integrating our data sets, we determined that BC progenitors undergo PRC-driven epigenetic reprogramming toward a convergent transcriptomic state. Specifically, PRC2 directs BC DNA hypermethylation, which in turn silences key genes involved in myeloid differentiation and tumor suppressor function via so-called epigenetic switching, whereas PRC1 represses an overlapping and distinct set of genes, including novel BC tumor suppressors. On the basis of these observations, we developed an integrated model of BC that facilitated the identification of combinatorial therapies capable of reversing BC reprogramming (decitabine+PRC1 inhibitors), novel PRC-silenced tumor suppressor genes (NR4A2), and gene expression signatures predictive of disease progression and drug resistance in CP.

Observations that suggest a contribution of altered dermal papilla mitochondrial function to androgenetic alopecia.

Androgenetic alopecia (AGA) is a prevalent hair loss condition in males that develops due to the influence of androgens and genetic predisposition. With the aim of elucidating genes involved in AGA pathogenesis, we modelled AGA with three-dimensional culture of keratinocyte-surrounded dermal papilla (DP) cells. We co-cultured immortalised balding and non-balding human DP cells (DPCs) derived from male AGA patients with epidermal keratinocyte (NHEK) using multi-interfacial polyelectrolyte complexation technique. We observed up-regulated mitochondria-related gene expression in balding compared with non-balding DP aggregates which indicated altered mitochondria metabolism. Further observation of significantly reduced electron transport chain complex activity (complexes I, IV and V), ATP levels and ability to uptake metabolites for ATP generation demonstrated compromised mitochondria function in balding DPC. Balding DP was also found to be under significantly higher oxidative stress than non-balding DP. Our experiments suggest that application of antioxidants lowers oxidative stress levels and improves metabolite uptake in balding DPC. We postulate that the observed up-regulation of mitochondria-related genes in balding DP aggregates resulted from an over-compensatory effort to rescue decreased mitochondrial function in balding DP through the attempted production of new functional mitochondria. In all, our three-dimensional co-culturing revealed mitochondrial dysfunction in balding DPC, suggesting a metabolic component in the aetiology of AGA.

Loss of the SWI/SNF-ATPase subunit members SMARCF1 (ARID1A), SMARCA2 (BRM), SMARCA4 (BRG1) and SMARCB1 (INI1) in oesophageal adenocarcinoma.

BACKGROUND: The SWI/SNF complex is an important chromatin remodeler, commonly dysregulated in cancer, with an estimated mutation frequency of 20%. ARID1A is the most frequently mutated subunit gene. Almost nothing is known about the other familiar members of the SWI/SNF complexes, SMARCA2 (BRM), SMARCA4 (BRG1) and SMARCB1 (INI1), in oesophageal adenocarcinoma (EAC). METHODS: We analysed a large cohort of 685 patients with EAC. We used four different antibodies to detect a loss-of-protein of ARID1A BRM, BRG1 and INI1 by immunohistochemistry and correlated these findings with molecular and clinical data. RESULTS: Loss of ARID1A, BRG1, BRM and INI1 was observed in 10.4, 3.4, 9.9 and 2% of EAC. We found a co-existing protein loss of ARID1A and BRM in 9.9% and of ARID1A and BRG1 in 2.2%. Patients with loss of ARID1A and TP53 wildtype EACs showed a shortened overall survival compared with AIRDA1A-positive tumours [median overall survival was 60.1 months (95%CI 1.2-139.9 months)] in patients with ARIDA-1A expression and 26.2 months (95%CI 3.7-19.1 months) in cases of ARIDA-1A loss (p = 0.044). Tumours with loss or expression of ARID1A and TP53 loss were not associated with a difference in survival. Only one tumour revealed high microsatellite instability (MSI-H) with concomitant ARID1A loss. All other ARID1A loss-EACs were microsatellite-stable (MSS). No predictive relevance was seen for SWI/SNF-complex alterations and simultaneous amplification of different genes (PIK3CA, KRAS, c-MYC, MET, GATA6, ERBB2). CONCLUSION: Our work describes, for the first time, loss of one of the SWI/SNF ATPase subunit proteins in a large number of adenocarcinomas of the oesophagus. Several papers discuss possible therapeutic interventions for tumours showing a loss of function of the SWI/SNF complex, such as PARP inhibitors or PI3K and AKT inhibitors. Future studies will be needed to show whether SWI/SNF complex-deficient EACs may benefit from personalized therapy.

Amplification of KRAS and its heterogeneity in non-Asian gastric adenocarcinomas.

BACKGROUND: Gastric cancer is one of the deadliest cancer entities worldwide. While surgery is the only curative treatment option in early tumors, for locally advanced and metastatic patients further therapeutic targets are needed. Several studies not only reported mutations but also amplifications of the KRAS locus in different cancer entities. More recently, KRAS amplification was discussed as a new therapeutic target. Little is known about the (prognostic) relevance and (heterogenic) distribution of KRAS amplification in gastric adenocarcinomas, especially in Non-Asian patients. METHODS: Amplification of the KRAS locus and corresponding protein expression was analyzed in 582 gastric adenocarcinomas employing fluorescence in-situ hybridization (FISH) and immunohistochemistry. Amplification status was correlated with clinico-pathological features, clinical outcome and molecular tumor data including a correlation to the TCGA subtypes of gastric carcinoma. RESULTS: KRAS amplification was detected in 27 out of 470 analysable tumors (5.7%) and correlated with protein expression of KRAS in all amplified tumors. Within the KRAS amplified gastric tumors 14/27 (51.9%) showed a heterogeneous distribution with also KRAS non-amplified tumor parts. According to TCGA 24 tumors (88.8%) were related to chromosomal instable tumors (CIN). The survival analysis of the entire patient cohort did not show any difference in overall survival in dependence on the KRAS status. However, a significant survival difference with a worse outcome for patients with KRAS amplified tumors was identified when analysing patients without neoadjuvant pre-treatment. CONCLUSIONS: We confirm the unfavorable prognosis of KRAS amplified tumors reported by other studies in (Asian) patient groups, at least in patients without neoadjuvant pre-treatment. Within KRAS amplified tumors we revealed intratumoral heterogeneity that may define a (more aggressive) tumor cell population which is more frequently observed in patients with lymph node metastases. Despite the heterogeneous distribution of KRAS amplified tumor clones, KRAS amplified locally advanced or metastasized gastric adenocarcinomas represent a therapeutically highly relevant tumor subgroup.

Personalized genome sequencing coupled with iPSC technology identifies GTDC1 as a gene involved in neurodevelopmental disorders.

The cellular and molecular mechanisms underlying neurodevelopmental conditions such as autism spectrum disorders have been studied intensively for decades. The ability to generate patient-specific induced pluripotent stem cells (iPSCs) now offers a novel strategy for modelling human diseases. Recent studies have reported the derivation of iPSCs from patients with neurological disorders. The key challenge remains the demonstration of disease-related phenotypes and the ability to model the disease. Here we report a case study with signs of neurodevelopmental disorders (NDDs) harbouring chromosomal rearrangements that were sequenced using long-insert DNA paired-end tag (DNA-PET) sequencing approach. We identified the disruption of a specific gene, GTDC1. By deriving iPSCs from this patient and differentiating them into neural progenitor cells (NPCs) and neurons we dissected the disease process at the cellular level and observed defects in both NPCs and neuronal cells. We also showed that disruption of GTDC1 expression in wild type human NPCs and neurons showed a similar phenotype as patient's iPSCs. Finally, we utilized a zebrafish model to demonstrate a role for GTDC1 in the development of the central nervous system. Our findings highlight the importance of combining sequencing technologies with the iPSC technology for NDDs modelling that could be applied for personalized medicine.

ALPK1 missense pathogenic variant in five families leads to ROSAH syndrome, an ocular multisystem autosomal dominant disorder.

PURPOSE: To identify the molecular cause in five unrelated families with a distinct autosomal dominant ocular systemic disorder we called ROSAH syndrome due to clinical features of retinal dystrophy, optic nerve edema, splenomegaly, anhidrosis, and migraine headache. METHODS: Independent discovery exome and genome sequencing in families 1, 2, and 3, and confirmation in families 4 and 5. Expression of wild-type messenger RNA and protein in human and mouse tissues and cell lines. Ciliary assays in fibroblasts from affected and unaffected family members. RESULTS: We found the heterozygous missense variant in the ɑ-kinase gene, ALPK1, (c.710C>T, [p.Thr237Met]), segregated with disease in all five families. All patients shared the ROSAH phenotype with additional low-grade ocular inflammation, pancytopenia, recurrent infections, and mild renal impairment in some. ALPK1 was notably expressed in retina, retinal pigment epithelium, and optic nerve, with immunofluorescence indicating localization to the basal body of the connecting cilium of the photoreceptors, and presence in the sweat glands. Immunocytofluorescence revealed expression at the centrioles and spindle poles during metaphase, and at the base of the primary cilium. Affected family member fibroblasts demonstrated defective ciliogenesis. CONCLUSION: Heterozygosity for ALPK1, p.Thr237Met leads to ROSAH syndrome, an autosomal dominant ocular systemic disorder.

Concurrent Single-Cell RNA and Targeted DNA Sequencing on an Automated Platform for Comeasurement of Genomic and Transcriptomic Signatures.

BACKGROUND: The comeasurement of both genomic and transcriptomic signatures in single cells is of fundamental importance to accurately assess how the genetic information correlates with the transcriptomic phenotype. However, existing technologies have low throughput and laborious work flows. METHODS: We developed a new method for concurrent sequencing of the transcriptome and targeted genomic regions (CORTAD-seq) within the same single cell on an automated microfluidic platform. The method was compatible with the downstream library preparation, allowing easy integration into existing next-generation sequencing work flows. We incorporated a single-cell bioinformatics pipeline for transcriptome and mutation analysis. RESULTS: As proof of principle, we applied CORTAD-seq to lung cancer cell lines to dissect the cellular consequences of mutations that result in resistance to targeted therapy. We obtained a mean detection of 6000 expressed genes and an exonic rate of 50%. The targeted DNA-sequencing data achieved a 97.8% detection rate for mutations and allowed for the identification of copy number variations and haplotype construction. We detected expression signatures of tyrosine kinase inhibitor (TKI) resistance, epidermal growth factor receptor (EGFR) amplification, and expansion of the T790M mutation among resistant cells. We also identified characteristics for TKI resistance that were independent of EGFR T790M, indicating that other alterations are required for resistance in this context. CONCLUSIONS: CORTAD-seq allows assessment of the interconnection between genetic and transcriptomic changes in single cells. It is operated on an automated, commercially available single-cell isolation platform, making its implementation straightforward.

Systems consequences of amplicon formation in human breast cancer.

Chromosomal structural variations play an important role in determining the transcriptional landscape of human breast cancers. To assess the nature of these structural variations, we analyzed eight breast tumor samples with a focus on regions of gene amplification using mate-pair sequencing of long-insert genomic DNA with matched transcriptome profiling. We found that tandem duplications appear to be early events in tumor evolution, especially in the genesis of amplicons. In a detailed reconstruction of events on chromosome 17, we found large unpaired inversions and deletions connect a tandemly duplicated ERBB2 with neighboring 17q21.3 amplicons while simultaneously deleting the intervening BRCA1 tumor suppressor locus. This series of events appeared to be unusually common when examined in larger genomic data sets of breast cancers albeit using approaches with lesser resolution. Using siRNAs in breast cancer cell lines, we showed that the 17q21.3 amplicon harbored a significant number of weak oncogenes that appeared consistently coamplified in primary tumors. Down-regulation of BRCA1 expression augmented the cell proliferation in ERBB2-transfected human normal mammary epithelial cells. Coamplification of other functionally tested oncogenic elements in other breast tumors examined, such as RIPK2 and MYC on chromosome 8, also parallel these findings. Our analyses suggest that structural variations efficiently orchestrate the gain and loss of cancer gene cassettes that engage many oncogenic pathways simultaneously and that such oncogenic cassettes are favored during the evolution of a cancer.

Novel therapeutic targets on the horizon for lung cancer.

Lung cancer is a leading cause of cancer-related mortality worldwide, and is classically divided into two major histological subtypes: non-small-cell lung cancer (NSCLC) and small-cell lung cancer (SCLC). Although NSCLC and SCLC are considered distinct entities with different genomic landscapes, emerging evidence highlights a convergence in therapeutically relevant targets for both histologies. In adenocarcinomas with defined alterations such as EGFR mutations and ALK translocations, targeted therapies are now first-line standard of care. By contrast, many experimental and targeted agents remain largely unsuccessful for SCLC. Intense preclinical research and clinical trials are underway to exploit unique traits of lung cancer, such as oncogene dependency, DNA damage response, angiogenesis, and cellular plasticity arising from presence of cancer stem cell lineages. In addition, the promising clinical activity observed in NSCLC in response to immune checkpoint blockade has spurred great interest in the field of immunooncology, with the scope to develop a diverse repertoire of synergistic and personalised immunotherapeutics. In this Review, we discuss novel therapeutic agents for lung cancer that are in early-stage development, and how prospective clinical trials and drug development may be shaped by a deeper understanding of this heterogeneous disease.

Six novel susceptibility Loci for early-onset androgenetic alopecia and their unexpected association with common diseases.

Androgenetic alopecia (AGA) is a highly heritable condition and the most common form of hair loss in humans. Susceptibility loci have been described on the X chromosome and chromosome 20, but these loci explain a minority of its heritable variance. We conducted a large-scale meta-analysis of seven genome-wide association studies for early-onset AGA in 12,806 individuals of European ancestry. While replicating the two AGA loci on the X chromosome and chromosome 20, six novel susceptibility loci reached genome-wide significance (p = 2.62×10⁻9-1.01×10⁻¹²). Unexpectedly, we identified a risk allele at 17q21.31 that was recently associated with Parkinson's disease (PD) at a genome-wide significant level. We then tested the association between early-onset AGA and the risk of PD in a cross-sectional analysis of 568 PD cases and 7,664 controls. Early-onset AGA cases had significantly increased odds of subsequent PD (OR = 1.28, 95% confidence interval: 1.06-1.55, p = 8.9×10⁻³). Further, the AGA susceptibility alleles at the 17q21.31 locus are on the H1 haplotype, which is under negative selection in Europeans and has been linked to decreased fertility. Combining the risk alleles of six novel and two established susceptibility loci, we created a genotype risk score and tested its association with AGA in an additional sample. Individuals in the highest risk quartile of a genotype score had an approximately six-fold increased risk of early-onset AGA [odds ratio (OR) = 5.78, p = 1.4×10⁻88]. Our results highlight unexpected associations between early-onset AGA, Parkinson's disease, and decreased fertility, providing important insights into the pathophysiology of these conditions.

Familial hypomagnesemia with hypercalciuria and nephrocalcinosis: variable phenotypic expression in three affected sisters from Mexican ancestry.

Familial hypomagnesemia with hypercalciuria and nephrocalcinosis is a rare autosomal recessive renal disease caused by mutations in genes for the tight junction transmembrane proteins Claudin-16 (CLDN16) and Claudin-19 (CLDN19). We present the first case report of a Mexican family with three affected sisters carrying a p.Gly20Asp mutation in CLDN19 whose heterozygous mother showed evident hypercalciuria and normal low magnesemia without any other clinical, laboratory, and radiological symptoms of renal disease making of her an unsuitable donor. The affected sisters showed variable phenotypic expression including age of first symptoms, renal urinary tract infections, nephrolithiasis, nephrocalcinosis, and eye symptoms consisting in retinochoroiditis, strabismus, macular scars, bilateral anisocoria, and severe myopia and astigmatism. End stage renal disease due to renal failure needed kidney transplantation in the three of them. Interesting findings were a heterozygous mother with asymptomatic hypercalciuria warning on the need of carefully explore clinical, laboratory, kidney ultrasonograpy, and mutation status in first degree asymptomatic relatives to avoid inappropriate kidney donors; an evident variable phenotypic expression among patients; the identification of a mutation almost confined to Spanish cases and a 3.5 Mb block of genomic homozygosis strongly suggesting a common remote parental ancestor for the gene mutation reported.

Hereditary parkinsonism with dementia is caused by mutations in ATP13A2, encoding a lysosomal type 5 P-type ATPase.

Neurodegenerative disorders such as Parkinson and Alzheimer disease cause motor and cognitive dysfunction and belong to a heterogeneous group of common and disabling disorders. Although the complex molecular pathophysiology of neurodegeneration is largely unknown, major advances have been achieved by elucidating the genetic defects underlying mendelian forms of these diseases. This has led to the discovery of common pathophysiological pathways such as enhanced oxidative stress, protein misfolding and aggregation and dysfunction of the ubiquitin-proteasome system. Here, we describe loss-of-function mutations in a previously uncharacterized, predominantly neuronal P-type ATPase gene, ATP13A2, underlying an autosomal recessive form of early-onset parkinsonism with pyramidal degeneration and dementia (PARK9, Kufor-Rakeb syndrome). Whereas the wild-type protein was located in the lysosome of transiently transfected cells, the unstable truncated mutants were retained in the endoplasmic reticulum and degraded by the proteasome. Our findings link a class of proteins with unknown function and substrate specificity to the protein networks implicated in neurodegeneration and parkinsonism.

Impaired development of neural-crest cell-derived organs and intellectual disability caused by MED13L haploinsufficiency.

MED13L is a component subunit of the Mediator complex, an important regulator of transcription that is highly conserved across eukaryotes. Here, we report MED13L disruption in a translocation t(12;19) breakpoint of a patient with Pierre-Robin syndrome, moderate intellectual disability, craniofacial anomalies, and muscular defects. The phenotype is similar to previously described patients with MED13L haploinsufficiency. Knockdown of MED13L orthologue in zebrafish, med13b, showed early defective migration of cranial neural crest cells (NCCs) that contributed to cartilage structure deformities in the later stage, recapitulating craniofacial anomalies seen in human patients. Notably, we observed abnormal distribution of developing neurons in different brain regions of med13b morphant embryos, which could be rescued upon introduction of full-length human MED13L mRNA. To compare with mammalian system, we suppressed MED13L expression by short-hairpin RNA in ES-derived human neural progenitors, and differentiated them into neurons. Transcriptome analysis revealed differential expression of components of Wnt and FGF signaling pathways in MED13L-deficient neurons. Our finding provides a novel insight into the mechanism of overlapping phenotypic outcome targeting NCCs derivatives organs in patients with MED13L haploinsufficiency, and emphasizes a clinically recognizable syndromic phenotype in these patients.

Functional chromatin features are associated with structural mutations in cancer.

BACKGROUND: Structural mutations (SMs) play a major role in cancer development. In some cancers, such as breast and ovarian, DNA double-strand breaks (DSBs) occur more frequently in transcribed regions, while in other cancer types such as prostate, there is a consistent depletion of breakpoints in transcribed regions. Despite such regularity, little is understood about the mechanisms driving these effects. A few works have suggested that protein binding may be relevant, e.g. in studies of androgen receptor binding and active chromatin in specific cell types. We hypothesized that this behavior might be general, i.e. that correlation between protein-DNA binding (and open chromatin) and breakpoint locations is common across divergent cancers. RESULTS: We investigated this hypothesis by comprehensively analyzing the relationship among 457 ENCODE protein binding ChIP-seq experiments, 125 DnaseI and 24 FAIRE experiments, and 14,600 SMs from 8 diverse cancer datasets covering 147 samples. In most cancers, including breast and ovarian, we found enrichment of protein binding and open chromatin in the vicinity of SM breakpoints at distances up to 200 kb. Furthermore, for all cancer types we observed an enhanced enrichment in regions distant from genes when compared to regions proximal to genes, suggesting that the SM-induction mechanism is independent from the bias of DSBs to occur near transcribed regions. We also observed a stronger effect for sites with more than one protein bound. CONCLUSIONS: Protein binding and open chromatin state are associated with nearby SM breakpoints in many cancer datasets. These observations suggest a consistent mechanism underlying SM locations across different cancers.

Transcriptional consequences of genomic structural aberrations in breast cancer.

Using a long-span, paired-end deep sequencing strategy, we have comprehensively identified cancer genome rearrangements in eight breast cancer genomes. Herein, we show that 40%-54% of these structural genomic rearrangements result in different forms of fusion transcripts and that 44% are potentially translated. We find that single segmental tandem duplication spanning several genes is a major source of the fusion gene transcripts in both cell lines and primary tumors involving adjacent genes placed in the reverse-order position by the duplication event. Certain other structural mutations, however, tend to attenuate gene expression. From these candidate gene fusions, we have found a fusion transcript (RPS6KB1-VMP1) recurrently expressed in ∼30% of breast cancers associated with potential clinical consequences. This gene fusion is caused by tandem duplication on 17q23 and appears to be an indicator of local genomic instability altering the expression of oncogenic components such as MIR21 and RPS6KB1.

Comprehensive long-span paired-end-tag mapping reveals characteristic patterns of structural variations in epithelial cancer genomes.

Somatic genome rearrangements are thought to play important roles in cancer development. We optimized a long-span paired-end-tag (PET) sequencing approach using 10-Kb genomic DNA inserts to study human genome structural variations (SVs). The use of a 10-Kb insert size allows the identification of breakpoints within repetitive or homology-containing regions of a few kilobases in size and results in a higher physical coverage compared with small insert libraries with the same sequencing effort. We have applied this approach to comprehensively characterize the SVs of 15 cancer and two noncancer genomes and used a filtering approach to strongly enrich for somatic SVs in the cancer genomes. Our analyses revealed that most inversions, deletions, and insertions are germ-line SVs, whereas tandem duplications, unpaired inversions, interchromosomal translocations, and complex rearrangements are over-represented among somatic rearrangements in cancer genomes. We demonstrate that the quantitative and connective nature of DNA-PET data is precise in delineating the genealogy of complex rearrangement events, we observe signatures that are compatible with breakage-fusion-bridge cycles, and we discover that large duplications are among the initial rearrangements that trigger genome instability for extensive amplification in epithelial cancers.

Frequency and functional characterization of fusion genes in squamous cell carcinoma of the lung.

INTRODUCTION: In contrast to lung adenocarcinoma (LUAD), targetable genetic alterations are less frequently detected in squamous cell carcinoma of the lung (LUSC). Over the last years, gene fusions have become promising targets in many solid cancers. Here, we analysed a cohort of LUSC, identified recurrent fusion genes and functionally characterised these tumour genomes. METHODS: A subset of 1608 squamous cell carcinomas of the lung was analysed by means of the FusionPlex® Lung Panel to identify potentially targetable gene fusions using targeted next-generation sequencing. Cases harbouring recurrent gene fusions were further analysed using FISH, Cytoscan HD arrays and cell culture experiments. RESULTS: We found both, known and novel gene fusions in about 3 % of the cases. Known fusions occurring in lung cancer included ALK::EML4, EGFRvIII, EZR::ROS1 and FGFR3::TACC. We further identified recurrent gene fusions of currently unknown biological function, involving EGFR::VSTM2A and NSD3::FGFR1 and showed that the occurrence of the EGFR::VSTM2A fusion is accompanied by high-level amplification of EGFR. Our analyses further revealed that the genomes of these LUSC patients are chromosomally unstable, which leads us to believe that such non-actionable genomic rearrangements may be a result of "chromosomal chaos" most probably not representing exclusive cancer-driving genes in this cancer entity. CONCLUSIONS: We emphasise that caution should be taken when novel fusions are found and that the appearance of new gene fusions should always be interpreted in the molecular context of the respective disease.

Tumor phylogeography reveals block-shaped spatial heterogeneity and the mode of evolution in Hepatocellular Carcinoma.

Solid tumors are complex ecosystems with heterogeneous 3D structures, but the spatial intra-tumor heterogeneity (sITH) at the macroscopic (i.e., whole tumor) level is under-explored. Using a phylogeographic approach, we sequence genomes and transcriptomes from 235 spatially informed sectors across 13 hepatocellular carcinomas (HCC), generating one of the largest datasets for studying sITH. We find that tumor heterogeneity in HCC segregates into spatially variegated blocks with large genotypic and phenotypic differences. By dissecting the transcriptomic heterogeneity, we discover that 30% of patients had a "spatially competing distribution" (SCD), where different spatial blocks have distinct transcriptomic subtypes co-existing within a tumor, capturing the critical transition period in disease progression. Interestingly, the tumor regions with more advanced transcriptomic subtypes (e.g., higher cell cycle) often take clonal dominance with a wider geographic range, rejecting neutral evolution for SCD patients. Extending the statistical tests for detecting natural selection to many non-SCD patients reveal varying levels of selective signal across different tumors, implying that many evolutionary forces including natural selection and geographic isolation can influence the overall pattern of sITH. Taken together, tumor phylogeography unravels a dynamic landscape of sITH, pinpointing important evolutionary and clinical consequences of spatial heterogeneity in cancer.