Core myopathy in two siblings with a biallelic variant in the CACNA1S gene-A case series study

Homozygous variants of Calcium Voltage-Gated Channel Subunit Alpha1 S (CACNA1S) gene mutation were previously identified as causes of periodic paralysis and congenital early-onset myopathy, while it could be manifested as a late-onset congenital core myopathy.

Abstract

Calcium Voltage-Gated Channel Subunit Alpha1 S (CACNA1S) gene mutation has been linked to various neuromuscular conditions in recent years. Congenital myopathy with core-like features is one of the cardinal associations reported previously, causing severe respiratory insufficiency and death in neonates. Informed consent was received from the patients. Subsequently, peripheral blood leukocytes were utilized to extract genomic DNA. Moreover, exome enrichment was implemented through the Twist Human Core Exome Kit (Twist Bioscience) and exome sequenced using Illumina NovaSeq 6000 platform (Illumina, San Diego, CA, USA). Sanger sequencing using BIG Dye Terminators confirmed the presence of the final variant. Finally, the candidate variants were classified based on the American College of Medical Genetics and Genomics (ACMG) guidelines. In this report, we describe two siblings, who presented with childhood and late-onset progressive muscle weakness, and had a homozygous variant in exon 2 of the CACNA1S gene defined as c.188C > A (p.Ala63Asp) (NM_000069.3). The SIFT, Polyphen2, CADD PHRED, and Mutation Taster analysis tools classified the variant as pathogenic/damaging. The muscle biopsy of the younger brother revealed intermyofibrillar network pattern disruption as cytoplasmic core-like lesions. The muscle magnetic resonance imaging (MRI) reported grade IIa and IIb fatty changes. Finally, the electromyography (EMG) findings suggested a myopathic change pattern. This report illustrates the clinical variability in CACNA1S-related myopathy by reviewing prior reports and adding newly found aspects, additionally expanding the gene defects associated with core myopathy.

HFE and Non-HFE Hereditary Hemochromatosis Based on Screening of 854 Individuals: 12 Years of an Iranian Experience

Introduction: The genetics of hereditary hemochromatosis (HH) is understudied in Iran. Here, we report the result of genetic screening of 854 individuals, referred as "suspected cases of HH," to a diagnostic laboratory in Iran over a 12-year period. Materials and Methods: From 2011 to 2012, 121 cases were screened for HH using Sanger sequencing of HFE exons. After 2012, this method was replaced by a commercial reverse hybridization assay (RHA) targeting 18 variants in the HFE, TFR2, and FPN1(SLC40A1) genes and 733 cases were screened using this method. Results: From the total studied population, HH was confirmed by genetic diagnosis in only seven cases (0.82%): two homozygotes for HFE:C282Y and five homozygotes for TFR2:AVAQ 594-597 deletion. In 254 cases (29.7%), H63D, C282Y, S65C, and four other HFE variants not targeted by RHA were identified. Although the resulting genotypes in the latter cases did not confirm HH, some of them were known modifying factors of iron overload or could cause HH in combination with a possibly undetected variant. No variant was detected in 593 cases (69.4%). Conclusion: This study showed that the spectrum of genetic variants of HH in the Iranian population includes HFE and TFR2 variants. However, HH was not confirmed in the majority (99.2%) of suspected cases. This could be explained by limitations of our genetic diagnostics and possible inaccuracies in clinical suspicion of HH. A cooperative clinical and genetic investigation is proposed as a solution to this issue.

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.