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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Limbic System Associated Membrane Protein Mutation in an Iranian Family Diagnosed with Ménière's Disease.

BACKGROUND: Ménière's disease (MD) is a common inner ear disorder which is characterized by recurrent attacks of vertigo, fluctuating sensorineural hearing loss (SNHL), tinnitus, and a sense of fullness in the affected ear. MD is a complex disorder; although six genes have been linked to familial autosomal dominant form of the disease, in many cases, the exact genetic etiology remains elusive. METHODS: To elucidate the genetic causes of MD in an Iranian family, we performed exome sequencing on all members of the family: consanguineous parents and four children (two affected and two unaffected). Variant filtering was completed using a customized workflow keeping variants based on segregation with MD in autosomal recessive (AR) inheritance pattern, minor allele frequency (MAF), and in-silico prediction of pathogenicity. RESULTS: Analysis revealed that in this family, 970 variants co-segregated with MD in AR pattern, out of which eight variants (one intergenic, four intronic, and three exonic) were extremely rare. The exonic variants included a synonymous substitution in USP3 gene, an in-frame deletion in ZBED2 gene, and a rare, highly conserved deleterious missense alteration in LSAMP gene. CONCLUSION: The phenotype observed in the proband described here, i.e. vertigo, poor sense of smell, tinnitus, and borderline hearing ability, may originate from aberrant changes in the cerebellum and limbic system due to a deleterious mutation in the LSAMP gene; hence, LSAMP mutation is a possible candidate for the etiology of MD in this family.

When transcripts matter: delineating between non-syndromic hearing loss DFNB32 and hearing impairment infertile male syndrome (HIIMS).

Mutations in the CDC14A (Cell Division-Cycle 14A) gene, which encodes a conserved dual-specificity protein tyrosine phosphatase, have been identified as a cause of autosomal recessive non-syndromic hearing loss (DFNB32) and hearing impairment infertility male syndrome (HIIMS). We used next-generation sequencing to screen six deaf probands from six families segregating sensorineural moderate-to-profound hearing loss. Data analysis and variant prioritization were completed using a custom bioinformatics pipeline. We identified three homozygous loss of function variants (p.Arg345Ter, p.Arg376Ter, and p.Ala451Thrfs*43) in the CDC14A gene, segregating with deafness in each family. Of the six families, four segregated the p.Arg376Ter mutation, one family segregated the p.Arg345Ter mutation and one family segregated a novel frameshift (p.Ala451Thrfs*43) mutation. In-depth phenotyping of affected individuals ruled out secondary syndromic findings. This study implicates the p.Arg376Ter mutation might be as a founder mutation in the Iranian population. It also provides the first semen analysis for deaf males carrying mutations in exon 11 of CDC14A and reveals a genotype-phenotype correlation that delineates between DFNB32 and HIIMS. The clinical results from affected males suggest the NM_033313.2 transcript alone is sufficient for proper male fertility, but not for proper auditory function. We conclude that DFNB32 is a distinct phenotypic entity in males.

Bi-allelic Pathogenic Variants in TUBGCP2 Cause Microcephaly and Lissencephaly Spectrum Disorders.

Lissencephaly comprises a spectrum of malformations of cortical development. This spectrum includes agyria, pachygyria, and subcortical band heterotopia; each represents anatomical malformations of brain cortical development caused by neuronal migration defects. The molecular etiologies of neuronal migration anomalies are highly enriched for genes encoding microtubules and microtubule-associated proteins, and this enrichment highlights the critical role for these genes in cortical growth and gyrification. Using exome sequencing and family based rare variant analyses, we identified a homozygous variant (c.997C>T [p.Arg333Cys]) in TUBGCP2, encoding gamma-tubulin complex protein 2 (GCP2), in two individuals from a consanguineous family; both individuals presented with microcephaly and developmental delay. GCP2 forms the multiprotein γ-tubulin ring complex (γ-TuRC) together with γ-tubulin and other GCPs to regulate the assembly of microtubules. By querying clinical exome sequencing cases and through GeneMatcher-facilitated collaborations, we found three additional families with bi-allelic variation and similarly affected phenotypes including a homozygous variant (c.1843G>C [p.Ala615Pro]) in two families and compound heterozygous variants consisting of one missense variant (c.889C>T [p.Arg297Cys]) and one splice variant (c.2025-2A>G) in another family. Brain imaging from all five affected individuals revealed varying degrees of cortical malformations including pachygyria and subcortical band heterotopia, presumably caused by disruption of neuronal migration. Our data demonstrate that pathogenic variants in TUBGCP2 cause an autosomal recessive neurodevelopmental trait consisting of a neuronal migration disorder, and our data implicate GCP2 as a core component of γ-TuRC in neuronal migrating cells.

Distinct genetic variation and heterogeneity of the Iranian population.

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

G130V de novo mutation in an Iranian pedigree with nonsyndromic hearing loss without palmoplantar keratoderma.

Mutations in the GJB2 gene encoding connexin 26 (Cx26) cause autosomal recessive and rarely dominant nonsyndromic sensorineural hearing loss as well as asyndromic hearing impairment with skin problems. A dominant GJB2 mutation, c.389G > T (p.G130V), has been reported previously in association with hearing impairment and palmoplantar keratoderm. Here we report the first de novo G130V mutation of GJB2 gene in a sporadic case of hearing loss in a consanguineous Iranian family which is not associated with skin disorder.

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

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

Novel Mutations in KCNQ4, LHFPL5 and COCH Genes in Iranian Families with Hearing Impairment.

BACKGROUND: Hearing loss (HL) is the most common sensory deficit in humans, and genetic factors contribute to about half of the cases. With 112 causative genes identified so far and a disproportionate share of the genes within different ethnic groups, HL has proven to be quite heterogeneous. METHODS: Twenty Iranian families having at least 2 children with hereditary HL were initially verified to be GJB2-negative and were then subjected to whole exome sequencing (WES). Sanger sequencing was used to confirm segregation of the variant identified in each family. RESULTS: In 3 families, WES revealed 3 novel variants in KCNQ4, LHFPL5 and COCH genes. The KCNQ4 gene (DFNA2A) encodes a potassium channel (KV7.4) and the heterozygous variant identified (c.1647C>G, p.F549L) resulted in the substitution of Phe549 residing in the KV7.4 cytoplasmic region. The homozygous variant (c.34A>T, p.K12X) was identified in the LHFPL5 gene (DFNB67) which encodes a transmembrane protein, and another variant in a homozygous state (c.116T>A, p.L39X) was identified in the COCH gene which encodes a secretory protein. Pathogenic variants in the COCH gene are associated with late onset autosomal dominant hearing loss (DFNA9) but the affected individuals displayed early onset HL with a recessive mode of inheritance. CONCLUSION: The 16% contribution of GJB2 to HL in the Iranian population necessitates the discovery of the remaining causal factors. This study is the first to report KCNQ4 and COCH related HL in the Iranian population and the second study, globally, to report HL due to biallelic inactivation of the COCH gene.

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

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

Genetics of intellectual disability in consanguineous families.

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