Genome sequencing reveals novel causative structural and single nucleotide variants in Pakistani families with congenital hypogonadotropic hypogonadism.

BACKGROUND/OBJECTIVES: This study aims to elucidate the genetic causes of congenital hypogonadotropic hypogonadism (CHH), a rare genetic disorder resulting in GnRH deficiency, in six families from Pakistan. METHODS: Eighteen DNA samples from six families underwent genome sequencing followed by standard evaluation for pathogenic single nucleotide variants (SNVs) and small indels. All families were subsequently analyzed for pathogenic copy number variants (CNVs) using CoverageMaster. RESULTS: Novel pathogenic homozygous SNVs in known CHH genes were identified in four families: two families with variants in GNRHR, and two others harboring KISS1R variants. Subsequent investigation of CNVs in the remaining two families identified novel unique large deletions in ANOS1. CONCLUSION: A combined, systematic analysis of single nucleotide and CNVs helps to improve the diagnostic yield for variants in patients with CHH.

Monoallelic and bi-allelic variants in NCDN cause neurodevelopmental delay, intellectual disability, and epilepsy.

Neurochondrin (NCDN) is a cytoplasmatic neural protein of importance for neural growth, glutamate receptor (mGluR) signaling, and synaptic plasticity. Conditional loss of Ncdn in mice neural tissue causes depressive-like behaviors, impaired spatial learning, and epileptic seizures. We report on NCDN missense variants in six affected individuals with variable degrees of developmental delay, intellectual disability (ID), and seizures. Three siblings were found homozygous for a NCDN missense variant, whereas another three unrelated individuals carried different de novo missense variants in NCDN. We assayed the missense variants for their capability to rescue impaired neurite formation in human neuroblastoma (SH-SY5Y) cells depleted of NCDN. Overexpression of wild-type NCDN rescued the neurite-phenotype in contrast to expression of NCDN containing the variants of affected individuals. Two missense variants, associated with severe neurodevelopmental features and epilepsy, were unable to restore mGluR5-induced ERK phosphorylation. Electrophysiological analysis of SH-SY5Y cells depleted of NCDN exhibited altered membrane potential and impaired action potentials at repolarization, suggesting NCDN to be required for normal biophysical properties. Using available transcriptome data from human fetal cortex, we show that NCDN is highly expressed in maturing excitatory neurons. In combination, our data provide evidence that bi-allelic and de novo variants in NCDN cause a clinically variable form of neurodevelopmental delay and epilepsy, highlighting a critical role for NCDN in human brain development.

Exome sequencing in four families with neurodevelopmental disorders: genotype-phenotype correlation and identification of novel disease-causing variants in VPS13B and RELN.

Neurodevelopmental disorders (NDDs) are a clinically and genetically heterogeneous group of early-onset pediatric disorders that affect the structure and/or function of the central or peripheral nervous system. Achieving a precise molecular diagnosis for NDDs may be challenging due to the diverse genetic underpinnings and clinical variability. In the current study, we investigated the underlying genetic cause(s) of NDDs in four unrelated Pakistani families. Using exome sequencing (ES) as a diagnostic approach, we identified disease-causing variants in established NDD-associated genes in all families, including one hitherto unreported variant in RELN and three recurrent variants in VPS13B, DEGS1, and SPG11. Overall, our study highlights the potential of ES as a tool for clinical diagnosis.

A homozygous founder variant in PDE2A causes paroxysmal dyskinesia with intellectual disability.

Intellectual developmental disorder with paroxysmal dyskinesia or seizures (IDDPADS, OMIM#619150) is an ultra-rare childhood-onset autosomal recessive movement disorder manifesting paroxysmal dyskinesia, global developmental delay, impaired cognition, progressive psychomotor deterioration and/or drug-refractory seizures. We investigated three consanguineous Pakistani families with six affected individuals presenting overlapping phenotypes partially consistent with the reported characteristics of IDDPADS. Whole exome sequencing identified a novel missense variant in Phosphodiesterase 2A (PDE2A): NM_002599.4: c.1514T > C p.(Phe505Ser) that segregated with the disease status of individuals in these families. Retrospectively, we performed haplotype analysis that revealed a 3.16 Mb shared haplotype at 11q13.4 among three families suggesting a founder effect in this region. Moreover, we also observed abnormal mitochondrial morphology in patient fibroblasts compared to controls. Belonging to diverse age groups (13 years-60 years), patients presented paroxysmal dyskinesia, developmental delay, cognitive abnormalities, speech impairment, and drug-refractory seizures with variable onset of disease (as early as 3 months of age to 7 years). Together with the previous reports, we observed that intellectual disability, progressive psychomotor deterioration, and drug-refractory seizures are consistent outcomes of the disease. However, permanent choreodystonia showed variability. We also noticed that the later onset of paroxysmal dyskinesia manifests severe attacks in terms of duration. Being the first report from Pakistan, we add to the clinical and mutation spectrum of PDE2A-related recessive disease raising the total number of patients from six to 12 and variants from five to six. Together, with our findings, the role of PDE2A is strengthened in critical physio-neurological processes.

Altered paracellular cation permeability due to a rare CLDN10B variant causes anhidrosis and kidney damage.

Claudins constitute the major component of tight junctions and regulate paracellular permeability of epithelia. Claudin-10 occurs in two major isoforms that form paracellular channels with ion selectivity. We report on two families segregating an autosomal recessive disorder characterized by generalized anhidrosis, severe heat intolerance and mild kidney failure. All affected individuals carry a rare homozygous missense mutation c.144C>G, p.(N48K) specific for the claudin-10b isoform. Immunostaining of sweat glands from patients suggested that the disease is associated with reduced levels of claudin-10b in the plasma membranes and in canaliculi of the secretory portion. Expression of claudin-10b N48K in a 3D cell model of sweat secretion indicated perturbed paracellular Na+ transport. Analysis of paracellular permeability revealed that claudin-10b N48K maintained cation over anion selectivity but with a reduced general ion conductance. Furthermore, freeze fracture electron microscopy showed that claudin-10b N48K was associated with impaired tight junction strand formation and altered cis-oligomer formation. These data suggest that claudin-10b N48K causes anhidrosis and our findings are consistent with a combined effect from perturbed TJ function and increased degradation of claudin-10b N48K in the sweat glands. Furthermore, affected individuals present with Mg2+ retention, secondary hyperparathyroidism and mild kidney failure that suggest a disturbed reabsorption of cations in the kidneys. These renal-derived features recapitulate several phenotypic aspects detected in mice with kidney specific loss of both claudin-10 isoforms. Our study adds to the spectrum of phenotypes caused by tight junction proteins and demonstrates a pivotal role for claudin-10b in maintaining paracellular Na+ permeability for sweat production and kidney function.

Primary microcephaly, primordial dwarfism, and brachydactyly in adult cases with biallelic skipping of RTTN exon 42.

Biallelic and pathogenic variants in the RTTN gene, encoding the centrosomal protein Rotatin, are associated with variable degrees of neurodevelopmental abnormalities, microcephaly, and extracranial malformations. To date, no reported case has reached their third decade. Herein, we report on a consanguineous family with three adult members, age 43, 57, and 60 years respectively, with primary microcephaly, developmental delay, primordial dwarfism, and brachydactyly segregating a homozygous splice site variant NM_173630.3:c.5648-5T>A in RTTN. The variant RTTN allele results in a nonhypomorphic skipping of exon 42 and a frameshift [(NP_775901.3:p.Ala1883Glyfs*6)]. Brain MRI of one affected individual showed markedly reduced volume of cerebral lobes and enlarged sulci but without signs of neural migration defects. Our assessment of three adult cases with a biallelic RTTN variant shows that a predicted shortened Rotatin, lacking the C-terminal end, are associated with stationary clinical features into the seventh decade. Furthermore, our report adds brachydactyly to the phenotypic spectrum in this pleiotropic entity.

Mutations of KIF14 cause primary microcephaly by impairing cytokinesis.

OBJECTIVE: Autosomal recessive primary microcephaly (MCPH) is a rare condition characterized by a reduced cerebral cortex accompanied with intellectual disability. Mutations in 17 genes have been shown to cause this phenotype. Recently, mutations in CIT, encoding CRIK (citron rho-interacting kinase)-a component of the central spindle matrix-were added. We aimed at identifying novel MCPH-associated genes and exploring their functional role in pathogenesis. METHODS: Linkage analysis and whole exome sequencing were performed in consanguineous and nonconsanguineous MCPH families to identify disease-causing variants. Functional consequences were investigated by RNA studies and on the cellular level using immunofluorescence and microscopy. RESULTS: We identified homozygous mutations in KIF14 (NM_014875.2;c.263T>A;pLeu88*, c.2480_2482delTTG; p.Val827del, and c.4071G>A;p.Gln1357=) as the likely cause in 3 MCPH families. Furthermore, in a patient presenting with a severe form of primary microcephaly and short stature, we identified compound heterozygous missense mutations in KIF14 (NM_014875.2;c.2545C>G;p.His849Asp and c.3662G>T;p.Gly1221Val). Three of the 5 identified mutations impaired splicing, and 2 resulted in a truncated protein. Intriguingly, Kif14 knockout mice also showed primary microcephaly. Human kinesin-like protein KIF14, a microtubule motor protein, localizes at the midbody to finalize cytokinesis by interacting with CRIK. We found impaired localization of both KIF14 and CRIK at the midbody in patient-derived fibroblasts. Furthermore, we observed a large number of binucleated and apoptotic cells-signs of failed cytokinesis that we also observed in experimentally KIF14-depleted cells. INTERPRETATION: Our data corroborate the role of an impaired cytokinesis in the etiology of primary and syndromic microcephaly, as has been proposed by recent findings on CIT mutations. Ann Neurol 2017;82:562-577.

Homozygous mutation in the NPHP3 gene causing foetal nephronophthisis.

We present a case of a foetal sonographic finding of hyper-echogenic kidneys, which led to a strategic series of genetic tests and identified a homozygous mutation (c.424C > T, p. R142*) in the NPHP3 gene. Our study provides a rare presentation of NPHP3-related ciliopathy and adds to the mutation spectrum of the gene, being the first one from Pakistani population. With a thorough literature review, it also advocates for molecular assessment of ciliopathies to improve risk estimate for future pregnancies, and identify predisposed asymptomatic carriers.

CDK6 associates with the centrosome during mitosis and is mutated in a large Pakistani family with primary microcephaly.

Autosomal recessive primary microcephaly (MCPH) is characterized by reduced head circumference, reduction in the size of the cerebral cortex with otherwise grossly normal brain structure and variable intellectual disability. MCPH is caused by mutations of 11 different genes which code for proteins implicated in cell division and cell cycle regulation. We studied a consanguineous eight-generation family from Pakistan with ten microcephalic children using homozygosity mapping and found a new MCPH locus at HSA 7q21.11-q21.3. Sanger sequencing of the most relevant candidate genes in this region revealed a homozygous single nucleotide substitution c.589G>A in CDK6, which encodes cyclin-dependent kinase 6. The mutation changes a highly conserved alanine at position 197 into threonine (p.Ala197Thr). Post hoc whole-exome sequencing corroborated this mutation's identification as the causal variant. CDK6 is an important protein for the control of the cell cycle and differentiation of various cell types. We show here for the first time that CDK6 associates with the centrosome during mitosis; however, this was not observed in patient fibroblasts. Moreover, the mutant primary fibroblasts exhibited supernumerary centrosomes, disorganized microtubules and mitotic spindles, an increased centrosome nucleus distance, reduced cell proliferation and impaired cell motility and polarity. Upon ectopic expression of the mutant protein and knockdown of CDK6 through shRNA, we noted similar effects. We propose that the identified CDK6 mutation leads to reduced cell proliferation and impairs the correct functioning of the centrosome in microtubule organization and its positioning near the nucleus which are key determinants during neurogenesis.

Genetic association of vitamin D receptor gene with female infertility.

BACKGROUND: Infertility is defined as failure to achieve a clinical pregnancy after 12 months of unprotected intercourse. It affects 15% of couples globally and 22% of couples within Pakistan. Female infertility can be caused by numerous genetic or environmental factors including hormone imbalances and exposure to chemicals or radiation. The prevalence of vitamin D deficiency among the adult population was reported to be 14-59% with a higher prevalence in Asian countries. Furthermore, the expression of Vitamin D receptor (VDR) can play a vital role in the reproductive organs of females. Hence, the aim of our present study was to check the association of VDR polymorphisms with infertile females. For this purpose, blood samples were collected for genotyping of four known VDR mutations [FokI (rs2228570), TaqI (rs731236), ApaI (rs7975232), and BsmI (rs1544410)] via PCR-based RFLP assay. RESULTS: Genotyping indicated that FokI, TaqI, and ApaI are associated with infertility (p = 0.004*, p = 0.013*, and p = 0.033*, respectively). However, BsmI did not show any significance. Multinomial regression analysis indicated that FokI heterozygous genotypes increase the risk of infertility by 2.5 times (hetero: OR = 2.5, 95%, p = 0.001*) as compared to wild type. Heterozygous genotypes of TaqI and ApaI were found to play a protective role and reduce the risk of infertility by 58 and 52%, respectively [TaqI: OR = 0.42, 95%, p = 0.004*, ApaI: OR = 0.48, 95%, p = 0.01*, respectively] as compared to wild type. Multinomial logistic regression analysis was also performed for allelic data as well. CONCLUSION: Thus, it could be summarized that among the studied polymorphisms of VDR, FokI SNP greatly increased the risk of infertility, while TaqI and ApaI genotypes protect from infertility. However, BsmI does not influence the risk of infertility in Pakistani females.

Generation of a ZEB2 deficient human iPSC line (KICRi002A-4).

The transcription factor ZEB2 is essential for early embryonic development. Using CRISPR/Cas9, we generated a ZEB2 deficient human iPSC cell line (KICRi002A-4), carrying a homozygous 790 bp deletion in ZEB2 that involves part of exon 5, intron 5 and part of exon 6. The deletion leads to markedly reduced levels of a truncated ZEB2 transcript. No ZEB2 protein was detected by immunopreciptation. The iPSC line expressed pluripotency markers and showed a capacity to differentiate into the three germ layers in vitro. Assessment of genomic integrity revealed a normal karyotype without detectable OFF-target editing. The iPSC line KICRi002A-4 thus offers a valuable resource to study the role of ZEB2 for the commitment and differentiation of various human cell lineages.

Genetic Investigation of Consanguineous Pakistani Families Segregating Rare Spinocerebellar Disorders.

Spinocerebellar disorders are a vast group of rare neurogenetic conditions, generally characterized by overlapping clinical symptoms including progressive cerebellar ataxia, spastic paraparesis, cognitive deficiencies, skeletal/muscular and ocular abnormalities. The objective of the present study is to identify the underlying genetic causes of the rare spinocerebellar disorders in the Pakistani population. Herein, nine consanguineous families presenting different spinocerebellar phenotypes have been investigated using whole exome sequencing. Sanger sequencing was performed for segregation analysis in all the available individuals of each family. The molecular analysis of these families identified six novel pathogenic/likely pathogenic variants; ZFYVE26: c.1093del, SACS: c.1201C>T, BICD2: c.2156A>T, ALS2: c.2171-3T>G, ALS2: c.3145T>A, and B4GALNT1: c.334_335dup, and three already reported pathogenic variants; FA2H: c.159_176del, APTX: c.689T>G, and SETX: c.5308_5311del. The clinical features of all patients in each family are concurrent with the already reported cases. Hence, the current study expands the mutation spectrum of rare spinocerebellar disorders and implies the usefulness of next-generation sequencing in combination with clinical investigation for better diagnosis of these overlapping phenotypes.

A Novel Nonsense Variant in GRM1 Causes Autosomal Recessive Spinocerebellar Ataxia 13 in a Consanguineous Pakistani Family.

Background and objectives: Autosomal recessive spinocerebellar ataxia-13 (SCAR13) is an ultra-rare disorder characterized by slowly progressive cerebellar ataxia, cognitive deficiencies, and skeletal and oculomotor abnormalities. The objective of this case report is to expand the clinical and molecular spectrum of SCAR13. Methods: We investigated a consanguineous Pakistani family with four patients partially presenting with clinical features of SCAR13 using whole exome sequencing. Segregation analysis was performed by Sanger sequencing in all the available individuals of the family. Results: Patients presented with quadrupedal gait, delayed developmental milestones, non-progressive peripheral neuropathy, and cognitive impairment. Whole exome sequencing identified a novel pathogenic nonsense homozygous variant, Gly240*, in the gene GRM1 as a cause of SCAR13 that segregates with the recessive disease. Discussion: We report a novel homozygous nonsense variant in the GRM1 gene in four Pakistani patients presenting with clinical features that partially overlap with the already reported phenotype of SCAR13. In addition, the family presented quadrupedal gait and non-progressive symptoms, manifestations which have not been recognized previously. So far, only four variants in GRM1 have been reported, in families of Roma, Iranian, and Tunisian origins. The current study adds to the mutation spectrum of GRM1 and provides a rare presentation of SCAR13, the first from the Pakistani population.

MCPH1: A Novel Case Report and a Review of the Literature.

Microcephaly primary hereditary (MCPH) is a congenital disease characterized by nonsyndromic reduction in brain size due to impaired neurogenesis, often associated with a variable degree of intellectual disability (ID). The genetic etiology of MCPH is heterogeneous and comprises more than 20 loci, nearly all following a recessive inheritance pattern. The first causative gene identified, MCPH1 or Microcephalin, encodes a centrosomal protein that modulates chromosome condensation and cell cycle progression. It is also involved in DNA damage response and telomere maintenance in the nucleus. Despite numerous studies on MCPH1 function, MCPH1-affected individuals are rare and the available clinical reports are not sufficient to define the natural history of the disease. Here, we present a novel patient with congenital microcephaly, ID, language delay, short stature, and other minor features such as strabismus. magnetic resonance imaging revealed ventriculomegaly, simplified gyral pattern in the frontal lobes, and a neuronal migration defect. Genetic testing detected a homozygous deletion of exons 1-8 of MCPH1. We compare the patients' characteristics with a list of features from MCPH1 cases described in the literature, in an effort to provide additional clues for a comprehensive definition of disease presentation and evolution.

The SARS-CoV-2 differential genomic adaptation in response to varying UVindex reveals potential genomic resources for better COVID-19 diagnosis and prevention.

Coronavirus disease 2019 (COVID-19) has been a pandemic disease reported in almost every country and causes life-threatening, severe respiratory symptoms. Recent studies showed that various environmental selection pressures challenge the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infectivity and, in response, the virus engenders new mutations, leading to the emergence of more virulent strains of WHO concern. Advance prediction of the forthcoming virulent SARS-CoV-2 strains in response to the principal environmental selection pressures like temperature and solar UV radiation is indispensable to overcome COVID-19. To discover the UV-solar radiation-driven genomic adaption of SARS-CoV-2, a curated dataset of 2,500 full-grade genomes from five different UVindex regions (25 countries) was subjected to in-depth downstream genome-wide analysis. The recurrent variants that best respond to UV-solar radiations were extracted and extensively annotated to determine their possible effects and impacts on gene functions. This study revealed 515 recurrent single nucleotide variants (rcntSNVs) as SARS-CoV-2 genomic responses to UV-solar radiation, of which 380 were found to be distinct. For all discovered rcntSNVs, 596 functional effects (rcntEffs) were detected, containing 290 missense, 194 synonymous, 81 regulatory, and 31 in the intergenic region. The highest counts of missense rcntSNVs in spike (27) and nucleocapsid (26) genes explain the SARS-CoV-2 genomic adjustment to escape immunity and prevent UV-induced DNA damage, respectively. Among all, the most commonly observed rcntEffs were four missenses (RdRp-Pro327Leu, N-Arg203Lys, N-Gly204Arg, and Spike-Asp614Gly) and one synonymous (ORF1ab-Phe924Phe) functional effects. The highest number of rcntSNVs found distinct and were uniquely attributed to the specific UVindex regions, proposing solar-UV radiation as one of the driving forces for SARS-CoV-2 differential genomic adaptation. The phylogenetic relationship indicated the high UVindex region populating SARS-CoV-2 as the recent progenitor of all included samples. Altogether, these results provide baseline genomic data that may need to be included for preparing UVindex region-specific future diagnostic and vaccine formulations.

ZEB2 haploinsufficient Mowat-Wilson syndrome induced pluripotent stem cells show disrupted GABAergic transcriptional regulation and function.

Mowat-Wilson syndrome (MWS) is a severe neurodevelopmental disorder caused by heterozygous variants in the gene encoding transcription factor ZEB2. Affected individuals present with structural brain abnormalities, speech delay and epilepsy. In mice, conditional loss of Zeb2 causes hippocampal degeneration, altered migration and differentiation of GABAergic interneurons, a heterogeneous population of mainly inhibitory neurons of importance for maintaining normal excitability. To get insights into GABAergic development and function in MWS we investigated ZEB2 haploinsufficient induced pluripotent stem cells (iPSC) of MWS subjects together with iPSC of healthy donors. Analysis of RNA-sequencing data at two time points of GABAergic development revealed an attenuated interneuronal identity in MWS subject derived iPSC with enrichment of differentially expressed genes required for transcriptional regulation, cell fate transition and forebrain patterning. The ZEB2 haploinsufficient neural stem cells (NSCs) showed downregulation of genes required for ventral telencephalon specification, such as FOXG1, accompanied by an impaired migratory capacity. Further differentiation into GABAergic interneuronal cells uncovered upregulation of transcription factors promoting pallial and excitatory neurons whereas cortical markers were downregulated. The differentially expressed genes formed a neural protein-protein network with extensive connections to well-established epilepsy genes. Analysis of electrophysiological properties in ZEB2 haploinsufficient GABAergic cells revealed overt perturbations manifested as impaired firing of repeated action potentials. Our iPSC model of ZEB2 haploinsufficient GABAergic development thus uncovers a dysregulated gene network leading to immature interneurons with mixed identity and altered electrophysiological properties, suggesting mechanisms contributing to the neuropathogenesis and seizures in MWS.

Syndromic RNA polymerase II insufficiency: Generation of a human induced pluripotent stem cell line (UUIGPi002A-5) with a heterozygous disruption of POLR2A.

Heterozygous variants in POLR2A, encoding the largest subunit of RNA polymerase II, cause severe neurodevelopmental and multisystem abnormalities in humans. Using CRISPR/Cas9 we generated the human iPSC line KICRi002A-5 with a heterozygous truncating 4 bp insertion in exon 5 of the POLR2A gene. Analysis using qRT-PCR confirmed reduced POLR2A mRNA in KICRi002A-5 vs. the isogenic WT iPSC line. The edited iPSC line expressed pluripotency markers and exhibited differentiation capacity into the three germ layers. Assessment of genomic integrity revealed a normal karyotype and OFF-target editing was excluded. The iPSC line KICRi002A-5 provides a useful resource to study mechanisms underlying developmental defects caused by RBP1 insufficiency.

Rare Pathogenic Variants in Genes Implicated in Glutamatergic Neurotransmission Pathway Segregate with Schizophrenia in Pakistani Families.

Schizophrenia is a disabling neuropsychiatric disorder of adulthood onset with high heritability. Worldwide collaborations have identified an association of ~270 common loci, with small individual effects and hence weak clinical implications. The recent technological feasibility of exome sequencing enables the identification of rare variants of high penetrance that refine previous findings and improve risk assessment and prognosis. We recruited two multiplex Pakistani families, having 11 patients and 19 unaffected individuals in three generations. We performed genome-wide SNP genotyping, next-generation mate pairing and whole-exome sequencing of selected members to unveil genetic components. Candidate variants were screened in unrelated cohorts of 508 cases, 300 controls and fifteen families (with 51 affected and 47 unaffected individuals) of Pakistani origin. The structural impact of substituted residues was assessed through in silico modeling using iTASSER. In one family, we identified a rare novel microduplication (5q14.1_q14.2) encompassing critical genes involved in glutamate signaling, such as CMYA5, HOMER and RasGRF2. The second family segregates two ultra-rare, predicted pathogenic variants in the GRIN2A (NM_001134407.3: c.3505C>T, (p.R1169W) and in the NRG3 NM_001010848.4: c.1951G>A, (p.E651K). These genes encode for parts of AMPA and NMDA receptors of glutamatergic neurotransmission, respectively, and the variants are predicted to compromise protein function by destabilizing their structures. The variants were absent in the aforementioned cohorts. Our findings suggest that rare, highly penetrant variants of genes involved in glutamatergic neurotransmission are contributing to the etiology of schizophrenia in these families. It also highlights that genetic investigations of multiplex, multigenerational families could be a powerful approach to identify rare genetic variants involved in complex disorders.