Ultraviolet radiation shapes dendritic cell leukaemia transformation in the skin.

Tumours most often arise from progression of precursor clones within a single anatomical niche. In the bone marrow, clonal progenitors can undergo malignant transformation to acute leukaemia, or differentiate into immune cells that contribute to disease pathology in peripheral tissues1-4. Outside the marrow, these clones are potentially exposed to a variety of tissue-specific mutational processes, although the consequences of this are unclear. Here we investigate the development of blastic plasmacytoid dendritic cell neoplasm (BPDCN)-an unusual form of acute leukaemia that often presents with malignant cells isolated to the skin5. Using tumour phylogenomics and single-cell transcriptomics with genotyping, we find that BPDCN arises from clonal (premalignant) haematopoietic precursors in the bone marrow. We observe that BPDCN skin tumours first develop at sun-exposed anatomical sites and are distinguished by clonally expanded mutations induced by ultraviolet (UV) radiation. A reconstruction of tumour phylogenies reveals that UV damage can precede the acquisition of alterations associated with malignant transformation, implicating sun exposure of plasmacytoid dendritic cells or committed precursors during BPDCN pathogenesis. Functionally, we find that loss-of-function mutations in Tet2, the most common premalignant alteration in BPDCN, confer resistance to UV-induced cell death in plasmacytoid, but not conventional, dendritic cells, suggesting a context-dependent tumour-suppressive role for TET2. These findings demonstrate how tissue-specific environmental exposures at distant anatomical sites can shape the evolution of premalignant clones to disseminated cancer.

Sex-Biased ZRSR2 Mutations in Myeloid Malignancies Impair Plasmacytoid Dendritic Cell Activation and Apoptosis.

Blastic plasmacytoid dendritic cell neoplasm (BPDCN) is an aggressive leukemia of plasmacytoid dendritic cells (pDC). BPDCN occurs at least three times more frequently in men than in women, but the reasons for this sex bias are unknown. Here, studying genomics of primary BPDCN and modeling disease-associated mutations, we link acquired alterations in RNA splicing to abnormal pDC development and inflammatory response through Toll-like receptors. Loss-of-function mutations in ZRSR2, an X chromosome gene encoding a splicing factor, are enriched in BPDCN, and nearly all mutations occur in males. ZRSR2 mutation impairs pDC activation and apoptosis after inflammatory stimuli, associated with intron retention and inability to upregulate the transcription factor IRF7. In vivo, BPDCN-associated mutations promote pDC expansion and signatures of decreased activation. These data support a model in which male-biased mutations in hematopoietic progenitors alter pDC function and confer protection from apoptosis, which may impair immunity and predispose to leukemic transformation. SIGNIFICANCE: Sex bias in cancer is well recognized, but the underlying mechanisms are incompletely defined. We connect X chromosome mutations in ZRSR2 to an extremely male-predominant leukemia. Aberrant RNA splicing induced by ZRSR2 mutation impairs dendritic cell inflammatory signaling, interferon production, and apoptosis, revealing a sex- and lineage-related tumor suppressor pathway.This article is highlighted in the In This Issue feature, p. 275.

Short loop functional commonality identified in leukaemia proteome highlights crucial protein sub-networks.

Direct drug targeting of mutated proteins in cancer is not always possible and efficacy can be nullified by compensating protein-protein interactions (PPIs). Here, we establish an in silico pipeline to identify specific PPI sub-networks containing mutated proteins as potential targets, which we apply to mutation data of four different leukaemias. Our method is based on extracting cyclic interactions of a small number of proteins topologically and functionally linked in the Protein-Protein Interaction Network (PPIN), which we call short loop network motifs (SLM). We uncover a new property of PPINs named 'short loop commonality' to measure indirect PPIs occurring via common SLM interactions. This detects 'modules' of PPI networks enriched with annotated biological functions of proteins containing mutation hotspots, exemplified by FLT3 and other receptor tyrosine kinase proteins. We further identify functional dependency or mutual exclusivity of short loop commonality pairs in large-scale cellular CRISPR-Cas9 knockout screening data. Our pipeline provides a new strategy for identifying new therapeutic targets for drug discovery.

Comprehensive metagenomic analysis of blastic plasmacytoid dendritic cell neoplasm.

Blastic plasmacytoid dendritic cell neoplasm (BPDCN) is a hematologic malignancy believed to originate from plasmacytoid dendritic cells (pDCs), the immune cells responsible for producing type 1 interferons during infection. Nearly all patients with BPDCN have prominent skin involvement, with cutaneous infiltration occupying the dermis and subcutis. One half of patients present with BPDCN cells only in the skin, with no evidence of disease elsewhere. Because normal pDCs are rare or absent in cutaneous sites, and they only traffic to the skin after activation by pathogen or inflammation, our aim was to determine if a microorganism is associated with BPDCN. We performed RNA sequencing in BPDCN skin and bone marrow, with cutaneous T-cell lymphoma (CTCL) and normal skin as controls. GATK-PathSeq was used to identify known microbial sequences. Bacterial reads in BPDCN skin were components of normal flora and did not distinguish BPDCN from controls. We then developed a new computational tool, virID (Viral Identification and Discovery; https://github.com/jnoms/virID), for identification of microbial-associated reads remaining unassigned after GATK-PathSeq. We found no evidence for a known or novel virus in BPDCN skin or bone marrow, despite confirming that virID could identify Merkel cell polyomavirus in Merkel cell carcinoma, human papillomavirus in head and neck squamous cell carcinoma, and Kaposi's sarcoma herpesvirus in Kaposi's sarcoma in a blinded fashion. Thus, at the level of sensitivity used here, we found no clear pathogen linked to BPDCN.

Genetic variants and protein-protein interactions: a multidimensional network-centric view.

We review recent progress in the mapping of genetic variants to proteins, in the context of their interactions, as measured from experiments and/or computational predictions. Such variants can impact on the molecular mechanisms underlying an interaction and its stability. We highlight recent work which relies on the effective use of protein-protein interaction networks (PPINs), integrated with 3D structural information, for evaluating disease-associated variants. Furthermore, we discuss how the integration of multiple layers of biological information, in the context of PPINs, can improve the interpretation of genetic variants and inspire new therapeutic strategies.

Deep phenotyping of Tregs identifies an immune signature for idiopathic aplastic anemia and predicts response to treatment.

Idiopathic aplastic anemia (AA) is an immune-mediated and serious form of bone marrow failure. Akin to other autoimmune diseases, we have previously shown that in AA regulatory T cells (Tregs) are reduced in number and function. The aim of this study was to further characterize Treg subpopulations in AA and investigate the potential correlation between specific Treg subsets and response to immunosuppressive therapy (IST) as well as their in vitro expandability for potential clinical use. Using mass cytometry and an unbiased multidimensional analytical approach, we identified 2 specific human Treg subpopulations (Treg A and Treg B) with distinct phenotypes, gene expression, expandability, and function. Treg B predominates in IST responder patients, has a memory/activated phenotype (with higher expression of CD95, CCR4, and CD45RO within FOXP3(hi), CD127(lo) Tregs), expresses the interleukin-2 (IL-2)/STAT5 pathway and cell-cycle commitment genes. Furthermore, in vitro-expanded Tregs become functional and take on the characteristics of Treg B. Collectively, this study identifies human Treg subpopulations that can be used as predictive biomarkers for response to IST in AA and potentially other autoimmune diseases. We also show that Tregs from AA patients are IL-2-sensitive and expandable in vitro, suggesting novel therapeutic approaches such as low-dose IL-2 therapy and/or expanded autologous Tregs and meriting further exploration.

Anatomy of protein disorder, flexibility and disease-related mutations.

Integration of protein structural information with human genetic variation and pathogenic mutations is essential to understand molecular mechanisms associated with the effects of polymorphisms on protein interactions and cellular processes. We investigate occurrences of non-synonymous SNPs in ordered and disordered protein regions by systematic mapping of common variants and disease-related SNPs onto these regions. We show that common variants accumulate in disordered regions; conversely pathogenic variants are significantly depleted in disordered regions. These different occurrences of pathogenic and common SNPs can be attributed to a negative selection on random mutations in structurally highly constrained regions. New approaches in the study of quantitative effects of pathogenic-related mutations should effectively account for all the possible contexts and relative functional constraints in which the sequence variation occurs.

Bridging topological and functional information in protein interaction networks by short loops profiling.

Protein-protein interaction networks (PPINs) have been employed to identify potential novel interconnections between proteins as well as crucial cellular functions. In this study we identify fundamental principles of PPIN topologies by analysing network motifs of short loops, which are small cyclic interactions of between 3 and 6 proteins. We compared 30 PPINs with corresponding randomised null models and examined the occurrence of common biological functions in loops extracted from a cross-validated high-confidence dataset of 622 human protein complexes. We demonstrate that loops are an intrinsic feature of PPINs and that specific cell functions are predominantly performed by loops of different lengths. Topologically, we find that loops are strongly related to the accuracy of PPINs and define a core of interactions with high resilience. The identification of this core and the analysis of loop composition are promising tools to assess PPIN quality and to uncover possible biases from experimental detection methods. More than 96% of loops share at least one biological function, with enrichment of cellular functions related to mRNA metabolic processing and the cell cycle. Our analyses suggest that these motifs can be used in the design of targeted experiments for functional phenotype detection.