Severe manifestation of Rauch-Azzarello syndrome associated with biallelic deletion of CTNND2

CTNND2 encodes δ-catenin, a component of an adherens junction complex, and plays an important role in neuronal structure and function. To date, only heterozygous loss-of-function CTNND2 variants have been associated with mild neurodevelopmental delay and behavioral anomalies, a condition, which we named Rauch-Azzarello syndrome. Here, we report three siblings of a consanguineous family of Syrian descent with a homozygous deletion encompassing the last 19 exons of CTNND2 predicted to disrupt the transcript. All presented with severe neurodevelopmental delay with absent speech, profound motor delay, stereotypic behavior, microcephaly, short stature, muscular hypotonia with lower limb hypertonia, and variable eye anomalies. The parents and the fourth sibling were heterozygous carriers of the deletion and exhibited mild neurodevelopmental impairment resembling that of the previously described heterozygous individuals. The present study unveils a severe manifestation of CTNND2-associated Rauch-Azzarello syndrome attributed to biallelic loss-of-function aberrations, clinically distinct from the already described mild presentation of heterozygous individuals. Furthermore, we demonstrate novel clinical features in homozygous individuals that have not been reported in heterozygous cases to date.

P2RX7 gene variants associate with altered inflammasome assembly and reduced pyroptosis in chronic nonbacterial osteomyelitis (CNO)

Chronic nonbacterial osteomyelitis (CNO), an autoinflammatory bone disease primarily affecting children, can cause pain, hyperostosis and fractures, affecting quality-of-life and psychomotor development. This study investigated CNO-associated variants in P2RX7, encoding for the ATP-dependent trans-membrane K+ channel P2X7, and their effects on NLRP3 inflammasome assembly. Whole exome sequencing in two related transgenerational CNO patients, and target sequencing of P2RX7 in a large CNO cohort (N = 190) were conducted. Results were compared with publicly available datasets and regional controls (N = 1873). Findings were integrated with demographic and clinical data. Patient-derived monocytes and genetically modified THP-1 cells were used to investigate potassium flux, inflammasome assembly, pyroptosis, and cytokine release. Rare presumably damaging P2RX7 variants were identified in two related CNO patients. Targeted P2RX7 sequencing identified 62 CNO patients with rare variants (32.4%), 11 of which (5.8%) carried presumably damaging variants (MAF <1%, SIFT "deleterious", Polyphen "probably damaging", CADD >20). This compared to 83 of 1873 controls (4.4%), 36 with rare and presumably damaging variants (1.9%). Across the CNO cohort, rare variants unique to one (Median: 42 versus 3.7) or more (≤11 patients) participants were over-represented when compared to 190 randomly selected controls. Patients with rare damaging variants more frequently experienced gastrointestinal symptoms and lymphadenopathy while having less spinal, joint and skin involvement (psoriasis). Monocyte-derived macrophages from patients, and genetically modified THP-1-derived macrophages reconstituted with CNO-associated P2RX7 variants exhibited altered potassium flux, inflammasome assembly, IL-1β and IL-18 release, and pyroptosis. Damaging P2RX7 variants occur in a small subset of CNO patients, and rare P2RX7 variants may represent a CNO risk factor. Observations argue for inflammasome inhibition and/or cytokine blockade and may allow future patient stratification and individualized care.

GWAS meta-analysis of psoriasis identifies new susceptibility alleles impacting disease mechanisms and therapeutic targets

Psoriasis is a common, debilitating immune-mediated skin disease. Genetic studies have identified biological mechanisms of psoriasis risk, including those targeted by effective therapies. However, the genetic liability to psoriasis is not fully explained by variation at robustly identified risk loci. To move towards a saturation map of psoriasis susceptibility we meta-analysed 18 GWAS comprising 36,466 cases and 458,078 controls and identified 109 distinct psoriasis susceptibility loci, including 45 that have not been previously reported. These include susceptibility variants at loci in which the therapeutic targets IL17RA and AHR are encoded, and deleterious coding variants supporting potential new drug targets (including in STAP2, CPVL and POU2F3). We conducted a transcriptome-wide association study to identify regulatory effects of psoriasis susceptibility variants and cross-referenced these against single cell expression profiles in psoriasis-affected skin, highlighting roles for the transcriptional regulation of haematopoietic cell development and epigenetic modulation of interferon signalling in psoriasis pathobiology.

Abstract 4587: Vascular-derived signature imprinted by tumor microenvironment-dependent transcriptional memory predicts colon cancer prognosis

Background: Antiangiogenic therapy is part of the guideline therapy of colorectal cancer (CRC). Surprisingly, the impact of the tumor microenvironment (TME) on tumor vessel endothelial cells (TECs) is largely unclear. The aim of this study was to investigate the presence of a TME-dependent transcriptional memory in TECs isolated from patients and to exploit it to retrieve signatures that characterize TECs in different intratumoral TMEs with impact on patient outcome.

Methods: ECs from tumor and normal colon tissues, PBMCs and tumor cells were isolated from CRC patients with different prognostic TMEs. Cells were analyzed by qPCR, immunocytochemistry and multi-omics (transcriptomics, EPICmethylation chips, exome sequencing). Integrative bioinformatics was used to identify TME-dependent memory genes predicting prognosis and scRNASeq to validate endothelial gene expression in CRC tissues.

Results: Ultrapure TECs were isolated from CRC with different prognostic TMEs (Th1 vs. non-Th1) and systematically compared by multi-omics. A transcriptional memory differentiating the respective TEC groups was identified. This in vivo imprinted transcriptional memory was preferentially regulated by epigenetic DNA methylation but not by genomic alterations and was different from an in vitro primed transcriptional memory to IFN-γ. Moreover, it was specific for TECs and not observed in CAFs. With integrative bioinformatics a TME-dependent memory signature was extracted and its expression in TECs in CRC tissues was confirmed by scRNASeq. Notably, the identified signature predicted the prognosis of CRC patients.

Conclusion: We identified a tumor vessel-derived TME-dependent transcriptional memory signature that was manifested by epigenetic mechanisms and allowed tumor vessel-based prediction of CRC patients prognosis.

Citation Format: Elisabeth Naschberger, Maximilian Fuchs, Nicholas Dickel, Meik Kunz, Charles G. Anchang, Richard Demmler, Bernt Popp, Arif B. Ekici, Steffen Uebe, Carol I. Geppert, Claudia Günther, Susanne Merkel, Vera S. Schellerer, Michael Stürzl. Vascular-derived signature imprinted by tumor microenvironment-dependent transcriptional memory predicts colon cancer prognosis. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 4587.

Interspecies Single-Cell RNA-Seq Analysis Reveals the Novel Trajectory of Osteoclast Differentiation and Therapeutic Targets

Bone turnover is finely tuned by cells in the bone milieu, including osteoblasts, osteoclasts, and osteocytes. Osteoclasts are multinucleated giant cells with a bone-resorbing function that play a critical role in regulating skeletal homeostasis. Osteoclast differentiation is characterized by dramatic changes in morphology and gene expression following receptor activator of nuclear factor-kappa-Β ligand (RANKL) stimulation. We performed single-cell RNA-sequencing analyses of human and murine osteoclast-lineage cells (OLCs) and found that OLCs in the mitotic phase do not differentiate into mature osteoclasts. We also identified a guanosine triphosphatase (GTPase) family member, RAB38, as a highly expressed molecule in both human and murine osteoclast clusters; RAB38 gene expression is associated with dynamic changes in histone modification and transcriptional regulation. Silencing Rab38 expression by using short hairpin RNA (shRNA) inhibited osteoclast differentiation and maturation. In summary, we established an integrated fate map of human and murine osteoclastogenesis; this will help identify therapeutic targets in bone diseases. © 2022 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Abstract 3802: Micromilieu-dependent transcriptional memory of tumor vessel endothelial cells in colorectal cancer

Purpose of the study: Tumor microenvironment (TME)-induced plasticity of fibroblasts and immune cells are known to contribute to cancer pathogenesis. In contrast, only little is known about the existence, stability, functions and prognostic potential of different plastic phenotypes of tumor vessel endothelial cells (TECs) in cancer. Recently, transcriptional memory effects have been observed in different cell lines after cytokine stimulation in vitro. From this we developed the hypothesis that isolated TEC may exhibit TME-dependent transcriptional memories, which may provide novel insights into the intratumoral effects on TEC plasticity in cancer. To investigate this, we isolated ultrapure TECs from colorectal cancer (CRC) with prognostically favorable (Th-1-like) or worse (non-Th1-like) immune-micro-environments and analyzed the prescence of TME-associated transcriptional memory patterns.

Methods: TECs and as controls normal colon endothelial cells (NECs), PBMCs and tumor cells (isolated by laser-microdissection) from the same CRC patients (n=75) were isolated by FACS-sorting. Ultra high purity (98%) of TECs and NECs was validated by qPCR and immunocytochemistry. Quality validated complete sample sets from 12 patients (Th-1-like, n=6, non-Th1-like, n=6) were compared using exome sequencing, Clariome S chips and EPICmethylation. Integrative bioinformatics was used to identify differential gene sets and to predict their association with prognosis.

Results: Differential transcriptional memory patterns were detected in TECs derived from CRC with different prognostic TMEs. Of note, the aggressiveness of the originating tumors was reflected by the functions of the vascular transcriptional memory genes. The transcriptional memory was imprinted by epigenetic DNA methylation and not due to acquired somatic variants in the genome. These results could be successfully confirmed by alternative methods. Most importantly, integrative bioinformatics analyses retrieved seven vascular memory genes (VMG), which were expressed in the TECs in a TME-dependent manner and predicted the outcome of patients in an independent cohort of CRC patients.

Conclusion: Our findings demonstrate for the first time a TME-dependent transcriptional memory of TECs, providing a clinically relevant indicator of intra-tumoral TME-dependent vascular plasticity with perspectives for patient stratification in vessel-directed tumor therapies.

Citation Format: Michael Stürzl, Maximilian Fuchs, Meik Kunz, Nicholas Dickel, Bernt Popp, Steffen Uebe, Arif B. Ekici, Carol I. Geppert, Susanne Merkel, Vera S. Schellerer, Elisabeth Naschberger. Micromilieu-dependent transcriptional memory of tumor vessel endothelial cells in colorectal cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 3802.

SRD5A3-CDG: Twins with an intragenic tandem duplication

Steroid 5α-reductase type 3 congenital disorder of glycosylation (SRD5A3-CDG) is a rare metabolic disease mainly characterized by psychomotor disability, visual impairment, and variable eye malformations caused by bi-allelic pathogenic variants in SRD5A3. So far, only 23 distinct mutations were described. Exome sequencing in 32-year old monozygotic male twins revealed only the heterozygous splice variant c.562+3delG in SRD5A3, but no second variant. The twins presented with psychomotor deficit and a complex eye disease including retinal dystrophy, pallor of the papilla, nystagmus, and strabismus suggestive of SRD5A3-CDG. Only when applying exome-based copy number analysis, we identified as a second compound heterozygous variant a previously not reported tandem duplication of exons 2-4 in SRD5A3. Next to the typical skeletal anomalies of SRD5A3-CDG such as kyphosis and scoliosis, extension deficits of the proximal interphalangeal (PIP) joints IV were observed. Since similar contractures were described once in a patient with SRD5A3-CDG, we suggest that this rare symptom is possibly associated with SRD5A3-CDG. Our findings further expand the mutational and clinical spectrum of SRD5A3-CDG and emphasize the importance of an intragenic copy number analysis in patients with strong clinical suspicion of SRD5A3-CDG and only one detectable sequence variant.

Transcriptomes of MPO-deficient patients with generalized pustular psoriasis reveals expansion of CD4+ cytotoxic T cells and an involvement of the complement system

Generalized pustular psoriasis (GPP) is a severe psoriatic subtype characterized by epidermal neutrophil infiltration. Although variants in IL36RN and MPO have been shown to affect immune cells, a systematic analysis of neutrophils and peripheral blood mononuclear cells (PBMCs) subsets and their differential gene expression dependent on MPO genotypes was not performed yet. We assessed transcriptomes of MPO-deficient patients using single cell RNA-sequencing (scRNAseq) of PBMCs and RNA-sequencing of neutrophils in stable disease state. Cell type annotation by multimodal reference mapping of scRNAseq data was verified by flow cytometry of surface and intracellular markers; proportions of CD4⁺ cytotoxic T-lymphocytes (CTLs) and other CD4⁺ effector cells were increased in GPP, while frequencies of naïve CD4⁺ T cells were significantly lower. The expression of FGFBP2 marking CD4⁺ CTLs and CD8⁺ effector memory T-cells (TEMs) was elevated in GPP patients with disease-contributing variants compared to non-carriers (p=0.0015). In neutrophils, differentially expressed genes (DEGs) were significantly enriched in genes of the classical complement activation pathway. Future studies assessing affected cell-types and pathways will show their contribution to GPP's pathogenesis, and indicate whether findings can be transferred to the acute epidermal situation and whether depletion or inactivation of CD4⁺ CTLs may be a reasonable therapeutic approach.

BDV Syndrome: An Emerging Syndrome With Profound Obesity and Neurodevelopmental Delay Resembling Prader-Willi Syndrome

CONTEXT

CPE encodes carboxypeptidase E, an enzyme that converts proneuropeptides and propeptide hormones to bioactive forms. It is widely expressed in the endocrine and central nervous system. To date, 4 individuals from 2 families with core clinical features including morbid obesity, neurodevelopmental delay, and hypogonadotropic hypogonadism, harboring biallelic loss-of-function (LoF) CPE variants, have been reported.

OBJECTIVE

We describe 4 affected individuals from 3 unrelated consanguineous families, 2 siblings of Syrian, 1 of Egyptian, and 1 of Pakistani descent, all harboring novel homozygous CPE LoF variants.

METHODS

After excluding Prader-Willi syndrome (PWS), exome sequencing was performed in both Syrian siblings. The variants identified in the other 2 individuals were reported as research variants in a large-scale exome study and in the ClinVar database. Computational modeling of all possible missense alterations allowed assessing CPE tolerance to missense variants.

RESULTS

All affected individuals were severely obese with neurodevelopmental delay and other endocrine anomalies. Three individuals from 2 families shared the same CPE homozygous truncating variant c.361C > T, p.(Arg121*), while the fourth carried the c.994del, p.(Ser333Alafs*22) variant. Comparison of clinical features with previously described cases and standardization according to the Human Phenotype Ontology terms indicated a recognizable clinical phenotype, which we termed Blakemore-Durmaz-Vasileiou (BDV) syndrome. Computational analysis indicated high conservation of CPE domains and intolerance to missense changes.

CONCLUSION

Biallelic truncating CPE variants are associated with BDV syndrome, a clinically recognizable monogenic recessive syndrome with childhood-onset obesity, neurodevelopmental delay, hypogonadotropic hypogonadism, and hypothyroidism. BDV syndrome resembles PWS. Our findings suggest missense variants may also be clinically relevant.

Experimental Epileptogenesis in a Cell Culture Model of Primary Neurons from Rat Brain: A Temporal Multi-Scale Study

Understanding seizure development requires an integrated knowledge of different scales of organization of epileptic networks. We developed a model of “epilepsy-in-a-dish” based on dissociated primary neuronal cells from neonatal rat hippocampus. We demonstrate how a single application of glutamate stimulated neurons to generate spontaneous synchronous spiking activity with further progression into spontaneous seizure-like events after a distinct latency period. By computational analysis, we compared the observed neuronal activity in vitro with intracranial electroencephalography (EEG) data recorded from epilepsy patients and identified strong similarities, including a related sequence of events with defined onset, progression, and termination. Next, a link between the neurophysiological changes with network composition and cellular structure down to molecular changes was established. Temporal development of epileptiform network activity correlated with increased neurite outgrowth and altered branching, increased ratio of glutamatergic over GABAergic synapses, and loss of calbindin-positive interneurons, as well as genome-wide alterations in DNA methylation. Differentially methylated genes were engaged in various cellular activities related to cellular structure, intracellular signaling, and regulation of gene expression. Our data provide evidence that a single short-term excess of glutamate is sufficient to induce a cascade of events covering different scales from molecule- to network-level, all of which jointly contribute to seizure development.

Manifestation of epilepsy in a patient with EED-related overgrowth (Cohen-Gibson syndrome)

Cohen-Gibson syndrome is a rare genetic disorder, characterized by fetal or early childhood overgrowth and mild to severe intellectual disability. It is caused by heterozygous aberrations in EED, which encodes an evolutionary conserved polycomb group (PcG) protein that forms the polycomb repressive complex-2 (PRC2) together with EZH2, SUZ12, and RBBP7/4. In total, 11 affected individuals with heterozygous pathogenic variants in EED were reported, so far. All variants affect a few key residues within the EED WD40 repeat domain. By trio exome sequencing, we identified the heterozygous missense variant c.581A > G, p.(Asn194Ser) in exon 6 of the EED-gene in an individual with moderate intellectual disability, overgrowth, and epilepsy. The same pathogenic variant was detected in 2 of the 11 previously reported cases. Epilepsy, however, was only diagnosed in one other individual with Cohen-Gibson syndrome before. Our findings further confirm that the WD40 repeat domain represents a mutational hotspot; they also expand the clinical spectrum of Cohen-Gibson syndrome and highlight the clinical variability even in individuals with the same pathogenic variant. Furthermore, they indicate a possible association between Cohen-Gibson syndrome and epilepsy.

The protective variant rs7173049 at LOXL1 locus impacts on retinoic acid signaling pathway in pseudoexfoliation syndrome

LOXL1 (lysyl oxidase-like 1) has been identified as the major effect locus in pseudoexfoliation (PEX) syndrome, a fibrotic disorder of the extracellular matrix and frequent cause of chronic open-angle glaucoma. However, all known PEX-associated common variants show allele effect reversal in populations of different ancestry, casting doubt on their biological significance. Based on extensive LOXL1 deep sequencing, we report here the identification of a common non-coding sequence variant, rs7173049A>G, located downstream of LOXL1, consistently associated with a decrease in PEX risk (odds ratio, OR = 0.63; P = 6.33 × 10⁻³¹) in nine different ethnic populations. We provide experimental evidence for a functional enhancer-like regulatory activity of the genomic region surrounding rs7173049 influencing expression levels of ISLR2 (immunoglobulin superfamily containing leucine-rich repeat protein 2) and STRA6 [stimulated by retinoic acid (RA) receptor 6], apparently mediated by allele-specific binding of the transcription factor thyroid hormone receptor beta. We further show that the protective rs7173049-G allele correlates with increased tissue expression levels of ISLR2 and STRA6 and that both genes are significantly downregulated in tissues of PEX patients together with other key components of the STRA6 receptor-driven RA signaling pathway. siRNA-mediated downregulation of RA signaling induces upregulation of LOXL1 and PEX-associated matrix genes in PEX-relevant cell types. These data indicate that dysregulation of STRA6 and impaired retinoid metabolism are involved in the pathophysiology of PEX syndrome and that the variant rs7173049-G, which represents the first common variant at the broad LOXL1 locus without allele effect reversal, mediates a protective effect through upregulation of STRA6 in ocular tissues.

EIF3F-related neurodevelopmental disorder: refining the phenotypic and expanding the molecular spectrum

BACKGROUND

An identical homozygous missense variant in EIF3F, identified through a large-scale genome-wide sequencing approach, was reported as causative in nine individuals with a neurodevelopmental disorder, characterized by variable intellectual disability, epilepsy, behavioral problems and sensorineural hearing-loss. To refine the phenotypic and molecular spectrum of EIF3F-related neurodevelopmental disorder, we examined independent patients.

RESULTS

21 patients were homozygous and one compound heterozygous for c.694T>G/p.(Phe232Val) in EIF3F. Haplotype analyses in 15 families suggested that c.694T>G/p.(Phe232Val) was a founder variant. All affected individuals had developmental delays including delayed speech development. About half of the affected individuals had behavioral problems, altered muscular tone, hearing loss, and short stature. Moreover, this study suggests that microcephaly, reduced sensitivity to pain, cleft lip/palate, gastrointestinal symptoms and ophthalmological symptoms are part of the phenotypic spectrum. Minor dysmorphic features were observed, although neither the individuals' facial nor general appearance were obviously distinctive. Symptoms in the compound heterozygous individual with an additional truncating variant were at the severe end of the spectrum in regard to motor milestones, speech delay, organic problems and pre- and postnatal growth of body and head, suggesting some genotype-phenotype correlation.

CONCLUSIONS

Our study refines the phenotypic and expands the molecular spectrum of EIF3F-related syndromic neurodevelopmental disorder.

Genome-wide meta-analysis identifies 127 open-angle glaucoma loci with consistent effect across ancestries

Primary open-angle glaucoma (POAG), is a heritable common cause of blindness world-wide. To identify risk loci, we conduct a large multi-ethnic meta-analysis of genome-wide association studies on a total of 34,179 cases and 349,321 controls, identifying 44 previously unreported risk loci and confirming 83 loci that were previously known. The majority of loci have broadly consistent effects across European, Asian and African ancestries. Cross-ancestry data improve fine-mapping of causal variants for several loci. Integration of multiple lines of genetic evidence support the functional relevance of the identified POAG risk loci and highlight potential contributions of several genes to POAG pathogenesis, including SVEP1, RERE, VCAM1, ZNF638, CLIC5, SLC2A12, YAP1, MXRA5, and SMAD6. Several drug compounds targeting POAG risk genes may be potential glaucoma therapeutic candidates.

Network- and systems-based re-engineering of dendritic cells with non-coding RNAs for cancer immunotherapy

Dendritic cells (DCs) are professional antigen-presenting cells that induce and regulate adaptive immunity by presenting antigens to T cells. Due to their coordinative role in adaptive immune responses, DCs have been used as cell-based therapeutic vaccination against cancer. The capacity of DCs to induce a therapeutic immune response can be enhanced by re-wiring of cellular signalling pathways with microRNAs (miRNAs). Methods: Since the activation and maturation of DCs is controlled by an interconnected signalling network, we deploy an approach that combines RNA sequencing data and systems biology methods to delineate miRNA-based strategies that enhance DC-elicited immune responses. Results: Through RNA sequencing of IKKβ-matured DCs that are currently being tested in a clinical trial on therapeutic anti-cancer vaccination, we identified 44 differentially expressed miRNAs. According to a network analysis, most of these miRNAs regulate targets that are linked to immune pathways, such as cytokine and interleukin signalling. We employed a network topology-oriented scoring model to rank the miRNAs, analysed their impact on immunogenic potency of DCs, and identified dozens of promising miRNA candidates, with miR-15a and miR-16 as the top ones. The results of our analysis are presented in a database that constitutes a tool to identify DC-relevant miRNA-gene interactions with therapeutic potential (https://www.synmirapy.net/dc-optimization). Conclusions: Our approach enables the systematic analysis and identification of functional miRNA-gene interactions that can be experimentally tested for improving DC immunogenic potency.

Myeloperoxidase Modulates Inflammation in Generalized Pustular Psoriasis and Additional Rare Pustular Skin Diseases

Generalized pustular psoriasis (GPP) is a severe multi-systemic inflammatory disease characterized by neutrophilic pustulosis and triggered by pro-inflammatory IL-36 cytokines in skin. While 19%–41% of affected individuals harbor bi-allelic mutations in IL36RN, the genetic cause is not known in most cases. To identify and characterize new pathways involved in the pathogenesis of GPP, we performed whole-exome sequencing in 31 individuals with GPP and demonstrated effects of mutations in MPO encoding the neutrophilic enzyme myeloperoxidase (MPO). We discovered eight MPO mutations resulting in MPO -deficiency in neutrophils and monocytes. MPO mutations, primarily those resulting in complete MPO deficiency, cumulatively associated with GPP (p = 1.85E−08; OR = 6.47). The number of mutant MPO alleles significantly differed between 82 affected individuals and >4,900 control subjects (p = 1.04E−09); this effect was stronger when including IL36RN mutations (1.48E−13) and correlated with a younger age of onset (p = 0.0018). The activity of four proteases, previously implicated as activating enzymes of IL-36 precursors, correlated with MPO deficiency. Phorbol-myristate-acetate-induced formation of neutrophil extracellular traps (NETs) was reduced in affected cells (p = 0.015), and phagocytosis assays in MPO-deficient mice and human cells revealed altered neutrophil function and impaired clearance of neutrophils by monocytes (efferocytosis) allowing prolonged neutrophil persistence in inflammatory skin. MPO mutations contribute significantly to GPP’s pathogenesis. We implicate MPO as an inflammatory modulator in humans that regulates protease activity and NET formation and modifies efferocytosis. Our findings indicate possible implications for the application of MPO inhibitors in cardiovascular diseases. MPO and affected pathways represent attractive targets for inducing resolution of inflammation in neutrophil-mediated skin diseases.

A case of severe autosomal recessive spinocerebellar ataxia type 18 with a novel nonsense variant in GRID2

Autosomal-recessive spinocerebellar ataxia type 18 (SCAR18) is a rare neurologic disorder. It is caused by bi-allelic aberrations in the GRID2 gene, encoding an ionotropic glutamate receptor. In total, 20 affected individuals with mainly homozygous/compound heterozygous intragenic deletions/duplications, two different missense variants and one nonsense variant in GRID2 have been reported, so far. SCAR18 is characterized by delayed psychomotor development, intellectual disability, severely impaired gait due to cerebellar ataxia, ocular movement abnormalities, and cerebellar atrophy in brain imaging. By trio exome sequencing, we now identified a novel homozygous nonsense variant (c.568C > T; p.Gln190*) in GRID2 in a four year old female from a consanguineous family who presented with a particularly severe manifestation of SCAR18. The girl was born after an uneventful pregnancy and showed early-onset, profoundly delayed psychomotor development with no achieved psychomotor milestones at age 4 years. Additionally, she presented with severe muscular hypotonia, progressive truncal and appendicular ataxia, binocular vertical nystagmus, central hearing loss and incomplete loss of sight. She was dystrophic, interacted only very little and had behavioral anomalies such as eating hair and bruxism. Brain imaging showed cerebellar hypoplasia, extended cerebrospinal fluid spaces and beginning reduction of cerebral volume. Our findings further delineate the mutational and clinical spectrum of GRID2-associated spinocerebellar ataxia type 18 and indicate that homozygous nonsense variants are possibly associated with the severe end of the SCAR18 phenotypic spectrum.

Prenatal diagnosis of HNF1B-associated renal cysts: Is there a need to differentiate intragenic variants from 17q12 microdeletion syndrome?

OBJECTIVE

17q12 microdeletions containing HNF1B and intragenic variants within this gene are associated with variable developmental, endocrine, and renal anomalies, often already noted prenatally as hyperechogenic/cystic kidneys. Here, we describe prenatal and postnatal phenotypes of seven individuals with HNF1B aberrations and compare their clinical and genetic data to those of previous studies.

METHODS

Prenatal sequencing and postnatal chromosomal microarray analysis were performed in seven individuals with renal and/or neurodevelopmental phenotypes. We evaluated HNF1B-related clinical features from 82 studies and reclassified 192 reported intragenic HNF1B variants.

RESULTS

In a prenatal case, we identified a novel in-frame deletion p.(Gly239del) within the HNF1B DNA-binding domain, a mutational hot spot as demonstrated by spatial clustering analysis and high computational prediction scores. The six postnatally diagnosed individuals harbored 17q12 microdeletions. Literature screening revealed variable reporting of HNF1B-associated clinical traits. Overall, both mutation groups showed a high phenotypic heterogeneity. The reclassification of all previously reported intragenic HNF1B variants provided an up-to-date overview of the mutational spectrum.

CONCLUSIONS

We highlight the value of prenatal HNF1B screening in renal developmental diseases. Standardized clinical reporting and systematic classification of HNF1B variants are necessary for a more accurate risk quantification of prenatal and postnatal clinical features, improving genetic counseling and prenatal decision making.

Dissecting TSC2-mutated renal and hepatic angiomyolipomas in an individual with ARID1B-associated intellectual disability

Background

Several subunits of the SWI/SNF chromatin remodeling complex are implicated in both cancer and neurodevelopmental disorders (NDD). Though there is no clinical evidence for an increased tumor risk in individuals with NDDs due to germline mutations in most of these genes so far, this has been repeatedly proposed and discussed. A young woman with NDD due to a de novo mutation in ARID1B now presented with a large renal (> 19 cm in diameter) and multiple hepatic angiomyolipomas (AMLs) but no other signs of tuberous sclerosis complex.

Methods

We analyzed tumor and healthy tissue samples with exome and panel sequencing.

Results

Additionally to the previously known, germline ARID1B variant we identified a post-zygotic truncating TSC2 variant in both renal and hepatic AMLs but not in any of the healthy tissues. We did not detect any further, obvious tumor driver events. The identification of a passenger variant in SIPA1L3 in both AMLs points to a common clonal origin. Metastasis of the renal AML into the liver is unlikely on the basis of discordant histopathological features. Our findings therefore point to very low-grade mosaicism for the TSC2 variant, possibly in a yet unknown mesenchymal precursor cell that expanded clonally during tumor development. A possible contribution of the germline ARID1B variant to the tumorigenesis remains unclear but cannot be excluded given the absence of any other evident tumor drivers in the AMLs.

Conclusion

This unique case highlights the blurred line between tumor genetics and post-zygotic events that can complicate exact molecular diagnoses in patients with rare manifestations. It also demonstrates the relevance of multiple disorders in a single individual, the challenges of detecting low-grade mosaicisms, and the importance of proper diagnosis for treatment and surveillance.

Electronic supplementary material

The online version of this article (10.1186/s12885-019-5633-1) contains supplementary material, which is available to authorized users.

Transcription factor Fra-1 targets arginase-1 to enhance macrophage-mediated inflammation in arthritis

The polarization of macrophages is regulated by transcription factors such as nuclear factor kappa B (NF-κB) and activator protein 1 (AP-1). In this manuscript, we delineated the role of the transcription factor Fos-related antigen 1 (Fra-1) during macrophage activation and development of arthritis. Network level interaction analysis of microarray data derived from Fra-1- or Fra-2-deficient macrophages revealed a central role of Fra-1, but not of Fra-2 in orchestrating the expression of genes related to wound response, toll-like receptor activation and interleukin signaling. Chromatin-immunoprecipitation (ChIP)-sequencing and standard ChIP analyses of macrophages identified arginase 1 (Arg1) as a target of Fra-1. Luciferase reporter assays revealed that Fra-1 down-regulated Arg1 expression by direct binding to the promoter region. Using macrophage-specific Fra-1- or Fra-2- deficient mice, we observed an enhanced expression and activity of Arg1 and a reduction of arthritis in the absence of Fra-1, but not of Fra-2. This phenotype was reversed by treatment with the arginase inhibitor Nω-hydroxy-nor-L-arginine, while ʟ-arginine supplementation increased arginase activity and alleviated arthritis, supporting the notion that reduced arthritis in macrophage-specific Fra-1-deficient mice resulted from enhanced Arg1 expression and activity. Moreover, patients with active RA showed increased Fra-1 expression in the peripheral blood and elevated Fra-1 protein in synovial macrophages compared to RA patients in remission. In addition, the Fra-1/ARG1 ratio in synovial macrophages was related to RA disease activity. In conclusion, these data suggest that Fra-1 orchestrates the inflammatory state of macrophages by inhibition of Arg1 expression and thereby impedes the resolution of inflammation.

Evolutionary conserved networks of human height identify multiple Mendelian causes of short stature

Height is a heritable and highly heterogeneous trait. Short stature affects 3% of the population and in most cases is genetic in origin. After excluding known causes, 67% of affected individuals remain without diagnosis. To identify novel candidate genes for short stature, we performed exome sequencing in 254 unrelated families with short stature of unknown cause and identified variants in 63 candidate genes in 92 (36%) independent families. Based on systematic characterization of variants and functional analysis including expression in chondrocytes, we classified 13 genes as strong candidates. Whereas variants in at least two families were detected for all 13 candidates, two genes had variants in 6 (UBR4) and 8 (LAMA5) families, respectively. To facilitate their characterization, we established a clustered network of 1025 known growth and short stature genes, which yielded 29 significantly enriched clusters, including skeletal system development, appendage development, metabolic processes, and ciliopathy. Eleven of the candidate genes mapped to 21 of these clusters, including CPZ, EDEM3, FBRS, IFT81, KCND1, PLXNA3, RASA3, SLC7A8, UBR4, USP45, and ZFHX3. Fifty additional growth-related candidates we identified await confirmation in other affected families. Our study identifies Mendelian forms of growth retardation as an important component of idiopathic short stature.

TRIM28 haploinsufficiency predisposes to Wilms tumor

Two percent of patients with Wilms tumors have a positive family history. In many of these cases the genetic cause remains unresolved. By applying germline exome sequencing in two families with two affected individuals with Wilms tumors, we identified truncating mutations in TRIM28. Subsequent mutational screening of germline and tumor DNA of 269 children affected by Wilms tumor was performed, and revealed seven additional individuals with germline truncating mutations, and one individual with a somatic truncating mutation in TRIM28. TRIM28 encodes a complex scaffold protein involved in many different processes, including gene silencing, DNA repair and maintenance of genomic integrity. Expression studies on mRNA and protein level showed reduction of TRIM28, confirming a loss-of-function effect of the mutations identified. The tumors showed an epithelial-type histology that stained negative for TRIM28 by immunohistochemistry. The tumors were bilateral in six patients, and 10/11 tumors are accompanied by perilobar nephrogenic rests. Exome sequencing on eight tumor DNA samples from six individuals showed loss-of-heterozygosity (LOH) of the TRIM28-locus by mitotic recombination in seven tumors, suggesting that TRIM28 functions as a tumor suppressor gene in Wilms tumor development. Additionally, the tumors showed very few mutations in known Wilms tumor driver genes, suggesting that loss of TRIM28 is the main driver of tumorigenesis. In conclusion, we identified heterozygous germline truncating mutations in TRIM28 in 11 children with mainly epithelial-type Wilms tumors, which become homozygous in tumor tissue. These data establish TRIM28 as a novel Wilms tumor predisposition gene, acting as a tumor suppressor gene by LOH.

The mutational and phenotypic spectrum of TUBA1A-associated tubulinopathy

BACKGROUND

The TUBA1A-associated tubulinopathy is clinically heterogeneous with brain malformations, microcephaly, developmental delay and epilepsy being the main clinical features. It is an autosomal dominant disorder mostly caused by de novo variants in TUBA1A.

RESULTS

In three individuals with developmental delay we identified heterozygous de novo missense variants in TUBA1A using exome sequencing. While the c.1307G > A, p.(Gly436Asp) variant was novel, the two variants c.518C > T, p.(Pro173Leu) and c.641G > A, p.(Arg214His) were previously described. We compared the variable phenotype observed in these individuals with a carefully conducted review of the current literature and identified 166 individuals, 146 born and 20 fetuses with a TUBA1A variant. In 107 cases with available clinical information we standardized the reported phenotypes according to the Human Phenotype Ontology. The most commonly reported features were developmental delay (98%), anomalies of the corpus callosum (96%), microcephaly (76%) and lissencephaly (agyria-pachygyria) (70%), although reporting was incomplete in the different studies. We identified a total of 121 specific variants, including 15 recurrent ones. Missense variants cluster in the C-terminal region around the most commonly affected amino acid position Arg402 (13.3%). In a three-dimensional protein model, 38.6% of all disease-causing variants including those in the C-terminal region are predicted to affect the binding of microtubule-associated proteins or motor proteins. Genotype-phenotype analysis for recurrent variants showed an overrepresentation of certain clinical features. However, individuals with these variants are often reported in the same publication.

CONCLUSIONS

With 166 individuals, we present the most comprehensive phenotypic and genotypic standardized synopsis for clinical interpretation of TUBA1A variants. Despite this considerable number, a detailed genotype-phenotype characterization is limited by large inter-study variability in reporting.

PU.1 controls fibroblast polarization and tissue fibrosis

Fibroblasts are polymorphic cells with pleiotropic roles in organ morphogenesis, tissue homeostasis and immune responses. In fibrotic diseases, fibroblasts synthesize abundant amounts of extracellular matrix which lead to scaring and organ failure. In contrast, the hallmark feature of fibroblasts in arthritis is matrix degradation by the release of metalloproteinases and degrading enzymes, and subsequent tissue destruction. The mechanisms driving these functionally opposing pro-fibrotic and pro-inflammatory phenotypes of fibroblasts are enigmatic. We identified the transcription factor PU.1 as an essential orchestrator of the pro-fibrotic gene expression program. The interplay between transcriptional and post-transcriptional mechanisms which normally control PU.1 expression is perturbed in various fibrotic diseases, resulting in upregulation of PU.1, induction of fibrosis-associated gene sets, and a phenotypic switch in matrix-producing pro-fibrotic fibroblasts. In contrast, pharmacological and genetic inactivation of PU.1 disrupts the fibrotic network and enables re-programming of fibrotic fibroblasts into resting fibroblasts with regression of fibrosis in different organs.

A biallelic truncating AEBP1 variant causes connective tissue disorder in two siblings

Biallelic variants in the AEBP1 gene cause a novel autosomal-recessive connective tissue disorder (CTD) reminiscent of Ehlers-Danlos Syndrome (EDS). The four previously reported individuals show considerable clinical variability. Unbiased high-throughput sequencing enables the rapid identification of additional cases for such rare entities. We identified the homozygous nonsense variant c.917dup, p.Tyr306* in AEBP1 using clinical exome sequencing in a female individual with previously unsolved CTD. Segregation testing confirmed homozygosity in the clinically affected brother and heterozygous carrier status in the healthy mother. Chromosomal microarray showed that the variant lies in a run of homozygosity, suggesting a common origin of this genomic segment. RT-PCR analysis in the mother revealed a monoallelic expression of the normal transcript supporting a nonsense-mediated mRNA decay and functional nullizygosity as disease mechanism. We describe two individuals from a fourth family with AEBP1-associated CTD. Our results further verify that autosomal-recessive inherited LOF variants in the AEBP1 gene cause clinical features of different EDS subtypes, but also of the marfanoid spectrum. As identification of further individuals is necessary to inform the clinical characterization, we stress the added value of exome sequencing for such rare diseases.

Need for high-resolution Genetic Analysis in iPSC: Results and Lessons from the ForIPS Consortium

Genetic integrity of induced pluripotent stem cells (iPSCs) is essential for their validity as disease models and for potential therapeutic use. We describe the comprehensive analysis in the ForIPS consortium: an iPSC collection from donors with neurological diseases and healthy controls. Characterization included pluripotency confirmation, fingerprinting, conventional and molecular karyotyping in all lines. In the majority, somatic copy number variants (CNVs) were identified. A subset with available matched donor DNA was selected for comparative exome sequencing. We identified single nucleotide variants (SNVs) at different allelic frequencies in each clone with high variability in mutational load. Low frequencies of variants in parental fibroblasts highlight the importance of germline samples. Somatic variant number was independent from reprogramming, cell type and passage. Comparison with disease genes and prediction scores suggest biological relevance for some variants. We show that high-throughput sequencing has value beyond SNV detection and the requirement to individually evaluate each clone.

Addition of triple negativity of breast cancer as an indicator for germline mutations in predisposing genes increases sensitivity of clinical selection criteria

Background

Breast cancer is the most common cancer in women. 12–15% of all tumors are triple-negative breast cancers (TNBC). So far, TNBC has been mainly associated with mutations in BRCA1. The presence of other predisposing genes seems likely since DNA damage repair is a complex process that involves several genes. Therefore we investigated if mutations in other genes are involved in cancer development and whether TNBC is an additional indicator of mutational status besides family history and age of onset.

Methods

We performed a germline panel-based screening of 10 high and low-moderate penetrance breast cancer susceptibility genes (BRCA1, BRCA2, ATM, CDH1, CHEK2, NBN, PALB2, RAD51C, RAD51D and TP53) in 229 consecutive individuals affected with TNBC unselected for age, family history or bilateral disease. Within this cohort we compared the number of mutation carriers fulfilling clinical selection criteria with the total number of carriers identified.

Results

Age at diagnosis ranged from 23 to 80 years with an average age of 50.2 years. In 57 women (24.9%) we detected a pathogenic mutation, with a higher frequency (29.7%) in the group manifesting cancer before 60 years. Deleterious BRCA1 mutations occurred in 14.8% of TNBC patients. These were predominantly recurrent frameshift mutations (24/34, 70.6%). Deleterious BRCA2 mutations occurred in 5.7% of patients, all but one (c.1813dupA) being unique.

While no mutations were found in CDH1 and TP53, 10 mutations were detected in one of the six other predisposition genes. Remarkably, neither of the ATM, RAD51D, CHEK2 and PALB2 mutation carriers had a family history. Furthermore, patients with non-BRCA1/2 mutations were not significantly younger than mutation negative women (p = 0.3341). Most importantly, among the 57 mutation carriers, ten (17.5%) would be missed using current clinical testing criteria including five (8%) with BRCA1/2 mutations.

Conclusions

In summary, our data confirm and expand previous studies of a high frequency of germline mutations in genes associated with ineffective repair of DNA damage in women with TNBCs. Neither age of onset, contralateral disease nor family history were able to discern all mutation positive individuals. Therefore, TNBC should be considered as an additional criterion for panel based genetic testing.

The polynucleotide kinase 3'-phosphatase gene (PNKP) is involved in Charcot-Marie-Tooth disease (CMT2B2) previously related to MED25

Charcot-Marie-Tooth disease (CMT) represents a heterogeneous group of hereditary peripheral neuropathies. We previously reported a CMT locus on chromosome 19q13.3 segregating with the disease in a large Costa Rican family with axonal neuropathy and autosomal recessive pattern of inheritance (CMT2B2). We proposed a homozygous missense variant in the Mediator complex 25 (MED25) gene as causative of the disease. Nevertheless, the fact that no other CMT individuals with MED25 variants were reported to date led us to reevaluate the original family. Using exome sequencing, we now identified a homozygous nonsense variant (p.Gln517ter) in the last exon of an adjacent gene, the polynucleotide kinase 3'-phosphatase (PNKP) gene. It encodes a DNA repair protein recently associated with recessive ataxia with oculomotor apraxia type 4 (AOA4) and microcephaly, seizures, and developmental delay (MCSZ). Subsequently, five unrelated Costa Rican CMT2 subjects initially identified as being heterozygous for the same MED25 variant were found to be also compound heterozygote for PNKP. All were heterozygous for the same variant found homozygous in the large family and a second one previously associated with ataxia (p.Thr408del). Detailed clinical reassessment of the initial family and the new individuals revealed in all an adult-onset slowly progressive CMT2 associated with signs of cerebellar dysfunction such as slurred speech and oculomotor involvement, but neither microcephaly, seizures, nor developmental delay. We propose that PKNP variants are the major causative variant for the CMT2 phenotype in these individuals and that the milder clinical manifestation is due to an allelic effect.

Systematic approach demonstrates enrichment of multiple interactions between non-HLA risk variants and HLA-DRB1 risk alleles in rheumatoid arthritis

OBJECTIVE

In anti-citrullinated protein antibody positive rheumatoid arthritis (ACPA-positive RA), a particular subset of HLA-DRB1 alleles, called shared epitope (SE) alleles, is a highly influential genetic risk factor. Here, we investigated whether non-HLA single nucleotide polymorphisms (SNP), conferring low disease risk on their own, interact with SE alleles more frequently than expected by chance and if such genetic interactions influence the HLA-DRB1 SE effect concerning risk to ACPA-positive RA.

METHODS

We computed the attributable proportion (AP) due to additive interaction at genome-wide level for two independent ACPA-positive RA cohorts: the Swedish epidemiological investigation of rheumatoid arthritis (EIRA) and the North American rheumatoid arthritis consortium (NARAC). Then, we tested for differences in the AP p value distributions observed for two groups of SNPs, non-associated and associated with disease. We also evaluated whether the SNPs in interaction with HLA-DRB1 were cis-eQTLs in the SE alleles context in peripheral blood mononuclear cells from patients with ACPA-positive RA (SE-eQTLs).

RESULTS

We found a strong enrichment of significant interactions (AP p<0.05) between the HLA-DRB1 SE alleles and the group of SNPs associated with ACPA-positive RA in both cohorts (Kolmogorov-Smirnov test D=0.35 for EIRA and D=0.25 for NARAC, p<2.2e-16 for both). Interestingly, 564 out of 1492 SNPs in consistent interaction for both cohorts were significant SE-eQTLs. Finally, we observed that the effect size of HLA-DRB1 SE alleles for disease decreases from 5.2 to 2.5 after removal of the risk alleles of the two top interacting SNPs (rs2476601 and rs10739581).

CONCLUSION

Our data demonstrate that there are massive genetic interactions between the HLA-DRB1 SE alleles and non-HLA genetic variants in ACPA-positive RA.

Integrative bioinformatics analysis characterizing the role of EDC3 in mRNA decay and its association to intellectual disability

BACKGROUND

Decapping of mRNA is an important step in the regulation of mRNA turnover and therefore of gene expression, which is a key process controlling development and homeostasis of all organisms. It has been shown that EDC3 plays a role in mRNA decapping, however its function is not well understood. Previously, we have associated a homozygous variant in EDC3 with autosomal recessive intellectual disability. Here, we investigate the functional role of EDC3.

METHODS

We performed transcriptome analyses in patients' samples. In addition, we established an EDC3 loss-of-function model using siRNA-based knockdown in the human neuroblastoma cell line SKNBE and carried out RNA sequencing. Integrative bioinformatics analyses were performed to identify EDC3-dependent candidate genes and/or pathways.

RESULTS

Our analyses revealed that 235 genes were differentially expressed in patients versus controls. In addition, AU-rich element (ARE)-containing mRNAs, whose degradation in humans has been suggested to involve EDC3, had higher fold changes than non-ARE-containing genes. The analysis of RNA sequencing data from the EDC3 in vitro loss-of-function model confirmed the higher fold changes of ARE-containing mRNAs compared to non-ARE-containing mRNAs and further showed an upregulation of long non-coding and coding RNAs. In total, 764 genes were differentially expressed. Integrative bioinformatics analyses of these genes identified dysregulated candidate pathways, including pathways related to synapses/coated vesicles and DNA replication/cell cycle.

CONCLUSION

Our data support the involvement of EDC3 in mRNA decay, including ARE-containing mRNAs, and suggest that EDC3 might be preferentially involved in the degradation of long coding and non-coding RNAs. Furthermore, our results associate ECD3 loss-of-function with synapses-related pathways. Collectively, our data provide novel information that might help elucidate the molecular mechanisms underlying the association of intellectual disability with the dysregulation of mRNA degradation.

T cells are influenced by a long non-coding RNA in the autoimmune associated PTPN2 locus

Non-coding SNPs in the protein tyrosine phosphatase non-receptor type 2 (PTPN2) locus have been linked with several autoimmune diseases, including rheumatoid arthritis, type I diabetes, and inflammatory bowel disease. However, the functional consequences of these SNPs are poorly characterized. Herein, we show in blood cells that SNPs in the PTPN2 locus are highly correlated with DNA methylation levels at four CpG sites downstream of PTPN2 and expression levels of the long non-coding RNA (lncRNA) LINC01882 downstream of these CpG sites. We observed that LINC01882 is mainly expressed in T cells and that anti-CD3/CD28 activated naïve CD4⁺ T cells downregulate the expression of LINC01882. RNA sequencing analysis of LINC01882 knockdown in Jurkat T cells, using a combination of antisense oligonucleotides and RNA interference, revealed the upregulation of the transcription factor ZEB1 and kinase MAP2K4, both involved in IL-2 regulation. Overall, our data suggests the involvement of LINC01882 in T cell activation and hints towards an auxiliary role of these non-coding SNPs in autoimmunity associated with the PTPN2 locus.

The genetic basis for most patients with pustular skin disease remains elusive

BACKGROUND

Rare variants in the genes IL36RN, CARD14 and AP1S3 have been identified to cause or contribute to pustular skin diseases, primarily generalized pustular psoriasis (GPP).

OBJECTIVES

To better understand the disease relevance of these genes, we screened our cohorts of patients with pustular skin diseases [primarily GPP and palmoplantar pustular psoriasis (PPP)] for coding changes in these three genes. Carriers of single heterozygous IL36RN mutations were screened for a second mutation in IL36RN.

METHODS

Coding exons of IL36RN, CARD14 and AP1S3 were sequenced in 67 patients - 61 with GPP, two with acute generalized exanthematous pustulosis and four with acrodermatitis continua of Hallopeau. We screened IL36RN and AP1S3 for intragenic copy-number variants and 258 patients with PPP for coding changes in AP1S3. Eleven heterozygous IL36RN mutations carriers were analysed for a second noncoding IL36RN mutation. Genotype-phenotype correlations in carriers/noncarriers of IL36RN mutations were assessed within the GPP cohort.

RESULTS

The majority of patients (GPP, 64%) did not carry rare variants in any of the three genes. Biallelic and monoallelic IL36RN mutations were identified in 15 and five patients with GPP, respectively. Noncoding rare IL36RN variants were not identified in heterozygous carriers. The only significant genotype-phenotype correlation observed for IL36RN mutation carriers was early age at disease onset. Additional rare CARD14 or AP1S3 variants were identified in 15% of IL36RN mutation carriers.

CONCLUSIONS

The identification of IL36RN mutation carriers harbouring additional rare variants in CARD14 or AP1S3 indicates a more complex mode of inheritance of pustular psoriasis. Our results suggest that, in heterozygous IL36RN mutation carriers, there are additional disease-causing genetic factors outside IL36RN.

High resolution chromosomal microarray analysis in paediatric obsessive-compulsive disorder

Background

Obsessive-Compulsive Disorder (OCD) is a common and chronic disorder in which a person has uncontrollable, reoccurring thoughts and behaviours. It is a complex genetic condition and, in case of early onset (EO), the patients manifest a more severe phenotype, and an increased heritability. Large (>500 kb) copy number variations (CNVs) previously associated with autism and schizophrenia have been reported in OCD. Recently, rare CNVs smaller than 500 kb overlapping risk loci for other neurodevelopmental conditions have also been reported in OCD, stressing the importance of examining CNVs of any size range. The aim of this study was to further investigate the role of rare and small CNVs in the aetiology of EO-OCD.

Methods

We performed high-resolution chromosomal microarray analysis in 121 paediatric OCD patients and in 124 random controls to identify rare CNVs (>50 kb) which might contribute to EO-OCD.

Results

The frequencies and the size of the observed rare CNVs in the patients did not differ from the controls. However, we observed a significantly higher frequency of rare CNVs affecting brain related genes, especially deletions, in the patients (OR = 1.98, 95% CI 1.02–3.84; OR = 3.61, 95% CI 1.14–11.41, respectively). Similarly, enrichment-analysis of CNVs gene content, performed with three independent methods, confirmed significant clustering of predefined genes involved in synaptic/brain related functional pathways in the patients but not in the controls. In two patients we detected de-novo CNVs encompassing genes previously associated with different neurodevelopmental disorders (NRXN1, ANKS1B, UHRF1BP1).

Conclusions

Our results further strengthen the role of small rare CNVs, particularly deletions, as susceptibility factors for paediatric OCD.

Electronic supplementary material

The online version of this article (10.1186/s12920-017-0299-5) contains supplementary material, which is available to authorized users.

Genetic screening confirms heterozygous mutations in ACAN as a major cause of idiopathic short stature

Short stature is a common pediatric disorder affecting 3% of the population. However, the clinical variability and genetic heterogeneity prevents the identification of the underlying cause in about 80% of the patients. Recently, heterozygous mutations in the ACAN gene coding for the proteoglycan aggrecan, a main component of the cartilage matrix, were associated with idiopathic short stature. To ascertain the prevalence of ACAN mutations and broaden the phenotypic spectrum in patients with idiopathic short stature we performed sequence analyses in 428 families. We identified heterozygous nonsense mutations in four and potentially disease-causing missense variants in two families (1.4%). These patients presented with a mean of −3.2 SDS and some suggestive clinical characteristics. The results suggest heterozygous mutations in ACAN as a common cause of isolated as well as inherited idiopathic short stature.

Genome-wide association and targeted analysis of copy number variants with psoriatic arthritis in German patients

Background

Psoriatic Arthritis (PsA) is a chronic inflammatory disease of the joints. PsA is etiologically complex, and 11 susceptibility loci have been identified so far. Most of these overlap with loci associated with psoriasis vulgaris (PsV), the most common psoriatic skin manifestation which is also frequently seen in PsA patients. In addition, two copy number variants (CNVs) are associated with PsV, one of which, located within the LCE3 gene cluster, is also associated with PsA. Finally, an intergenic deletion has been reported as a PsA-specific CNV.

Methods

We performed a genome-wide association study (GWAS) of CNVs in PsA and assessed the contribution to disease risk by CNVs at known psoriasis susceptibility loci.

Results

After stringent quality assessment and validation of CNVs of the GWAS with an alternative quantitative method, two significantly associated CNVs remained, one near UXS1, the other one at the TRB locus. However, MLPA analysis did not confirm the CN state in ~1/3 of individuals, and an analysis of an independent case-control-study failed to confirm the initial associations. Furthermore, detailed PCR-based analysis of the sequence at TRB revealed the existence of a more complex genomic sequence most accurately represented by freeze hg18 which accordingly failed to confirm the hg19 sequence.

Only rare CNVs were detected at psoriasis susceptibility loci. At three of 12 susceptibility loci with CNVs (CSMD1, IL12B, RYR2), CN variability was confirmed independently by MLPA. Overall, the rate of CNV confirmation by MLPA was strongly dependent upon CNV type, CNV size and the number of array markers involved in a CNV.

Conclusion

Although we identified PsA associations at several loci and confirmed that the common CNVs at these sites were real, ~1/3 of the common CNV states could not be reproduced. Furthermore, replication analysis failed to confirm the original association. Furthermore, SNP array-based analyses of CNVs were found to be more reliable for deletions than duplications, independent of the respective CNV allele frequency. CNVs are thus good candidate disease variants, while the methods to detect them should be applied cautiously and reproduced by an independent method.

Electronic supplementary material

The online version of this article (doi:10.1186/s12881-017-0447-y) contains supplementary material, which is available to authorized users.

Fra-2 regulates B cell development by enhancing IRF4 and Foxo1 transcription

The role of AP-1 transcription factors in early B cell development and function is still incompletely characterized. Ubieta at al. describe the function of the Fra-2/AP-1 transcription factor as a regulator of Foxo1 and Irf4 expression in B cells. Fra-2 affects B cell proliferation and maintains their number in bone marrow.

The role of AP-1 transcription factors in early B cell development and function is still incompletely characterized. Here we address the role of Fra-2 in B cell differentiation. Deletion of Fra-2 leads to impaired B cell proliferation in the bone marrow. In addition, IL-7–stimulated pro–B cell cultures revealed a reduced differentiation from large pre–B cells to small B cells and immature B cells. Gene profiling and chromatin immunoprecipitation sequencing analyses unraveled a transcriptional reduction of the transcription factors Foxo1, Irf4, Ikaros, and Aiolos in Fra-2–deficient B cells. Moreover, expression of IL7Rα and Rag 1/2, downstream targets of Irf4 and Foxo1, were also reduced in the absence of Fra-2. Pro–B cell proliferation and small pre–B cell differentiation were fully rescued by expression of Foxo1 and Irf4 in Fra-2–deficient pro–B cells. Hence, Fra-2 is a key upstream regulator of Foxo1 and Irf4 expression and influences proliferation and differentiation of B cells at multiple stages.

Genetic association study of exfoliation syndrome identifies a protective rare variant at LOXL1 and five new susceptibility loci

Exfoliation syndrome (XFS) is the most common known risk factor for secondary glaucoma and a major cause of blindness worldwide. Variants in two genes, LOXL1 and CACNA1A, have previously been associated with XFS. To further elucidate the genetic basis of XFS, we collected a global sample of XFS cases to refine the association at LOXL1, which previously showed inconsistent results across populations, and to identify new variants associated with XFS. We identified a rare protective allele at LOXL1 (p.Phe407, odds ratio (OR) = 25, P = 2.9 × 10^-14) through deep resequencing of XFS cases and controls from nine countries. A genome-wide association study (GWAS) of XFS cases and controls from 24 countries followed by replication in 18 countries identified seven genome-wide significant loci (P < 5 × 10^-8). We identified association signals at 13q12 (POMP), 11q23.3 (TMEM136), 6p21 (AGPAT1), 3p24 (RBMS3) and 5q23 (near SEMA6A). These findings provide biological insights into the pathology of XFS and highlight a potential role for naturally occurring rare LOXL1 variants in disease biology.

Pseudoexfoliation syndrome-associated genetic variants affect transcription factor binding and alternative splicing of LOXL1

Although lysyl oxidase-like 1 (LOXL1) is known as the principal genetic risk factor for pseudoexfoliation (PEX) syndrome, a major cause of glaucoma and cardiovascular complications, no functional variants have been identified to date. Here, we conduct a genome-wide association scan on 771 German PEX patients and 1,350 controls, followed by independent testing of associated variants in Italian and Japanese data sets. We focus on a 3.5-kb four-component polymorphic locus positioned spanning introns 1 and 2 of LOXL1 with enhancer-like chromatin features. We find that the rs11638944:C>G transversion exerts a cis-acting effect on the expression levels of LOXL1, mediated by differential binding of the transcription factor RXRα (retinoid X receptor alpha) and by modulating alternative splicing of LOXL1, eventually leading to reduced levels of LOXL1 mRNA in cells and tissues of risk allele carriers. These findings uncover a functional mechanism by which common noncoding variants influence LOXL1 expression.

LOXL1 is a genetic risk factor for pseudoexfoliation syndrome of the eye but a causal variant has not been identified. Here, Pasutto et al., find intronic LOXL1 risk variants influence transcription factor binding and alternative splicing of LOXL1 in affected tissues reducing levels of LOXL1 mRNA.

Association Between Telomere Length and Risk of Cancer and Non-Neoplastic Diseases: A Mendelian Randomization Study

IMPORTANCE

The causal direction and magnitude of the association between telomere length and incidence of cancer and non-neoplastic diseases is uncertain owing to the susceptibility of observational studies to confounding and reverse causation.

OBJECTIVE

To conduct a Mendelian randomization study, using germline genetic variants as instrumental variables, to appraise the causal relevance of telomere length for risk of cancer and non-neoplastic diseases.

DATA SOURCES

Genomewide association studies (GWAS) published up to January 15, 2015.

STUDY SELECTION

GWAS of noncommunicable diseases that assayed germline genetic variation and did not select cohort or control participants on the basis of preexisting diseases. Of 163 GWAS of noncommunicable diseases identified, summary data from 103 were available.

DATA EXTRACTION AND SYNTHESIS

Summary association statistics for single nucleotide polymorphisms (SNPs) that are strongly associated with telomere length in the general population.

MAIN OUTCOMES AND MEASURES

Odds ratios (ORs) and 95% confidence intervals (CIs) for disease per standard deviation (SD) higher telomere length due to germline genetic variation.

RESULTS

Summary data were available for 35 cancers and 48 non-neoplastic diseases, corresponding to 420 081 cases (median cases, 2526 per disease) and 1 093 105 controls (median, 6789 per disease). Increased telomere length due to germline genetic variation was generally associated with increased risk for site-specific cancers. The strongest associations (ORs [95% CIs] per 1-SD change in genetically increased telomere length) were observed for glioma, 5.27 (3.15-8.81); serous low-malignant-potential ovarian cancer, 4.35 (2.39-7.94); lung adenocarcinoma, 3.19 (2.40-4.22); neuroblastoma, 2.98 (1.92-4.62); bladder cancer, 2.19 (1.32-3.66); melanoma, 1.87 (1.55-2.26); testicular cancer, 1.76 (1.02-3.04); kidney cancer, 1.55 (1.08-2.23); and endometrial cancer, 1.31 (1.07-1.61). Associations were stronger for rarer cancers and at tissue sites with lower rates of stem cell division. There was generally little evidence of association between genetically increased telomere length and risk of psychiatric, autoimmune, inflammatory, diabetic, and other non-neoplastic diseases, except for coronary heart disease (OR, 0.78 [95% CI, 0.67-0.90]), abdominal aortic aneurysm (OR, 0.63 [95% CI, 0.49-0.81]), celiac disease (OR, 0.42 [95% CI, 0.28-0.61]) and interstitial lung disease (OR, 0.09 [95% CI, 0.05-0.15]).

CONCLUSIONS AND RELEVANCE

It is likely that longer telomeres increase risk for several cancers but reduce risk for some non-neoplastic diseases, including cardiovascular diseases.

DNA methylation mediates genotype and smoking interaction in the development of anti-citrullinated peptide antibody-positive rheumatoid arthritis

BACKGROUND

Multiple factors, including interactions between genetic and environmental risks, are important in susceptibility to rheumatoid arthritis (RA). However, the underlying mechanism is not fully understood. This study was undertaken to evaluate whether DNA methylation can mediate the interaction between genotype and smoking in the development of anti-citrullinated peptide antibody (ACPA)-positive RA.

METHODS

We investigated the gene-smoking interactions in DNA methylation using 393 individuals from the Epidemiological Investigation of Rheumatoid Arthritis (EIRA). The interaction between rs6933349 and smoking in the risk of developing ACPA-positive RA was further evaluated in a larger portion of the EIRA (1119 controls and 944 ACPA-positive patients with RA), and in the Malaysian Epidemiological Investigation of Rheumatoid Arthritis (MyEIRA) (1556 controls and 792 ACPA-positive patients with RA). Finally, mediation analysis was performed to investigate whether DNA methylation of cg21325723 mediates this gene-environment interaction on the risk of developing of ACPA-positive RA.

RESULTS

We identified and replicated one significant gene-environment interaction between rs6933349 and smoking in DNA methylation of cg21325723. This gene-smoking interaction is a novel interaction in the risk of developing ACPA-positive in both Caucasian (multiplicative P value = 0.056; additive P value = 0.016) and Asian populations (multiplicative P value = 0.035; additive P value = 0.00027), and it is mediated through DNA methylation of cg21325723.

CONCLUSIONS

We showed that DNA methylation of cg21325723 can mediate the gene-environment interaction between rs6933349 and smoking, impacting the risk of developing ACPA-positive RA, thus being a potential regulator that integrates both internal genetic and external environmental risk factors.

Hypomorphic Pathogenic Variants in TAF13 Are Associated with Autosomal-Recessive Intellectual Disability and Microcephaly

In two independent consanguineous families each with two children affected by mild intellectual disability and microcephaly, we identified two homozygous missense variants (c.119T>A [p.Met40Lys] and c.92T>A [p.Leu31His]) in TATA-box-binding-protein-associated factor 13 (TAF13). Molecular modeling suggested a pathogenic effect of both variants through disruption of the interaction between TAF13 and TAF11. These two proteins form a histone-like heterodimer that is essential for their recruitment into the general RNA polymerase II transcription factor IID (TFIID) complex. Co-immunoprecipitation in HeLa cells transfected with plasmids encoding TAF11 and TAF13 revealed that both variants indeed impaired formation of the TAF13-TAF11 heterodimer, thus confirming the protein modeling analysis. To further understand the functional role of TAF13, we performed RNA sequencing of neuroblastoma cell lines upon TAF13 knockdown. The transcriptional profile showed significant deregulation of gene expression patterns with an emphasis on genes related to neuronal and skeletal functions and those containing E-box motives in their promoters. Here, we expand the spectrum of TAF-associated phenotypes and highlight the importance of TAF13 in neuronal functions.

Diagnostic Yield and Novel Candidate Genes by Exome Sequencing in 152 Consanguineous Families With Neurodevelopmental Disorders

IMPORTANCE

Autosomal recessive inherited neurodevelopmental disorders are highly heterogeneous, and many, possibly most, of the disease genes are still unknown.

OBJECTIVES

To promote the identification of disease genes through confirmation of previously described genes and presentation of novel candidates and provide an overview of the diagnostic yield of exome sequencing in consanguineous families.

DESIGN, SETTING, AND PARTICIPANTS

Autozygosity mapping in families and exome sequencing of index patients were performed in 152 consanguineous families (the parents descended from a same ancestor) with at least 1 offspring with intellectual disability (ID). The study was conducted from July 1, 2008, to June 30, 2015, and data analysis was conducted from July 1, 2015, to August 31, 2016.

RESULTS

Of the 152 consanguineous families enrolled, 1 child (in 45 families [29.6%]) or multiple children (107 families [70.4%]) had ID; additional features were present in 140 of the families (92.1%). The mean (SD) age of the children was 10.3 (9.0) years, and 171 of 297 (57.6%) were male. In 109 families (71.7%), potentially protein-disrupting and clinically relevant variants were identified. Of these, a clear clinical genetic diagnosis was made in 56 families (36.8%) owing to 57 (likely) pathogenic variants in 50 genes already established in neurodevelopmental disorders (46 autosomal recessive, 2 X-linked, and 2 de novo) or in 7 previously proposed recessive candidates. In 5 of these families, potentially treatable disorders were diagnosed (mutations in PAH, CBS, MTHFR, CYP27A1, and HIBCH), and in 1 family, 2 disease-causing homozygous variants in different genes were identified. In another 48 families (31.6%), 52 convincing recessive variants in candidate genes that were not previously reported in regard to neurodevelopmental disorders were identified. Of these, 14 were homozygous and truncating in GRM7, STX1A, CCAR2, EEF1D, GALNT2, SLC44A1, LRRIQ3, AMZ2, CLMN, SEC23IP, INIP, NARG2, FAM234B, and TRAP1. The diagnostic yield was higher in individuals with severe ID (35 of 77 [45.5%]), in multiplex families (42 of 107 [39.3%]), in patients with additional features (30 of 70 [42.9%]), and in those with remotely related parents (15 of 34 [44.1%]).

CONCLUSIONS AND RELEVANCE

Because of the high diagnostic yield of 36.8% and the possibility of identifying treatable diseases or the coexistence of several disease-causing variants, using exome sequencing as a first-line diagnostic approach in consanguineous families with neurodevelopmental disorders is recommended. Furthermore, the literature is enriched with 52 convincing candidate genes that are awaiting confirmation in independent families.