Ultrarare Variants in DNA Damage Repair Genes in Pediatric Acute-Onset Neuropsychiatric Syndrome or Acute Behavioral Regression in Neurodevelopmental Disorders.

Acute onset of severe psychiatric symptoms or regression may occur in children with premorbid neurodevelopmental disorders, although typically developing children can also be affected. Infections or other stressors are likely triggers. The underlying causes are unclear, but a current hypothesis suggests the convergence of genes that influence neuronal and immunological function. We previously identified 11 genes in Pediatric Acute-Onset Neuropsychiatry Syndrome (PANS), in which two classes of genes related to either synaptic function or the immune system were found. Among the latter, three affect the DNA damage response (DDR): PPM1D, CHK2, and RAG1. We now report an additional 17 cases with mutations in PPM1D and other DDR genes in patients with acute onset of psychiatric symptoms and/or regression that their clinicians classified as PANS or another inflammatory brain condition. The genes include clusters affecting p53 DNA repair (PPM1D, ATM, ATR, 53BP1, and RMRP), and the Fanconi Anemia Complex (FANCE, SLX4/FANCP, FANCA, FANCI, and FANCC). We hypothesize that defects in DNA repair genes, in the context of infection or other stressors, could contribute to decompensated states through an increase in genomic instability with a concomitant accumulation of cytosolic DNA in immune cells triggering DNA sensors, such as cGAS-STING and AIM2 inflammasomes, as well as central deficits on neuroplasticity. In addition, increased senescence and defective apoptosis affecting immunological responses could be playing a role. These compelling preliminary findings motivate further genetic and functional characterization as the downstream impact of DDR deficits may point to novel treatment strategies.

A serum B-lymphocyte activation signature is a key distinguishing feature of the immune response in sarcoidosis compared to tuberculosis.

Sarcoidosis and tuberculosis (TB) are two granulomatous diseases that often share overlapping clinical features, including uveitis. We measured 368 inflammation-related proteins in serum in both diseases, with and without uveitis from two distinct geographically separated cohorts: sarcoidosis from the Netherlands and TB from Indonesia. A total of 192 and 102 differentially expressed proteins were found in sarcoidosis and active pulmonary TB compared to their geographical healthy controls, respectively. While substantial overlap exists in the immune-related pathways involved in both diseases, activation of B cell activating factor (BAFF) signaling and proliferation-inducing ligand (APRIL) mediated signaling pathways was specifically associated with sarcoidosis. We identified a B-lymphocyte activation signature consisting of BAFF, TNFRSF13B/TACI, TRAF2, IKBKG, MAPK9, NFATC1, and DAPP1 that was associated with sarcoidosis, regardless of the presence of uveitis. In summary, a difference in B-lymphocyte activation is a key discriminative immunological feature between sarcoidosis/ocular sarcoidosis (OS) and TB/ocular TB (OTB).

Bioinformatic meta-analysis reveals novel differentially expressed genes and pathways in sarcoidosis.

Introduction: Sarcoidosis is a multi-system inflammatory disease of unknown origin with heterogeneous clinical manifestations varying from a single organ non-caseating granuloma site to chronic systemic inflammation and fibrosis. Gene expression studies have suggested several genes and pathways implicated in the pathogenesis of sarcoidosis, however, due to differences in study design and variable statistical approaches, results were frequently not reproducible or concordant. Therefore, meta-analysis of sarcoidosis gene-expression datasets is of great importance to robustly establish differentially expressed genes and signalling pathways. Methods: We performed meta-analysis on 22 published gene-expression studies on sarcoidosis. Datasets were analysed systematically using same statistical cut-offs. Differentially expressed genes were identified by pooling of p-values using Edgington's method and analysed for pathways using Ingenuity Pathway Analysis software. Results: A consistent and significant signature of novel and well-known genes was identified, those collectively implicated both type I and type II interferon mediated signalling pathways in sarcoidosis. In silico functional analysis showed consistent downregulation of eukaryotic initiation factor 2 signalling, whereas cytokines like interferons and transcription factor STAT1 were upregulated. Furthermore, we analysed affected tissues to detect differentially expressed genes likely to be involved in granuloma biology. This revealed that matrix metallopeptidase 12 was exclusively upregulated in affected tissues, suggesting a crucial role in disease pathogenesis. Discussion: Our analysis provides a concise gene signature in sarcoidosis and expands our knowledge about the pathogenesis. Our results are of importance to improve current diagnostic approaches and monitoring strategies as well as in the development of targeted therapeutics.

Plasticity mechanisms of genetically distinct Purkinje cells.

Despite its uniform appearance, the cerebellar cortex is highly heterogeneous in terms of structure, genetics and physiology. Purkinje cells (PCs), the principal and sole output neurons of the cerebellar cortex, can be categorized into multiple populations that differentially express molecular markers and display distinctive physiological features. Such features include action potential rate, but also their propensity for synaptic and intrinsic plasticity. However, the precise molecular and genetic factors that correlate with the differential physiological properties of PCs remain elusive. In this article, we provide a detailed overview of the cellular mechanisms that regulate PC activity and plasticity. We further perform a pathway analysis to highlight how molecular characteristics of specific PC populations may influence their physiology and plasticity mechanisms.

Whole exome sequencing of patients with varicella-zoster virus and herpes simplex virus induced acute retinal necrosis reveals rare disease-associated genetic variants.

Purpose: Herpes simplex virus (HSV) and varicella-zoster virus (VZV) are neurotropic human alphaherpesviruses endemic worldwide. Upon primary infection, both viruses establish lifelong latency in neurons and reactivate intermittently to cause a variety of mild to severe diseases. Acute retinal necrosis (ARN) is a rare, sight-threatening eye disease induced by ocular VZV or HSV infection. The virus and host factors involved in ARN pathogenesis remain incompletely described. We hypothesize an underlying genetic defect in at least part of ARN cases. Methods: We collected blood from 17 patients with HSV-or VZV-induced ARN, isolated DNA and performed Whole Exome Sequencing by Illumina followed by analysis in Varseq with criteria of CADD score > 15 and frequency in GnomAD < 0.1% combined with biological filters. Gene modifications relative to healthy control genomes were filtered according to high quality and read-depth, low frequency, high deleteriousness predictions and biological relevance. Results: We identified a total of 50 potentially disease-causing genetic variants, including missense, frameshift and splice site variants and on in-frame deletion in 16 of the 17 patients. The vast majority of these genes are involved in innate immunity, followed by adaptive immunity, autophagy, and apoptosis; in several instances variants within a given gene or pathway was identified in several patients. Discussion: We propose that the identified variants may contribute to insufficient viral control and increased necrosis ocular disease presentation in the patients and serve as a knowledge base and starting point for the development of improved diagnostic, prophylactic, and therapeutic applications.

Pathogenic variants in the paired-related homeobox 1 gene (PRRX1) cause craniosynostosis with incomplete penetrance.

PURPOSE: Studies have previously implicated PRRX1 in craniofacial development, including demonstration of murine Prrx1 expression in the preosteogenic cells of the cranial sutures. We investigated the role of heterozygous missense and loss-of-function (LoF) variants in PRRX1 associated with craniosynostosis. METHODS: Trio-based genome, exome, or targeted sequencing were used to screen PRRX1 in patients with craniosynostosis; immunofluorescence analyses were used to assess nuclear localization of wild-type and mutant proteins. RESULTS: Genome sequencing identified 2 of 9 sporadically affected individuals with syndromic/multisuture craniosynostosis, who were heterozygous for rare/undescribed variants in PRRX1. Exome or targeted sequencing of PRRX1 revealed a further 9 of 1449 patients with craniosynostosis harboring deletions or rare heterozygous variants within the homeodomain. By collaboration, 7 additional individuals (4 families) were identified with putatively pathogenic PRRX1 variants. Immunofluorescence analyses showed that missense variants within the PRRX1 homeodomain cause abnormal nuclear localization. Of patients with variants considered likely pathogenic, bicoronal or other multisuture synostosis was present in 11 of 17 cases (65%). Pathogenic variants were inherited from unaffected relatives in many instances, yielding a 12.5% penetrance estimate for craniosynostosis. CONCLUSION: This work supports a key role for PRRX1 in cranial suture development and shows that haploinsufficiency of PRRX1 is a relatively frequent cause of craniosynostosis.

A germline STAT6 gain-of-function variant is associated with early-onset allergies.

BACKGROUND: The signal transducer and activator of transcription 6 (STAT6) signaling pathway plays a central role in allergic inflammation. To date, however, there have been no descriptions of STAT6 gain-of-function variants leading to allergies in humans. OBJECTIVE: We report a STAT6 gain-of-function variant associated with early-onset multiorgan allergies in a family with 3 affected members. METHODS: Exome sequencing and immunophenotyping of T-helper cell subsets were conducted. The function of the STAT6 protein was analyzed by Western blot, immunofluorescence, electrophoretic mobility shift assays, and luciferase assays. Gastric organoids obtained from the index patient were used to study downstream effector cytokines. RESULTS: We identified a heterozygous missense variant (c.1129G>A;p.Glu377Lys) in the DNA binding domain of STAT6 that was de novo in the index patient's father and was inherited by 2 of his 3 children. Severe atopic dermatitis and food allergy were key presentations. Clinical heterogeneity was observed among the affected individuals. Higher levels of peripheral blood TH2 lymphocytes were detected. The mutant STAT6 displayed a strong preference for nuclear localization, increased DNA binding affinity, and spontaneous transcriptional activity. Moreover, gastric organoids showed constitutive activation of STAT6 downstream signaling molecules. CONCLUSIONS: A germline STAT6 gain-of-function variant results in spontaneous activation of the STAT6 signaling pathway and is associated with an early-onset and severe allergic phenotype in humans. These observations enhance our knowledge of the molecular mechanisms underlying allergic diseases and will potentially contribute to novel therapeutic interventions.

Missense variants in ANKRD11 cause KBG syndrome by impairment of stability or transcriptional activity of the encoded protein.

PURPOSE: Although haploinsufficiency of ANKRD11 is among the most common genetic causes of neurodevelopmental disorders, the role of rare ANKRD11 missense variation remains unclear. We characterized clinical, molecular, and functional spectra of ANKRD11 missense variants. METHODS: We collected clinical information of individuals with ANKRD11 missense variants and evaluated phenotypic fit to KBG syndrome. We assessed pathogenicity of variants through in silico analyses and cell-based experiments. RESULTS: We identified 20 unique, mostly de novo, ANKRD11 missense variants in 29 individuals, presenting with syndromic neurodevelopmental disorders similar to KBG syndrome caused by ANKRD11 protein truncating variants or 16q24.3 microdeletions. Missense variants significantly clustered in repression domain 2 at the ANKRD11 C-terminus. Of the 10 functionally studied missense variants, 6 reduced ANKRD11 stability. One variant caused decreased proteasome degradation and loss of ANKRD11 transcriptional activity. CONCLUSION: Our study indicates that pathogenic heterozygous ANKRD11 missense variants cause the clinically recognizable KBG syndrome. Disrupted transrepression capacity and reduced protein stability each independently lead to ANKRD11 loss-of-function, consistent with haploinsufficiency. This highlights the diagnostic relevance of ANKRD11 missense variants, but also poses diagnostic challenges because the KBG-associated phenotype may be mild and inherited pathogenic ANKRD11 (missense) variants are increasingly observed, warranting stringent variant classification and careful phenotyping.

CRISPRs in the human genome are differentially expressed between malignant and normal adjacent to tumor tissue.

Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs) have been identified in bacteria, archaea and mitochondria of plants, but not in eukaryotes. Here, we report the discovery of 12,572 putative CRISPRs randomly distributed across the human chromosomes, which we termed hCRISPRs. By using available transcriptome datasets, we demonstrate that hCRISPRs are distinctively expressed as small non-coding RNAs (sncRNAs) in cell lines and human tissues. Moreover, expression patterns thereof enabled us to distinguish normal from malignant tissues. In prostate cancer, we confirmed the differential hCRISPR expression between normal adjacent and malignant primary prostate tissue by RT-qPCR and demonstrate that the SHERLOCK and DETECTR dipstick tools are suitable to detect these sncRNAs. We anticipate that the discovery of CRISPRs in the human genome can be further exploited for diagnostic purposes in cancer and other medical conditions, which certainly will lead to the development of point-of-care tests based on the differential expression of the hCRISPRs.

Identification of a novel variant of the ciliopathic gene FUZZY associated with craniosynostosis.

Craniosynostosis is a birth defect occurring in approximately one in 2000 live births, where premature fusion of the cranial bones inhibits growth of the skull during critical periods of brain development. The resulting changes in skull shape can lead to compression of the brain, causing severe complications. While we have some understanding of the molecular pathology of craniosynostosis, a large proportion of cases are of unknown genetic aetiology. Based on studies in mouse, we previously proposed that the ciliopathy gene Fuz should be considered a candidate craniosynostosis gene. Here, we report a novel variant of FUZ (c.851 G > C, p.(Arg284Pro)) found in monozygotic twins presenting with craniosynostosis. To investigate whether Fuz has a direct role in regulating osteogenic fate and mineralisation, we cultured primary osteoblasts and mouse embryonic fibroblasts (MEFs) from Fuz mutant mice. Loss of Fuz resulted in increased osteoblastic mineralisation. This suggests that FUZ protein normally acts as a negative regulator of osteogenesis. We then used Fuz mutant MEFs, which lose functional primary cilia, to test whether the FUZ p.(Arg284Pro) variant could restore FUZ function during ciliogenesis. We found that expression of the FUZ p.(Arg284Pro) variant was sufficient to partially restore cilia numbers, but did not mediate a comparable response to Hedgehog pathway activation. Together, this suggests the osteogenic effects of FUZ p.(Arg284Pro) do not depend upon initiation of ciliogenesis.

Vitreous proteomics, a gateway to improved understanding and stratification of diverse uveitis aetiologies.

PURPOSE: The vitreous proteome might provide an attractive gateway to discriminate between various uveitis aetiologies and gain novel insights into the underlying pathophysiological processes. Here, we investigated 180 vitreous proteins to discover novel biomarkers and broaden disease insights by comparing (1). primary vitreoretinal lymphoma ((P)VRL) versus other aetiologies, (2). sarcoid uveitis versus tuberculosis (TB)-associated uveitis and (3). granulomatous (sarcoid and TB) uveitis versus other aetiologies. METHODS: Vitreous protein levels were determined by proximity extension assay in 47 patients with intraocular inflammation and a prestudy diagnosis (cohort 1; training) and 22 patients with a blinded diagnosis (cohort 2; validation). Differentially expressed proteins identified by t-tests on cohort 1 were used to calculate Youden's indices. Pathway and network analysis was performed by ingenuity pathway analysis. A random forest classifier was trained to predict the diagnosis of blinded patients. RESULTS: For (P)VRL stratification, the previously reported combined diagnostic value of IL-10 and IL-6 was confirmed. Additionally, CD70 was identified as potential novel marker for (P)VRL. However, the classifier trained on the entire cohort (cohort 1 and 2) relied primarily on the interleukin score for intraocular lymphoma diagnosis (ISOLD) or IL-10/IL-6 ratio and only showed a supportive role for CD70. Furthermore, sarcoid uveitis displayed increased levels of vitreous CCL17 as compared to TB-associated uveitis. CONCLUSION: We underline the previously reported value of the ISOLD and the IL-10/IL-6 ratio for (P)VRL identification and present CD70 as a potentially valuable target for (P)VRL stratification. Finally, we also show that increased CCL17 levels might help to distinguish sarcoid uveitis from TB-associated uveitis.

Exploring Differentially Methylated Genes in Vulvar Squamous Cell Carcinoma.

DNA methylation is the most widely studied mechanism of epigenetic modification, which can influence gene expression without alterations in DNA sequences. Aberrations in DNA methylation are known to play a role in carcinogenesis, and methylation profiling has enabled the identification of biomarkers of potential clinical interest for several cancers. For vulvar squamous cell carcinoma (VSCC), however, methylation profiling remains an under-studied area. We sought to identify differentially methylated genes (DMGs) in VSCC, by performing Infinium MethylationEPIC BeadChip (Illumina) array sequencing, on a set of primary VSCC (n = 18), and normal vulvar tissue from women with no history of vulvar (pre)malignancies (n = 6). Using a false-discovery rate of 0.05, beta-difference (Δß) of ±0.5, and CpG-island probes as cut-offs, 199 DMGs (195 hyper-methylated, 4 hypo-methylated) were identified for VSCC. Most of the hyper-methylated genes were found to be involved in transcription regulator activity, indicating that disruption of this process plays a vital role in VSCC development. The majority of VSCCs harbored amplifications of chromosomes 3, 8, and 9. We identified a set of DMGs in this exploratory, hypothesis-generating study, which we hope will facilitate epigenetic profiling of VSCCs. Prognostic relevance of these DMGs deserves further exploration in larger cohorts of VSCC and its precursor lesions.

Nuclear factor IB is downregulated in vulvar squamous cell carcinoma (VSCC): Unravelling differentially expressed genes in VSCC through gene expression dataset analysis.

Vulvar squamous cell carcinoma (VSCC) comprises two distinct etiopathological subtypes: i) Human papilloma virus (HPV)-related VSCC, which arises via the precursor high grade squamous intraepithelial lesion (HSIL); and ii) HPV-independent VSCC, which arises via precursor, differentiated vulvar intraepithelial neoplasia (dVIN), driven by TP53 mutations. However, the mechanism of carcinogenesis of VSCC is poorly understood. The current study aimed to gain insight into VSCC carcinogenesis by identifying differentially expressed genes (DEGs) for each VSCC subtype. The expression of certain DEGs was then further assessed by performing immunohistochemistry (IHC) on whole tissue sections of VSCC and its precursors. Statistical analysis of microarrays was performed on two independent gene expression datasets (GSE38228 and a study from Erasmus MC) on VSCC and normal vulva. DEGs were identified that were similarly (up/down) regulated with statistical significance in both datasets. For HPV-related VSCCs, this constituted 88 DEGs, and for HPV-independent VSCCs, this comprised 46 DEGs. IHC was performed on VSCC (n=11), dVIN (n=6), HSIL (n=6) and normal vulvar tissue (n=7) with i) signal transducer and activator of transcription 1 (STAT1; an upregulated DEGs); ii) nuclear factor IB (NFIB; a downregulated DEG); iii) p16 (to determine the HPV status of tissues); and iv) p53 (to confirm the histological diagnoses). Strong and diffuse NFIB expression was observed in the basal and para-basal layers of normal vulvar tissue, whereas NFIB expression was minimal or completely negative in dVIN and in both subtypes of VSCC. In contrast, no discernable difference was observed in STAT1 expression among normal vulvar tissue, dVIN, HSIL or VSCC. By leveraging bioinformatics, the current study identified DEGs that can facilitate research into VSCC carcinogenesis. The results suggested that NFIB is downregulated in VSCC and its relevance as a diagnostic/prognostic biomarker deserves further exploration.

Evaluation of Immunohistochemical Markers, CK17 and SOX2, as Adjuncts to p53 for the Diagnosis of Differentiated Vulvar Intraepithelial Neoplasia (dVIN).

Histological diagnosis of differentiated vulvar intraepithelial neoplasia (dVIN), the precursor of human papillomavirus (HPV)-independent vulvar squamous cell carcinoma (VSCC), can be challenging, as features of dVIN may mimic those of non-dysplastic dermatoses. To aid the diagnosis, p53-immunohistochemistry (IHC) is commonly used, and mutant expression patterns are used to support a histological diagnosis of dVIN. However, a proportion of dVIN can show wild-type p53-expression, which is characteristic of non-dysplastic dermatoses. Furthermore, recent research has identified a novel precursor of HPV-independent VSCC-the p53-wild-type differentiated exophytic vulvar intraepithelial lesion (de-VIL). Currently, there are no established diagnostic IHC-markers for p53-wild-type dVIN or de-VIL. We evaluated IHC-markers, cytokeratin 17 (CK17), and SRY-box 2 (SOX2), as diagnostic adjuncts for dVIN. For this, IHC-expression of CK17, SOX2, and p53 was studied in dVIN (n = 56), de-VIL (n = 8), and non-dysplastic vulvar tissues (n = 46). For CK17 and SOX2, the percentage of cells showing expression, and the intensity and distribution of expression were recorded. We also performed next generation targeted sequencing (NGTS) on a subset of dVIN (n = 8) and de-VIL (n = 8). With p53-IHC, 74% of dVIN showed mutant patterns and 26% showed wild-type expression. Median percentage of cells expressing CK17 or SOX2 was significantly higher in dVIN (p53-mutant or p53-wild-type) and de-VIL than in non-dysplastic tissues (p < 0.01). Diffuse, moderate-to-strong, full epithelial expression of CK17 or SOX2 was highly specific for dVIN and de-VIL. With NGTS, TP53 mutations were detected in both dVIN and de-VIL. We infer that immunohistochemical markers CK17 and SOX2, when used along with p53, may help support the histological diagnosis of dVIN.

Protein instability associated with AARS1 and MARS1 mutations causes trichothiodystrophy.

Trichothiodystrophy (TTD) is a rare hereditary neurodevelopmental disorder defined by sulfur-deficient brittle hair and nails and scaly skin, but with otherwise remarkably variable clinical features. The photosensitive TTD (PS-TTD) forms exhibits in addition to progressive neuropathy and other features of segmental accelerated aging and is associated with impaired genome maintenance and transcription. New factors involved in various steps of gene expression have been identified for the different non-photosensitive forms of TTD (NPS-TTD), which do not appear to show features of premature aging. Here, we identify alanyl-tRNA synthetase 1 and methionyl-tRNA synthetase 1 variants as new gene defects that cause NPS-TTD. These variants result in the instability of the respective gene products alanyl- and methionyl-tRNA synthetase. These findings extend our previous observations that TTD mutations affect the stability of the corresponding proteins and emphasize this phenomenon as a common feature of TTD. Functional studies in skin fibroblasts from affected individuals demonstrate that these new variants also impact on the rate of tRNA charging, which is the first step in protein translation. The extension of reduced abundance of TTD factors to translation as well as transcription redefines TTD as a syndrome in which proteins involved in gene expression are unstable.

Sex difference in the mouse BAT transcriptome reveals a role of progesterone.

Brown adipose tissue (BAT) is a metabolically active organ that exhibits sex-differential features, that is, being generally more abundant and active in females than in males. Although sex steroids, particularly estrogens, have been shown to regulate BAT thermogenic function, the underlying molecular mechanisms contributing to sexual dimorphism in basal BAT activity have not been elucidated. Therefore, we assessed the transcriptome of interscapular BAT of male and female C57BL/6J mice by RNA sequencing and identified 295 genes showing ≥2-fold differential expression (adjusted P < 0.05). In silico functional annotation clustering suggested an enrichment of genes encoding proteins involved in cell-cell contact, interaction, and adhesion. Ovariectomy reduced the expression of these genes in female BAT toward a male pattern whereas orchiectomy had marginal effects on the transcriptional pattern, indicating a prominent role of female gonadal hormones in this sex-differential expression pattern. Progesterone was identified as a possible upstream regulator of the sex-differentially expressed genes. Studying the direct effects of progesterone in vitro in primary adipocytes showed that progesterone significantly altered the transcription of several of the identified genes, possibly via the glucocorticoid receptor. In conclusion, this study reveals a sexually dimorphic transcription profile in murine BAT at general housing conditions and demonstrates a role for progesterone in the regulation of the interscapular BAT transcriptome.

Correction: SMAD6 variants in craniosynostosis: genotype and phenotype evaluation.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

SMAD6 variants in craniosynostosis: genotype and phenotype evaluation.

PURPOSE: Enrichment of heterozygous missense and truncating SMAD6 variants was previously reported in nonsyndromic sagittal and metopic synostosis, and interaction of SMAD6 variants with a common polymorphism nearBMP2 (rs1884302) was proposed to contribute to inconsistent penetrance. We determined the occurrence of SMAD6 variants in all types of craniosynostosis, evaluated the impact of different missense variants on SMAD6 function, and tested independently whether rs1884302 genotype significantly modifies the phenotype. METHODS: We performed resequencing of SMAD6 in 795 unsolved patients with any type of craniosynostosis and genotyped rs1884302 in SMAD6-positive individuals and relatives. We examined the inhibitory activity and stability of SMAD6 missense variants. RESULTS: We found 18 (2.3%) different rare damaging SMAD6 variants, with the highest prevalence in metopic synostosis (5.8%) and an 18.3-fold enrichment of loss-of-function variants comparedwith gnomAD data (P < 10-7). Combined with eight additional variants, ≥20/26 were transmitted from an unaffected parent but rs1884302 genotype did not predict phenotype. CONCLUSION: Pathogenic SMAD6 variants substantially increase the risk of both nonsyndromic and syndromic presentations of craniosynostosis, especially metopic synostosis. Functional analysis is important to evaluate missense variants. Genotyping of rs1884302 is not clinically useful. Mechanisms to explain the remarkable diversity of phenotypes associated with SMAD6 variants remain obscure.