Enrichment of Rare Variants of Hemophagocytic Lymphohistiocytosis Genes in Systemic Juvenile Idiopathic Arthritis.

Objective: To evaluate whether there is an enrichment of rare variants in familial hemophagocytic lymphohistiocytosis (HLH) genes and systemic juvenile idiopathic arthritis (sJIA) with or without macrophage activation syndrome (MAS). Methods: Targeted sequencing of HLH genes (LYST, PRF1, RAB27A, STX11, STXBP2, UNC13D) was performed in sJIA subjects from an established cohort. Sequence data from control subjects were obtained in silico (dbGaP:phs000280.v8.p2). Rare variant association testing (RVT) was performed with sequence kernel association test (SKAT) package. Significance was defined as p<0.05 after 100,000 permutations. Results: Sequencing data from 524 sJIA cases were jointly called and harmonized with exome-derived target data from 3000 controls. Quality control operations produced a set of 481 cases and 2924 ancestrally-matched control subjects. RVT of sJIA cases and controls revealed a significant association with rare protein-altering variants (minor allele frequency [MAF]<0.01) of STXBP2 (p=0.020), and ultra-rare variants (MAF<0.001) of STXBP2 (p=0.007) and UNC13D (p=0.045). A subanalysis of 32 cases with known MAS and 90 without revealed significant association of rare UNC13D variants (p=0.0047). Additionally, sJIA patients more often carried ≥2 HLH variants than did controls (p=0.007), driven largely by digenic combinations involving LYST. Conclusion: We identified an enrichment of rare HLH variants in sJIA patients compared with healthy controls, driven by STXBP2 and UNC13D. Biallelic variation in HLH genes was associated with sJIA, driven by LYST. Only UNC13D displayed enrichment in patients with MAS. This suggests that HLH variants may contribute to the pathophysiology of sJIA, even without MAS.

Ultra-rare genetic variation in relapsing polychondritis: a whole-exome sequencing study.

OBJECTIVE: Relapsing polychondritis (RP) is a systemic inflammatory disease of unknown aetiology. The objective of this study was to examine the contribution of rare genetic variations to RP. METHODS: We performed a case-control exome-wide rare variant association analysis that included 66 unrelated European American cases with RP and 2923 healthy controls (HC). Gene-level collapsing analysis was performed using Firth's logistics regression. Exploratory pathway analysis was performed using three different methods: Gene Set Enrichment Analysis, sequence kernel association test and higher criticism test. Plasma DCBLD2 levels were measured in patients with RP and HC using ELISA. RESULTS: In the collapsing analysis, RP was associated with a significantly higher burden of ultra-rare damaging variants in the DCBLD2 gene (7.6% vs 0.1%, unadjusted OR=79.8, p=2.93×10-7). Plasma DCBLD2 protein levels were significantly higher in RP than in HC (median 4.06 ng/µL vs 0.05 ng/µL, p<0.001). The pathway analysis revealed a statistically significant enrichment of genes in the tumour necrosis factor signalling pathway driven by rare damaging variants in RELB, RELA and REL using higher criticism test weighted by eigenvector centrality. CONCLUSIONS: This study identified specific rare variants in the DCBLD2 gene as a putative genetic risk factor for RP. These findings should be validated in additional patients with RP and supported by future functional experiments.

Implications of combined NOD2 and other gene mutations in autoinflammatory diseases.

NOD-like receptors (NLRs) are intracellular sensors associated with systemic autoinflammatory diseases (SAIDs). We investigated the largest monocentric cohort of patients with adult-onset SAIDs for coinheritance of low frequency and rare mutations in NOD2 and other autoinflammatory genes. Sixty-three patients underwent molecular testing for SAID gene panels after extensive clinical workups. Whole exome sequencing data from the large Atherosclerosis Risk in Communities (ARIC) study of individuals of European-American ancestry were used as control. Of 63 patients, 44 (69.8%) were found to carry combined gene variants in NOD2 and another gene (Group 1), and 19 (30.2%) were carriers only for NOD2 variants (Group 2). The genetic variant combinations in SAID patients were digenic in 66% (NOD2/MEFV, NOD2/NLRP12, NOD2/NLRP3, and NOD2/TNFRSF1A) and oligogenic in 34% of cases. These variant combinations were either absent or significantly less frequent in the control population. By phenotype-genotype correlation, approximately 40% of patients met diagnostic criteria for a specific SAID, and 60% had mixed diagnoses. There were no statistically significant differences in clinical manifestations between the two patient groups except for chest pain. Due to overlapping phenotypes and mixed genotypes, we have suggested a new term, "Mixed NLR-associated Autoinflammatory Disease ", to describe this disease scenario. Gene variant combinations are significant in patients with SAIDs primarily presenting with mixed clinical phenotypes. Our data support the proposition that immunological disease expression is modified by genetic background and environmental exposure. We provide a preliminary framework in diagnosis, management, and interpretation of the clinical scenario.

Transmission disequilibrium analysis of whole genome data in childhood-onset systemic lupus erythematosus.

Childhood-onset systemic lupus erythematosus (cSLE) patients are unique, with hallmarks of Mendelian disorders (early-onset and severe disease) and thus are an ideal population for genetic investigation of SLE. In this study, we use the transmission disequilibrium test (TDT), a family-based genetic association analysis that employs robust methodology, to analyze whole genome sequencing data. We aim to identify novel genetic associations in an ancestrally diverse, international cSLE cohort. Forty-two cSLE patients and 84 unaffected parents from 3 countries underwent whole genome sequencing. First, we performed TDT with single nucleotide variant (SNV)-based (common variants) using PLINK 1.9, and gene-based (rare variants) analyses using Efficient and Parallelizable Association Container Toolbox (EPACTS) and rare variant TDT (rvTDT), which applies multiple gene-based burden tests adapted for TDT, including the burden of rare variants test. Applying the GWAS standard threshold (5.0 × 10-8) to common variants, our SNV-based analysis did not return any genome-wide significant SNVs. The rare variant gene-based TDT analysis identified many novel genes significantly enriched in cSLE patients, including HNRNPUL2, a DNA repair protein, and DNAH11, a ciliary movement protein, among others. Our approach identifies several novel SLE susceptibility genes in an ancestrally diverse childhood-onset lupus cohort.

Ultra-Rare Genetic Variation in Relapsing Polychondritis: A Whole-Exome Sequencing Study.

Objective: Relapsing polychondritis (RP) is a systemic inflammatory disease of unknown etiology. The study objective was to examine the contribution of rare genetic variations in RP. Methods: We performed a case-control exome-wide rare variant association analysis including 66 unrelated European American RP cases and 2923 healthy controls. Gene-level collapsing analysis was performed using Firth's logistics regression. Pathway analysis was performed on an exploratory basis with three different methods: Gene Set Enrichment Analysis (GSEA), sequence kernel association test (SKAT) and higher criticism test. Plasma DCBLD2 levels were measured in patients with RP and healthy controls using enzyme-linked immunosorbent assay (ELISA). Results: In the collapsing analysis, RP was associated with higher burden of ultra-rare damaging variants in the DCBLD2 gene (7.6% vs 0.1%, unadjusted odds ratio = 79.8, p = 2.93 × 10-7). Patients with RP and ultra-rare damaging variants in DCBLD2 had a higher prevalence of cardiovascular manifestations. Plasma DCBLD2 protein levels were significantly higher in RP than healthy controls (5.9 vs 2.3, p < 0.001). Pathway analysis showed statistically significant enrichment of genes in the tumor necrosis factor (TNF) signaling pathway driven by rare damaging variants in RELB, RELA and REL using higher criticism test weighted by degree and eigenvector centrality. Conclusions: This study identified specific rare variants in DCBLD2 as putative genetic risk factors for RP. Genetic variation within the TNF pathway is also potentially associated with development of RP. These findings should be validated in additional patients with RP and supported by future functional experiments.

Risk of relapse of ANCA-associated vasculitis among patients homozygous for the proteinase 3 gene Val119Ile polymorphism.

BACKGROUND: The frequency of proteinase 3 gene (PRTN3) polymorphisms in patients with antineutrophil cytoplasmic antibody (ANCA)-associated vasculitis (AAV) is not fully characterised. We hypothesise that the presence of a PRTN3 gene polymorphism (single nucleotide polymorphism (SNP) rs351111) is relevant for clinical outcomes. METHODS: DNA variant calling for SNP rs351111 (chr.19:844020, c.355G>A) in PRTN3 gene assessed the allelic frequency in patients with PR3-AAV included in the Rituximab in ANCA-Associated Vasculitis trial. This was followed by RNA-seq variant calling to characterise the mRNA expression. We compared clinical outcomes between patients homozygous for PRTN3-Ile119 or PRTN3-Val119. RESULTS: Whole blood samples for DNA calling were available in 188 patients. 75 patients with PR3-AAV had the allelic variant: 62 heterozygous PRTN3-Val119Ile and 13 homozygous for PRTN3-Ile119. RNA-seq was available for 89 patients and mRNA corresponding to the allelic variant was found in 32 patients with PR3-AAV: 25 heterozygous PRTN3-Val119Ile and 7 homozygous for PRTN3-Ile119. The agreement between the DNA calling results and mRNA expression of the 86 patients analysed by both methods was 100%. We compared the clinical outcomes of 64 patients with PR3-AAV: 51 homozygous for PRTN3-Val119 and 13 homozygous for PRTN3-Ile119. The frequency of severe flares at 18 months in homozygous PRTN3-Ile119 was significantly higher when compared with homozygous PRTN3-Val119 (46.2% vs 19.6%, p=0.048). Multivariate analysis identified homozygous PR3-Ile119 as main predictor of severe relapse (HR 4.67, 95% CI 1.16 to 18.86, p=0.030). CONCLUSION: In patients with PR3-AAV, homozygosity for PRTN3-Val119Ile polymorphism appears associated with higher frequency of severe relapse. Further studies are necessary to better understand the association of this observation with the risk of severe relapse.

Constitutively active Lyn kinase causes a cutaneous small vessel vasculitis and liver fibrosis syndrome.

Neutrophilic inflammation is a hallmark of many monogenic autoinflammatory diseases; pathomechanisms that regulate extravasation of damaging immune cells into surrounding tissues are poorly understood. Here we identified three unrelated boys with perinatal-onset of neutrophilic cutaneous small vessel vasculitis and systemic inflammation. Two patients developed liver fibrosis in their first year of life. Next-generation sequencing identified two de novo truncating variants in the Src-family tyrosine kinase, LYN, p.Y508*, p.Q507* and a de novo missense variant, p.Y508F, that result in constitutive activation of Lyn kinase. Functional studies revealed increased expression of ICAM-1 on induced patient-derived endothelial cells (iECs) and of ß2-integrins on patient neutrophils that increase neutrophil adhesion and vascular transendothelial migration (TEM). Treatment with TNF inhibition improved systemic inflammation; and liver fibrosis resolved on treatment with the Src kinase inhibitor dasatinib. Our findings reveal a critical role for Lyn kinase in modulating inflammatory signals, regulating microvascular permeability and neutrophil recruitment, and in promoting hepatic fibrosis.

Genetic variations in NLRP3 and NLRP12 genes in adult-onset patients with autoinflammatory diseases: a comparative study.

Objectives: Cryopyrin-associated periodic syndrome or NLRP3-associated autoinflammatory disease (NLRP3-AID) and NLRP12-AID are both Mendelian disorders with autosomal dominant inheritance. Both diseases are rare, primarily reported in the pediatric population, and are thought to be phenotypically indistinguishable. We provide the largest cohort of adult-onset patients and compared these diseases and the gene variant frequency to population controls. Methods: A cohort of adult patients with AIDs were retrospectively studied. All underwent molecular testing for periodic fever syndrome gene panels after extensive and negative workups for systemic autoimmune and other related diseases. Patients were divided into Group 1- NLRP3-AID patients with NLRP3 variants (N=15), Group 2- NLRP12-AID with NLRP12 variants (N=14) and Group 3- both NLRP3 and NLRP12 (N=9) variants. Exome sequence data of two large control populations including the ARIC study were used to compare gene variant distribution and frequency. Results: All 38 patients were Caucasian with women accounting for 82%. Median age at diagnosis was 41 ± 23 years and the disease duration at diagnosis was 14 ± 13 years. We identified statistically significant differences between the groups, notably that gastrointestinal symptoms as well as evaluations for same were significantly more frequent in patients with NLRP12 variants, and headaches/dizziness were less common among the NLRP12 patients. Livedo reticularis was noted in four patients, exclusively among NLRP12 carriers. Over 50% of patients in Groups 1 and 2 carry low-frequency disease-associated variants, while the remaining carry rare variants. We unprecedently identified digenic variants, i.e., the coexistence of NLRP3 and NLRP12, which were either both low frequency or low frequency/rare. Allele frequencies of all variants identified in our cohort were either absent or significantly lower in the control populations, further strengthening the evidence of susceptibility of these variants to SAID phenotypes. Conclusion: Our comparative study shows that both NLRP3-AID and NLRP12-AID share similar clinical phenotypes, yet there are significant differences between them with regard to gastrointestinal and neurological symptoms. A spectrum of high to low genetic variations in both genes can contribute to SAID individually or in combination.

Genetically programmed alternative splicing of NEMO mediates an autoinflammatory disease phenotype.

Host defense and inflammation are regulated by the NF-κB essential modulator (NEMO), a scaffolding protein with a broad immune cell and tissue expression profile. Hypomorphic mutations in inhibitor of NF-κB kinase regulatory subunit gamma (IKBKG) encoding NEMO typically present with immunodeficiency. Here, we characterized a pediatric autoinflammatory syndrome in 3 unrelated male patients with distinct X-linked IKBKG germline mutations that led to overexpression of a NEMO protein isoform lacking the domain encoded by exon 5 (NEMO-Δex5). This isoform failed to associate with TANK binding kinase 1 (TBK1), and dermal fibroblasts from affected patients activated NF-κB in response to TNF but not TLR3 or RIG-I-like receptor (RLR) stimulation when isoform levels were high. By contrast, T cells, monocytes, and macrophages that expressed NEMO-Δex5 exhibited increased NF-κB activation and IFN production, and blood cells from these patients expressed a strong IFN and NF-κB transcriptional signature. Immune cells and TNF-stimulated dermal fibroblasts upregulated the inducible IKK protein (IKKi) that was stabilized by NEMO-Δex5, promoting type I IFN induction and antiviral responses. These data revealed how IKBKG mutations that lead to alternative splicing of skipping exon 5 cause a clinical phenotype we have named NEMO deleted exon 5 autoinflammatory syndrome (NDAS), distinct from the immune deficiency syndrome resulting from loss-of-function IKBKG mutations.

Somatic Mutations in UBA1 Define a Distinct Subset of Relapsing Polychondritis Patients With VEXAS.

OBJECTIVE: Somatic mutations in UBA1 cause a newly defined syndrome known as VEXAS (vacuoles, E1 enzyme, X-linked, autoinflammatory, somatic syndrome). More than 50% of patients currently identified as having VEXAS met diagnostic criteria for relapsing polychondritis (RP), but clinical features that characterize VEXAS within a cohort of patients with RP have not been defined. We undertook this study to define the prevalence of somatic mutations in UBA1 in patients with RP and to create an algorithm to identify patients with genetically confirmed VEXAS among those with RP. METHODS: Exome and targeted sequencing of UBA1 was performed in a prospective observational cohort of patients with RP. Clinical and immunologic characteristics of patients with RP were compared based on the presence or absence of UBA1 mutations. The random forest method was used to derive a clinical algorithm to identify patients with UBA1 mutations. RESULTS: Seven of 92 patients with RP (7.6%) had UBA1 mutations (referred to here as VEXAS-RP). Patients with VEXAS-RP were all male, were on average ≥45 years of age at disease onset, and commonly had fever, ear chondritis, skin involvement, deep vein thrombosis, and pulmonary infiltrates. No patient with VEXAS-RP had chondritis of the airways or costochondritis. Mortality was greater in VEXAS-RP than in RP (23% versus 4%; P = 0.029). Elevated acute-phase reactants and hematologic abnormalities (e.g., macrocytic anemia, thrombocytopenia, lymphopenia, multiple myeloma, myelodysplastic syndrome) were prevalent in VEXAS-RP. A decision tree algorithm based on male sex, a mean corpuscular volume >100 fl, and a platelet count <200 ×103 /µl differentiated VEXAS-RP from RP with 100% sensitivity and 96% specificity. CONCLUSION: Mutations in UBA1 were causal for disease in a subset of patients with RP. This subset of patients was defined by disease onset in the fifth decade of life or later, male sex, ear/nose chondritis, and hematologic abnormalities. Early identification is important in VEXAS given the associated high mortality rate.

Somatic Mutations in UBA1 and Severe Adult-Onset Autoinflammatory Disease.

BACKGROUND: Adult-onset inflammatory syndromes often manifest with overlapping clinical features. Variants in ubiquitin-related genes, previously implicated in autoinflammatory disease, may define new disorders. METHODS: We analyzed peripheral-blood exome sequence data independent of clinical phenotype and inheritance pattern to identify deleterious mutations in ubiquitin-related genes. Sanger sequencing, immunoblotting, immunohistochemical testing, flow cytometry, and transcriptome and cytokine profiling were performed. CRISPR-Cas9-edited zebrafish were used as an in vivo model to assess gene function. RESULTS: We identified 25 men with somatic mutations affecting methionine-41 (p.Met41) in UBA1, the major E1 enzyme that initiates ubiquitylation. (The gene UBA1 lies on the X chromosome.) In such patients, an often fatal, treatment-refractory inflammatory syndrome develops in late adulthood, with fevers, cytopenias, characteristic vacuoles in myeloid and erythroid precursor cells, dysplastic bone marrow, neutrophilic cutaneous and pulmonary inflammation, chondritis, and vasculitis. Most of these 25 patients met clinical criteria for an inflammatory syndrome (relapsing polychondritis, Sweet's syndrome, polyarteritis nodosa, or giant-cell arteritis) or a hematologic condition (myelodysplastic syndrome or multiple myeloma) or both. Mutations were found in more than half the hematopoietic stem cells, including peripheral-blood myeloid cells but not lymphocytes or fibroblasts. Mutations affecting p.Met41 resulted in loss of the canonical cytoplasmic isoform of UBA1 and in expression of a novel, catalytically impaired isoform initiated at p.Met67. Mutant peripheral-blood cells showed decreased ubiquitylation and activated innate immune pathways. Knockout of the cytoplasmic UBA1 isoform homologue in zebrafish caused systemic inflammation. CONCLUSIONS: Using a genotype-driven approach, we identified a disorder that connects seemingly unrelated adult-onset inflammatory syndromes. We named this disorder the VEXAS (vacuoles, E1 enzyme, X-linked, autoinflammatory, somatic) syndrome. (Funded by the NIH Intramural Research Programs and the EU Horizon 2020 Research and Innovation Program.).

Somatic SMAD3-activating mutations cause melorheostosis by up-regulating the TGF-ß/SMAD pathway.

Melorheostosis is a rare sclerosing dysostosis characterized by asymmetric exuberant bone formation. Recently, we reported that somatic mosaicism for MAP2K1-activating mutations causes radiographical "dripping candle wax" melorheostosis. We now report somatic SMAD3 mutations in bone lesions of four unrelated patients with endosteal pattern melorheostosis. In vitro, the SMAD3 mutations stimulated the TGF-ß pathway in osteoblasts, enhanced nuclear translocation and target gene expression, and inhibited proliferation. Osteoblast differentiation and mineralization were stimulated by the SMAD3 mutation, consistent with higher mineralization in affected than in unaffected bone, but differing from MAP2K1 mutation-positive melorheostosis. Conversely, osteoblast differentiation and mineralization were inhibited when osteogenesis of affected osteoblasts was driven in the presence of BMP2. Transcriptome profiling displayed that TGF-ß pathway activation and ossification-related processes were significantly influenced by the SMAD3 mutation. Co-expression clustering illuminated melorheostosis pathophysiology, including alterations in ECM organization, cell growth, and interferon signaling. These data reveal antagonism of TGF-ß/SMAD3 activation by BMP signaling in SMAD3 mutation-positive endosteal melorheostosis, which may guide future therapies.

Distinct interferon signatures and cytokine patterns define additional systemic autoinflammatory diseases.

BACKGROUNDUndifferentiated systemic autoinflammatory diseases (USAIDs) present diagnostic and therapeutic challenges. Chronic interferon (IFN) signaling and cytokine dysregulation may identify diseases with available targeted treatments.METHODSSixty-six consecutively referred USAID patients underwent underwent screening for the presence of an interferon signature using a standardized type-I IFN-response-gene score (IRG-S), cytokine profiling, and genetic evaluation by next-generation sequencing.RESULTSThirty-six USAID patients (55%) had elevated IRG-S. Neutrophilic panniculitis (40% vs. 0%), basal ganglia calcifications (46% vs. 0%), interstitial lung disease (47% vs. 5%), and myositis (60% vs. 10%) were more prevalent in patients with elevated IRG-S. Moderate IRG-S elevation and highly elevated serum IL-18 distinguished 8 patients with pulmonary alveolar proteinosis (PAP) and recurrent macrophage activation syndrome (MAS). Among patients with panniculitis and progressive cytopenias, 2 patients were compound heterozygous for potentially novel LRBA mutations, 4 patients harbored potentially novel splice variants in IKBKG (which encodes NF-κB essential modulator [NEMO]), and 6 patients had de novo frameshift mutations in SAMD9L. Of additional 12 patients with elevated IRG-S and CANDLE-, SAVI- or Aicardi-Goutières syndrome-like (AGS-like) phenotypes, 5 patients carried mutations in either SAMHD1, TREX1, PSMB8, or PSMG2. Two patients had anti-MDA5 autoantibody-positive juvenile dermatomyositis, and 7 could not be classified. Patients with LRBA, IKBKG, and SAMD9L mutations showed a pattern of IRG elevation that suggests prominent NF-κB activation different from the canonical interferonopathies CANDLE, SAVI, and AGS.CONCLUSIONSIn patients with elevated IRG-S, we identified characteristic clinical features and 3 additional autoinflammatory diseases: IL-18-mediated PAP and recurrent MAS (IL-18PAP-MAS), NEMO deleted exon 5-autoinflammatory syndrome (NEMO-NDAS), and SAMD9L-associated autoinflammatory disease (SAMD9L-SAAD). The IRG-S expands the diagnostic armamentarium in evaluating USAIDs and points to different pathways regulating IRG expression.TRIAL REGISTRATIONClinicalTrials.gov NCT02974595.FUNDINGThe Intramural Research Program of the NIH, NIAID, NIAMS, and the Clinical Center.

Inactivating Mutation in IRF8 Promotes Osteoclast Transcriptional Programs and Increases Susceptibility to Tooth Root Resorption.

This is the first study to our knowledge to report a novel mutation in the interferon regulatory factor 8 gene (IRF8G388S ) associated with multiple idiopathic tooth root resorption, a form of periodontal disease. The IRF8G388S variant in the highly conserved C-terminal motif is predicted to alter the protein structure, likely impairing IRF8 function. Functional assays demonstrated that the IRF8G388S mutant promoted osteoclastogenesis and failed to inhibit NFATc1-dependent transcriptional activation when compared with IRF8WT control. Further, similar to subjects with heterozygous IRF8G388S mutation, Irf8+/- mice exhibited increased osteoclast activity in the mandibular alveolar bone surrounding molar teeth. Immunohistochemistry illustrated increased NFATc1 expression in the dentoalveolar region of Irf8-/- and Irf8+/- mice when compared with Irf8+/+ controls. Genomewide analyses revealed that IRF8 constitutively bound to regulatory regions of several thousand genes in osteoclast precursors, and genetic aberration of IRF8 significantly enhanced many osteoclast-specific transcripts. Collectively, this study delineates the critical role of IRF8 in defining osteoclast lineage and osteoclast transcriptional program, which may help in better understanding of various osteoclast-mediated disorders, including periodontal disease. © 2019 American Society for Bone and Mineral Research.

Melorheostotic Bone Lesions Caused by Somatic Mutations in MAP2K1 Have Deteriorated Microarchitecture and Periosteal Reaction.

Melorheostosis is a rare non-hereditary condition characterized by dense hyperostotic lesions with radiographic "dripping candle wax" appearance. Somatic activating mutations in MAP2K1 have recently been identified as a cause of melorheostosis. However, little is known about the development, composition, structure, and mechanical properties of the bone lesions. We performed a multi-method phenotype characterization of material properties in affected and unaffected bone biopsy samples from six melorheostosis patients with MAP2K1 mutations. On standard histology, lesions show a zone with intensively remodeled osteonal-like structure and prominent osteoid accumulation, covered by a shell formed through bone apposition, consisting of compact multi-layered lamellae oriented parallel to the periosteal surface and devoid of osteoid. Compared with unaffected bone, melorheostotic bone has lower average mineralization density measured by quantitative backscattered electron imaging (CaMean: -4.5%, p = 0.04). The lamellar portion of the lesion is even less mineralized, possibly because the newly deposited material has younger tissue age. Affected bone has higher porosity by micro-CT, due to increased tissue vascularity and elevated 2D-microporosity (osteocyte lacunar porosity: +39%, p = 0.01) determined on quantitative backscattered electron images. Furthermore, nano-indentation modulus characterizing material hardness and stiffness was strictly dependent on tissue mineralization (correlation with typical calcium concentration, CaPeak: r = 0.8984, p = 0.0150, and r = 0.9788, p = 0.0007, respectively) in both affected and unaffected bone, indicating that the surgical hardness of melorheostotic lesions results from their lamellar structure. The results suggest a model for pathophysiology of melorheostosis caused by somatic activating mutations in MAP2K1, in which the genetically induced gradual deterioration of bone microarchitecture triggers a periosteal reaction, similar to the process found to occur after bone infection or local trauma, and leads to an overall cortical outgrowth. The micromechanical properties of the lesions reflect their structural heterogeneity and correlate with local variations in mineral content, tissue age, and remodeling rates, in the same way as normal bone. © 2018 American Society for Bone and Mineral Research.

Distinct Clinical and Pathological Features of Melorheostosis Associated With Somatic MAP2K1 Mutations.

Melorheostosis is a rare hyperostotic disease of the long bones classically characterized by a "dripping candle-wax" radiographic appearance. We recently described somatic activating mutations in MAP2K1 as a cause of melorheostosis. Here, we report distinguishing characteristics of patients with MAP2K1-positive melorheostosis. Fifteen unrelated patients with radiographic appearance of melorheostosis underwent paired biopsies of affected and unaffected bone for whole-exome sequencing, histology, and cell culture. Eight patients with mutations in MAP2K1 in affected bone were compared to the seven MAP2K1-negative patients to identify distinguishing characteristics. Patients with MAP2K1-positive melorheostosis had a distinct phenotype with classic "dripping candle-wax" appearance on radiographs (p = 0.01), characteristic vascular lesions on skin overlying affected bone (p = 0.01), and higher prevalence of extraosseous mineralization and joint involvement (p = 0.04 for both). Melorheostotic bone from both MAP2K1-positive and MAP2K1-negative patients showed two zones of distinct morphology-an outer segment of parallel layers of primary lamellar bone and a deeper zone of intensely remodeled highly porous osteonal-like bone. Affected bone from MAP2K1-positive patients showed excessive osteoid (p = 0.0012), increased number of osteoblasts (p = 0.012) and osteoclasts (p = 0.04), and increased vascularity on histology in comparison to paired unaffected bone which was not seen in affected bone in most MAP2K1-negative patients. The identification of a distinct phenotype of patients with MAP2K1-positive melorheostosis demonstrates clinical and genetic heterogeneity among patients with the disease. Further studies are needed to better understand the underlying pathophysiology and associated skin findings. © 2018 American Society for Bone and Mineral Research.

Somatic activating mutations in MAP2K1 cause melorheostosis.

Melorheostosis is a sporadic disease of uncertain etiology characterized by asymmetric bone overgrowth and functional impairment. Using whole exome sequencing, we identify somatic mosaic MAP2K1 mutations in affected, but not unaffected, bone of eight unrelated patients with melorheostosis. The activating mutations (Q56P, K57E and K57N) cluster tightly in the MEK1 negative regulatory domain. Affected bone displays a mosaic pattern of increased p-ERK1/2 in osteoblast immunohistochemistry. Osteoblasts cultured from affected bone comprise two populations with distinct p-ERK1/2 levels by flow cytometry, enhanced ERK1/2 activation, and increased cell proliferation. However, these MAP2K1 mutations inhibit BMP2-mediated osteoblast mineralization and differentiation in vitro, underlying the markedly increased osteoid detected in affected bone histology. Mosaicism is also detected in the skin overlying bone lesions in four of five patients tested. Our data show that the MAP2K1 oncogene is important in human bone formation and implicate MEK1 inhibition as a potential treatment avenue for melorheostosis.