A novel SMARCC1 BAFopathy implicates neural progenitor epigenetic dysregulation in human hydrocephalus

Hydrocephalus, characterized by cerebral ventriculomegaly, is the most common disorder requiring brain surgery in children. Recent studies have implicated SMARCC1, a component of the BRG1-associated factor (BAF) chromatin remodelling complex, as a candidate congenital hydrocephalus gene. However, SMARCC1 variants have not been systematically examined in a large patient cohort or conclusively linked with a human syndrome. Moreover, congenital hydrocephalus-associated SMARCC1 variants have not been functionally validated or mechanistically studied in vivo. Here, we aimed to assess the prevalence of SMARCC1 variants in an expanded patient cohort, describe associated clinical and radiographic phenotypes, and assess the impact of Smarcc1 depletion in a novel Xenopus tropicalis model of congenital hydrocephalus. To do this, we performed a genetic association study using whole-exome sequencing from a cohort consisting of 2697 total ventriculomegalic trios, including patients with neurosurgically-treated congenital hydrocephalus, that total 8091 exomes collected over 7 years (2016-23). A comparison control cohort consisted of 1798 exomes from unaffected siblings of patients with autism spectrum disorder and their unaffected parents were sourced from the Simons Simplex Collection. Enrichment and impact on protein structure were assessed in identified variants. Effects on the human fetal brain transcriptome were examined with RNA-sequencing and Smarcc1 knockdowns were generated in Xenopus and studied using optical coherence tomography imaging, in situ hybridization and immunofluorescence. SMARCC1 surpassed genome-wide significance thresholds, yielding six rare, protein-altering de novo variants localized to highly conserved residues in key functional domains. Patients exhibited hydrocephalus with aqueductal stenosis; corpus callosum abnormalities, developmental delay, and cardiac defects were also common. Xenopus knockdowns recapitulated both aqueductal stenosis and cardiac defects and were rescued by wild-type but not patient-specific variant SMARCC1. Hydrocephalic SMARCC1-variant human fetal brain and Smarcc1-variant Xenopus brain exhibited a similarly altered expression of key genes linked to midgestational neurogenesis, including the transcription factors NEUROD2 and MAB21L2. These results suggest de novo variants in SMARCC1 cause a novel human BAFopathy we term 'SMARCC1-associated developmental dysgenesis syndrome', characterized by variable presence of cerebral ventriculomegaly, aqueductal stenosis, developmental delay and a variety of structural brain or cardiac defects. These data underscore the importance of SMARCC1 and the BAF chromatin remodelling complex for human brain morphogenesis and provide evidence for a 'neural stem cell' paradigm of congenital hydrocephalus pathogenesis. These results highlight utility of trio-based whole-exome sequencing for identifying pathogenic variants in sporadic congenital structural brain disorders and suggest whole-exome sequencing may be a valuable adjunct in clinical management of congenital hydrocephalus patients.

GIMAP5 deficiency reveals a mammalian ceramide-driven longevity assurance pathway

Preserving cells in a functional, non-senescent state is a major goal for extending human healthspans. Model organisms reveal that longevity and senescence are genetically controlled, but how genes control longevity in different mammalian tissues is unknown. Here, we report a new human genetic disease that causes cell senescence, liver and immune dysfunction, and early mortality that results from deficiency of GIMAP5, an evolutionarily conserved GTPase selectively expressed in lymphocytes and endothelial cells. We show that GIMAP5 restricts the pathological accumulation of long-chain ceramides (CERs), thereby regulating longevity. GIMAP5 controls CER abundance by interacting with protein kinase CK2 (CK2), attenuating its ability to activate CER synthases. Inhibition of CK2 and CER synthase rescues GIMAP5-deficient T cells by preventing CER overaccumulation and cell deterioration. Thus, GIMAP5 controls longevity assurance pathways crucial for immune function and healthspan in mammals.

LRRC23 truncation impairs radial spoke 3 head assembly and sperm motility underlying male infertility

Radial spokes (RS) are T-shaped multiprotein complexes on the axonemal microtubules. Repeated RS1, RS2, and RS3 couple the central pair to modulate ciliary and flagellar motility. Despite the cell type specificity of RS3 substructures, their molecular components remain largely unknown. Here, we report that a leucine-rich repeat-containing protein, LRRC23, is an RS3 head component essential for its head assembly and flagellar motility in mammalian spermatozoa. From infertile male patients with defective sperm motility, we identified a splice site variant of LRRC23. A mutant mouse model mimicking this variant produces a truncated LRRC23 at the C-terminus that fails to localize to the sperm tail, causing male infertility due to defective sperm motility. LRRC23 was previously proposed to be an ortholog of the RS stalk protein RSP15. However, we found that purified recombinant LRRC23 interacts with an RS head protein RSPH9, which is abolished by the C-terminal truncation. Evolutionary and structural comparison also shows that LRRC34, not LRRC23, is the RSP15 ortholog. Cryo-electron tomography clearly revealed that the absence of the RS3 head and the sperm-specific RS2-RS3 bridge structure in LRRC23 mutant spermatozoa. Our study provides new insights into the structure and function of RS3 in mammalian spermatozoa and the molecular pathogenicity of LRRC23 underlying reduced sperm motility in infertile human males.

Mutation of key signaling regulators of cerebrovascular development in vein of Galen malformations

To elucidate the pathogenesis of vein of Galen malformations (VOGMs), the most common and most severe of congenital brain arteriovenous malformations, we performed an integrated analysis of 310 VOGM proband-family exomes and 336,326 human cerebrovasculature single-cell transcriptomes. We found the Ras suppressor p120 RasGAP (RASA1) harbored a genome-wide significant burden of loss-of-function de novo variants (2042.5-fold, p = 4.79 x 10-7). Rare, damaging transmitted variants were enriched in Ephrin receptor-B4 (EPHB4) (17.5-fold, p = 1.22 x 10-5), which cooperates with p120 RasGAP to regulate vascular development. Additional probands had damaging variants in ACVRL1, NOTCH1, ITGB1, and PTPN11. ACVRL1 variants were also identified in a multi-generational VOGM pedigree. Integrative genomic analysis defined developing endothelial cells as a likely spatio-temporal locus of VOGM pathophysiology. Mice expressing a VOGM-specific EPHB4 kinase-domain missense variant (Phe867Leu) exhibited disrupted developmental angiogenesis and impaired hierarchical development of arterial-capillary-venous networks, but only in the presence of a "second-hit" allele. These results illuminate human arterio-venous development and VOGM pathobiology and have implications for patients and their families.

Isradipine therapy in Cacna1dIle772Met/+ mice ameliorates primary aldosteronism and neurologic abnormalities

Somatic gain-of-function mutations in the L-type calcium channel CaV1.3 (CACNA1D gene) cause adrenal aldosterone-producing adenomas and micronodules. De novo germline mutations are found in a syndrome of primary aldosteronism, seizures, and neurologic abnormalities (PASNA) as well as in autism spectrum disorder. Using CRISPR/Cas9, we here generated mice with a Cacna1d gain-of-function mutation found in both adenomas and PASNA syndrome (Cacna1dIle772Met/+). These mice show reduced body weight and increased mortality from weaning to approximately 100 days of age. Male mice do not breed, likely due to neuromotor impairment, and the offspring of female mice die perinatally, likely due to lack of maternal care. Mice generated by in vitro fertilization showed elevated intracellular calcium in the aldosterone-producing zona glomerulosa, an elevated aldosterone/renin ratio, and persistently elevated serum aldosterone on a high-salt diet as signs of primary aldosteronism. Anesthesia with ketamine and xylazine induced tonic-clonic seizures. Neurologic abnormalities included hyperlocomotion, impaired performance in the rotarod test, impaired nest building, and slight changes in social behavior. Intracellular calcium in the zona glomerulosa, aldosterone levels, and rotarod performance responded to treatment with the calcium channel blocker isradipine, with implications for the therapy of patients with aldosterone-producing lesions and with PASNA syndrome.

LRRC23 truncation impairs radial spoke 3 head assembly and sperm motility underlying male infertility

Radial spokes (RS) are T-shaped multiprotein complexes on the axonemal microtubules. Repeated RS1, RS2, and RS3 couple the central pair to modulate ciliary and flagellar motility. Despite the cell type specificity of RS3 substructures, their molecular components remain largely unknown. Here, we report that a leucine-rich repeat-containing protein, LRRC23, is an RS3 head component essential for its head assembly and flagellar motility in mammalian spermatozoa. From infertile male patients with defective sperm motility, we identified a splice site variant of LRRC23. A mutant mouse model mimicking this variant produces a truncated LRRC23 at the C-terminus that fails to localize to the sperm tail, causing male infertility due to defective sperm motility. LRRC23 was previously proposed to be an ortholog of the RS stalk protein RSP15. However, we found that purified recombinant LRRC23 interacts with an RS head protein RSPH9, which is abolished by the C-terminal truncation. Evolutionary and structural comparison also shows that LRRC34, not LRRC23, is the RSP15 ortholog. Cryo-electron tomography clearly revealed that the absence of the RS3 head and the sperm-specific RS2-RS3 bridge structure in LRRC23 mutant spermatozoa. Our study provides new insights into the structure and function of RS3 in mammalian spermatozoa and the molecular pathogenicity of LRRC23 underlying reduced sperm motility in infertile human males.

Genome-wide association analyses define pathogenic signaling pathways and prioritize drug targets for IgA nephropathy

IgA nephropathy (IgAN) is a progressive form of kidney disease defined by glomerular deposition of IgA. Here we performed a genome-wide association study of 10,146 kidney-biopsy-diagnosed IgAN cases and 28,751 controls across 17 international cohorts. We defined 30 genome-wide significant risk loci explaining 11% of disease risk. A total of 16 loci were new, including TNFSF4/TNFSF18, REL, CD28, PF4V1, LY86, LYN, ANXA3, TNFSF8/TNFSF15, REEP3, ZMIZ1, OVOL1/RELA, ETS1, IGH, IRF8, TNFRSF13B and FCAR. The risk loci were enriched in gene orthologs causing abnormal IgA levels when genetically manipulated in mice. We also observed a positive genetic correlation between IgAN and serum IgA levels. High polygenic score for IgAN was associated with earlier onset of kidney failure. In a comprehensive functional annotation analysis of candidate causal genes, we observed convergence of biological candidates on a common set of inflammatory signaling pathways and cytokine ligand-receptor pairs, prioritizing potential new drug targets.

Contribution of Somatic Ras/Raf/Mitogen-Activated Protein Kinase Variants in the Hippocampus in Drug-Resistant Mesial Temporal Lobe Epilepsy

Importance

Mesial temporal lobe epilepsy (MTLE) is the most common focal epilepsy subtype and is often refractory to antiseizure medications. While most patients with MTLE do not have pathogenic germline genetic variants, the contribution of postzygotic (ie, somatic) variants in the brain is unknown.

Objective

To test the association between pathogenic somatic variants in the hippocampus and MTLE.

Design, Setting, and Participants

This case-control genetic association study analyzed the DNA derived from hippocampal tissue of neurosurgically treated patients with MTLE and age-matched and sex-matched neurotypical controls. Participants treated at level 4 epilepsy centers were enrolled from 1988 through 2019, and clinical data were collected retrospectively. Whole-exome and gene-panel sequencing (each genomic region sequenced more than 500 times on average) were used to identify candidate pathogenic somatic variants. A subset of novel variants was functionally evaluated using cellular and molecular assays. Patients with nonlesional and lesional (mesial temporal sclerosis, focal cortical dysplasia, and low-grade epilepsy-associated tumors) drug-resistant MTLE who underwent anterior medial temporal lobectomy were eligible. All patients with available frozen tissue and appropriate consents were included. Control brain tissue was obtained from neurotypical donors at brain banks. Data were analyzed from June 2020 to August 2022.

Exposures

Drug-resistant MTLE.

Main Outcomes and Measures

Presence and abundance of pathogenic somatic variants in the hippocampus vs the unaffected temporal neocortex.

Results

Of 105 included patients with MTLE, 53 (50.5%) were female, and the median (IQR) age was 32 (26-44) years; of 30 neurotypical controls, 11 (36.7%) were female, and the median (IQR) age was 37 (18-53) years. Eleven pathogenic somatic variants enriched in the hippocampus relative to the unaffected temporal neocortex (median [IQR] variant allele frequency, 1.92 [1.5-2.7] vs 0.3 [0-0.9]; P = .01) were detected in patients with MTLE but not in controls. Ten of these variants were in PTPN11, SOS1, KRAS, BRAF, and NF1, all predicted to constitutively activate Ras/Raf/mitogen-activated protein kinase (MAPK) signaling. Immunohistochemical studies of variant-positive hippocampal tissue demonstrated increased Erk1/2 phosphorylation, indicative of Ras/Raf/MAPK activation, predominantly in glial cells. Molecular assays showed abnormal liquid-liquid phase separation for the PTPN11 variants as a possible dominant gain-of-function mechanism.

Conclusions and Relevance

Hippocampal somatic variants, particularly those activating Ras/Raf/MAPK signaling, may contribute to the pathogenesis of sporadic, drug-resistant MTLE. These findings may provide a novel genetic mechanism and highlight new therapeutic targets for this common indication for epilepsy surgery.

Rare Single Nucleotide and Copy Number Variants and the Etiology of Congenital Obstructive Uropathy: Implications for Genetic Diagnosis

SIGNIFICANCE STATEMENT

Congenital obstructive uropathy (COU) is a prevalent human developmental defect with highly heterogeneous clinical presentations and outcomes. Genetics may refine diagnosis, prognosis, and treatment, but the genomic architecture of COU is largely unknown. Comprehensive genomic screening study of 733 cases with three distinct COU subphenotypes revealed disease etiology in 10.0% of them. We detected no significant differences in the overall diagnostic yield among COU subphenotypes, with characteristic variable expressivity of several mutant genes. Our findings therefore may legitimize a genetic first diagnostic approach for COU, especially when burdening clinical and imaging characterization is not complete or available.

BACKGROUND

Congenital obstructive uropathy (COU) is a common cause of developmental defects of the urinary tract, with heterogeneous clinical presentation and outcome. Genetic analysis has the potential to elucidate the underlying diagnosis and help risk stratification.

METHODS

We performed a comprehensive genomic screen of 733 independent COU cases, which consisted of individuals with ureteropelvic junction obstruction ( n =321), ureterovesical junction obstruction/congenital megaureter ( n =178), and COU not otherwise specified (COU-NOS; n =234).

RESULTS

We identified pathogenic single nucleotide variants (SNVs) in 53 (7.2%) cases and genomic disorders (GDs) in 23 (3.1%) cases. We detected no significant differences in the overall diagnostic yield between COU sub-phenotypes, and pathogenic SNVs in several genes were associated to any of the three categories. Hence, although COU may appear phenotypically heterogeneous, COU phenotypes are likely to share common molecular bases. On the other hand, mutations in TNXB were more often identified in COU-NOS cases, demonstrating the diagnostic challenge in discriminating COU from hydronephrosis secondary to vesicoureteral reflux, particularly when diagnostic imaging is incomplete. Pathogenic SNVs in only six genes were found in more than one individual, supporting high genetic heterogeneity. Finally, convergence between data on SNVs and GDs suggest MYH11 as a dosage-sensitive gene possibly correlating with severity of COU.

CONCLUSIONS

We established a genomic diagnosis in 10.0% of COU individuals. The findings underscore the urgent need to identify novel genetic susceptibility factors to COU to better define the natural history of the remaining 90% of cases without a molecular diagnosis.

De novo variants implicate chromatin modification, transcriptional regulation, and retinoic acid signaling in syndromic craniosynostosis

Craniosynostosis (CS) is the most common congenital cranial anomaly. Several Mendelian forms of syndromic CS are well described, but a genetic etiology remains elusive in a substantial fraction of probands. Analysis of exome sequence data from 526 proband-parent trios with syndromic CS identified a marked excess (observed 98, expected 33, p = 4.83 × 10-20) of damaging de novo variants (DNVs) in genes highly intolerant to loss-of-function variation (probability of LoF intolerance > 0.9). 30 probands harbored damaging DNVs in 21 genes that were not previously implicated in CS but are involved in chromatin modification and remodeling (4.7-fold enrichment, p = 1.1 × 10-11). 17 genes had multiple damaging DNVs, and 13 genes (CDK13, NFIX, ADNP, KMT5B, SON, ARID1B, CASK, CHD7, MED13L, PSMD12, POLR2A, CHD3, and SETBP1) surpassed thresholds for genome-wide significance. A recurrent gain-of-function DNV in the retinoic acid receptor alpha (RARA; c.865G>A [p.Gly289Arg]) was identified in two probands with similar CS phenotypes. CS risk genes overlap with those identified for autism and other neurodevelopmental disorders, are highly expressed in cranial neural crest cells, and converge in networks that regulate chromatin modification, gene transcription, and osteoblast differentiation. Our results identify several CS loci and have major implications for genetic testing and counseling.

Pleiotropic role of TRAF7 in skull-base meningiomas and congenital heart disease

While somatic variants of TRAF7 (Tumor necrosis factor receptor-associated factor 7) underlie anterior skull-base meningiomas, here we report the inherited mutations of TRAF7 that cause congenital heart defects. We show that TRAF7 mutants operate in a dominant manner, inhibiting protein function via heterodimerization with wild-type protein. Further, the shared genetics of the two disparate pathologies can be traced to the common origin of forebrain meninges and cardiac outflow tract from the TRAF7-expressing neural crest. Somatic and inherited mutations disrupt TRAF7-IFT57 interactions leading to cilia degradation. TRAF7-mutant meningioma primary cultures lack cilia, and TRAF7 knockdown causes cardiac, craniofacial, and ciliary defects in Xenopus and zebrafish, suggesting a mechanistic convergence for TRAF7-driven meningiomas and developmental heart defects.

Rare predicted loss-of-function variants of type I IFN immunity genes are associated with life-threatening COVID-19

BACKGROUND

We previously reported that impaired type I IFN activity, due to inborn errors of TLR3- and TLR7-dependent type I interferon (IFN) immunity or to autoantibodies against type I IFN, account for 15-20% of cases of life-threatening COVID-19 in unvaccinated patients. Therefore, the determinants of life-threatening COVID-19 remain to be identified in ~ 80% of cases.

METHODS

We report here a genome-wide rare variant burden association analysis in 3269 unvaccinated patients with life-threatening COVID-19, and 1373 unvaccinated SARS-CoV-2-infected individuals without pneumonia. Among the 928 patients tested for autoantibodies against type I IFN, a quarter (234) were positive and were excluded.

RESULTS

No gene reached genome-wide significance. Under a recessive model, the most significant gene with at-risk variants was TLR7, with an OR of 27.68 (95%CI 1.5-528.7, P = 1.1 × 10-4) for biochemically loss-of-function (bLOF) variants. We replicated the enrichment in rare predicted LOF (pLOF) variants at 13 influenza susceptibility loci involved in TLR3-dependent type I IFN immunity (OR = 3.70[95%CI 1.3-8.2], P = 2.1 × 10-4). This enrichment was further strengthened by (1) adding the recently reported TYK2 and TLR7 COVID-19 loci, particularly under a recessive model (OR = 19.65[95%CI 2.1-2635.4], P = 3.4 × 10-3), and (2) considering as pLOF branchpoint variants with potentially strong impacts on splicing among the 15 loci (OR = 4.40[9%CI 2.3-8.4], P = 7.7 × 10-8). Finally, the patients with pLOF/bLOF variants at these 15 loci were significantly younger (mean age [SD] = 43.3 [20.3] years) than the other patients (56.0 [17.3] years; P = 1.68 × 10-5).

CONCLUSIONS

Rare variants of TLR3- and TLR7-dependent type I IFN immunity genes can underlie life-threatening COVID-19, particularly with recessive inheritance, in patients under 60 years old.

Genetic dysregulation of an endothelial Ras signaling network in vein of Galen malformations

To elucidate the pathogenesis of vein of Galen malformations (VOGMs), the most common and severe congenital brain arteriovenous malformation, we performed an integrated analysis of 310 VOGM proband-family exomes and 336,326 human cerebrovasculature single-cell transcriptomes. We found the Ras suppressor p120 RasGAP ( RASA1 ) harbored a genome-wide significant burden of loss-of-function de novo variants (p=4.79×10 ^-7 ). Rare, damaging transmitted variants were enriched in Ephrin receptor-B4 ( EPHB4 ) (p=1.22×10 ^-5 ), which cooperates with p120 RasGAP to limit Ras activation. Other probands had pathogenic variants in ACVRL1 , NOTCH1 , ITGB1 , and PTPN11 . ACVRL1 variants were also identified in a multi-generational VOGM pedigree. Integrative genomics defined developing endothelial cells as a key spatio-temporal locus of VOGM pathophysiology. Mice expressing a VOGM-specific EPHB4 kinase-domain missense variant exhibited constitutive endothelial Ras/ERK/MAPK activation and impaired hierarchical development of angiogenesis-regulated arterial-capillary-venous networks, but only when carrying a "second-hit" allele. These results illuminate human arterio-venous development and VOGM pathobiology and have clinical implications.

A novel SMARCC1 -mutant BAFopathy implicates epigenetic dysregulation of neural progenitors in hydrocephalus

Importance

Hydrocephalus, characterized by cerebral ventriculomegaly, is the most common disorder requiring brain surgery. A few familial forms of congenital hydrocephalus (CH) have been identified, but the cause of most sporadic cases of CH remains elusive. Recent studies have implicated SMARCC1 , a component of the B RG1- a ssociated factor (BAF) chromatin remodeling complex, as a candidate CH gene. However, SMARCC1 variants have not been systematically examined in a large patient cohort or conclusively linked with a human syndrome. Moreover, CH-associated SMARCC1 variants have not been functionally validated or mechanistically studied in vivo .

Objectives

The aims of this study are to (i) assess the extent to which rare, damaging de novo mutations (DNMs) in SMARCC1 are associated with cerebral ventriculomegaly; (ii) describe the clinical and radiographic phenotypes of SMARCC1 -mutated patients; and (iii) assess the pathogenicity and mechanisms of CH-associated SMARCC1 mutations in vivo .

Design setting and participants

A genetic association study was conducted using whole-exome sequencing from a cohort consisting of 2,697 ventriculomegalic trios, including patients with neurosurgically-treated CH, totaling 8,091 exomes collected over 5 years (2016-2021). Data were analyzed in 2023. A comparison control cohort consisted of 1,798 exomes from unaffected siblings of patients with autism spectrum disorder and their unaffected parents sourced from the Simons simplex consortium.

Main outcomes and measures

Gene variants were identified and filtered using stringent, validated criteria. Enrichment tests assessed gene-level variant burden. In silico biophysical modeling estimated the likelihood and extent of the variant impact on protein structure. The effect of a CH-associated SMARCC1 mutation on the human fetal brain transcriptome was assessed by analyzing RNA-sequencing data. Smarcc1 knockdowns and a patient-specific Smarcc1 variant were tested in Xenopus and studied using optical coherence tomography imaging, in situ hybridization, and immunofluorescence microscopy.

Results

SMARCC1 surpassed genome-wide significance thresholds in DNM enrichment tests. Six rare protein-altering DNMs, including four loss-of-function mutations and one recurrent canonical splice site mutation (c.1571+1G>A) were detected in unrelated patients. DNMs localized to the highly conserved DNA-interacting SWIRM, Myb-DNA binding, Glu-rich, and Chromo domains of SMARCC1 . Patients exhibited developmental delay (DD), aqueductal stenosis, and other structural brain and heart defects. G0 and G1 Smarcc1 Xenopus mutants exhibited aqueductal stenosis and cardiac defects and were rescued by human wild-type SMARCC1 but not a patient-specific SMARCC1 mutant. Hydrocephalic SMARCC1 -mutant human fetal brain and Smarcc1 -mutant Xenopus brain exhibited a similarly altered expression of key genes linked to midgestational neurogenesis, including the transcription factors NEUROD2 and MAB21L2 .

Conclusions

SMARCC1 is a bona fide CH risk gene. DNMs in SMARCC1 cause a novel human BAFopathy we term " S MARCC1- a ssociated D evelopmental D ysgenesis S yndrome (SaDDS)", characterized by cerebral ventriculomegaly, aqueductal stenosis, DD, and a variety of structural brain or cardiac defects. These data underscore the importance of SMARCC1 and the BAF chromatin remodeling complex for human brain morphogenesis and provide evidence for a "neural stem cell" paradigm of human CH pathogenesis. These results highlight the utility of trio-based WES for identifying risk genes for congenital structural brain disorders and suggest WES may be a valuable adjunct in the clinical management of CH patients.

KEY POINTS

Question: What is the role of SMARCC1 , a core component of the B RG1- a ssociated factor (BAF) chromatin remodeling complex, in brain morphogenesis and congenital hydrocephalus (CH)? Findings: SMARCC1 harbored an exome-wide significant burden of rare, protein-damaging de novo mutations (DNMs) (p = 5.83 × 10 -9 ) in the largest ascertained cohort to date of patients with cerebral ventriculomegaly, including treated CH (2,697 parent-proband trios). SMARCC1 contained four loss-of-function DNMs and two identical canonical splice site DNMs in a total of six unrelated patients. Patients exhibited developmental delay, aqueductal stenosis, and other structural brain and cardiac defects. Xenopus Smarcc1 mutants recapitulated core human phenotypes and were rescued by the expression of human wild-type but not patient-mutant SMARCC1 . Hydrocephalic SMARCC1 -mutant human brain and Smarcc1 -mutant Xenopus brain exhibited similar alterationsin the expression of key transcription factors that regulate neural progenitor cell proliferation. Meaning: SMARCC1 is essential for human brain morphogenesis and is a bona fide CH risk gene. SMARCC1 mutations cause a novel human BAFopathy we term " S MARCC1- a ssociated D evelopmental D ysgenesis S yndrome (SaDDS)". These data implicate epigenetic dysregulation of fetal neural progenitors in the pathogenesis of hydrocephalus, with diagnostic and prognostic implications for patients and caregivers.

Molecular genetics of human developmental neurocranial anomalies: towards "precision surgery"

Recent trio-based whole-exome sequencing studies of congenital hydrocephalus and nonsyndromic craniosynostosis have identified multiple novel disease genes that have illuminated the pathogenesis of these disorders and shed new insight into the genetic regulation of human brain and skull development. Continued study of these and other historically understudied developmental anomalies has the potential to replace the current antiquated, anatomically based disease classification systems with a molecular nomenclature that may increase precision for genetic counseling, prognostication, and surgical treatment stratification-including when not to operate. Data will also inform future clinical trials, catalyze the development of targeted therapies, and generate infrastructure and publicly available data sets relevant for other related nonsurgical neurodevelopmental and neuropsychiatric diseases.

Multiomic analyses implicate a neurodevelopmental program in the pathogenesis of cerebral arachnoid cysts

Cerebral arachnoid cysts (ACs) are one of the most common and poorly understood types of developmental brain lesion. To begin to elucidate AC pathogenesis, we performed an integrated analysis of 617 patient-parent (trio) exomes, 152,898 human brain and mouse meningeal single-cell RNA sequencing transcriptomes and natural language processing data of patient medical records. We found that damaging de novo variants (DNVs) were highly enriched in patients with ACs compared with healthy individuals (P = 1.57 × 10^-33). Seven genes harbored an exome-wide significant DNV burden. AC-associated genes were enriched for chromatin modifiers and converged in midgestational transcription networks essential for neural and meningeal development. Unsupervised clustering of patient phenotypes identified four AC subtypes and clinical severity correlated with the presence of a damaging DNV. These data provide insights into the coordinated regulation of brain and meningeal development and implicate epigenomic dysregulation due to DNVs in AC pathogenesis. Our results provide a preliminary indication that, in the appropriate clinical context, ACs may be considered radiographic harbingers of neurodevelopmental pathology warranting genetic testing and neurobehavioral follow-up. These data highlight the utility of a systems-level, multiomics approach to elucidate sporadic structural brain disease.

The choroid plexus links innate immunity to CSF dysregulation in hydrocephalus

The choroid plexus (ChP) is the blood-cerebrospinal fluid (CSF) barrier and the primary source of CSF. Acquired hydrocephalus, caused by brain infection or hemorrhage, lacks drug treatments due to obscure pathobiology. Our integrated, multi-omic investigation of post-infectious hydrocephalus (PIH) and post-hemorrhagic hydrocephalus (PHH) models revealed that lipopolysaccharide and blood breakdown products trigger highly similar TLR4-dependent immune responses at the ChP-CSF interface. The resulting CSF "cytokine storm", elicited from peripherally derived and border-associated ChP macrophages, causes increased CSF production from ChP epithelial cells via phospho-activation of the TNF-receptor-associated kinase SPAK, which serves as a regulatory scaffold of a multi-ion transporter protein complex. Genetic or pharmacological immunomodulation prevents PIH and PHH by antagonizing SPAK-dependent CSF hypersecretion. These results reveal the ChP as a dynamic, cellularly heterogeneous tissue with highly regulated immune-secretory capacity, expand our understanding of ChP immune-epithelial cell cross talk, and reframe PIH and PHH as related neuroimmune disorders vulnerable to small molecule pharmacotherapy.

Inborn errors of OAS-RNase L in SARS-CoV-2-related multisystem inflammatory syndrome in children

Multisystem inflammatory syndrome in children (MIS-C) is a rare and severe condition that follows benign COVID-19. We report autosomal recessive deficiencies of OAS1, OAS2, or RNASEL in five unrelated children with MIS-C. The cytosolic double-stranded RNA (dsRNA)-sensing OAS1 and OAS2 generate 2'-5'-linked oligoadenylates (2-5A) that activate the single-stranded RNA-degrading ribonuclease L (RNase L). Monocytic cell lines and primary myeloid cells with OAS1, OAS2, or RNase L deficiencies produce excessive amounts of inflammatory cytokines upon dsRNA or severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) stimulation. Exogenous 2-5A suppresses cytokine production in OAS1-deficient but not RNase L-deficient cells. Cytokine production in RNase L-deficient cells is impaired by MDA5 or RIG-I deficiency and abolished by mitochondrial antiviral-signaling protein (MAVS) deficiency. Recessive OAS-RNase L deficiencies in these patients unleash the production of SARS-CoV-2-triggered, MAVS-mediated inflammatory cytokines by mononuclear phagocytes, thereby underlying MIS-C.

Genetic Variants in ARHGEF6 Cause Congenital Anomalies of the Kidneys and Urinary Tract in Humans, Mice, and Frogs

BACKGROUND

About 40 disease genes have been described to date for isolated CAKUT, the most common cause of childhood CKD. However, these genes account for only 20% of cases. ARHGEF6, a guanine nucleotide exchange factor that is implicated in biologic processes such as cell migration and focal adhesion, acts downstream of integrin-linked kinase (ILK) and parvin proteins. A genetic variant of ILK that causes murine renal agenesis abrogates the interaction of ILK with a murine focal adhesion protein encoded by Parva , leading to CAKUT in mice with this variant.

METHODS

To identify novel genes that, when mutated, result in CAKUT, we performed exome sequencing in an international cohort of 1265 families with CAKUT. We also assessed the effects in vitro of wild-type and mutant ARHGEF6 proteins, and the effects of Arhgef6 deficiency in mouse and frog models.

RESULTS

We detected six different hemizygous variants in the gene ARHGEF6 (which is located on the X chromosome in humans) in eight individuals from six families with CAKUT. In kidney cells, overexpression of wild-type ARHGEF6 -but not proband-derived mutant ARHGEF6 -increased active levels of CDC42/RAC1, induced lamellipodia formation, and stimulated PARVA-dependent cell spreading. ARHGEF6-mutant proteins showed loss of interaction with PARVA. Three-dimensional Madin-Darby canine kidney cell cultures expressing ARHGEF6-mutant proteins exhibited reduced lumen formation and polarity defects. Arhgef6 deficiency in mouse and frog models recapitulated features of human CAKUT.

CONCLUSIONS

Deleterious variants in ARHGEF6 may cause dysregulation of integrin-parvin-RAC1/CDC42 signaling, thereby leading to X-linked CAKUT.

De novo mutations in the BMP signaling pathway in lambdoid craniosynostosis

Lambdoid craniosynostosis (CS) is a congenital anomaly resulting from premature fusion of the cranial suture between the parietal and occipital bones. Predominantly sporadic, it is the rarest form of CS and its genetic etiology is largely unexplored. Exome sequencing of 25 kindreds, including 18 parent-offspring trios with sporadic lambdoid CS, revealed a marked excess of damaging (predominantly missense) de novo mutations that account for ~ 40% of sporadic cases. These mutations clustered in the BMP signaling cascade (P = 1.6 × 10^-7), including mutations in genes encoding BMP receptors (ACVRL1 and ACVR2A), transcription factors (SOX11, FOXO1) and a transcriptional co-repressor (IFRD1), none of which have been implicated in other forms of CS. These missense mutations are at residues critical for substrate or target sequence recognition and many are inferred to cause genetic gain-of-function. Additionally, mutations in transcription factor NFIX were implicated in syndromic craniosynostosis affecting diverse sutures. Single cell RNA sequencing analysis of the mouse lambdoid suture identified enrichment of mutations in osteoblast precursors (P = 1.6 × 10^-6), implicating perturbations in the balance between proliferation and differentiation of osteoprogenitor cells in lambdoid CS. The results contribute to the growing knowledge of the genetics of CS, have implications for genetic counseling, and further elucidate the molecular etiology of premature suture fusion.

Rare predicted loss-of-function variants of type I IFN immunity genes are associated with life-threatening COVID-19

Background

We previously reported inborn errors of TLR3- and TLR7-dependent type I interferon (IFN) immunity in 1-5% of unvaccinated patients with life-threatening COVID-19, and auto-antibodies against type I IFN in another 15-20% of cases.

Methods

We report here a genome-wide rare variant burden association analysis in 3,269 unvaccinated patients with life-threatening COVID-19 (1,301 previously reported and 1,968 new patients), and 1,373 unvaccinated SARS-CoV-2-infected individuals without pneumonia. A quarter of the patients tested had antibodies against type I IFN (234 of 928) and were excluded from the analysis.

Results

No gene reached genome-wide significance. Under a recessive model, the most significant gene with at-risk variants was TLR7 , with an OR of 27.68 (95%CI:1.5-528.7, P= 1.1×10 -4 ), in analyses restricted to biochemically loss-of-function (bLOF) variants. We replicated the enrichment in rare predicted LOF (pLOF) variants at 13 influenza susceptibility loci involved in TLR3-dependent type I IFN immunity (OR=3.70 [95%CI:1.3-8.2], P= 2.1×10 -4 ). Adding the recently reported TYK2 COVID-19 locus strengthened this enrichment, particularly under a recessive model (OR=19.65 [95%CI:2.1-2635.4]; P= 3.4×10 -3 ). When these 14 loci and TLR7 were considered, all individuals hemizygous ( n =20) or homozygous ( n =5) for pLOF or bLOF variants were patients (OR=39.19 [95%CI:5.2-5037.0], P =4.7×10 -7 ), who also showed an enrichment in heterozygous variants (OR=2.36 [95%CI:1.0-5.9], P =0.02). Finally, the patients with pLOF or bLOF variants at these 15 loci were significantly younger (mean age [SD]=43.3 [20.3] years) than the other patients (56.0 [17.3] years; P= 1.68×10 -5 ).

Conclusions

Rare variants of TLR3- and TLR7-dependent type I IFN immunity genes can underlie life-threatening COVID-19, particularly with recessive inheritance, in patients under 60 years old.

Recessive inborn errors of type I IFN immunity in children with COVID-19 pneumonia

Recessive or dominant inborn errors of type I interferon (IFN) immunity can underlie critical COVID-19 pneumonia in unvaccinated adults. The risk of COVID-19 pneumonia in unvaccinated children, which is much lower than in unvaccinated adults, remains unexplained. In an international cohort of 112 children (<16 yr old) hospitalized for COVID-19 pneumonia, we report 12 children (10.7%) aged 1.5-13 yr with critical (7 children), severe (3), and moderate (2) pneumonia and 4 of the 15 known clinically recessive and biochemically complete inborn errors of type I IFN immunity: X-linked recessive TLR7 deficiency (7 children) and autosomal recessive IFNAR1 (1), STAT2 (1), or TYK2 (3) deficiencies. Fibroblasts deficient for IFNAR1, STAT2, or TYK2 are highly vulnerable to SARS-CoV-2. These 15 deficiencies were not found in 1,224 children and adults with benign SARS-CoV-2 infection without pneumonia (P = 1.2 × 10-11) and with overlapping age, sex, consanguinity, and ethnicity characteristics. Recessive complete deficiencies of type I IFN immunity may underlie ∼10% of hospitalizations for COVID-19 pneumonia in children.

Quantifying concordant genetic effects of de novo mutations on multiple disorders

Exome sequencing on tens of thousands of parent-proband trios has identified numerous deleterious de novo mutations (DNMs) and implicated risk genes for many disorders. Recent studies have suggested shared genes and pathways are enriched for DNMs across multiple disorders. However, existing analytic strategies only focus on genes that reach statistical significance for multiple disorders and require large trio samples in each study. As a result, these methods are not able to characterize the full landscape of genetic sharing due to polygenicity and incomplete penetrance. In this work, we introduce EncoreDNM, a novel statistical framework to quantify shared genetic effects between two disorders characterized by concordant enrichment of DNMs in the exome. EncoreDNM makes use of exome-wide, summary-level DNM data, including genes that do not reach statistical significance in single-disorder analysis, to evaluate the overall and annotation-partitioned genetic sharing between two disorders. Applying EncoreDNM to DNM data of nine disorders, we identified abundant pairwise enrichment correlations, especially in genes intolerant to pathogenic mutations and genes highly expressed in fetal tissues. These results suggest that EncoreDNM improves current analytic approaches and may have broad applications in DNM studies.

Mutation spectrum of congenital heart disease in a consanguineous Turkish population

BACKGROUNDS

While many studies agree that consanguinity increases the rate of congenital heart disease (CHD), few genome analyses have been conducted with consanguineous CHD cohorts.

METHODS

We recruited 73 CHD probands from consanguineous families in Turkey and used whole-exome sequencing (WES) to identify genetic lesions in these patients.

RESULTS

On average, each patient had 6.95 rare damaging homozygous variants, 0.68 of which are loss-of-function (LoF) variants. Seven patients (9.6%) carried damaging homozygous variants in five causal CHD genes. Six of those patients exhibited laterality defects (six HTX and one D-TGA). Three additional patients (4.1%) harbored other types of CHD-associated genomic alterations, which overall explained 13.7% (10/73) of the cohort. The contribution from recessive variants in our cohort is higher than 1.8% reported from a cohort of 2871 CHD subjects where 5.6% of subjects met the criteria for consanguinity.

CONCLUSIONS

Our WES screen of a Turkish consanguineous population with structural CHD revealed its unique genetic architecture. Six of seven damaging homozygous variants in CHD causal genes occur in the setting of laterality defects implies a strong contribution from consanguinity to these defects specifically. Our study thus provided valuable information about the genetic landscape of CHD in consanguineous families in Turkey.

Whole-exome sequencing reveals damaging gene variants associated with hypoalphalipoproteinemia

Low levels of high density lipoprotein-cholesterol (HDL-C) are associated with an elevated risk of arteriosclerotic coronary heart disease. Heritability of HDL-C levels is high. In this research discovery study, we used whole-exome sequencing to identify damaging gene variants that may play significant roles in determining HDL-C levels. We studied 204 individuals with a mean HDL-C level of 27.8 ± 6.4 mg/dl (range: 4-36 mg/dl). Data were analyzed by statistical gene burden testing and by filtering against candidate gene lists. We found 120 occurrences of probably damaging variants (116 heterozygous; four homozygous) among 45 of 104 recognized HDL candidate genes. Those with the highest prevalence of damaging variants were ABCA1 (n = 20), STAB1 (n = 9), OSBPL1A (n = 8), CPS1 (n = 8), CD36 (n = 7), LRP1 (n = 6), ABCA8 (n = 6), GOT2 (n = 5), AMPD3 (n = 5), WWOX (n = 4), and IRS1 (n = 4). Binomial analysis for damaging missense or loss-of-function variants identified the ABCA1 and LDLR genes at genome-wide significance. In conclusion, whole-exome sequencing of individuals with low HDL-C showed the burden of damaging rare variants in the ABCA1 and LDLR genes is particularly high and revealed numerous occurrences in HDL candidate genes, including many genes identified in genome-wide association study reports. Many of these genes are involved in cancer biology, which accords with epidemiologic findings of the association of HDL deficiency with increased risk of cancer, thus presenting a new area of interest in HDL genomics.

The risk of COVID-19 death is much greater and age dependent with type I IFN autoantibodies

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection fatality rate (IFR) doubles with every 5 y of age from childhood onward. Circulating autoantibodies neutralizing IFN-α, IFN-ω, and/or IFN-β are found in ∼20% of deceased patients across age groups, and in ∼1% of individuals aged <70 y and in >4% of those >70 y old in the general population. With a sample of 1,261 unvaccinated deceased patients and 34,159 individuals of the general population sampled before the pandemic, we estimated both IFR and relative risk of death (RRD) across age groups for individuals carrying autoantibodies neutralizing type I IFNs, relative to noncarriers. The RRD associated with any combination of autoantibodies was higher in subjects under 70 y old. For autoantibodies neutralizing IFN-α2 or IFN-ω, the RRDs were 17.0 (95% CI: 11.7 to 24.7) and 5.8 (4.5 to 7.4) for individuals <70 y and ≥70 y old, respectively, whereas, for autoantibodies neutralizing both molecules, the RRDs were 188.3 (44.8 to 774.4) and 7.2 (5.0 to 10.3), respectively. In contrast, IFRs increased with age, ranging from 0.17% (0.12 to 0.31) for individuals <40 y old to 26.7% (20.3 to 35.2) for those ≥80 y old for autoantibodies neutralizing IFN-α2 or IFN-ω, and from 0.84% (0.31 to 8.28) to 40.5% (27.82 to 61.20) for autoantibodies neutralizing both. Autoantibodies against type I IFNs increase IFRs, and are associated with high RRDs, especially when neutralizing both IFN-α2 and IFN-ω. Remarkably, IFRs increase with age, whereas RRDs decrease with age. Autoimmunity to type I IFNs is a strong and common predictor of COVID-19 death.

Centers for Mendelian Genomics: A decade of facilitating gene discovery

PURPOSE

Mendelian disease genomic research has undergone a massive transformation over the past decade. With increasing availability of exome and genome sequencing, the role of Mendelian research has expanded beyond data collection, sequencing, and analysis to worldwide data sharing and collaboration.

METHODS

Over the past 10 years, the National Institutes of Health-supported Centers for Mendelian Genomics (CMGs) have played a major role in this research and clinical evolution.

RESULTS

We highlight the cumulative gene discoveries facilitated by the program, biomedical research leveraged by the approach, and the larger impact on the research community. Beyond generating a list of gene-phenotype relationships and participating in widespread data sharing, the CMGs have created resources, tools, and training for the larger community to foster understanding of genes and genome variation. The CMGs have participated in a wide range of data sharing activities, including deposition of all eligible CMG data into the Analysis, Visualization, and Informatics Lab-space (AnVIL), sharing candidate genes through the Matchmaker Exchange and the CMG website, and sharing variants in Genotypes to Mendelian Phenotypes (Geno2MP) and VariantMatcher.

CONCLUSION

The work is far from complete; strengthening communication between research and clinical realms, continued development and sharing of knowledge and tools, and improving access to richly characterized data sets are all required to diagnose the remaining molecularly undiagnosed patients.

Genome-Wide De Novo Variants in Congenital Heart Disease Are Not Associated With Maternal Diabetes or Obesity

BACKGROUND

Congenital heart disease (CHD) is the most common anomaly at birth, with a prevalence of ≈1%. While infants born to mothers with diabetes or obesity have a 2- to 3-fold increased incidence of CHD, the cause of the increase is unknown. Damaging de novo variants (DNV) in coding regions are more common among patients with CHD, but genome-wide rates of coding and noncoding DNVs associated with these prenatal exposures have not been studied in patients with CHD.

METHODS

DNV frequencies were determined for 1812 patients with CHD who had whole-genome sequencing and prenatal history data available from the Pediatric Cardiac Genomics Consortium's CHD GENES study (Genetic Network). The frequency of DNVs was compared between subgroups using t test or linear model.

RESULTS

Among 1812 patients with CHD, the number of DNVs per patient was higher with maternal diabetes (76.5 versus 72.1, t test P=3.03×10^-11), but the difference was no longer significant after including parental ages in a linear model (paternal and maternal correction P=0.42). No interaction was observed between diabetes risk and parental age (paternal and maternal interaction P=0.80 and 0.68, respectively). No difference was seen in DNV count per patient based on maternal obesity (72.0 versus 72.2 for maternal body mass index <25 versus maternal body mass index >30, t test P=0.86).

CONCLUSIONS

After accounting for parental age, the offspring of diabetic or obese mothers have no increase in DNVs compared with other children with CHD. These results emphasize the role for other mechanisms in the cause of CHD associated with these prenatal exposures.

REGISTRATION

URL: https://clinicaltrials.gov; NCT01196182.

Impaired neurogenesis alters brain biomechanics in a neuroprogenitor-based genetic subtype of congenital hydrocephalus

Hydrocephalus, characterized by cerebral ventricular dilatation, is routinely attributed to primary defects in cerebrospinal fluid (CSF) homeostasis. This fosters CSF shunting as the leading reason for brain surgery in children despite considerable disease heterogeneity. In this study, by integrating human brain transcriptomics with whole-exome sequencing of 483 patients with congenital hydrocephalus (CH), we found convergence of CH risk genes in embryonic neuroepithelial stem cells. Of all CH risk genes, TRIM71/lin-41 harbors the most de novo mutations and is most specifically expressed in neuroepithelial cells. Mice harboring neuroepithelial cell-specific Trim71 deletion or CH-specific Trim71 mutation exhibit prenatal hydrocephalus. CH mutations disrupt TRIM71 binding to its RNA targets, causing premature neuroepithelial cell differentiation and reduced neurogenesis. Cortical hypoplasia leads to a hypercompliant cortex and secondary ventricular enlargement without primary defects in CSF circulation. These data highlight the importance of precisely regulated neuroepithelial cell fate for normal brain-CSF biomechanics and support a clinically relevant neuroprogenitor-based paradigm of CH.

Unexpected role of SIX1 variants in craniosynostosis: expanding the phenotype of SIX1-related disorders

BACKGROUND

Pathogenic heterozygous SIX1 variants (predominantly missense) occur in branchio-otic syndrome (BOS), but an association with craniosynostosis has not been reported.

METHODS

We investigated probands with craniosynostosis of unknown cause using whole exome/genome (n=628) or RNA (n=386) sequencing, and performed targeted resequencing of SIX1 in 615 additional patients. Expression of SIX1 protein in embryonic cranial sutures was examined in the Six1nLacZ/+ reporter mouse.

RESULTS

From 1629 unrelated cases with craniosynostosis we identified seven different SIX1 variants (three missense, including two de novo mutations, and four nonsense, one of which was also present in an affected twin). Compared with population data, enrichment of SIX1 loss-of-function variants was highly significant (p=0.00003). All individuals with craniosynostosis had sagittal suture fusion; additionally four had bilambdoid synostosis. Associated BOS features were often attenuated; some carrier relatives appeared non-penetrant. SIX1 is expressed in a layer basal to the calvaria, likely corresponding to the dura mater, and in the mid-sagittal mesenchyme.

CONCLUSION

Craniosynostosis is associated with heterozygous SIX1 variants, with possible enrichment of loss-of-function variants compared with classical BOS. We recommend screening of SIX1 in craniosynostosis, particularly when sagittal±lambdoid synostosis and/or any BOS phenotypes are present. These findings highlight the role of SIX1 in cranial suture homeostasis.

Whole exome sequencing identifies potential candidate genes for spina bifida derived from mouse models

Spina bifida (SB) is the second most common nonlethal congenital malformation. The existence of monogenic SB mouse models and human monogenic syndromes with SB features indicate that human SB may be caused by monogenic genes. We hypothesized that whole exome sequencing (WES) allows identification of potential candidate genes by (i) generating a list of 136 candidate genes for SB, and (ii) by unbiased exome-wide analysis. We generated a list of 136 potential candidate genes from three categories and evaluated WES data of 50 unrelated SB cases for likely deleterious variants in 136 potential candidate genes, and for potential SB candidate genes exome-wide. We identified 6 likely deleterious variants in 6 of the 136 potential SB candidate genes in 6 of the 50 SB cases, whereof 4 genes were derived from mouse models, 1 gene was derived from human nonsyndromic SB, and 1 gene was derived from candidate genes known to cause human syndromic SB. In addition, by unbiased exome-wide analysis, we identified 12 genes as potential candidates for SB. Identification of these 18 potential candidate genes in larger SB cohorts will help decide which ones can be considered as novel monogenic causes of human SB.

The Genomic and Phenotypic Landscape of Ichthyosis: An Analysis of 1000 Kindreds

Importance

Ichthyoses are clinically and genetically heterogeneous disorders characterized by scaly skin. Despite decades of investigation identifying pathogenic variants in more than 50 genes, clear genotype-phenotype associations have been difficult to establish.

Objective

To expand the genotypic and phenotypic spectra of ichthyosis and delineate genotype-phenotype associations.

Design, Setting, and Participants

This cohort study recruited an international group of individuals with ichthyosis and describes characteristic and distinguishing features of common genotypes, including genotype-phenotype associations, during a 10-year period from June 2011 to July 2021. Participants of all ages, races, and ethnicities were included and were enrolled worldwide from referral centers and patient advocacy groups. A questionnaire to assess clinical manifestations was completed by those with a genetic diagnosis.

Main Outcomes and Measures

Genetic analysis of saliva or blood DNA, a phenotyping questionnaire, and standardized clinical photographs. Descriptive statistics, such as frequency counts, were used to describe the cases in the cohort. Fisher exact tests identified significant genotype-phenotype associations.

Results

Results were reported for 1000 unrelated individuals enrolled from around the world (mean [SD] age, 50.0 [34.0] years; 524 [52.4%] were female, 427 [42.7%] were male, and 49 [4.9%] were not classified); 75% were from the US, 12% from Latin America, 4% from Canada, 3% from Europe, 3% from Asia, 2% from Africa, 1% from the Middle East, and 1% from Australia and New Zealand. A total of 266 novel disease-associated variants in 32 genes were identified among 869 kindreds. Of these, 241 (91%) pathogenic variants were found through multiplex amplicon sequencing and 25 (9%) through exome sequencing. Among the 869 participants with a genetic diagnosis, 304 participants (35%) completed the phenotyping questionnaire. Analysis of clinical manifestations in these 304 individuals revealed that pruritus, hypohydrosis, skin pain, eye problems, skin odor, and skin infections were the most prevalent self-reported features. Genotype-phenotype association analysis revealed that the presence of a collodion membrane at birth (odds ratio [OR], 6.7; 95% CI, 3.0-16.7; P < .001), skin odor (OR, 2.8; 95% CI, 1.1-6.8; P = .02), hearing problems (OR, 2.9; 95% CI, 1.6-5.5; P < .001), eye problems (OR, 3.0; 95% CI, 1.5-6.0; P < .001), and alopecia (OR, 4.6; 95% CI, 2.4-9.0; P < .001) were significantly associated with TGM1 variants compared with other ichthyosis genotypes studied. Skin pain (OR, 6.8; 95% CI, 1.6-61.2; P = .002), odor (OR, 5.7; 95% CI, 2.0-19.7; P < .001), and infections (OR, 3.1; 95% CI, 1.4-7.7; P = .03) were significantly associated with KRT10 pathogenic variants compared with disease-associated variants in other genes that cause ichthyosis. Pathogenic variants were identified in 869 (86.9%) participants. Most of the remaining individuals had unique phenotypes, enabling further genetic discovery.

Conclusions and Relevance

This cohort study expands the genotypic and phenotypic spectrum of ichthyosis, establishing associations between clinical manifestations and genotypes. Collectively, the findings may help improve clinical assessment, assist with developing customized management plans, and improve clinical course prognostication.

PPIL4 is essential for brain angiogenesis and implicated in intracranial aneurysms in humans

Intracranial aneurysm (IA) rupture leads to subarachnoid hemorrhage, a sudden-onset disease that often causes death or severe disability. Although genome-wide association studies have identified common genetic variants that increase IA risk moderately, the contribution of variants with large effect remains poorly defined. Using whole-exome sequencing, we identified significant enrichment of rare, deleterious mutations in PPIL4, encoding peptidyl-prolyl cis-trans isomerase-like 4, in both familial and index IA cases. Ppil4 depletion in vertebrate models causes intracerebral hemorrhage, defects in cerebrovascular morphology and impaired Wnt signaling. Wild-type, but not IA-mutant, PPIL4 potentiates Wnt signaling by binding JMJD6, a known angiogenesis regulator and Wnt activator. These findings identify a novel PPIL4-dependent Wnt signaling mechanism involved in brain-specific angiogenesis and maintenance of cerebrovascular integrity and implicate PPIL4 gene mutations in the pathogenesis of IA.

Autoantibodies neutralizing type I IFNs are present in 4% of uninfected individuals over 70 years old and account for 20% of COVID-19 deaths

Circulating autoantibodies (auto-Abs) neutralizing high concentrations (10 ng/mL, in plasma diluted 1 to 10) of IFN-α and/or -ω are found in about 10% of patients with critical COVID-19 pneumonia, but not in subjects with asymptomatic infections. We detect auto-Abs neutralizing 100-fold lower, more physiological, concentrations of IFN-α and/or -ω (100 pg/mL, in 1/10 dilutions of plasma) in 13.6% of 3,595 patients with critical COVID-19, including 21% of 374 patients > 80 years, and 6.5% of 522 patients with severe COVID-19. These antibodies are also detected in 18% of the 1,124 deceased patients (aged 20 days-99 years; mean: 70 years). Moreover, another 1.3% of patients with critical COVID-19 and 0.9% of the deceased patients have auto-Abs neutralizing high concentrations of IFN-β. We also show, in a sample of 34,159 uninfected subjects from the general population, that auto-Abs neutralizing high concentrations of IFN-α and/or -ω are present in 0.18% of individuals between 18 and 69 years, 1.1% between 70 and 79 years, and 3.4% >80 years. Moreover, the proportion of subjects carrying auto-Abs neutralizing lower concentrations is greater in a subsample of 10,778 uninfected individuals: 1% of individuals <70 years, 2.3% between 70 and 80 years, and 6.3% >80 years. By contrast, auto-Abs neutralizing IFN-β do not become more frequent with age. Auto-Abs neutralizing type I IFNs predate SARS-CoV-2 infection and sharply increase in prevalence after the age of 70 years. They account for about 20% of both critical COVID-19 cases in the over-80s, and total fatal COVID-19 cases.

X-linked recessive TLR7 deficiency in ~1% of men under 60 years old with life-threatening COVID-19

Autosomal inborn errors of type I IFN immunity and autoantibodies against these cytokines underlie at least 10% of critical COVID-19 pneumonia cases. We report very rare, biochemically deleterious X-linked TLR7 variants in 16 unrelated male individuals aged 7 to 71 years (mean: 36.7 years) from a cohort of 1,202 male patients aged 0.5 to 99 years (mean: 52.9 years) with unexplained critical COVID-19 pneumonia. None of the 331 asymptomatically or mildly infected male individuals aged 1.3 to 102 years (mean: 38.7 years) tested carry such TLR7 variants (p = 3.5 × 10-5). The phenotypes of five hemizygous relatives of index cases infected with SARS-CoV-2 include asymptomatic or mild infection (n=2, 5 and 38 years), or moderate (n=1, 5 years), severe (n=1, 27 years), or critical (n=1, 29 years) pneumonia. Two boys (aged 7 and 12 years) from a cohort of 262 male patients with severe COVID-19 pneumonia (mean: 51.0 years) are hemizygous for a deleterious TLR7 variant. The cumulative allele frequency for deleterious TLR7 variants in the male general population is < 6.5x10-4 We also show that blood B cell lines and myeloid cell subsets from the patients do not respond to TLR7 stimulation, a phenotype rescued by wild-type TLR7 The patients' blood plasmacytoid dendritic cells (pDCs) produce low levels of type I IFNs in response to SARS-CoV-2. Overall, X-linked recessive TLR7 deficiency is a highly penetrant genetic etiology of critical COVID-19 pneumonia, in about 1.8% of male patients below the age of 60 years. Human TLR7 and pDCs are essential for protective type I IFN immunity against SARS-CoV-2 in the respiratory tract.

Exome survey of individuals affected by VATER/VACTERL with renal phenotypes identifies phenocopies and novel candidate genes

The acronym VATER/VACTERL refers to the rare nonrandom association of the following component features (CFs): vertebral defects (V), anorectal malformations (ARM) (A), cardiac anomalies (C), tracheoesophageal fistula with or without esophageal atresia (TE), renal malformations (R), and limb anomalies (L). For the clinical diagnosis, the presence of at least three CFs is required, individuals presenting with only two CFs have been categorized as VATER/VACTERL-like. The majority of VATER/VACTERL individuals displays a renal phenotype. Hitherto, variants in FGF8, FOXF1, HOXD13, LPP, TRAP1, PTEN, and ZIC3 have been associated with the VATER/VACTERL association; however, large-scale re-sequencing could only confirm TRAP1 and ZIC3 as VATER/VACTERL disease genes, both associated with a renal phenotype. In this study, we performed exome sequencing in 21 individuals and their families with a renal VATER/VACTERL or VATER/VACTERL-like phenotype to identify potentially novel genetic causes. Exome analysis identified biallelic and X-chromosomal hemizygous potentially pathogenic variants in six individuals (29%) in B9D1, FREM1, ZNF157, SP8, ACOT9, and TTLL11, respectively. The online tool GeneMatcher revealed another individual with a variant in ZNF157. Our study suggests six biallelic and X-chromosomal hemizygous VATER/VACTERL disease genes implicating all six genes in the expression of human renal malformations.

DIAPH1 Variants in Non-East Asian Patients With Sporadic Moyamoya Disease

Importance

Moyamoya disease (MMD), a progressive vasculopathy leading to narrowing and ultimate occlusion of the intracranial internal carotid arteries, is a cause of childhood stroke. The cause of MMD is poorly understood, but genetic factors play a role. Several familial forms of MMD have been identified, but the cause of most cases remains elusive, especially among non-East Asian individuals.

Objective

To assess whether ultrarare de novo and rare, damaging transmitted variants with large effect sizes are associated with MMD risk.

Design, Setting, and Participants

A genetic association study was conducted using whole-exome sequencing case-parent MMD trios in a small discovery cohort collected over 3.5 years (2016-2019); data were analyzed in 2020. Medical records from US hospitals spanning a range of 1 month to 1.5 years were reviewed for phenotyping. Exomes from a larger validation cohort were analyzed to identify additional rare, large-effect variants in the top candidate gene. Participants included patients with MMD and, when available, their parents. All participants who met criteria and were presented with the option to join the study agreed to do so; none were excluded. Twenty-four probands (22 trios and 2 singletons) composed the discovery cohort, and 84 probands (29 trios and 55 singletons) composed the validation cohort.

Main Outcomes and Measures

Gene variants were identified and filtered using stringent criteria. Enrichment and case-control tests assessed gene-level variant burden. In silico modeling estimated the probability of variant association with protein structure. Integrative genomics assessed expression patterns of MMD risk genes derived from single-cell RNA sequencing data of human and mouse brain tissue.

Results

Of the 24 patients in the discovery cohort, 14 (58.3%) were men and 18 (75.0%) were of European ancestry. Three of 24 discovery cohort probands contained 2 do novo (1-tailed Poisson P = 1.1 × 10-6) and 1 rare, transmitted damaging variant (12.5% of cases) in DIAPH1 (mammalian diaphanous-1), a key regulator of actin remodeling in vascular cells and platelets. Four additional ultrarare damaging heterozygous DIAPH1 variants (3 unphased) were identified in 3 other patients in an 84-proband validation cohort (73.8% female, 77.4% European). All 6 patients were non-East Asian. Compound heterozygous variants were identified in ena/vasodilator-stimulated phosphoproteinlike protein EVL, a mammalian diaphanous-1 interactor that regulates actin polymerization. DIAPH1 and EVL mutant probands had severe, bilateral MMD associated with transfusion-dependent thrombocytopenia. DIAPH1 and other MMD risk genes are enriched in mural cells of midgestational human brain. The DIAPH1 coexpression network converges in vascular cell actin cytoskeleton regulatory pathways.

Conclusions and Relevance

These findings provide the largest collection to date of non-East Asian individuals with sporadic MMD harboring pathogenic variants in the same gene. The results suggest that DIAPH1 is a novel MMD risk gene and impaired vascular cell actin remodeling in MMD pathogenesis, with diagnostic and therapeutic ramifications.

GIMAP5 maintains liver endothelial cell homeostasis and prevents portal hypertension

Portal hypertension is a major contributor to decompensation and death from liver disease, a global health problem. Here, we demonstrate homozygous damaging mutations in GIMAP5, a small organellar GTPase, in four families with unexplained portal hypertension. We show that GIMAP5 is expressed in hepatic endothelial cells and that its loss in both humans and mice results in capillarization of liver sinusoidal endothelial cells (LSECs); this effect is also seen when GIMAP5 is selectively deleted in endothelial cells. Single-cell RNA-sequencing analysis in a GIMAP5-deficient mouse model reveals replacement of LSECs with capillarized endothelial cells, a reduction of macrovascular hepatic endothelial cells, and places GIMAP5 upstream of GATA4, a transcription factor required for LSEC specification. Thus, GIMAP5 is a critical regulator of liver endothelial cell homeostasis and, when absent, produces portal hypertension. These findings provide new insight into the pathogenesis of portal hypertension, a major contributor to morbidity and mortality from liver disease.

Genetic Defects in DNAH2 Underlie Male Infertility With Multiple Morphological Abnormalities of the Sperm Flagella in Humans and Mice

Asthenozoospermia accounts for over 80% of primary male infertility cases. Reduced sperm motility in asthenozoospermic patients are often accompanied by teratozoospermia, or defective sperm morphology, with varying severity. Multiple morphological abnormalities of the flagella (MMAF) is one of the most severe forms of asthenoteratozoospermia, characterized by heterogeneous flagellar abnormalities. Among various genetic factors known to cause MMAF, multiple variants in the DNAH2 gene are reported to underlie MMAF in humans. However, the pathogenicity by DNAH2 mutations remains largely unknown. In this study, we identified a novel recessive variant (NM_020877:c.12720G > T;p.W4240C) in DNAH2 by whole-exome sequencing, which fully co-segregated with the infertile male members in a consanguineous Pakistani family diagnosed with asthenozoospermia. 80–90% of the sperm from the patients are morphologically abnormal, and in silico analysis models reveal that the non-synonymous variant substitutes a residue in dynein heavy chain domain and destabilizes DNAH2. To better understand the pathogenicity of various DNAH2 variants underlying MMAF in general, we functionally characterized Dnah2-mutant mice generated by CRISPR/Cas9 genome editing. Dnah2-null males, but not females, are infertile. Dnah2-null sperm cells display absent, short, bent, coiled, and/or irregular flagella consistent with the MMAF phenotype. We found misexpression of centriolar proteins and delocalization of annulus proteins in Dnah2-null spermatids and sperm, suggesting dysregulated flagella development in spermiogenesis. Scanning and transmission electron microscopy analyses revealed that flagella ultrastructure is severely disorganized in Dnah2-null sperm. Absence of DNAH2 compromises the expression of other axonemal components such as DNAH1 and RSPH3. Our results demonstrate that DNAH2 is essential for multiple steps in sperm flagella formation and provide insights into molecular and cellular mechanisms of MMAF pathogenesis.

Association of Damaging Variants in Genes With Increased Cancer Risk Among Patients With Congenital Heart Disease

Importance

Patients with congenital heart disease (CHD), the most common birth defect, have increased risks for cancer. Identification of the variables that contribute to cancer risk is essential for recognizing patients with CHD who warrant longitudinal surveillance and early interventions.

Objective

To compare the frequency of damaging variants in cancer risk genes among patients with CHD and control participants and identify associated clinical variables in patients with CHD who have cancer risk variants.

Design, Setting, and Participants

This multicenter case-control study included participants with CHD who had previously been recruited to the Pediatric Cardiac Genomics Consortium based on presence of structural cardiac anomaly without genetic diagnosis at the time of enrollment. Permission to use published sequencing data from unaffected adult participants was obtained from 2 parent studies. Data were collected for this study from December 2010 to April 2019.

Exposures

Presence of rare (allele frequency, <1 × 10-5) loss-of-function (LoF) variants in cancer risk genes.

Main Outcomes and Measures

Frequency of LoF variants in cancer risk genes (defined in the Catalogue of Somatic Mutations in Cancer-Cancer Gene Consensus database), were statistically assessed by binomial tests in patients with CHD and control participants.

Results

A total of 4443 individuals with CHD (mean [range] age, 13.0 [0-84] years; 2225 of 3771 with reported sex [59.0%] male) and 9808 control participants (mean [range] age, 52.1 [1-92] years; 4967 of 9808 [50.6%] male) were included. The frequency of LoF variants in regulatory cancer risk genes was significantly higher in patients with CHD than control participants (143 of 4443 [3.2%] vs 166 of 9808 [1.7%]; odds ratio [OR], 1.93 [95% CI, 1.54-2.42]; P = 1.38 × 10-12), and among CHD genes previously associated with cancer risk (58 of 4443 [1.3%] vs 18 of 9808 [0.18%]; OR, 7.2 [95% CI, 4.2-12.2]; P < 2.2 × 10-16). The LoF variants were also nominally increased in 14 constrained cancer risk genes with high expression in the developing heart. Seven of these genes (ARHGEF12, CTNNB1, LPP, MLLT4, PTEN, TCF12, and TFRC) harbored LoF variants in multiple patients with unexplained CHD. The highest rates for LoF variants in cancer risk genes occurred in patients with CHD and extracardiac anomalies (248 of 1482 individuals [16.7%]; control: 1099 of 9808 individuals [11.2%]; OR, 1.59 [95% CI, 1.37-1.85]; P = 1.3 × 10-10) and/or neurodevelopmental delay (209 of 1393 individuals [15.0%]; control: 1099 of 9808 individuals [11.2%]; OR, 1.40 [95% CI, 1.19-1.64]; P = 9.6 × 10-6).

Conclusions and Relevance

Genotypes of CHD may account for increased cancer risks. In this cohort, damaging variants were prominent in the 216 genes that predominantly encode regulatory proteins. Consistent with their fundamental developmental functions, patients with CHD and damaging variants in these genes often had extracardiac manifestations. These data may also implicate cancer risk genes that are repeatedly varied in patients with unexplained CHD as CHD genes.

Mutations in PRDM15 Are a Novel Cause of Galloway-Mowat Syndrome

BACKGROUND

Galloway-Mowat syndrome (GAMOS) is characterized by neurodevelopmental defects and a progressive nephropathy, which typically manifests as steroid-resistant nephrotic syndrome. The prognosis of GAMOS is poor, and the majority of children progress to renal failure. The discovery of monogenic causes of GAMOS has uncovered molecular pathways involved in the pathogenesis of disease.

METHODS

Homozygosity mapping, whole-exome sequencing, and linkage analysis were used to identify mutations in four families with a GAMOS-like phenotype, and high-throughput PCR technology was applied to 91 individuals with GAMOS and 816 individuals with isolated nephrotic syndrome. In vitro and in vivo studies determined the functional significance of the mutations identified.

RESULTS

Three biallelic variants of the transcriptional regulator PRDM15 were detected in six families with proteinuric kidney disease. Four families with a variant in the protein's zinc-finger (ZNF) domain have additional GAMOS-like features, including brain anomalies, cardiac defects, and skeletal defects. All variants destabilize the PRDM15 protein, and the ZNF variant additionally interferes with transcriptional activation. Morpholino oligonucleotide-mediated knockdown of Prdm15 in Xenopus embryos disrupted pronephric development. Human wild-type PRDM15 RNA rescued the disruption, but the three PRDM15 variants did not. Finally, CRISPR-mediated knockout of PRDM15 in human podocytes led to dysregulation of several renal developmental genes.

CONCLUSIONS

Variants in PRDM15 can cause either isolated nephrotic syndrome or a GAMOS-type syndrome on an allelic basis. PRDM15 regulates multiple developmental kidney genes, and is likely to play an essential role in renal development in humans.

Mechanisms of Congenital Heart Disease Caused by NAA15 Haploinsufficiency

[Figure: see text].

De novo TRIM8 variants impair its protein localization to nuclear bodies and cause developmental delay, epilepsy, and focal segmental glomerulosclerosis

Focal segmental glomerulosclerosis (FSGS) is the main pathology underlying steroid-resistant nephrotic syndrome (SRNS) and a leading cause of chronic kidney disease. Monogenic forms of pediatric SRNS are predominantly caused by recessive mutations, while the contribution of de novo variants (DNVs) to this trait is poorly understood. Using exome sequencing (ES) in a proband with FSGS/SRNS, developmental delay, and epilepsy, we discovered a nonsense DNV in TRIM8, which encodes the E3 ubiquitin ligase tripartite motif containing 8. To establish whether TRIM8 variants represent a cause of FSGS, we aggregated exome/genome-sequencing data for 2,501 pediatric FSGS/SRNS-affected individuals and 48,556 control subjects, detecting eight heterozygous TRIM8 truncating variants in affected subjects but none in control subjects (p = 3.28 × 10-11). In all six cases with available parental DNA, we demonstrated de novo inheritance (p = 2.21 × 10-15). Reverse phenotyping revealed neurodevelopmental disease in all eight families. We next analyzed ES from 9,067 individuals with epilepsy, yielding three additional families with truncating TRIM8 variants. Clinical review revealed FSGS in all. All TRIM8 variants cause protein truncation clustering within the last exon between residues 390 and 487 of the 551 amino acid protein, indicating a correlation between this syndrome and loss of the TRIM8 C-terminal region. Wild-type TRIM8 overexpressed in immortalized human podocytes and neuronal cells localized to nuclear bodies, while constructs harboring patient-specific variants mislocalized diffusely to the nucleoplasm. Co-localization studies demonstrated that Gemini and Cajal bodies frequently abut a TRIM8 nuclear body. Truncating TRIM8 DNVs cause a neuro-renal syndrome via aberrant TRIM8 localization, implicating nuclear bodies in FSGS and developmental brain disease.

Recessive NOS1AP variants impair actin remodeling and cause glomerulopathy in humans and mice

Nephrotic syndrome (NS) is a leading cause of chronic kidney disease. We found recessive NOS1AP variants in two families with early-onset NS by exome sequencing. Overexpression of wild-type (WT) NOS1AP, but not cDNA constructs bearing patient variants, increased active CDC42 and promoted filopodia and podosome formation. Pharmacologic inhibition of CDC42 or its effectors, formin proteins, reduced NOS1AP-induced filopodia formation. NOS1AP knockdown reduced podocyte migration rate (PMR), which was rescued by overexpression of WT Nos1ap but not by constructs bearing patient variants. PMR in NOS1AP knockdown podocytes was also rescued by constitutively active CDC42Q61L or the formin DIAPH3 Modeling a NOS1AP patient variant in knock-in human kidney organoids revealed malformed glomeruli with increased apoptosis. Nos1apEx3-/Ex3- mice recapitulated the human phenotype, exhibiting proteinuria, foot process effacement, and glomerulosclerosis. These findings demonstrate that recessive NOS1AP variants impair CDC42/DIAPH-dependent actin remodeling, cause aberrant organoid glomerulogenesis, and lead to a glomerulopathy in humans and mice.

DAAM2 Variants Cause Nephrotic Syndrome via Actin Dysregulation

The discovery of >60 monogenic causes of nephrotic syndrome (NS) has revealed a central role for the actin regulators RhoA/Rac1/Cdc42 and their effectors, including the formin INF2. By whole-exome sequencing (WES), we here discovered bi-allelic variants in the formin DAAM2 in four unrelated families with steroid-resistant NS. We show that DAAM2 localizes to the cytoplasm in podocytes and in kidney sections. Further, the variants impair DAAM2-dependent actin remodeling processes: wild-type DAAM2 cDNA, but not cDNA representing missense variants found in individuals with NS, rescued reduced podocyte migration rate (PMR) and restored reduced filopodia formation in shRNA-induced DAAM2-knockdown podocytes. Filopodia restoration was also induced by the formin-activating molecule IMM-01. DAAM2 also co-localizes and co-immunoprecipitates with INF2, which is intriguing since variants in both formins cause NS. Using in vitro bulk and TIRF microscopy assays, we find that DAAM2 variants alter actin assembly activities of the formin. In a Xenopus daam2-CRISPR knockout model, we demonstrate actin dysregulation in vivo and glomerular maldevelopment that is rescued by WT-DAAM2 mRNA. We conclude that DAAM2 variants are a likely cause of monogenic human SRNS due to actin dysregulation in podocytes. Further, we provide evidence that DAAM2-associated SRNS may be amenable to treatment using actin regulating compounds.