Swarm Learning for decentralized and confidential clinical machine learning.

Fast and reliable detection of patients with severe and heterogeneous illnesses is a major goal of precision medicine1,2. Patients with leukaemia can be identified using machine learning on the basis of their blood transcriptomes3. However, there is an increasing divide between what is technically possible and what is allowed, because of privacy legislation4,5. Here, to facilitate the integration of any medical data from any data owner worldwide without violating privacy laws, we introduce Swarm Learning-a decentralized machine-learning approach that unites edge computing, blockchain-based peer-to-peer networking and coordination while maintaining confidentiality without the need for a central coordinator, thereby going beyond federated learning. To illustrate the feasibility of using Swarm Learning to develop disease classifiers using distributed data, we chose four use cases of heterogeneous diseases (COVID-19, tuberculosis, leukaemia and lung pathologies). With more than 16,400 blood transcriptomes derived from 127 clinical studies with non-uniform distributions of cases and controls and substantial study biases, as well as more than 95,000 chest X-ray images, we show that Swarm Learning classifiers outperform those developed at individual sites. In addition, Swarm Learning completely fulfils local confidentiality regulations by design. We believe that this approach will notably accelerate the introduction of precision medicine.

HMGB2 Loss upon Senescence Entry Disrupts Genomic Organization and Induces CTCF Clustering across Cell Types.

Processes like cellular senescence are characterized by complex events giving rise to heterogeneous cell populations. However, the early molecular events driving this cascade remain elusive. We hypothesized that senescence entry is triggered by an early disruption of the cells' three-dimensional (3D) genome organization. To test this, we combined Hi-C, single-cell and population transcriptomics, imaging, and in silico modeling of three distinct cells types entering senescence. Genes involved in DNA conformation maintenance are suppressed upon senescence entry across all cell types. We show that nuclear depletion of the abundant HMGB2 protein occurs early on the path to senescence and coincides with the dramatic spatial clustering of CTCF. Knocking down HMGB2 suffices for senescence-induced CTCF clustering and for loop reshuffling, while ectopically expressing HMGB2 rescues these effects. Our data suggest that HMGB2-mediated genomic reorganization constitutes a primer for the ensuing senescent program.

Biallelic loss-of-function variants of ZFTRAF1 cause neurodevelopmental disorder with microcephaly and hypotonia.

PURPOSE: Neurodevelopmental disorders exhibit clinical and genetic heterogeneity, ergo manifest dysfunction in components of diverse cellular pathways; the precise pathomechanism for the majority remains elusive. METHODS: We studied 5 affected individuals from 3 unrelated families manifesting global developmental delay, postnatal microcephaly, and hypotonia. We used exome sequencing and prioritized variants that were subsequently characterized using immunofluorescence, immunoblotting, pulldown assays, and RNA sequencing. RESULTS: We identified biallelic variants in ZFTRAF1, encoding a protein of yet unknown function. Four affected individuals from 2 unrelated families segregated 2 homozygous frameshift variants in ZFTRAF1, whereas, in the third family, an intronic splice site variant was detected. We investigated ZFTRAF1 at the cellular level and signified it as a nucleocytoplasmic protein in different human cell lines. ZFTRAF1 was completely absent in the fibroblasts of 2 affected individuals. We also identified 110 interacting proteins enriched in mRNA processing and autophagy-related pathways. Based on profiling of autophagy markers, patient-derived fibroblasts show irregularities in the protein degradation process. CONCLUSION: Thus, our findings suggest that biallelic variants of ZFTRAF1 cause a severe neurodevelopmental disorder.

The genomic landscape across 474 surgically accessible epileptogenic human brain lesions.

Understanding the exact molecular mechanisms involved in the aetiology of epileptogenic pathologies with or without tumour activity is essential for improving treatment of drug-resistant focal epilepsy. Here, we characterize the landscape of somatic genetic variants in resected brain specimens from 474 individuals with drug-resistant focal epilepsy using deep whole-exome sequencing (>350×) and whole-genome genotyping. Across the exome, we observe a greater number of somatic single-nucleotide variants in low-grade epilepsy-associated tumours (7.92 ± 5.65 single-nucleotide variants) than in brain tissue from malformations of cortical development (6.11 ± 4 single-nucleotide variants) or hippocampal sclerosis (5.1 ± 3.04 single-nucleotide variants). Tumour tissues also had the largest number of likely pathogenic variant carrying cells. low-grade epilepsy-associated tumours had the highest proportion of samples with one or more somatic copy-number variants (24.7%), followed by malformations of cortical development (5.4%) and hippocampal sclerosis (4.1%). Recurring somatic whole chromosome duplications affecting Chromosome 7 (16.8%), chromosome 5 (10.9%), and chromosome 20 (9.9%) were observed among low-grade epilepsy-associated tumours. For germline variant-associated malformations of cortical development genes such as TSC2, DEPDC5 and PTEN, germline single-nucleotide variants were frequently identified within large loss of heterozygosity regions, supporting the recently proposed 'second hit' disease mechanism in these genes. We detect somatic variants in 12 established lesional epilepsy genes and demonstrate exome-wide statistical support for three of these in the aetiology of low-grade epilepsy-associated tumours (e.g. BRAF) and malformations of cortical development (e.g. SLC35A2 and MTOR). We also identify novel significant associations for PTPN11 with low-grade epilepsy-associated tumours and NRAS Q61 mutated protein with a complex malformation of cortical development characterized by polymicrogyria and nodular heterotopia. The variants identified in NRAS are known from cancer studies to lead to hyperactivation of NRAS, which can be targeted pharmacologically. We identify large recurrent 1q21-q44 duplication including AKT3 in association with focal cortical dysplasia type 2a with hyaline astrocytic inclusions, another rare and possibly under-recognized brain lesion. The clinical-genetic analyses showed that the numbers of somatic single-nucleotide variant across the exome and the fraction of affected cells were positively correlated with the age at seizure onset and surgery in individuals with low-grade epilepsy-associated tumours. In summary, our comprehensive genetic screen sheds light on the genome-scale landscape of genetic variants in epileptic brain lesions, informs the design of gene panels for clinical diagnostic screening and guides future directions for clinical implementation of epilepsy surgery genetics.

Genomic ALK alterations in primary and relapsed neuroblastoma.

BACKGROUND: Genomic alterations of the anaplastic lymphoma kinase gene (ALK) occur recurrently in neuroblastoma, a pediatric malignancy of the sympathetic nervous system. However, information on their development over time has remained sparse. METHODS: ALK alterations were assessed in neuroblastomas at diagnosis and/or relapse from a total of 943 patients, covering all stages of disease. Longitudinal information on diagnostic and relapsed samples from individual patients was available in 101 and 102 cases for mutation and amplification status, respectively. RESULTS: At diagnosis, ALK point mutations occurred in 10.5% of all cases, with highest frequencies in stage 4 patients <18 months. At relapse, ALK alteration frequency increased by 70%, both in high-risk and non-high-risk cases. The increase was most likely due to de novo mutations, frequently leading to R1275Q substitutions, which are sensitive to pharmacological ALK inhibition. By contrast, the frequency of ALK amplifications did not change over the course of the disease. ALK amplifications, but not mutations, were associated with poor patient outcome. CONCLUSIONS: The considerably increased frequency of ALK mutations at relapse and their high prevalence in young stage 4 patients suggest surveying the genomic ALK status regularly in these patient cohorts, and to evaluate ALK-targeted treatment also in intermediate-risk patients.

Homozygous loss-of-function variants in FILIP1 cause autosomal recessive arthrogryposis multiplex congenita with microcephaly.

Arthrogryposis multiplex congenita forms a broad group of clinically and etiologically heterogeneous disorders characterized by congenital joint contractures that involve at least two different parts of the body. Neurological and muscular disorders are commonly underlying arthrogryposis. Here, we report five affected individuals from three independent families sharing an overlapping phenotype with congenital contractures affecting shoulder, elbow, hand, hip, knee and foot as well as scoliosis, reduced palmar and plantar skin folds, microcephaly and facial dysmorphism. Using exome sequencing, we identified homozygous truncating variants in FILIP1 in all patients. FILIP1 is a regulator of filamin homeostasis required for the initiation of cortical cell migration in the developing neocortex and essential for the differentiation process of cross-striated muscle cells during myogenesis. In summary, our data indicate that bi-allelic truncating variants in FILIP1 are causative of a novel autosomal recessive disorder and expand the spectrum of genetic factors causative of arthrogryposis multiplex congenita.

Identification of pathogenic variant enriched regions across genes and gene families.

Missense variant interpretation is challenging. Essential regions for protein function are conserved among gene-family members, and genetic variants within these regions are potentially more likely to confer risk to disease. Here, we generated 2871 gene-family protein sequence alignments involving 9990 genes and performed missense variant burden analyses to identify novel essential protein regions. We mapped 2,219,811 variants from the general population into these alignments and compared their distribution with 76,153 missense variants from patients. With this gene-family approach, we identified 465 regions enriched for patient variants spanning 41,463 amino acids in 1252 genes. As a comparison, by testing the same genes individually, we identified fewer patient variant enriched regions, involving only 2639 amino acids and 215 genes. Next, we selected de novo variants from 6753 patients with neurodevelopmental disorders and 1911 unaffected siblings and observed an 8.33-fold enrichment of patient variants in our identified regions (95% C.I. = 3.90-Inf, P-value = 2.72 × 10-11). By using the complete ClinVar variant set, we found that missense variants inside the identified regions are 106-fold more likely to be classified as pathogenic in comparison to benign classification (OR = 106.15, 95% C.I = 70.66-Inf, P-value < 2.2 × 10-16). All pathogenic variant enriched regions (PERs) identified are available online through "PER viewer," a user-friendly online platform for interactive data mining, visualization, and download. In summary, our gene-family burden analysis approach identified novel PERs in protein sequences. This annotation can empower variant interpretation.

Expanding the phenotypic spectrum and clinical severity associated with WLS gene.

WLS (Wnt ligand secretion mediator or Wntless) orchestrates the secretion of all Wnt proteins, a family of evolutionary conserved proteins, involved in Wnt signaling pathway that has many essential biological functions including the regulation of development, cell proliferation, migration and apoptosis. Biallelic variants in WLS have recently been described in 10 patients with pleiotropic multiple congenital anomalies (MCA) known as Zaki syndrome. We identified a likely disease-causing variant in WLS (c.1579G>A, p.Gly527Arg) in a boy presented with a broad range of MCA including microcephaly, facial dysmorphism, alopecia, ophthalmologic anomalies, and complete soft tissue syndactyly. These features were reminiscent of Zaki syndrome although variable clinical severity was observed. In a detailed clinical assessment, our patient also displayed microphthalmia, dental anomalies, skeletal dysplasia with spontaneous fractures and Dandy-Walker malformation. As such, we extend the phenotype linked to Zaki syndrome. This study further highlights the importance of a thorough clinical evaluation to delineate the phenotypic spectrum associated with WLS variants and suggests that genotype-phenotype correlations due to variant localization seems likely. However, future work on additional patients and more functional studies may give further insights into genotype-phenotype correlations and the complex function of WLS.

Homozygous frameshift variant in desmoglein 2 gene causes biventricular arrhythmogenic right ventricular cardiomyopathy.

Arrhythmogenic right ventricular cardiomyopathy caused by a homozygous frameshift variant (c.2358delA) in DSG2.

Expanding the spectrum of EEF1D neurodevelopmental disorders: Biallelic variants in the guanine exchange domain.

Protein translation is an essential cellular process and dysfunctional protein translation causes various neurodevelopmental disorders. The eukaryotic translation elongation factor 1A (eEF1A) delivers aminoacyl-tRNA to the ribosome, while the eEF1B complex acts as a guanine exchange factor (GEF) of GTP for GDP indirectly catalyzing the release of eEF1A from the ribosome. The gene EEF1D encodes the eEF1Bδ subunit of the eEF1B complex. EEF1D is alternatively spliced giving rise to one long and three short isoforms. Two different homozygous, truncating variants in EEF1D had been associated with severe intellectual disability and microcephaly in two families. The published variants only affect the long isoform of EEF1D that acts as a transcription factor of heat shock element proteins. By exome sequencing, we identified two different homozygous variants in EEF1D in two families with severe developmental delay, severe microcephaly, spasticity, and failure to thrive with optic atrophy, poor feeding, and recurrent aspiration pneumonia. The EEF1D variants reported in this study are localized in the C-terminal GEF domain, suggesting that a disturbed protein translation machinery might contribute to the neurodevelopmental phenotype. Pathogenic variants localized in both the alternatively spliced domain or the GEF domain of EEF1D cause a severe neurodevelopmental disorder with microcephaly and spasticity.

Ganglioglioma with adverse clinical outcome and atypical histopathological features were defined by alterations in PTPN11/KRAS/NF1 and other RAS-/MAP-Kinase pathway genes.

Exome-wide sequencing studies recently described PTPN11 as a novel brain somatic epilepsy gene. In contrast, germline mutations of PTPN11 are known to cause Noonan syndrome, a multisystem disorder characterized by abnormal facial features, developmental delay, and sporadically, also brain tumors. Herein, we performed a deep phenotype-genotype analysis of a comprehensive series of ganglioglioma (GG) with brain somatic alterations of the PTPN11/KRAS/NF1 genes compared to GG with common MAP-Kinase signaling pathway alterations, i.e., BRAFV600E. Seventy-two GG were submitted to whole exome sequencing and genotyping and 84 low grade epilepsy associated tumors (LEAT) to DNA-methylation analysis. In 28 tumours, both analyses were available from the same sample. Clinical data were retrieved from hospital files including disease onset, age at surgery, brain localization, and seizure outcome. A comprehensive histopathology staining panel was available in all cases. We identified eight GG with PTPN11 alterations, copy number variant (CNV) gains of chromosome 12, and the commonality of additional CNV gains in NF1, KRAS, FGFR4 and RHEB, as well as BRAFV600E alterations. Histopathology revealed an atypical glio-neuronal phenotype with subarachnoidal tumor spread and large, pleomorphic, and multinuclear cellular features. Only three out of eight patients with GG and PTPN11/KRAS/NF1 alterations were free of disabling-seizures 2 years after surgery (38% had Engel I). This was remarkably different from our series of GG with only BRAFV600E mutations (85% had Engel I). Unsupervised cluster analysis of DNA methylation arrays separated these tumours from well-established LEAT categories. Our data point to a subgroup of GG with cellular atypia in glial and neuronal cell components, adverse postsurgical outcome, and genetically characterized by complex alterations in PTPN11 and other RAS-/MAP-Kinase and/or mTOR signaling pathways. These findings need prospective validation in clinical practice as they argue for an adaptation of the WHO grading system in developmental, glio-neuronal tumors associated with early onset focal epilepsy.

Mutations in TAF8 cause a neurodegenerative disorder.

TAF8 is part of the transcription factor II D complex, composed of the TATA-binding protein and 13 TATA-binding protein-associated factors (TAFs). Transcription factor II D is the first general transcription factor recruited at promoters to assemble the RNA polymerase II preinitiation complex. So far disorders related to variants in 5 of the 13 subunits of human transcription factor II D have been described. Recently, a child with a homozygous c.781-1G>A mutation in TAF8 has been reported. Here we describe seven further patients with mutations in TAF8 and thereby confirm the TAF8 related disorder. In two sibling patients, we identified two novel compound heterozygous TAF8 splice site mutations, c.45+4A > G and c.489G>A, which cause aberrant splicing as well as reduced expression and mislocalization of TAF8. In five further patients, the previously described c.781-1G > A mutation was present on both alleles. The clinical phenotype associated with the different TAF8 mutations is characterized by severe psychomotor retardation with almost absent development, feeding problems, microcephaly, growth retardation, spasticity and epilepsy. Cerebral imaging showed hypomyelination, a thin corpus callosum and brain atrophy. Moreover, repeated imaging in the sibling pair demonstrated progressive cerebral and cerebellar atrophy. Consistently, reduced N-acetylaspartate, a marker of neuronal viability, was observed on magnetic resonance spectroscopy. Further review of the literature shows that mutations causing a reduced expression of transcription factor II D subunits have an overlapping phenotype of microcephaly, developmental delay and intellectual disability. Although transcription factor II D plays an important role in RNA polymerase II transcription in all cells and tissues, the symptoms associated with such defects are almost exclusively neurological. This might indicate a specific vulnerability of neuronal tissue to widespread deregulation of gene expression as also seen in Rett syndrome or Cornelia de Lange syndrome.

Loss-of-function variants in DNM1 cause a specific form of developmental and epileptic encephalopathy only in biallelic state.

BACKGROUND: Developmental and epileptic encephalopathies (DEEs) represent a group of severe neurological disorders characterised by an onset of refractory seizures during infancy or early childhood accompanied by psychomotor developmental delay or regression. DEEs are genetically heterogeneous with, to date, more than 80 different genetic subtypes including DEE31 caused by heterozygous missense variants in DNM1. METHODS: We performed a detailed clinical characterisation of two unrelated patients with DEE and used whole-exome sequencing to identify causative variants in these individuals. The identified variants were tested for cosegregation in the respective families. RESULTS: We excluded pathogenic variants in known, DEE-associated genes. We identified homozygous nonsense variants, c.97C>T; p.(Gln33*) in family 1 and c.850C>T; p.(Gln284*) in family 2, in the DNM1 gene, indicating that biallelic, loss-of-function pathogenic variants in DNM1 cause DEE. CONCLUSION: Our finding that homozygous, loss-of-function variants in DNM1 cause DEE expands the spectrum of pathogenic variants in DNM1. All parents who were heterozygous carriers of the identified loss-of-function variants were healthy and did not show any clinical symptoms, indicating that the type of mutation in DNM1 determines the pattern of inheritance.

Sensitivity and specificity of loss of heterozygosity analysis for the classification of rare germline variants in BRCA1/2: results of the observational AGO-TR1 study (NCT02222883).

Variant-specific loss of heterozygosity (LOH) analyses may be useful to classify BRCA1/2 germline variants of unknown significance (VUS). The sensitivity and specificity of this approach, however, remains unknown. We performed comparative next-generation sequencing analyses of the BRCA1/2 genes using blood-derived and tumour-derived DNA of 488 patients with ovarian cancer enrolled in the observational AGO-TR1 trial (NCT02222883). Overall, 94 pathogenic, 90 benign and 24 VUS were identified in the germline. A significantly increased variant fraction (VF) of a germline variant in the tumour indicates loss of the wild-type allele; a decreased VF indicates loss of the variant allele. We demonstrate that significantly increased VFs predict pathogenicity with high sensitivity (0.84, 95% CI 0.77 to 0.91), poor specificity (0.63, 95% CI 0.53 to 0.73) and poor positive predictive value (PPV; 0.71, 95% CI 0.62 to 0.79). Significantly decreased VFs predict benignity with low sensitivity (0.26, 95% CI 0.17 to 0.35), high specificity (1.0, 95% CI 0.96 to 1.00) and PPV (1.0, 95% CI 0.85 to 1.00). Variant classification based on significantly increased VFs results in an unacceptable proportion of false-positive results. A significantly decreased VF in the tumour may be exploited as a reliable predictor for benignity, with no false-negative result observed. When applying the latter approach, VUS identified in four patients can now be considered benign. Trial registration number NCT02222883.

Pseudouridylation defect due to DKC1 and NOP10 mutations causes nephrotic syndrome with cataracts, hearing impairment, and enterocolitis.

RNA modifications play a fundamental role in cellular function. Pseudouridylation, the most abundant RNA modification, is catalyzed by the H/ACA small ribonucleoprotein (snoRNP) complex that shares four core proteins, dyskerin (DKC1), NOP10, NHP2, and GAR1. Mutations in DKC1, NOP10, or NHP2 cause dyskeratosis congenita (DC), a disorder characterized by telomere attrition. Here, we report a phenotype comprising nephrotic syndrome, cataracts, sensorineural deafness, enterocolitis, and early lethality in two pedigrees: males with DKC1 p.Glu206Lys and two children with homozygous NOP10 p.Thr16Met. Females with heterozygous DKC1 p.Glu206Lys developed cataracts and sensorineural deafness, but nephrotic syndrome in only one case of skewed X-inactivation. We found telomere attrition in both pedigrees, but no mucocutaneous abnormalities suggestive of DC. Both mutations fall at the dyskerin-NOP10 binding interface in a region distinct from those implicated in DC, impair the dyskerin-NOP10 interaction, and disrupt the catalytic pseudouridylation site. Accordingly, we found reduced pseudouridine levels in the ribosomal RNA (rRNA) of the patients. Zebrafish dkc1 mutants recapitulate the human phenotype and show reduced 18S pseudouridylation, ribosomal dysregulation, and a cell-cycle defect in the absence of telomere attrition. We therefore propose that this human disorder is the consequence of defective snoRNP pseudouridylation and ribosomal dysfunction.

WARS1 and SARS1: Two tRNA synthetases implicated in autosomal recessive microcephaly.

Aminoacylation of transfer RNA (tRNA) is a key step in protein biosynthesis, carried out by highly specific aminoacyl-tRNA synthetases (ARSs). ARSs have been implicated in autosomal dominant and autosomal recessive human disorders. Autosomal dominant variants in tryptophanyl-tRNA synthetase 1 (WARS1) are known to cause distal hereditary motor neuropathy and Charcot-Marie-Tooth disease, but a recessively inherited phenotype is yet to be clearly defined. Seryl-tRNA synthetase 1 (SARS1) has rarely been implicated in an autosomal recessive developmental disorder. Here, we report five individuals with biallelic missense variants in WARS1 or SARS1, who presented with an overlapping phenotype of microcephaly, developmental delay, intellectual disability, and brain anomalies. Structural mapping showed that the SARS1 variant is located directly within the enzyme's active site, most likely diminishing activity, while the WARS1 variant is located in the N-terminal domain. We further characterize the identified WARS1 variant by showing that it negatively impacts protein abundance and is unable to rescue the phenotype of a CRISPR/Cas9 wars1 knockout zebrafish model. In summary, we describe two overlapping autosomal recessive syndromes caused by variants in WARS1 and SARS1, present functional insights into the pathogenesis of the WARS1-related syndrome and define an emerging disease spectrum: ARS-related developmental disorders with or without microcephaly.

Transmission ratio distortion of mutations in the master regulator of centriole biogenesis PLK4.

The evolutionary conserved Polo-like kinase 4 (PLK4) is essential for centriole duplication, spindle assembly, and de novo centriole formation. In man, homozygous mutations in PLK4 lead to primary microcephaly, altered PLK4 expression is associated with aneuploidy in human embryos. Here, we report on a consanguineous four-generation family with 8 affected individuals compound heterozygous for a novel missense variant, c.881 T > G, and a deletion of the PLK4 gene. The clinical phenotype of the adult patients is mild compared to individuals with previously described PLK4 mutations. One individual was homozygous for the variant c.881G and phenotypically unaffected. The deletion was inherited by 14 of 16 offspring and thus exhibits transmission ratio distortion (TRD). Moreover, based on the already published families with PLK4 mutations, it could be shown that due to the preferential transmission of the mutant alleles, the number of affected offspring is significantly increased. It is assumed that reduced expression of PLK4 decreases the intrinsically high error rate of the first cell divisions after fertilization, increases the number of viable embryos and thus leads to preferential transmission of the deleted/mutated alleles.

Autosomal-Recessive Mutations in MESD Cause Osteogenesis Imperfecta.

Osteogenesis imperfecta (OI) comprises a genetically heterogeneous group of skeletal fragility diseases. Here, we report on five independent families with a progressively deforming type of OI, in whom we identified four homozygous truncation or frameshift mutations in MESD. Affected individuals had recurrent fractures and at least one had oligodontia. MESD encodes an endoplasmic reticulum (ER) chaperone protein for the canonical Wingless-related integration site (WNT) signaling receptors LRP5 and LRP6. Because complete absence of MESD causes embryonic lethality in mice, we hypothesized that the OI-associated mutations are hypomorphic alleles since these mutations occur downstream of the chaperone activity domain but upstream of ER-retention domain. This would be consistent with the clinical phenotypes of skeletal fragility and oligodontia in persons deficient for LRP5 and LRP6, respectively. When we expressed wild-type (WT) and mutant MESD in HEK293T cells, we detected WT MESD in cell lysate but not in conditioned medium, whereas the converse was true for mutant MESD. We observed that both WT and mutant MESD retained the ability to chaperone LRP5. Thus, OI-associated MESD mutations produce hypomorphic alleles whose failure to remain within the ER significantly reduces but does not completely eliminate LRP5 and LRP6 trafficking. Since these individuals have no eye abnormalities (which occur in individuals completely lacking LRP5) and have neither limb nor brain patterning defects (both of which occur in mice completely lacking LRP6), we infer that bone mass accrual and dental patterning are more sensitive to reduced canonical WNT signaling than are other developmental processes. Biologic agents that can increase LRP5 and LRP6-mediated WNT signaling could benefit individuals with MESD-associated OI.

Monoallelic and biallelic variants in LEF1 are associated with a new syndrome combining ectodermal dysplasia and limb malformations caused by altered WNT signaling.

PURPOSE: LEF1 encodes a transcription factor acting downstream of the WNT-ß-catenin signaling pathway. It was recently suspected as a candidate for ectodermal dysplasia in 2 individuals carrying 4q35 microdeletions. We report on 12 individuals harboring LEF1 variants. METHODS: High-throughput sequencing was employed to delineate the genetic underpinnings of the disease. Cellular consequences were characterized by immunofluorescence, immunoblotting, pulldown assays, and/or RNA sequencing. RESULTS: Monoallelic variants in LEF1 were detected in 11 affected individuals from 4 unrelated families, and a biallelic variant was detected in an affected individual from a consanguineous family. The phenotypic spectrum includes various limb malformations, such as radial ray defects, polydactyly or split hand/foot, and ectodermal dysplasia. Depending on the type and location of LEF1 variants, the inheritance of this novel Mendelian condition can be either autosomal dominant or recessive. Our functional data indicate that 2 molecular mechanisms are at play: haploinsufficiency or loss of DNA binding are responsible for a mild to moderate phenotype, whereas loss of ß-catenin binding caused by biallelic variants is associated with a severe phenotype. Transcriptomic studies reveal an alteration of WNT signaling. CONCLUSION: Our findings establish mono- and biallelic variants in LEF1 as a cause for a novel syndrome comprising limb malformations and ectodermal dysplasia.