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.

Copy Number Variations in Hereditary Spastic Paraplegia-Related Genes: Evaluation of an Iranian Hereditary Spastic Paraplegia Cohort and Literature Review.

Introduction: In human genetic disorders, copy number variations (CNVs) are considered a considerable underlying cause. CNVs are generally detected by array-based methods but can also be discovered by read-depth analysis of whole-exome sequencing (WES) data. We performed WES-based CNV identification in a cohort of 35 Iranian families with hereditary spastic paraplegia (HSP) patients. Methods: Thirty-five patients whose routine single-nucleotide variants (SNVs) and insertion/deletion analyses from exome data were unrevealing underwent a pipeline of CNV analysis using the read-depth detection method. Subsequently, a comprehensive search about the existence of CNVs in all 84 known HSP-causing genes was carried out in all reported HSP cases, so far. Results and Discussion: CNV analysis of exome data indicated that 1 patient harbored a heterozygous deletion in exon 17 of the SPAST gene. Multiplex ligation-dependent probe amplification analysis confirmed this deletion in the proband and his affected father. Literature review demonstrated that, to date, pathogenic CNVs have been identified in 30 out of 84 HSP-causing genes (∼36%). However, CNVs in only 17 of these genes were specifically associated with the HSP phenotype. Among them, CNVs were more common in L1CAM, PLP1, SPAST, SPG7, SPG11, and REEP1 genes. The identification of the CNV in 1 of our patients suggests that WES allows the detection of both SNVs and CNVs from a single method without additional costs and execution time. However, because of intrinsic issues of WES in the detection of large rearrangements, it may not yet be exploited to replace the CNV detection methods in standard clinical practice.

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.

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.

Genetic etiology of hearing loss in Iran.

Hearing loss (HL) is an etiologically heterogeneous disorder that affects around 5% of the world's population. There has been an exponential increase in the identification of genes and variants responsible for hereditary HL over recent years. Iran, a country located in the Middle East, has a high prevalence of consanguineous marriages, so heterogeneous diseases such as HL are more common. Comprehensive studies using different strategies from linkage analysis to next-generation sequencing, especially exome-sequencing, have achieved significant success in identifying possible pathogens in deaf Iranian families. About 12% of non-syndromic autosomal recessive HL genes investigated to date, were first identified in families from Iran. Variations of 56 genes have been observed in families with NSHL in Iran. Variants in GJB2, SLC26A4, MYO15A, MYO7A, CDH23, and TMC1 account for 16.5%, 16.25%, 13.5%, 9.35%, 6.9% and 4.92%, cases of NSHL, respectively. In summary, there are also different diagnostic rates between studies conducted in Iran. In the comprehensive investigations conducted by the Genetic Research Center of the University of Social Welfare and Rehabilitation Sciences over the past 20 years, the overall diagnosis rate is about 80% while there are other studies with lower diagnostic rates which could reflect differences in project designs, sampling, and accuracy and validity of the methods used. Furthermore, there are several syndromic HHLs in Iran including, Waardenburg syndrome, BOR syndrome, Brown-Vialetto-Van Laere syndrome, Wolfram syndrome, among which Pendred and Usher syndromes are well-studied. These results are of importance for further investigation and elucidation of the molecular basis of HHL in Iran.

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.

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.

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.

Chromosomal aberrations in pregnancy and fetal loss: Insight on the effect of consanguinity, review of 1625 cases.

BACKGROUND: Pregnancy loss affects 10%-15% of pregnancies and is caused by several factors, maternal and fetal. Most common cause is chromosomal aneuploidy and has traditionally been detected by karyotyping product of conception and/or fetal tissue. In recent years, array comparative genomic hybridization (a-CGH) has been used because of its higher detection and lower failure rates. METHODS: DNA was extracted from 1625 products of abortion or fetal tissue. In 1,104 cases both quantitative fluorescent-polymerase chain reaction (QF-PCR) and a-CGH, and in 521 cases only a-CGH, was performed. RESULTS: The detection rate using QF-PCR and a-CGH is 20% compared to 12.7%, overall, and 15.7%, excluding failed samples, by karyotypes in our center. QF-PCR and a-CGH failed in 1.9% of cases, while the failure rate for karyotypes was 20.1%. The difference of detection and failure rates is significant (p-value < 0.001 and p-value < 0.001 respectively). Unexpectedly we also found a significant difference in frequency of imbalances in related versus unrelated couples. (χ2  = 11.4926, p-value < 0.001). CONCLUSION: It is highly likely that the pregnancy loss in consanguineous couples is caused by other genetic and immune mechanisms. It is plausible that, through the same mechanism by which single gene disorders have a higher prevalence of manifesting disease in consanguineous couples, they can cause lethal genetic disorders leading to pregnancy loss and intra-uterine fetal death (IUFD) in these couples. Our findings suggest that this is a matter for further study as it will greatly influence the approach to counseling and managing consanguineous couples with pregnancy loss.

SLC52A2 mutations cause SCABD2 phenotype: A second report.

INTRODUCTION: Autosomal recessive cerebellar ataxias (ARCAs) are a large group of neurodegenerative disorders that manifest mainly in children and young adults. Most ARCAs are heterogeneous with respect to age at onset, severity of disease progression, and frequency of extracerebellar and systemic signs. METHODS: The phenotype of a consanguineous Iranian family was characterized using clinical testing and pedigree analysis. Whole-exome sequencing was used to identify the disease-causing gene in this family. RESULTS AND CONCLUSION: Using whole exome sequencing (WES), a novel missense mutation in SLC52A2 gene is reported in a consanguineous Iranian family with progressive severe hearing loss, optic atrophy and ataxia. This is the second report of the genotype-phenotype correlation between this syndrome named spinocerebellar ataxia with blindness and deafness type 2 (SCABD2) and SLC52A2 gene.

C-reactive protein and complement factor H polymorphism interaction in advanced exudative age-related macular degeneration.

PURPOSE: To determine the association of C-reactive protein (CRP) and complement factor H (CFH) gene with exudative age-related macular degeneration (AMD) and any possible interaction among these factors. METHODS: In this case-control study, 139 unrelated patients with exudative AMD and 123 non-AMD controls were recruited. Blood sample was taken for analysis of the CRP levels and DNA testing. DNA fragments of CFH gene variants containing 4 single nucleotide polymorphisms including rs800292, rs1061170, rs2274700, and rs3753395 were assessed. A CRP level of ≥3 mg/L was considered as elevated. The association of elevated CRP and CFH gene variants polymorphism with exudative AMD was compared between the groups. RESULTS: Mean age was 72.6 ± 6.4 for controls and 74.9 ± 7.4 for case group (P = 0.006). The difference between CRP levels in cases and controls was not statistically significant (P = 0.055). However, Y402H variant of CFH in both homozygous and heterozygous carriers C allele was significantly more frequent among exudative AMD patients than controls, 32.1 versus 6.5 % (P < 0.001). Evaluating various CRP levels in patients with CC and non-CC genotypes disclosed that in CC genotype group, higher CRP level (>3 mg/L) was associated with higher risk of developing exudative AMD (OR = 12.0, CI: 1.5-98.8) compared with the control group. CONCLUSION: This study disclosed no difference in CRP levels per se between exudative AMD patients with control group. However, higher levels of CRP in the presence of C allele of Y402H might confer more risk for the development of exudative AMD.

Heterogeneity of Hereditary Hearing Loss in Iran: a Comprehensive Review.

A significant contribution to the causes of hereditary hearing impairment comes from genetic factors. More than 120 genes and 160 loci have been identified to be involved in hearing impairment. Given that consanguine populations are more vulnerable to most inherited diseases, such as hereditary hearing loss (HHL), the genetic picture of HHL among the Iranian population, which consists of at least eight ethnic subgroups with a high rate of intermarriage, is expected to be highly heterogeneous. Using an electronic literature review through various databases such as PubMed, MEDLINE, and Scopus, we review the current picture of HHL in Iran. In this review, we present more than 39 deafness genes reported to cause non-syndromic HHL in Iran, of which the most prevalent causative genes include GJB2, SLC26A4, MYO15A, and MYO7A. In addition, we highlight some of the more common genetic causes of syndromic HHL in Iran. These results are of importance for further investigation and elucidation of the molecular basis of HHL in Iran and also for developing a national diagnostic tool tailored to the Iranian context enabling early and efficient diagnosis of hereditary hearing impairment.

Investigating the CFH Gene Polymorphisms as a Risk Factor for Age-related Macular Degeneration in an Iranian Population.

BACKGROUND: Age-related macular degeneration (AMD) is a complex disorder which results in irreversible vision loss and progressive impairment of central vision. Disease susceptibility is influenced by multiple genetic and environmental factors. Single nucleotide polymorphisms (SNP) in the complement factor H gene are the most important genetic risk factors. We conducted a case-control study to investigate the association four SNPs (dbSNP ID: rs800292, rs1061170, rs2274700 and rs3753395) of CFH gene with AMD in the Iranian population. MATERIALS AND METHODS: We recruited 100 AMD patients and 100 age- and sex-matched normal controls. Direct sequencing for three SNPs (rs800292, rs2274700 and rs3753395) and restriction fragment length polymorphism utilized for rs1061170. Allele and genotype frequencies of SNPs were calculated and tested for departure from Hardy-Weinberg equilibrium using the Chi-square test. An allelic and genotypic association was compared by logistic regression analysis using the SNPassoc. RESULTS: According to our results, the frequencies of risk allele for all SNPs (G, G, A, and C alleles of rs800292, rs2274700, rs3753395 and rs1061170, respectively) were significantly higher in AMD patients (p value < 0.001). AMD individuals who had at least one copy of the C allele of rs1061170 had an increased risk of disease compared with cases with the T allele. Other studied polymorphisms showed the same association. CONCLUSION: Our results suggest the contribution of all four predicted CFH polymorphisms in AMD susceptibility among the Iranian population. This association with CFH may lead to early detection and new strategies for prevention and treatment of AMD.

Impact of whole exome sequencing among Iranian patients with autosomal recessive retinitis pigmentosa.

BACKGROUND: Non-syndromic autosomal recessive Retinitis Pigmentosa (arRP) is a highly heterogeneous genetic visual disorder with a large number of causative genes. We aimed to determine the power of Whole Exome Sequencing (WES) in the identification of the genes responsible for non-syndromic arRP among Iranian patients. METHODS: We used WES, followed by the Sanger sequencing to identify the underlying gene mutations causing non-syndromic arRP. RESULTS: Our study revealed disease-causing mutations in known arRP genes for 10 of the 13 families studied (76.9%). These mutations included two-frameshift insertion/deletion in CRB1 and ABCA4, one splicing mutation in PDE6B, four missense mutations in RP1, CRB1, PANK2 and IFT140, as well as three stop codon mutations in RDH12, PRCD, and C2orf71. Three remaining families harbored no mutation in previously known RP genes. Of the 10 diseases causing mutations identified among the investigated Iranian patients with non-syndromic arRP, eight variants had not been reported previously. We confirmed segregation of all 10 mutations with disease phenotypes in our studied population. CONCLUSION: This study supports the genetic heterogeneity of non-syndromic arRP in Iranian patients, and provides an opportunity to show the effectiveness of WES in the identification of pathogenic mutations among patients with non-syndromic arRP born to consanguineous parents.

Characterising the spectrum of autosomal recessive hereditary hearing loss in Iran.

BACKGROUND: Countries with culturally accepted consanguinity provide a unique resource for the study of rare recessively inherited genetic diseases. Although hereditary hearing loss (HHL) is not uncommon, it is genetically heterogeneous, with over 85 genes causally implicated in non-syndromic hearing loss (NSHL). This heterogeneity makes many gene-specific types of NSHL exceedingly rare. We sought to define the spectrum of autosomal recessive HHL in Iran by investigating both common and rarely diagnosed deafness-causing genes. DESIGN: Using a custom targeted genomic enrichment (TGE) panel, we simultaneously interrogated all known genetic causes of NSHL in a cohort of 302 GJB2-negative Iranian families. RESULTS: We established a genetic diagnosis for 67% of probands and their families, with over half of all diagnoses attributable to variants in five genes: SLC26A4, MYO15A, MYO7A, CDH23 and PCDH15. As a reflection of the power of consanguinity mapping, 26 genes were identified as causative for NSHL in the Iranian population for the first time. In total, 179 deafness-causing variants were identified in 40 genes in 201 probands, including 110 novel single nucleotide or small insertion-deletion variants and three novel CNV. Several variants represent founder mutations. CONCLUSION: This study attests to the power of TGE and massively parallel sequencing as a diagnostic tool for the evaluation of hearing loss in Iran, and expands on our understanding of the genetics of HHL in this country. Families negative for variants in the genes represented on this panel represent an excellent cohort for novel gene discovery.

PDZD7 and hearing loss: More than just a modifier.

Deafness is the most frequent sensory disorder. With over 90 genes and 110 loci causally implicated in non-syndromic hearing loss, it is phenotypically and genetically heterogeneous. Here, we investigate the genetic etiology of deafness in four families of Iranian origin segregating autosomal recessive non-syndromic hearing loss (ARNSHL). We used a combination of linkage analysis, homozygosity mapping, and a targeted genomic enrichment platform to simultaneously screen 90 known deafness-causing genes for pathogenic variants. Variant segregation was confirmed by Sanger sequencing. Linkage analysis and homozygosity mapping showed segregation with the DFNB57 locus on chromosome 10 in two families. Targeted genomic enrichment with massively parallel sequencing identified causal variants in PDZD7: a homozygous missense variant (p.Gly103Arg) in one family and compound heterozygosity for missense (p.Met285Arg) and nonsense (p.Tyr500Ter) variants in the second family. Screening of two additional families identified two more variants: (p.Gly228Arg) and (p.Gln526Ter). Variant segregation with the hearing loss phenotype was confirmed in all families by Sanger sequencing. The missense variants are predicted to be deleterious, and the two nonsense mutations produce null alleles. This report is the first to show that mutations in PDZD7 cause ARNSHL, a finding that offers addition insight into the USH2 interactome. We also describe a novel likely disease-causing mutation in CIB2 and illustrate the complexity associated with gene identification in diseases that exhibit large genetic and phenotypic heterogeneity.

Association of polymorphisms at LDLR locus with coronary artery disease independently from lipid profile.

Coronary artery disease (CAD) is the leading cause of mortality in many parts of the world. Genome-wide association studies (GWAS) have identified several genetic variants associated with CAD in Low-density lipoprotein receptor (LDLR) locus. This study was evaluated the possible association of genetic markers at LDLR locus with CAD irrespective to lipid profile and as well as the association of these SNPs with severity of CAD in Iranian population. Sequencing of 2 exons in LDLR gene (Exon 2, 12) and part of intron 30 of SMARCA4 gene include rs1122608, was performed in 170 Iranian patients angiographically confirmed CAD and 104 healthy controls by direct sequencing. Sullivan's scoring system was used for determining the severity of CAD in cases. Our results showed that homozygote genotypes of rs1122608 (P<0.0001), rs4300767 (P<0.005) and rs10417578 (p<0.007) SNPs have strong protective effects on the CAD. In addition, we found that rs1122608 (GT or TT) was at higher risk of three vessel involvement compared to single vessels affecting (P=0.01).

Investigation of the asporin gene polymorphism as a risk factor for knee osteoarthritis in Iran.

Osteoarthritis (OA) is a degenerative disease of the joints characterized by degradation of the hyaline articular cartilage and remodeling of the subchondral bone with sclerosis. The asporin (ASPN) gene encodes a cartilage extracellular protein belonging to the small leucine-rich proteoglycan family. Polymorphisms in the aspartic acid (D) repeat region in the second exon of this gene, which consist of GAT repeats, are associated with OA susceptibility. The D14 allele, which contains 14 D-repeats, is associated with increased OA susceptibility in the Japanese and the Han Chinese but is not an important factor in OA etiology among Caucasians, though the D15 allele is a risk allele for the Greek population. To examine the possibility of this controversial association, we explored the effect of ASPN on Iranians with knee OA. The D-repeat polymorphism was genotyped in 100 knee OA patients and 100 controls, and the allelic association of the D-repeat was examined. Our data suggest that the D15 allele could be considered a risk allele significant only for women (P = .045, odds ratio = 1.73, 95% confidence interval [CI] = 1.01-2.94) in the Iranian population. This association is in part similar to that found for the Greek population.

A comprehensive study to determine heterogeneity of autosomal recessive nonsyndromic hearing loss in Iran.

Hearing loss is the most common sensory disorder worldwide and affects 1 of every 500 newborns. In developed countries, at least 50% of cases are genetic, most often resulting in nonsyndromic deafness (70%), which is usually autosomal recessive (∼80%). Although the cause of hearing loss is heterogeneous, mutations in GJB2 gene at DFNB1 locus are the major cause of autosomal recessive nonsyndromic hearing loss (ARNSHL) in many populations. Our previous study showed that mutations of GJB2 gene do not contribute to the major genetic load of deafness in the Iranian population (∼16%). Therefore, to define the importance of other genes in contributing to an ARNSHL phenotype in the Iranian population, we used homozygosity mapping to identify regions of autozygosity-by-descent in 144 families which two or more progeny had ARNSHL but were negative for GJB2 gene mutations. Using flanking or intragenic short-tandem repeat markers for 33 loci we identified 33 different homozygous variations in 10 genes, of which 9 are novel. In aggregate, these data explain ∼40% of genetic background of ARNHSL in the Iranian population.