Exome sequencing in four families with neurodevelopmental disorders: genotype-phenotype correlation and identification of novel disease-causing variants in VPS13B and RELN.

Neurodevelopmental disorders (NDDs) are a clinically and genetically heterogeneous group of early-onset pediatric disorders that affect the structure and/or function of the central or peripheral nervous system. Achieving a precise molecular diagnosis for NDDs may be challenging due to the diverse genetic underpinnings and clinical variability. In the current study, we investigated the underlying genetic cause(s) of NDDs in four unrelated Pakistani families. Using exome sequencing (ES) as a diagnostic approach, we identified disease-causing variants in established NDD-associated genes in all families, including one hitherto unreported variant in RELN and three recurrent variants in VPS13B, DEGS1, and SPG11. Overall, our study highlights the potential of ES as a tool for clinical diagnosis.

TNF promoter hypomethylation is associated with mucosal inflammation in IBD and anti-TNF response.

Background and aims: Inflammatory Bowel Diseases (IBD) are chronic inflammatory conditions influenced heavily by environmental factors. DNA methylation is a form of epigenetic regulation linking environmental stimuli to gene expression changes and inflammation. Here, we investigated how DNA methylation of the TNF promoter differs between inflamed and uninflamed mucosa of IBD patients, including anti-TNF responders and non-responders. Methods: We obtained mucosal biopsies from 200 participants (133 IBD and 67 controls) and analyzed TNF promoter methylation using bisulfite sequencing, comparing inflamed with uninflamed segments, in addition to paired inflamed/uninflamed samples from individual patients. We conducted similar analyses on purified intestinal epithelial cells from bowel resections. We also compared TNF methylation levels of inflamed and uninflamed mucosa from a separate cohort of 15 anti-TNF responders and 17 non-responders. Finally, we sequenced DNA methyltransferase genes to identify rare variants in IBD patients and functionally tested them using rescue experiments in a zebrafish genetic model of DNA methylation deficiency. Results: TNF promoter methylation levels were decreased in inflamed mucosa of IBD patients and correlated with disease severity. Isolated IECs from inflamed tissue showed proportional decreases in TNF methylation. Anti-TNF non-responders showed lower levels of TNF methylation than responders in uninflamed mucosa. Our sequencing analysis revealed two missense variants in DNMT1, one of which had reduced function in vivo. Conclusions: Our study reveals an association of TNF promoter hypomethylation with mucosal inflammation, suggesting that IBD patients may be particularly sensitive to inflammatory environmental insults affecting DNA methylation. Together, our analyses indicate that TNF promoter methylation analysis may aid in the characterization of IBD status and evaluation of anti-TNF therapy response.

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

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

Exome sequencing in a Romanian Bardet-Biedl syndrome cohort revealed an overabundance of causal BBS12 variants.

Bardet-Biedl syndrome (BBS), is an emblematic ciliopathy hallmarked by pleiotropy, phenotype variability, and extensive genetic heterogeneity. BBS is a rare (~1/140,000 to ~1/160,000 in Europe) autosomal recessive pediatric disorder characterized by retinal degeneration, truncal obesity, polydactyly, cognitive impairment, renal dysfunction, and hypogonadism. Twenty-eight genes involved in ciliary structure or function have been implicated in BBS, and explain the molecular basis for ~75%-80% of individuals. To investigate the mutational spectrum of BBS in Romania, we ascertained a cohort of 24 individuals in 23 families. Following informed consent, we performed proband exome sequencing (ES). We detected 17 different putative disease-causing single nucleotide variants or small insertion-deletions and two pathogenic exon disruptive copy number variants in known BBS genes in 17 pedigrees. The most frequently impacted genes were BBS12 (35%), followed by BBS4, BBS7, and BBS10 (9% each) and BBS1, BBS2, and BBS5 (4% each). Homozygous BBS12 p.Arg355* variants were present in seven pedigrees of both Eastern European and Romani origin. Our data show that although the diagnostic rate of BBS in Romania is likely consistent with other worldwide cohorts (74%), we observed a unique distribution of causal BBS genes, including overrepresentation of BBS12 due to a recurrent nonsense variant, that has implications for regional diagnostics.

Association of deep phenotyping with diagnostic yield of prenatal exome sequencing for fetal brain abnormalities.

PURPOSE: To evaluate whether deep prenatal phenotyping of fetal brain abnormalities (FBAs) increases diagnostic yield of trio-exome sequencing (ES) compared with standard phenotyping. METHODS: Retrospective exploratory analysis of a multicenter prenatal ES study. Participants were eligible if an FBA was diagnosed and subsequently found to have a normal microarray. Deep phenotyping was defined as phenotype based on targeted ultrasound plus prenatal/postnatal magnetic resonance imaging, autopsy, and/or known phenotypes of other affected family members. Standard phenotyping was based on targeted ultrasound alone. FBAs were categorized by major brain findings on prenatal ultrasound. Cases with positive ES results were compared with those that have negative results by available phenotyping, as well as diagnosed FBAs. RESULTS: A total of 76 trios with FBAs were identified, of which 25 (33%) cases had positive ES results and 51 (67%) had negative results. Individual modalities of deep phenotyping were not associated with diagnostic ES results. The most common FBAs identified were posterior fossa anomalies and midline defects. Neural tube defects were significantly associated with receipt of a negative ES result (0% vs 22%, P = .01). CONCLUSION: Deep phenotyping was not associated with increased diagnostic yield of ES for FBA in this small cohort. Neural tube defects were associated with negative ES results.

Biallelic variants in ribonuclease inhibitor (RNH1), an inflammasome modulator, are associated with a distinctive subtype of acute, necrotizing encephalopathy.

PURPOSE: Mendelian etiologies for acute encephalopathies in previously healthy children are poorly understood, with the exception of RAN binding protein 2 (RANBP2)-associated acute necrotizing encephalopathy subtype 1 (ANE1). We provide clinical, genetic, and neuroradiological evidence that biallelic variants in ribonuclease inhibitor (RNH1) confer susceptibility to a distinctive ANE subtype. METHODS: This study aimed to evaluate clinical data, neuroradiological studies, genomic sequencing, and protein immunoblotting results in 8 children from 4 families who experienced acute febrile encephalopathy. RESULTS: All 8 healthy children became acutely encephalopathic during a viral/febrile illness and received a variety of immune modulation treatments. Long-term outcomes varied from death to severe neurologic deficits to normal outcomes. The neuroradiological findings overlapped with ANE but had distinguishing features. All affected children had biallelic predicted damaging variants in RNH1: a subset that was studied had undetectable RNH1 protein. Incomplete penetrance of the RNH1 variants was evident in 1 family. CONCLUSION: Biallelic variants in RNH1 confer susceptibility to a subtype of ANE (ANE2) in previously healthy children. Intensive immunological treatments may alter outcomes. Genomic sequencing in children with unexplained acute febrile encephalopathy can detect underlying genetic etiologies, such as RNH1, and improve outcomes in the probands and at-risk siblings.

Abrogation of MAP4K4 protein function causes congenital anomalies in humans and zebrafish.

We report 21 families displaying neurodevelopmental differences and multiple congenital anomalies while bearing a series of rare variants in mitogen-activated protein kinase kinase kinase kinase 4 (MAP4K4). MAP4K4 has been implicated in many signaling pathways including c-Jun N-terminal and RAS kinases and is currently under investigation as a druggable target for multiple disorders. Using several zebrafish models, we demonstrate that these human variants are either loss-of-function or dominant-negative alleles and show that decreasing Map4k4 activity causes developmental defects. Furthermore, MAP4K4 can restrain hyperactive RAS signaling in early embryonic stages. Together, our data demonstrate that MAP4K4 negatively regulates RAS signaling in the early embryo and that variants identified in affected humans abrogate its function, establishing MAP4K4 as a causal locus for individuals with syndromic neurodevelopmental differences.

Diagnosis of TBC1D32-associated conditions: Expanding the phenotypic spectrum of a complex ciliopathy.

Exome sequencing is a powerful tool in prenatal and postnatal genetics and can help identify novel candidate genes critical to human development. We describe seven unpublished probands with rare likely pathogenic variants or variants of uncertain significance that segregate with recessive disease in TBC1D32, including four fetal probands in three unrelated pedigrees and three pediatric probands in unrelated pedigrees. We also report clinical comparisons with seven previously published patients. Index probands were identified through an ongoing prenatal exome sequencing study and through an online data sharing platform (Gene Matcher™). A literature review was also completed. TBC1D32 is involved in the development and function of cilia and is expressed in the developing hypothalamus and pituitary gland. We provide additional data to expand the phenotype correlated with TBC1D32 variants, including a severe prenatal phenotype associated with life-limiting congenital anomalies.

Pathogenic variants in SLF2 and SMC5 cause segmented chromosomes and mosaic variegated hyperploidy.

Embryonic development is dictated by tight regulation of DNA replication, cell division and differentiation. Mutations in DNA repair and replication genes disrupt this equilibrium, giving rise to neurodevelopmental disease characterized by microcephaly, short stature and chromosomal breakage. Here, we identify biallelic variants in two components of the RAD18-SLF1/2-SMC5/6 genome stability pathway, SLF2 and SMC5, in 11 patients with microcephaly, short stature, cardiac abnormalities and anemia. Patient-derived cells exhibit a unique chromosomal instability phenotype consisting of segmented and dicentric chromosomes with mosaic variegated hyperploidy. To signify the importance of these segmented chromosomes, we have named this disorder Atelís (meaning - incomplete) Syndrome. Analysis of Atelís Syndrome cells reveals elevated levels of replication stress, partly due to a reduced ability to replicate through G-quadruplex DNA structures, and also loss of sister chromatid cohesion. Together, these data strengthen the functional link between SLF2 and the SMC5/6 complex, highlighting a distinct role for this pathway in maintaining genome stability.

The p190 RhoGAPs, ARHGAP35, and ARHGAP5 are implicated in GnRH neuronal development: Evidence from patients with idiopathic hypogonadotropic hypogonadism, zebrafish, and in vitro GAP activity assay.

PURPOSE: The study aimed to identify novel genes for idiopathic hypogonadotropic hypogonadism (IHH). METHODS: A cohort of 1387 probands with IHH underwent exome sequencing and de novo, familial, and cohort-wide investigations. Functional studies were performed on 2 p190 Rho GTPase-activating proteins (p190 RhoGAP), ARHGAP35 and ARHGAP5, which involved in vivo modeling in larval zebrafish and an in vitro p190A-GAP activity assay. RESULTS: Rare protein-truncating variants (PTVs; n = 5) and missense variants in the RhoGAP domain (n = 7) in ARHGAP35 were identified in IHH cases (rare variant enrichment: PTV [unadjusted P = 3.1E-06] and missense [adjusted P = 4.9E-03] vs controls). Zebrafish modeling using gnrh3:egfp phenotype assessment showed that mutant larvae with deficient arhgap35a, the predominant ARHGAP35 paralog in the zebrafish brain, display decreased GnRH3-GFP+ neuronal area, a readout for IHH. In vitro GAP activity studies showed that 1 rare missense variant [ARHGAP35 p.(Arg1284Trp)] had decreased GAP activity. Rare PTVs (n = 2) also were discovered in ARHGAP5, a paralog of ARHGAP35; however, arhgap5 zebrafish mutants did not display significant GnRH3-GFP+ abnormalities. CONCLUSION: This study identified ARHGAP35 as a new autosomal dominant genetic driver for IHH and ARHGAP5 as a candidate gene for IHH. These observations suggest a novel role for the p190 RhoGAP proteins in GnRH neuronal development and integrity.

A cross-disorder dosage sensitivity map of the human genome.

Rare copy-number variants (rCNVs) include deletions and duplications that occur infrequently in the global human population and can confer substantial risk for disease. In this study, we aimed to quantify the properties of haploinsufficiency (i.e., deletion intolerance) and triplosensitivity (i.e., duplication intolerance) throughout the human genome. We harmonized and meta-analyzed rCNVs from nearly one million individuals to construct a genome-wide catalog of dosage sensitivity across 54 disorders, which defined 163 dosage sensitive segments associated with at least one disorder. These segments were typically gene dense and often harbored dominant dosage sensitive driver genes, which we were able to prioritize using statistical fine-mapping. Finally, we designed an ensemble machine-learning model to predict probabilities of dosage sensitivity (pHaplo & pTriplo) for all autosomal genes, which identified 2,987 haploinsufficient and 1,559 triplosensitive genes, including 648 that were uniquely triplosensitive. This dosage sensitivity resource will provide broad utility for human disease research and clinical genetics.

Testing for rare genetic causes of obesity: findings and experiences from a pediatric weight management program.

BACKGROUND: Genetic screening for youth with obesity in the absence of syndromic findings has not been part of obesity management. For children with early onset obesity, genetic screening is recommended for those having clinical features of genetic obesity syndromes (including hyperphagia). OBJECTIVES: The overarching goal of this work is to report the findings and experiences from one pediatric weight management program that implemented targeted sequencing analysis for genes known to cause rare genetic disorders of obesity. SUBJECTS/METHODS: This exploratory study evaluated youth tested over an 18-month period using a panel of 40-genes in the melanocortin 4 receptor pathway. Medical records were reviewed for demographic and visit information, including body mass index (BMI) percent of 95th percentile (%BMIp95) and two eating behaviors. RESULTS: Of 117 subjects: 51.3% were male; 53.8% Hispanic; mean age 10.2 years (SD 3.8); mean %BMIp95 157% (SD 29%). Most subjects were self- or caregiver-reported to have overeating to excess or binge eating (80.3%) and sneaking food or eating in secret (59.0%). Among analyzed genes, 72 subjects (61.5%) had at least one variant reported; 50 (42.7%) had a single variant reported; 22 (18.8%) had 2-4 variants reported; most variants were rare (<0.05% minor allele frequency [MAF]), and of uncertain significance; all variants were heterozygous. Nine subjects (7.7%) had a variant reported as PSCK1 "risk" or MC4R "likely pathogenic"; 39 (33.3%) had a Bardet-Biedl Syndrome (BBS) gene variant (4 with "pathogenic" or "likely pathogenic" variants). Therefore, 9 youth (7.7%) had gene variants previously identified as increasing risk for obesity and 4 youth (3.4%) had BBS carrier status. CONCLUSIONS: Panel testing identified rare variants of uncertain significance in most youth tested, and infrequently identified variants previously reported to increase the risk for obesity. Further research in larger cohorts is needed to understand how genetic variants influence the expression of non-syndromic obesity.

A recurrent rare intronic variant in CAPN3 alters mRNA splicing and causes autosomal recessive limb-girdle muscular dystrophy-1 in three Pakistani pedigrees.

Autosomal recessive limb-girdle muscular dystrophy-1 (LGMDR1) is an autosomal recessive disorder characterized by progressive weakness of the proximal limb and girdle muscles. Biallelic mutations in CAPN3 are reported frequently to cause LGMDR1. Here, we describe 11 individuals from three unrelated consanguineous families that present with typical features of LGMDR1 that include proximal muscle wasting, weakness of the upper and lower limbs, and elevated serum creatine kinase. Whole-exome sequencing identified a rare homozygous CAPN3 variant near the exon 2 splice donor site that segregates with disease in all three families. mRNA splicing studies showed partial retention of intronic sequence and subsequent introduction of a premature stop codon (NM_000070.3: c.379 + 3A>G; p.Asp128Glyfs*15). Furthermore, we observe reduced CAPN3 expression in primary dermal fibroblasts derived from an affected individual, suggesting instability and/or nonsense-mediated decay of mutation-bearing mRNA. Genome-wide homozygosity mapping and single-nucleotide polymorphism analysis identified a shared haplotype and supports a possible founder effect for the CAPN3 variant. Together, our data extend the mutational spectrum of LGMDR1 and have implications for improved diagnostics for individuals of Pakistani origin.

Multidisciplinary approaches for elucidating genetics and molecular pathogenesis of urinary tract malformations.

Advances in clinical diagnostics and molecular tools have improved our understanding of the genetically heterogeneous causes underlying congenital anomalies of kidney and urinary tract (CAKUT). However, despite a sharp incline of CAKUT reports in the literature within the past 2 decades, there remains a plateau in the genetic diagnostic yield that is disproportionate to the accelerated ability to generate robust genome-wide data. Explanations for this observation include (i) diverse inheritance patterns with incomplete penetrance and variable expressivity, (ii) rarity of single-gene drivers such that large sample sizes are required to meet the burden of proof, and (iii) multigene interactions that might produce either intra- (e.g., copy number variants) or inter- (e.g., effects in trans) locus effects. These challenges present an opportunity for the community to implement innovative genetic and molecular avenues to explain the missing heritability and to better elucidate the mechanisms that underscore CAKUT. Here, we review recent multidisciplinary approaches at the intersection of genetics, genomics, in vivo modeling, and in vitro systems toward refining a blueprint for overcoming the diagnostic hurdles that are pervasive in urinary tract malformation cohorts. These approaches will not only benefit clinical management by reducing age at molecular diagnosis and prompting early evaluation for comorbid features but will also serve as a springboard for therapeutic development.

A case of Bardet-Biedl syndrome caused by a recurrent variant in BBS12: A case report.

Bardet-Biedl syndrome (BBS) is a clinically and genetically heterogenous disorder that manifests as a result of primary cilia impairment. Cilia are present on most cell types, thus BBS is a multisystemic condition involving the majority of organ systems. The core features of the syndrome include retinal degeneration, obesity, polydactyly, cognitive impairment, renal anomalies and urogenital malformations. To date, pathogenic variants in 26 genes have been shown to be involved in the molecular basis of this rare ciliopathy. Of these causal loci, BBS12 accounts for ~8% of all cases. In this case report, an individual with BBS caused by a rare recurrent variant in BBS12 (NM_152618.3: c.1063C>T; p.Arg355*) is described and compared with others with the same DNA variant, placing this finding in the context of the current literature.

A recessive variant in TFAM causes mtDNA depletion associated with primary ovarian insufficiency, seizures, intellectual disability and hearing loss.

Mitochondrial disorders are collectively common, genetically heterogeneous disorders in both pediatric and adult populations. They are caused by molecular defects in oxidative phosphorylation, failure of essential bioenergetic supply to mitochondria, and apoptosis. Here, we present three affected individuals from a consanguineous family of Pakistani origin with variable seizures and intellectual disability. Both females display primary ovarian insufficiency (POI), while the male shows abnormal sex hormone levels. We performed whole exome sequencing and identified a recessive missense variant c.694C > T, p.Arg232Cys in TFAM that segregates with disease. TFAM (mitochondrial transcription factor A) is a component of the mitochondrial replisome machinery that maintains mtDNA transcription and replication. In primary dermal fibroblasts, we show depletion of mtDNA and significantly altered mitochondrial function and morphology. Moreover, we observed reduced nucleoid numbers with significant changes in nucleoid size or shape in fibroblasts from an affected individual compared to controls. We also investigated the effect of tfam impairment in zebrafish; homozygous tfam mutants carrying an in-frame c.141_149 deletion recapitulate the mtDNA depletion and ovarian dysgenesis phenotypes observed in affected humans. Together, our genetic and functional data confirm that TFAM plays a pivotal role in gonad development and expands the repertoire of mitochondrial disease phenotypes.

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

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

Acoustofluidic rotational tweezing enables high-speed contactless morphological phenotyping of zebrafish larvae.

Modern biomedical research and preclinical pharmaceutical development rely heavily on the phenotyping of small vertebrate models for various diseases prior to human testing. In this article, we demonstrate an acoustofluidic rotational tweezing platform that enables contactless, high-speed, 3D multispectral imaging and digital reconstruction of zebrafish larvae for quantitative phenotypic analysis. The acoustic-induced polarized vortex streaming achieves contactless and rapid (~1 s/rotation) rotation of zebrafish larvae. This enables multispectral imaging of the zebrafish body and internal organs from different viewing perspectives. Moreover, we develop a 3D reconstruction pipeline that yields accurate 3D models based on the multi-view images for quantitative evaluation of basic morphological characteristics and advanced combinations of metrics. With its contactless nature and advantages in speed and automation, our acoustofluidic rotational tweezing system has the potential to be a valuable asset in numerous fields, especially for developmental biology, small molecule screening in biochemistry, and pre-clinical drug development in pharmacology.

HEK293T Cells with TFAM Disruption by CRISPR-Cas9 as a Model for Mitochondrial Regulation.

The mitochondrial transcription factor A (TFAM) is considered a key factor in mitochondrial DNA (mtDNA) copy number. Given that the regulation of active copies of mtDNA is still not fully understood, we investigated the effects of CRISPR-Cas9 gene editing of TFAM in human embryonic kidney (HEK) 293T cells on mtDNA copy number. The aim of this study was to generate a new in vitro model by CRISPR-Cas9 system by editing the TFAM locus in HEK293T cells. Among the resulting single-cell clones, seven had high mutation rates (67-96%) and showed a decrease in mtDNA copy number compared to control. Cell staining with Mitotracker Red showed a reduction in fluorescence in the edited cells compared to the non-edited cells. Our findings suggest that the mtDNA copy number is directly related to TFAM control and its disruption results in interference with mitochondrial stability and maintenance.

A BBS1 SVA F retrotransposon insertion is a frequent cause of Bardet-Biedl syndrome.

Bardet-Biedl syndrome (BBS) is a ciliopathy characterized by retinitis pigmentosa, obesity, polydactyly, cognitive impairment and renal failure. Pathogenic variants in 24 genes account for the molecular basis of >80% of cases. Toward saturated discovery of the mutational basis of the disorder, we carefully explored our cohorts and identified a hominid-specific SINE-R/VNTR/Alu type F (SVA-F) insertion in exon 13 of BBS1 in eight families. In six families, the repeat insertion was found in trans with c.1169 T > G, p.Met390Arg and in two families the insertion was found in addition to other recessive BBS loci. Whole genome sequencing, de novo assembly and SNP array analysis were performed to characterize the genomic event. This insertion is extremely rare in the general population (found in 8 alleles of 8 BBS cases but not in >10 800 control individuals from gnomAD-SV) and due to a founder effect. Its 2435 bp sequence contains hallmarks of LINE1 mediated retrotransposition. Functional studies with patient-derived cell lines confirmed that the BBS1 SVA-F is deleterious as evidenced by a significant depletion of both mRNA and protein levels. Such findings highlight the importance of dedicated bioinformatics pipelines to identify all types of variation.