Biallelic variants identified in 36 Pakistani families and trios with autism spectrum disorder.

With its high rate of consanguineous marriages and diverse ethnic population, little is currently understood about the genetic architecture of autism spectrum disorder (ASD) in Pakistan. Pakistan has a highly ethnically diverse population, yet with a high proportion of endogamous marriages, and is therefore anticipated to be enriched for biallelic disease-relate variants. Here, we attempt to determine the underlying genetic abnormalities causing ASD in thirty-six small simplex or multiplex families from Pakistan. Microarray genotyping followed by homozygosity mapping, copy number variation analysis, and whole exome sequencing were used to identify candidate. Given the high levels of consanguineous marriages among these families, autosomal recessively inherited variants were prioritized, however de novo/dominant and X-linked variants were also identified. The selected variants were validated using Sanger sequencing. Here we report the identification of sixteen rare or novel coding variants in fifteen genes (ARAP1, CDKL5, CSMD2, EFCAB12, EIF3H, GML, NEDD4, PDZD4, POLR3G, SLC35A2, TMEM214, TMEM232, TRANK1, TTC19, and ZNF292) in affected members in eight of the families, including ten homozygous variants in four families (nine missense, one loss of function). Three heterozygous de novo mutations were also identified (in ARAP1, CSMD2, and NEDD4), and variants in known X-linked neurodevelopmental disorder genes CDKL5 and SLC35A2. The current study offers information on the genetic variability associated with ASD in Pakistan, and demonstrates a marked enrichment for biallelic variants over that reported in outbreeding populations. This information will be useful for improving approaches for studying ASD in populations where endogamy is commonly practiced.

Biallelic Loss of Function Mutation in Sodium Channel Gene SCN10A in an Autism Spectrum Disorder Trio from Pakistan.

The genetic dissection of autism spectrum disorders (ASD) has uncovered the contribution of de novo mutations in many single genes as well as de novo copy number variants. More recent work also suggests a strong contribution from recessively inherited variants, particularly in populations in which consanguineous marriages are common. What is also becoming more apparent is the degree of pleiotropy, whereby mutations in the same gene may have quite different phenotypic and clinical consequences. We performed whole exome sequencing in a group of 115 trios from countries with a high level of consanguineous marriages. In this paper we report genetic and clinical findings on a proband with ASD, who inherited a biallelic truncating pathogenic/likely pathogenic variant in the gene encoding voltage-gated sodium channel X alpha subunit, SCN10A (NM_006514.2:c.937G>T:(p.Gly313*)). The biallelic pathogenic/likely pathogenic variant in this study have different clinical features than heterozygous mutations in the same gene. The study of consanguineous families for autism spectrum disorder is highly valuable.

Genetic studies of multiple consanguineous Pakistani families segregating oculocutaneous albinism identified novel and reported mutations.

Oculocutaneous albinism (OCA) is an autosomal-recessive disorder of a defective melanin pathway. The condition is characterized by hypopigmentation of hair, dermis, and ocular tissue. Genetic studies have reported seven nonsyndromic OCA genes, among which Pakistani OCA families mostly segregate TYR and OCA2 gene mutations. Here in the present study, we investigate the genetic factors of eight consanguineous OCA families from Pakistan. Genetic analysis was performed through single-nucleotide polymorphism (SNP) genotyping (for homozygosity mapping), whole exome sequencing (for mutation identification), Sanger sequencing (for validation and segregation analysis), and quantitative PCR (qPCR) (for copy number variant [CNV] validation). Genetic mapping in one family identified a novel homozygous deletion mutation of the entire TYRP1 gene, and a novel deletion of exon 19 in the OCA2 gene in two apparently unrelated families. In three further families, we identified homozygous mutations in TYR (NM_000372.4:c.1424G > A; p.Trp475*), NM_000372.4:c.895C > T; p.Arg299Cys), and SLC45A2 (NM_016180:c.1532C > T; p.Ala511Val). For the remaining two families, G and H, compound heterozygous TYR variants NM_000372.4:c.1037-7T > A, NM_000372.4:c.1255G > A (p.Gly419Arg), and NM_000372.4:c.1255G > A (p.Gly419Arg) and novel variant NM_000372.4:c.248T > G; (p.Val83Gly), respectively, were found. Our study further extends the evidence of TYR and OCA2 as genetic mutation hot spots in Pakistani families. Genetic screening of additional OCA cases may also contribute toward the development of Pakistani specific molecular diagnostic tests, genetic counseling, and personalized healthcare.

Three Mutations in the Bilateral Frontoparietal Polymicrogyria Gene GPR56 in Pakistani Intellectual Disability Families.

Bilateral frontoparietal polymicrogyria (BFPP, MIM 606854) is a heterogeneous autosomal recessive disorder of abnormal cortical lamination, leading to moderate-to-severe intellectual disability (ID), seizure disorder, and motor difficulties, and caused by mutations in the G protein-coupled receptor 56 ( GPR56 ) gene. Twenty-eight mutations in 40 different families have been reported in the literature. The clinical and neuroimaging phenotype is consistent in these cases. The BFPP cortex consists of numerous small gyral cells, with scalloping of the cortical-white matter junction. There are also associated white matter, brain stem, and cerebellar changes. GPR56 is a member of an adhesion G protein-coupled receptor family with a very long N-terminal stalk and seven transmembrane domains. In this study, we identified three families from Pakistan, ascertained primarily for ID, with overlapping approximately 1 Mb region (chr16:56,973,335-57,942,866) of homozygosity by descent, including 24 RefSeq genes. We found three GPR56 homozygous mutations, using next-generation sequencing. These mutations include a substitutional variant, c.1460T > C; p.L487P, (chr16:57693480 T > C), a 13-bp insertion causing the frameshift and truncating mutation, p.Leu269Hisfs*21 (NM_005682.6:c.803_804insCCATGGAGGTGCT; Chr16: 57689345_57689346insCCATGGAGGTGCT), and a truncating mutation c.1426C > T; p.Arg476* (Chr16:57693446C > T). These mutations fully segregated with ID in these families and were absent in the Exome Aggregation Consortium database that has approximately 8,000 control samples of South Asian origin. Two of these mutations have been reported in ClinVar database, and the third one has not been reported before. Three families from Pakistan with GPR56 mutations have been reported before. With the addition of our findings, the total number of mutations reported in Pakistani patients now is six. These results increase our knowledge regarding the mutational spectrum of the GPR56 gene causing BFPP/ID.

Identification of a homozygous missense mutation in LRP2 and a hemizygous missense mutation in TSPYL2 in a family with mild intellectual disability.

Non-syndromic autosomal recessive intellectual disability (ID) is a genetically heterogeneous disorder with more than 50 mutated genes to date. ID is characterized by deficits in memory skills and language development with difficulty in learning, problem solving, and adaptive behaviors, and affects ∼ 1% of the population. For detection of disease-causing mutations in such a heterogeneous disorder, homozygosity mapping together with exome sequencing is a powerful approach, as almost all known genes can be assessed simultaneously in a high-throughput manner. In this study, a hemizygous c.786C>G:p.Ile262Met in the testis specific protein Y-encoded-like 2 (TSPYL2) gene and a homozygous c.11335G>A:p.Asp3779Asn in the low-density lipoprotein receptor-related protein 2 (LRP2) gene were detected after genome-wide genotyping and exome sequencing in a consanguineous Pakistani family with two boys with mild ID. Mutations in the LRP2 gene have previously been reported in patients with Donnai-Barrow and Stickler syndromes. LRP2 has also been associated with a 2q locus for autism (AUTS5). The TSPYL2 variant is not listed in any single-nucleotide polymorphism databases, and the LRP2 variant was absent in 400 ethnically matched healthy control chromosomes, and is not listed in single-nucleotide polymorphism databases as a common polymorphism. The LRP2 mutation identified here is located in one of the low-density lipoprotein-receptor class A domains, which is a cysteine-rich repeat that plays a central role in mammalian cholesterol metabolism, suggesting that alteration of cholesterol processing pathway can contribute to ID.

Mutations in the histamine N-methyltransferase gene, HNMT, are associated with nonsyndromic autosomal recessive intellectual disability.

Histamine (HA) acts as a neurotransmitter in the brain, which participates in the regulation of many biological processes including inflammation, gastric acid secretion and neuromodulation. The enzyme histamine N-methyltransferase (HNMT) inactivates HA by transferring a methyl group from S-adenosyl-l-methionine to HA, and is the only well-known pathway for termination of neurotransmission actions of HA in mammalian central nervous system. We performed autozygosity mapping followed by targeted exome sequencing and identified two homozygous HNMT alterations, p.Gly60Asp and p.Leu208Pro, in patients affected with nonsyndromic autosomal recessive intellectual disability from two unrelated consanguineous families of Turkish and Kurdish ancestry, respectively. We verified the complete absence of a functional HNMT in patients using in vitro toxicology assay. Using mutant and wild-type DNA constructs as well as in silico protein modeling, we confirmed that p.Gly60Asp disrupts the enzymatic activity of the protein, and that p.Leu208Pro results in reduced protein stability, resulting in decreased HA inactivation. Our results highlight the importance of inclusion of HNMT for genetic testing of individuals presenting with intellectual disability.

New recessive truncating mutation in LTBP3 in a family with oligodontia, short stature, and mitral valve prolapse.

Latent TGFB-binding protein 3 (LTBP3) is known to increase bio-availability of TGFB. A homozygous mutation in this gene has previously been associated with oligodontia and short stature in a single family. We report on two sisters with homozygous truncating mutations in LTBP3. In addition to oligodontia and short stature, both sisters have mitral valve prolapse, suggesting a link between truncating LTBP3 mutations and mitral valve disease mediated through the TGFB pathway.

Mutations in DCPS and EDC3 in autosomal recessive intellectual disability indicate a crucial role for mRNA decapping in neurodevelopment.

There are two known mRNA degradation pathways, 3' to 5' and 5' to 3'. We identified likely pathogenic variants in two genes involved in these two pathways in individuals with intellectual disability. In a large family with multiple branches, we identified biallelic variants in DCPS in three affected individuals; a splice site variant (c.636+1G>A) that results in an in-frame insertion of 45 nucleotides and a missense variant (c.947C>T; p.Thr316Met). DCPS decaps the cap structure generated by 3' to 5' exonucleolytic degradation of mRNA. In vitro decapping assays showed an ablation of decapping function for both variants in DCPS. In another family, we identified a homozygous mutation (c.161T>C; p.Phe54Ser) in EDC3 in two affected children. EDC3 stimulates DCP2, which decaps mRNAs at the beginning of the 5' to 3' degradation pathway. In vitro decapping assays showed that altered EDC3 is unable to enhance DCP2 decapping at low concentrations and even inhibits DCP2 decapping at high concentration. We show that individuals with biallelic mutations in these genes of seemingly central functions are viable and that these possibly lead to impairment of neurological functions linking mRNA decapping to normal cognition. Our results further affirm an emerging theme linking aberrant mRNA metabolism to neurological defects.

Biallelic truncating mutations in FMN2, encoding the actin-regulatory protein Formin 2, cause nonsyndromic autosomal-recessive intellectual disability.

Dendritic spines represent the major site of neuronal activity in the brain; they serve as the receiving point for neurotransmitters and undergo rapid activity-dependent morphological changes that correlate with learning and memory. Using a combination of homozygosity mapping and next-generation sequencing in two consanguineous families affected by nonsyndromic autosomal-recessive intellectual disability, we identified truncating mutations in formin 2 (FMN2), encoding a protein that belongs to the formin family of actin cytoskeleton nucleation factors and is highly expressed in the maturing brain. We found that FMN2 localizes to punctae along dendrites and that germline inactivation of mouse Fmn2 resulted in animals with decreased spine density; such mice were previously demonstrated to have a conditioned fear-learning defect. Furthermore, patient neural cells derived from induced pluripotent stem cells showed correlated decreased synaptic density. Thus, FMN2 mutations link intellectual disability either directly or indirectly to the regulation of actin-mediated synaptic spine density.

Identification of a homozygous splice site mutation in the dynein axonemal light chain 4 gene on 22q13.1 in a large consanguineous family from Pakistan with congenital mirror movement disorder.

Mirror movements (MRMV) are involuntary movements on one side of the body that mirror voluntary movements on the opposite side. Congenital mirror movement disorder is a rare, typically autosomal-dominant disorder, although it has been suspected that some sporadic cases may be due to recessive inheritance. Using a linkage analysis and a candidate gene approach, two genes have been implicated in congenital MRMV disorder to date: DCC on 18q21.2 (MRMV1), which encodes a netrin receptor, and RAD51 on 15q15.1 (MRMV2), which is involved in the maintenance of genomic integrity. Here, we describe a large consanguineous Pakistani family with 11 cases of congenital MRMV disorder reported across five generations, with autosomal recessive inheritance likely. Sanger sequencing of DCC and RAD51 did not identify a mutation. We then employed microarray genotyping and autozygosity mapping to identify a shared region of homozygosity-by-descent among the affected individuals. We identified a large autozygous region of ~3.3 Mb on chromosome 22q13.1 (Chr22:36605976-39904648). We used Sanger sequencing to exclude several candidate genes within this region, including DMC1 and NPTXR. Whole exome sequencing was employed, and identified a splice site mutation in the dynein axonemal light chain 4 gene, DNAL4. This splice site change leads to skipping of exon 3, and omission of 28 amino acids from DNAL4 protein. Linkage analysis using Simwalk2 gives a maximum Lod score of 6.197 at this locus. Whether or how DNAL4 function may relate to the function of DCC or RAD51 is not known. Also, there is no suggestion of primary ciliary dyskinesis, situs inversus, or defective sperm in affected family members, which might be anticipated given a putative role for DNAL4 in axonemal-based dynein complexes. We suggest that DNAL4 plays a role in the cytoplasmic dynein complex for netrin-1-directed retrograde transport, and in commissural neurons of the corpus callosum in particular. This, in turn, could lead to faulty cross-brain wiring, resulting in MRMV.

Over-expression of either MECP2_e1 or MECP2_e2 in neuronally differentiated cells results in different patterns of gene expression.

Mutations in MECP2 are responsible for the majority of Rett syndrome cases. MECP2 is a regulator of transcription, and has two isoforms, MECP2_e1 and MECP2_e2. There is accumulating evidence that MECP2_e1 is the etiologically relevant variant for Rett. In this study we aim to detect genes that are differentially transcribed in neuronal cells over-expressing either of these two MECP2 isoforms. The human neuroblastoma cell line SK-N-SH was stably infected by lentiviral vectors over-expressing MECP2_e1, MECP2_e2, or eGFP, and were then differentiated into neurons. The same lentiviral constructs were also used to infect mouse Mecp2 knockout (Mecp2(tm1.1Bird)) fibroblasts. RNA from these cells was used for microarray gene expression analysis. For the human neuronal cells, ∼ 800 genes showed >three-fold change in expression level with the MECP2_e1 construct, and ∼ 230 with MECP2_e2 (unpaired t-test, uncorrected p value <0.05). We used quantitative RT-PCR to verify microarray results for 41 of these genes. We found significant up-regulation of several genes resulting from over-expression of MECP2_e1 including SRPX2, NAV3, NPY1R, SYN3, and SEMA3D. DOCK8 was shown via microarray and qRT-PCR to be upregulated in both SK-N-SH cells and mouse fibroblasts. Both isoforms up-regulated GABRA2, KCNA1, FOXG1 and FOXP2. Down-regulation of expression in the presence of MECP2_e1 was seen with UNC5C and RPH3A. Understanding the biology of these differentially transcribed genes and their role in neurodevelopment may help us to understand the relative functions of the two MECP2 isoforms, and ultimately develop a better understanding of RTT etiology and determine the clinical relevance of isoform-specific mutations.

Truncation of the E3 ubiquitin ligase component FBXO31 causes non-syndromic autosomal recessive intellectual disability in a Pakistani family.

In this study, we have performed autozygosity mapping on a large consanguineous Pakistani family segregating with intellectual disability. We identified two large regions of homozygosity-by-descent (HBD) on 16q12.2-q21 and 16q24.1-q24.3. Whole exome sequencing (WES) was performed on an affected individual from the family, but initially, no obvious mutation was detected. However, three genes within the HBD regions that were not fully captured during the WES were Sanger sequenced and we identified a five base pair deletion (actually six base pairs deleted plus one base pair inserted) in exon 7 of the gene FBXO31. The variant segregated completely in the family, in recessive fashion giving a LOD score of 3.95. This variant leads to a frameshift and a premature stop codon and truncation of the FBXO31 protein, p.(Cys283Asnfs*81). Quantification of mRNA and protein expression suggests that nonsense-mediated mRNA decay also contributes to the loss of FBXO31 protein in affected individuals. FBXO31 functions as a centrosomal E3 ubiquitin ligase, in association with SKP1 and Cullin-1, involved in ubiquitination of proteins targeted for degradation. The FBXO31/SKP1/Cullin1 complex is important for neuronal morphogenesis and axonal identity. FBXO31 also plays a role in dendrite growth and neuronal migration in developing cerebellar cortex. Our finding adds further evidence of the involvement of disruption of the protein ubiquitination pathway in intellectual disability.

Identification of risk genes for autism spectrum disorder through copy number variation analysis in Austrian families.

Autism or autism spectrum disorder (ASD) is a range of neurodevelopmental disorders starting in early childhood and is characterized by impairments in communication and reciprocal social interaction and presence of restricted and repetitive patterns of behavior. The contribution of genetic factors to autism is clear in twin and family studies. It is apparent that, overall, ASD is a complex non-Mendelian disorder. Recent studies suggest that copy number variations (CNVs) play a significant role in the etiology of ASD. For the current work, we recruited 245 family members from 73 ASD families from Styria, Austria. The DNA from probands was genotyped with Affymetrix single nucleotide polymorphism (SNP) 6.0 microarrays to screen for CNVs in their genomes. Analysis of the microarray data was performed using three different algorithms, and a list of stringent calls was compared to existing CNV data from over 2,357 controls of European ancestry. For stringent calls not present in controls, quantitative real-time PCR (qRT-PCR) was used to validate the CNVs in the probands and in their family members. Twenty-two CNVs were validated from this set (five of which are apparently de novo), many of which appear likely to disrupt genes that may be considered as good candidates for neuropsychiatric disorders, including DLG2, S100B, ARX, DIP2A, HPCAL1, and GPHN. Several others disrupt genes that have previously been implicated in autism, such as BDNF, AUTS2, DPP6, and C18orf22, and our data add to the growing evidence of their involvement in ASD.

CC2D2A, encoding a coiled-coil and C2 domain protein, causes autosomal-recessive mental retardation with retinitis pigmentosa.

Autosomal-recessive inheritance is believed to be relatively common in mental retardation (MR), although only four genes for nonsyndromic autosomal-recessive mental retardation (ARMR) have been reported. In this study, we ascertained a consanguineous Pakistani family with ARMR in four living individuals from three branches of the family, plus an additional affected individual later identified as a phenocopy. Retinitis pigmentosa was present in affected individuals, but no other features suggestive of a syndromic form of MR were found. We used Affymetrix 500K microarrays to perform homozygosity mapping and identified a homozygous and haploidentical region of 11.2 Mb on chromosome 4p15.33-p15.2. Linkage analysis across this region produced a maximum two-point LOD score of 3.59. We sequenced genes within the critical region and identified a homozygous splice-site mutation segregating in the family, within a coiled-coil and C2 domain-containing gene, CC2D2A. This mutation leads to the skipping of exon 19, resulting in a frameshift and a truncated protein lacking the C2 domain. Conservation analysis for CC2D2A suggests a functional domain near the C terminus as well as the C2 domain. Preliminary functional studies of CC2D2A suggest a possible role in Ca(2+)-dependent signal transduction. Identifying the function of CC2D2A, and a possible common pathway with CC2D1A, in correct neuronal development and functioning may help identify possible therapeutic targets for MR.