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.

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.

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.

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.

The inner ear has fluid-filled compartments of different ionic compositions, including the endolymphatic and perilymphatic spaces of the organ of Corti; the separation from one another by epithelial barriers is required for normal hearing. TRIC encodes tricellulin, a recently discovered tight-junction (TJ) protein that contributes to the structure and function of tricellular contacts of neighboring cells in many epithelial tissues. We show that, in humans, four different recessive mutations of TRIC cause nonsyndromic deafness (DFNB49), a surprisingly limited phenotype, given the widespread tissue distribution of tricellulin in epithelial cells. In the inner ear, tricellulin is concentrated at the tricellular TJs in cochlear and vestibular epithelia, including the structurally complex and extensive junctions between supporting and hair cells. We also demonstrate that there are multiple alternatively spliced isoforms of TRIC in various tissues and that mutations of TRIC associated with hearing loss remove all or most of a conserved region in the cytosolic domain that binds to the cytosolic scaffolding protein ZO-1. A wild-type isoform of tricellulin, which lacks this conserved region, is unaffected by the mutant alleles and is hypothesized to be sufficient for structural and functional integrity of epithelial barriers outside the inner ear.

Primary congenital glaucoma (PCG) is a clinically and genetically heterogeneous disease. The present study was undertaken to find the genetic causes of PCG segregating in 36 large consanguineous Pakistani families.

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.

Intellectual disability (ID) is a genetically and clinically heterogeneous disorder, characterized by limited cognitive abilities and impaired adaptive behaviors. In recent years, exome sequencing (ES) has been instrumental in deciphering the genetic etiology of ID. Here, through ES of a large cohort of individuals with ID, we identified two bi-allelic frameshift variants in METTL5, c.344_345delGA (p.Arg115Asnfs∗19) and c.571_572delAA (p.Lys191Valfs∗10), in families of Pakistani and Yemenite origin. Both of these variants were segregating with moderate to severe ID, microcephaly, and various facial dysmorphisms, in an autosomal-recessive fashion. METTL5 is a member of the methyltransferase-like protein family, which encompasses proteins with a seven-beta-strand methyltransferase domain. We found METTL5 expression in various substructures of rodent and human brains and METTL5 protein to be enriched in the nucleus and synapses of the hippocampal neurons. Functional studies of these truncating variants in transiently transfected orthologous cells and cultured hippocampal rat neurons revealed no effect on the localization of METTL5 but alter its level of expression. Our in silico analysis and 3D modeling simulation predict disruption of METTL5 function by both variants. Finally, mettl5 knockdown in zebrafish resulted in microcephaly, recapitulating the human phenotype. This study provides evidence that biallelic variants in METTL5 cause ID and microcephaly in humans and highlights the essential role of METTL5 in brain development and neuronal function.

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.

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.

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.

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.

Non-secretor status due to homozygosity for the common FUT2 variant c.461G>A (p.Trp154∗) is associated with either risk for autoimmune diseases or protection against viral diarrhea and HIV. We determined the role of FUT2 in otitis media susceptibility by obtaining DNA samples from 609 multi-ethnic families and simplex case subjects with otitis media. Exome and Sanger sequencing, linkage analysis, and Fisher exact and transmission disequilibrium tests (TDT) were performed. The common FUT2 c.604C>T (p.Arg202∗) variant co-segregates with otitis media in a Filipino pedigree (LOD = 4.0). Additionally, a rare variant, c.412C>T (p.Arg138Cys), is associated with recurrent/chronic otitis media in European-American children (p = 1.2 × 10-5) and US trios (TDT p = 0.01). The c.461G>A (p.Trp154∗) variant was also over-transmitted in US trios (TDT p = 0.01) and was associated with shifts in middle ear microbiota composition (PERMANOVA p < 10-7) and increased biodiversity. When all missense and nonsense variants identified in multi-ethnic US trios with CADD > 20 were combined, FUT2 variants were over-transmitted in trios (TDT p = 0.001). Fut2 is transiently upregulated in mouse middle ear after inoculation with non-typeable Haemophilus influenzae. Four FUT2 variants-namely p.Ala104Val, p.Arg138Cys, p.Trp154∗, and p.Arg202∗-reduced A antigen in mutant-transfected COS-7 cells, while the nonsense variants also reduced FUT2 protein levels. Common and rare FUT2 variants confer susceptibility to otitis media, likely by modifying the middle ear microbiome through regulation of A antigen levels in epithelial cells. Our families demonstrate marked intra-familial genetic heterogeneity, suggesting that multiple combinations of common and rare variants plus environmental factors influence the individual otitis media phenotype as a complex trait.

Juvenile open angle glaucoma (JOAG), which is an uncommon form of primary open angle glaucoma, is a clinically and genetically heterogeneous disorder. We report on a family with a recessively inherited form of JOAG. The proband has a superior and an inferior never fiber layer thinning in both the eyes and the nasal visual field (VF) defects in the left eye, which are clinical findings consistent with glaucomatous optic neuropathy. Whole exome sequencing revealed two novel compound heterozygous variants [c.2966C>G, p.(Pro989Arg); c.5235T>G, p.(Asn1745Lys)] in latent transforming growth factor-beta-binding protein 2 (LTBP2) segregating with the phenotype. Both these variants are predicted to replace evolutionary conserved amino acids, have a pathogenic effect on the encode protein, and have very low frequencies in the control databases. Mutations in LTBP2 are known to cause the Weill-Marchesani syndrome and a Weill-Marchesani-like syndrome, which include glaucoma in their clinical presentation. However, to our knowledge, this is the first published case of a JOAG subject associated with recessively inherited variants of LTPB2 and, thus, expands the repertoire of the known genetic causes of JOAG and the phenotypic spectrum of LTBP2 alleles.

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.

Epidermolysis bullosa (EB) includes a group of rare gesnodermatoses that result in blistering and erosions of the skin and mucous membranes. Genetically, pathogenic variants in around 20 genes are known to alter the structural and functional integrity of intraepidermal adhesion and dermo-epidermal anchorage, leading to four different types of EB. Here we report the underlying genetic causes of EB phenotypes segregating in seven large consanguineous families, recruited from different regions of Pakistan. Whole exome sequencing, followed by segregation analysis of candidate variants through Sanger sequencing, identified eight pathogenic variants, including three novel (ITGB4: c.1285G>T, and c.3373G>A; PLEC: c.1828A>G) and five previously reported variants (COL7A1: c.6209G>A, and c.1573C>T; FERMT1: c.676insC; LAMA3: c.151insG; LAMB3: c.1705C>T). All identified variants were either absent or had very low frequencies in the control databases. Our in-silico analyses and 3-dimensional (3D) molecular modeling support the deleterious impact of these variants on the encoded proteins. Intriguingly, we report the first case of a recessively inherited form of rare EBS-Ogna associated with a homozygous variant in the PLEC gene. Our study highlights the clinical and genetic diversity of EB in the Pakistani population and expands the mutation spectrum of EB; it could also be useful for prenatal diagnosis and genetic counseling of the affected families.

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.

The genetic etiologies of more than half of rare diseases remain unknown. Standardized genome sequencing and phenotyping of large patient cohorts provide an opportunity for discovering the unknown etiologies, but this depends on efficient and powerful analytical methods. We built a compact database, the 'Rareservoir', containing the rare variant genotypes and phenotypes of 77,539 participants sequenced by the 100,000 Genomes Project. We then used the Bayesian genetic association method BeviMed to infer associations between genes and each of 269 rare disease classes assigned by clinicians to the participants. We identified 241 known and 19 previously unidentified associations. We validated associations with ERG, PMEPA1 and GPR156 by searching for pedigrees in other cohorts and using bioinformatic and experimental approaches. We provide evidence that (1) loss-of-function variants in the Erythroblast Transformation Specific (ETS)-family transcription factor encoding gene ERG lead to primary lymphoedema, (2) truncating variants in the last exon of transforming growth factor-β regulator PMEPA1 result in Loeys-Dietz syndrome and (3) loss-of-function variants in GPR156 give rise to recessive congenital hearing impairment. The Rareservoir provides a lightweight, flexible and portable system for synthesizing the genetic and phenotypic data required to study rare disease cohorts with tens of thousands of participants.

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.