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.

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.

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.

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.

Correlative analysis of wastewater trends with clinical cases and hospitalizations through five dominant variant waves of COVID-19.

Wastewater-based epidemiology (WBE) has been utilized to track community infections of Coronavirus Disease 2019 (COVID-19) by detecting RNA of the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), within samples collected from wastewater. The correlations between community infections and wastewater measurements of the RNA can potentially change as SARS-CoV-2 evolves into new variations by mutating. This study analyzed SARS-CoV-2 RNA, and indicators of human waste in wastewater from two sewersheds of different scales (University of Miami (UM) campus and Miami-Dade County Central District wastewater treatment plant (CDWWTP)) during five internally defined COVID-19 variant dominant periods (Initial, Pre-Delta, Delta, Omicron and Post-Omicron wave). SARS-CoV-2 RNA quantities were compared against COVID-19 clinical cases and hospitalizations to evaluate correlations with wastewater SARS-CoV-2 RNA. Although correlations between documented clinical cases and hospitalizations were high, prevalence for a given wastewater SARS-CoV-2 level varied depending upon the variant analyzed. The correlative relationship was significantly steeper (more cases per level found in wastewater) for the Omicron-dominated period. For hospitalization, the relationships were steepest for the Initial wave, followed by the Delta wave with flatter slopes during all other waves. Overall results were interpreted in the context of SARS-CoV-2 virulence and vaccination rates among the community.

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.

Relationships between SARS-CoV-2 in Wastewater and COVID-19 Clinical Cases and Hospitalizations, with and without Normalization against Indicators of Human Waste.

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) in wastewater has been used to track community infections of coronavirus disease-2019 (COVID-19), providing critical information for public health interventions. Since levels in wastewater are dependent upon human inputs, we hypothesize that tracking infections can be improved by normalizing wastewater concentrations against indicators of human waste [Pepper Mild Mottle Virus (PMMoV), ß-2 Microglobulin (B2M), and fecal coliform]. In this study, we analyzed SARS-CoV-2 and indicators of human waste in wastewater from two sewersheds of different scales: a University campus and a wastewater treatment plant. Wastewater data were combined with complementary COVID-19 case tracking to evaluate the efficiency of wastewater surveillance for forecasting new COVID-19 cases and, for the larger scale, hospitalizations. Results show that the normalization of SARS-CoV-2 levels by PMMoV and B2M resulted in improved correlations with COVID-19 cases for campus data using volcano second generation (V2G)-qPCR chemistry (r s = 0.69 without normalization, r s = 0.73 with normalization). Mixed results were obtained for normalization by PMMoV for samples collected at the community scale. Overall benefits from normalizing with measures of human waste depend upon qPCR chemistry and improves with smaller sewershed scale. We recommend further studies that evaluate the efficacy of additional normalization targets.

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.

Comparison of Electronegative Filtration to Magnetic Bead-Based Concentration and V2G-qPCR to RT-qPCR for Quantifying Viral SARS-CoV-2 RNA from Wastewater.

Methods of wastewater concentration (electronegative filtration (ENF) versus magnetic bead-based concentration (MBC)) were compared for the analysis of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), beta-2 microglobulin, and human-coronavirus OC43. Using ENF as the concentration method, two quantitative Polymerase Chain Reaction (qPCR) analytical methods were also compared: Volcano 2nd Generation (V2G)-qPCR and reverse transcriptase (RT)-qPCR measuring three different targets of the virus responsible for the COVID-19 illness (N1, modified N3, and ORF1ab). Correlations between concentration methods were strong and statistically significant for SARS-CoV-2 (r=0.77, p<0.001) and B2M (r=0.77, p<0.001). Comparison of qPCR analytical methods indicate that, on average, each method provided equivalent results with average ratios of 0.96, 0.96 and 1.02 for N3 to N1, N3 to ORF1ab, and N1 to ORF1ab and were supported by significant (p<0.001) correlation coefficients (r =0.67 for V2G (N3) to RT (N1), r =0.74 for V2G (N3) to RT (ORF1ab), r = 0.81 for RT (N1) to RT (ORF1ab)). Overall results suggest that the two concentration methods and qPCR methods provide equivalent results, although variability is observed for individual measurements. Given the equivalency of results, additional advantages and disadvantages, as described in the discussion, are to be considered when choosing an appropriate method.

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.

A Rapid, Isothermal, and Point-of-Care System for COVID-19 Diagnostics.

The COVID-19 pandemic has had a profound, detrimental effect on economies and societies worldwide. Where the pandemic has been controlled, extremely high rates of diagnostic testing for the SARS-CoV-2 virus have proven critical, enabling isolation of cases and contact tracing. Recently, diagnostic testing has been supplemented with wastewater measures to evaluate the degree to which communities have infections. Whereas much testing has been done through traditional, centralized, clinical, or environmental laboratory methods, point-of-care testing has proven successful in reducing time to result. As the pandemic progresses and becomes more broadly distributed, further decentralization of diagnostic testing will be helpful to mitigate its spread. This will be particularly both challenging and critical in settings with limited resources due to lack of medical infrastructure and expertise as well as requirements to return results quickly. In this article, we validate the tiny isothermal nucleic acid quantification system (TINY) and a novel loop-mediated isothermal amplification (LAMP)-based assay for the point-of-care diagnosis of SARS-CoV-2 infection in humans and also for in-the-field, point-of-collection surveillance of wastewater. The TINY system is portable and designed for use in settings with limited resources. It can be powered by electrical, solar, or thermal energy and is robust against interruptions in services. These applied testing examples demonstrate that this novel detection platform is a simpler procedure than reverse-transcription quantitative polymerase chain reaction, and moreover, this TINY instrument and LAMP assay combination has the potential to effectively provide both point-of-care diagnosis of individuals and point-of-collection environmental surveillance using wastewater.

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.

Next Generation Sequencing of Prenatal Structural Chromosomal Rearrangements Using Large-Insert Libraries.

Precise tests for genomic structural variation (SV) are essential for accurate diagnosis of prenatal genome abnormalities. The two most ubiquitous traditional methods for prenatal SV assessment, karyotyping and chromosomal microarrays, do not provide sufficient resolution for some clinically actionable SVs. Standard whole-genome sequencing (WGS) overcomes shortcomings of traditional techniques by providing base-pair resolution of the entire accessible genome. However, while sequencing costs have continued to decline in recent years, conventional WGS costs remain high for most routine clinical applications. Here, we describe a specialized WGS technique using large inserts (liWGS; also known as "jumping libraries") to resolve large (>5000-10,000 nucleotides) SVs at kilobase-resolution in prenatal samples, and at a fraction of the cost of standard WGS. We explicate the protocols for generating liWGS libraries and supplement with an overview for processing and analyzing liWGS data.