SUMMARY While alternative splicing is known to diversify the functional characteristics of some genes, the extent to which protein isoforms globally contribute to functional complexity on a proteomic scale remains unknown. To address this systematically, we cloned full-length open reading frames of alternatively spliced transcripts for a large number of human genes and used protein-protein interaction profiling to functionally compare hundreds of protein isoform pairs. The majority of isoform pairs share less than 50% of their interactions. In the global context of interactome network maps, alternative isoforms tend to behave like distinct proteins rather than minor variants of each other. Interaction partners specific to alternative isoforms tend to be expressed in a highly tissue-specific manner and belong to distinct functional modules. Our strategy, applicable to other functional characteristics, reveals a widespread expansion of protein interaction capabilities through alternative splicing and suggests that many alternative ‘isoforms’ are functionally divergent (i.e., ‘functional alloforms’).

Genotypic differences greatly influence susceptibility and resistance to disease. Understanding genotype-phenotype relationships requires that phenotypes be viewed as manifestations of network properties, rather than simply as the result of individual genomic variations1. Genome sequencing efforts have identified numerous germline mutations associated with cancer predisposition and large numbers of somatic genomic alterations2. However, it remains challenging to distinguish between background, or “passenger” and causal, or “driver” cancer mutations in these datasets. Human viruses intrinsically depend on their host cell during the course of infection and can elicit pathological phenotypes similar to those arising from mutations3. To test the hypothesis that genomic variations and tumour viruses may cause cancer via related mechanisms, we systematically examined host interactome and transcriptome network perturbations caused by DNA tumour virus proteins. The resulting integrated viral perturbation data reflects rewiring of the host cell networks, and highlights pathways that go awry in cancer, such as Notch signalling and apoptosis. We show that systematic analyses of host targets of viral proteins can identify cancer genes with a success rate on par with their identification through functional genomics and large-scale cataloguing of tumour mutations. Together, these complementary approaches result in increased specificity for cancer gene identification. Combining systems-level studies of pathogen-encoded gene products with genomic approaches will facilitate prioritization of cancer-causing driver genes so as to advance understanding of the genetic basis of human cancer.

Crucial transitions in cancer-including tumor initiation, local expansion, metastasis, and therapeutic resistance-involve complex interactions between cells within the dynamic tumor ecosystem. Transformative single-cell genomics technologies and spatial multiplex in situ methods now provide an opportunity to interrogate this complexity at unprecedented resolution. The Human Tumor Atlas Network (HTAN), part of the National Cancer Institute (NCI) Cancer Moonshot Initiative, will establish a clinical, experimental, computational, and organizational framework to generate informative and accessible three-dimensional atlases of cancer transitions for a diverse set of tumor types. This effort complements both ongoing efforts to map healthy organs and previous largescale cancer genomics approaches focused on bulk sequencing at a single point in time. Generating single-cell, multiparametric, longitudinal atlases and integrating them with clinical outcomes should help identify novel predictive biomarkers and features as well as therapeutically relevant cell types, cell states, and cellular interactions across transitions. The resulting tumor atlases should have a profound impact on our understanding of cancer biology and have the potential to improve cancer detection, prevention, and therapeutic discovery for better precision-medicine treatments of cancer patients and those at risk for cancer.Cancer forms and progresses through a series of critical transitions-from pre-malignant to malignant states, from locally contained to metastatic disease, and from treatment-responsive to treatment-resistant tumors (Figure 1). Although specifics differ across tumor types and patients, all transitions involve complex dynamic interactions between diverse pre-malignant, malignant, and non-malignant cells (e.g., stroma cells and immune cells), often organized in specific patterns within the tumor

Many human diseases, arising from mutations of disease susceptibility genes (genetic diseases), are also associated with viral infections (virally implicated diseases), either in a directly causal manner or by indirect associations. Here we examine whether viral perturbations of host interactome may underlie such virally implicated disease relationships. Using as models two different human viruses, Epstein-Barr virus (EBV) and human papillomavirus (HPV), we find that host targets of viral proteins reside in network proximity to products of disease susceptibility genes. Expression changes in virally implicated disease tissues and comorbidity patterns cluster significantly in the network vicinity of viral targets. The topological proximity found between cellular targets of viral proteins and disease genes was exploited to uncover a novel pathway linking HPV to Fanconi anemia.

Human papillomavirus 18 (HPV18) is the second most carcinogenic HPV type, after HPV16, and it accounts for approximately 12% of squamous cell carcinoma (SCC) as well as 37% of adenocarcinoma (ADC) of the cervix worldwide. We aimed to evaluate the worldwide diversity and carcinogenicity of HPV18 genetic variants by sequencing the entire long control region ( In region-stratified analyses, there were no significant differences in the distribution of HPV18 variant lineages between cervical cancer cases and controls or between ADC and SCC. In conclusion, our findings do not support the role of HPV18 (sub)lineages for discriminating cancer risk or explaining why HPV18 is more strongly linked with ADC than SCC. IMPORTANCEThis is the largest and most geographically/ethnically diverse study of the genetic variation of HPV18 to date, providing a comprehensive reference for phylogenetic classification of HPV18 sublineages for epidemiological and biological studies.

Neutrophil extracellular trap (NET) formation can generate short-term, functional anucleate cytoplasts and trigger loss of cell viability. We demonstrated that the necroptotic cell death effector mixed lineage kinase domain-like (MLKL) translocated from the cytoplasm to the plasma membrane and stimulated downstream NADPH oxidase-independent ROS production, loss of cytoplasmic granules, breakdown of the nuclear membrane, chromatin decondensation, histone hypercitrullination, and extrusion of bacteriostatic NETs. This process was coordinated by receptor-interacting protein kinase-1 (RIPK1), which activated the caspase-8-dependent apoptotic or RIPK3/MLKL-dependent necroptotic death of mouse and human neutrophils. Genetic deficiency of RIPK3 and MLKL prevented NET formation but did not prevent cell death, which was because of residual caspase-8-dependent activity. Peptidylarginine deiminase 4 (PAD4) was activated downstream of RIPK1/RIPK3/MLKL and was required for maximal histone hypercitrullination and NET extrusion. This work defines a distinct signaling network that activates PAD4-dependent NET release for the control of methicillin-resistant (MRSA) infection.

Human papillomavirus (HPV) type 16 infection is an etiologic factor in a subset of head and neck squamous cell carcinomas (HNSCC). It is unknown if host genetic susceptibility modifies the HPV16-HNSCC association. DNA samples collected as part of a Boston area case-control study of HNSCC were genotyped for single-nucleotide polymorphisms (SNPs) from the National Cancer Institute's SNP500Cancer database. Analysis of demographic, phenotypic and genotypic data for 319 HNSCC cases and 495 frequency-matched controls was performed using unconditional logistic regression. All reported P-values are two sided. We identified a polymorphism in the sodium-dependent vitamin C transporter SLC23A2 that modifies the risk of HNSCC associated with HPV16 infection. Among those with a wild-type allele at SLC23A2, the risk of HNSCC associated with HPV16-positive serology was 5.0 (95% confidence interval (CI) = 3.2-7.8). However, among those with a homozygous variant genotype, the risk of HNSCC associated with HPV16 was attenuated [odds ratio (OR) = 2.8; 95% CI = 1.2-6.2]. Further, when we tested whether genotype modified the interaction between citrus exposure, HPV16, and HNSCC, we found a dramatically increased risk of HNSCC for those with a wild-type SLC23A2 allele, HPV16-positive serology and high citrus intake (OR = 7.4; 95% CI = 3.6-15.1). These results suggest that SLC23A2 genetic variation alters HPV16-associated HNSCC while also highlighting the important role of citrus exposure in this disease.

Cutaneous melanoma is a highly immunogenic disease, surgically curable at early stages, but life-threatening when metastatic. Here we integrate high-plex imaging, 3D high-resolution microscopy, and spatially-resolved micro-region transcriptomics to study immune evasion and immunoediting in primary melanoma. We find that recurrent cellular neighborhoods involving tumor, immune, and stromal cells change significantly along a progression axis from precursor states to melanoma in situ to invasive tumor. Hallmarks of immunosuppression are detectable by the precursor stage, and when tumors become locally invasive, a consolidated and spatially restricted suppressive environment forms along the tumor-stromal boundary. This environment is established by cytokine gradients that promote expression of MHC-II and IDO1 and by PDL1-expressing macrophages and dendritic cells engaging activated T cells. However, a few mm away, T cells synapse with melanoma cells in fields of tumor regression. Thus, invasion and immunoediting can co-exist within a few millimeters of each other in a single specimen.