Rare variant enrichment analysis supports GREB1L as a contributory driver gene in the etiology of Mayer-Rokitansky-Küster-Hauser syndrome.

Mayer-Rokitansky-Küster-Hauser (MRKH) syndrome is characterized by aplasia of the female reproductive tract; the syndrome can include renal anomalies, absence or dysgenesis, and skeletal anomalies. While functional models have elucidated several candidate genes, only WNT4 (MIM: 603490) variants have been definitively associated with a subtype of MRKH with hyperandrogenism (MIM: 158330). DNA from 148 clinically diagnosed MRKH probands across 144 unrelated families and available family members from North America, Europe, and South America were exome sequenced (ES) and by family-based genomics analyzed for rare likely deleterious variants. A replication cohort consisting of 442 Han Chinese individuals with MRKH was used to further reproduce GREB1L findings in diverse genetic backgrounds. Proband and OMIM phenotypes annotated using the Human Phenotype Ontology were analyzed to quantitatively delineate the phenotypic spectrum associated with GREB1L variant alleles found in our MRKH cohort and those previously published. This study reports 18 novel GREB1L variant alleles, 16 within a multiethnic MRKH cohort and two within a congenital scoliosis cohort. Cohort-wide analyses for a burden of rare variants within a single gene identified likely damaging variants in GREB1L (MIM: 617782), a known disease gene for renal hypoplasia and uterine abnormalities (MIM: 617805), in 16 of 590 MRKH probands. GREB1L variant alleles, including a CNV null allele, were found in 8 MRKH type 1 probands and 8 MRKH type II probands. This study used quantitative phenotypic analyses in a worldwide multiethnic cohort to identify and strengthen the association of GREB1L to isolated uterine agenesis (MRKH type I) and syndromic MRKH type II.

Novel LSS variants in alopecia and intellectual disability syndrome: New case report and clinical spectrum of LSS-related rare disease traits.

Pathogenic biallelic variants in LSS are associated with three Mendelian rare disease traits including congenital cataract type 44, autosomal recessive hypotrichosis type 14, and alopecia-intellectual disability syndrome type 4 (APMR4). We performed trio research exome sequencing on a family with a four-year-old male with global developmental delay, epilepsy and striking alopecia, and identified novel compound heterozygous LSS splice site (c.14+2T>C) and missense (c.1357 G>A; p.V453L) variant alleles. Rare features associated with APMR4 such as cryptorchidism, micropenis, mild cortical brain atrophy and thin corpus callosum were detected. Previously unreported APMR4 findings including cerebellar involvement in the form of unsteady ataxic gait, small vermis with prominent folia, were noted. A review of all reported variants to date in 29 families with LSS-related phenotypes showed an emerging genotype-phenotype correlation. Our report potentially expands LSS-related phenotypic spectrum and highlights the importance of performing brain imaging in LSS-related conditions.

SRSF1 haploinsufficiency is responsible for a syndromic developmental disorder associated with intellectual disability.

SRSF1 (also known as ASF/SF2) is a non-small nuclear ribonucleoprotein (non-snRNP) that belongs to the arginine/serine (R/S) domain family. It recognizes and binds to mRNA, regulating both constitutive and alternative splicing. The complete loss of this proto-oncogene in mice is embryonically lethal. Through international data sharing, we identified 17 individuals (10 females and 7 males) with a neurodevelopmental disorder (NDD) with heterozygous germline SRSF1 variants, mostly de novo, including three frameshift variants, three nonsense variants, seven missense variants, and two microdeletions within region 17q22 encompassing SRSF1. Only in one family, the de novo origin could not be established. All individuals featured a recurrent phenotype including developmental delay and intellectual disability (DD/ID), hypotonia, neurobehavioral problems, with variable skeletal (66.7%) and cardiac (46%) anomalies. To investigate the functional consequences of SRSF1 variants, we performed in silico structural modeling, developed an in vivo splicing assay in Drosophila, and carried out episignature analysis in blood-derived DNA from affected individuals. We found that all loss-of-function and 5 out of 7 missense variants were pathogenic, leading to a loss of SRSF1 splicing activity in Drosophila, correlating with a detectable and specific DNA methylation episignature. In addition, our orthogonal in silico, in vivo, and epigenetics analyses enabled the separation of clearly pathogenic missense variants from those with uncertain significance. Overall, these results indicated that haploinsufficiency of SRSF1 is responsible for a syndromic NDD with ID due to a partial loss of SRSF1-mediated splicing activity.

Metabolic individuality: Limitations, challenges, and potential for clinical utility.

In Nature Medicine, Surendran and colleagues recently reported the analysis of human plasma metabolomic data for 913 metabolites in ∼20,000 individuals, identifying 2,599 metabolite-genetic variant associations and >400 metabolite signatures comprised of jointly regulated metabolites. This extensive atlas of variant-metabolite relationships reveals novel genomic mechanisms driving metabolic phenotypes.

A biallelic frameshift indel in PPP1R35 as a cause of primary microcephaly.

Protein phosphatase 1 regulatory subunit 35 (PPP1R35) encodes a centrosomal protein required for recruiting microtubule-binding elongation machinery. Several proteins in this centriole biogenesis pathway correspond to established primary microcephaly (MCPH) genes, and multiple model organism studies hypothesize PPP1R35 as a candidate MCPH gene. Here, using exome sequencing (ES) and family-based rare variant analyses, we report a homozygous, frameshifting indel deleting the canonical stop codon in the last exon of PPP1R35 [Chr7: c.753_*3delGGAAGCGTAGACCinsCG (p.Trp251Cysfs*22)]; the variant allele maps in a 3.7 Mb block of absence of heterozygosity (AOH) in a proband with severe MCPH (-4.3 SD at birth, -6.1 SD by 42 months), pachygyria, and global developmental delay from a consanguineous Turkish kindred. Droplet digital PCR (ddPCR) confirmed mutant mRNA expression in fibroblasts. In silico prediction of the translation of mutant PPP1R35 is expected to be elongated by 18 amino acids before encountering a downstream stop codon. This complex indel allele is absent in public databases (ClinVar, gnomAD, ARIC, 1000 genomes) and our in-house database of 14,000+ exomes including 1800+ Turkish exomes supporting predicted pathogenicity. Comprehensive literature searches for PPP1R35 variants yielded two probands affected with severe microcephaly (-15 SD and -12 SD) with the same homozygous indel from a single, consanguineous, Iranian family from a cohort of 404 predominantly Iranian families. The lack of heterozygous cases in two large cohorts representative of the genetic background of these two families decreased our suspicion of a founder allele and supports the contention of a recurrent mutation. We propose two potential secondary structure mutagenesis models for the origin of this variant allele mediated by hairpin formation between complementary GC rich segments flanking the stop codon via secondary structure mutagenesis.

An ELF4 hypomorphic variant results in NK cell deficiency.

NK cell deficiencies (NKD) are a type of primary immune deficiency in which the major immunologic abnormality affects NK cell number, maturity, or function. Since NK cells contribute to immune defense against virally infected cells, patients with NKD experience higher susceptibility to chronic, recurrent, and fatal viral infections. An individual with recurrent viral infections and mild hypogammaglobulinemia was identified to have an X-linked damaging variant in the transcription factor gene ELF4. The variant does not decrease expression but disrupts ELF4 protein interactions and DNA binding, reducing transcriptional activation of target genes and selectively impairing ELF4 function. Corroborating previous murine models of ELF4 deficiency (Elf4-/-) and using a knockdown human NK cell line, we determined that ELF4 is necessary for normal NK cell development, terminal maturation, and function. Through characterization of the NK cells of the proband, expression of the proband's variant in Elf4-/- mouse hematopoietic precursor cells, and a human in vitro NK cell maturation model, we established this ELF4 variant as a potentially novel cause of NKD.

TCEAL1 loss-of-function results in an X-linked dominant neurodevelopmental syndrome and drives the neurological disease trait in Xq22.2 deletions.

An Xq22.2 region upstream of PLP1 has been proposed to underly a neurological disease trait when deleted in 46,XX females. Deletion mapping revealed that heterozygous deletions encompassing the smallest region of overlap (SRO) spanning six Xq22.2 genes (BEX3, RAB40A, TCEAL4, TCEAL3, TCEAL1, and MORF4L2) associate with an early-onset neurological disease trait (EONDT) consisting of hypotonia, intellectual disability, neurobehavioral abnormalities, and dysmorphic facial features. None of the genes within the SRO have been associated with monogenic disease in OMIM. Through local and international collaborations facilitated by GeneMatcher and Matchmaker Exchange, we have identified and herein report seven de novo variants involving TCEAL1 in seven unrelated families: three hemizygous truncating alleles; one hemizygous missense allele; one heterozygous TCEAL1 full gene deletion; one heterozygous contiguous deletion of TCEAL1, TCEAL3, and TCEAL4; and one heterozygous frameshift variant allele. Variants were identified through exome or genome sequencing with trio analysis or through chromosomal microarray. Comparison with previously reported Xq22 deletions encompassing TCEAL1 identified a more-defined syndrome consisting of hypotonia, abnormal gait, developmental delay/intellectual disability especially affecting expressive language, autistic-like behavior, and mildly dysmorphic facial features. Additional features include strabismus, refractive errors, variable nystagmus, gastroesophageal reflux, constipation, dysmotility, recurrent infections, seizures, and structural brain anomalies. An additional maternally inherited hemizygous missense allele of uncertain significance was identified in a male with hypertonia and spasticity without syndromic features. These data provide evidence that TCEAL1 loss of function causes a neurological rare disease trait involving significant neurological impairment with features overlapping the EONDT phenotype in females with the Xq22 deletion.

The multiple de novo copy number variant (MdnCNV) phenomenon presents with peri-zygotic DNA mutational signatures and multilocus pathogenic variation.

BACKGROUND: The multiple de novo copy number variant (MdnCNV) phenotype is described by having four or more constitutional de novo CNVs (dnCNVs) arising independently throughout the human genome within one generation. It is a rare peri-zygotic mutational event, previously reported to be seen once in every 12,000 individuals referred for genome-wide chromosomal microarray analysis due to congenital abnormalities. These rare families provide a unique opportunity to understand the genetic factors of peri-zygotic genome instability and the impact of dnCNV on human diseases. METHODS: Chromosomal microarray analysis (CMA), array-based comparative genomic hybridization, short- and long-read genome sequencing (GS) were performed on the newly identified MdnCNV family to identify de novo mutations including dnCNVs, de novo single-nucleotide variants (dnSNVs), and indels. Short-read GS was performed on four previously published MdnCNV families for dnSNV analysis. Trio-based rare variant analysis was performed on the newly identified individual and four previously published MdnCNV families to identify potential genetic etiologies contributing to the peri-zygotic genomic instability. Lin semantic similarity scores informed quantitative human phenotype ontology analysis on three MdnCNV families to identify gene(s) driving or contributing to the clinical phenotype. RESULTS: In the newly identified MdnCNV case, we revealed eight de novo tandem duplications, each ~ 1 Mb, with microhomology at 6/8 breakpoint junctions. Enrichment of de novo single-nucleotide variants (SNV; 6/79) and de novo indels (1/12) was found within 4 Mb of the dnCNV genomic regions. An elevated post-zygotic SNV mutation rate was observed in MdnCNV families. Maternal rare variant analyses identified three genes in distinct families that may contribute to the MdnCNV phenomenon. Phenotype analysis suggests that gene(s) within dnCNV regions contribute to the observed proband phenotype in 3/3 cases. CNVs in two cases, a contiguous gene duplication encompassing PMP22 and RAI1 and another duplication affecting NSD1 and SMARCC2, contribute to the clinically observed phenotypic manifestations. CONCLUSIONS: Characteristic features of dnCNVs reported here are consistent with a microhomology-mediated break-induced replication (MMBIR)-driven mechanism during the peri-zygotic period. Maternal genetic variants in DNA repair genes potentially contribute to peri-zygotic genomic instability. Variable phenotypic features were observed across a cohort of three MdnCNV probands, and computational quantitative phenotyping revealed that two out of three had evidence for the contribution of more than one genetic locus to the proband's phenotype supporting the hypothesis of de novo multilocus pathogenic variation (MPV) in those families.

Bi-allelic CAMSAP1 variants cause a clinically recognizable neuronal migration disorder.

Non-centrosomal microtubules are essential cytoskeletal filaments that are important for neurite formation, axonal transport, and neuronal migration. They require stabilization by microtubule minus-end-targeting proteins including the CAMSAP family of molecules. Using exome sequencing on samples from five unrelated families, we show that bi-allelic CAMSAP1 loss-of-function variants cause a clinically recognizable, syndromic neuronal migration disorder. The cardinal clinical features of the syndrome include a characteristic craniofacial appearance, primary microcephaly, severe neurodevelopmental delay, cortical visual impairment, and seizures. The neuroradiological phenotype comprises a highly recognizable combination of classic lissencephaly with a posterior more severe than anterior gradient similar to PAFAH1B1(LIS1)-related lissencephaly and severe hypoplasia or absence of the corpus callosum; dysplasia of the basal ganglia, hippocampus, and midbrain; and cerebellar hypodysplasia, similar to the tubulinopathies, a group of monogenic tubulin-associated disorders of cortical dysgenesis. Neural cell rosette lineages derived from affected individuals displayed findings consistent with these phenotypes, including abnormal morphology, decreased cell proliferation, and neuronal differentiation. Camsap1-null mice displayed increased perinatal mortality, and RNAScope studies identified high expression levels in the brain throughout neurogenesis and in facial structures, consistent with the mouse and human neurodevelopmental and craniofacial phenotypes. Together our findings confirm a fundamental role of CAMSAP1 in neuronal migration and brain development and define bi-allelic variants as a cause of a clinically distinct neurodevelopmental disorder in humans and mice.

Loss of non-motor kinesin KIF26A causes congenital brain malformations via dysregulated neuronal migration and axonal growth as well as apoptosis.

Kinesins are canonical molecular motors but can also function as modulators of intracellular signaling. KIF26A, an unconventional kinesin that lacks motor activity, inhibits growth-factor-receptor-bound protein 2 (GRB2)- and focal adhesion kinase (FAK)-dependent signal transduction, but its functions in the brain have not been characterized. We report a patient cohort with biallelic loss-of-function variants in KIF26A, exhibiting a spectrum of congenital brain malformations. In the developing brain, KIF26A is preferentially expressed during early- and mid-gestation in excitatory neurons. Combining mice and human iPSC-derived organoid models, we discovered that loss of KIF26A causes excitatory neuron-specific defects in radial migration, localization, dendritic and axonal growth, and apoptosis, offering a convincing explanation of the disease etiology in patients. Single-cell RNA sequencing in KIF26A knockout organoids revealed transcriptional changes in MAPK, MYC, and E2F pathways. Our findings illustrate the pathogenesis of KIF26A loss-of-function variants and identify the surprising versatility of this non-motor kinesin.

Early-Onset and Severe Complex Hereditary Spastic Paraplegia Caused by De Novo Variants in SPAST.

BACKGROUND: Familial hereditary spastic paraplegia (HSP)-SPAST (SPG4) typically presents with a pure HSP phenotype. OBJECTIVE: The aim of this study was to delineate the genotypic and phenotypic spectrum of children with de novo HSP-SPAST. METHODS: This study used a systematic cross-sectional analysis of clinical and molecular features. RESULTS: We report the clinical and molecular spectrum of 40 patients with heterozygous pathogenic de novo variants in SPAST (age range: 2.2-27.7 years). We identified 19 unique variants (16/40 carried the same recurrent variant, p.Arg499His). Symptom onset was in early childhood (median: 11.0 months, interquartile range: 6.0 months) with significant motor and speech delay, followed by progressive ascending spasticity, dystonia, neurogenic bladder dysfunction, gastrointestinal dysmotility, and epilepsy. The mean Spastic Paraplegia Rating Scale score was 32.8 ± 9.7 (standard deviation). CONCLUSIONS: These results confirm that de novo variants in SPAST lead to a severe and complex form of HSP that differs from classic familial pure HSP-SPAST. Clinicians should be aware of this syndrome in the differential diagnosis for cerebral palsy. © 2022 International Parkinson and Movement Disorder Society.

Long read sequencing and expression studies of AHDC1 deletions in Xia-Gibbs syndrome reveal a novel genetic regulatory mechanism.

Xia-Gibbs syndrome is a rare mendelian disorder characterized by Development Delay (DD), intellectual disability (ID), and hypotonia. Individuals with XGS typically harbor de novo protein-truncating mutations in the AT-Hook DNA binding motif containing 1 (AHDC1) gene, although some missense mutations can also cause XGS. Large de novo heterozygous deletions that encompass the AHDC1 gene have also been ascribed as diagnostic for the disorder, without substantial evidence to support their pathogenicity. We analyzed 19 individuals with large contiguous deletions involving AHDC1, along with other genes. One individual bore the smallest known contiguous AHDC1 deletion (∼350 Kb), encompassing eight other genes within chr1p36.11 (Feline Gardner-Rasheed, IFI6, FAM76A, STX12, PPP1R8, THEMIS2, RPA2, SMPDL3B) and terminating within the first intron of AHDC1. The breakpoint junctions and phase of the deletion were identified using both short and long read sequencing (Oxford Nanopore). Quantification of RNA expression patterns in whole blood revealed that AHDC1 exhibited a mono-allelic expression pattern with no deficiency in overall AHDC1 expression levels, in contrast to the other deleted genes, which exhibited a 50% reduction in mRNA expression. These results suggest that AHDC1 expression in this individual is compensated by a novel regulatory mechanism and advances understanding of mutational and regulatory mechanisms in neurodevelopmental disorders.

Functional characteristics of a broad spectrum of TBX6 variants in Mayer-Rokitansky-Küster-Hauser syndrome.

PURPOSE: Mayer-Rokitansky-Küster-Hauser syndrome (MRKHS) is characterized by congenital absence of the uterus, cervix, and the upper part of the vagina in females. Whole-gene deletion and loss-of-function variants in TBX6 have been identified in association with MRKHS. We aimed to expand the spectrum of TBX6 variants in MRKHS and explore the biological effect of the variant alleles. METHODS: Rare variants in TBX6 were called from a combined multiethnic cohort of 622 probands with MRKHS who underwent exome sequencing or genome sequencing. Multiple in vitro functional experiments were performed, including messenger RNA analysis, western blotting, transcriptional activity assay, and immunofluorescence staining. RESULTS: We identified 16 rare variants in TBX6 from the combined cohort, including 1 protein-truncating variant reported in our previous study and 15 variants with unknown effects. By comparing the prevalence of TBX6 variants in the Chinese MRKHS cohort vs 1038 female controls, we observed a significant mutational burden of TBX6 in affected individuals (P = .0004, odds ratio = 5.25), suggesting a causal role of TBX6 variants in MRKHS. Of the 15 variants with uncertain effects, 7 were shown to induce a loss-of-function effect through various mechanisms. The c.423G>A (p.Leu141=) and c.839+5G>A variants impaired the normal splicing of TBX6 messenger RNA, c.422T>C (p.Leu141Pro) and c.745G>A (p.Val249Met) led to decreased protein expression, c.10C>T (p.Pro4Ser) and c.400G>A (p.Glu134Lys) resulted in perturbed transcriptional activity, and c.356G>A (p.Arg119His) caused protein mislocalization. We observed incomplete penetrance and variable expressivity in families carrying deleterious variants, which indicates a more complex genetic mechanism than classical Mendelian inheritance. CONCLUSION: Our study expands the mutational spectrum of TBX6 in MRKHS and delineates the molecular pathogenesis of TBX6 variants, supporting the association between deleterious variants in TBX6 and MRKHS.

A reverse genetics and genomics approach to gene paralog function and disease: Myokymia and the juxtaparanode.

The leucine-rich glioma-inactivated (LGI) family consists of four highly conserved paralogous genes, LGI1-4, that are highly expressed in mammalian central and/or peripheral nervous systems. LGI1 antibodies are detected in subjects with autoimmune limbic encephalitis and peripheral nerve hyperexcitability syndromes (PNHSs) such as Isaacs and Morvan syndromes. Pathogenic variations of LGI1 and LGI4 are associated with neurological disorders as disease traits including familial temporal lobe epilepsy and neurogenic arthrogryposis multiplex congenita 1 with myelin defects, respectively. No human disease has been reported associated with either LGI2 or LGI3. We implemented exome sequencing and family-based genomics to identify individuals with deleterious variants in LGI3 and utilized GeneMatcher to connect practitioners and researchers worldwide to investigate the clinical and electrophysiological phenotype in affected subjects. We also generated Lgi3-null mice and performed peripheral nerve dissection and immunohistochemistry to examine the juxtaparanode LGI3 microarchitecture. As a result, we identified 16 individuals from eight unrelated families with loss-of-function (LoF) bi-allelic variants in LGI3. Deep phenotypic characterization showed LGI3 LoF causes a potentially clinically recognizable PNHS trait characterized by global developmental delay, intellectual disability, distal deformities with diminished reflexes, visible facial myokymia, and distinctive electromyographic features suggestive of motor nerve instability. Lgi3-null mice showed reduced and mis-localized Kv1 channel complexes in myelinated peripheral axons. Our data demonstrate bi-allelic LoF variants in LGI3 cause a clinically distinguishable disease trait of PNHS, most likely caused by disturbed Kv1 channel distribution in the absence of LGI3.

Biallelic variants in ADAMTS15 cause a novel form of distal arthrogryposis.

PURPOSE: We aimed to identify the underlying genetic cause for a novel form of distal arthrogryposis. METHODS: Rare variant family-based genomics, exome sequencing, and disease-specific panel sequencing were used to detect ADAMTS15 variants in affected individuals. Adamts15 expression was analyzed at the single-cell level during murine embryogenesis. Expression patterns were characterized using in situ hybridization and RNAscope. RESULTS: We identified homozygous rare variant alleles of ADAMTS15 in 5 affected individuals from 4 unrelated consanguineous families presenting with congenital flexion contractures of the interphalangeal joints and hypoplastic or absent palmar creases. Radiographic investigations showed physiological interphalangeal joint morphology. Additional features included knee, Achilles tendon, and toe contractures, spinal stiffness, scoliosis, and orthodontic abnormalities. Analysis of mouse whole-embryo single-cell sequencing data revealed a tightly regulated Adamts15 expression in the limb mesenchyme between embryonic stages E11.5 and E15.0. A perimuscular and peritendinous expression was evident in in situ hybridization in the developing mouse limb. In accordance, RNAscope analysis detected a significant coexpression with Osr1, but not with markers for skeletal muscle or joint formation. CONCLUSION: In aggregate, our findings provide evidence that rare biallelic recessive trait variants in ADAMTS15 cause a novel autosomal recessive connective tissue disorder, resulting in a distal arthrogryposis syndrome.

Developmental genomics of limb malformations: Allelic series in association with gene dosage effects contribute to the clinical variability.

Genetic heterogeneity, reduced penetrance, and variable expressivity, the latter including asymmetric body axis plane presentations, have all been described in families with congenital limb malformations (CLMs). Interfamilial and intrafamilial heterogeneity highlight the complexity of the underlying genetic pathogenesis of these developmental anomalies. Family-based genomics by exome sequencing (ES) and rare variant analyses combined with whole-genome array-based comparative genomic hybridization were implemented to investigate 18 families with limb birth defects. Eleven of 18 (61%) families revealed explanatory variants, including 7 single-nucleotide variant alleles and 3 copy number variants (CNVs), at previously reported "disease trait associated loci": BHLHA9, GLI3, HOXD cluster, HOXD13, NPR2, and WNT10B. Breakpoint junction analyses for all three CNV alleles revealed mutational signatures consistent with microhomology-mediated break-induced replication, a mechanism facilitated by Alu/Alu-mediated rearrangement. Homozygous duplication of BHLHA9 was observed in one Turkish kindred and represents a novel contributory genetic mechanism to Gollop-Wolfgang Complex (MIM: 228250), where triplication of the locus has been reported in one family from Japan (i.e., 4n = 2n + 2n versus 4n = 3n + 1n allelic configurations). Genes acting on limb patterning are sensitive to a gene dosage effect and are often associated with an allelic series. We extend an allele-specific gene dosage model to potentially assist, in an adjuvant way, interpretations of interconnections among an allelic series, clinical severity, and reduced penetrance of the BHLHA9-related CLM spectrum.

De novo heterozygous variants in SLC30A7 are a candidate cause for Joubert syndrome.

Joubert syndrome (JS), a well-established ciliopathy, is characterized by the distinctive molar tooth sign on brain MRI, ataxia, and neurodevelopmental features. Other manifestations can include polydactyly, accessory frenula, renal, or liver disease. Here, we report individuals meeting criteria for JS with de novo heterozygous variants in SLC30A7 (Chr1p21.2). The first individual is a female with history of unilateral postaxial polydactyly, classic molar tooth sign on MRI, macrocephaly, ataxia, ocular motor apraxia, neurodevelopmental delay, and precocious puberty. Exome sequencing detected a de novo heterozygous missense variant in SLC30A7: NM_133496.5: c.407 T > C, (p.Val136Ala). The second individual had bilateral postaxial polydactyly, molar tooth sign, macrocephaly, developmental delay, and an extra oral frenulum. A de novo deletion-insertion variant in SLC30A7, c.490_491delinsAG (p.His164Ser) was found. Both de novo variants affect highly conserved residues. Variants were not identified in known Joubert genes for either case. SLC30A7 has not yet been associated with a human phenotype. The SLC30 family of zinc transporters, like SLC30A7, permit cellular efflux of zinc, and although it is expressed in the brain its functions remain unknown. Published data from proteomic studies support SLC30A7 interaction with TCTN3, another protein associated with JS. The potential involvement of such genes in primary cilia suggest a role in Sonic Hedgehog signaling. SLC30A7 is a candidate JS-associated gene. Future work could be directed toward further characterization of SLC30A7 variants and understanding its function.

Biallelic pathogenic variants in roundabout guidance receptor 1 associate with syndromic congenital anomalies of the kidney and urinary tract.

Congenital anomalies of the kidney and urinary tract (CAKUT) represent the most common cause of chronic kidney failure in children. Despite growing knowledge of the genetic causes of CAKUT, the majority of cases remain etiologically unsolved. Genetic alterations in roundabout guidance receptor 1 (ROBO1) have been associated with neuronal and cardiac developmental defects in living individuals. Although Slit-Robo signaling is pivotal for kidney development, diagnostic ROBO1 variants have not been reported in viable CAKUT to date. By next-generation-sequencing methods, we identified six unrelated individuals and two non-viable fetuses with biallelic truncating or combined missense and truncating variants in ROBO1. Kidney and genitourinary manifestation included unilateral or bilateral kidney agenesis, vesicoureteral junction obstruction, vesicoureteral reflux, posterior urethral valve, genital malformation, and increased kidney echogenicity. Further clinical characteristics were remarkably heterogeneous, including neurodevelopmental defects, intellectual impairment, cerebral malformations, eye anomalies, and cardiac defects. By in silico analysis, we determined the functional significance of identified missense variants and observed absence of kidney ROBO1 expression in both human and murine mutant tissues. While its expression in multiple tissues may explain heterogeneous organ involvement, variability of the kidney disease suggests gene dosage effects due to a combination of null alleles with mild hypomorphic alleles. Thus, comprehensive genetic analysis in CAKUT should include ROBO1 as a new cause of recessively inherited disease. Hence, in patients with already established ROBO1-associated cardiac or neuronal disorders, screening for kidney involvement is indicated.