Multiomics profiling of primary lung cancers and distant metastases reveals immunosuppression as a common characteristic of tumor cells with metastatic plasticity.

BACKGROUND: Metastasis is the primary cause of cancer mortality accounting for 90% of cancer deaths. Our understanding of the molecular mechanisms driving metastasis is rudimentary. RESULTS: We perform whole exome sequencing (WES), RNA sequencing, methylation microarray, and immunohistochemistry (IHC) on 8 pairs of non-small cell lung cancer (NSCLC) primary tumors and matched distant metastases. Furthermore, we analyze published WES data from 35 primary NSCLC and metastasis pairs, and transcriptomic data from 4 autopsy cases with metastatic NSCLC and one metastatic lung cancer mouse model. The majority of somatic mutations are shared between primary tumors and paired distant metastases although mutational signatures suggest different mutagenesis processes in play before and after metastatic spread. Subclonal analysis reveals evidence of monoclonal seeding in 41 of 42 patients. Pathway analysis of transcriptomic data reveals that downregulated pathways in metastases are mainly immune-related. Further deconvolution analysis reveals significantly lower infiltration of various immune cell types in metastases with the exception of CD4+ T cells and M2 macrophages. These results are in line with lower densities of immune cells and higher CD4/CD8 ratios in metastases shown by IHC. Analysis of transcriptomic data from autopsy cases and animal models confirms that immunosuppression is also present in extracranial metastases. Significantly higher somatic copy number aberration and allelic imbalance burdens are identified in metastases. CONCLUSIONS: Metastasis is a molecularly late event, and immunosuppression driven by different molecular events, including somatic copy number aberration, may be a common characteristic of tumors with metastatic plasticity.

Integrative genomic analysis of adult mixed phenotype acute leukemia delineates lineage associated molecular subtypes.

Mixed phenotype acute leukemia (MPAL) is a rare subtype of acute leukemia characterized by leukemic blasts presenting myeloid and lymphoid markers. Here we report data from integrated genomic analysis on 31 MPAL samples and compare molecular profiling with that from acute myeloid leukemia (AML), B cell acute lymphoblastic leukemia (B-ALL), and T cell acute lymphoblastic leukemia (T-ALL). Consistent with the mixed immunophenotype, both AML-type and ALL-type mutations are detected in MPAL. Myeloid-B and myeloid-T MPAL show distinct mutation and methylation signatures that are associated with differences in lineage-commitment gene expressions. Genome-wide methylation comparison among MPAL, AML, B-ALL, and T-ALL sub-classifies MPAL into AML-type and ALL-type MPAL, which is associated with better clinical response when lineage-matched therapy is given. These results elucidate the genetic and epigenetic heterogeneity of MPAL and its genetic distinction from AML, B-ALL, and T-ALL and further provide proof of concept for a molecularly guided precision therapy approach in MPAL.

Clearance of Somatic Mutations at Remission and the Risk of Relapse in Acute Myeloid Leukemia.

Purpose The aim of the current study was to determine whether the degree of mutation clearance at remission predicts the risk of relapse in patients with acute myeloid leukemia (AML). Patients and Methods One hundred thirty-one previously untreated patients with AML who received intensive induction chemotherapy and attained morphologic complete remission (CR) at day 30 were studied. Pretreatment and CR bone marrow were analyzed using targeted capture DNA sequencing. We analyzed the association between mutation clearance (MC) on the basis of variant allele frequency (VAF) at CR (MC2.5: if the VAF of residual mutations was < 2.5%; MC1.0: if the VAF was < 1%; and complete MC [CMC]: if no detectable residual mutations) and event-free survival, overall survival (OS), and cumulative incidence of relapse (CIR). Results MC1.0 and CMC were associated with significantly better OS (2-year OS: 75% v 61% in MC1.0 v non-MC1.0; P = .0465; 2-year OS: 77% v 60% in CMC v non-CMC; P = .0303) and lower CIR (2-year CIR: 26% v 46% in MC1.0 v non-MC 1.0; P = .0349; 2 year-CIR: 24% v 46% in CMC v non-CMC; P = .03), whereas there was no significant difference in any of the above outcomes by MC2.5. Multivariable analysis adjusting for age, cytogenetic risk, allogeneic stem-cell transplantation, and flow cytometry-based minimal residual disease revealed that patients with CMC had significantly better event-free survival (hazard ratio [HR], 0.43; P = .0083), OS (HR, 0.47; P = .04), and CIR (HR, 0.27; P < .001) than did patients without CMC. These prognostic associations were stronger when preleukemic mutations, such as DNMT3A, TET2, and ASXL1, were removed from the analysis. Conclusion Clearance of somatic mutation at CR, particularly in nonpreleukemic genes, was associated with significantly better survival and less risk of relapse. Somatic mutations in nonpreleukemic genes may function as a molecular minimal residual disease marker in AML.

A massive expansion of effector genes underlies gall-formation in the wheat pest Mayetiola destructor.

Gall-forming arthropods are highly specialized herbivores that, in combination with their hosts, produce extended phenotypes with unique morphologies [1]. Many are economically important, and others have improved our understanding of ecology and adaptive radiation [2]. However, the mechanisms that these arthropods use to induce plant galls are poorly understood. We sequenced the genome of the Hessian fly (Mayetiola destructor; Diptera: Cecidomyiidae), a plant parasitic gall midge and a pest of wheat (Triticum spp.), with the aim of identifying genic modifications that contribute to its plant-parasitic lifestyle. Among several adaptive modifications, we discovered an expansive reservoir of potential effector proteins. Nearly 5% of the 20,163 predicted gene models matched putative effector gene transcripts present in the M. destructor larval salivary gland. Another 466 putative effectors were discovered among the genes that have no sequence similarities in other organisms. The largest known arthropod gene family (family SSGP-71) was also discovered within the effector reservoir. SSGP-71 proteins lack sequence homologies to other proteins, but their structures resemble both ubiquitin E3 ligases in plants and E3-ligase-mimicking effectors in plant pathogenic bacteria. SSGP-71 proteins and wheat Skp proteins interact in vivo. Mutations in different SSGP-71 genes avoid the effector-triggered immunity that is directed by the wheat resistance genes H6 and H9. Results point to effectors as the agents responsible for arthropod-induced plant gall formation.

The first myriapod genome sequence reveals conservative arthropod gene content and genome organisation in the centipede Strigamia maritima.

Myriapods (e.g., centipedes and millipedes) display a simple homonomous body plan relative to other arthropods. All members of the class are terrestrial, but they attained terrestriality independently of insects. Myriapoda is the only arthropod class not represented by a sequenced genome. We present an analysis of the genome of the centipede Strigamia maritima. It retains a compact genome that has undergone less gene loss and shuffling than previously sequenced arthropods, and many orthologues of genes conserved from the bilaterian ancestor that have been lost in insects. Our analysis locates many genes in conserved macro-synteny contexts, and many small-scale examples of gene clustering. We describe several examples where S. maritima shows different solutions from insects to similar problems. The insect olfactory receptor gene family is absent from S. maritima, and olfaction in air is likely effected by expansion of other receptor gene families. For some genes S. maritima has evolved paralogues to generate coding sequence diversity, where insects use alternate splicing. This is most striking for the Dscam gene, which in Drosophila generates more than 100,000 alternate splice forms, but in S. maritima is encoded by over 100 paralogues. We see an intriguing linkage between the absence of any known photosensory proteins in a blind organism and the additional absence of canonical circadian clock genes. The phylogenetic position of myriapods allows us to identify where in arthropod phylogeny several particular molecular mechanisms and traits emerged. For example, we conclude that juvenile hormone signalling evolved with the emergence of the exoskeleton in the arthropods and that RR-1 containing cuticle proteins evolved in the lineage leading to Mandibulata. We also identify when various gene expansions and losses occurred. The genome of S. maritima offers us a unique glimpse into the ancestral arthropod genome, while also displaying many adaptations to its specific life history.