We have sequenced 463 presenting cases of myeloma entered into the UK Myeloma XI study using whole exome sequencing. Here we identify mutations induced as a consequence of misdirected AID in the partner oncogenes of IGH translocations, which are activating and associated with impaired clinical outcome. An APOBEC mutational signature is seen in 3.8% of cases and is linked to the translocation-mediated deregulation of MAF and MAFB, a known poor prognostic factor. Patients with this signature have an increased mutational load and a poor prognosis. Loss of MAF or MAFB expression results in decreased APOBEC3B and APOBEC4 expression, indicating a transcriptional control mechanism. Kataegis, a further mutational pattern associated with APOBEC deregulation, is seen at the sites of the MYC translocation. The APOBEC mutational signature seen in myeloma is, therefore, associated with poor prognosis primary and secondary translocations and the molecular mechanisms involved in generating them.

To identify risk variants for multiple myeloma, we conducted a genome-wide association study of 1,675 individuals with multiple myeloma and 5,903 control subjects. We identified risk loci for multiple myeloma at 3p22.1 (rs1052501 in ULK4; odds ratio (OR) = 1.32; P = 7.47 × 10(-9)) and 7p15.3 (rs4487645, OR = 1.38; P = 3.33 × 10(-15)). In addition, we observed a promising association at 2p23.3 (rs6746082, OR = 1.29; P = 1.22 × 10(-7)). Our study identifies new genomic regions associated with multiple myeloma risk that may lead to new etiological insights.

Genomic regions of homozygosity (ROH), detectable in outbred populations, have been implicated as determinants of inherited risk. To examine whether ROH is associated with risk of multiple myeloma (MM), we performed whole-genome homozygosity analysis using single-nucleotide polymorphism genotype data from 2,282 MM cases and 5,197 controls, with replication in an additional 878 MM cases and 7,083 controls. Globally, the distribution of ROH between cases and controls was not significantly different. However, one ROH at chromosome 9q21, harboring the B-cell transcription factor gene KLF9, showed evidence of a consistent association and may therefore warrant further investigation as a candidate risk factor for MM. Overall, our analysis provides little support for a homozygosity signature being a significant factor in MM risk.

VCP/p97 regulates numerous cellular functions by mediating protein degradation through its segregase activity. Its key role in governing protein homoeostasis has made VCP/p97 an appealing anticancer drug target. Here, we provide evidence that VCP/p97 acts as a regulator of cellular metabolism. We found that VCP/p97 was tied to multiple metabolic processes on the gene expression level in a diverse range of cancer cell lines and in patient-derived multiple myeloma cells. Cellular VCP/p97 dependency to maintain proteostasis was increased under conditions of glucose and glutamine limitation in a range of cancer cell lines from different tissues. Moreover, glutamine depletion led to increased VCP/p97 expression, whereas VCP/p97 inhibition perturbed metabolic processes and intracellular amino acid turnover. GCN2, an amino acid-sensing kinase, attenuated stress signalling and cell death triggered by VCP/p97 inhibition and nutrient shortages and modulated ERK activation, autophagy, and glycolytic metabolite turnover. Together, our data point to an interconnected role of VCP/p97 and GCN2 in maintaining cancer cell metabolic and protein homoeostasis.

The clustering of different types of B-cell malignancies in families raises the possibility of shared aetiology. To examine this, we performed cross-trait linkage disequilibrium (LD)-score regression of multiple myeloma (MM) and chronic lymphocytic leukaemia (CLL) genome-wide association study (GWAS) data sets, totalling 11,734 cases and 29,468 controls. A significant genetic correlation between these two B-cell malignancies was shown (Rg = 0.4, P = 0.0046). Furthermore, four of the 45 known CLL risk loci were shown to associate with MM risk and five of the 23 known MM risk loci associate with CLL risk. By integrating eQTL, Hi-C and ChIP-seq data, we show that these pleiotropic risk loci are enriched for B-cell regulatory elements and implicate B-cell developmental genes. These data identify shared biological pathways influencing the development of CLL and, MM and further our understanding of the aetiological basis of these B-cell malignancies.

B-cell malignancies (BCM) originate from the same cell of origin, but at different maturation stages and have distinct clinical phenotypes. Although genetic risk variants for individual BCMs have been identified, an agnostic, genome-wide search for shared genetic susceptibility has not been performed. We explored genome-wide association studies of chronic lymphocytic leukaemia (CLL, N = 1,842), Hodgkin lymphoma (HL, N = 1,465) and multiple myeloma (MM, N = 3,790). We identified a novel pleiotropic risk locus at 3q22.2 (NCK1, rs11715604, P = 1.60 × 10-9) with opposing effects between CLL (P = 1.97 × 10-8) and HL (P = 3.31 × 10-3). Eight established non-HLA risk loci showed pleiotropic associations. Within the HLA region, Ser37 + Phe37 in HLA-DRB1 (P = 1.84 × 10-12) was associated with increased CLL and HL risk (P = 4.68 × 10-12), and reduced MM risk (P = 1.12 × 10-2), and Gly70 in HLA-DQB1 (P = 3.15 × 10-10) showed opposing effects between CLL (P = 3.52 × 10-3) and HL (P = 3.41 × 10-9). By integrating eQTL, Hi-C and ChIP-seq data, we show that the pleiotropic risk loci are enriched for B-cell regulatory elements, as well as an over-representation of binding of key B-cell transcription factors. These data identify shared biological pathways influencing the development of CLL, HL and MM. The identification of these risk loci furthers our understanding of the aetiological basis of BCMs.

Genome-wide association studies (GWAS) have transformed our understanding of susceptibility to multiple myeloma (MM), but much of the heritability remains unexplained. We report a new GWAS, a meta-analysis with previous GWAS and a replication series, totalling 9974 MM cases and 247,556 controls of European ancestry. Collectively, these data provide evidence for six new MM risk loci, bringing the total number to 23. Integration of information from gene expression, epigenetic profiling and in situ Hi-C data for the 23 risk loci implicate disruption of developmental transcriptional regulators as a basis of MM susceptibility, compatible with altered B-cell differentiation as a key mechanism. Dysregulation of autophagy/apoptosis and cell cycle signalling feature as recurrently perturbed pathways. Our findings provide further insight into the biological basis of MM.

Myeloma is characterized by a highly variable clinical outcome. Despite the effectiveness of high-dose therapy, 15% of patients relapse within 1 year. We show that these cases also have a significantly shorter post-relapse survival compared to the others (median 14.9 months vs. 40 months, p = 8.03 × 10(- 14)). There are no effective approaches to define this potentially distinct biological group such that treatment could be altered. In this work a series of uniformly treated patients with myeloma were used to develop a gene expression profiling (GEP)-based signature to identify this high risk clinical behavior. Gene enrichment analyses applied to the top differentially expressed genes showed a significant enrichment of epigenetic regulators as well as "stem cell" myeloma genes. A derived 17-gene signature effectively identifies patients at high risk of early relapse as well as impaired overall survival. Integrative genomic analyses showed that epigenetic mechanisms may play an important role on transcription of these genes.

Human cytomegalovirus (HCMV) is a ubiquitous virus that is a major pathogen in newborns and immunocompromised or immunosuppressed patients. HCMV infects a wide variety of cell types using distinct entry pathways that involve different forms of the gH/gL glycoprotein: gH/gL/gO and gH/gL/UL128-131 as well as the viral fusion glycoprotein, gB. However, the minimal or core fusion machinery (sufficient for cell-cell fusion) is just gH/gL and gB. Here, we demonstrate that HCMV gB and gH/gL form a stable complex early after their synthesis and in the absence of other viral proteins. gH/gL can interact with gB mutants that are unable to mediate cell-cell fusion. gB-gH/gL complexes included as much as 16-50% of the total gH/gL in HCMV virus particles. In contrast, only small amounts of gH/gL/gO and gH/gL/UL128-131 complexes were found associated with gB. All herpesviruses express gB and gH/gL molecules and most models describing herpesvirus entry suggest that gH/gL interacts with gB to mediate membrane fusion, although there is no direct evidence for this. For herpes simplex virus (HSV-1) it has been suggested that after receptor binding gH/gL binds to gB either just before, or coincident with membrane fusion. Therefore, our results have major implications for these models, demonstrating that HCMV gB and gH/gL forms stable gB-gH/gL complexes that are incorporated virions without receptor binding or membrane fusion. Moreover, our data is the best support to date for the proposal that gH/gL interacts with gB.

A number of specific chromosomal abnormalities define the subgroups of multiple myeloma. In a meta-analysis of two genome-wide association studies of multiple myeloma including a total of 1,661 affected individuals, we investigated risk for developing a specific tumor karyotype. The t(11;14)(q13;q32) translocation in which CCND1 is placed under the control of the immunoglobulin heavy chain enhancer was strongly associated with the CCND1 c.870G>A polymorphism (P = 7.96 × 10(-11)). These results provide a model in which a constitutive genetic factor is associated with risk of a specific chromosomal translocation.

The mechanisms involved in the progression from monoclonal gammopathy of undetermined significance (MGUS) and smoldering myeloma (SMM) to malignant multiple myeloma (MM) and plasma cell leukemia (PCL) are poorly understood but believed to involve the sequential acquisition of genetic hits. We performed exome and whole-genome sequencing on a series of MGUS (n=4), high-risk (HR)SMM (n=4), MM (n=26) and PCL (n=2) samples, including four cases who transformed from HR-SMM to MM, to determine the genetic factors that drive progression of disease. The pattern and number of non-synonymous mutations show that the MGUS disease stage is less genetically complex than MM, and HR-SMM is similar to presenting MM. Intraclonal heterogeneity is present at all stages and using cases of HR-SMM, which transformed to MM, we show that intraclonal heterogeneity is a typical feature of the disease. At the HR-SMM stage of disease, the majority of the genetic changes necessary to give rise to MM are already present. These data suggest that clonal progression is the key feature of transformation of HR-SMM to MM and as such the invasive clinically predominant clone typical of MM is already present at the SMM stage and would be amenable to therapeutic intervention at that stage.

Survival following a diagnosis of multiple myeloma (MM) varies between patients and some of these differences may be a consequence of inherited genetic variation. In this study, to identify genetic markers associated with MM overall survival (MM-OS), we conduct a meta-analysis of four patient series of European ancestry, totalling 3,256 patients with 1,200 MM-associated deaths. Each series is genotyped for ∼600,000 single nucleotide polymorphisms across the genome; genotypes for six million common variants are imputed using 1000 Genomes Project and UK10K as the reference. The association between genotype and OS is assessed by Cox proportional hazards model adjusting for age, sex, International staging system and treatment. We identify a locus at 6q25.1 marked by rs12374648 associated with MM-OS (hazard ratio=1.34, 95% confidence interval=1.22-1.48, P=4.69 × 10(-9)). Our findings have potential clinical implications since they demonstrate that inherited genotypes can provide prognostic information in addition to conventional tumor acquired prognostic factors.

Human cytomegalovirus (HCMV) entry involves trimer (gH/gL/gO) that interacts with PDGFRα in fibroblasts. Entry into epithelial and endothelial cells requires trimer, which binds unidentified receptors, and pentamer (gH/gL/UL128-131), which binds neuropilin-2. To identify functionally important domains in trimer, we screened an overlapping 20-mer gO peptide library and identified two sets of peptides: 19/20 (a.a. 235-267) and 32/33 (a.a. 404-436) that could block virus entry. Soluble trimer containing wild type gO blocked HCMV entry, whereas soluble trimers with the 19/20 or 32/33 sequences mutated did not block entry. Interestingly, the mutant trimers retained the capacity to bind to cellular receptors including PDGFRα. Peptide 19/20 and 32/33 sequences formed a lobe extending from the surface of gO and an adjacent concave structure, respectively. Neither of these sets of sequences contacted PDGFRα. Instead, our data support a model in which the 19/20 and 32/33 trimer sequences function downstream of receptor binding, e.g. trafficking of HCMV into endosomes or binding to gB for entry fusion. We also screened for peptides that bound antibodies (Abs) in human sera, observing that peptides 20 and 26 bound Abs. These peptides engendered neutralizing Abs (NAbs) after immunization of rabbits and could pull out NAbs from human sera. Peptides 20 and 26 sequences represent the first NAb epitopes identified in trimer. These studies describe two important surfaces on gO defined by: i) peptides 19/20 and 32/33, which apparently act downstream of receptor binding and ii) peptide 26 that interacts with PDGFRα. Both these surfaces are targets of NAbs.

Cross-species studies enable rapid translational discovery and produce the broadest impact when both mechanism and phenotype are consistent across organisms. We developed a knock-in mouse that biologically recapitulates a common human mutation in the gene for fatty acid amide hydrolase (FAAH) (C385A; rs324420), the primary catabolic enzyme for the endocannabinoid anandamide. This common polymorphism impacts the expression and activity of FAAH, thereby increasing anandamide levels. Here, we show that the genetic knock-in mouse and human variant allele carriers exhibit parallel alterations in biochemisty, neurocircuitry and behaviour. Specifically, there is reduced FAAH expression associated with the variant allele that selectively enhances fronto-amygdala connectivity and fear extinction learning, and decreases anxiety-like behaviours. These results suggest a gain of function in fear regulation and may indicate for whom and for what anxiety symptoms FAAH inhibitors or exposure-based therapies will be most efficacious, bridging an important translational gap between the mouse and human.

Anxiety disorders peak in incidence during adolescence, a developmental window that is marked by dynamic changes in gene expression, endocannabinoid signaling, and frontolimbic circuitry. We tested whether genetic alterations in endocannabinoid signaling related to a common polymorphism in fatty acid amide hydrolase (FAAH), which alters endocannabinoid anandamide (AEA) levels, would impact the development of frontolimbic circuitry implicated in anxiety disorders. In a pediatric imaging sample of over 1,000 3- to 21-y-olds, we show effects of the FAAH genotype specific to frontolimbic connectivity that emerge by ∼12 y of age and are paralleled by changes in anxiety-related behavior. Using a knock-in mouse model of the FAAH polymorphism that controls for genetic and environmental backgrounds, we confirm phenotypic differences in frontoamygdala circuitry and anxiety-related behavior by postnatal day 45 (P45), when AEA levels begin to decrease, and also, at P75 but not before. These results, which converge across species and level of analysis, highlight the importance of underlying developmental neurobiology in the emergence of genetic effects on brain circuitry and function. Moreover, the results have important implications for the identification of risk for disease and precise targeting of treatments to the biological state of the developing brain as a function of developmental changes in gene expression and neural circuit maturation.

The activating receptor, NKG2D, is expressed on a variety of immune effector cells and recognizes divergent families of major histocompatibility complex (MHC) class I-related ligands, including the MIC and ULBP proteins. Infection, stress, or transformation can induce NKG2D ligand expression, resulting in effector cell activation and killing of the ligand-expressing target cell. The human cytomegalovirus (HCMV) membrane glycoprotein, UL16, binds to three of the five known ligands for human NKG2D. UL16 is retained in the endoplasmic reticulum and cis-Golgi apparatus of cells and causes MICB to be similarly retained and stabilized within cells. Coexpression of UL16 markedly reduces cell surface levels of MICB, ULBP1, and ULBP2, and decreases susceptibility to natural killer cell-mediated cytotoxicity. Domain swapping experiments demonstrate that the transmembrane and cytoplasmic domains of UL16 are important for intracellular retention of UL16, whereas the ectodomain of UL16 participates in down-regulation of NKG2D ligands. The intracellular sequestration of NKG2D ligands by UL16 represents a novel HCMV immune evasion mechanism to add to the well-documented viral strategies directed against antigen presentation by classical MHC molecules.

Multiple myeloma (MM) is a malignancy of plasma cells. Genome-wide association studies have shown that variation at 5q15 influences MM risk. Here, we have sought to decipher the causal variant at 5q15 and the mechanism by which it influences tumorigenesis. We show that rs6877329 G > C resides in a predicted enhancer element that physically interacts with the transcription start site of ELL2. The rs6877329-C risk allele is associated with reduced enhancer activity and lowered ELL2 expression. Since ELL2 is critical to the B cell differentiation process, reduced ELL2 expression is consistent with inherited genetic variation contributing to arrest of plasma cell development, facilitating MM clonal expansion. These data provide evidence for a biological mechanism underlying a hereditary risk of MM at 5q15.

Obsessive compulsive disorder (OCD) is substantially heritable, but few molecular genetic risk factors have been identified. Knockout mice lacking SLIT and NTRK-Like Family, Member 5 (SLITRK5) display OCD-like phenotypes including serotonin reuptake inhibitor-sensitive pathologic grooming, and corticostriatal dysfunction. Thus, mutations that impair SLITRK5 function may contribute to the genetic risk for OCD. We re-sequenced the protein-coding sequence of the human SLITRK5 gene (SLITRK5) in three hundred and seventy seven OCD subjects and compared rare non-synonymous mutations (RNMs) in that sample with similar mutations in the 1000 Genomes database. We also performed in silico assessments and in vitro functional synaptogenesis assays on the Slitrk5 mutations identified. We identified four RNM's among these OCD subjects. There were no significant differences in the prevalence or in silico effects of rare non-synonymous mutations in the OCD sample versus controls. Direct functional testing of recombinant SLITRK5 proteins found that all mutations identified in OCD subjects impaired synaptogenic activity whereas none of the pseudo-matched mutations identified in 1000 Genomes controls had significant effects on SLITRK5 function (Fisher's exact test P = 0.028). These results demonstrate that rare functional mutations in SLITRK5 contribute to the genetic risk for OCD in human populations. They also highlight the importance of biological characterization of allelic effects in understanding genotype-phenotype relationships as there were no statistical differences in overall prevalence or bioinformatically predicted effects of OCD case versus control mutations. Finally, these results converge with others to highlight the role of aberrant synaptic function in corticostriatal neurons in the pathophysiology of OCD.

Human cytomegalovirus (HCMV) replicates in many diverse cell types in vivo, and entry into different cells involves distinct entry mechanisms and different envelope glycoproteins. HCMV glycoprotein gB is thought to act as the virus fusogen, apparently after being triggered by different gH/gL proteins that bind distinct cellular receptors or entry mediators. A trimer of gH/gL/gO is required for entry into all cell types, and entry into fibroblasts involves trimer binding to platelet-derived growth factor receptor alpha (PDGFRα). HCMV entry into biologically relevant epithelial and endothelial cells and monocyte-macrophages also requires a pentamer, gH/gL complexed with UL128, UL130, and UL131, and there is evidence that the pentamer binds unidentified receptors. We screened an epithelial cell cDNA library and identified the cell surface protein CD147, which increased entry of pentamer-expressing HCMV into HeLa cells but not entry of HCMV that lacked the pentamer. A panel of CD147-specific monoclonal antibodies inhibited HCMV entry into epithelial and endothelial cells, but not entry into fibroblasts. shRNA silencing of CD147 in endothelial cells inhibited HCMV entry but not entry into fibroblasts. CD147 colocalized with HCMV particles on cell surfaces and in endosomes. CD147 also promoted cell-cell fusion induced by expression of pentamer and gB in epithelial cells. However, soluble CD147 did not block HCMV entry and trimer and pentamer did not bind directly to CD147, supporting the hypothesis that CD147 acts indirectly through other proteins. CD147 represents the first HCMV entry mediator that specifically functions to promote entry of pentamer-expressing HCMV into epithelial and endothelial cells.IMPORTANCE Human cytomegalovirus infects nearly 80% of the world's population and causes significant morbidity and mortality. The current method of treatment involves the use of antiviral agents that are prone to resistance and can be highly toxic to patients; currently, there is no vaccine against HCMV available. HCMV infections involve virus dissemination throughout the body, infecting a wide variety of tissues; however, the mechanism of spread is not well understood, particularly with regard to which cellular proteins are utilized by HCMV to establish infection. This report describes the characterization of a newly identified cellular molecule that affects HCMV entry into epithelial and endothelial cells. These results will lead to a better understanding of HCMV pathogenesis and have implications for the development of future therapeutics.

Human cytomegalovirus (HCMV) establishes lifelong infection with recurrent episodes of virus production and shedding despite the presence of adaptive immunological memory responses including HCMV immune immunoglobulin G (IgG). Very little is known how HCMV evades from humoral and cellular IgG-dependent immune responses, the latter being executed by cells expressing surface receptors for the Fc domain of IgG (FcγRs). Remarkably, HCMV expresses the RL11-encoded gp34 and UL119-118-encoded gp68 type I transmembrane glycoproteins which bind Fcγ with nanomolar affinity. Using a newly developed FcγR activation assay, we tested if the HCMV-encoded Fcγ binding proteins (HCMV FcγRs) interfere with individual host FcγRs. In absence of gp34 or/and gp68, HCMV elicited a much stronger activation of FcγRIIIA/CD16, FcγRIIA/CD32A and FcγRI/CD64 by polyclonal HCMV-immune IgG as compared to wildtype HCMV. gp34 and gp68 co-expression culminates in the late phase of HCMV replication coinciding with the emergence of surface HCMV antigens triggering FcγRIII/CD16 responses by polyclonal HCMV-immune IgG. The gp34- and gp68-dependent inhibition of HCMV immune IgG was fully reproduced when testing the activation of primary human NK cells. Their broad antagonistic function towards FcγRIIIA, FcγRIIA and FcγRI activation was also recapitulated in a gain-of-function approach based on humanized monoclonal antibodies (trastuzumab, rituximab) and isotypes of different IgG subclasses. Surface immune-precipitation showed that both HCMV-encoded Fcγ binding proteins have the capacity to bind trastuzumab antibody-HER2 antigen complexes demonstrating simultaneous linkage of immune IgG with antigen and the HCMV inhibitors on the plasma membrane. Our studies reveal a novel strategy by which viral FcγRs can compete for immune complexes against various Fc receptors on immune cells, dampening their activation and antiviral immunity.