Contribution of genetic variants in the development of familial premature coronary artery disease in a cohort of cardiac patients

Coronary artery disease (CAD), the most prevalent cardiovascular disease, is the leading cause of death worldwide. Heritable factors play a significant role in the pathogenesis of CAD. It has been proposed that approximately one-third of patients with CAD have a positive family history, and individuals with such history are at ~1.5-fold increased risk of CAD in their lifespans. Accordingly, the long-recognized familial clustering of CAD is a strong risk factor for this disease. Our study aimed to identify candidate genetic variants contributing to CAD by studying a cohort of 60 large Iranian families with at least two members in different generations afflicted with premature CAD (PCAD), defined as established disease at ≤45 years in men and ≤55 years in women. Exome sequencing was performed for a subset of the affected individuals, followed by prioritization and Sanger sequencing of candidate variants in all available family members. Subsequently, apparently healthy carriers of potential risk variants underwent coronary computed tomography angiography (CCTA), followed by co-segregation analysis of the combined data. Putative causal variants were identified in seven genes, ABCG8, CD36, CYP27A1, PIK3C2G, RASSF9, RYR2, and ZFYVE21, co-segregating with familial PCAD in seven unrelated families. Among these, PIK3C2G, RASSF9, and ZFYVE21 are novel candidate CAD susceptibility genes. Our findings indicate that rare variants in genes identified in this study are involved in CAD development.

The First Iranian Case of Unstable Hemoglobin Santa Ana

In this report, we describe a 6-year-old girl with a medical history of pallor, mild icterus, anemia, blood transfusion and abnormal hemoglobin variant analysis on capillary electrophoresis. She was referred for further analysis. DNA sequencing of the proband revealed a de novo mutation in Codon 88 (CTG > CCG) of the β-globin gene (HBB: c.266T > C) in a heterozygous state compatible with hemoglobin Santa Ana, an unstable hemoglobin. This is the first case of Hb Santa Ana from Iran associated with moderate to severe anemia who underwent splenectomy with clinical improvement.

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.

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.

Identification of a Novel Variant in CC2D1A Gene Linked to Autosomal Recessive Intellectual Disability 3 in an Iranian Family and Investigating the Structure and Pleiotropic Effects of this Gene

Objectives

Intellectual disability (ID) represents a significant health challenge due to its diverse and intricate nature. A multitude of genes play a role in brain development and function, with defects in these genes potentially leading to ID. Considering that many of these genes have yet to be identified, and those identified have only been found in a small number of patients, no complete description of the phenotype created by these genes is available. CC2D1A is one of the genes whose loss-of-function mutation leads to a rare form of non-syndromic ID-3(OMIM*610055), and four pathogenic variants have been reported in this gene so far.

Materials & Methods

n the current study, two affected females were included with an initial diagnosis of ID who were from an Iranian family with consanguineous marriage. Whole-exome sequencing was used to identify the probable genetic defects. The Genotypic and phenotypic characteristics of the patients were compared with a mutation in the CC2D1A gene, and then the structure of the gene and its reported variants were investigated.

Results

The patients carried a novel homozygous splicing variant (NM_017721, c.1641+1G>A) in intron 14, which is pathogenic according to the ACMG guideline. Loss-of-function mutations in CC2D1A have severe phenotypic consequences such as ID, autism spectrum disorder (ASD), and seizures. However, missense mutations lead to ASD with or without ID, and in some patients, they cause ciliopathy.

Conclusion

This study reports the fifth novel, probably pathogenic variant in the CC2D1A gene. Comparing the clinical and molecular genetic features of the patients with loss-of-function mutation helped to describe the phenotype caused by this gene more precisely. Investigating the CC2D1A gene's mutations and structure revealed that it performs multiple functions. The DM14 domain appears more pivotal in triggering severe clinical symptoms, including ID, than the C2 domain.

The liver-derived exosomes stimulate insulin gene expression in pancreatic beta cells under condition of insulin resistance

Introduction

An insufficient functional beta cell mass is a core pathological hallmark of type 2 diabetes (T2D). Despite the availability of several effective pharmaceuticals for diabetes management, there is an urgent need for novel medications to protect pancreatic beta cells under diabetic conditions. Integrative organ cross-communication controls the energy balance and glucose homeostasis. The liver and pancreatic islets have dynamic cross-communications where the liver can trigger a compensatory beta cell mass expansion and enhanced hormonal secretion in insulin-resistant conditions. However, the indispensable element(s) that foster beta cell proliferation and insulin secretion have yet to be completely identified. Exosomes are important extracellular vehicles (EVs) released by most cell types that transfer biological signal(s), including metabolic messengers such as miRNA and peptides, between cells and organs.

Methods

We investigated whether beta cells can take up liver-derived exosomes and examined their impact on beta cell functional genes and insulin expression. Exosomes isolated from human liver HepG2 cells were characterized using various methods, including Transmission Electron Microscopy (TEM), dynamic light scattering (DLS), and Western blot analysis of exosomal markers. Exosome labeling and cell uptake were assessed using CM-Dil dye. The effect of liver cell-derived exosomes on Min6 beta cells was determined through gene expression analyses of beta cell markers and insulin using qPCR, as well as Akt signaling using Western blotting.

Results

Treatment of Min6 beta cells with exosomes isolated from human liver HepG2 cells treated with insulin receptor antagonist S961 significantly increased the expression of beta cell markers Pdx1, NeuroD1, and Ins1 compared to the exosomes isolated from untreated cells. In line with this, the activity of AKT kinase, an integral component of the insulin receptor pathway, is elevated in pancreatic beta cells, as represented by an increase in AKT's downstream substrate, FoxO1 phosphorylation.

Discussions

This study suggests that liver-derived exosomes may carry a specific molecular cargo that can affect insulin expression in pancreatic beta cells, ultimately affecting glucose homeostasis.

The Spectrum of HBB Mutations among 2315 Beta Thalassemia Patients of a Reference Clinic in Tehran-Iran

Beta Thalassemia is the most prevalent and well-studied single gene disorder in Iran. Here, we investigated the spectrum of HBB gene mutations, identified among 2315 patients, referred to a reference thalassemia clinic in Tehran, on the basis of suspicion to thalassemia major or intermedia. The patients were homozygous or compound heterozygous for HBB mutations, and were referred from various Iranian provinces, during 15 years (2001- 2016). The HBB mutations were classified based on their frequency, and the result was compared to a meta-analysis of 14,293 beta thalassemia cases in the Iranian population, within the same time period. The mutation spectrum in this study contained 43 HBB mutations, compared to the 90, presented by the meta-analysis. Similar to the meta-analysis, IVSII-1 (G > A) and IVSI-5 (G > C) were the most common mutations in this study. These two comprised 62.40% of the total HBB mutant alleles in the studied population, comparable to 51.92% of that in the meta-analysis. IVSII-1 (G > A) and IVSI-5 (G > C), followed by 17 other mutations that had frequencies ranging from 0.15% to 5.44%, were among the 20 common HBB mutations in Iran and neighboring countries, according to the meta-analysis. This study provided further evidence to support the spectrum of the most common HBB mutations in the Iranian population.

The prevalence and phenotypic range associated with biallelic PKDCC variants

PKDCC encodes a component of Hedgehog signalling required for normal chondrogenesis and skeletal development. Although biallelic PKDCC variants have been implicated in rhizomelic shortening of limbs with variable dysmorphic features, this association was based on just two patients. In this study, data from the 100 000 Genomes Project was used in conjunction with exome sequencing and panel-testing results accessed via international collaboration to assemble a cohort of eight individuals from seven independent families with biallelic PKDCC variants. The allelic series included six frameshifts, a previously described splice-donor site variant and a likely pathogenic missense variant observed in two families that was supported by in silico structural modelling. Database queries suggested that the prevalence of this condition is between 1 of 127 and 1 of 721 in clinical cohorts with skeletal dysplasia of unknown aetiology. Clinical assessments, combined with data from previously published cases, indicate a predominantly upper limb involvement. Micrognathia, hypertelorism and hearing loss appear to be commonly co-occurring features. In conclusion, this study strengthens the link between biallelic inactivation of PKDCC and rhizomelic limb-shortening and will enable clinical testing laboratories to better interpret variants in this gene.

Expanding the Molecular Spectrum of HK1-Related Charcot-Marie-Tooth Disease, Type 4G; the First Report in Iran

Charcot-Marie-Tooth disease type 4G (CMT4G) was first reported in Balkan Gypsies as a myelinopathy starting with progressive distal lower limb weakness, followed by upper limb involvement and prominent distal sensory impairment later in the patient's life. So far, CMT4G has been only reported in European Roma communities with two founder homozygous variants; g.9712G>C and g.11027G>A, located in the 5'-UTR of the HK1 gene. Here, we present the first Iranian CMT4G patient manifesting progressive distal lower limb weakness from 11 years of age and diagnosed with chronic demyelinating sensorimotor polyneuropathy. Whole-exome sequencing for this patient revealed a homozygous c.19C>T (p. Arg7*) variant in the HK1 gene. This report expands the mutational spectrum of the HK1-related CMT disorder and provides supporting evidence for the observation of CMT4G outside the Roma population. Interestingly, the same Arg7* variant is recently observed in another unrelated Pakistani CMT patient, proposing a possible prevalence of this variant in the Middle Eastern populations.

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.

Identification of a homozygous frameshift mutation in the FGF3 gene in a consanguineous Iranian family: First report of labyrinthine aplasia, microtia, and microdontia syndrome in Iran and literature review

BACKGROUND

To date, over 400 syndromes with hearing impairment have been identified which altogether constitute almost 30% of hereditary hearing loss (HL) cases around the globe. Manifested as complete or partial labyrinthine aplasia (severe malformations of the inner ear structure), type I microtia (smaller outer ear with shortened auricles), and microdontia (small and widely spaced teeth), labyrinthine aplasia, microtia, and microdontia (LAMM) syndrome (OMIM 610706) is an extremely rare autosomal recessive condition caused by bi-allelic mutations in the FGF3 gene.

METHODS

Using the whole-exome sequencing (WES) data of the proband, we analyzed a consanguineous Iranian family with three affected members presenting with congenital bilateral HL, type I microtia, and microdontia.

RESULTS

We discovered the homozygous deletion c.45delC in the first exon of the FGF3 gene, overlapping a 38.72 Mb homozygosity region in chromosome 11. Further investigations using Sanger sequencing revealed that this variant co-segregated with the phenotype observed in the family.

CONCLUSION

Here, we report the first identified case of LAMM syndrome in Iran, and by identifying a frameshift variant in the first exon of the FGF3 gene, our result will help better clarify the phenotype-genotype relation of LAMM syndrome.

An Extended Iranian Family with Autosomal Dominant Non-syndromic Hearing Loss Associated with A Nonsense Mutation in the DIAPH1 Gene

Genetic analysis of non-syndromic hearing loss (NSHL) has been challenged due to marked clinical and genetic heterogeneity. Today, advanced next-generation sequencing (NGS) technologies, such as exome sequencing (ES), have drastically increased the efficacy of gene identification in heterogeneous Mendelian disorders. Here, we present the utility of ES and re-evaluate the phenotypic data for identifying candidate causal variants for previously unexplained progressive moderate to severe NSHL in an extended Iranian family. Using this method, we identified a known heterozygous nonsense variant in exon 26 of the DIAPH1 gene (MIM: 602121), which led to "Deafness, autosomal dominant 1, with or without thrombocytopenia; DFNA1" (MIM: 124900) in this large family in the absence of GJB2 disease-causing variants and also OtoSCOPE-negative results. To the best of our knowledge, this nonsense variant (NM_001079812.3):c.3610C>T (p.Arg1204Ter) is the first report of the DIAPH1 gene variant for autosomal dominant non-syndromic hearing loss (ADNSHL) in Iran.

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.

Implementation of an In-House Platform for Rapid Screening of SARS-CoV-2 Genome Variations

BACKGROUND

Global real-time monitoring of SARS-CoV-2 variants is crucial to controlling the COVID-19 outbreak. The purpose of this study was to set up a Sanger-based platform for massive SARS-CoV-2 variant tracking in laboratories in low-resource settings.

METHODS

We used nested RT-PCR assay, Sanger sequencing and lineage assignment for 930-bp of the SARS-CoV-2 spike gene, which harbors specific variants of concern (VOCs) mutations. We set up our platform by comparing its results with whole genome sequencing (WGS) data on 137 SARS-CoV-2 positive samples. Then, we applied it on 1028 samples from March-September 2021.

RESULTS

In total, 125 out of 137 samples showed 91.24% concordance in mutation detection. In lineage assignment, 123 out of 137 samples demonstrated 89.78% concordance, 65 of which were assigned as VOCs and showed 100% concordance. Of 1028 samples screened by our in-house method, 78 distinct mutations were detected. The most common mutations were: S:D614G (21.91%), S:P681R (12.19%), S:L452R (12.15%), S:T478K (12.15%), S:N501Y (8.91%), S:A570D (8.89%), S:P681H (8.89%), S:T716I (8.74%), S:L699I (3.50%) and S:S477N (0.28%). Of 1028 samples, 980 were attributed as VOCs, which include the Delta (B.1.617.2) and Alpha (B.1.1.7) variants.

CONCLUSION

Our proposed in-house Sanger-based assay for SARS-CoV-2 lineage assignment is an accessible strategy in countries with poor infrastructure facilities. It can be applied in the rapid tracking of SARS-CoV-2 VOCs in the SARS-CoV-2 pandemic.

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.

Aberrant phase separation and nucleolar dysfunction in rare genetic diseases

Thousands of genetic variants in protein-coding genes have been linked to disease. However, the functional impact of most variants is unknown as they occur within intrinsically disordered protein regions that have poorly defined functions^1-3. Intrinsically disordered regions can mediate phase separation and the formation of biomolecular condensates, such as the nucleolus^4,5. This suggests that mutations in disordered proteins may alter condensate properties and function^6-8. Here we show that a subset of disease-associated variants in disordered regions alter phase separation, cause mispartitioning into the nucleolus and disrupt nucleolar function. We discover de novo frameshift variants in HMGB1 that cause brachyphalangy, polydactyly and tibial aplasia syndrome, a rare complex malformation syndrome. The frameshifts replace the intrinsically disordered acidic tail of HMGB1 with an arginine-rich basic tail. The mutant tail alters HMGB1 phase separation, enhances its partitioning into the nucleolus and causes nucleolar dysfunction. We built a catalogue of more than 200,000 variants in disordered carboxy-terminal tails and identified more than 600 frameshifts that create arginine-rich basic tails in transcription factors and other proteins. For 12 out of the 13 disease-associated variants tested, the mutation enhanced partitioning into the nucleolus, and several variants altered rRNA biogenesis. These data identify the cause of a rare complex syndrome and suggest that a large number of genetic variants may dysregulate nucleoli and other biomolecular condensates in humans.

Characterizing Genotypes and Phenotypes Associated with Dysfunction of Channel-Encoding Genes in a Cohort of Patients with Intellectual Disability

BACKGROUND

Ion channel dysfunction in the brain can lead to impairment of neuronal membranes and generate several neurological diseases, especially neurodevelopmental disorders.

METHODS

In this study, we set out to delineate the genotype and phenotype spectrums of 14 Iranian patients from 7 families with intellectual disability (ID) and/or developmental delay (DD) in whom genetic mutations were identified by next-generation sequencing (NGS) in 7 channel-encoding genes: KCNJ10, KCNQ3, KCNK6, CACNA1C, CACNA1G, SCN8A, and GRIN2B. Moreover, the data of 340 previously fully reported ID and/or DD cases with a mutation in any of these seven genes were combined with our patients to clarify the genotype and phenotype spectrum in this group.

RESULTS

In total, the most common phenotypes in 354 cases with ID/DD in whom mutation in any of these 7 channel-encoding genes was identified were as follows: ID (77.4%), seizure (69.8%), DD (59.8%), behavioral abnormality (29.9%), hypotonia (21.7%), speech disorder (21.5%), gait disturbance (20.9%), and ataxia (20.3%). Electroencephalography abnormality (33.9%) was the major brain imaging abnormality.

CONCLUSION

The results of this study broaden the molecular spectrum of channel pathogenic variants associated with different clinical presentations in individuals with ID and/or DD.

Genetic Diagnosis of Pyruvate Kinase Deficiency in Undiagnosed Iranian Patients with Severe Hemolytic Anemia, using Whole Exome Sequencing

BACKGROUND

After ruling out the most common causes of severe hemolytic anemia by routine diagnostic tests, certain patients remain without a diagnosis. The aim of this study was to elucidate the genetic cause of the disease in these patients using next generation sequencing (NGS).

METHODS

Four unrelated Iranian families including six blood transfusion dependent cases and their parents were referred to us from a specialist center in Tehran. There was no previous history of anemia in the families and the parents had no abnormal hematological presentations. All probands presented severe congenital hemolytic anemia, neonatal jaundice and splenomegaly. Common causes of hemolytic anemia were ruled out prior to this investigation in these patients and they had no diagnosis. Whole exome sequencing (WES) was performed in the probands and the results were confirmed by Sanger sequencing and subsequent family studies.

RESULTS

We identified five variants in the PKLR gene, including a novel unpublished frameshift in these families. These variants were predicted as pathogenic according to the ACMG guidelines by Intervar and/or Varsome prediction tools. Subsequent family studies by Sanger sequencing supported the diagnosis of pyruvate kinase deficiency (PKD) in six affected individuals and the carrier status of disease in their parents.

CONCLUSION

These findings show that PKD is among the rare blood disorders that could remain undiagnosed or even ruled out in Iranian population without performing NGS. This could be due to pitfalls in clinical, hematological or biochemical approaches in diagnosing PKD. Furthermore, genotyping PKD patients in Iran could reveal novel mutations in the PKLR gene.

Targeted Next Generation Sequencing Revealed Novel Variants in the PKD1 and PKD2 Genes of Iranian Patients with Autosomal Dominant Polycystic Kidney Disease

BACKGROUND

Autosomal dominant polycystic kidney disease (ADPKD), one of the common inherited disorders in humans, is characterized by the development and enlargement of renal cysts, often leading to end-stage renal disease (ESRD). In this study, Iranian ADPKD families were subjected to high-throughput DNA sequencing to find potential causative variants facilitating the way toward risk assessment and targeted therapy.

METHODS

Our protocol was based on the targeted next generation sequencing (NGS) panel previously developed in our center comprising 12 genes involved in PKD. This panel has been applied to investigate the genetic causes of 32 patients with a clinical suspicion of ADPKD.

RESULTS

We identified a total of 31 variants for 32 individuals, two of which were each detected in two individuals. Twenty-seven out of 31 detected variants were interpreted as pathogenic/likely pathogenic and the remaining 4 of uncertain significance with a molecular diagnostic success rate of 87.5%. Among these variants, 25 PKD1/2 pathogenic/likely pathogenic variants were detected in 32 index patients (78.1%), and variants of uncertain significance in four individuals (12.5% in PKD1/2). The majority of variants was identified in PKD1 (74.2%). Autosomal recessive PKD was identified in one patient, indicating the similarities between recessive and dominant PKD. In concordance with earlier studies, this biallelic PKD1 variant, p.Arg3277Cys, leads to rapidly progressive and severe disease with very early-onset ADPKD.

CONCLUSION

Our findings suggest that targeted gene panel sequencing is expected to be the method of choice to improve diagnostic and prognostic accuracy in PKD patients with heterogeneity in genetic background.

Disease Waves of SARS-CoV-2 in Iran Closely Mirror Global Pandemic Trends

BACKGROUND

Complete SARS-CoV-2 genome sequencing in the early phase of the outbreak in Iran showed two independent viral entries. Subsequently, as part of a genome surveillance project, we aimed to characterize the genetic diversity of SARS-CoV-2 in Iran over one year after emerging.

METHODS

We provided 319 SARS-CoV-2 whole-genome sequences used to monitor circulating lineages in March 2020-May 2021 time interval.

RESULTS

The temporal dynamics of major SARS-CoV-2 clades/lineages circulating in Iran is comparable to the global perspective and represent the 19A clade (B.4) dominating the first disease wave, followed by 20A (B.1.36), 20B (B.1.1.413), 20I (B.1.1.7), leading the second, third and fourth waves, respectively. We observed a mixture of circulating B.1.36, B.1.1.413, B.1.1.7 lineages in winter 2021, paralleled in a fading manner for B.1.36/B.1.1.413 and a growing rise for B.1.1.7, prompting the fourth outbreak. Entry of the Delta variant, leading to the fifth disease wave in summer 2021, was detected in April 2021. This study highlights three lineages as hallmarks of the SARS-CoV-2 outbreak in Iran; B4, dominating early periods of the epidemic, B.1.1.413 (B.1.1 with the combination of [D138Y-S477N-D614G] spike mutations) as a characterizing lineage in Iran, and the co-occurrence of [I100T-L699I] spike mutations in half of B.1.1.7 sequences mediating the fourth peak. It also designates the renowned combination of G and GR clades' mutations as the top recurrent mutations.

CONCLUSION

In brief, we provided a real-time and comprehensive picture of the SARS-CoV-2 genetic diversity in Iran and shed light on the SARS-CoV-2 transmission and circulation on the regional scale.

Intellectual disability associated with craniofacial dysmorphism due to POLR3B mutation and defect in spliceosomal machinery

Background

Intellectual disability (ID) is a clinically important disease and a most prevalent neurodevelopmental disorder. The etiology and pathogenesis of ID are poorly recognized. Exome sequencing revealed a homozygous missense mutation in the POLR3B gene in a consanguineous family with three Intellectual disability with craniofacial anomalies patients. POLR3B gene encoding the second largest subunit of RNA polymerase III.

Methods

We performed RNA sequencing on blood samples to obtain insights into the biological pathways influenced by POLR3B mutation. We applied the results of our RNA-Seq analysis to several gene ontology programs such as ToppGene, Enrichr, KEGG.

Results

A significant decrease in expression of several spliceosomal RNAs, ribosomal proteins, and transcription factors was detected in the affected, compared to unaffected, family members.

Conclusions

We hypothesize that POLR3B mutation dysregulates the expression of some important transcription factors, ribosomal and spliceosomal genes, and impairments in protein synthesis and splicing mediated in part by transcription factors such as FOXC2 and GATA1 contribute to impaired neuronal function and concurrence of intellectual disability and craniofacial anomalies in our patients. Our study highlights the emerging role of the spliceosome and ribosomal proteins in intellectual disability.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01237-5.

Identification of microRNAs associated with human fragile X syndrome using next-generation sequencing

Fragile X syndrome (FXS) is caused by a mutation in the FMR1 gene which can lead to a loss or shortage of the FMR1 protein. This protein interacts with specific miRNAs and can cause a range of neurological disorders. Therefore, miRNAs could act as a novel class of biomarkers for common CNS diseases. This study aimed to test this theory by exploring the expression profiles of various miRNAs in Iranian using deep sequencing-based technologies and validating the miRNAs affecting the expression of the FMR1 gene. Blood samples were taken from 15 patients with FXS (9 males, 6 females) and 12 controls. 25 miRNAs were differentially expressed in individuals with FXS compared to controls. Levels of 9 miRNAs were found to be significantly changed (3 upregulated and 6 downregulated). In Patients, the levels of hsa-miR-532-5p, hsa-miR-652-3p and hsa-miR-4797-3p were significantly upregulated while levels of hsa-miR-191-5p, hsa-miR-181-5p, hsa-miR-26a-5p, hsa-miR-30e-5p, hsa-miR-186-5p, and hsa-miR-4797-5p exhibited significant downregulation; and these dysregulations were confirmed by RT-qPCR. This study presents among the first evidence of altered miRNA expression in blood samples from patients with FXS, which could be used for diagnostic, prognostic, and treatment purposes. Larger studies are required to confirm these preliminary results.

Phenotype and genotype spectrum of variants in guanine nucleotide exchange factor genes in a broad cohort of Iranian patients

Background

Guanine nucleotide exchange factors (GEFs) play pivotal roles in neuronal cell functions by exchanging GDP to GTP nucleotide and activation of GTPases. We aimed to determine the genotype and phenotype spectrum of GEF mutations by collecting data from a large Iranian cohort with intellectual disability (ID) and/or developmental delay (DD).

Methods

We collected data from nine families with 20 patients extracted from Iranian cohort of 640 families with ID and/or DD. Next‐generation sequencing (NGS) was used to identify the causing variants in recruited families. We also compared our clinical and molecular findings with previously reported patients carrying mutations in these GEF genes in the literature published until mid‐2021.

Results

We identified disease‐causing variants in eight GEF genes including ALS2, IQSEC2, MADD, RAB3GAP1, RAB3GAP2, TRIO, ITSN1, and DENND2A. The major clinical manifestations in 203 previously reported cases along with our 20 patients with disease causing variants in eight GEF genes were as follow; speech disorder (85.2%), ID (81.6%), DD (81.1%), inability to walk (71.3%), facial dysmorphisms features (52.4%), abnormalities in skull morphology (55.6%), hypotonia and muscle weakness (47%), and brain MRI abnormalities (43.4%).

Conclusion

Our study provides new insights into the genotype and phenotype spectrum of mutations in GEF genes.

DNA Banking to Assess Genetic Influences on Schizophrenia

Background: Schizophrenia is among the most prevalent psychiatric disorders globally, with a lifetime prevalence rate of 0.3% to 0.7%, characterized by the heterogeneous presence of positive, negative, and cognitive symptoms that affect all aspects of mental activity. We aimed to describe the genetics of schizophrenia to widening our understanding of the inheritance of this illness.

Methods: This quasi-experimental study was conducted in Razi psychiatric hospital in Tehran province, Iran. Recruitment of the study samples was conducted in Tehran, Iran, among patients with schizophrenia and their families. For this purpose, individuals with schizophrenia in 40 families with at least 1 to 2 affected members were identified and selected based on a clinical interview conducted by a psychiatrist and according to the Diagnostic and Statistical Manual of Mental Disorders, 5th Edition. The clinical and paraclinical data, drug and substance usage, and medical treatments were collected through a standardized clinical questionnaire. Besides, the Global Assessment Scale and the Positive and Negative Syndrome Scale were completed for all study participants.

Results: A total of 22 families had a negative family history, and 1 affected member and the rest of the studied families had a positive family history and at least 2 affected members. In addition, genealogical data (family tree) and lymphoblastic cell categories were developed to examine genes, and subsequent research results will be reported in the future.

Conclusion: As the research continues, the approach to sampling must be modified to ensure that the deoxyribonucleic acid bank is as extensively representative as possible of all schizophrenia cases.

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.

The Effect of Xmn -1 Polymorphism and Coinheritance of Alpha Mutations on Age at First Blood Transfusion in Iranian Patients with Homozygote IVSI-5 Mutation

Background: Thalassemia syndromes are the most prevalent hereditary hemoglobinopathies in the world. Iran is located on the thalassemia belt.  In this study, the effect of Xmn -1 polymorphism and coinheritance of alpha mutations on age at first transfusion and also transfusion interval in Iranian thalassemic patients with homozygous IVSI-5 mutation were assessed.

Materials and Methods: In this retrospective cross-sectional study 154 transfusion dependent thalassemia (TDT) patients (140 patients with β-thalassemia major and 14 cases with β-thalassemia intermedia) who were homozygote of IVSI-5 mutation have been participated. Blood samples were collected from participants using EDTA containers for genomic DNA analysis. DNA extraction and amplification-refractory mutation to determine the Xmn -1 polymorphism were performed. Multiplex PCR was performed to identify alpha globin deletions. 

Results: The mean age of participants was 29±7, 58 of them were male and 96 were female. A significant relation between presence of Xmn -1 polymorphism and age at receiving first transfusion was detected. Coinheritance of alpha thalassemia mutation does not have significant effect on age at first transfusion or transfusion interval.

Conclusion: Presence of Xmn -1 polymorphism can delay the onset of transfusion in patients with homozygote IVSI-5 mutation.

A novel variant of C12orf4 linked to autosomal recessive intellectual disability type 66 with phenotype expansion

BACKGROUND

Intellectual disability (ID) is a hallmark of many rare disorders that are highly heterogeneous and complex. A large number of specific genes are involved in development of this heterogeneity, and each of these genes is only found in a small number of patients. This weakens the definition of the predominant genotype and the phenotypic characteristics associated with that gene. Autosomal recessive ID type 66 (OMIM # 618221) is one of these very rare diseases created by defects in C12orf4 gene.

METHODS

The researchers in the current study included two patients from an Iranian family with initial diagnosis of non-syndromic ID to identify the possible genetic cause(s), and whole-exome sequencing (WES) was performed for the proband. The obtained variant was confirmed by Sanger sequencing and co-segregated in the family.

RESULT

The patients carried a novel pathogenic splicing variant called c.1441-1G>A in exon 12 of C12orf4 gene (NM_001304811). They predominantly manifested ID, behavioral problems, speech impairment, and dysmorphic facial features, some of which had not been reported in previous studies.

CONCLUSION

A novel pathogenic splicing variant was identified named c.1441-1G>A in C12orf4 gene. So far, only seven families have been reported with defects in this gene. Previous studies have not highlighted the exact clinical manifestations of these patients, thus, this study could contribute to better delineation of the genotype-phenotype correlation and interpretation of very rare variants of the gene.

The recurrent missense mutation p.(Arg367Trp) in YARS1 causes a distinct neurodevelopmental phenotype

Pathogenic variants in aminoacyl-tRNA synthetases (ARS1) cause a diverse spectrum of autosomal recessive disorders. Tyrosyl tRNA synthetase (TyrRS) is encoded by YARS1 (cytosolic, OMIM*603,623) and is responsible of coupling tyrosine to its specific tRNA. Next to the enzymatic domain, TyrRS has two additional functional domains (N-Terminal TyrRSMini and C-terminal EMAP-II-like domain) which confer cytokine-like functions. Mutations in YARS1 have been associated with autosomal-dominant Charcot-Marie-Tooth (CMT) neuropathy type C and a heterogenous group of autosomal recessive, multisystem diseases. We identified 12 individuals from 6 families with the recurrent homozygous missense variant c.1099C > T;p.(Arg367Trp) (NM_003680.3) in YARS1. This variant causes a multisystem disorder with developmental delay, microcephaly, failure to thrive, short stature, muscular hypotonia, ataxia, brain anomalies, microcytic anemia, hepatomegaly, and hypothyroidism. In silico analyses show that the p.(Arg367Trp) does not affect the catalytic domain responsible of enzymatic coupling, but destabilizes the cytokine-like C-terminal domain. The phenotype associated with p.(Arg367Trp) is distinct from the other biallelic pathogenic variants that reside in different functional domains of TyrRS which all show some common, but also divergent clinical signs [(e.g., p.(Phe269Ser)-retinal anomalies, p.(Pro213Leu)/p.(Gly525Arg)-mild ID, p.(Pro167Thr)-high fatality)]. The diverse clinical spectrum of ARS1-associated disorders is related to mutations affecting the various non-canonical domains of ARS1, and impaired protein translation is likely not the exclusive disease-causing mechanism of YARS1- and ARS1-associated neurodevelopmental disorders. KEY MESSAGES: The missense variant p.(Arg367Trp) in YARS1 causes a distinct multisystem disorder. p.(Arg367Trp) affects a non-canonical domain with cytokine-like functions. Phenotypic heterogeneity associates with the different affected YARS1 domains. Impaired protein translation is likely not the exclusive mechanism of ARS1-associated disorders.

The PTRHD1 Mutation in Intellectual Disability

BACKGROUND

Intellectual disability (ID) is a heterogonous disorder with complex etiology. The frequency of autosomal recessive inheritance defects was elevated in a consanguineous family.

METHODS

In this study, high-throughput DNA sequencing was performed in an Iranian consanguineous family with two affected individuals to find potential causative variants. Whole-exome sequencing was carried out on the proband and Sanger sequencing was implemented for validation of the likely causative variant in the family members.

RESULTS

A novel homozygous missense mutation (p.Arg122Trp) was detected in the PTRHD1 gene.

CONCLUSION

PTRHD1 has been recently introduced as a candidate ID and Parkinsonism causing gene. Our findings are in agreement with the clinical spectrum of PTRHD1 mutations; however, our affected individuals suffer from ID manifestations.

Exome sequencing utility in defining the genetic landscape of hearing loss and novel-gene discovery in Iran

Hearing loss (HL) is one of the most common sensory defects affecting more than 466 million individuals worldwide. It is clinically and genetically heterogeneous with over 120 genes causing non-syndromic HL identified to date. Here, we performed exome sequencing (ES) on a cohort of Iranian families with no disease-causing variants in known deafness-associated genes after screening with a targeted gene panel. We identified likely causal variants in 20 out of 71 families screened. Fifteen families segregated variants in known deafness-associated genes. Eight families segregated variants in novel candidate genes for HL: DBH, TOP3A, COX18, USP31, TCF19, SCP2, TENM1, and CARMIL1. In the three of these families, intrafamilial locus heterogeneity was observed with variants in both known and novel candidate genes. In aggregate, we were able to identify the underlying genetic cause of HL in nearly 30% of our study cohort using ES. This study corroborates the observation that high-throughput DNA sequencing in populations with high rates of consanguineous marriages represents a more appropriate strategy to elucidate the genetic etiology of heterogeneous conditions such as HL.

A Novel β0-Thalassemia Mutation, HBB: c.356_357delTT [Codon 118 (-TT)] in an Iraqi Kurd

We report a novel frameshift β-thalassemia (β-thal) mutation due to a two-nucleotide deletion at codon 118 of the β-globin gene (HBB: c.356_357delTT) in a 4-year-old Iraqi Kurd female presenting as transfusion-dependent β-thal. This frameshift mutation, unlike many others involving the third exon, behaved as a recessive β⁰ defect and not as dominant β-thal mutation.

CEP104 and CEP290; Genes with Ciliary Functions Cause Intellectual Disability in Multiple Families

BACKGROUND

Neurodevelopmental and intellectual impairments are extremely heterogeneous disorders caused by a diverse variety of genes involved in different molecular pathways and networks. Genetic alterations in cilia, highly-conserved organelles with sensorineural and signal transduction roles can compromise their proper functions and lead to so-called "ciliopathies" featuring intellectual disability (ID) or neurodevelopmental disorders as frequent clinical manifestations. Here, we report several Iranian families affected with ID and other ciliopathy-associated features carrying known and novel variants in two ciliary genes; CEP104 and CEP290.

METHODS

Whole exome and Targeted exome sequencing were carried out on affected individuals. Lymphoblastoid cell lines (LCLs) derived from the members of affected families were established for two families carrying CEP104 mutations. RNA and protein expression studies were carried out on these cells using qPCR and Western blot, respectively.

RESULTS

A novel homozygous variant; NM_025114.3:c.7341_7344dupACTT p.(Ser2449Thrfs*8) and four previously reported homozygous variants; NM_025114.3:c.322C>T p.(Arg108*), NM_025114.3:c.4393C>T p.(Arg1465*), NM_025114.3:c.5668G>T p.(Gly1890*) and NM_025114.3:c.1666dupA p.(Ile556Asnfs*20) were identified in CEP290. In two other families, two novel homozygous variants; NM_014704:c.2356_2357insTT p.(Cys786Phefs*11) and NM_014704:c.1901_1902insT p.(Leu634Phefs*33) were identified in CEP104, another ciliary gene. qPCR and Western blot analyses showed significantly lower levels of CEP104 transcripts and protein in patients compared to heterozygous or normal family members.

CONCLUSION

We emphasize on the clinical variability and pleiotropic phenotypes due to variants of these genes. In conclusion, our findings support the pivotal role of these genes resulting in cognitive and neurodevelopmental features.

SARS-CoV-2 outbreak in Iran: The dynamics of the epidemic and evidence on two independent introductions

The SARS‐CoV‐2 virus has been rapidly spreading globally since December 2019, triggering a pandemic, soon after its emergence. While Iran was among the first countries confronted with rapid spread of virus in February 2020, no real‐time SARS‐CoV‐2 whole‐genome tracking in early phase of outbreak was performed in the country. To address this issue, we provided 50 whole‐genome sequences of viral isolates ascertained from different geographical locations in Iran during March–July 2020. The corresponding analysis on origins, transmission dynamics and genetic diversity of SARS‐CoV‐2 virus, represented at least two introductions of the virus into the country, constructing two major clusters defined as B.4 and B.1*. The first entry of the virus might have occurred around very late 2019/early 2020, as suggested by the time to the most recent common ancestor, followed by a rapid community transmission that led to dominancy of B.4 lineage in early epidemic till the end of June. Gradually, reduction in dominancy of B.4 occurred possibly as a result of other entries of the virus, followed by surge of B.1* lineages, as of mid‐May. Remarkably, variation tracking of the virus indicated the increase in frequency of D614G mutation, along with B.1* lineages, which showed continuity till October 2020. The increase in frequency of D614G mutation and B.1* lineages from mid‐May onwards predicts a rapid viral transmission that may push the country into a critical health situation followed by a considerable change in composition of viral lineages circulating in the country.

Discovery of a neuromuscular syndrome caused by biallelic variants in ASCC3

Activating Signal Cointegrator 1 Complex, Subunit 3 (ASCC3) is part of the four-part ASC-1 transcriptional cointegrator complex. This complex includes ASCC1 (associated with spinal muscular atrophy with congenital bone fractures 2), TRIP4 (associated with spinal muscular atrophy with congenital bone fractures 1), and ASCC2 (not yet associated with human disease.) ASCC3 encodes a DNA helicase responsible for generating single-stranded DNA as part of the DNA damage response. Interestingly, ASCC3 expresses coding and non-coding isoforms, which act in opposition to balance the recovery of gene transcription after UV-induced DNA damage. Here we report the discovery of ASCC3 as the cause of a neuromuscular syndrome in seven unreported individuals from six unrelated families and updates on the one previously reported family. All the individuals share a neurologic phenotype that ranges from severe developmental delay to muscle fatigue. There appears to be genotype-phenotype correlation, as the most mildly affected individual is homozygous for a rare missense variant, while the more severely affected individuals are compound heterozygotes for a missense and a presumed loss-of-function (LOF) variant. There are no individuals with biallelic presumed LOF variants in our cohort or in gnomAD, as this genotype may not be compatible with life. In summary we report a syndrome in these eleven individuals from seven families with biallelic variants in ASCC3.

Identifying the causes of recurrent pregnancy loss in consanguineous couples using whole exome sequencing on the products of miscarriage with no chromosomal abnormalities

Recurrent miscarriages occur in about 5% of couples trying to conceive. In the past decade, the products of miscarriage have been studied using array comparative genomic hybridization (a-CGH). Within the last decade, an association has been proposed between miscarriages and single or multigenic changes, introducing the possibility of detecting other underlying genetic factors by whole exome sequencing (WES). We performed a-CGH on the products of miscarriage from 1625 Iranian women in consanguineous or non-consanguineous marriages. WES was carried out on DNA extracted from the products of miscarriage from 20 Iranian women in consanguineous marriages and with earlier normal genetic testing. Using a-CGH, a statistically significant difference was detected between the frequency of imbalances in related vs. unrelated couples (P < 0.001). WES positively identified relevant alterations in 11 genes in 65% of cases. In 45% of cases, we were able to classify these variants as pathogenic or likely pathogenic, according to the American College of Medical Genetics and Genomics guidelines, while in the remainder, the variants were classified as of unknown significance. To the best of our knowledge, our study is the first to employ WES on the products of miscarriage in consanguineous families with recurrent miscarriages regardless of the presence of fetal abnormalities. We propose that WES can be helpful in making a diagnosis of lethal disorders in consanguineous couples after prior genetic testing.

Clinical and Genetic Characteristics of Splicing Variant in CYP27A1 in an Iranian Family with Cerebrotendinous Xanthomatosis

Background:

CTX is a rare congenital lipid-storage disorder, leading to a progressive multisystem disease. CTX with autosomal recessive inheritance is caused by a defect in the CYP27A1 gene. Chronic diarrhea, tendon xanthomas, neurologic impairment, and bilateral cataracts are common symptoms of the disease.

Methods:

Three affected siblings with an initial diagnosis of non-syndromic intellectual disability were recruited for further molecular investigations. To identify the possible genetic cause(s), WES was performed on the proband. Sanger sequencing was applied to confirm the final variant. The clinical and molecular genetic features of the three siblings from the new CTX family and other patients with the same mutations, as previously reported, were analyzed. The CYP27A1 gene was also studied for the number of pathogenic variants and their location.

Results:

We found a homozygous splicing mutation, NM_000784: exon6: c.1184+1G>A, in CYP27A1 gene, which was confirmed by Sanger sequencing. Among the detected pathogenic variants, the splice site mutation had the highest prevalence, and the mutations were mostly found in exon 4.

Conclusion:

This study is the first to report the c.1184+1G>A mutation in Iran. Our findings highlight the other feature of the disease, which is the lack of relationship between phenotype and genotype. Due to nonspecific symptoms and delay in diagnosis, CYP27A1 genetic analysis should be the definitive method for CTX diagnosis.

YIF1B mutations cause a post-natal neurodevelopmental syndrome associated with Golgi and primary cilium alterations

Human post-natal neurodevelopmental delay is often associated with cerebral alterations that can lead, by themselves or associated with peripheral deficits, to premature death. Here, we report the clinical features of 10 patients from six independent families with mutations in the autosomal YIF1B gene encoding a ubiquitous protein involved in anterograde traffic from the endoplasmic reticulum to the cell membrane, and in Golgi apparatus morphology. The patients displayed global developmental delay, motor delay, visual deficits with brain MRI evidence of ventricle enlargement, myelination alterations and cerebellar atrophy. A similar profile was observed in the Yif1b knockout (KO) mouse model developed to identify the cellular alterations involved in the clinical defects. In the CNS, mice lacking Yif1b displayed neuronal reduction, altered myelination of the motor cortex, cerebellar atrophy, enlargement of the ventricles, and subcellular alterations of endoplasmic reticulum and Golgi apparatus compartments. Remarkably, although YIF1B was not detected in primary cilia, biallelic YIF1B mutations caused primary cilia abnormalities in skin fibroblasts from both patients and Yif1b-KO mice, and in ciliary architectural components in the Yif1b-KO brain. Consequently, our findings identify YIF1B as an essential gene in early post-natal development in human, and provide a new genetic target that should be tested in patients developing a neurodevelopmental delay during the first year of life. Thus, our work is the first description of a functional deficit linking Golgipathies and ciliopathies, diseases so far associated exclusively to mutations in genes coding for proteins expressed within the primary cilium or related ultrastructures. We therefore propose that these pathologies should be considered as belonging to a larger class of neurodevelopmental diseases depending on proteins involved in the trafficking of proteins towards specific cell membrane compartments.

Subcellular relocalization and nuclear redistribution of the RNA methyltransferases TRMT1 and TRMT1L upon neuronal activation

RNA modifications are dynamic chemical entities that expand the RNA lexicon and regulate RNA fate. The most abundant modification present in mRNAs, N6-methyladenosine (m⁶A), has been implicated in neurogenesis and memory formation. However, whether additional RNA modifications may be playing a role in neuronal functions and in response to environmental queues is largely unknown. Here we characterize the biochemical function and cellular dynamics of two human RNA methyltransferases previously associated with neurological dysfunction, TRMT1 and its homolog, TRMT1-like (TRMT1L). Using a combination of next-generation sequencing, LC-MS/MS, patient-derived cell lines and knockout mouse models, we confirm the previously reported dimethylguanosine (m^2,2G) activity of TRMT1 in tRNAs, as well as reveal that TRMT1L, whose activity was unknown, is responsible for methylating a subset of cytosolic tRNA^Ala(AGC) isodecoders at position 26. Using a cellular in vitro model that mimics neuronal activation and long term potentiation, we find that both TRMT1 and TRMT1L change their subcellular localization upon neuronal activation. Specifically, we observe a major subcellular relocalization from mitochondria and other cytoplasmic domains (TRMT1) and nucleoli (TRMT1L) to different small punctate compartments in the nucleus, which are as yet uncharacterized. This phenomenon does not occur upon heat shock, suggesting that the relocalization of TRMT1 and TRMT1L is not a general reaction to stress, but rather a specific response to neuronal activation. Our results suggest that subcellular relocalization of RNA modification enzymes may play a role in neuronal plasticity and transmission of information, presumably by addressing new targets.

NGLY1 deficiency: Novel variants and literature review

NGLY1 deficiency is a recently described autosomal recessive disorder, involved in deglycosylation of proteins, and for that reason grouped as the congenital disorders of deglycosylation together with the lysosomal storage disorders. The typical phenotype is characterized by intellectual disability, liver malfunctioning, muscular hypotonia, involuntary movements, and decreased or absent tear production. Liver biopsy demonstrates vacuolar amorphous cytoplasmic storage material. NGLY1 deficiency is caused by bi-allelic variants in NGLY1 which catalyzes protein deglycosylation. We describe five patients from two families with NGLY1 deficiency due to homozygosity for two novel NGLY1 variants, and compare their findings to those of earlier reported patients. The typical features of the disorder are present in a limited way, and there is intra-familial variability. In addition in one of the families the muscle atrophy and posture abnormalities are marked. These can be explained either as variability of the phenotype or as sign of slowly progression of features as the present affected individuals are older than earlier reported patients.

Biallelic mutations in the death domain of PIDD1 impair caspase-2 activation and are associated with intellectual disability

PIDD1 encodes p53-Induced Death Domain protein 1, which acts as a sensor surveilling centrosome numbers and p53 activity in mammalian cells. Early results also suggest a role in DNA damage response where PIDD1 may act as a cell-fate switch, through interaction with RIP1 and NEMO/IKKg, activating NF-κB signaling for survival, or as an apoptosis-inducing protein by activating caspase-2. Biallelic truncating mutations in CRADD—the protein bridging PIDD1 and caspase-2—have been reported in intellectual disability (ID), and in a form of lissencephaly. Here, we identified five families with ID from Iran, Pakistan, and India, with four different biallelic mutations in PIDD1, all disrupting the Death Domain (DD), through which PIDD1 interacts with CRADD or RIP1. Nonsense mutations Gln863* and Arg637* directly disrupt the DD, as does a missense mutation, Arg815Trp. A homozygous splice mutation in the fifth family is predicted to disrupt splicing upstream of the DD, as confirmed using an exon trap. In HEK293 cells, we show that both Gln863* and Arg815Trp mutants fail to co-localize with CRADD, leading to its aggregation and mis-localization, and fail to co-precipitate CRADD. Using genome-edited cell lines, we show that these three PIDD1 mutations all cause loss of PIDDosome function. Pidd1 null mice show decreased anxiety, but no motor abnormalities. Together this indicates that PIDD1 mutations in humans may cause ID (and possibly lissencephaly) either through gain of function or secondarily, due to altered scaffolding properties, while complete loss of PIDD1, as modeled in mice, may be well tolerated or is compensated for.

Biallelic mutations in the death domain of PIDD1 impair caspase-2 activation and are associated with intellectual disability

PIDD1 encodes p53-Induced Death Domain protein 1, which acts as a sensor surveilling centrosome numbers and p53 activity in mammalian cells. Early results also suggest a role in DNA damage response where PIDD1 may act as a cell-fate switch, through interaction with RIP1 and NEMO/IKKg, activating NF-κB signaling for survival, or as an apoptosis-inducing protein by activating caspase-2. Biallelic truncating mutations in CRADD-the protein bridging PIDD1 and caspase-2-have been reported in intellectual disability (ID), and in a form of lissencephaly. Here, we identified five families with ID from Iran, Pakistan, and India, with four different biallelic mutations in PIDD1, all disrupting the Death Domain (DD), through which PIDD1 interacts with CRADD or RIP1. Nonsense mutations Gln863* and Arg637* directly disrupt the DD, as does a missense mutation, Arg815Trp. A homozygous splice mutation in the fifth family is predicted to disrupt splicing upstream of the DD, as confirmed using an exon trap. In HEK293 cells, we show that both Gln863* and Arg815Trp mutants fail to co-localize with CRADD, leading to its aggregation and mis-localization, and fail to co-precipitate CRADD. Using genome-edited cell lines, we show that these three PIDD1 mutations all cause loss of PIDDosome function. Pidd1 null mice show decreased anxiety, but no motor abnormalities. Together this indicates that PIDD1 mutations in humans may cause ID (and possibly lissencephaly) either through gain of function or secondarily, due to altered scaffolding properties, while complete loss of PIDD1, as modeled in mice, may be well tolerated or is compensated for.

Novel Mutation in LARP7 in Two Iranian Consanguineous Families with Syndromic Intellectual Disability and Facial Dysmorphism

BACKGROUND

Recently, we have reported mutations in LARP7 gene, leading to neurodevelopmental disorders (NDDs), the most frequent cause of disability in children with a broad phenotype spectrum and diverse genetic landscape.

METHODS

Here, we present two Iranian patients from consanguineous families with syndromic intellectual disability, facial dysmorphism, and short stature.

RESULTS

Whole-exome sequencing (WES) revealed a novel homozygous stop-gain (c.C925T, p.R309X) variant and a previously known homozygous acceptor splice-site (c.1669-1_1671del) variant in LARP7 gene, indicating the diagnosis of Alazami syndrome.

CONCLUSION

These identified variants in patients with Alazami syndrome were consistent with previously reported loss of function variants in LARP7 and provide further evidence that loss of function of LARP7 is the disease mechanism.

Heterozygosity for the Novel HBA2: c.*91_*92delTA Polyadenylation Site Variant on the α2-Globin Gene Expanding the Genetic Spectrum of α-Thalassemia in Iran

There are four copy numbers of α-globin genes (16p13.3) in the human genome and the number of defective α-globin genes dictates the severity of α-thalassemia (α-thal). Mutations that occur in the 3' untranslated region (3'UTR), and especially at the polyadenylation (polyA) sites, affect the translation, stability and export of mRNA. A patient with hypochromic microcytic anemia was referred to the Kariminejad-Najmabadi Pathology & Genetics Center, Tehran, Iran by the health network. Molecular analysis of genomic DNA for the evaluation of mutations on the α- and β-globin genes was performed. Direct sequencing of the hemoglobin (Hb) subunit α2 (HBA2) gene revealed a two nucleotide deletion between +816 and +817 in the 3'UTR, located at the polyA site, which seems to be a novel pathogenic variant. This novel variant expands the genetic spectrum of α-thal in the 3'UTR of the HBA2 gene.

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.

Whole genome sequencing identifies a duplicated region encompassing Xq13.2q13.3 in a large Iranian family with intellectual disability

Background

The X chromosome has historically been one of the most thoroughly investigated chromosomes regarding intellectual disability (ID), whose etiology is attributed to many factors including copy number variations (CNVs). Duplications of the long arm of the X chromosome have been reported in patients with ID, short stature, facial anomalies, and in many cases hypoplastic genitalia and/or behavioral abnormalities.

Methods

Here, we report on a large Iranian family with X‐linked ID caused by a duplication on the X chromosome identified by whole genome sequencing in combination with linkage analysis.

Results

Seven affected males in different branches of the family presented with ID, short stature, seizures, facial anomalies, behavioral abnormalities (aggressiveness, self‐injury, anxiety, impaired social interactions, and shyness), speech impairment, and micropenis. The duplication of the region Xq13.2q13.3, which is ~1.8 Mb in size, includes seven protein‐coding OMIM genes. Three of these genes, namely SLC16A2, RLIM, and NEXMIF, if impaired, can lead to syndromes presenting with ID. Of note, this duplicated region was located within a linkage interval with a LOD score >3.

Conclusion

Our report indicates that CNVs should be considered in multi‐affected families where no candidate gene defect has been identified in sequencing data analysis.

Adult-onset very-long-chain acyl-CoA dehydrogenase deficiency (VLCADD)

BACKGROUND AND PURPOSE

Very-long-chain acyl-CoA dehydrogenase deficiency (VLCADD) is a hereditary disorder of mitochondrial long-chain fatty acid oxidation that has variable presentations, including exercise intolerance, cardiomyopathy and liver disease. The aim of this study was to describe the clinical and genetic manifestations of six patients with adult-onset VLCADD.

METHODS

In this study, the clinical, pathological and genetic findings of six adult patients (four from Iran and two from Serbia) with VLCADD and their response to treatment are described.

RESULTS

The median (range) age of patients at first visit was 31 (27-38) years, and the median (range) age of onset was 26.5 (19-33) years. Parental consanguinity was present for four patients. Four patients had a history of rhabdomyolysis, and the recorded CK level ranged between 67 and 90 000 IU/l. Three patients had a history of exertional myalgia, and one patient had a non-fluctuating weakness. Through next-generation sequencing analysis, we identified six cases with variants in the ACADVL gene and a confirmed diagnosis of VLCADD. Of the total six variants identified, five were missense, and one was a novel frameshift mutation identified in two unrelated individuals. Two variants were novel, and three were previously reported. We treated the patients with a combination of L-carnitine, Coenzyme Q10 and riboflavin. Three patients responded favorably to the treatment.

CONCLUSION

Adult-onset VLCADD is a rare entity with various presentations. Patients may respond favorably to a cocktail of L-carnitine, Coenzyme Q10, and riboflavin.

Molecular Diagnosis of Hereditary Neuropathies by Whole Exome Sequencing and Expanding the Phenotype Spectrum

BACKGROUND

Inherited peripheral neuropathies (IPNs) are a group of neuropathies affecting peripheral motor and sensory neurons. Charcot-Marie-Tooth (CMT) disease is the most common disease in this group. With recent advances in next-generation sequencing (NGS) technologies, more than 100 genes have been implicated for different types of CMT and other clinically and genetically inherited neuropathies. There are also a number of genes where neuropathy is a major feature of the disease such as spinocerebellar ataxia (SCA) and hereditary spastic paraplegia (HSP). We aimed to determine the genetic causes underlying IPNs in Iranian families.

METHODS

We performed whole exome sequencing (WES) for 58 PMP22 deletion-/duplication-negative unrelated Iranian patients with a spectrum of phenotypes and with a preliminary diagnosis of hereditary neuropathies.

RESULTS

Twenty-seven (46.6%) of the cases were genetically diagnosed with pathogenic or likely pathogenic variants. In this study, we identified genetically strong variants within genes not previously linked to any established disease phenotype in five (8.6%) patients.

CONCLUSION

Our results highlight the advantage of using WES for genetic diagnosis in highly heterogeneous diseases such as IPNs. Moreover, functional analysis is required for novel and uncertain variants.