Ultra-Rare Variants Identify Biological Pathways and Candidate Genes in the Pathobiology of Non-Syndromic Cleft Palate Only.

In recent decades, many efforts have been made to elucidate the genetic causes of non-syndromic cleft palate (nsCPO), a complex congenital disease caused by the interaction of several genetic and environmental factors. Since genome-wide association studies have evidenced a minor contribution of common polymorphisms in nsCPO inheritance, we used whole exome sequencing data to explore the role of ultra-rare variants in this study. In a cohort of 35 nsCPO cases and 38 controls, we performed a gene set enrichment analysis (GSEA) and a hypergeometric test for assessing significant overlap between genes implicated in nsCPO pathobiology and genes enriched in ultra-rare variants in our cohort. GSEA highlighted an enrichment of ultra-rare variants in genes principally belonging to cytoskeletal protein binding pathway (Probability Density Function corrected p-value = 1.57 × 10-4); protein-containing complex binding pathway (p-value = 1.06 × 10-2); cell adhesion molecule binding pathway (p-value = 1.24 × 10-2); ECM-receptor interaction pathway (p-value = 1.69 × 10-2); and in the Integrin signaling pathway (p-value = 1.28 × 10-2). Two genes implicated in nsCPO pathobiology, namely COL2A1 and GLI3, ranked among the genes (n = 34) with nominal enrichment in the ultra-rare variant collapsing analysis (Fisher's exact test p-value < 0.05). These genes were also part of an independent list of genes highly relevant to nsCPO biology (n = 25). Significant overlap between the two sets of genes (hypergeometric test p-value = 5.86 × 10-3) indicated that enriched genes are likely to be implicated in physiological palate development and/or the pathological processes of oral clefting. In conclusion, ultra-rare variants collectively impinge on biological pathways crucial to nsCPO pathobiology and point to candidate genes that may contribute to the individual risk of disease. Sequencing can be an effective approach to identify candidate genes and pathways for nsCPO.

Biallelic loss of LDB3 leads to a lethal pediatric dilated cardiomyopathy.

Autosomal dominant variants in LDB3 (also known as ZASP), encoding the PDZ-LIM domain-binding factor, have been linked to a late onset phenotype of cardiomyopathy and myofibrillar myopathy in humans. However, despite knockout mice displaying a much more severe phenotype with premature death, bi-allelic variants in LDB3 have not yet been reported. Here we identify biallelic loss-of-function variants in five unrelated cardiomyopathy families by next-generation sequencing. In the first family, we identified compound heterozygous LOF variants in LDB3 in a fetus with bilateral talipes and mild left cardiac ventricular enlargement. Ultra-structural examination revealed highly irregular Z-disc formation, and RNA analysis demonstrated little/no expression of LDB3 protein with a functional C-terminal LIM domain in muscle tissue from the affected fetus. In a second family, a homozygous LDB3 nonsense variant was identified in a young girl with severe early-onset dilated cardiomyopathy with left ventricular non-compaction; the same homozygous nonsense variant was identified in a third unrelated female infant with dilated cardiomyopathy. We further identified homozygous LDB3 frameshift variants in two unrelated probands diagnosed with cardiomegaly and severely reduced left ventricular ejection fraction. Our findings demonstrate that recessive LDB3 variants can lead to an early-onset severe human phenotype of cardiomyopathy and myopathy, reminiscent of the knockout mouse phenotype, and supporting a loss of function mechanism.

Clinical and genetic profile of children with unexplained intellectual disability/developmental delay and epilepsy.

OBJECTIVE: This study was conducted to evaluate the validity of performing whole exome sequencing in children with unexplained intellectual disability (ID), developmental delay (DD), and epilepsy. METHODS: We enrolled 61 Iranian children with unexplained DD/ID, and epilepsy with no etiologic diagnosis. 64 % of cases were male and 36 % were female, with a mean age of 6.2 years (range, 38 days to 15 years). Approximately 79 % of patients were born to consanguineous parents or had non-related parents from a highly inbred local region. Whole-exome sequencing analysis followed by Sanger sequencing was performed in all patients. RESULTS: Pathogenic/likely pathogenic variants were identified in 59% (36/61) of patients, consisting of 26 novel and 14 known alterations. Variants of unknown significance were observed in 6.5 % (4/61) of patients. Variants in 28 genes have not been previously reported in Iranian patients with ID. Several additional phenotypes, mostly microcephaly, were common in 57.4 % of cases. Additionally, epilepsy was refractory in 40 % of patients. Three groups of brain anomalies consisting of brain dysgenesis, brain atrophy, and leukodystrophy were identified in our cohort. Mutations in genes implicated in cellular metabolic pathways were the most common, followed by ion channel/ion transporter and transcription pathways. DISCUSSION: High-throughput DNA sequencing of the Iranian population with a high rate of parental consanguinity is a valuable strategy for identifying genetic etiology in children with unexplained ID/DD and epilepsy. Determining the genetic basis and most commonly involved pathways may help to identify novel genes and targeted antiepileptic treatments.

Autozygosity-driven genetic diagnosis in consanguineous families from Italy and the Greater Middle East.

Autozygosity-driven exome analysis has been shown effective for identification of genes underlying recessive diseases especially in countries of the so-called Greater Middle East (GME), where high consanguinity unravels the phenotypic effects of recessive alleles and large family sizes facilitate homozygosity mapping. In Italy, as in most European countries, consanguinity is estimated low. Nonetheless, consanguineous Italian families are not uncommon in publications of genetic findings and are often key to new associations of genes with rare diseases. We collected 52 patients from 47 consanguineous families with suspected recessive diseases, 29 originated in GME countries and 18 of Italian descent. We performed autozygosity-driven exome analysis by detecting long runs of homozygosity (ROHs > 1.5 Mb) and by prioritizing candidate clinical variants within. We identified a pathogenic synonymous variant that had been previously missed in NARS2 and we increased an initial high diagnostic rate (47%) to 55% by matchmaking our candidate genes and including in the analysis shorter ROHs that may also happen to be autozygous. GME and Italian families contributed to diagnostic yield comparably. We found no significant difference either in the extension of the autozygous genome, or in the distribution of candidate clinical variants between GME and Italian families, while we showed that the average autozygous genome was larger and the mean number of candidate clinical variants was significantly higher (p = 0.003) in mutation-positive than in mutation-negative individuals, suggesting that these features influence the likelihood that the disease is autozygosity-related. We highlight the utility of autozygosity-driven genomic analysis also in countries and/or communities, where consanguinity is not widespread cultural tradition.

ROCK1 is associated with non-syndromic cleft palate.

BACKGROUND: Craniofacial morphogenesis is the result of an intricate multistep network of tightly controlled spatial and temporal signalling that involves several molecules and transcription factors organized into highly coordinated pathways. Any alteration in even one step of this delicate process can lead to congenital malformations such as cleft palate. One of the first steps in embryonal orofacial development is the migration of cells from the neural crests to the branchial arches. Next, the cells have to proliferate, differentiate, move and connect to each other in order to correctly form the palate. Cell contraction, promoted by the interaction of non-muscle myosin II and actin A, is a crucial step in morphogenesis and is regulated by ROCK1 protein. METHODS: A family-based association study was carried out in order to verify whether or not genetic variants of ROCK1 were associated with non-syndromic cleft palate (nsCP). Two cohorts from Italy and Iran, a total of 189 nsCP cases and their parents were enrolled. RESULTS: The rs35996865-G allele was under-transmitted in cases of nsCP [P = .006, odds ratio (OR) = 0.63 (95% CI 0.45-0.88)]. CONCLUSION: This investigation reveals for the first time data supporting a role for ROCK1 in nsCP aetiology.

Non-syndromic cleft palate: Association analysis on three gene polymorphisms of the folate pathway in Asian and Italian populations.

Periconceptional folic acid supplementation can reduce the risk of inborn malformations, including orofacial clefts. Polymorphisms of MTHFR, TCN2, and CBS folate-related genes seem to modulate the risk of cleft lip with or without cleft palate (CL/P) in some populations. CL/P and cleft palate only (CPO) are different malformations that share several features and possibly etiological causes. In the present investigation, we conducted a family-based, candidate gene association study of non-syndromic CPO. Three single nucleotide polymorphisms, namely, rs1801133 of MTHFR, rs1801198 of TCN2, and rs4920037 of CBS, were investigated in a sample that included 129 Italian and 65 Asian families. No evidence of association between the three genotyped polymorphisms and CPO was found in the Italian and Asian cases, indeed the transmission disequilibrium test did not detect any asymmetry of transmission of alleles. This investigation, although with some limitation, further supports that CL/P and CPO diverge in their genetic background.

Possible effect of SNAIL family transcriptional repressor 1 polymorphisms in non-syndromic cleft lip with or without cleft palate.

OBJECTIVE: Orofacial development is a complex process subjected to failure impairing. Indeed, the cleft of the lip and/or of the palate is among the most frequent inborn malformations. The JARID2 gene has been suggested to be involved in non-syndromic cleft lip with or without cleft palate (nsCL/P) etiology. JARID2 interacts with the polycomb repressive complex 2 (PRC2) in regulating the expression patterns of developmental genes by modifying the chromatin state. MATERIALS AND METHODS: Genes coding for the PRC2 components, as well as other genes active in cell differentiation and embryonic development, were selected for a family-based association study to verify their involvement in nsCL/P. A total of 632 families from Italy and Asia participated to the study. RESULTS: Evidence of allelic association was found with polymorphisms of SNAI1; in particular, the rs16995010-G allele was undertransmitted to the nsCL/P cases [P = 0.004, odds ratio = 0.69 (95% C.I. 0.54-0.89)]. However, the adjusted significance value corrected for all the performed tests was P = 0.051. CONCLUSIONS: The findings emerging by the present study suggest for the first time an involvement of SNAI1 in the nsCL/P onset. CLINICAL RELEVANCE: Interestingly, SNAI1 is known to promote epithelial to mesenchymal transition by repressing E-cadherin expression, but it needs an intact PRC2 to act this function. Alterations of this process could contribute to the complex etiology of nsCL/P.

Consanguineous marriages in the genetic counseling centers of Isfahan and the ethical issues of clinical consultations.

Consanguineous marriage, which is common in many regions in the world, has absorbed much attention as a causative factor in raising the incidence of genetic diseases. The adverse effects may be attributed to the expression of the genes received from common ancestors and mortality and morbidity of the offspring. Iran has a high rate of consanguineous marriages. In recent years genetic counseling has come to be considered in health care services. This cross-sectional study was conducted in order to determine the prevalence and types of consanguineous marriages in the genetic clinics in Isfahan. We aimed to define the different types of marriages, specific categories of genetic disorders associated with consanguineous marriages, and mode of inheritance in the family tree. We also narratively reviewed the ethical aspects of the issue. The data were collected using a simple questionnaire. A total number of 1535 couples from urban and rural areas formed the study population. The marriages were classified according to the degree of the relationship between couples, including: double cousin, first cousin, first cousin once removed, second cousin and beyond second cousin. The SPSS software version 16 was used for data analysis. Data obtained through genetic counseling offered during a 5-year period revealed that 74.3% had consanguineous relationships, 62.3% were first cousins, 1% were double cousins and 7.8% were second cousins. In addition, 76% of the couples had at least one genetic disease in their family tree. Related ethical issues were also considered in this study, including autonomy and informed decision making, benefit and harm assessment, confidentiality, ethics in research, justice in access to counseling services, financial problems ethics, and the intellectual property of scientific success.

Reciprocal 22q11.2 Deletion and Duplication in Siblings with Karyotypically Normal Parents.

The 22q11.2 locus is known to harbor a high risk for structural variation caused by non-allelic homologous recombination, resulting in deletions and duplications. Here, we describe the first family with one sibling carrying the 22q11 deletion and the other carrying the reciprocal duplication. FISH and SNP array analysis of the parents show a maternal origin for both deletion and duplication, without indications of balanced deletions/duplications or mosaicism. We hypothesize that germline mosaicism in the mother underlies the deletion and duplication, which would implicate a high recurrence risk for her offspring.

Prevalence of 22q11.2 microdeletion syndrome in Iranian patients with cleft palate.

BACKGROUND: 22q11.2 microdeletion syndrome is the most common multiple genetic disorder associated with learning disabilities, developmental delays, immune deficiency, hypocalcemia, and cleft palate. Finding some valid criteria for screening of 22q11.2 deletion syndromes in infants would be very helpful in early diagnosis and treatment. MATERIALS AND METHODS: Since 69% of individuals with 22q11.2 deletion have a palatal abnormality, we studied the prevalence of 22q11.2 deletion syndrome in 378 Iranian patients during a 5-year period, including 291 patients affected with cleft palate only without cleft lip (CPO) and 87 patients affected with velopharyngeal incompetence (VPI) and/or submucous cleft palate (SMCP). DNA copy number was analyzed with multiplex ligation-dependent probe amplification (MLPA) technique. RESULTS: In our study, 15/378 (3.97%) patients with palatal anomalies showed 22q11.2 deletion. Interestingly, this prevalence between syndromic patients was 15/104 (14.42%). CONCLUSION: It seems that SMCP or VPI, in addition to one or more another features of 22q11.2 deletions, especially developmental delay, may be good criteria for molecular investigation of 22q11.2 region.

Non-syndromic cleft lip with or without cleft palate in Asian populations: Association analysis on three gene polymorphisms of the folate pathway.

OBJECTIVE: Orofacial clefts (OFCs) are one of the most common birth defects in humans. They are the subject of a number of investigations aimed at elucidating the bases of their complex mode of inheritance involving both genetic and environmental factors. Genes belonging to the folate pathway have been among the most studied. The aim of the investigation was to replicate previous studies reporting evidence of association between polymorphisms of folate related genes and the occurrence of non-syndromic cleft lip with or without cleft palate (NSCL/P), using three independent samples of different ancestry: from Tibet, Bangladesh and Iran, respectively. DESIGN: Specifically, the polymorphisms rs1801133 of MTHFR, rs1801198 of TCN2, and rs4920037 of CBS, were tested. RESULTS: A decreased risk of NSCL/P was observed in patients presenting the C677T variant at MTHFR gene (relative risk for heterozygotes=0.53; 95% confidence interval [C.I.]=0.32-0.87). The investigated polymorphisms mapping at TCN2 and CBS genes did not provide any evidence of association. CONCLUSION: Overall, these results indicate that NSCL/P risk factors differ among populations and confirm the importance of testing putative susceptibility variants in different genetic backgrounds.

Identification of Proximal and Distal 22q11.2 Microduplications among Patients with Cleft Lip and/or Palate: A Novel Inherited Atypical 0.6 Mb Duplication.

Misalignments of low-copy repeats (LCRs) located in chromosome 22, particularly band 22q11.2, predispose to rearrangements. A variety of phenotypic features are associated with 22q11.2 microduplication syndrome which makes it challenging for the genetic counselors to recommend appropriate genetic assessment and counseling for the patients. In this study, multiplex ligation probe dependent amplification (MLPA) analysis was performed on 378 patients with cleft lip and/or palate to characterize rearrangements in patients suspected of 22q11.2 microduplication and microdeletion syndromes. Of 378 cases, 15 were diagnosed with a microdeletion with various sizes and 3 with duplications. For the first time in this study an atypical 0.6 Mb duplication is reported. Illustration of the phenotypes associated with the microduplications increases the knowledge of phenotypes reported in the literature.

Family-based association analysis between nonsyndromic cleft lip with or without cleft palate and IRF6 polymorphism in an Iranian population.

OBJECTIVES: Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is a common birth defect which is strongly associated with genetic factors. Previous studies in several populations showed a significant correlation between IRF6 rs642961 polymorphism and NSCL/P. The aim of this study is to indicate the correlation of IRF6 rs642961 polymorphism and NSCL/P in Iranian NSCL/P families. MATERIAL AND METHODS: In this study, we analyzed IRF6 rs642961 genotype in 352 individuals from 102 Iranian nuclear families affected by NSCL/P using iPlex assay on a Sequenom MassARRAY platform. Hardy-Weinberg equilibrium and Mendelian error checking were performed by Haploview 4.2. Allelic association analysis was conducted with family-based association tests implemented in FBAT program v2.03. RESULTS: The family-based association analysis revealed no significant association between IRF6 rs642961 genotypes and an increased NSCL/P risk. CONCLUSIONS: In contrast to other Asian populations, our study indicates that the IRF6 rs642961 polymorphism cannot be a risk factor for NSCL/P in an Iranian population. CLINICAL RELEVANCE: Genetic factors have an important role in NSCL/P, among which interferon regulatory factor 6 (IRF6) has been reported as a risk factor for NSCL/P in several populations; however, our data indicated no significant association between IRF6 polymorphism and NSCL/P in an Iranian population.

Detection of intragenic SMN1 mutations in spinal muscular atrophy patients with a single copy of SMN1.

Proximal spinal muscular atrophy is an autosomal recessive disorder characterized by symmetrical muscle weakness due to degeneration of alpha motor neurons in the spinal cord. Homozygous deletions in the SMN1 have been reported in more than 90% of spinal muscular atrophy cases. Compound heterozygous patients account for approximately 4% of spinal muscular atrophy cases. In this study, we performed a quantitative test in 20 of 87 spinal muscular atrophy patients who did not have homozygous deletion of SMN1. Mutation screening of SMN1 gene was performed in 4 patients who have only 1 copy of SMN1 to identify intragenic mutations. In addition to a previously described missense mutation in exon 4 (p.A188S/ c.562G>T), we identified 2 novel mutations including a single nucleotide insertion in exon 7 (c.861_862insT/p.R288X) and a deletion of nucleotide G in exon 3 (c.286delG/p.D96Tfs*53). Our results suggested that about 4% of spinal muscular atrophy patients have subtle mutations and might be considered in laboratory examination.

Evaluation of multiplex ligation-dependent probe amplification analysis versus multiplex polymerase chain reaction assays in the detection of dystrophin gene rearrangements in an Iranian population subset.

BACKGROUND: The Duchenne muscular dystrophy (DMD) gene is located in the short arm of the X chromosome (Xp21). It spans 2.4 Mb of the human genomic DNA and is composed of 79 exons. Mutations in the Dystrophin gene result in DMD and Becker muscular dystrophy. In this study, the efficiency of multiplex ligation-dependent probe amplification (MLPA) over multiplex polymerase chain reaction (PCR) assays in an Iranian population was investigated. MATERIALS AND METHODS: Multiplex PCR assays and MLPA analysis were carried out in 74 patients affected with DMD. RESULTS: Multiplex PCR detected deletions in 51% of the patients with DMD. MLPA analysis could determine all the deletions detected by the multiplex PCR. Additionally, MLPA was able to identify one more deletion and duplication in patients without detectable mutations by multiplex PCR. Moreover, MLPA precisely determined the exact size of the deletions. CONCLUSION: Although MLPA analysis is more sensitive for detection of deletions and duplications in the dystrophin gene, multiplex PCR might be used for the initial analysis of the boys affected with DMD in the Iranian population as it was able to detect 95% of the rearrangements in patients with DMD.

Genotype-phenotype correlation of survival motor neuron and neuronal apoptosis inhibitory protein genes in spinal muscular atrophy patients from Iran.

BACKGROUND: Proximal spinal muscular atrophy (SMA) is an autosomal recessive neuromuscular disease characterized by symmetrical proximal muscle weakness and atrophy. According to the severity of the disease and the age of onset, SMA can be divided into three groups. The survival motor neuron (SMN) gene that is located on 5q13 is identified as the disease determining gene. Another gene in this region is neuronal apoptosis inhibitory protein (NAIP), and its functional role in the pathogenesis of SMA has not been fully elucidated. Here, we investigated the correlation between deletions in SMN and NAIP genes with clinical features of SMA patients. MATERIALS AND METHODS: In the current study, 71 unrelated Iranian patients were investigated for the detection of deletions in SMN1 and NAIP genes. Polymerase chain reaction (PCR) was used to detect the deletions of exon 4 and 5 of the NAIP gene. Deletions in exon 7 and 8 of SMN1 gene were detected by RFLP-PCR with DraI and DdeI, respectively. RESULTS: Our results showed that 51 patients have homozygous deletions in SMN1 and/or NAIP genes. Among these 51 patients, deletion in NAIP gene were found in 35 patients (65.7% of type I, 22.5% type II and 11.42% type III). CONCLUSION: Defect in SMN1 gene plays a major role in manifesting of the disease and NAIP (4 and 5) gene acts as a modifying factor in severity of symptoms. Correlation between NAIP gene defect and severity of the disease is confirmed. However, the exact role of NAIP gene in SMA has yet to be fully clarified.

Use of in silico tools for classification of novel missense mutations identified in dystrophin gene in developing countries.

DMD gene which is composed of 79 exons is the largest known gene located on X chromosome (Xp21). Point mutations in the dystrophin gene are responsible for 30-35% of cases with DMD/BMD. Mutation analysis of all the exons of the DMD gene is costly in developing countries, therefore, a few of the exons are selected to be analyzed routinely in clinical laboratories. In this study, direct sequencing was used for detection of point mutations in 10 exons of dystrophin gene in patients affected with DMD without detectable large rearrangements. Freely available programs were used to predict the damaging effects of the mutations. Point mutations were successfully detected in three patients. Three novel mutations, two missense mutations located on nonconservative domains and a single nucleotide deletion, were detected. Missense mutations were predicted to change splicing efficiency. Detection of point mutations by DNA analysis followed by prediction of the pathogenecity by using bioinformatic tool might be an asset to provide proper diagnosis or genetic counseling to patients and their family.

An insight into genetics of non-syndromic cleft palate.

Our proband is a 29-year-old man, who is affected with soft cleft palate and hypernasality. A study of about six generations of this family pedigree shows that cleft palate has repeatedly occurred in males, with probably a X-linked recessive pattern of inheritance. Interestingly, the sister of the proband is affected with hypernasality and she has an affected son. This is the first report of X-linked inheritance pattern of cleft palate in Iran.

Mutation analysis of CACNA1A gene in Iranian migrainous and review literatures.

BACKGROUND: There are contrary results about the role of CACNA1A gene in the causation of common migraine in different populations. However, migraine may be genetically heterogeneous and more studies in different families and populations are required for a definite conclusion. The aim of this study was to surveyed leukocyte genomic DNA mutation of CACNA1A in Iranian migraine patients with [MA] and without aura [MO] who has family history of migraine and we performed a narrative review of all studies that evaluated CACNA1A gene, non-hemiplegic migraine [MA and MO] and FHM [familial hemiplegic migraine]. MATERIALS AND METHODS: The 30 patients with family history of migraine were selected for mutations analysis for CACNA1A gene by PCR method. For review, we searched MEDLINE-PUBMED, ISI, Scopus and Cochrane databases up to December 2012. RESULTS: Mutation analysis of the 4 exons of the CACNA1A gene in these patients revealed no mutations in this gene. Direct sequencing revealed a polymorphism previously reported G to A transition in the exon 16 [nt2369, G→A] in 9 patients. In review, the correlation of FHM loci [CACNA1A gene] with MA and MO has been showed in different population and only small population from Caucasians presented this correlation. CONCLUSION: CACNA1A is most likely not a major susceptibility gene for common migraine in Iranian maigrainous. It's essential to study more on larger series and covering all 47 exons of the CACNA1A gene to confirm this hypothesis.

Cutis laxa type II with mutation in the pyrroline-5-carboxylate reductase 1 gene.

A 14-year-old Iranian boy with congenital cutis laxa and several other typical autosomal recessive type II features was examined. Mutation analysis of the pyrroline-5-carboxylate reductase 1 gene revealed a single-base deletion (c.345delC) in exon 4 leading to frame shift and premature termination of translation.