Characterizing Homozygous Variants in Bardet-Biedl Syndrome-Associated Genes Within Iranian Families: Unveiling a Founder Variant in BBS2, c.471G>A.

Bardet-Biedl syndrome (BBS) is a rare inherited ciliopathy disorder characterized by a broad spectrum of clinical symptoms such as retinal dystrophy, obesity, polydactyly, genitourinary and kidney anomalies, learning disability, and hypogonadism. The understanding of the variants involved in BBS-causing genes remains incomplete, highlighting the need for further research to develop a molecular diagnostic strategy for this syndrome. Singleton whole-exome sequencing (WES) was performed on sixteen patients. Our study revealed (1) nine patients carried eight homozygous pathogenic variants with four of them being novel (2) Specifically, a synonymous splicing variant (c.471G > A) in BBS2 gene in six patients with Baloch ethnicity. The identification of runs of homozygosity (ROH) calling was performed using the BCFtools/RoH software on WES data of patients harboring c.471G > A variant. The presence of shared homozygous regions containing the identified variant was confirmed in these patients. In-silico analysis predicted the effect of the c.471G > A variants on BBS2 mRNA splicing. This variant results in disrupted wild-type donor site and intron retention in the mature mRNA. (3) And a deletion of exons 14 to 17 in the BBS1 gene was identified in one patient by Copy-Number Variation (CNV) analysis using the ExomeDepth pipeline. Our results identified the founder variant c.471G > A in the BBS2 gene in the Baloch ethnicity of the Iranian population. This finding can guide the diagnostic approach of this syndrome in future studies.

ZBTB11 dysfunction: spectrum of brain abnormalities, biochemical signature and cellular consequences.

Bi-allelic pathogenic variants in ZBTB11 have been associated with intellectual developmental disorder, autosomal recessive 69 (MRT69; OMIM 618383). We report five patients from three families with novel, bi-allelic variants in ZBTB11. We have expanded the clinical phenotype of MRT69, documenting varied severity of atrophy affecting different brain regions and described combined malonic and methylmalonic aciduria as a biochemical manifestation. As ZBTB11 encodes for a transcriptional regulator, we performeded chromatin immunoprecipitation-sequencing targeting ZBTB11 in fibroblasts from patients and controls. Chromatin immunoprecipitation-sequencing revealed binding of wild-type ZBTB11 to promoters in 238 genes, among which genes encoding proteins involved in mitochondrial functions and RNA processing are over-represented. Mutated ZBTB11 showed reduced binding to 61 of the targeted genes, indicating that the variants act as loss of function. Most of these genes are related to mitochondrial functions. Transcriptome analysis of the patient fibroblasts revealed dysregulation of mitochondrial functions. In addition, we uncovered that reduced binding of the mutated ZBTB11 to ACSF3 leads to decreased ACSF3 transcript level, explaining combined malonic and methylmalonic aciduria. Collectively, these results expand the clinical spectrum of ZBTB11-related neurological disease and give insight into the pathophysiology in which the dysfunctional ZBTB11 affect mitochondrial functions and RNA processing contributing to the neurological and biochemical phenotypes.

Distinct genetic variation and heterogeneity of the Iranian population.

Iran, despite its size, geographic location and past cultural influence, has largely been a blind spot for human population genetic studies. With only sparse genetic information on the Iranian population available, we pursued its genome-wide and geographic characterization based on 1021 samples from eleven ethnic groups. We show that Iranians, while close to neighboring populations, present distinct genetic variation consistent with long-standing genetic continuity, harbor high heterogeneity and different levels of consanguinity, fall apart into a cluster of similar groups and several admixed ones and have experienced numerous language adoption events in the past. Our findings render Iran an important source for human genetic variation in Western and Central Asia, will guide adequate study sampling and assist the interpretation of putative disease-implicated genetic variation. Given Iran's internal genetic heterogeneity, future studies will have to consider ethnic affiliations and possible admixture.

Comprehensive genotype-phenotype correlation in AP-4 deficiency syndrome; Adding data from a large cohort of Iranian patients.

Mutations in adaptor protein complex-4 (AP-4) genes have first been identified in 2009, causing a phenotype termed as AP-4 deficiency syndrome. Since then several patients with overlapping phenotypes, comprised of intellectual disability (ID) and spastic tetraplegia have been reported. To delineate the genotype-phenotype correlation of the AP-4 deficiency syndrome, we add the data from 30 affected individuals from 12 out of 640 Iranian families with ID in whom we detected disease-causing variants in AP-4 complex subunits, using next-generation sequencing. Furthermore, by comparing genotype-phenotype findings of those affected individuals with previously reported patients, we further refine the genotype-phenotype correlation in this syndrome. The most frequent reported clinical findings in the 101 cases consist of ID and/or global developmental delay (97%), speech disorders (92.1%), inability to walk (90.1%), spasticity (77.2%), and microcephaly (75.2%). Spastic tetraplegia has been reported in 72.3% of the investigated patients. The major brain imaging findings are abnormal corpus callosum morphology (63.4%) followed by ventriculomegaly (44.5%). Our result might suggest the AP-4 deficiency syndrome as a major differential diagnostic for unknown hereditary neurodegenerative disorders.

Biallelic missense variants in ZBTB11 can cause intellectual disability in humans.

Exploring genes and pathways underlying intellectual disability (ID) provides insight into brain development and function, clarifying the complex puzzle of how cognition develops. As part of ongoing systematic studies to identify candidate ID genes, linkage analysis and next-generation sequencing revealed Zinc Finger and BTB Domain Containing 11 (ZBTB11) as a novel candidate ID gene. ZBTB11 encodes a little-studied transcription regulator, and the two identified missense variants in this study are predicted to disrupt canonical Zn2+-binding residues of its C2H2 zinc finger domain, leading to possible altered DNA binding. Using HEK293T cells transfected with wild-type and mutant GFP-ZBTB11 constructs, we found the ZBTB11 mutants being excluded from the nucleolus, where the wild-type recombinant protein is predominantly localized. Pathway analysis applied to ChIP-seq data deposited in the ENCODE database supports the localization of ZBTB11 in nucleoli, highlighting associated pathways such as ribosomal RNA synthesis, ribosomal assembly, RNA modification and stress sensing, and provides a direct link between subcellular ZBTB11 location and its function. Furthermore, given the report of prominent brain and spinal cord degeneration in a zebrafish Zbtb11 mutant, we investigated ZBTB11-ortholog knockdown in Drosophila melanogaster brain by targeting RNAi using the UAS/Gal4 system. The observed approximate reduction to a third of the mushroom body size-possibly through neuronal reduction or degeneration-may affect neuronal circuits in the brain that are required for adaptive behavior, specifying the role of this gene in the nervous system. In conclusion, we report two ID families segregating ZBTB11 biallelic mutations disrupting Zn2+-binding motifs and provide functional evidence linking ZBTB11 dysfunction to this phenotype.

Genetics of intellectual disability in consanguineous families.

Autosomal recessive (AR) gene defects are the leading genetic cause of intellectual disability (ID) in countries with frequent parental consanguinity, which account for about 1/7th of the world population. Yet, compared to autosomal dominant de novo mutations, which are the predominant cause of ID in Western countries, the identification of AR-ID genes has lagged behind. Here, we report on whole exome and whole genome sequencing in 404 consanguineous predominantly Iranian families with two or more affected offspring. In 219 of these, we found likely causative variants, involving 77 known and 77 novel AR-ID (candidate) genes, 21 X-linked genes, as well as 9 genes previously implicated in diseases other than ID. This study, the largest of its kind published to date, illustrates that high-throughput DNA sequencing in consanguineous families is a superior strategy for elucidating the thousands of hitherto unknown gene defects underlying AR-ID, and it sheds light on their prevalence.

Iranome: A catalog of genomic variations in the Iranian population.

Considering the application of human genome variation databases in precision medicine, population-specific genome projects are continuously being developed. However, the Middle Eastern population is underrepresented in current databases. Accordingly, we established Iranome database (www.iranome.com) by performing whole exome sequencing on 800 individuals from eight major Iranian ethnic groups representing the second largest population of Middle East. We identified 1,575,702 variants of which 308,311 were novel (19.6%). Also, by presenting higher frequency for 37,384 novel or known rare variants, Iranome database can improve the power of molecular diagnosis. Moreover, attainable clinical information makes this database a good resource for classifying pathogenicity of rare variants. Principal components analysis indicated that, apart from Iranian-Baluchs, Iranian-Turkmen, and Iranian-Persian Gulf Islanders, who form their own clusters, rest of the population were genetically linked, forming a super-population. Furthermore, only 0.6% of novel variants showed counterparts in "Greater Middle East Variome Project", emphasizing the value of Iranome at national level by releasing a comprehensive catalog of Iranian genomic variations and also filling another gap in the catalog of human genome variations at international level. We introduce Iranome as a resource which may also be applicable in other countries located in neighboring regions historically called Greater Iran (Persia).

Identification of disease-causing variants in the EXOSC gene family underlying autosomal recessive intellectual disability in Iranian families.

Neurodevelopmental delay and intellectual disability (ID) can arise from numerous genetic defects. To date, variants in the EXOSC gene family have been associated with such disorders. Using next-generation sequencing (NGS), known and novel variants in this gene family causing autosomal recessive ID (ARID) have been identified in five Iranian families. By collecting clinical information on these families and comparing their phenotypes with previously reported patients, we further describe the clinical variability of ARID resulting from alterations in the EXOSC gene family, and emphasize the role of RNA processing dysregulation in ARID.

Effect of inbreeding on intellectual disability revisited by trio sequencing.

In outbred Western populations, most individuals with intellectual disability (ID) are sporadic cases, dominant de novo mutations (DNM) are frequent, and autosomal recessive ID (ARID) is very rare. Because of the high rate of parental consanguinity, which raises the risk for ARID and other recessive disorders, the prevalence of ID is significantly higher in near- and middle-east countries. Indeed, homozygosity mapping and sequencing in consanguineous families have already identified a plethora of ARID genes, but because of the design of these studies, DNMs could not be systematically assessed, and the proportion of cases that are potentially preventable by avoiding consanguineous marriages or through carrier testing is hitherto unknown. This prompted us to perform whole-exome sequencing in 100 sporadic ID patients from Iran and their healthy consanguineous parents. In 61 patients, we identified apparently causative changes in known ID genes. Of these, 44 were homozygous recessive and 17 dominant DNMs. Assuming that the DNM rate is stable, these results suggest that parental consanguinity raises the ID risk about 3.6-fold, and about 4.1 to 4.25-fold for children of first-cousin unions. These results do not rhyme with recent opinions that consanguinity-related health risks are generally small and have been "overstated" in the past.

GPR126: A novel candidate gene implicated in autosomal recessive intellectual disability.

Intellectual disability (ID), a genetically and clinically heterogeneous disorder, affects 1%-3% of the general population and is a major health problem, especially in developing countries and in populations with a high frequency of consanguineous marriage. Using whole exome sequencing, a homozygous missense variation (c.3264G>C, p.W1088C) in a plausible disease causing gene, GPR126, was identified in two patients presenting with profound ID, severe speech impairment, microcephaly, seizures during infancy, and spasticity accompanied by cerebellar hypoplasia. The role of GPR126 in radial sorting and myelination in Schwann cells suggests a mechanism of pathogenesis for ID. Involvement of GPR126 in lethal congenital contracture syndrome 9 has been identified previously, but this is the first report of a plausible candidate gene, GPR126, in ID.

CNKSR1 gene defect can cause syndromic autosomal recessive intellectual disability.

The advent of high-throughput sequencing technologies has led to an exponential increase in the identification of novel disease-causing genes in highly heterogeneous diseases. A novel frameshift mutation in CNKSR1 gene was detected by Next-Generation Sequencing (NGS) in an Iranian family with syndromic autosomal recessive intellectual disability (ARID). CNKSR1 encodes a connector enhancer of kinase suppressor of Ras 1, which acts as a scaffold component for receptor tyrosine kinase in mitogen-activated protein kinase (MAPK) cascades. CNKSR1 interacts with proteins which have already been shown to be associated with intellectual disability (ID) in the MAPK signaling pathway and promotes cell migration through RhoA-mediated c-Jun N-terminal kinase (JNK) activation. Lack of CNKSR1 transcripts and protein was observed in lymphoblastoid cells derived from affected patients using qRT-PCR and western blot analysis, respectively. Furthermore, RNAi-mediated knockdown of cnk, the CNKSR1 orthologue in Drosophila melanogaster brain, led to defects in eye and mushroom body (MB) structures. In conclusion, our findings support the possible role of CNKSR1 in brain development which can lead to cognitive impairment.

Deep sequencing reveals 50 novel genes for recessive cognitive disorders.

Common diseases are often complex because they are genetically heterogeneous, with many different genetic defects giving rise to clinically indistinguishable phenotypes. This has been amply documented for early-onset cognitive impairment, or intellectual disability, one of the most complex disorders known and a very important health care problem worldwide. More than 90 different gene defects have been identified for X-chromosome-linked intellectual disability alone, but research into the more frequent autosomal forms of intellectual disability is still in its infancy. To expedite the molecular elucidation of autosomal-recessive intellectual disability, we have now performed homozygosity mapping, exon enrichment and next-generation sequencing in 136 consanguineous families with autosomal-recessive intellectual disability from Iran and elsewhere. This study, the largest published so far, has revealed additional mutations in 23 genes previously implicated in intellectual disability or related neurological disorders, as well as single, probably disease-causing variants in 50 novel candidate genes. Proteins encoded by several of these genes interact directly with products of known intellectual disability genes, and many are involved in fundamental cellular processes such as transcription and translation, cell-cycle control, energy metabolism and fatty-acid synthesis, which seem to be pivotal for normal brain development and function.

Mutation profile of BBS genes in Iranian patients with Bardet-Biedl syndrome: genetic characterization and report of nine novel mutations in five BBS genes.

Bardet-Biedl syndrome (BBS) is a rare ciliopathy disorder that is clinically and genetically heterogeneous with 18 known genes. This study was performed to characterize responsible genes and mutation spectrum in a cohort of 14 Iranian families with BBS. Sanger sequencing of the most commonly mutated genes (BBS1, BBS2 and BBS10) accounting for ∼50% of BBS patients determined mutations only in BBS2, including three novel mutations. Next, three of the remaining patients were subjected to whole exome sequencing with 96% at 20 × depth of coverage that revealed novel BBS4 mutation. Observation of no mutation in the other patients represents the possible presence of novel genes. Screening of the remaining patients for six other genes (BBS3, BBS4, BBS6, BBS7, BBS9 and BBS12) revealed five novel mutations. This result represents another indication for the genetic heterogeneity of BBS and extends the mutational spectrum of the disease by introducing nine novel mutations in five BBS genes. In conclusion, although BBS1 and BBS10 are among the most commonly mutated genes in other populations like Caucasian, these two seem not to have an important role in Iranian patients. This suggests that a different strategy in molecular genetics diagnostic approaches in Middle Eastern countries such as Iran should be considered.

Drp1 depletion protects against ferroptotic cell death by preserving mitochondrial integrity and redox homeostasis.

Mitochondria are highly dynamic organelles which undergo constant fusion and fission as part of the mitochondrial quality control. In genetic diseases and age-related neurodegenerative disorders, altered mitochondrial fission-fusion dynamics have been linked to impaired mitochondrial quality control, disrupted organelle integrity and function, thereby promoting neural dysfunction and death. The key enzyme regulating mitochondrial fission is the GTPase Dynamin-related Protein 1 (Drp1), which is also considered as a key player in mitochondrial pathways of regulated cell death. In particular, increasing evidence suggests a role for impaired mitochondrial dynamics and integrity in ferroptosis, which is an iron-dependent oxidative cell death pathway with relevance in neurodegeneration. In this study, we demonstrate that CRISPR/Cas9-mediated genetic depletion of Drp1 exerted protective effects against oxidative cell death by ferroptosis through preserved mitochondrial integrity and maintained redox homeostasis. Knockout of Drp1 resulted in mitochondrial elongation, attenuated ferroptosis-mediated impairment of mitochondrial membrane potential, and stabilized iron trafficking and intracellular iron storage. In addition, Drp1 deficiency exerted metabolic effects, with reduced basal and maximal mitochondrial respiration and a metabolic shift towards glycolysis. These metabolic effects further alleviated the mitochondrial contribution to detrimental ROS production thereby significantly enhancing neural cell resilience against ferroptosis. Taken together, this study highlights the key role of Drp1 in mitochondrial pathways of ferroptosis and expose the regulator of mitochondrial dynamics as a potential therapeutic target in neurological diseases involving oxidative dysregulation.

Genetic Analysis of 27 Y-STR Haplotypes in 11 Iranian Ethnic Groups.

BACKGROUND: The study of Y-chromosomal variations provides valuable insights into male susceptibility in certain diseases like cardiovascular disease (CVD). In this study, we analyzed paternal lineage in different Iranian ethnic groups, not only to identify developing medical etiology, but also to pave the way for gender-specific targeted strategies and personalized medicine in medical genetic research studies. METHODS: The diversity of eleven Iranian ethnic groups was studied using 27 Y-chromosomal short tandem repeat (Y-STR) haplotypes from Y-filer® Plus kit. Analysis of molecular variance (AMOVA) based on pair-wise RST along with multidimensional scaling (MDS) calculation and Network phylogenic analysis was employed to quantify the differences between 503 unrelated individuals from each ethnicity. RESULTS: Results from AMOVA calculation confirmed that Gilaks and Azeris showed the largest genetic distance (RST=0.35434); however, Sistanis and Lurs had the smallest considerable genetic distance (RST=0.00483) compared to other ethnicities. Although Azeris had a considerable distance from other ethnicities, they were still close to Turkmens. MDS analysis of ethnic groups gave the indication of lack of similarity between different ethnicities. Besides, network phylogenic analysis demonstrated insignificant clustering between samples. CONCLUSION: The AMOVA analysis results explain that the close distance of Azeris and Turkmens may be the effect of male-dominant expansions across Central Asia that contributed to historical and demographics of populations in the region. Insignificant differences in network analysis could be the consequence of high mutation events that happened in the Y-STR regions over the years. Considering the ethnic group affiliations in medical research, our results provided an understanding and characterization of Iranian male population for future medical and population genetics studies.

Clinical application of next generation sequencing for Mendelian disease diagnosis in the Iranian population.

Next-generation sequencing (NGS) has been proven to be one of the most powerful diagnostic tools for rare Mendelian disorders. Several studies on the clinical application of NGS in unselected cohorts of Middle Eastern patients have reported a high diagnostic yield of up to 48%, correlated with a high level of consanguinity in these populations. We evaluated the diagnostic utility of NGS-based testing across different clinical indications in 1436 patients from Iran, representing the first study of its kind in this highly consanguineous population. A total of 1075 exome sequencing and 361 targeted gene panel sequencing were performed over 8 years at a single clinical genetics laboratory, with the majority of cases tested as proband-only (91.6%). The overall diagnostic rate was 46.7%, ranging from 24% in patients with an abnormality of prenatal development to over 67% in patients with an abnormality of the skin. We identified 660 pathogenic or likely pathogenic variants, including 241 novel variants, associated with over 342 known genetic conditions. The highly consanguineous nature of this cohort led to the diagnosis of autosomal recessive disorders in the majority of patients (79.1%) and allowed us to determine the shared carrier status of couples for suspected recessive phenotypes in their deceased child(ren) when direct testing was not possible. We also highlight the observations of recessive inheritance of genes previously associated only with dominant disorders and provide an expanded genotype-phenotype spectrum for multiple less-characterized genes. We present the largest mutational spectrum of known Mendelian disease, including possible founder variants, throughout the Iranian population, which can serve as a unique resource for clinical genomic studies locally and beyond.

SNUPN deficiency causes a recessive muscular dystrophy due to RNA mis-splicing and ECM dysregulation.

SNURPORTIN-1, encoded by SNUPN, plays a central role in the nuclear import of spliceosomal small nuclear ribonucleoproteins. However, its physiological function remains unexplored. In this study, we investigate 18 children from 15 unrelated families who present with atypical muscular dystrophy and neurological defects. Nine hypomorphic SNUPN biallelic variants, predominantly clustered in the last coding exon, are ascertained to segregate with the disease. We demonstrate that mutant SPN1 failed to oligomerize leading to cytoplasmic aggregation in patients' primary fibroblasts and CRISPR/Cas9-mediated mutant cell lines. Additionally, mutant nuclei exhibit defective spliceosomal maturation and breakdown of Cajal bodies. Transcriptome analyses reveal splicing and mRNA expression dysregulation, particularly in sarcolemmal components, causing disruption of cytoskeletal organization in mutant cells and patient muscle tissues. Our findings establish SNUPN deficiency as the genetic etiology of a previously unrecognized subtype of muscular dystrophy and provide robust evidence of the role of SPN1 for muscle homeostasis.

Bi-allelic loss of function variant in the NRCAM gene is associated with motor-predominant axonal polyneuropathy; the second report.

BACKGROUND: The role of biallelic variants in the NRCAM gene underlying a neurodevelopmental disorder has been defined recently. The phenotype is mainly recognized by varying severity of global developmental delay/intellectual disability, hypotonia, spasticity, and peripheral neuropathy. METHODS: Here, we describe a patient with an initial diagnosis of motor-predominant axonal polyneuropathy or a form of distal SMA. Whole-exome sequencing (WES), in parallel with WES-based CNV detection and assessment of homozygosity runs, was performed to identify this patient's possible genetic cause. RESULTS: Whole exome sequencing revealed a homozygous variant, c.73C > T (p.Gln25*), in the NRCAM gene, while the patient manifests a mild range of phenotypes compared to NRCAM-related disorder. He presented only motor-predominant axonal polyneuropathy with no other signs of central nervous system involvement. CONCLUSIONS: This study is the second report of an association between biallelic NRCAM gene variants and a Mendelian disorder. The obtained clinical data, together with the molecular findings in this patient, expands the clinical and molecular spectrum of NRCAM-related disorder and highlights its phenotypic complexity. Although patients with loss of function variants in this gene have previously presented severe clinical features, we show that type of the pathogenic variant does not necessarily determine the severity of this phenotype.

Emerging Epidemiological Data on Rare Intellectual Disability Syndromes from Analyzing the Data of a Large Iranian Cohort.

BACKGROUND: Intellectual disability (ID) is a genetically heterogeneous condition, and so far, 1679 human genes have been identified for this phenotype. Countries with a high rate of parental consanguinity, such as Iran, provide an excellent opportunity to identify the remaining novel ID genes, especially those with an autosomal recessive (AR) mode of inheritance. This study aimed to investigate the most prevalent ID genes identified via next-generation sequencing (NGS) in a large ID cohort at the Genetics Research Center (GRC) of the University of Social Welfare and Rehabilitation Sciences. METHODS: First, we surveyed the epidemiological data of 619 of 1295 families in our ID cohort, who referred to the Genetics Research Center from all over the country between 2004 and 2021 for genetic investigation via the NGS pipeline. We then compared our data with those of several prominent studies conducted in consanguineous countries. Data analysis, including cohort data extraction, categorization, and comparison, was performed using the R program version 4.1.2. RESULTS: We categorized the most common ID genes that were mutated in more than two families into 17 categories. The most common syndromic ID in our cohort was AP4 deficiency syndrome, and the most common non-syndromic autosomal recessive intellectual disability (ARID) gene was ASPM. We identified two unrelated families for the 36 ID genes. We found 14 genes in common between our cohort and the Arab and Pakistani groups, of which three genes (AP4M1, AP4S1, and ADGRG1) were repeated more than once. CONCLUSION: To date, there has been no comprehensive targeted NGS platform for the detection of ID genes in our country. Due to the large sample size of our study, our data may provide the initial step toward designing an indigenously targeted NGS platform for the diagnosis of ID, especially common ARID in our population.