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.

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.

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.

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.

Lapachone is a natural naphthoquinone-derived compound found in Tabebuia avellanedae. It is well-known for its analgesic, anti-inflammatory, anti-microbial, diuretic, and anti-cancerous effects. However, the wound-healing effects of this compound are not known yet. The aim of this study was to investigate the wound healing activity of naphthoquinones (α-lapachone and β-lapachone) from Handroanthus impetiginosus.

Cargo transport along the cytoplasmic microtubular network is essential for neuronal function, and cytoplasmic dynein-1 is an established molecular motor that is critical for neurogenesis and homeostasis. We performed whole-exome sequencing, homozygosity mapping, and chromosomal microarray studies in five individuals from three independent pedigrees and identified likely-pathogenic variants in DYNC1I2 (Dynein Cytoplasmic 1 Intermediate Chain 2), encoding a component of the cytoplasmic dynein 1 complex. In a consanguineous Pakistani family with three affected individuals presenting with microcephaly, severe intellectual disability, simplification of cerebral gyration, corpus callosum hypoplasia, and dysmorphic facial features, we identified a homozygous splice donor site variant (GenBank: NM_001378.2:c.607+1G>A). We report two additional individuals who have similar neurodevelopmental deficits and craniofacial features and harbor deleterious variants; one individual bears a c.740A>G (p.Tyr247Cys) change in trans with a 374 kb deletion encompassing DYNC1I2, and an unrelated individual harbors the compound-heterozygous variants c.868C>T (p.Gln290∗) and c.740A>G (p.Tyr247Cys). Zebrafish larvae subjected to CRISPR-Cas9 gene disruption or transient suppression of dync1i2a displayed significantly altered craniofacial patterning with concomitant reduction in head size. We monitored cell death and cell cycle progression in dync1i2a zebrafish models and observed significantly increased apoptosis, likely due to prolonged mitosis caused by abnormal spindle morphology, and this finding offers initial insights into the cellular basis of microcephaly. Additionally, complementation studies in zebrafish demonstrate that p.Tyr247Cys attenuates gene function, consistent with protein structural analysis. Our genetic and functional data indicate that DYNC1I2 dysfunction probably causes an autosomal-recessive microcephaly syndrome and highlight further the critical roles of the dynein-1 complex in neurodevelopment.

Germline pathogenic variants in chromatin-modifying enzymes are a common cause of pediatric developmental disorders. These enzymes catalyze reactions that regulate epigenetic inheritance via histone post-translational modifications and DNA methylation. Cytosine methylation (5-methylcytosine [5mC]) of DNA is the quintessential epigenetic mark, yet no human Mendelian disorder of DNA demethylation has yet been delineated. Here, we describe in detail a Mendelian disorder caused by the disruption of DNA demethylation. TET3 is a methylcytosine dioxygenase that initiates DNA demethylation during early zygote formation, embryogenesis, and neuronal differentiation and is intolerant to haploinsufficiency in mice and humans. We identify and characterize 11 cases of human TET3 deficiency in eight families with the common phenotypic features of intellectual disability and/or global developmental delay; hypotonia; autistic traits; movement disorders; growth abnormalities; and facial dysmorphism. Mono-allelic frameshift and nonsense variants in TET3 occur throughout the coding region. Mono-allelic and bi-allelic missense variants localize to conserved residues; all but one such variant occur within the catalytic domain, and most display hypomorphic function in an assay of catalytic activity. TET3 deficiency and other Mendelian disorders of the epigenetic machinery show substantial phenotypic overlap, including features of intellectual disability and abnormal growth, underscoring shared disease mechanisms.

TET3 encodes an essential dioxygenase involved in epigenetic regulation through DNA demethylation. TET3 deficiency, or Beck-Fahrner syndrome (BEFAHRS; MIM: 618798), is a recently described neurodevelopmental disorder of the DNA demethylation machinery with a nonspecific phenotype resembling other chromatin-modifying disorders, but inconsistent variant types and inheritance patterns pose diagnostic challenges. Given TET3's direct role in regulating 5-methylcytosine and recent identification of syndrome-specific DNA methylation profiles, we analyzed genome-wide DNA methylation in whole blood of TET3-deficient individuals and identified an episignature that distinguishes affected and unaffected individuals and those with mono-allelic and bi-allelic pathogenic variants. Validation and testing of the episignature correctly categorized known TET3 variants and determined pathogenicity of variants of uncertain significance. Clinical utility was demonstrated when the episignature alone identified an affected individual from over 1000 undiagnosed cases and was confirmed upon distinguishing TET3-deficient individuals from those with 46 other disorders. The TET3-deficient signature - and the signature resulting from activating mutations in DNMT1 which normally opposes TET3 - are characterized by hypermethylation, which for BEFAHRS involves CpG sites that may be biologically relevant. This work expands the role of epi-phenotyping in molecular diagnosis and reveals genome-wide DNA methylation profiling as a quantitative, functional readout for characterization of this new biochemical category of disease.

Hearing impairment (HI) is a common disorder of sensorineural function with a highly heterogeneous genetic background. Although substantial progress has been made in the understanding of the genetic etiology of hereditary HI, many genes implicated in HI remain undiscovered. Via exome and Sanger sequencing of DNA samples obtained from consanguineous Pakistani families that segregate profound prelingual sensorineural HI, we identified rare homozygous missense variants in four genes (ADAMTS1, MPDZ, MVD, and SEZ6) that are likely the underlying cause of HI. Linkage analysis provided statistical evidence that these variants are associated with autosomal recessive nonsyndromic HI. In silico analysis of the mutant proteins encoded by these genes predicted structural, conformational or interaction changes. RNAseq data analysis revealed expression of these genes in the sensory epithelium of the mouse inner ear during embryonic, postnatal, and adult stages. Immunohistochemistry of the mouse cochlear tissue, further confirmed the expression of ADAMTS1, SEZ6, and MPDZ in the neurosensory hair cells of the organ of Corti, while MVD expression was more prominent in the spiral ganglion cells. Overall, supported by in silico mutant protein analysis, animal models, linkage analysis, and spatiotemporal expression profiling in the mouse inner ear, we propose four new candidate genes for HI and expand our understanding of the etiology of HI.

We have used genome-wide genotyping to identify an overlapping homozygosity-by-descent locus on chromosome 9q34.3 (MRT15) in four consanguineous families affected by nonsyndromic autosomal-recessive intellectual disability (NS-ARID) and one in which the patients show additional clinical features. Four of the families are from Pakistan, and one is from Iran. Using a combination of next-generation sequencing and Sanger sequencing, we have identified mutations in the gene MAN1B1, encoding a mannosyl oligosaccharide, alpha 1,2-mannosidase. In one Pakistani family, MR43, a homozygous nonsense mutation (RefSeq number NM_016219.3: c.1418G>A [p.Trp473*]), segregated with intellectual disability and additional dysmorphic features. We also identified the missense mutation c. 1189G>A (p.Glu397Lys; RefSeq number NM_016219.3), which segregates with NS-ARID in three families who come from the same village and probably have shared inheritance. In the Iranian family, the missense mutation c.1000C>T (p.Arg334Cys; RefSeq number NM_016219.3) also segregates with NS-ARID. Both missense mutations are at amino acid residues that are conserved across the animal kingdom, and they either reduce k(cat) by ∼1300-fold or disrupt stable protein expression in mammalian cells. MAN1B1 is one of the few NS-ARID genes with an elevated mutation frequency in patients with NS-ARID from different populations.

To investigate the molecular basis of retinitis pigmentosa in two consanguineous families of Pakistani origin with multiple affected members.

The DFNB74 locus for autosomal-recessive, nonsyndromic deafness segregating in three families was previously mapped to a 5.36 Mb interval on chromosome 12q14.2-q15. Subsequently, we ascertained five additional consanguineous families in which deafness segregated with markers at this locus and refined the critical interval to 2.31 Mb. We then sequenced the protein-coding exons of 18 genes in this interval. The affected individuals of six apparently unrelated families were homozygous for the same transversion (c.265T>G) in MSRB3, which encodes a zinc-containing methionine sulfoxide reductase B3. c.265T>G results in a substitution of glycine for cysteine (p.Cys89Gly), and this substitution cosegregates with deafness in the six DFNB74 families. This cysteine residue of MSRB3 is conserved in orthologs from yeast to humans and is involved in binding structural zinc. In vitro, p.Cys89Gly abolished zinc binding and MSRB3 enzymatic activity, indicating that p.Cys89Gly is a loss-of-function allele. The affected individuals in two other families were homozygous for a transition mutation (c.55T>C), which results in a nonsense mutation (p.Arg19X) in alternatively spliced exon 3, encoding a mitochondrial localization signal. This finding suggests that DFNB74 deafness is due to a mitochondrial dysfunction. In a cohort of 1,040 individuals (aged 53-67 years) of European ancestry, we found no association between 17 tagSNPs for MSRB3 and age-related hearing loss. Mouse Msrb3 is expressed widely. In the inner ear, it is found in the sensory epithelium of the organ of Corti and vestibular end organs as well as in cells of the spiral ganglion. Taken together, MSRB3-catalyzed reduction of methionine sulfoxides to methionine is essential for hearing.

The RDX gene anchors cytoskeletal actin of stereocilia to hair cell transmembrane and is responsible for autosomal recessive nonsyndromic hearing impairment (ARNSHI) due to DFNB24. A genome scan was performed using DNA samples from a consanguineous Pakistani family with ARNSHI. A significant maximum two-point LOD score of 4.5 (θ = 0) and multipoint LOD score of 5.8 were achieved at marker D11S1998 (chr11 : 117.20 Mb). The region of homozygosity is bounded by markers D11S2000 (105.06 Mb) and D11S4464 (123.13 Mb) and contains the NSHI genes TECTA and RDX. Although no potentially causal variants were identified in the TECTA gene, within the RDX gene a novel deletion c.1076_1079delTTAA (p.Ile359Lysfs∗6) was identified. The RDX deletion segregates with ARNSHI within the family and was not observed in 500 control chromosomes. It is predicted to cause premature truncation of radixin at the α-helical domain and to result in nonfunctional transcripts within the cochlea. RDX isoforms which encode the coiled-coil region of the α-helical domain are deemed necessary for proper function of hair cell stereocilia.

Mutations in the estrogen-related receptor beta (ESRRB) gene is the underlying cause of autosomal recessive nonsyndromic hearing impairment (ARNSHI) due to the DFNB35 locus which maps to 14q24.3. A genome scan of a large consanguineous Pakistani pedigree with ARNSHI established linkage with a maximum multipoint LOD score of 4.2 to the 14q24 region and the region of homozygosity contained the ESRRB gene. Sequencing of the ESRRB gene using DNA samples from hearing-impaired family members uncovered a novel three-nucleotide deletion c.1018_1020delGAG (p.Glu340del). The deletion segregates with hearing impairment in the pedigree and was not observed in 500 control chromosomes. The deletion of glutamic acid residue occurs in the ligand-binding domain of ESRRB protein. It is expected that the deletion affects the ligand-binding activity of the domain in ESRRB, which leads to the ARNSHI.

A novel ARNSHI (autosomal recessive non-syndromic hearing impairment) locus, DFNB71, was localized to 8p22-21.3. To map the locus, a whole-genome scan was carried out using DNA samples from a consanguineous seven-generational Pakistani family with profound prelingual ARNSHI. A maximum multipoint LOD score of 4.2 occurred at marker D8S261. The DFNB71 region of homozygosity and 3-unit support interval is flanked by markers D8S1130 and D8S1786. This region has a genetic distance of 19.1 cM and contains 10.6 Mb of sequence.

Cohen syndrome (CS) is a rare genetic disorder caused by variants of the VPS13B gene. CS patients are affected with a severe form of retinal dystrophy, and in several cases cataracts also develop. The purpose of this study was to investigate the mechanisms and risk factors for cataract in CS, as well as to report on cataract surgeries in CS patients.

The recent availability of whole-exome sequencing has opened new possibilities for the evaluation of individuals with genetically undiagnosed intellectual disability.

We report two consanguineous families with probands that exhibit intellectual disability, developmental delay, short stature, aphasia, and hypotonia in which homozygous non-synonymous variants were identified in IQSEC1 (GenBank: NM_001134382.3). In a Pakistani family, the IQSEC1 segregating variant is c.1028C>T (p.Thr343Met), while in a Saudi Arabian family the variant is c.962G>A (p.Arg321Gln). IQSEC1-3 encode guanine nucleotide exchange factors for the small GTPase ARF6 and their loss affects a variety of actin-dependent cellular processes, including AMPA receptor trafficking at synapses. The ortholog of IQSECs in the fly is schizo and its loss affects growth cone guidance at the midline in the CNS, also an actin-dependent process. Overexpression of the reference IQSEC1 cDNA in wild-type flies is lethal, but overexpression of the two variant IQSEC1 cDNAs did not affect viability. Loss of schizo caused embryonic lethality that could be rescued to 2nd instar larvae by moderate expression of the human reference cDNA. However, the p.Arg321Gln and p.Thr343Met variants failed to rescue embryonic lethality. These data indicate that the variants behave as loss-of-function mutations. We also show that schizo in photoreceptors is required for phototransduction. Finally, mice with a conditional Iqsec1 deletion in cortical neurons exhibited an increased density of dendritic spines with an immature morphology. The phenotypic similarity of the affecteds and the functional experiments in flies and mice indicate that IQSEC1 variants are the cause of a recessive disease with intellectual disability, developmental delay, and short stature, and that axonal guidance and dendritic projection defects as well as dendritic spine dysgenesis may underlie disease pathogenesis.

LHFPL5, the gene for DFNB67, underlies autosomal recessive nonsyndromic hearing impairment. We identified seven Pakistani families that mapped to 6p21.31, which includes the LHFPL5 gene. Sanger sequencing of LHFPL5 using DNA samples from hearing impaired and unaffected members of these seven families identified four variants. Among the identified variants, two were novel: one missense c.452 G > T (p.Gly151Val) and one splice site variant (c.*16 + 1 G > A) were each identified in two families. Two known variants: c.250delC (p.Leu84*) and c.380 A > G (p.Tyr127Cys) were also observed in two families and a single family, respectively. Nucleotides c.452G and c.*16 + 1G and amino-acid residue p.Gly151 are under strong evolutionary conservation. In silico bioinformatics analyses predicted these variants to be damaging. The splice site variant (c.*16 + 1 G > A) is predicted to affect pre-mRNA splicing and a loss of the 5' donor splice site in the 3'-untranslated region (3'-UTR). Further analysis supports the activation of a cryptic splice site approximately 357-bp downstream, leading to an extended 3'-UTR with additional regulatory motifs. In conclusion, we identified two novel variants in LHFPL5, including a unique 3'-UTR splice site variant that is predicted to impact pre-mRNA splicing and regulation through an extended 3'-UTR.