PSMD11 loss-of-function variants correlate with a neurobehavioral phenotype, obesity, and increased interferon response.

Primary proteasomopathies have recently emerged as a new class of rare early-onset neurodevelopmental disorders (NDDs) caused by pathogenic variants in the PSMB1, PSMC1, PSMC3, or PSMD12 proteasome genes. Proteasomes are large multi-subunit protein complexes that maintain cellular protein homeostasis by clearing ubiquitin-tagged damaged, misfolded, or unnecessary proteins. In this study, we have identified PSMD11 as an additional proteasome gene in which pathogenic variation is associated with an NDD-causing proteasomopathy. PSMD11 loss-of-function variants caused early-onset syndromic intellectual disability and neurodevelopmental delay with recurrent obesity in 10 unrelated children. Our findings demonstrate that the cognitive impairment observed in these individuals could be recapitulated in Drosophila melanogaster with depletion of the PMSD11 ortholog Rpn6, which compromised reversal learning. Our investigations in subject samples further revealed that PSMD11 loss of function resulted in impaired 26S proteasome assembly and the acquisition of a persistent type I interferon (IFN) gene signature, mediated by the integrated stress response (ISR) protein kinase R (PKR). In summary, these data identify PSMD11 as an additional member of the growing family of genes associated with neurodevelopmental proteasomopathies and provide insights into proteasomal biology in human health.

Role of CAMK2D in neurodevelopment and associated conditions.

The calcium/calmodulin-dependent protein kinase type 2 (CAMK2) family consists of four different isozymes, encoded by four different genes-CAMK2A, CAMK2B, CAMK2G, and CAMK2D-of which the first three have been associated recently with neurodevelopmental disorders. CAMK2D is one of the major CAMK2 proteins expressed in the heart and has been associated with cardiac anomalies. Although this CAMK2 isoform is also known to be one of the major CAMK2 subtypes expressed during early brain development, it has never been linked with neurodevelopmental disorders until now. Here we show that CAMK2D plays an important role in neurodevelopment not only in mice but also in humans. We identified eight individuals harboring heterozygous variants in CAMK2D who display symptoms of intellectual disability, delayed speech, behavioral problems, and dilated cardiomyopathy. The majority of the variants tested lead to a gain of function (GoF), which appears to cause both neurological problems and dilated cardiomyopathy. In contrast, loss-of-function (LoF) variants appear to induce only neurological symptoms. Together, we describe a cohort of individuals with neurodevelopmental disorders and cardiac anomalies, harboring pathogenic variants in CAMK2D, confirming an important role for the CAMK2D isozyme in both heart and brain function.

Bi-allelic variants in DOHH, catalyzing the last step of hypusine biosynthesis, are associated with a neurodevelopmental disorder.

Deoxyhypusine hydroxylase (DOHH) is the enzyme catalyzing the second step in the post-translational synthesis of hypusine [Nε-(4-amino-2-hydroxybutyl)lysine] in the eukaryotic initiation factor 5A (eIF5A). Hypusine is formed exclusively in eIF5A by two sequential enzymatic steps catalyzed by deoxyhypusine synthase (DHPS) and deoxyhypusine hydroxylase (DOHH). Hypusinated eIF5A is essential for translation and cell proliferation in eukaryotes, and all three genes encoding eIF5A, DHPS, and DOHH are highly conserved throughout eukaryotes. Pathogenic variants affecting either DHPS or EIF5A have been previously associated with neurodevelopmental disorders. Using trio exome sequencing, we identified rare bi-allelic pathogenic missense and truncating DOHH variants segregating with disease in five affected individuals from four unrelated families. The DOHH variants are associated with a neurodevelopmental phenotype that is similar to phenotypes caused by DHPS or EIF5A variants and includes global developmental delay, intellectual disability, facial dysmorphism, and microcephaly. A two-dimensional gel analyses revealed the accumulation of deoxyhypusine-containing eIF5A [eIF5A(Dhp)] and a reduction in the hypusinated eIF5A in fibroblasts derived from affected individuals, providing biochemical evidence for deficiency of DOHH activity in cells carrying the bi-allelic DOHH variants. Our data suggest that rare bi-allelic variants in DOHH result in reduced enzyme activity, limit the hypusination of eIF5A, and thereby lead to a neurodevelopmental disorder.

THUMPD1 bi-allelic variants cause loss of tRNA acetylation and a syndromic neurodevelopmental disorder.

Covalent tRNA modifications play multi-faceted roles in tRNA stability, folding, and recognition, as well as the rate and fidelity of translation, and other cellular processes such as growth, development, and stress responses. Mutations in genes that are known to regulate tRNA modifications lead to a wide array of phenotypes and diseases including numerous cognitive and neurodevelopmental disorders, highlighting the critical role of tRNA modification in human disease. One such gene, THUMPD1, is involved in regulating tRNA N4-acetylcytidine modification (ac4C), and recently was proposed as a candidate gene for autosomal-recessive intellectual disability. Here, we present 13 individuals from 8 families who harbor rare loss-of-function variants in THUMPD1. Common phenotypic findings included global developmental delay, speech delay, moderate to severe intellectual deficiency, behavioral abnormalities such as angry outbursts, facial dysmorphism, and ophthalmological abnormalities. We demonstrate that the bi-allelic variants identified cause loss of function of THUMPD1 and that this defect results in a loss of ac4C modification in small RNAs, and of individually purified tRNA-Ser-CGA. We further corroborate this effect by showing a loss of tRNA acetylation in two CRISPR-Cas9-generated THUMPD1 KO cell lines. In addition, we also show the resultant amino acid substitution that occurs in a missense THUMPD1 allele identified in an individual with compound heterozygous variants results in a marked decrease in THUMPD1 stability and RNA-binding capacity. Taken together, these results suggest that the lack of tRNA acetylation due to THUMPD1 loss of function results in a syndromic form of intellectual disability associated with developmental delay, behavioral abnormalities, hearing loss, and facial dysmorphism.

Circulating white blood cell traits and colorectal cancer risk: A Mendelian randomisation study.

Observational studies have suggested a protective role for eosinophils in colorectal cancer (CRC) development and implicated neutrophils, but the causal relationships remain unclear. Here, we aimed to estimate the causal effect of circulating white blood cell (WBC) counts (N = ~550 000) for basophils, eosinophils, monocytes, lymphocytes and neutrophils on CRC risk (N = 52 775 cases and 45 940 controls) using Mendelian randomisation (MR). For comparison, we also examined this relationship using individual-level data from UK Biobank (4043 incident CRC cases and 332 773 controls) in a longitudinal cohort analysis. The inverse-variance weighted (IVW) MR analysis suggested a protective effect of increased basophil count and eosinophil count on CRC risk [OR per 1-SD increase: 0.88, 95% CI: 0.78-0.99, P = .04; OR: 0.93, 95% CI: 0.88-0.98, P = .01]. The protective effect of eosinophils remained [OR per 1-SD increase: 0.88, 95% CI: 0.80-0.97, P = .01] following adjustments for all other WBC subtypes, to account for genetic correlation between the traits, using multivariable MR. A protective effect of increased lymphocyte count on CRC risk was also found [OR: 0.84, 95% CI: 0.76-0.93, P = 6.70e-4] following adjustment. Consistent with MR results, a protective effect for eosinophils in the cohort analysis in the fully adjusted model [RR per 1-SD increase: 0.96, 95% CI: 0.93-0.99, P = .02] and following adjustment for the other WBC subtypes [RR: 0.96, 95% CI: 0.93-0.99, P = .001] was observed. Our study implicates peripheral blood immune cells, in particular eosinophils and lymphocytes, in CRC development, highlighting a need for mechanistic studies to interrogate these relationships.

SEMA6B variants cause intellectual disability and alter dendritic spine density and axon guidance.

Intellectual disability (ID) is a neurodevelopmental disorder frequently caused by monogenic defects. In this study, we collected 14 SEMA6B heterozygous variants in 16 unrelated patients referred for ID to different centers. Whereas, until now, SEMA6B variants have mainly been reported in patients with progressive myoclonic epilepsy, our study indicates that the clinical spectrum is wider and also includes non-syndromic ID without epilepsy or myoclonus. To assess the pathogenicity of these variants, selected mutated forms of Sema6b were overexpressed in Human Embryonic Kidney 293T (HEK293T) cells and in primary neuronal cultures. shRNAs targeting Sema6b were also used in neuronal cultures to measure the impact of the decreased Sema6b expression on morphogenesis and synaptogenesis. The overexpression of some variants leads to a subcellular mislocalization of SEMA6B protein in HEK293T cells and to a reduced spine density owing to loss of mature spines in neuronal cultures. Sema6b knockdown also impairs spine density and spine maturation. In addition, we conducted in vivo rescue experiments in chicken embryos with the selected mutated forms of Sema6b expressed in commissural neurons after knockdown of endogenous SEMA6B. We observed that expression of these variants in commissural neurons fails to rescue the normal axon pathway. In conclusion, identification of SEMA6B variants in patients presenting with an overlapping phenotype with ID and functional studies highlight the important role of SEMA6B in neuronal development, notably in spine formation and maturation and in axon guidance. This study adds SEMA6B to the list of ID-related genes.

Missense variants in DPYSL5 cause a neurodevelopmental disorder with corpus callosum agenesis and cerebellar abnormalities.

The collapsin response mediator protein (CRMP) family proteins are intracellular mediators of neurotrophic factors regulating neurite structure/spine formation and are essential for dendrite patterning and directional axonal pathfinding during brain developmental processes. Among this family, CRMP5/DPYSL5 plays a significant role in neuronal migration, axonal guidance, dendrite outgrowth, and synapse formation by interacting with microtubules. Here, we report the identification of missense mutations in DPYSL5 in nine individuals with brain malformations, including corpus callosum agenesis and/or posterior fossa abnormalities, associated with variable degrees of intellectual disability. A recurrent de novo p.Glu41Lys variant was found in eight unrelated patients, and a p.Gly47Arg variant was identified in one individual from the first family reported with Ritscher-Schinzel syndrome. Functional analyses of the two missense mutations revealed impaired dendritic outgrowth processes in young developing hippocampal primary neuronal cultures. We further demonstrated that these mutations, both located in the same loop on the surface of DPYSL5 monomers and oligomers, reduced the interaction of DPYSL5 with neuronal cytoskeleton-associated proteins MAP2 and ßIII-tubulin. Our findings collectively indicate that the p.Glu41Lys and p.Gly47Arg variants impair DPYSL5 function on dendritic outgrowth regulation by preventing the formation of the ternary complex with MAP2 and ßIII-tubulin, ultimately leading to abnormal brain development. This study adds DPYSL5 to the list of genes implicated in brain malformation and in neurodevelopmental disorders.

Classification of 101 BRCA1 and BRCA2 variants of uncertain significance by cosegregation study: A powerful approach.

Up to 80% of BRCA1 and BRCA2 genetic variants remain of uncertain clinical significance (VUSs). Only variants classified as pathogenic or likely pathogenic can guide breast and ovarian cancer prevention measures and treatment by PARP inhibitors. We report the first results of the ongoing French national COVAR (cosegregation variant) study, the aim of which is to classify BRCA1/2 VUSs. The classification method was a multifactorial model combining different associations between VUSs and cancer, including cosegregation data. At this time, among the 653 variants selected, 101 (15%) distinct variants shared by 1,624 families were classified as pathogenic/likely pathogenic or benign/likely benign by the COVAR study. Sixty-six of the 101 (65%) variants classified by COVAR would have remained VUSs without cosegregation data. Of note, among the 34 variants classified as pathogenic by COVAR, 16 remained VUSs or likely pathogenic when following the ACMG/AMP variant classification guidelines. Although the initiation and organization of cosegregation analyses require a considerable effort, the growing number of available genetic tests results in an increasing number of families sharing a particular variant, and thereby increases the power of such analyses. Here we demonstrate that variant cosegregation analyses are a powerful tool for the classification of variants in the BRCA1/2 breast-ovarian cancer predisposition genes.

Haploinsufficiency of the Sin3/HDAC corepressor complex member SIN3B causes a syndromic intellectual disability/autism spectrum disorder.

Proteins involved in transcriptional regulation harbor a demonstrated enrichment of mutations in neurodevelopmental disorders. The Sin3 (Swi-independent 3)/histone deacetylase (HDAC) complex plays a central role in histone deacetylation and transcriptional repression. Among the two vertebrate paralogs encoding the Sin3 complex, SIN3A variants cause syndromic intellectual disability, but the clinical consequences of SIN3B haploinsufficiency in humans are uncharacterized. Here, we describe a syndrome hallmarked by intellectual disability, developmental delay, and dysmorphic facial features with variably penetrant autism spectrum disorder, congenital malformations, corpus callosum defects, and impaired growth caused by disruptive SIN3B variants. Using chromosomal microarray or exome sequencing, and through international data sharing efforts, we identified nine individuals with heterozygous SIN3B deletion or single-nucleotide variants. Five individuals harbor heterozygous deletions encompassing SIN3B that reside within a ∼230 kb minimal region of overlap on 19p13.11, two individuals have a rare nonsynonymous substitution, and two individuals have a single-nucleotide deletion that results in a frameshift and predicted premature termination codon. To test the relevance of SIN3B impairment to measurable aspects of the human phenotype, we disrupted the orthologous zebrafish locus by genome editing and transient suppression. The mutant and morphant larvae display altered craniofacial patterning, commissural axon defects, and reduced body length supportive of an essential role for Sin3 function in growth and patterning of anterior structures. To investigate further the molecular consequences of SIN3B variants, we quantified genome-wide enhancer and promoter activity states by using H3K27ac ChIP-seq. We show that, similar to SIN3A mutations, SIN3B disruption causes hyperacetylation of a subset of enhancers and promoters in peripheral blood mononuclear cells. Together, these data demonstrate that SIN3B haploinsufficiency leads to a hitherto unknown intellectual disability/autism syndrome, uncover a crucial role of SIN3B in the central nervous system, and define the epigenetic landscape associated with Sin3 complex impairment.

Biallelic USP14 variants cause a syndromic neurodevelopmental disorder.

PURPOSE: Imbalances in protein homeostasis affect human brain development, with the ubiquitin-proteasome system (UPS) and autophagy playing crucial roles in neurodevelopmental disorders (NDD). This study explores the impact of biallelic USP14 variants on neurodevelopment, focusing on its role as a key hub connecting UPS and autophagy. METHODS: Here, we identified biallelic USP14 variants in 4 individuals from 3 unrelated families: 1 fetus, a newborn with a syndromic NDD and 2 siblings affected by a progressive neurological disease. Specifically, the 2 siblings from the latter family carried 2 compound heterozygous variants c.8T>C p.(Leu3Pro) and c.988C>T p.(Arg330∗), whereas the fetus had a homozygous frameshift c.899_902del p.(Lys300Serfs∗24) variant, and the newborn patient harbored a homozygous frameshift c.233_236del p.(Leu78Glnfs∗11) variant. Functional studies were conducted using sodium dodecyl-sulfate polyacrylamide gel electrophoresis, western blotting, and mass spectrometry analyses in both patient-derived and CRISPR-Cas9-generated cells. RESULTS: Our investigations indicated that the USP14 variants correlated with reduced N-terminal methionine excision, along with profound alterations in proteasome, autophagy, and mitophagy activities. CONCLUSION: Biallelic USP14 variants in NDD patients perturbed protein degradation pathways, potentially contributing to disorder etiology. Altered UPS, autophagy, and mitophagy activities underscore the intricate interplay, elucidating their significance in maintaining proper protein homeostasis during brain development.

Biallelic loss-of-function variants in CACHD1 cause a novel neurodevelopmental syndrome with facial dysmorphism and multisystem congenital abnormalities.

PURPOSE: We established the genetic etiology of a syndromic neurodevelopmental condition characterized by variable cognitive impairment, recognizable facial dysmorphism, and a constellation of extra-neurological manifestations. METHODS: We performed phenotypic characterization of 6 participants from 4 unrelated families presenting with a neurodevelopmental syndrome and used exome sequencing to investigate the underlying genetic cause. To probe relevance to the neurodevelopmental phenotype and craniofacial dysmorphism, we established two- and three-dimensional human stem cell-derived neural models and generated a stable cachd1 zebrafish mutant on a transgenic cartilage reporter line. RESULTS: Affected individuals showed mild cognitive impairment, dysmorphism featuring oculo-auriculo abnormalities, and developmental defects involving genitourinary and digestive tracts. Exome sequencing revealed biallelic putative loss-of-function variants in CACHD1 segregating with disease in all pedigrees. RNA sequencing in CACHD1-depleted neural progenitors revealed abnormal expression of genes with key roles in Wnt signaling, neurodevelopment, and organ morphogenesis. CACHD1 depletion in neural progenitors resulted in reduced percentages of post-mitotic neurons and enlargement of 3D neurospheres. Homozygous cachd1 mutant larvae showed mandibular patterning defects mimicking human facial dysmorphism. CONCLUSION: Our findings support the role of loss-of-function variants in CACHD1 as the cause of a rare neurodevelopmental syndrome with facial dysmorphism and multisystem abnormalities.

Growth charts in DYRK1A syndrome.

DYRK1A Syndrome (OMIM #614104) is caused by pathogenic variations in the DYRK1A gene located on 21q22. Haploinsufficiency of DYRK1A causes a syndrome with global psychomotor delay and intellectual disability. Low birth weight, growth restriction with feeding difficulties, stature insufficiency, and microcephaly are frequently reported. This study aims to create specific growth charts for individuals with DYRK1A Syndrome and identify parameters for size prognosis. Growth parameters were obtained for 92 individuals with DYRK1A Syndrome (49 males vs. 43 females). The data were obtained from pediatric records, parent reporting, and scientific literature. Growth charts for height, weight, body mass index (BMI), and occipitofrontal circumference (OFC) were generated using generalized additive models through R package gamlss. The growth curves include height, weight, and OFC measurements for patients aged 0-5 years. In accordance with the literature, the charts show that individuals are more likely to present intrauterine growth restriction with low birth weight and microcephaly. The growth is then characterized by severe microcephaly, low weight, and short stature. This study proposes growth charts for widespread use in the management of patients with DYRK1A syndrome.

Integrating RNA-Seq into genome sequencing workflow enhances the analysis of structural variants causing neurodevelopmental disorders.

BACKGROUND: Molecular diagnosis of neurodevelopmental disorders (NDDs) is mainly based on exome sequencing (ES), with a diagnostic yield of 31% for isolated and 53% for syndromic NDD. As sequencing costs decrease, genome sequencing (GS) is gradually replacing ES for genome-wide molecular testing. As many variants detected by GS only are in deep intronic or non-coding regions, the interpretation of their impact may be difficult. Here, we showed that integrating RNA-Seq into the GS workflow can enhance the analysis of the molecular causes of NDD, especially structural variants (SVs), by providing valuable complementary information such as aberrant splicing, aberrant expression and monoallelic expression. METHODS: We performed trio-GS on a cohort of 33 individuals with NDD for whom ES was inconclusive. RNA-Seq on skin fibroblasts was then performed in nine individuals for whom GS was inconclusive and optical genome mapping (OGM) was performed in two individuals with an SV of unknown significance. RESULTS: We identified pathogenic or likely pathogenic variants in 16 individuals (48%) and six variants of uncertain significance. RNA-Seq contributed to the interpretation in three individuals, and OGM helped to characterise two SVs. CONCLUSION: Our study confirmed that GS significantly improves the diagnostic performance of NDDs. However, most variants detectable by GS alone are structural or located in non-coding regions, which can pose challenges for interpretation. Integration of RNA-Seq data overcame this limitation by confirming the impact of variants at the transcriptional or regulatory level. This result paves the way for new routinely applicable diagnostic protocols.

Probing the diabetes and colorectal cancer relationship using gene - environment interaction analyses.

BACKGROUND: Diabetes is an established risk factor for colorectal cancer. However, the mechanisms underlying this relationship still require investigation and it is not known if the association is modified by genetic variants. To address these questions, we undertook a genome-wide gene-environment interaction analysis. METHODS: We used data from 3 genetic consortia (CCFR, CORECT, GECCO; 31,318 colorectal cancer cases/41,499 controls) and undertook genome-wide gene-environment interaction analyses with colorectal cancer risk, including interaction tests of genetics(G)xdiabetes (1-degree of freedom; d.f.) and joint testing of Gxdiabetes, G-colorectal cancer association (2-d.f. joint test) and G-diabetes correlation (3-d.f. joint test). RESULTS: Based on the joint tests, we found that the association of diabetes with colorectal cancer risk is modified by loci on chromosomes 8q24.11 (rs3802177, SLC30A8 - ORAA: 1.62, 95% CI: 1.34-1.96; ORAG: 1.41, 95% CI: 1.30-1.54; ORGG: 1.22, 95% CI: 1.13-1.31; p-value3-d.f.: 5.46 × 10-11) and 13q14.13 (rs9526201, LRCH1 - ORGG: 2.11, 95% CI: 1.56-2.83; ORGA: 1.52, 95% CI: 1.38-1.68; ORAA: 1.13, 95% CI: 1.06-1.21; p-value2-d.f.: 7.84 × 10-09). DISCUSSION: These results suggest that variation in genes related to insulin signaling (SLC30A8) and immune function (LRCH1) may modify the association of diabetes with colorectal cancer risk and provide novel insights into the biology underlying the diabetes and colorectal cancer relationship.

Mutations in the Kinesin-2 Motor KIF3B Cause an Autosomal-Dominant Ciliopathy.

Kinesin-2 enables ciliary assembly and maintenance as an anterograde intraflagellar transport (IFT) motor. Molecular motor activity is driven by a heterotrimeric complex comprised of KIF3A and KIF3B or KIF3C plus one non-motor subunit, KIFAP3. Using exome sequencing, we identified heterozygous KIF3B variants in two unrelated families with hallmark ciliopathy phenotypes. In the first family, the proband presents with hepatic fibrosis, retinitis pigmentosa, and postaxial polydactyly; he harbors a de novo c.748G>C (p.Glu250Gln) variant affecting the kinesin motor domain encoded by KIF3B. The second family is a six-generation pedigree affected predominantly by retinitis pigmentosa. Affected individuals carry a heterozygous c.1568T>C (p.Leu523Pro) KIF3B variant segregating in an autosomal-dominant pattern. We observed a significant increase in primary cilia length in vitro in the context of either of the two mutations while variant KIF3B proteins retained stability indistinguishable from wild type. Furthermore, we tested the effects of KIF3B mutant mRNA expression in the developing zebrafish retina. In the presence of either missense variant, rhodopsin was sequestered to the photoreceptor rod inner segment layer with a concomitant increase in photoreceptor cilia length. Notably, impaired rhodopsin trafficking is also characteristic of recessive KIF3B models as exemplified by an early-onset, autosomal-recessive, progressive retinal degeneration in Bengal cats; we identified a c.1000G>A (p.Ala334Thr) KIF3B variant by genome-wide association study and whole-genome sequencing. Together, our genetic, cell-based, and in vivo modeling data delineate an autosomal-dominant syndromic retinal ciliopathy in humans and suggest that multiple KIF3B pathomechanisms can impair kinesin-driven ciliary transport in the photoreceptor.

Deep Learning for Detecting BRCA Mutations in High-Grade Ovarian Cancer Based on an Innovative Tumor Segmentation Method From Whole Slide Images.

BRCA1 and BRCA2 genes play a crucial role in repairing DNA double-strand breaks through homologous recombination. Their mutations represent a significant proportion of homologous recombination deficiency and are a reliable effective predictor of sensitivity of high-grade ovarian cancer (HGOC) to poly(ADP-ribose) polymerase inhibitors. However, their testing by next-generation sequencing is costly and time-consuming and can be affected by various preanalytical factors. In this study, we present a deep learning classifier for BRCA mutational status prediction from hematoxylin-eosin-safran-stained whole slide images (WSI) of HGOC. We constituted the OvarIA cohort composed of 867 patients with HGOC with known BRCA somatic mutational status from 2 different pathology departments. We first developed a tumor segmentation model according to dynamic sampling and then trained a visual representation encoder with momentum contrastive learning on the predicted tumor tiles. We finally trained a BRCA classifier on more than a million tumor tiles in multiple instance learning with an attention-based mechanism. The tumor segmentation model trained on 8 WSI obtained a dice score of 0.915 and an intersection-over-union score of 0.847 on a test set of 50 WSI, while the BRCA classifier achieved the state-of-the-art area under the receiver operating characteristic curve of 0.739 in 5-fold cross-validation and 0.681 on the testing set. An additional multiscale approach indicates that the relevant information for predicting BRCA mutations is located more in the tumor context than in the cell morphology. Our results suggest that BRCA somatic mutations have a discernible phenotypic effect that could be detected by deep learning and could be used as a prescreening tool in the future.

Comprehensive in silico and functional studies for classification of EPAS1/HIF2A genetic variants identified in patients with erythrocytosis.

Gain-of-function mutations in the EPAS1/HIF2A gene have been identified in patients with hereditary erythrocytosis that can be associated with the development of paraganglioma, pheochromocytoma and somatostatinoma. In the present study, we describe a unique European collection of 41 patients and 28 relatives diagnosed with an erythrocytosis associated with a germline genetic variant in EPAS1. In addition we identified two infants with severe erythrocytosis associated with a mosaic mutation present in less than 2% of the blood, one of whom later developed a paraganglioma. The aim of this study was to determine the causal role of these genetic variants, to establish pathogenicity, and to identify potential candidates eligible for the new hypoxia-inducible factor-2 α (HIF-2α) inhibitor treatment. Pathogenicity was predicted with in silico tools and the impact of 13 HIF-2b variants has been studied by using canonical and real-time reporter luciferase assays. These functional assays consisted of a novel edited vector containing an expanded region of the erythropoietin promoter combined with distal regulatory elements which substantially enhanced the HIF-2α-dependent induction. Altogether, our studies allowed the classification of 11 mutations as pathogenic in 17 patients and 23 relatives. We described four new mutations (D525G, L526F, G527K, A530S) close to the key proline P531, which broadens the spectrum of mutations involved in erythrocytosis. Notably, we identified patients with only erythrocytosis associated with germline mutations A530S and Y532C previously identified at somatic state in tumors, thereby raising the complexity of the genotype/phenotype correlations. Altogether, this study allows accurate clinical follow-up of patients and opens the possibility of benefiting from HIF-2α inhibitor treatment, so far the only targeted treatment in hypoxia-related erythrocytosis disease.

Characterization of genetic variants in the EGLN1/PHD2 gene identified in a European collection of patients with erythrocytosis.

Hereditary erythrocytosis is a rare hematologic disorder characterized by an excess of red blood cell production. Here we describe a European collaborative study involving a collection of 2,160 patients with erythrocytosis sequenced in ten different laboratories. We focused our study on the EGLN1 gene and identified 39 germline missense variants including one gene deletion in 47 probands. EGLN1 encodes the PHD2 prolyl 4-hydroxylase, a major inhibitor of hypoxia-inducible factor. We performed a comprehensive study to evaluate the causal role of the identified PHD2 variants: (i) in silico studies of localization, conservation, and deleterious effects; (ii) analysis of hematologic parameters of carriers identified in the UK Biobank; (iii) functional studies of the protein activity and stability; and (iv) a comprehensive study of PHD2 splicing. Altogether, these studies allowed the classification of 16 pathogenic or likely pathogenic mutants in a total of 48 patients and relatives. The in silico studies extended to the variants described in the literature showed that a minority of PHD2 variants can be classified as pathogenic (36/96), without any differences from the variants of unknown significance regarding the severity of the developed disease (hematologic parameters and complications). Here, we demonstrated the great value of federating laboratories working on such rare disorders in order to implement the criteria required for genetic classification, a strategy that should be extended to all hereditary hematologic diseases.

Burden of rare variants in arrhythmogenic cardiomyopathy with right dominant form-associated genes provides new insights for molecular diagnosis and clinical management.

Arrhythmogenic cardiomyopathy with right dominant form (ACR) is a rare heritable cardiac cardiomyopathy disorder associated with sudden cardiac death. Pathogenic variants (PVs) in desmosomal genes have been causally related to ACR in 40% of cases. Other genes encoding nondesmosomal proteins have been described in ACR, but their contribution in this pathology is still debated. A panel of 71 genes associated with inherited cardiopathies was screened in an ACR population of 172 probands and 856 individuals from the general population. PVs and uncertain significance variants (VUS) have been identified in 36% and 18.6% of patients, respectively. Among the cardiopathy-associated genes, burden tests show a significant enrichment in PV and VUS only for desmosomal genes PKP2 (plakophilin-2), DSP (desmoplakin), DSC2 (desmocollin-2), and DSG2 (desmoglein-2). Importantly, VUS may account for 15% of ACR cases and should then be considered for molecular diagnosis. Among the other genes, no evidence of enrichment was detected, suggesting an extreme caution in the interpretation of these genetic variations without associated functional or segregation data. Genotype-phenotype correlation points to (1) a more severe and earlier onset of the disease in PV and VUS carriers, underlying the importance to carry out presymptomatic diagnosis in relatives and (2) to a more prevalent left ventricular dysfunction in DSP variant carriers.

RPL13 Variants Cause Spondyloepimetaphyseal Dysplasia with Severe Short Stature.

Variants in genes encoding ribosomal proteins have thus far been associated with Diamond-Blackfan anemia, a rare inherited bone marrow failure, and isolated congenital asplenia. Here, we report one de novo missense variant and three de novo splice variants in RPL13, which encodes ribosomal protein RPL13 (also called eL13), in four unrelated individuals with a rare bone dysplasia causing severe short stature. The three splice variants (c.477+1G>T, c.477+1G>A, and c.477+2 T>C) result in partial intron retention, which leads to an 18-amino acid insertion. In contrast to observations from Diamond-Blackfan anemia, we detected no evidence of significant pre-rRNA processing disturbance in cells derived from two affected individuals. Consistently, we showed that the insertion-containing protein is stably expressed and incorporated into 60S subunits similar to the wild-type protein. Erythroid proliferation in culture and ribosome profile on sucrose gradient are modified, suggesting a change in translation dynamics. We also provide evidence that RPL13 is present at high levels in chondrocytes and osteoblasts in mouse growth plates. Taken together, we show that the identified RPL13 variants cause a human ribosomopathy defined by a rare skeletal dysplasia, and we highlight the role of this ribosomal protein in bone development.