Optical genome mapping unveils hidden structural variants in neurodevelopmental disorders.

While short-read sequencing currently dominates genetic research and diagnostics, it frequently falls short of capturing certain structural variants (SVs), which are often implicated in the etiology of neurodevelopmental disorders (NDDs). Optical genome mapping (OGM) is an innovative technique capable of capturing SVs that are undetectable or challenging-to-detect via short-read methods. This study aimed to investigate NDDs using OGM, specifically focusing on cases that remained unsolved after standard exome sequencing. OGM was performed in 47 families using ultra-high molecular weight DNA. Single-molecule maps were assembled de novo, followed by SV and copy number variant calling. We identified 7 variants of interest, of which 5 (10.6%) were classified as likely pathogenic or pathogenic, located in BCL11A, OPHN1, PHF8, SON, and NFIA. We also identified an inversion disrupting NAALADL2, a gene which previously was found to harbor complex rearrangements in two NDD cases. Variants in known NDD genes or candidate variants of interest missed by exome sequencing mainly consisted of larger insertions (> 1kbp), inversions, and deletions/duplications of a low number of exons (1-4 exons). In conclusion, in addition to improving molecular diagnosis in NDDs, this technique may also reveal novel NDD genes which may harbor complex SVs often missed by standard sequencing techniques.

Heterogeneous genetic patterns in bilateral perisylvian polymicrogyria: insights from a Finnish family cohort.

Bilateral perisylvian polymicrogyria is the most common form of regional polymicrogyria within malformations of cortical development, constituting 20% of all malformations of cortical development. Bilateral perisylvian polymicrogyria is characterized by an excessive folding of the cerebral cortex and abnormal cortical layering. Notable clinical features include upper motoneuron dysfunction, dysarthria and asymmetric quadriparesis. Cognitive impairment and epilepsy are frequently observed. To identify genetic variants underlying bilateral perisylvian polymicrogyria in Finland, we examined 21 families using standard exome sequencing, complemented by optical genome mapping and/or deep exome sequencing. Pathogenic or likely pathogenic variants were identified in 5/21 (24%) of families, of which all were confirmed as de novo. These variants were identified in five genes, i.e. DDX23, NUS1, SCN3A, TUBA1A and TUBB2B, with NUS1 and DDX23 being associated with bilateral perisylvian polymicrogyria for the first time. In conclusion, our results confirm the previously reported genetic heterogeneity of bilateral perisylvian polymicrogyria and underscore the necessity of more advanced methods to elucidate the genetic background of bilateral perisylvian polymicrogyria.

A novel variant in CYFIP2 in a girl with severe disabilities and bilateral perisylvian polymicrogyria.

Bilateral perisylvian polymicrogyria (BPP) is a structural malformation of the cerebral cortex that can be caused by several genetic abnormalities. The most common clinical manifestations of BPP include intellectual disability and epilepsy. Cytoplasmic FMRP-interacting protein 2 (CYFIP2) is a protein that interacts with the fragile X mental retardation protein (FMRP). CYFIP2 variants can cause various brain structural abnormalities with the most common clinical manifestations of intellectual disability, epileptic encephalopathy and dysmorphic features. We present a girl with multiple disabilities and BPP caused by a heterozygous, novel, likely pathogenic variant (c.1651G>C: p.(Val551Leu) in the CYFIP2 gene. Our case report broadens the spectrum of genetic diversity associated with BPP by incorporating CYFIP2.

A loss of function variant in AGPAT3 underlies intellectual disability and retinitis pigmentosa (IDRP) syndrome.

Intellectual disability (ID) and retinal dystrophy (RD) are the frequently found features of multiple syndromes involving additional systemic manifestations. Here, we studied a family with four members presenting severe ID and retinitis pigmentosa (RP). Using genome wide genotyping and exome sequencing, we identified a nonsense variant c.747 C > A (p.Tyr249Ter) in exon 7 of AGPAT3 which co-segregates with the disease phenotype. Western blot analysis of overexpressed WT and mutant AGPAT3 in HEK293T cells showed the absence of AGPAT3, suggesting instability of the truncated protein. Knockdown of Agpat3 in the embryonic mouse brain caused marked deficits in neuronal migration, strongly suggesting that reduced expression of AGPAT3 affects neuronal function. Altogether, our data indicates that AGPAT3 activity is essential for neuronal functioning and loss of its activity probably causes intellectual disability and retinitis pigmentosa (IDRP) syndrome.

Variants in EFCAB7 underlie nonsyndromic postaxial polydactyly.

Polydactyly is the most common limb malformation that occurs in 1.6-10.6 per one thousand live births, with incidence varying with ancestry. The underlying gene has been identified for many of the ~100 syndromes that include polydactyly. While for the more common form, nonsydromic polydactyly, eleven candidate genes have been reported. We investigated the underlying genetic cause of autosomal recessive nonsyndromic postaxial polydactyly in four consanguineous Pakistani families. Some family members with postaxial polydactyly also present with syndactyly, camptodactyly, or clinodactyly. Analysis of the exome sequence data revealed two novel homozygous frameshift deletions in EFCAB7: [c.830delG;p.(Gly277Valfs*5)]; in three families and [c.1350_1351delGA;p.(Asn451Phefs*2)] in one family. Sanger sequencing confirmed that these variants segregated with postaxial polydactyly, i.e., family members with postaxial polydactyly were found to be homozygous while unaffected members were heterozygous or wild type. EFCAB7 displays expressions in the skeletal muscle and on the cellular level in cilia. IQCE-EFCAB7 and EVC-EVC2 are part of the heterotetramer EvC complex, which is a positive regulator of the Hedgehog (Hh) pathway, that plays a key role in limb formation. Depletion of either EFCAB7 or IQCE inhibits induction of Gli1, a direct Hh target gene. Variants in IQCE and GLI1 have been shown to cause nonsyndromic postaxial polydactyly, while variants in EVC and EVC2 underlie Ellis van Creveld and Weyers syndromes, which include postaxial polydactyly as a phenotype. This is the first report of the involvement of EFCAB7 in human disease etiology.

Comparison of multiple imputation and other methods for the analysis of imputed genotypes.

BACKGROUND: Analysis of imputed genotypes is an important and routine component of genome-wide association studies and the increasing size of imputation reference panels has facilitated the ability to impute and test low-frequency variants for associations. In the context of genotype imputation, the true genotype is unknown and genotypes are inferred with uncertainty using statistical models. Here, we present a novel method for integrating imputation uncertainty into statistical association tests using a fully conditional multiple imputation (MI) approach which is implemented using the Substantive Model Compatible Fully Conditional Specification (SMCFCS). We compared the performance of this method to an unconditional MI and two additional approaches that have been shown to demonstrate excellent performance: regression with dosages and a mixture of regression models (MRM). RESULTS: Our simulations considered a range of allele frequencies and imputation qualities based on data from the UK Biobank. We found that the unconditional MI was computationally costly and overly conservative across a wide range of settings. Analyzing data with Dosage, MRM, or MI SMCFCS resulted in greater power, including for low frequency variants, compared to unconditional MI while effectively controlling type I error rates. MRM andl MI SMCFCS are both more computationally intensive then using Dosage. CONCLUSIONS: The unconditional MI approach for association testing is overly conservative and we do not recommend its use in the context of imputed genotypes. Given its performance, speed, and ease of implementation, we recommend using Dosage for imputed genotypes with MAF [Formula: see text] 0.001 and Rsq [Formula: see text] 0.3.

Variants in JAZF1 are associated with asthma, type 2 diabetes, and height in the United Kingdom biobank population.

Background: Asthma, type 2 diabetes (T2D), and anthropometric measures are correlated complex traits that all have a major genetic component. Objective: To investigate the overlap in genetic variants associated with these complex traits. Methods: Using United Kingdom Biobank data, we performed univariate association analysis, fine-mapping, and mediation analysis to identify and dissect shared genomic regions associated with asthma, T2D, height, weight, body mass index (BMI), and waist circumference (WC). Results: We found several genome-wide significant variants in and around the JAZF1 gene that are associated with asthma, T2D, or height with two of these variants shared by the three phenotypes. We also observed an association in this region with WC when adjusted for BMI. However, there was no association with WC when it was not adjusted for BMI or weight. Additionally, only suggestive associations between variants in this region and BMI were observed. Fine-mapping analyses suggested that within JAZF1 there are non-overlapping regions harboring causal susceptibility variants for asthma, T2D, and height. Mediation analyses supported the conclusion that these are independent associations. Conclusion: Our findings indicate that variants in the JAZF1 are associated with asthma, T2D, and height, but the associated causal variant(s) are different for each of the three phenotypes.

Nonsense variant in a consanguineous family expands the phenotype of KPTN gene-related syndrome to include hearing impairment.

A short report with two affected siblings from consanguineous family born with intellectual disability, motor disability, language deficit, and hearing impairment and found to carry biallelic nonsense variant in KPTN gene known to be associated with KPTN gene related syndrome.

A Novel Variant in VPS13B Underlying Cohen Syndrome.

Pathogenic variants in vacuolar protein sorting 13 homolog B (VPS13B) cause Cohen syndrome (CS), a clinically diverse neurodevelopmental disorder. We used whole exome and Sanger sequencing to identify disease-causing variants in a Pakistani family with intellectual disability, microcephaly, facial dysmorphism, neutropenia, truncal obesity, speech delay, motor delay, and insomnia. We identified a novel homozygous nonsense variant c.8841G > A: p.(W2947∗) in VPS13B (NM_017890.5) which segregated with the disease. Sleep disturbances are commonly seen in neurodevelopmental disorders and can exacerbate medical issues if left untreated. We demonstrate that individuals with Cohen syndrome may also be affected by sleep disturbances. In conclusion, we expand the genetic and phenotypic features of Cohen syndrome in the Pakistani population.

The genetic contribution of the X chromosome in age-related hearing loss.

Age-related (AR) hearing loss (HL) is the most common sensory impairment with heritability of 55%. The aim of this study was to identify genetic variants on chromosome X associated with ARHL through the analysis of data obtained from the UK Biobank. We performed association analysis between self-reported measures of HL and genotyped and imputed variants on chromosome X from ∼460,000 white Europeans. We identified three loci associated with ARHL with a genome-wide significance level (p < 5 × 10-8), ZNF185 (rs186256023, p = 4.9 × 10-10) and MAP7D2 (rs4370706, p = 2.3 × 10-8) in combined analysis of males and females, and LOC101928437 (rs138497700, p = 8.9 × 10-9) in the sex-stratified analysis of males. In-silico mRNA expression analysis showed MAP7D2 and ZNF185 are expressed in mice and adult human inner ear tissues, particularly in the inner hair cells. We estimated that only a small amount of variation of ARHL, 0.4%, is explained by variants on the X chromosome. This study suggests that although there are likely a few genes contributing to ARHL on the X chromosome, the role that the X chromosome plays in the etiology of ARHL may be limited.

An FBN1 deep intronic variant is associated with pseudoexon formation and a variable Marfan phenotype in a five generation family.

Exome sequencing of genes associated with heritable thoracic aortic disease (HTAD) failed to identify a pathogenic variant in a large family with Marfan syndrome (MFS). A genome-wide linkage analysis for thoracic aortic disease identified a peak at 15q21.1, and genome sequencing identified a novel deep intronic FBN1 variant that segregated with thoracic aortic disease in the family (LOD score 2.7) and was predicted to alter splicing. RT-PCR and bulk RNA sequencing of RNA harvested from fibroblasts explanted from the affected proband revealed an insertion of a pseudoexon between exons 13 and 14 of the FBN1 transcript, predicted to lead to nonsense mediated decay (NMD). Treating the fibroblasts with an NMD inhibitor, cycloheximide, greatly improved the detection of the pseudoexon-containing transcript. Family members with the FBN1 variant had later onset aortic events and fewer MFS systemic features than typical for individuals with haploinsufficiency of FBN1. Variable penetrance of the phenotype and negative genetic testing in MFS families should raise the possibility of deep intronic FBN1 variants and the need for additional molecular studies.

Rare-variant association analysis reveals known and new age-related hearing loss genes.

Age-related (AR) hearing loss (HL) is a prevalent sensory deficit in the elderly population. Several studies showed that common variants increase ARHL susceptibility. Here, we demonstrate that rare-variants play a crucial role in ARHL etiology. We analyzed exome and imputed data from white-European UK Biobank volunteers, performing both single-variant and rare-variant aggregate association analyses using self-reported ARHL phenotypes. We identified and replicated associations between ARHL and rare-variants in KLHDC7B, PDCD6, MYO6, SYNJ2, and TECTA. PUS7L and EYA4 also revealed rare-variant associations with ARHL. EYA4, MYO6, and TECTA are all known to underline Mendelian nonsyndromic HL. PDCD6, a new HL gene, plays an important role in apoptosis and has widespread inner ear expression, particularly in the inner hair cells. An unreplicated common variant association was previously observed for KHLDC7B, here we demonstrate that rare-variants in this gene also play a role in ARHL etiology. Additionally, the first replicated association between SYNJ2 and ARHL was detected. Analysis of common variants revealed several previously reported, i.e., ARHGEF28, and new, i.e., PIK3R3, ARHL associations, as well as ones we replicate here for the first time, i.e., BAIAP2L2, CRIP3, KLHDC7B, MAST2, and SLC22A7. It was also observed that the odds ratios for rare-variant ARHL associations, were higher than those for common variants. In conclusion, we demonstrate the vital role rare-variants, including those in Mendelian nonsyndromic HL genes, play in the etiology of ARHL.

Variants of LRP2, encoding a multifunctional cell-surface endocytic receptor, associated with hearing loss and retinal dystrophy.

Hereditary deafness and retinal dystrophy are each genetically heterogenous and clinically variable. Three small unrelated families segregating the combination of deafness and retinal dystrophy were studied by exome sequencing (ES). The proband of Family 1 was found to be compound heterozygous for NM_004525.3: LRP2: c.5005A > G, p.(Asn1669Asp) and c.149C > G, p.(Thr50Ser). In Family 2, two sisters were found to be compound heterozygous for LRP2 variants, p.(Tyr3933Cys) and an experimentally confirmed c.7715 + 3A > T consensus splice-altering variant. In Family 3, the proband is compound heterozygous for a consensus donor splice site variant LRP2: c.8452_8452 + 1del and p.(Cys3150Tyr). In mouse cochlea, Lrp2 is expressed abundantly in the stria vascularis marginal cells demonstrated by smFISH, single-cell and single-nucleus RNAseq, suggesting that a deficiency of LRP2 may compromise the endocochlear potential, which is required for hearing. LRP2 variants have been associated with Donnai-Barrow syndrome and other multisystem pleiotropic phenotypes different from the phenotypes of the four cases reported herein. Our data expand the phenotypic spectrum associated with pathogenic variants in LRP2 warranting their consideration in individuals with a combination of hereditary hearing loss and retinal dystrophy.

Syntaxin 4 is essential for hearing in human and zebrafish.

Congenital hearing impairment (HI) is a genetically highly heterogeneous disorder in which prompt recognition and intervention are crucial to optimize outcomes. In this study, we used exome sequencing to investigate a large consanguineous Pakistani family with eight affected individuals showing bilateral severe-to-profound HI. This identified a homozygous splice region variant in STX4 (c.232 + 6T>C), which causes exon skipping and a frameshift, that segregated with HI (two-point logarithm of odds (LOD) score = 5.9). STX4, a member of the syntaxin family, is a component of the SNARE machinery involved in several vesicle transport and recycling pathways. In silico analysis showed that murine orthologue Stx4a is highly and widespread expressed in the developing and adult inner ear. Immunofluorescent imaging revealed localization of STX4A in the cell body, cell membrane and stereocilia of inner and outer hair cells. Furthermore, a morpholino-based knockdown of stx4 in zebrafish showed an abnormal startle response, morphological and developmental defects, and a disrupted mechanotransduction function in neuromast hair cells measured via FM1-43 uptake. Our findings indicate that STX4 dysfunction leads to HI in humans and zebrafish and supports the evolutionary conserved role of STX4 in inner ear development and hair cell functioning.

Clinical, muscle imaging, and genetic characteristics of dystrophinopathies with deep-intronic DMD variants.

BACKGROUND: Phenotypic heterogeneity within or between families with a same deep-intronic splice-altering variant in the DMD gene has never been systematically analyzed. This study aimed to determine the phenotypic and genetic characteristics of patients with deep-intronic DMD variants. METHODS: Of 1338 male patients with a suspected dystrophinopathy, 38 were confirmed to have atypical pathogenic DMD variants via our comprehensive genetic testing approach. Of the 38 patients, 30 patients from 22 unrelated families with deep-intronic DMD variants underwent a detailed clinical and imaging assessment. RESULTS: Nineteen different deep-intronic DMD variants were identified in the 30 patients, including 15 with Duchenne muscular dystrophy (DMD), 14 with Becker muscular dystrophy (BMD), and one with X-linked dilated cardiomyopathy. Of the 19 variants, 15 were single-nucleotide variants, 2 were structural variants (SVs), and 2 were pure-intronic large-scale SVs causing aberrant inclusion of other protein-coding genes sequences into the mature DMD transcripts. The trefoil with single fruit sign was observed in 18 patients and the concentric fatty infiltration pattern was observed in 2 patients. Remarkable phenotypic heterogeneity was observed not only in skeletal but also cardiac muscle involvement in 2 families harboring a same deep-intronic variant. Different skeletal muscle involvement between families with a same variant was observed in 4 families. High inter-individual phenotypic heterogeneity was observed within two BMD families and one DMD family. CONCLUSIONS: Our study first highlights the variable phenotypic expressivity of deep-intronic DMD variants and demonstrates a new class of deep-intronic DMD variants, i.e., pure-intronic SVs involving other protein-coding genes.

A novel autosomal dominant GREB1L variant associated with non-syndromic hearing impairment in Ghana.

BACKGROUND: Childhood hearing impairment (HI) is genetically heterogeneous with many implicated genes, however, only a few of these genes are reported in African populations. METHODS: This study used exome and Sanger sequencing to resolve the possible genetic cause of non-syndromic HI in a Ghanaian family. RESULTS: We identified a novel variant c.3041G > A: p.(Gly1014Glu) in GREB1L (DFNA80) in the index case. The GREB1L: p.(Gly1014Glu) variant had a CADD score of 26.5 and was absent from human genomic databases such as TopMed and gnomAD. In silico homology protein modeling approaches displayed major structural differences between the wildtype and mutant proteins. Additionally, the variant was predicted to probably affect the secondary protein structure that may impact its function. Publicly available expression data shows a higher expression of Greb1L in the inner ear of mice during development and a reduced expression in adulthood, underscoring its importance in the development of the inner ear structures. CONCLUSION: This report on an African individual supports the association of GREB1L variant with non-syndromic HI and extended the evidence of the implication of GREB1L variants in HI in diverse populations.

Whole genome sequencing identifies candidate genes for familial essential tremor and reveals biological pathways implicated in essential tremor aetiology.

BACKGROUND: Essential tremor (ET), one of the most common neurological disorders, has a phenotypically heterogeneous presentation characterized by bilateral kinetic tremor of the arms and, in some patients, tremor involving other body regions (e.g., head, voice). Genetic studies suggest that ET is genetically heterogeneous. METHODS: We analyzed whole genome sequence data (WGS) generated on 104 multi-generational white families with European ancestry affected by ET. Genome-wide parametric linkage and association scans were analyzed using adjusted logistic regression models through the application of the Pseudomarker software. To investigate the additional contribution of rare variants in familial ET, we also performed an aggregate variant non-parametric linkage (NPL) analysis using the collapsed haplotype method implemented in CHP-NPL software. FINDINGS: Parametric linkage analysis of common variants identified several loci with significant evidence of linkage (HLOD ≥3.6). Among the gene regions within the strongest ET linkage peaks were BTC (4q13.3, HLOD=4.53), N6AMT1 (21q21.3, HLOD=4.31), PCDH9 (13q21.32, HLOD=4.21), EYA1 (8q13.3, HLOD=4.04), RBFOX1 (16p13.3, HLOD=4.02), MAPT (17q21.31, HLOD=3.99) and SCARB2 (4q21.1, HLOD=3.65). CHP-NPL analysis identified fifteen additional genes with evidence of significant linkage (LOD ≥3.8). These genes include TUBB2A, VPS33B, STEAP1B, SPINK5, ZRANB1, TBC1D3C, PDPR, NPY4R, ETS2, ZNF736, SPATA21, ARL17A, PZP, BLK and CCDC94. In one ET family contributing to the linkage peak on chromosome 16p13.3, we identified a likely pathogenic heterozygous canonical splice acceptor variant in exon 2 of RBFOX1 (ENST00000547372; c.4-2A>G), that co-segregated with the ET phenotype in the family. INTERPRETATION: Linkage and association analyses of WGS identified several novel ET candidate genes, which are implicated in four major pathways that include 1) the epidermal growth factor receptor-phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha-AKT serine/threonine kinase 1 (EGFR-PI3K-AKT) and Mitogen-activated protein Kinase 1 (ERK) pathways, 2) Reactive oxygen species (ROS) and DNA repair, 3) gamma-aminobutyric acid-ergic (GABAergic) system and 4) RNA binding and regulation of RNA processes. Our study provides evidence for a possible overlap in the genetic architecture of ET, neurological disease, cancer and aging. The genes and pathways identified can be prioritized in future genetic and functional studies. FUNDING: National Institutes of Health, NINDS, NS073872 (USA) and NIA AG058131(USA).

A splice site variant in TCTN3 underlies an atypical form of orofaciodigital syndrome IV.

Orofaciodigital syndrome (OFD) is clinically heterogeneous and is characterized by abnormalities in the oral cavity, facial features, digits, and central nervous system. At least 18 subtypes of the condition have been described in the literature. OFD is caused by variants in several genes with overlapping phenotypes. We studied a consanguineous Pakistani family with two affected siblings with an atypical form of OFD type 4 (OFD4). In addition to the typical features of OFD4 that include limb defects and growth retardation, the siblings displayed rare features of scaphocephaly and seizures. Exome sequencing analysis revealed a novel homozygous splice site variant c.257-1G>A in TCTN3 that segregated with disease. This homozygous splice site variant in TCTN3 is most likely the underlying cause of the atypical form of OFD4 observed in this family. Our results contribute to the phenotypic spectrum of TCTN3 associated ciliopathies and will facilitate better clinical diagnosis.

A monoallelic variant in EYA1 is associated with Branchio-Otic syndrome in a Malian family.

BACKGROUND: Branchio-otic syndrome (BO) is one of the most common types of syndromic hearing impairment (HI) with an incidence of 1/40,000 globally. It is an autosomal dominant disorder typically characterized by the coexistence of branchial cysts or fistulae, malformations of the external, middle, and inner ears with preauricular pits or tags and a variable degree of HI. Most cases of BO have been reported in populations of European ancestry. To date, only few cases have been reported in people from African descent. METHODS: After a careful clinical examination, a pure tone audiometry was performed. DNA was extracted from peripheral blood and whole exome, and Sanger sequencing were performed for genetic analysis. RESULTS: Eight individuals from a large non-consanguineous Malian family, with autosomal dominant inheritance were enrolled. The ages at diagnosis ranged from 8 to 54 years. A high phenotypic variability was noted among the affected individuals. Four patients presented with a post-lingual and mixed type of HI, one individual had conductive HI while three had normal hearing but presented other BO features namely branchial fistulae and preauricular sinus. Serum creatinine level and renal ultrasonography were normal in three affected individuals who performed them. Genetic testing identified a monoallelic pathogenic variant in EYA1 (c.1286A > G; p.Asp429Gly) segregating with BO syndrome in the family. CONCLUSION: This is the first genetically confirmed case of BO syndrome caused by EYA1 variant in the sub-Saharan African population, expanding the genetic spectrum of the condition.