This study was undertaken to investigate the causal mutations responsible for autosomal recessive congenital stationary night blindness (CSNB) in consanguineous Pakistani families.

A challenge in clinical genomics is to predict whether copy number variation (CNV) affecting a gene or multiple genes will manifest as disease. Increasing recognition of gene dosage effects in neurodevelopmental disorders prompted us to develop a computational approach based on critical-exon (highly expressed in brain, highly conserved) examination for potential etiologic effects. Using a large CNV dataset, our updated analyses revealed significant (P < 1.64 × 10(-15)) enrichment of critical-exons within rare CNVs in cases compared to controls. Separately, we used a weighted gene co-expression network analysis (WGCNA) to construct an unbiased protein module from prenatal and adult tissues and found it significantly enriched for critical exons in prenatal (P < 1.15 × 10(-50), OR = 2.11) and adult (P < 6.03 × 10(-18), OR = 1.55) tissues. WGCNA yielded 1,206 proteins for which we prioritized the corresponding genes as likely to have a role in neurodevelopmental disorders. We compared the gene lists obtained from critical-exon and WGCNA analysis and found 438 candidate genes associated with CNVs annotated as pathogenic, or as variants of uncertain significance (VOUS), from among 10,619 developmental delay cases. We identified genes containing CNVs previously considered to be VOUS to be new candidate genes for neurodevelopmental disorders (GIT1, MVB12B and PPP1R9A) demonstrating the utility of this strategy to index the clinical effects of CNVs.

This study was performed to investigate the genetic determinants of autosomal recessive congenital cataracts in large consanguineous families.

Cohen Syndrome (COH1) is a rare autosomal recessive disorder, principally identified by ocular, neural and muscular deficits. We identified three large consanguineous Pakistani families with intellectual disability and in some cases with autistic traits.

This study was conducted to localize and identify causal mutations associated with autosomal recessive retinitis pigmentosa (RP) in consanguineous familial cases of Pakistani origin.

Next-generation sequencing (NGS) of exomes and genomes has accelerated the identification of genes involved in Mendelian phenotypes. However, many NGS studies fall short of identifying causal variants, with estimates for success rates as low as 25% for uncovering the pathological variant underlying disease etiology. An important reason for such failures is familial locus heterogeneity, where within a single pedigree causal variants in two or more genes underlie Mendelian trait etiology. As examples of intra- and inter-sibship familial locus heterogeneity, we present 10 consanguineous Pakistani families segregating hearing impairment due to homozygous variants in two different hearing impairment genes and a European-American pedigree in which hearing impairment is caused by four variants in three different genes. We have identified 41 additional pedigrees with syndromic and nonsyndromic hearing impairment for which a single previously reported hearing impairment gene has been identified but only segregates with the phenotype in a subset of affected pedigree members. We estimate that locus heterogeneity occurs in 15.3% (95% confidence interval: 11.9%, 19.9%) of the families in our collection. We demonstrate novel approaches to apply linkage analysis and homozygosity mapping (for autosomal recessive consanguineous pedigrees), which can be used to detect locus heterogeneity using either NGS or SNP array data. Results from linkage analysis and homozygosity mapping can also be used to group sibships or individuals most likely to be segregating the same causal variants and thereby increase the success rate of gene identification.

Ectodermal dysplasias (EDs) are a large heterogeneous group of inherited disorders exhibiting abnormalities in ectodermally derived appendages such as hair, nails, teeth and sweat glands. EDs associated with reticulated pigmentation phenotype are rare entities for which the genetic basis and pathophysiology are not well characterized. The present study describes a five generation consanguineous Pakistani family segregating an autosomal recessive form of a novel type of ectodermal dysplasia. The affected members present with sparse and woolly hair, severe nail dystrophy and reticulate skin pigmentation. After exclusion of known gene loci related with other skin disorders, genome-wide linkage analysis was performed using Illumina HumanOmniExpress beadchip SNP arrays. We linked this form of ED to human chromosome 18p11.32-p11.31 flanked by the SNPs rs9284390 (0.113Mb) and rs4797100 (3.14 Mb). A maximum two-point LOD score of 3.3 was obtained with several markers along the disease interval. The linkage interval of 3.03 Mb encompassed seventeen functional genes. However, sequence analysis of all these genes did not discover any potentially disease causing-variants. The identification of this novel locus provides additional information regarding the mapping of a rare form of ED. Further research, such as the use of whole-genome sequencing, would be expected to reveal any pathogenic mutation within the disease locus.

We describe six persons from three families with three homozygous protein truncating variants in PUS7: c.89_90del (p.Thr30Lysfs∗20), c.1348C>T (p.Arg450∗), and a deletion of the penultimate exon 15. All these individuals have intellectual disability with speech delay, short stature, microcephaly, and aggressive behavior. PUS7 encodes the RNA-independent pseudouridylate synthase 7. Pseudouridylation is the most abundant post-transcriptional modification in RNA, which is primarily thought to stabilize secondary structures of RNA. We show that the disease-related variants lead to abolishment of PUS7 activity on both tRNA and mRNA substrates. Moreover, pus7 knockout in Drosophila melanogaster results in a number of behavioral defects, including increased activity, disorientation, and aggressiveness supporting that neurological defects are caused by PUS7 variants. Our findings demonstrate that RNA pseudouridylation by PUS7 is essential for proper neuronal development and function.

Developmental eye defects often severely reduce vision. Despite extensive efforts, for a substantial fraction of these cases the molecular causes are unknown. Recessive eye disorders are frequent in consanguineous populations and such large families with multiple affected individuals provide an opportunity to identify recessive causative genes. We studied a Pakistani consanguineous family with three affected individuals with congenital vision loss and progressive eye degeneration. The family was analyzed by exome sequencing of one affected individual and genotyping of all family members. We have identified a non-synonymous homozygous variant (NM_001128918.2: c.1708C > G: p.Arg570Gly) in the MARK3 gene as the likely cause of the phenotype. Given that MARK3 is highly conserved in flies (I: 55%; S: 67%) we knocked down the MARK3 homologue, par-1, in the eye during development. This leads to a significant reduction in eye size, a severe loss of photoreceptors and loss of vision based on electroretinogram (ERG) recordings. Expression of the par-1 p.Arg792Gly mutation (equivalent to the MARK3 variant found in patients) in developing fly eyes also induces loss of eye tissue and reduces the ERG signals. The data in flies and human indicate that the MARK3 variant corresponds to a loss of function. We conclude that the identified mutation in MARK3 establishes a new gene-disease link, since it likely causes structural abnormalities during eye development and visual impairment in humans, and that the function of MARK3/par-1 is evolutionarily conserved in eye development.

Retinitis pigmentosa is an important cause of severe visual dysfunction. This study reports a novel splicing mutation in the lecithin retinol acyltransferase (LRAT) gene associated with early onset retinitis pigmentosa and characterizes the effects of this mutation on mRNA splicing and structure.

Inherited deafness is clinically and genetically heterogeneous. We recently mapped DFNB86, a locus associated with nonsyndromic deafness, to chromosome 16p. In this study, whole-exome sequencing was performed with genomic DNA from affected individuals from three large consanguineous families in which markers linked to DFNB86 segregate with profound deafness. Analyses of these data revealed homozygous mutation c.208G>T (p.Asp70Tyr) or c.878G>C (p.Arg293Pro) in TBC1D24 as the underlying cause of deafness in the three families. Sanger sequence analysis of TBC1D24 in an additional large family in which deafness segregates with DFNB86 identified the c.208G>T (p.Asp70Tyr) substitution. These mutations affect TBC1D24 amino acid residues that are conserved in orthologs ranging from fruit fly to human. Neither variant was observed in databases of single-nucleotide variants or in 634 chromosomes from ethnically matched control subjects. TBC1D24 in the mouse inner ear was immunolocalized predominantly to spiral ganglion neurons, indicating that DFNB86 deafness might be an auditory neuropathy spectrum disorder. Previously, six recessive mutations in TBC1D24 were reported to cause seizures (hearing loss was not reported) ranging in severity from epilepsy with otherwise normal development to epileptic encephalopathy resulting in childhood death. Two of our four families in which deafness segregates with mutant alleles of TBC1D24 were available for neurological examination. Cosegregation of epilepsy and deafness was not observed in these two families. Although the causal relationship between genotype and phenotype is not presently understood, our findings, combined with published data, indicate that recessive alleles of TBC1D24 can cause either epilepsy or nonsyndromic deafness.

Mesenchymal stem cells (MSCs) have the potential for treatment of diabetic cardiomyopathy; however, the repair capability of MSCs declines with age and disease. MSCs from diabetic animals exhibit impaired survival, proliferation, and differentiation and therefore require a strategy to improve their function. The aim of the study was to develop a preconditioning strategy to augment the ability of MSCs from diabetes patients to repair the diabetic heart.

Autism spectrum disorder (ASD) is a highly heritable disorder of complex and heterogeneous aetiology. It is primarily characterized by altered cognitive ability including impaired language and communication skills and fundamental deficits in social reciprocity. Despite some notable successes in neuropsychiatric genetics, overall, the high heritability of ASD (~90%) remains poorly explained by common genetic risk variants. However, recent studies suggest that rare genomic variation, in particular copy number variation, may account for a significant proportion of the genetic basis of ASD. We present a large scale analysis to identify candidate genes which may contain low-frequency recessive variation contributing to ASD while taking into account the potential contribution of population differences to the genetic heterogeneity of ASD. Our strategy, homozygous haplotype (HH) mapping, aims to detect homozygous segments of identical haplotype structure that are shared at a higher frequency amongst ASD patients compared to parental controls. The analysis was performed on 1,402 Autism Genome Project trios genotyped for 1 million single nucleotide polymorphisms (SNPs). We identified 25 known and 1,218 novel ASD candidate genes in the discovery analysis including CADM2, ABHD14A, CHRFAM7A, GRIK2, GRM3, EPHA3, FGF10, KCND2, PDZK1, IMMP2L and FOXP2. Furthermore, 10 of the previously reported ASD genes and 300 of the novel candidates identified in the discovery analysis were replicated in an independent sample of 1,182 trios. Our results demonstrate that regions of HH are significantly enriched for previously reported ASD candidate genes and the observed association is independent of gene size (odds ratio 2.10). Our findings highlight the applicability of HH mapping in complex disorders such as ASD and offer an alternative approach to the analysis of genome-wide association data.

Mesenchymal stem cell (MSC) transplantation has emerged as a promising therapy for liver fibrosis. Issues concerning poor MSC survival and engraftment in the fibrotic liver still persist and warrant development of a strategy to increase MSC potency for liver repair. The present study was designed to examine a synergistic role for Interleukin-6 (IL-6) and MSCs therapy in the recovery of carbon tetrachloride (CCl(4)) induced injured hepatocytes in vitro and in vivo.

Cosegregation of profound, congenital deafness with markers on chromosome 6q13 in three Pakistani families defines a new recessive deafness locus, DFNB37. Haplotype analyses reveal a 6-cM linkage region, flanked by markers D6S1282 and D6S1031, that includes the gene encoding unconventional myosin VI. In families with recessively inherited deafness, DFNB37, our sequence analyses of MYO6 reveal a frameshift mutation (36-37insT), a nonsense mutation (R1166X), and a missense mutation (E216V). These mutations, along with a previously published missense allele linked to autosomal dominant progressive hearing loss (DFNA22), provide an allelic spectrum that probes the relationship between myosin VI dysfunction and the resulting phenotype.

Previously, we have described the clinical and molecular characterization of a de novo 14q13.1-q21.1 microdeletion, less than 3.5 Mb in size, in a patient with severe microcephaly, psychomotor retardation, and other clinical anomalies. Here we report the characterization of the genomic structure of the human tuberin-like protein gene 1 (TULIP1; approved gene symbol GARNL1), a CpGisland-associated, brain-expressed candidate gene for the neurological findings in our patient, and its murine homologue. The human TULIP1 gene was mapped to chromosome band 14q13.2 by fluorescence in situ hybridization of BAC clone RP11-355C3 (GenBank Accession No. AL160231), containing the 3' region of the gene. TULIP1 spans about 271 kb of human genomic DNA and is divided into 41 exons. An untranscribed, processed pseudogene of TULIP1 was found on human chromosome band 9q31.1. The active locus TULIP1, encoding a predicted protein of 2036 amino acids, is expressed ubiquitously in pre- and postnatal human tissues. The murine homologue Tulip1 spans about 220 kb of mouse genomic DNA and is also divided into 41 exons, encoding a predicted protein of 2035 amino acids. No pseudogene could be found in the available mouse sequence data. Several splicing variants were found. Considering the location, expression profile, and predicted function, TULIP1 is a strong candidate for several neurological features seen in 14q deletion patients. Additionally we searched for mutations in the coding region of TULIP1 in subjects from a family with idiopathic basal ganglia calcification (IBGC; Fahr disease), previously linked to chromosome 14q. We identified two novel SNPs in the intron-exon boundaries; however, they did not segregate only with affected subjects in the predicted model of an autosomal dominant disease such as IBGC.

We identified overlapping homozygous regions within the DFNB25 locus in two Dutch and ten Pakistani families with sensorineural autosomal-recessive nonsyndromic hearing impairment (arNSHI). Only one of the families, W98-053, was not consanguineous, and its sibship pointed toward a reduced critical region of 0.9 Mb. This region contained the GRXCR1 gene, and the orthologous mouse gene was described to be mutated in the pirouette (pi) mutant with resulting hearing loss and circling behavior. Sequence analysis of the GRXCR1 gene in hearing-impaired family members revealed splice-site mutations in two Dutch families and a missense and nonsense mutation, respectively, in two Pakistani families. The splice-site mutations are predicted to cause frameshifts and premature stop codons. In family W98-053, this could be confirmed by cDNA analysis. GRXCR1 is predicted to contain a GRX-like domain. GRX domains are involved in reversible S-glutathionylation of proteins and thereby in the modulation of activity and/or localization of these proteins. The missense mutation is located in this domain, whereas the nonsense and splice-site mutations may result in complete or partial absence of the GRX-like domain or of the complete protein. Hearing loss in patients with GRXCR1 mutations is congenital and is moderate to profound. Progression of the hearing loss was observed in family W98-053. Vestibular dysfunction was observed in some but not all affected individuals. Quantitative analysis of GRXCR1 transcripts in fetal and adult human tissues revealed a preferential expression of the gene in fetal cochlea, which may explain the nonsyndromic nature of the hearing impairment.

Joubert syndrome and related disorders (JSRD) are primarily autosomal-recessive conditions characterized by hypotonia, ataxia, abnormal eye movements, and intellectual disability with a distinctive mid-hindbrain malformation. Variable features include retinal dystrophy, cystic kidney disease, and liver fibrosis. JSRD are included in the rapidly expanding group of disorders called ciliopathies, because all six gene products implicated in JSRD (NPHP1, AHI1, CEP290, RPGRIP1L, TMEM67, and ARL13B) function in the primary cilium/basal body organelle. By using homozygosity mapping in consanguineous families, we identify loss-of-function mutations in CC2D2A in JSRD patients with and without retinal, kidney, and liver disease. CC2D2A is expressed in all fetal and adult tissues tested. In ciliated cells, we observe localization of recombinant CC2D2A at the basal body and colocalization with CEP290, whose cognate gene is mutated in multiple hereditary ciliopathies. In addition, the proteins can physically interact in vitro, as shown by yeast two-hybrid and GST pull-down experiments. A nonsense mutation in the zebrafish CC2D2A ortholog (sentinel) results in pronephric cysts, a hallmark of ciliary dysfunction analogous to human cystic kidney disease. Knockdown of cep290 function in sentinel fish results in a synergistic pronephric cyst phenotype, revealing a genetic interaction between CC2D2A and CEP290 and implicating CC2D2A in cilium/basal body function. These observations extend the genetic spectrum of JSRD and provide a model system for studying extragenic modifiers in JSRD and other ciliopathies.

Tissue-specific alternative splicing is an important mechanism for providing spatiotemporal protein diversity. Here we show that an in-frame splice mutation in BBS8, one of the genes involved in pleiotropic Bardet-Biedl syndrome (BBS), is sufficient to cause nonsyndromic retinitis pigmentosa (RP). A genome-wide scan of a consanguineous RP pedigree mapped the trait to a 5.6 Mb region; subsequent systematic sequencing of candidate transcripts identified a homozygous splice-site mutation in a previously unknown BBS8 exon. The allele segregated with the disorder, was absent from controls, was completely invariant across evolution, and was predicted to lead to the elimination of a 10 amino acid sequence from the protein. Subsequent studies showed the exon to be expressed exclusively in the retina and enriched significantly in the photoreceptor layer. Importantly, we found this exon to represent the major BBS8 mRNA species in the mammalian photoreceptor, suggesting that the encoded 10 amino acids play a pivotal role in the function of BBS8 in this organ. Understanding the role of this additional sequence might therefore inform the mechanism of retinal degeneration in patients with syndromic BBS or other related ciliopathies.

We have identified a large multigenerational Austrian family displaying a novel form of X-linked recessive myopathy. Affected individuals develop an adult-onset scapulo-axio-peroneal myopathy with bent-spine syndrome characterized by specific atrophy of postural muscles along with pseudoathleticism or hypertrophy and cardiac involvement. Known X-linked myopathies were excluded by simple-tandem-repeat polymorphism (STRP) and single-nucleotide polymorphism (SNP) analysis, direct gene sequencing, and immunohistochemical analysis. STRP analysis revealed significant linkage at Xq25-q27.1. Haplotype analysis based on SNP microarray data from selected family members confirmed this linkage region on the distal arm of the X chromosome, thereby narrowing down the critical interval to 12 Mb. Sequencing of functional candidate genes led to the identification of a missense mutation within the four and a half LIM domain 1 gene (FHL1), which putatively disrupts the fourth LIM domain of the protein. Mutation screening of FHL1 in a myopathy family from the UK exhibiting an almost identical phenotype revealed a 3 bp insertion mutation within the second LIM domain. FHL1 on Xq26.3 is highly expressed in skeletal and cardiac muscles. Western-blot analysis of muscle biopsies showed a marked decrease in protein expression of FHL1 in patients, in concordance with the genetic data. In summary, we have to our knowledge characterized a new disorder, X-linked myopathy with postural muscle atrophy (XMPMA), and identified FHL1 as the causative gene. This is the first FHL protein to be identified in conjunction with a human genetic disorder and further supports the role of FHL proteins in the development and maintenance of muscle tissue. Mutation screening of FHL1 should be considered for patients with uncharacterized myopathies and cardiomyopathies.