Capitalizing on the heterogeneous effects of CFTR nonsense and frameshift variants to inform therapeutic strategy for cystic fibrosis.

CFTR modulators have revolutionized the treatment of individuals with cystic fibrosis (CF) by improving the function of existing protein. Unfortunately, almost half of the disease-causing variants in CFTR are predicted to introduce premature termination codons (PTC) thereby causing absence of full-length CFTR protein. We hypothesized that a subset of nonsense and frameshift variants in CFTR allow expression of truncated protein that might respond to FDA-approved CFTR modulators. To address this concept, we selected 26 PTC-generating variants from four regions of CFTR and determined their consequences on CFTR mRNA, protein and function using intron-containing minigenes expressed in 3 cell lines (HEK293, MDCK and CFBE41o-) and patient-derived conditionally reprogrammed primary nasal epithelial cells. The PTC-generating variants fell into five groups based on RNA and protein effects. Group A (reduced mRNA, immature (core glycosylated) protein, function <1% (n = 5)) and Group B (normal mRNA, immature protein, function <1% (n = 10)) variants were unresponsive to modulator treatment. However, Group C (normal mRNA, mature (fully glycosylated) protein, function >1% (n = 5)), Group D (reduced mRNA, mature protein, function >1% (n = 5)) and Group E (aberrant RNA splicing, mature protein, function > 1% (n = 1)) variants responded to modulators. Increasing mRNA level by inhibition of NMD led to a significant amplification of modulator effect upon a Group D variant while response of a Group A variant was unaltered. Our work shows that PTC-generating variants should not be generalized as genetic 'nulls' as some may allow generation of protein that can be targeted to achieve clinical benefit.

Systematic Computational Identification of Variants That Activate Exonic and Intronic Cryptic Splice Sites.

We developed a variant-annotation method that combines sequence-based machine-learning classification with a context-dependent algorithm for selecting splice variants. Our approach is distinctive in that it compares the splice potential of a sequence bearing a variant with the splice potential of the reference sequence. After training, classification accurately identified 168 of 180 (93.3%) canonical splice sites of five genes. The combined method, CryptSplice, identified and correctly predicted the effect of 18 of 21 (86%) known splice-altering variants in CFTR, a well-studied gene whose loss-of-function variants cause cystic fibrosis (CF). Among 1,423 unannotated CFTR disease-associated variants, the method identified 32 potential exonic cryptic splice variants, two of which were experimentally evaluated and confirmed. After complete CFTR sequencing, the method found three cryptic intronic splice variants (one known and two experimentally verified) that completed the molecular diagnosis of CF in 6 of 14 individuals. CryptSplice interrogation of sequence data from six individuals with X-linked dyskeratosis congenita caused by an unknown disease-causing variant in DKC1 identified two splice-altering variants that were experimentally verified. To assess the extent to which disease-associated variants might activate cryptic splicing, we selected 458 pathogenic variants and 348 variants of uncertain significance (VUSs) classified as high confidence from ClinVar. Splice-site activation was predicted for 129 (28%) of the pathogenic variants and 75 (22%) of the VUSs. Our findings suggest that cryptic splice-site activation is more common than previously thought and should be routinely considered for all variants within the transcribed regions of genes.

Integrated genome and transcriptome sequencing identifies a noncoding mutation in the genome replication factor DONSON as the cause of microcephaly-micromelia syndrome.

While next-generation sequencing has accelerated the discovery of human disease genes, progress has been largely limited to the "low hanging fruit" of mutations with obvious exonic coding or canonical splice site impact. In contrast, the lack of high-throughput, unbiased approaches for functional assessment of most noncoding variants has bottlenecked gene discovery. We report the integration of transcriptome sequencing (RNA-seq), which surveys all mRNAs to reveal functional impacts of variants at the transcription level, into the gene discovery framework for a unique human disease, microcephaly-micromelia syndrome (MMS). MMS is an autosomal recessive condition described thus far in only a single First Nations population and causes intrauterine growth restriction, severe microcephaly, craniofacial anomalies, skeletal dysplasia, and neonatal lethality. Linkage analysis of affected families, including a very large pedigree, identified a single locus on Chromosome 21 linked to the disease (LOD > 9). Comprehensive genome sequencing did not reveal any pathogenic coding or canonical splicing mutations within the linkage region but identified several nonconserved noncoding variants. RNA-seq analysis detected aberrant splicing in DONSON due to one of these noncoding variants, showing a causative role for DONSON disruption in MMS. We show that DONSON is expressed in progenitor cells of embryonic human brain and other proliferating tissues, is co-expressed with components of the DNA replication machinery, and that Donson is essential for early embryonic development in mice as well, suggesting an essential conserved role for DONSON in the cell cycle. Our results demonstrate the utility of integrating transcriptomics into the study of human genetic disease when DNA sequencing alone is not sufficient to reveal the underlying pathogenic mutation.

Mutations in mitochondrial enzyme GPT2 cause metabolic dysfunction and neurological disease with developmental and progressive features.

Mutations that cause neurological phenotypes are highly informative with regard to mechanisms governing human brain function and disease. We report autosomal recessive mutations in the enzyme glutamate pyruvate transaminase 2 (GPT2) in large kindreds initially ascertained for intellectual and developmental disability (IDD). GPT2 [also known as alanine transaminase 2 (ALT2)] is one of two related transaminases that catalyze the reversible addition of an amino group from glutamate to pyruvate, yielding alanine and α-ketoglutarate. In addition to IDD, all affected individuals show postnatal microcephaly and ∼80% of those followed over time show progressive motor symptoms, a spastic paraplegia. Homozygous nonsense p.Arg404* and missense p.Pro272Leu mutations are shown biochemically to be loss of function. The GPT2 gene demonstrates increasing expression in brain in the early postnatal period, and GPT2 protein localizes to mitochondria. Akin to the human phenotype, Gpt2-null mice exhibit reduced brain growth. Through metabolomics and direct isotope tracing experiments, we find a number of metabolic abnormalities associated with loss of Gpt2. These include defects in amino acid metabolism such as low alanine levels and elevated essential amino acids. Also, we find defects in anaplerosis, the metabolic process involved in replenishing TCA cycle intermediates. Finally, mutant brains demonstrate misregulated metabolites in pathways implicated in neuroprotective mechanisms previously associated with neurodegenerative disorders. Overall, our data reveal an important role for the GPT2 enzyme in mitochondrial metabolism with relevance to developmental as well as potentially to neurodegenerative mechanisms.

Mutations in QARS, encoding glutaminyl-tRNA synthetase, cause progressive microcephaly, cerebral-cerebellar atrophy, and intractable seizures.

Progressive microcephaly is a heterogeneous condition with causes including mutations in genes encoding regulators of neuronal survival. Here, we report the identification of mutations in QARS (encoding glutaminyl-tRNA synthetase [QARS]) as the causative variants in two unrelated families affected by progressive microcephaly, severe seizures in infancy, atrophy of the cerebral cortex and cerebellar vermis, and mild atrophy of the cerebellar hemispheres. Whole-exome sequencing of individuals from each family independently identified compound-heterozygous mutations in QARS as the only candidate causative variants. QARS was highly expressed in the developing fetal human cerebral cortex in many cell types. The four QARS mutations altered highly conserved amino acids, and the aminoacylation activity of QARS was significantly impaired in mutant cell lines. Variants p.Gly45Val and p.Tyr57His were located in the N-terminal domain required for QARS interaction with proteins in the multisynthetase complex and potentially with glutamine tRNA, and recombinant QARS proteins bearing either substitution showed an over 10-fold reduction in aminoacylation activity. Conversely, variants p.Arg403Trp and p.Arg515Trp, each occurring in a different family, were located in the catalytic core and completely disrupted QARS aminoacylation activity in vitro. Furthermore, p.Arg403Trp and p.Arg515Trp rendered QARS less soluble, and p.Arg403Trp disrupted QARS-RARS (arginyl-tRNA synthetase 1) interaction. In zebrafish, homozygous qars loss of function caused decreased brain and eye size and extensive cell death in the brain. Our results highlight the importance of QARS during brain development and that epilepsy due to impairment of QARS activity is unusually severe in comparison to other aminoacyl-tRNA synthetase disorders.

Somatic mutations in cerebral cortical malformations.

BACKGROUND: Although there is increasing recognition of the role of somatic mutations in genetic disorders, the prevalence of somatic mutations in neurodevelopmental disease and the optimal techniques to detect somatic mosaicism have not been systematically evaluated. METHODS: Using a customized panel of known and candidate genes associated with brain malformations, we applied targeted high-coverage sequencing (depth, ≥200×) to leukocyte-derived DNA samples from 158 persons with brain malformations, including the double-cortex syndrome (subcortical band heterotopia, 30 persons), polymicrogyria with megalencephaly (20), periventricular nodular heterotopia (61), and pachygyria (47). We validated candidate mutations with the use of Sanger sequencing and, for variants present at unequal read depths, subcloning followed by colony sequencing. RESULTS: Validated, causal mutations were found in 27 persons (17%; range, 10 to 30% for each phenotype). Mutations were somatic in 8 of the 27 (30%), predominantly in persons with the double-cortex syndrome (in whom we found mutations in DCX and LIS1), persons with periventricular nodular heterotopia (FLNA), and persons with pachygyria (TUBB2B). Of the somatic mutations we detected, 5 (63%) were undetectable with the use of traditional Sanger sequencing but were validated through subcloning and subsequent sequencing of the subcloned DNA. We found potentially causal mutations in the candidate genes DYNC1H1, KIF5C, and other kinesin genes in persons with pachygyria. CONCLUSIONS: Targeted sequencing was found to be useful for detecting somatic mutations in patients with brain malformations. High-coverage sequencing panels provide an important complement to whole-exome and whole-genome sequencing in the evaluation of somatic mutations in neuropsychiatric disease. (Funded by the National Institute of Neurological Disorders and Stroke and others.).

CRISPR Del/Rei: a simple, flexible, and efficient pipeline for scarless genome editing.

Scarless genome editing of induced pluripotent stem cells (iPSCs) is crucial for the precise modeling of genetic disease. Here we present CRISPR Del/Rei, a two-step deletion-reinsertion strategy with high editing efficiency and simple PCR-based screening that generates isogenic clones in ~ 2 months. We apply our strategy to edit iPSCs at 3 loci with only rare off target editing.

Transcriptomic data of Clozapine-treated Ngn2-induced Human Excitatory Neurons.

We generated human excitatory neurons using a protocol for rapid 21-day induction using neurogenin-2 overexpression (Zhang et al., 2013) in a publicly available control iPSC line. We validated the glutamatergic neuronal identity of the neurons by immunofluorescence and transcriptomics. We exposed 6 of the 12 replicate neuron cultures to therapeutic plasma levels of clozapine (300 ng/mL) for the last 3 days of culture, and the remaining 6 to replicates to the clozapine solvent alone (methanol) to be used as controls. We harvested the cultures and extracted total RNA, depleted ribosomal RNA and subjected them to RNA sequencing. Of the 6 control replicates 2 failed RNA quality control, and thus a total of 6 exposed and 4 control cultures were used for further analysis. Here, we provide that raw sequencing data as well as a list of all of the genes and their expression levels resulting from the RNA-sequencing. This dataset can be used as a reference data for future studies that access additional neuronal cell types, clozapine exposure conditions, and other antipsychotic medication. Related Research Article: Das, D., Peng, X., Lam, A.N., Bader, J.S., Avramopoulos, D., 2021. Transcriptome analysis of human induced excitatory neurons supports a strong effect of clozapine on cholesterol biosynthesis. Schizophr Res 228, 324-326. (Das et al., 2021).

Transcriptome analysis of human induced excitatory neurons supports a strong effect of clozapine on cholesterol biosynthesis.

Antipsychotics are known to modulate dopamine and other neurotransmitters which is often thought to be the mechanism underlying their therapeutic effects. Nevertheless, other less studied consequences of antipsychotics on neuronal function may contribute to their efficacy. Revealing the complete picture behind their action is of paramount importance for precision medicine and accurate drug selection. Progress in cell engineering allows the generation of induced pluripotent stem cells (iPSCs) and their differentiation to a variety of neuronal types, providing new tools to study antipsychotics. Here we use excitatory cortical neurons derived from iPSCs to explore their response to therapeutic levels of Clozapine as measured by their transcriptomic output, a proxy for neuronal homeostasis. To our surprise, but in agreement with the results of many investigators studying glial-like cells, Clozapine had a very strong effect on cholesterol metabolism. More than a quarter (12) of all annotated cholesterol genes (46) in the genome were significantly changed at FDR < 0.1, all upregulated. This is a 35-fold enrichment with an adjusted p = 8 × 10-11. Notably no other functional category showed evidence of enrichment. Cholesterol is a major component of the neuronal membrane and myelin but it does not cross the blood brain barrier, it is produced locally mostly by glia but also by neurons. By singling out increased expression of cholesterol metabolism genes as the main response of cortical excitatory neurons to antipsychotics, our work supports the hypothesis that cholesterol metabolism may be a contributing mechanism to the beneficial effects of Clozapine and possibly other antipsychotics.

SLC26A9 SNP rs7512462 is not associated with lung disease severity or lung function response to ivacaftor in cystic fibrosis patients with G551D-CFTR.

BACKGROUND: The CFTR modulator ivacaftor has been variably effective in treating individuals with cystic fibrosis (CF) who harbor CFTR gating variants such as G551D, as well as other classes of CFTR variants when used with other modulators. Because CFTR genotype does not fully explain this variability, defining genetic modifiers of response to modulator therapy is of particular interest to the field of individualized CF drug therapy. Previous studies have proposed that a variant in SLC26A9 (rs7512462) is associated with lung disease severity and with response to treatment with ivacaftor in individuals with CF who carry G551D or gating variants. METHODS: Given the implications for CF treatment, we re-examined the reported associations in three cohorts; patients enrolled in the Twin and Siblings study at Johns Hopkins University, the CF modifier study at the University of North Carolina at Chapel Hill, and the prospective G551D Observational (GOAL) study. The GOAL study was specifically designed to measure lung function response to ivacaftor. RESULTS: We find no association between SLC26A9 (rs7512462) genotype and lung disease severity (n = 272) or change in lung function at one-, three-, and six-month intervals following ivacaftor treatment(n = 141) in individuals with CF who carry at least one G551D variant. CONCLUSIONS: Our inability to replicate this association indicates that rs7512462 genotype should not be used in treatment decisions.

Increased expression of anion transporter SLC26A9 delays diabetes onset in cystic fibrosis.

Diabetes is a common complication of cystic fibrosis (CF) that affects approximately 20% of adolescents and 40%-50% of adults with CF. The age at onset of CF-related diabetes (CFRD) (marked by clinical diagnosis and treatment initiation) is an important measure of the disease process. DNA variants associated with age at onset of CFRD reside in and near SLC26A9. Deep sequencing of the SLC26A9 gene in 762 individuals with CF revealed that 2 common DNA haplotypes formed by the risk variants account for the association with diabetes. Single-cell RNA sequencing (scRNA-Seq) indicated that SLC26A9 is predominantly expressed in pancreatic ductal cells and frequently coexpressed with CF transmembrane conductance regulator (CFTR) along with transcription factors that have binding sites 5' of SLC26A9. These findings were replicated upon reanalysis of scRNA-Seq data from 4 independent studies. DNA fragments derived from the 5' region of SLC26A9-bearing variants from the low-risk haplotype generated 12%-20% higher levels of expression in PANC-1 and CFPAC-1 cells compared with the high- risk haplotype. Taken together, our findings indicate that an increase in SLC26A9 expression in ductal cells of the pancreas delays the age at onset of diabetes, suggesting a CFTR-agnostic treatment for a major complication of CF.

Homozygous deletions implicate non-coding epigenetic marks in Autism spectrum disorder.

More than 98% of the human genome is made up of non-coding DNA, but techniques to ascertain its contribution to human disease have lagged far behind our understanding of protein coding variations. Autism spectrum disorder (ASD) has been mostly associated with coding variations via de novo single nucleotide variants (SNVs), recessive/homozygous SNVs, or de novo copy number variants (CNVs); however, most ASD cases continue to lack a genetic diagnosis. We analyzed 187 consanguineous ASD families for biallelic CNVs. Recessive deletions were significantly enriched in affected individuals relative to their unaffected siblings (17% versus 4%, p < 0.001). Only a small subset of biallelic deletions were predicted to result in coding exon disruption. In contrast, biallelic deletions in individuals with ASD were enriched for overlap with regulatory regions, with 23/28 CNVs disrupting histone peaks in ENCODE (p < 0.009). Overlap with regulatory regions was further demonstrated by comparisons to the 127-epigenome dataset released by the Roadmap Epigenomics project, with enrichment for enhancers found in primary brain tissue and neuronal progenitor cells. Our results suggest a novel noncoding mechanism of ASD, describe a powerful method to identify important noncoding regions in the human genome, and emphasize the potential significance of gene activation and regulation in cognitive and social function.

Biallelic mutations in human DCC cause developmental split-brain syndrome.

Motor, sensory, and integrative activities of the brain are coordinated by a series of midline-bridging neuronal commissures whose development is tightly regulated. Here we report a new human syndrome in which these commissures are widely disrupted, thus causing clinical manifestations of horizontal gaze palsy, scoliosis, and intellectual disability. Affected individuals were found to possess biallelic loss-of-function mutations in the gene encoding the axon-guidance receptor 'deleted in colorectal carcinoma' (DCC), which has been implicated in congenital mirror movements when it is mutated in the heterozygous state but whose biallelic loss-of-function human phenotype has not been reported. Structural MRI and diffusion tractography demonstrated broad disorganization of white-matter tracts throughout the human central nervous system (CNS), including loss of all commissural tracts at multiple levels of the neuraxis. Combined with data from animal models, these findings show that DCC is a master regulator of midline crossing and development of white-matter projections throughout the human CNS.

Targeted DNA Sequencing from Autism Spectrum Disorder Brains Implicates Multiple Genetic Mechanisms.

Single nucleotide variants (SNVs), particularly loss-of-function mutations, are significant contributors to autism spectrum disorder (ASD) risk. Here we report the first systematic deep sequencing study of 55 postmortem ASD brains for SNVs in 78 known ASD candidate genes. Remarkably, even without parental samples, we find more ASD brains with mutations that are protein-altering (26/55 cases versus 12/50 controls, p = 0.015), deleterious (16/55 versus 5/50, p = 0.016), or loss-of-function (6/55 versus 0/50, p = 0.028) compared to controls, with recurrent deleterious mutations in ARID1B, SCN1A, SCN2A, and SETD2, suggesting these mutations contribute to ASD risk. In several cases, the identified mutations and medical records suggest syndromic ASD diagnoses. Two ASD and one Fragile X premutation case showed deleterious somatic mutations, providing evidence that somatic mutations occur in ASD cases, and supporting a model in which a combination of germline and/or somatic mutations may contribute to ASD risk on a case-by-case basis.

Deletions in GRID2 lead to a recessive syndrome of cerebellar ataxia and tonic upgaze in humans.

OBJECTIVE: To identify the genetic cause of a syndrome causing cerebellar ataxia and eye movement abnormalities. METHODS: We identified 2 families with cerebellar ataxia, eye movement abnormalities, and global developmental delay. We performed genetic analyses including single nucleotide polymorphism genotyping, linkage analysis, array comparative genomic hybridization, quantitative PCR, and Sanger sequencing. We obtained eye movement recordings of mutant mice deficient for the ortholog of the identified candidate gene, and performed immunohistochemistry using human and mouse brain specimens. RESULTS: All affected individuals had ataxia, eye movement abnormalities, most notably tonic upgaze, and delayed speech and cognitive development. Homozygosity mapping identified the disease locus on chromosome 4q. Within this region, a homozygous deletion of GRID2 exon 4 in the index family and compound heterozygous deletions involving GRID2 exon 2 in the second family were identified. Grid2-deficient mice showed larger spontaneous and random eye movements compared to wild-type mice. In developing mouse and human cerebella, GRID2 localized to the Purkinje cell dendritic spines. Brain MRI in 2 affected children showed progressive cerebellar atrophy, which was more severe than that of Grid2-deficient mice. CONCLUSIONS: Biallelic deletions of GRID2 lead to a syndrome of cerebellar ataxia and tonic upgaze in humans. The phenotypic resemblance and similarity in protein expression pattern between humans and mice suggest a conserved role for GRID2 in the synapse organization between parallel fibers and Purkinje cells. However, the progressive and severe cerebellar atrophy seen in the affected individuals could indicate an evolutionarily unique role for GRID2 in the human cerebellum.