Construction of copy number variation landscape and characterization of associated genes in a Bangladeshi cohort of neurodevelopmental disorders

Introduction: Copy number variations (CNVs) play a critical role in the pathogenesis of neurodevelopmental disorders (NDD) among children. In this study, we aim to identify clinically relevant CNVs, genes and their phenotypic characteristics in an ethnically underrepresented homogenous population of Bangladesh. Methods: We have conducted chromosomal microarray analysis (CMA) for 212 NDD patients with male to female ratio of 2.2:1.0 to identify rare CNVs. To identify candidate genes within the rare CNVs, gene constraint metrics [i.e., “Critical-Exon Genes (CEGs)”] were applied to the population data. Autism Diagnostic Observation Schedule-Second Edition (ADOS-2) was followed in a subset of 95 NDD patients to assess the severity of autism and all statistical tests were performed using the R package. Results: Of all the samples assayed, 12.26% (26/212) and 57.08% (121/212) patients carried pathogenic and variant of uncertain significance (VOUS) CNVs, respectively. While 2.83% (6/212) patients’ pathogenic CNVs were found to be located in the subtelomeric regions. Further burden test identified females are significant carriers of pathogenic CNVs compared to males (OR = 4.2; p = 0.0007). We have observed an increased number of Loss of heterozygosity (LOH) within cases with 23.85% (26/109) consanguineous parents. Our analyses on imprinting genes show, 36 LOH variants disrupting 69 unique imprinted genes and classified these variants as VOUS. ADOS-2 subset shows severe social communication deficit (p = 0.014) and overall ASD symptoms severity (p = 0.026) among the patients carrying duplication CNV compared to the CNV negative group. Candidate gene analysis identified 153 unique CEGs in pathogenic CNVs and 31 in VOUS. Of the unique genes, 18 genes were found to be in smaller (<1 MB) focal CNVs in our NDD cohort and we identified PSMC3 gene as a strong candidate gene for Autism Spectrum Disorder (ASD). Moreover, we hypothesized that KMT2B gene duplication might be associated with intellectual disability. Conclusion: Our results show the utility of CMA for precise genetic diagnosis and its integration into the diagnosis, therapy and management of NDD patients

Numerical data for Nassir et al, 2019

This file contains all numerical data for the plots described in the manuscript, Nassir et al., 2019. Numerical data for each figure is represented in each sheet

Data from: A telomerase with novel non-canonical roles: TERT controls cellular aggregation and tissue size in Dictyostelium

Telomerase, particularly its main subunit, the reverse transcriptase, TERT, prevents DNA erosion during eukaryotic chromosomal replication, but also has poorly understood non-canonical functions. Here, in the model social amoeba Dictyostelium discoideum, we show that the protein encoded by tert has telomerase-like motifs, and regulates, non-canonically, important developmental processes. Expression levels of wild-type (WT) tert were biphasic, peaking at 8 and 12 h post-starvation, aligning with developmental events, such as the initiation of streaming (~7 h) and mound formation (~10 h). In tert KO mutants, however, aggregation was delayed until 16 h. Large, irregular streams formed, then broke up, forming small mounds. The mound-size defect was not induced when a KO mutant of countin (a master size-regulating gene) was treated with TERT inhibitors, but anti-countin antibodies did rescue size in the tert KO. Although, conditioned medium (CM) from countin mutants failed to rescue size in the tert KO, tert KO CM rescued the countin KO phenotype. These and additional observations indicate that TERT acts upstream of smlA/countin: (i) the observed expression levels of smlA and countin, being respectively lower and higher (than WT) in the tert KO; (ii) the levels of known size-regulation intermediates, glucose (low) and adenosine (high), in the tert mutant, and the size defect’s rescue by supplemented glucose or the adenosine-antagonist, caffeine; (iii) the induction of the size defect in the WT by tert KO CM and TERT inhibitors. The tert KO’s other defects (delayed aggregation, irregular streaming) were associated with changes to cAMP-regulated processes (e.g. chemotaxis, cAMP pulsing) and their regulatory factors (e.g. cAMP; acaA, carA expression). Overexpression of WT tert in the tert KO rescued these defects (and size), and restored a single cAMP signaling centre. Our results indicate that TERT acts in novel, non-canonical and upstream ways, regulating key developmental events in Dictyostelium

Neuronal SNCA transcription during Lewy body formation

Abstract Misfolded α-synuclein (α-syn) is believed to contribute to neurodegeneration in Lewy body disease (LBD) based on considerable evidence including a gene-dosage effect observed in relation to point mutations and multiplication of SNCA in familial Parkinson’s disease. A contradictory concept proposes early loss of the physiological α-syn as the major driver of neurodegeneration. There is a paucity of data on SNCA transcripts in various α-syn immunoreactive cytopathologies. Here, the total cell body, nuclear, and cytoplasmic area density of SNCA transcripts in neurons without and with various α-syn immunoreactive cytopathologies in the substantia nigra and amygdala in autopsy cases of LBD (n = 5) were evaluated using RNAscope combined with immunofluorescence for disease-associated α-syn. Single-nucleus RNA sequencing was performed to elucidate cell-type specific SNCA expression in non-diseased frontal cortex (n = 3). SNCA transcripts were observed in the neuronal nucleus and cytoplasm in neurons without α-syn, those containing punctate α-syn immunoreactivity, irregular-shaped compact inclusion, and brainstem-type and cortical-type LBs. However, SNCA transcripts were only rarely found in the α-syn immunoreactive LB areas. The total cell body SNCA transcript area densities in neurons with punctate α-syn immunoreactivity were preserved but were significantly reduced in neurons with compact α-syn inclusions both in the substantia nigra and amygdala. This reduction was also observed in the cytoplasm but not in the nucleus. Only single SNCA transcripts were detected in astrocytes with or without disease-associated α-syn immunoreactivity in the amygdala. Single-nucleus RNA sequencing revealed that excitatory and inhibitory neurons, oligodendrocyte progenitor cells, oligodendrocytes, and homeostatic microglia expressed SNCA transcripts, while expression was largely absent in astrocytes and microglia. The preserved cellular SNCA expression in the more abundant non-Lewy body type α-syn cytopathologies might provide a pool for local protein production that can aggregate and serve as a seed for misfolded α-syn. Successful segregation of disease-associated α-syn is associated with the exhaustion of SNCA production in the terminal cytopathology, the Lewy body. Our observations inform therapy development focusing on targeting SNCA transcription in LBD

Mutational Landscape of Autism Spectrum Disorder Brain Tissue

Rare post-zygotic mutations in the brain are now known to contribute to several neurodevelopmental disorders, including autism spectrum disorder (ASD). However, due to the limited availability of brain tissue, most studies rely on estimates of mosaicism from peripheral samples. In this study, we undertook whole exome sequencing on brain tissue from 26 ASD brain donors from the Harvard Brain Tissue Resource Center (HBTRC) and ascertained the presence of post-zygotic and germline mutations categorized as pathological, including those impacting known ASD-implicated genes. Although quantification did not reveal enrichment for post-zygotic mutations compared with the controls (n = 15), a small number of pathogenic, potentially ASD-implicated mutations were identified, notably in TRAK1 and CLSTN3. Furthermore, germline mutations were identified in the same tissue samples in several key ASD genes, including PTEN, SC1A,&nbsp;CDH13, and CACNA1C. The establishment of tissue resources that are available to the scientific community will facilitate the discovery of new mutations for ASD and other neurodevelopmental disorders

Long-Read Sequencing Improves the Detection of Structural Variations Impacting Complex Non-Coding Elements of the Genome

The advent of long-read sequencing offers a new assessment method of detecting genomic structural variation (SV) in numerous rare genetic diseases. For autism spectrum disorders (ASD) cases where pathogenic variants fail to be found in the protein-coding genic regions along chromosomes, we proposed a scalable workflow to characterize the risk factor of SVs impacting non-coding elements of the genome. We applied whole-genome sequencing on an Emirati family having three children with ASD using long and short-read sequencing technology. A series of analytical pipelines were established to identify a set of SVs with high sensitivity and specificity. At 15-fold coverage, we observed that long-read sequencing technology (987 variants) detected a significantly higher number of SVs when compared to variants detected using short-read technology (509 variants) (p-value &lt; 1.1020 × 10−57). Further comparison showed 97.9% of long-read sequencing variants were spanning within the 1–100 kb size range (p-value &lt; 9.080 × 10−67) and impacting over 5000 genes. Moreover, long-read variants detected 604 non-coding RNAs (p-value &lt; 9.02 × 10−9), comprising 58% microRNA, 31.9% lncRNA, and 9.1% snoRNA. Even at low coverage, long-read sequencing has shown to be a reliable technology in detecting SVs impacting complex elements of the genome

Single-cell transcriptome identifies molecular subtype of autism spectrum disorder impacted by de novo loss-of-function variants regulating glial cells

Abstract Background In recent years, several hundred autism spectrum disorder (ASD) implicated genes have been discovered impacting a wide range of molecular pathways. However, the molecular underpinning of ASD, particularly from the point of view of ‘brain to behaviour’ pathogenic mechanisms, remains largely unknown. Methods We undertook a study to investigate patterns of spatiotemporal and cell type expression of ASD-implicated genes by integrating large-scale brain single-cell transcriptomes (> million cells) and de novo loss-of-function (LOF) ASD variants (impacting 852 genes from 40,122 cases). Results We identified multiple single-cell clusters from three distinct developmental human brain regions (anterior cingulate cortex, middle temporal gyrus and primary visual cortex) that evidenced high evolutionary constraint through enrichment for brain critical exons and high pLI genes. These clusters also showed significant enrichment with ASD loss-of-function variant genes (p < 5.23 × 10–11) that are transcriptionally highly active in prenatal brain regions (visual cortex and dorsolateral prefrontal cortex). Mapping ASD de novo LOF variant genes into large-scale human and mouse brain single-cell transcriptome analysis demonstrate enrichment of such genes into neuronal subtypes and are also enriched for subtype of non-neuronal glial cell types (astrocyte, p < 6.40 × 10–11, oligodendrocyte, p < 1.31 × 10–09). Conclusion Among the ASD genes enriched with pathogenic de novo LOF variants (i.e. KANK1, PLXNB1), a subgroup has restricted transcriptional regulation in non-neuronal cell types that are evolutionarily conserved. This association strongly suggests the involvement of subtype of non-neuronal glial cells in the pathogenesis of ASD and the need to explore other biological pathways for this disorder