DYT-THAP1: exploring gene expression in fibroblasts for potential biomarker discovery

Dystonia due to pathogenic variants in the THAP1 gene (DYT-THAP1) shows variable expressivity and reduced penetrance of ~ 50%. Since THAP1 encodes a transcription factor, modifiers influencing this variability likely operate at the gene expression level. This study aimed to assess the transferability of differentially expressed genes (DEGs) in neuronal cells related to pathogenic variants in the THAP1 gene, which were previously identified by transcriptome analyses. For this, we performed quantitative (qPCR) and Digital PCR (dPCR) in cultured fibroblasts. RNA was extracted from THAP1 manifesting (MMCs) and non-manifesting mutation carriers (NMCs) as well as from healthy controls. The expression profiles of ten of 14 known neuronal DEGs demonstrated differences in fibroblasts between these three groups. This included transcription factors and targets (ATF4, CLN3, EIF2A, RRM1, YY1), genes involved in G protein-coupled receptor signaling (BDKRB2, LPAR1), and a gene linked to apoptosis and DNA replication/repair (CRADD), which all showed higher expression levels in MMCs and NMCs than in controls. Moreover, the analysis of genes linked to neurological disorders (STXBP1, TOR1A) unveiled differences in expression patterns between MMCs and controls. Notably, the genes CUEDC2, DRD4, ECH1, and SIX2 were not statistically significantly differentially expressed in fibroblast cultures. With > 70% of the tested genes being DEGs also in fibroblasts, fibroblasts seem to be a suitable model for DYT-THAP1 research despite some restrictions. Furthermore, at least some of these DEGs may potentially also serve as biomarkers of DYT-THAP1 and influence its penetrance and expressivity.

Genetics and Pathogenesis of Dystonia

Dystonia is a clinically and genetically highly heterogeneous neurological disorder characterized by abnormal movements and postures caused by involuntary sustained or intermittent muscle contractions. A number of groundbreaking genetic and molecular insights have recently been gained. While they enable genetic testing and counseling, their translation into new therapies is still limited. However, we are beginning to understand shared pathophysiological pathways and molecular mechanisms. It has become clear that dystonia results from a dysfunctional network involving the basal ganglia, cerebellum, thalamus, and cortex. On the molecular level, more than a handful of, often intertwined, pathways have been linked to pathogenic variants in dystonia genes, including gene transcription during neurodevelopment (e.g., KMT2B, THAP1), calcium homeostasis (e.g., ANO3, HPCA), striatal dopamine signaling (e.g., GNAL), endoplasmic reticulum stress response (e.g., EIF2AK2, PRKRA, TOR1A), autophagy (e.g., VPS16), and others. Thus, different forms of dystonia can be molecularly grouped, which may facilitate treatment development in the future.

Transcriptional regulatory network for neuron-glia interactions and its implication for DYT6 dystonia

Advances in sequencing technologies have identified novel genes associated with inherited forms of dystonia, providing valuable insights into its genetic basis and revealing diverse genetic pathways and mechanisms involved in its pathophysiology. Since identifying genetic variation in the transcription factor coding THAP1 gene linked to isolated dystonia, numerous investigations have employed transcriptomic studies in DYT-THAP1 models to uncover pathogenic molecular mechanisms underlying dystonia. This review examines key findings from transcriptomic studies conducted on in vivo and in vitro DYT-THAP1 models, which demonstrate that the THAP1-regulated transcriptome is diverse and cell-specific, yet it is bound and co-regulated by a common set of proteins. Prominent among its functions, THAP1 and its co-regulatory network target molecular pathways critical for generating myelinating oligodendrocytes that ensheath axons and generate white matter in the central nervous system. Several lines of investigation have demonstrated the importance of myelination and oligodendrogenesis in motor function during development and in adults, emphasizing the non-cell autonomous contributions of glial cells to neural circuits involved in motor function. Further research on the role of myelin abnormalities in motor deficits in DYT6 models will enhance our understanding of axon-glia interactions in dystonia pathophysiology and provide potential therapeutic interventions targeting these pathways.

Emerging role of a systems biology approach to elucidate factors of reduced penetrance: transcriptional changes in THAP1-linked dystonia as an example

Pathogenic variants in THAP1 can cause dystonia with a penetrance of about 50 %. The underlying mechanisms are unknown and can be considered as means of endogenous disease protection. Since THAP1 encodes a transcription factor, drivers of this variability putatively act at the transcriptome level. Several transcriptome studies tried to elucidate THAP1 function in diverse cellular and mouse models, including mutation carrier-derived cells and iPSC-derived neurons, unveiling various differentially expressed genes and affected pathways. These include nervous system development, dopamine signalling, myelination, or cell-cell adhesion. A network diffusion analysis revealed mRNA splicing, mitochondria, DNA repair, and metabolism as significant pathways that may represent potential targets for therapeutic interventions.

A pathogenic DYT-THAP1 dystonia mutation causes hypomyelination and loss of YY1 binding

Dystonia is a disabling disease that manifests as prolonged involuntary twisting movements. DYT-THAP1 is an inherited form of isolated dystonia caused by mutations in THAP1 encoding the transcription factor THAP1. The phe81leu (F81L) missense mutation is representative of a category of poorly understood mutations that do not occur on residues critical for DNA binding. Here, we demonstrate that the F81L mutation (THAP1F81L) impairs THAP1 transcriptional activity and disrupts CNS myelination. Strikingly, THAP1F81L exhibits normal DNA binding but causes a significantly reduced DNA binding of YY1, its transcriptional partner that also has an established role in oligodendrocyte lineage progression. Our results suggest a model of molecular pathogenesis whereby THAP1F81L normally binds DNA but is unable to efficiently organize an active transcription complex.

Ying Yang 1 engagement in brain pathology

Herein, we discuss data concerning the involvement of transcription factor Yin Yang 1 (YY1) in the development of brain diseases, highlighting mechanisms of its pathological actions. YY1 plays an important role in the developmental and adult pathology of the nervous system. YY1 is essential for neurulation as well as maintenance and differentiation of neuronal progenitor cells and oligodendrocytes regulating both neural and glial tissues of the brain. Lack of a YY1 gene causes many developmental abnormalities and anatomical malformations of the central nervous system (CNS). Once dysregulated, YY1 exerts multiple neuropathological actions being involved in the induction of many brain disorders like stroke, epilepsy, Alzheimer's and Parkinson's diseases, autism spectrum disorder, dystonia, and brain tumors. Better understanding of YY1's dysfunction in the nervous system may lead to the development of novel therapeutic strategies related to YY1's actions.

The dystonia gene THAP1 controls DNA double-strand break repair choice

The Shieldin complex shields double-strand DNA breaks (DSBs) from nucleolytic resection. Curiously, the penultimate Shieldin component, SHLD1, is one of the least abundant mammalian proteins. Here, we report that the transcription factors THAP1, YY1, and HCF1 bind directly to the SHLD1 promoter, where they cooperatively maintain the low basal expression of SHLD1, thereby ensuring a proper balance between end protection and resection during DSB repair. The loss of THAP1-dependent SHLD1 expression confers cross-resistance to poly (ADP-ribose) polymerase (PARP) inhibitor and cisplatin in BRCA1-deficient cells and shorter progression-free survival in ovarian cancer patients. Moreover, the embryonic lethality and PARPi sensitivity of BRCA1-deficient mice is rescued by ablation of SHLD1. Our study uncovers a transcriptional network that directly controls DSB repair choice and suggests a potential link between DNA damage and pathogenic THAP1 mutations, found in patients with the neurodevelopmental movement disorder adult-onset torsion dystonia type 6.

Isolated dystonia: clinical and genetic updates

Four genes associated with isolated dystonia are currently well replicated and validated. DYT-THAP1 manifests as young-onset generalized dystonia with predominant craniocervical symptoms; and is associated with mostly deleterious missense variation in the THAP1 gene. De novo and inherited missense and protein truncating variation in GNAL as well as primarily missense variation in ANO3 cause isolated focal and/or segmental dystonia with preference for the upper half of the body and older ages at onset. The GAG deletion in TOR1A is associated with generalized dystonia with onset in childhood in the lower limbs. Rare variation in these genes causes monogenic sporadic and inherited forms of isolated dystonia; common variation may confer risk and imply that dystonia is a polygenic trait in a subset of cases. Although candidate gene screens have been successful in the past in detecting gene-disease associations, recent application of whole-genome and whole-exome sequencing methods enable unbiased capture of all genetic variation that may explain the phenotype. However, careful variant-level evaluation is necessary in every case, even in genes that have previously been associated with disease. We review the genetic architecture and phenotype of DYT-THAP1, DYT-GNAL, DYT-ANO3, and DYT-TOR1A by collecting case reports from the literature and performing variant classification using pathogenicity criteria.

HCF-1 promotes cell cycle progression by regulating the expression of CDC42

The eukaryotic cell cycle involves a highly orchestrated series of events in which the cellular genome is replicated during a synthesis (S) phase and each of the two resulting copies are segregated properly during mitosis (M). Host cell factor-1 (HCF-1) is a transcriptional co-regulator that is essential for and has been implicated in basic cellular processes, such as transcriptional regulation and cell cycle progression. Although a series of HCF-1 transcriptional targets have been identified, few functional clues have been provided, especially for chromosome segregation. Our results showed that HCF-1 activated CDC42 expression by binding to the −881 to −575 region upstream of the CDC42 transcription start site, and the regulation of CDC42 expression by HCF-1 was correlated with cell cycle progression. The overexpression of a spontaneously cycling and constitutively active CDC42 mutant (CDC42F28L) rescued G1 phase delay and multinucleate defects in mitosis upon the loss of HCF-1. Therefore, these results establish that HCF-1 ensures proper cell cycle progression by regulating the expression of CDC42, which indicates a possible mechanism of cell cycle coordination and the regulation mode of typical Rho GTPases.

DIVERSITY in binding, regulation, and evolution revealed from high-throughput ChIP

Genome-wide in vivo protein-DNA interactions are routinely mapped using high-throughput chromatin immunoprecipitation (ChIP). ChIP-reported regions are typically investigated for enriched sequence-motifs, which are likely to model the DNA-binding specificity of the profiled protein and/or of co-occurring proteins. However, simple enrichment analyses can miss insights into the binding-activity of the protein. Note that ChIP reports regions making direct contact with the protein as well as those binding through intermediaries. For example, consider a ChIP experiment targeting protein X, which binds DNA at its cognate sites, but simultaneously interacts with four other proteins. Each of these proteins also binds to its own specific cognate sites along distant parts of the genome, a scenario consistent with the current view of transcriptional hubs and chromatin loops. Since ChIP will pull down all X-associated regions, the final reported data will be a union of five distinct sets of regions, each containing binding sites of one of the five proteins, respectively. Characterizing all five different motifs and the corresponding sets is important to interpret the ChIP experiment and ultimately, the role of X in regulation. We present diversity which attempts exactly this: it partitions the data so that each partition can be characterized with its own de novo motif. Diversity uses a Bayesian approach to identify the optimal number of motifs and the associated partitions, which together explain the entire dataset. This is in contrast to standard motif finders, which report motifs individually enriched in the data, but do not necessarily explain all reported regions. We show that the different motifs and associated regions identified by diversity give insights into the various complexes that may be forming along the chromatin, something that has so far not been attempted from ChIP data. Webserver at http://diversity.ncl.res.in/; standalone (Mac OS X/Linux) from https://github.com/NarlikarLab/DIVERSITY/releases/tag/v1.0.0.

A high-throughput chromatin immunoprecipitation (ChIP) experiment identifies genomic regions bound by a protein in vivo. Current motif-discovery approaches seek an enriched motif signature in the reported regions, which they can attribute to the protein’s binding preferences. However, Diversity models the fact that since a ChIP experiment pulls down regions participating in all complexes involving the profiled protein, the reported regions are in all likelihood, a collection of different types of protein-DNA contacts. Diversity asks a different question: what sequence component caused a specific region to be reported in a ChIP experiment? The answer, in combination with additional data such as sequence conservation, SNPs, chromatin structure, downstream gene-expression, etc. can yield insights into the diverse regulatory mechanisms at play. The added benefits of a webserver and a standalone parallel version make diversity a practical tool for discovering new biology from ChIP experiments.

Mutations in THAP1/DYT6 reveal that diverse dystonia genes disrupt similar neuronal pathways and functions

Dystonia is characterized by involuntary muscle contractions. Its many forms are genetically, phenotypically and etiologically diverse and it is unknown whether their pathogenesis converges on shared pathways. Mutations in THAP1 [THAP (Thanatos-associated protein) domain containing, apoptosis associated protein 1], a ubiquitously expressed transcription factor with DNA binding and protein-interaction domains, cause dystonia, DYT6. There is a unique, neuronal 50-kDa Thap1-like immunoreactive species, and Thap1 levels are auto-regulated on the mRNA level. However, THAP1 downstream targets in neurons, and the mechanism via which it causes dystonia are largely unknown. We used RNA-Seq to assay the in vivo effect of a heterozygote Thap1 C54Y or ΔExon2 allele on the gene transcription signatures in neonatal mouse striatum and cerebellum. Enriched pathways and gene ontology terms include eIF2α Signaling, Mitochondrial Dysfunction, Neuron Projection Development, Axonal Guidance Signaling, and Synaptic LongTerm Depression, which are dysregulated in a genotype and tissue-dependent manner. Electrophysiological and neurite outgrowth assays were consistent with those enrichments, and the plasticity defects were partially corrected by salubrinal. Notably, several of these pathways were recently implicated in other forms of inherited dystonia, including DYT1. We conclude that dysfunction of these pathways may represent a point of convergence in the pathophysiology of several forms of inherited dystonia.

Dystonia is a brain disorder that causes disabling involuntary muscle contractions and abnormal postures. Mutations in THAP1, a zinc-finger transcription factor, cause DYT6, but its neuronal targets and functions are unknown. In this study, we sought to determine the effects of Thap1C54Y and ΔExon2 alleles on the gene transcription signatures at postnatal day 1 (P1) in the mouse striatum and cerebellum in order to correlate function with specific genes or pathways. Our unbiased transcriptomics approach showed that Thap1 mutants revealed multiple signaling pathways involved in neuronal plasticity, axonal guidance, and oxidative stress response, which are also present in other forms of dystonia, particularly DYT1. We conclude that dysfunction of these pathways may represent a point of convergence on the pathogenesis of unrelated forms of inherited dystonia.