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.

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.

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.

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.

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 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.

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.

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.

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.

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.

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.

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.

Modifier Genes in Microcephaly: A Report on WDR62, CEP63, RAD50 and PCNT Variants Exacerbating Disease Caused by Biallelic Mutations of ASPM and CENPJ.

Congenital microcephaly is the clinical presentation of significantly reduced head circumference at birth. It manifests as both non-syndromic-microcephaly primary hereditary (MCPH)-and syndromic forms and shows considerable inter- and intrafamilial variability. It has been hypothesized that additional genetic variants may be responsible for this variability, but data are sparse. We have conducted deep phenotyping and genotyping of five Pakistani multiplex families with either MCPH (n = 3) or Seckel syndrome (n = 2). In addition to homozygous causal variants in ASPM or CENPJ, we discovered additional heterozygous modifier variants in WDR62, CEP63, RAD50 and PCNT-genes already known to be associated with neurological disorders. MCPH patients carrying an additional heterozygous modifier variant showed more severe phenotypic features. Likewise, the phenotype of Seckel syndrome caused by a novel CENPJ variant was aggravated to microcephalic osteodysplastic primordial dwarfism type II (MOPDII) in conjunction with an additional PCNT variant. We show that the CENPJ missense variant impairs splicing and decreases protein expression. We also observed centrosome amplification errors in patient cells, which were twofold higher in MOPDII as compared to Seckel cells. Taken together, these observations advocate for consideration of additional variants in related genes for their role in modifying the expressivity of the phenotype and need to be considered in genetic counseling and risk assessment.

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.

RRP7A links primary microcephaly to dysfunction of ribosome biogenesis, resorption of primary cilia, and neurogenesis.

Primary microcephaly (MCPH) is characterized by reduced brain size and intellectual disability. The exact pathophysiological mechanism underlying MCPH remains to be elucidated, but dysfunction of neuronal progenitors in the developing neocortex plays a major role. We identified a homozygous missense mutation (p.W155C) in Ribosomal RNA Processing 7 Homolog A, RRP7A, segregating with MCPH in a consanguineous family with 10 affected individuals. RRP7A is highly expressed in neural stem cells in developing human forebrain, and targeted mutation of Rrp7a leads to defects in neurogenesis and proliferation in a mouse stem cell model. RRP7A localizes to centrosomes, cilia and nucleoli, and patient-derived fibroblasts display defects in ribosomal RNA processing, primary cilia resorption, and cell cycle progression. Analysis of zebrafish embryos supported that the patient mutation in RRP7A causes reduced brain size, impaired neurogenesis and cell proliferation, and defective ribosomal RNA processing. These findings provide novel insight into human brain development and MCPH.

Generation of two human iPSC lines (UUIGPi013-A and UUIPGi014-A) from cases with Down syndrome and full trisomy for chromosome 21 (T21).

Down syndrome (DS) is caused by trisomy for chromosome 21 (T21). We generated two induced pluripotent stem cell (iPSC) lines from skin fibroblasts of two males with DS using Sendai virus delivery of OCT4, SOX2, KLF4, and c-MYC. Characterization of the two iPSC lines, UUIGPi013-A and UUIPGi014-A, showed that they are genetically stable with a 47,XY,+21 karyotype. Both lines displayed expression of pluripotency markers and trilineage differentiation capacity. These two iPSC lines provide a useful resource for DS modeling and pharmacological interventions.

Aberrant splicing due to a novel RPS7 variant causes Diamond-Blackfan Anemia associated with spontaneous remission and meningocele.

Diamond-Blackfan Anemia (DBA) is a congenital pure red cell aplasia caused by heterozygous variants in ribosomal protein genes. The hematological features associated with DBA are highly variable and non-hematological abnormalities are common. We report herein on an affected mother and her daughter presenting with transfusion-dependent anemia. The mother showed mild physical abnormalities and entered spontaneous remission at age 13 years. Her daughter was born with occipital meningocele. Exome sequencing of DNA from the mother revealed a heterozygous novel splice site variant (NM_001011.4:c.508-3T > G) in the Ribosomal Protein S7 gene (RPS7) inherited by the daughter. Functional analysis of the RPS7 variant expressed from a mini-gene construct revealed that the exon 7 acceptor splice site was replaced by a cryptic splice resulting in a transcript missing 64 bp of exon 7 (p.Val170Serfs*8). Our study confirms a pathogenic effect of a novel RPS7 variant in DBA associated with spontaneous remission in the mother and meningocele in her daughter, thus adding to the genotype-phenotype correlations in DBA.

Incontinentia pigmenti: Generation of an IKBKG deficient human iPSC line (KICRi002-A-1) on a 46,XY background using CRISPR/Cas9.

Incontinentia pigmenti (IP) is an X-linked dominant neuroectodermal dysplasia caused by loss-of-function mutations in the IKBKG gene. Using CRISPR/Cas9 technology, we generated an IKBKG knock-out iPSC line (KICRi002-A-1) on a 46,XY background. The iPSC line showed a normal karyotype, expressed pluripotency markers and exhibited capability to differentiate into the three germ layers in vitro. Off-target editing was excluded and no IKBKG mRNA expression could be detected. Our line offers a useful resource to elucidate mechanisms caused by IKBKG deficiency that leads to disrupted male fetal development and for drug screening to improve treatment of female patients with IP.

Generation of a human Neurochondrin deficient iPSC line KICRi002-A-3 using CRISPR/Cas9.

The role of Neurochondrin (NCDN) in humans is not well understood. Mice with a conditional Ncdn knock-out show epileptic seizures, depressive-like behaviours and impaired spatial learning. Using CRISPR/Cas9, we generated a Neurochondrin deficient human iPSC line KICRi002-A-3 carrying a homozygous 752 bp deletion / 2 bp insertion in the NCDN gene. The iPSC line maintained a normal 46,XY karyotype, expressed pluripotency markers and exhibited capability to differentiate into the three germ layers in vitro. Off-target editing was excluded and Neurochondrin expression was not detectable. The iPSC line offers a valuable resource to study the role of Neurochondrin during human neurogenesis.