Dissecting the Causal Mechanism of X-Linked Dystonia-Parkinsonism by Integrating Genome and Transcriptome Assembly.

X-linked Dystonia-Parkinsonism (XDP) is a Mendelian neurodegenerative disease that is endemic to the Philippines and is associated with a founder haplotype. We integrated multiple genome and transcriptome assembly technologies to narrow the causal mutation to the TAF1 locus, which included a SINE-VNTR-Alu (SVA) retrotransposition into intron 32 of the gene. Transcriptome analyses identified decreased expression of the canonical cTAF1 transcript among XDP probands, and de novo assembly across multiple pluripotent stem-cell-derived neuronal lineages discovered aberrant TAF1 transcription that involved alternative splicing and intron retention (IR) in proximity to the SVA that was anti-correlated with overall TAF1 expression. CRISPR/Cas9 excision of the SVA rescued this XDP-specific transcriptional signature and normalized TAF1 expression in probands. These data suggest an SVA-mediated aberrant transcriptional mechanism associated with XDP and may provide a roadmap for layered technologies and integrated assembly-based analyses for other unsolved Mendelian disorders.

Systematic evaluation of genome sequencing for the diagnostic assessment of autism spectrum disorder and fetal structural anomalies.

Short-read genome sequencing (GS) holds the promise of becoming the primary diagnostic approach for the assessment of autism spectrum disorder (ASD) and fetal structural anomalies (FSAs). However, few studies have comprehensively evaluated its performance against current standard-of-care diagnostic tests: karyotype, chromosomal microarray (CMA), and exome sequencing (ES). To assess the clinical utility of GS, we compared its diagnostic yield against these three tests in 1,612 quartet families including an individual with ASD and in 295 prenatal families. Our GS analytic framework identified a diagnostic variant in 7.8% of ASD probands, almost 2-fold more than CMA (4.3%) and 3-fold more than ES (2.7%). However, when we systematically captured copy-number variants (CNVs) from the exome data, the diagnostic yield of ES (7.4%) was brought much closer to, but did not surpass, GS. Similarly, we estimated that GS could achieve an overall diagnostic yield of 46.1% in unselected FSAs, representing a 17.2% increased yield over karyotype, 14.1% over CMA, and 4.1% over ES with CNV calling or 36.1% increase without CNV discovery. Overall, GS provided an added diagnostic yield of 0.4% and 0.8% beyond the combination of all three standard-of-care tests in ASD and FSAs, respectively. This corresponded to nine GS unique diagnostic variants, including sequence variants in exons not captured by ES, structural variants (SVs) inaccessible to existing standard-of-care tests, and SVs where the resolution of GS changed variant classification. Overall, this large-scale evaluation demonstrated that GS significantly outperforms each individual standard-of-care test while also outperforming the combination of all three tests, thus warranting consideration as the first-tier diagnostic approach for the assessment of ASD and FSAs.

Transcriptional and functional consequences of alterations to MEF2C and its topological organization in neuronal models.

Point mutations and structural variants that directly disrupt the coding sequence of MEF2C have been associated with a spectrum of neurodevelopmental disorders (NDDs). However, the impact of MEF2C haploinsufficiency on neurodevelopmental pathways and synaptic processes is not well understood, nor are the complex mechanisms that govern its regulation. To explore the functional changes associated with structural variants that alter MEF2C expression and/or regulation, we generated an allelic series of 204 isogenic human induced pluripotent stem cell (hiPSC)-derived neural stem cells and glutamatergic induced neurons. These neuronal models harbored CRISPR-engineered mutations that involved direct deletion of MEF2C or deletion of the boundary points for topologically associating domains (TADs) and chromatin loops encompassing MEF2C. Systematic profiling of mutation-specific alterations, contrasted to unedited controls that were exposed to the same guide RNAs for each edit, revealed that deletion of MEF2C caused differential expression of genes associated with neurodevelopmental pathways and synaptic function. We also discovered significant reduction in synaptic activity measured by multielectrode arrays (MEAs) in neuronal cells. By contrast, we observed robust buffering against MEF2C regulatory disruption following deletion of a distal 5q14.3 TAD and loop boundary, whereas homozygous loss of a proximal loop boundary resulted in down-regulation of MEF2C expression and reduced electrophysiological activity on MEA that was comparable to direct gene disruption. Collectively, these studies highlight the considerable functional impact of MEF2C deletion in neuronal cells and systematically characterize the complex interactions that challenge a priori predictions of regulatory consequences from structural variants that disrupt three-dimensional genome organization.

Tissue- and cell-type-specific molecular and functional signatures of 16p11.2 reciprocal genomic disorder across mouse brain and human neuronal models.

Chromosome 16p11.2 reciprocal genomic disorder, resulting from recurrent copy-number variants (CNVs), involves intellectual disability, autism spectrum disorder (ASD), and schizophrenia, but the responsible mechanisms are not known. To systemically dissect molecular effects, we performed transcriptome profiling of 350 libraries from six tissues (cortex, cerebellum, striatum, liver, brown fat, and white fat) in mouse models harboring CNVs of the syntenic 7qF3 region, as well as cellular, transcriptional, and single-cell analyses in 54 isogenic neural stem cell, induced neuron, and cerebral organoid models of CRISPR-engineered 16p11.2 CNVs. Transcriptome-wide differentially expressed genes were largely tissue-, cell-type-, and dosage-specific, although more effects were shared between deletion and duplication and across tissue than expected by chance. The broadest effects were observed in the cerebellum (2,163 differentially expressed genes), and the greatest enrichments were associated with synaptic pathways in mouse cerebellum and human induced neurons. Pathway and co-expression analyses identified energy and RNA metabolism as shared processes and enrichment for ASD-associated, loss-of-function constraint, and fragile X messenger ribonucleoprotein target gene sets. Intriguingly, reciprocal 16p11.2 dosage changes resulted in consistent decrements in neurite and electrophysiological features, and single-cell profiling of organoids showed reciprocal alterations to the proportions of excitatory and inhibitory GABAergic neurons. Changes both in neuronal ratios and in gene expression in our organoid analyses point most directly to calretinin GABAergic inhibitory neurons and the excitatory/inhibitory balance as targets of disruption that might contribute to changes in neurodevelopmental and cognitive function in 16p11.2 carriers. Collectively, our data indicate the genomic disorder involves disruption of multiple contributing biological processes and that this disruption has relative impacts that are context specific.

Dystonia-specific mutations in THAP1 alter transcription of genes associated with neurodevelopment and myelin.

Dystonia is a neurologic disorder associated with an increasingly large number of genetic variants in many genes, resulting in characteristic disturbances in volitional movement. Dissecting the relationships between these mutations and their functional outcomes is critical in understanding the pathways that drive dystonia pathogenesis. Here we established a pipeline for characterizing an allelic series of dystonia-specific mutations. We used this strategy to investigate the molecular consequences of genetic variation in THAP1, which encodes a transcription factor linked to neural differentiation. Multiple pathogenic mutations associated with dystonia cluster within distinct THAP1 functional domains and are predicted to alter DNA-binding properties and/or protein interactions differently, yet the relative impact of these varied changes on molecular signatures and neural deficits is unclear. To determine the effects of these mutations on THAP1 transcriptional activity, we engineered an allelic series of eight alterations in a common induced pluripotent stem cell background and differentiated these lines into a panel of near-isogenic neural stem cells (n = 94 lines). Transcriptome profiling followed by joint analysis of the most robust signatures across mutations identified a convergent pattern of dysregulated genes functionally related to neurodevelopment, lysosomal lipid metabolism, and myelin. On the basis of these observations, we examined mice bearing Thap1-disruptive alleles and detected significant changes in myelin gene expression and reduction of myelin structural integrity relative to control mice. These results suggest that deficits in neurodevelopment and myelination are common consequences of dystonia-associated THAP1 mutations and highlight the potential role of neuron-glial interactions in the pathogenesis of dystonia.

Kctd13-deficient mice display short-term memory impairment and sex-dependent genetic interactions.

The 16p11.2 BP4-BP5 deletion and duplication syndromes are associated with a complex spectrum of neurodevelopmental phenotypes that includes developmental delay and autism spectrum disorder, with a reciprocal effect on head circumference, brain structure and body mass index. Mouse models of the 16p11.2 copy number variant have recapitulated some of the patient phenotypes, while studies in flies and zebrafish have uncovered several candidate contributory genes within the region, as well as complex genetic interactions. We evaluated one of these loci, KCTD13, by modeling haploinsufficiency and complete knockout in mice. In contrast to the zebrafish model, and in agreement with recent data, we found normal brain structure in heterozygous and homozygous mutants. However, recapitulating previously observed genetic interactions, we discovered sex-specific brain volumetric alterations in double heterozygous Kctd13xMvp and Kctd13xLat mice. Behavioral testing revealed a significant deficit in novel object recognition, novel location recognition and social transmission of food preference in Kctd13 mutants. These phenotypes were concomitant with a reduction in density of mature spines in the hippocampus, but potentially independent of RhoA abundance, which was unperturbed postnatally in our mutants. Furthermore, transcriptome analyses from cortex and hippocampus highlighted the dysregulation of pathways important in neurodevelopment, the most significant of which was synaptic formation. Together, these data suggest that KCTD13 contributes to the neurocognitive aspects of patients with the BP4-BP5 deletion, likely through genetic interactions with other loci.

TAF1 Transcripts and Neurofilament Light Chain as Biomarkers for X-linked Dystonia-Parkinsonism.

BACKGROUND: X-linked dystonia-parkinsonism is a rare neurological disease endemic to the Philippines. Dystonic symptoms appear in males at the mean age of 40 years and progress to parkinsonism with degenerative pathology in the striatum. A retrotransposon inserted in intron 32 of the TAF1 gene leads to alternative splicing in the region and a reduction of the full-length mRNA transcript. OBJECTIVES: The objective of this study was to discover cell-based and biofluid-based biomarkers for X-linked dystonia-parkinsonism. METHODS: RNA from patient-derived neural progenitor cells and their secreted extracellular vesicles were used to screen for dysregulation of TAF1 expression. Droplet-digital polymerase chain reaction was used to quantify the expression of TAF1 mRNA fragments 5' and 3' to the retrotransposon insertion and the disease-specific splice variant TAF1-32i in whole-blood RNA. Plasma levels of neurofilament light chain were measured using single-molecule array. RESULTS: In neural progenitor cells and their extracellular vesicles, we confirmed that the TAF1-3'/5' ratio was lower in patient samples, whereas TAF1-32i expression is higher relative to controls. In whole-blood RNA, both TAF1-3'/5' ratio and TAF1-32i expression can differentiate patient (n = 44) from control samples (n = 18) with high accuracy. Neurofilament light chain plasma levels were significantly elevated in patients (n = 43) compared with both carriers (n = 16) and controls (n = 21), with area under the curve of 0.79. CONCLUSIONS: TAF1 dysregulation in blood serves as a disease-specific biomarker that could be used as a readout for monitoring therapies targeting TAF1 splicing. Neurofilament light chain could be used in monitoring neurodegeneration and disease progression in patients. © 2020 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Loss of LDAH associated with prostate cancer and hearing loss.

Great strides in gene discovery have been made using a multitude of methods to associate phenotypes with genetic variants, but there still remains a substantial gap between observed symptoms and identified genetic defects. Herein, we use the convergence of various genetic and genomic techniques to investigate the underpinnings of a constellation of phenotypes that include prostate cancer (PCa) and sensorineural hearing loss (SNHL) in a human subject. Through interrogation of the subject's de novo, germline, balanced chromosomal translocation, we first identify a correlation between his disorders and a poorly annotated gene known as lipid droplet associated hydrolase (LDAH). Using data repositories of both germline and somatic variants, we identify convergent genomic evidence that substantiates a correlation between loss of LDAH and PCa. This correlation is validated through both in vitro and in vivo models that show loss of LDAH results in increased risk of PCa and, to a lesser extent, SNHL. By leveraging convergent evidence in emerging genomic data, we hypothesize that loss of LDAH is involved in PCa and other phenotypes observed in support of a genotype-phenotype association in an n-of-one human subject.

ROR1 is essential for proper innervation of auditory hair cells and hearing in humans and mice.

Hair cells of the inner ear, the mechanosensory receptors, convert sound waves into neural signals that are passed to the brain via the auditory nerve. Little is known about the molecular mechanisms that govern the development of hair cell-neuronal connections. We ascertained a family with autosomal recessive deafness associated with a common cavity inner ear malformation and auditory neuropathy. Via whole-exome sequencing, we identified a variant (c.2207G>C, p.R736T) in ROR1 (receptor tyrosine kinase-like orphan receptor 1), cosegregating with deafness in the family and absent in ethnicity-matched controls. ROR1 is a tyrosine kinase-like receptor localized at the plasma membrane. At the cellular level, the mutation prevents the protein from reaching the cellular membrane. In the presence of WNT5A, a known ROR1 ligand, the mutated ROR1 fails to activate NF-κB. Ror1 is expressed in the inner ear during development at embryonic and postnatal stages. We demonstrate that Ror1 mutant mice are severely deaf, with preserved otoacoustic emissions. Anatomically, mutant mice display malformed cochleae. Axons of spiral ganglion neurons show fasciculation defects. Type I neurons show impaired synapses with inner hair cells, and type II neurons display aberrant projections through the cochlear sensory epithelium. We conclude that Ror1 is crucial for spiral ganglion neurons to innervate auditory hair cells. Impairment of ROR1 function largely affects development of the inner ear and hearing in humans and mice.

Describing sequencing results of structural chromosome rearrangements with a suggested next-generation cytogenetic nomenclature.

With recent rapid advances in genomic technologies, precise delineation of structural chromosome rearrangements at the nucleotide level is becoming increasingly feasible. In this era of "next-generation cytogenetics" (i.e., an integration of traditional cytogenetic techniques and next-generation sequencing), a consensus nomenclature is essential for accurate communication and data sharing. Currently, nomenclature for describing the sequencing data of these aberrations is lacking. Herein, we present a system called Next-Gen Cytogenetic Nomenclature, which is concordant with the International System for Human Cytogenetic Nomenclature (2013). This system starts with the alignment of rearrangement sequences by BLAT or BLAST (alignment tools) and arrives at a concise and detailed description of chromosomal changes. To facilitate usage and implementation of this nomenclature, we are developing a program designated BLA(S)T Output Sequence Tool of Nomenclature (BOSToN), a demonstrative version of which is accessible online. A standardized characterization of structural chromosomal rearrangements is essential both for research analyses and for application in the clinical setting.

Parallelized engineering of mutational models using piggyBac transposon delivery of CRISPR libraries.

New technologies and large-cohort studies have enabled novel variant discovery and association at unprecedented scale, yet functional characterization of these variants remains paramount to deciphering disease mechanisms. Approaches that facilitate parallelized genome editing of cells of interest or induced pluripotent stem cells (iPSCs) have become critical tools toward this goal. Here, we developed an approach that incorporates libraries of CRISPR-Cas9 guide RNAs (gRNAs) together with inducible Cas9 into a piggyBac (PB) transposon system to engineer dozens to hundreds of genomic variants in parallel against isogenic cellular backgrounds. This method empowers loss-of-function (LoF) studies through the introduction of insertions or deletions (indels) and copy-number variants (CNVs), though generating specific nucleotide changes is possible with prime editing. The ability to rapidly establish high-quality mutational models at scale will facilitate the development of isogenic cellular collections and catalyze comparative functional genomic studies investigating the roles of hundreds of genes and mutations in development and disease.

Transcriptome analysis in a humanized mouse model of familial dysautonomia reveals tissue-specific gene expression disruption in the peripheral nervous system.

Familial dysautonomia (FD) is a rare recessive neurodevelopmental disease caused by a splice mutation in the Elongator acetyltransferase complex subunit 1 (ELP1) gene. This mutation results in a tissue-specific reduction of ELP1 protein, with the lowest levels in the central and peripheral nervous systems (CNS and PNS, respectively). FD patients exhibit complex neurological phenotypes due to the loss of sensory and autonomic neurons. Disease symptoms include decreased pain and temperature perception, impaired or absent myotatic reflexes, proprioceptive ataxia, and progressive retinal degeneration. While the involvement of the PNS in FD pathogenesis has been clearly recognized, the underlying mechanisms responsible for the preferential neuronal loss remain unknown. In this study, we aimed to elucidate the molecular mechanisms underlying FD by conducting a comprehensive transcriptome analysis of neuronal tissues from the phenotypic mouse model TgFD9; Elp1Δ20/flox. This mouse recapitulates the same tissue-specific ELP1 mis-splicing observed in patients while modeling many of the disease manifestations. Comparison of FD and control transcriptomes from dorsal root ganglion (DRG), trigeminal ganglion (TG), medulla (MED), cortex, and spinal cord (SC) showed significantly more differentially expressed genes (DEGs) in the PNS than the CNS. We then identified genes that were tightly co-expressed and functionally dependent on the level of full-length ELP1 transcript. These genes, defined as ELP1 dose-responsive genes, were combined with the DEGs to generate tissue-specific dysregulated FD signature genes and networks. Within the PNS networks, we observed direct connections between Elp1 and genes involved in tRNA synthesis and genes related to amine metabolism and synaptic signaling. Importantly, transcriptomic dysregulation in PNS tissues exhibited enrichment for neuronal subtype markers associated with peptidergic nociceptors and myelinated sensory neurons, which are known to be affected in FD. In summary, this study has identified critical tissue-specific gene networks underlying the etiology of FD and provides new insights into the molecular basis of the disease.

Wastewater based surveillance can be used to reduce clinical testing intensity on a university campus.

Clinical testing has been a vital part of the response to and suppression of the COVID-19 pandemic; however, testing imposes significant burdens on a population. College students had to contend with clinical testing while simultaneously dealing with health risks and the academic pressures brought on by quarantines, changes to virtual platforms, and other disruptions to daily life. The objective of this study was to analyze whether wastewater surveillance can be used to decrease the intensity of clinical testing while maintaining reliable measurements of diseases incidence on campus. Twelve months of human health and wastewater surveillance data for eight residential buildings on a university campus were analyzed to establish how SARS-CoV-2 levels in the wastewater can be used to minimize clinical testing burden on students. Wastewater SARS-CoV-2 levels were used to create multiple scenarios, each with differing levels of testing intensity, which were compared to the actual testing volumes implemented by the university. We found that scenarios in which testing intensity fluctuations matched rise and falls in SARS-CoV-2 wastewater levels had stronger correlations between SARS-CoV-2 levels and recorded clinical positives. In addition to stronger correlations, most scenarios resulted in overall fewer weekly clinical tests performed. We suggest the use of wastewater surveillance to guide COVID-19 testing as it can significantly increase the efficacy of COVID-19 surveillance while reducing the burden placed on college students during a pandemic. Future efforts should be made to integrate wastewater surveillance into clinical testing strategies implemented on college campuses.

The conserved tetrameric subunit stoichiometry of Slc26 proteins.

The Slc26 family proteins, with one possible exception, transport anions across membranes in a wide variety of tissues in vertebrates, invertebrates, and plants. Mutations in human members of the family are a significant cause of disease. Slc26 family proteins are thought to be oligomers, but their stoichiometry of association is in dispute. A recent study, using sequential bleaching of single fluorophore-coupled molecules in membrane fragments, demonstrated that mammalian Slc26a5 (prestin) is a tetramer. In this article, the stoichiometry of two nonmammalian prestins and three human SLC26 proteins has been analyzed by the same method, including the evolutionarily-distant SLC26A11. The analysis showed that tetramerization is common and likely to be ubiquitous among Slc26 proteins, at least in vertebrates. The implication of the findings is that tetramerization is present for functional reasons.

Transcriptome analysis in a humanized mouse model of familial dysautonomia reveals tissue-specific gene expression disruption in the peripheral nervous system.

Familial dysautonomia (FD) is a rare recessive neurodevelopmental disease caused by a splice mutation in the Elongator acetyltransferase complex subunit 1 ( ELP1 ) gene. This mutation results in a tissue-specific reduction of ELP1 protein, with the lowest levels in the central and peripheral nervous systems (CNS and PNS, respectively). FD patients exhibit complex neurological phenotypes due to the loss of sensory and autonomic neurons. Disease symptoms include decreased pain and temperature perception, impaired or absent myotatic reflexes, proprioceptive ataxia, and progressive retinal degeneration. While the involvement of the PNS in FD pathogenesis has been clearly recognized, the underlying mechanisms responsible for the preferential neuronal loss remain unknown. In this study, we aimed to elucidate the molecular mechanisms underlying FD by conducting a comprehensive transcriptome analysis of neuronal tissues from the phenotypic mouse model TgFD9 ; Elp1 Δ 20/flox . This mouse recapitulates the same tissue-specific ELP1 mis-splicing observed in patients while modeling many of the disease manifestations. Comparison of FD and control transcriptomes from dorsal root ganglion (DRG), trigeminal ganglion (TG), medulla (MED), cortex, and spinal cord (SC) showed significantly more differentially expressed genes (DEGs) in the PNS than the CNS. We then identified genes that were tightly co-expressed and functionally dependent on the level of full-length ELP1 transcript. These genes, defined as ELP1 dose-responsive genes, were combined with the DEGs to generate tissue-specific dysregulated FD signature genes and networks. Within the PNS networks, we observed direct connections between Elp1 and genes involved in tRNA synthesis and genes related to amine metabolism and synaptic signaling. Importantly, transcriptomic dysregulation in PNS tissues exhibited enrichment for neuronal subtype markers associated with peptidergic nociceptors and myelinated sensory neurons, which are known to be affected in FD. In summary, this study has identified critical tissue-specific gene networks underlying the etiology of FD and provides new insights into the molecular basis of the disease.

Expanding a Wastewater-Based Surveillance Methodology for DNA Isolation from a Workflow Optimized for SARS-CoV-2 RNA Quantification.

Wastewater-based surveillance (WBS) is a noninvasive, epidemiological strategy for assessing the spread of COVID-19 in communities. This strategy was based upon wastewater RNA measurements of the viral target, severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). The utility of WBS for assessing the spread of COVID-19 has motivated research to measure targets beyond SARS-CoV-2, including pathogens containing DNA. The objective of this study was to establish the necessary steps for isolating DNA from wastewater by modifying a long-standing RNA-specific extraction workflow optimized for SARS-CoV-2 detection. Modifications were made to the sample concentration process and included an evaluation of bead bashing prior to the extraction of either DNA or RNA. Results showed that bead bashing reduced detection of RNA from wastewater but improved recovery of DNA as assessed by quantitative polymerase chain reaction (qPCR). Bead bashing is therefore not recommended for the quantification of RNA viruses using qPCR. Whereas for Mycobacterium bacterial DNA isolation, bead bashing was necessary for improving qPCR quantification. Overall, we recommend 2 separate workflows, one for RNA viruses that does not include bead bashing and one for other microbes that use bead bashing for DNA isolation. The experimentation done here shows that current-standing WBS program methodologies optimized for SARS-CoV-2 need to be modified and reoptimized to allow for alternative pathogens to be readily detected and monitored, expanding its utility as a tool for public health assessment.

Geospatially-resolved public-health surveillance via wastewater sequencing.

Wastewater, which contains everything from pathogens to pollutants, is a geospatially-and temporally-linked microbial fingerprint of a given population. As a result, it can be leveraged for monitoring multiple dimensions of public health across locales and time. Here, we integrate targeted and bulk RNA sequencing (n=1,419 samples) to track the viral, bacterial, and functional content over geospatially distinct areas within Miami Dade County from 2020-2022. First, we used targeted amplicon sequencing (n=966) to track diverse SARS-CoV-2 variants across space and time, and we found a tight correspondence with clinical caseloads from University students (N = 1,503) and Miami-Dade County hospital patients (N = 3,939 patients), as well as an 8-day earlier detection of the Delta variant in wastewater vs. in patients. Additionally, in 453 metatranscriptomic samples, we demonstrate that different wastewater sampling locations have clinically and public-health-relevant microbiota that vary as a function of the size of the human population they represent. Through assembly, alignment-based, and phylogenetic approaches, we also detect multiple clinically important viruses (e.g., norovirus ) and describe geospatial and temporal variation in microbial functional genes that indicate the presence of pollutants. Moreover, we found distinct profiles of antimicrobial resistance (AMR) genes and virulence factors across campus buildings, dorms, and hospitals, with hospital wastewater containing a significant increase in AMR abundance. Overall, this effort lays the groundwork for systematic characterization of wastewater to improve public health decision making and a broad platform to detect emerging pathogens.