RESUMEN
Transcription Factors (TFs) regulate gene expression by facilitating or disrupting the formation of transcription initiation machinery at particular genomic loci. Since TF occupancy is driven in part by recognition of DNA sequence, genetic variation can influence TF-DNA associations and gene regulation. To identify variants that impact TF binding in human brain tissues, we assessed allele specific binding (ASB) at heterozygous variants for 94 TFs in 9 brain regions from two donors. Leveraging graph genomes constructed from phased genomic sequence data, we compared ChIP-seq signals between alleles at heterozygous variants within each brain region and identified thousands of variants exhibiting ASB for at least one TF. ASB reproducibility was measured by comparisons between independent experiments both within and between donors. We found that rarer alleles in the general population more frequently led to reduced TF binding, whereas common variation had an equal likelihood of increasing or decreasing binding. Motif analysis revealed TF-specific effects, with ASB variants for certain TFs displaying a greater incidence of motif alterations, as well as enrichments for variants under purifying selection. Notably, neuron-specific cis-regulatory elements (cCREs) showed depletion for ASB variants. We identified 2,670 ASB variants with prior evidence of allele-specific gene expression in the brain from GTEx data and observed increasing eQTL effect direction concordance as ASB significance increases. These results provide a valuable and unique resource for mechanistic analysis of cis-regulatory variation in human brain tissue.
RESUMEN
An expanded CAG repeat in the huntingtin gene (HTT) causes Huntington's disease (HD). Since the length of uninterrupted CAG repeat, not polyglutamine, determines the age-at-onset in HD, base editing strategies to convert CAG to CAA are anticipated to delay onset by shortening the uninterrupted CAG repeat. Here, we developed base editing strategies to convert CAG in the repeat to CAA and determined their molecular outcomes and effects on relevant disease phenotypes. Base editing strategies employing combinations of cytosine base editors and guide RNAs (gRNAs) efficiently converted CAG to CAA at various sites in the CAG repeat without generating significant indels, off-target edits, or transcriptome alterations, demonstrating their feasibility and specificity. Candidate BE strategies converted CAG to CAA on both expanded and non-expanded CAG repeats without altering HTT mRNA and protein levels. In addition, somatic CAG repeat expansion, which is the major disease driver in HD, was significantly decreased in the liver by a candidate BE strategy treatment in HD knock-in mice carrying canonical CAG repeats. Notably, CAG repeat expansion was abolished entirely in HD knock-in mice carrying CAA-interrupted repeats, supporting the therapeutic potential of CAG-to-CAA conversion strategies in HD and potentially other repeat expansion disorders.
Asunto(s)
Edición Génica , Proteína Huntingtina , Enfermedad de Huntington , Expansión de Repetición de Trinucleótido , Enfermedad de Huntington/genética , Enfermedad de Huntington/terapia , Animales , Edición Génica/métodos , Ratones , Proteína Huntingtina/genética , Proteína Huntingtina/metabolismo , Expansión de Repetición de Trinucleótido/genética , Modelos Animales de Enfermedad , Humanos , Mutación , Técnicas de Sustitución del GenRESUMEN
Transcription factors (TFs) orchestrate gene expression programs crucial for brain function, but we lack detailed information about TF binding in human brain tissue. We generated a multiomic resource (ChIP-seq, ATAC-seq, RNA-seq, DNA methylation) on bulk tissues and sorted nuclei from several postmortem brain regions, including binding maps for more than 100 TFs. We demonstrate improved measurements of TF activity, including motif recognition and gene expression modeling, upon identification and removal of high TF occupancy regions. Further, predictive TF binding models demonstrate a bias for these high-occupancy sites. Neuronal TFs SATB2 and TBR1 bind unique regions depleted for such sites and promote neuronal gene expression. Binding sites for TFs, including TBR1 and PKNOX1, are enriched for risk variants associated with neuropsychiatric disorders, predominantly in neurons. This work, titled BrainTF, is a powerful resource for future studies seeking to understand the roles of specific TFs in regulating gene expression in the human brain.
Asunto(s)
Encéfalo , Factores de Transcripción , Humanos , Factores de Transcripción/metabolismo , Factores de Transcripción/genética , Encéfalo/metabolismo , Metilación de ADN , Neuronas/metabolismo , Sitios de Unión , Unión Proteica , Secuenciación de Inmunoprecipitación de CromatinaRESUMEN
Transcription Factors (TFs) influence gene expression by facilitating or disrupting the formation of transcription initiation machinery at particular genomic loci. Because genomic localization of TFs is in part driven by TF recognition of DNA sequence, variation in TF binding sites can disrupt TF-DNA associations and affect gene regulation. To identify variants that impact TF binding in human brain tissues, we quantified allele bias for 93 TFs analyzed with ChIP-seq experiments of multiple structural brain regions from two donors. Using graph genomes constructed from phased genomic sequence data, we compared ChIP-seq signal between alleles at heterozygous variants within each tissue sample from each donor. Comparison of results from different brain regions within donors and the same regions between donors provided measures of allele bias reproducibility. We identified thousands of DNA variants that show reproducible bias in ChIP-seq for at least one TF. We found that alleles that are rarer in the general population were more likely than common alleles to exhibit large biases, and more frequently led to reduced TF binding. Combining ChIP-seq with RNA-seq, we identified TF-allele interaction biases with RNA bias in a phased allele linked to 6,709 eQTL variants identified in GTEx data, 3,309 of which were found in neural contexts. Our results provide insights into the effects of both common and rare variation on gene regulation in the brain. These findings can facilitate mechanistic understanding of cis-regulatory variation associated with biological traits, including disease.
RESUMEN
An expanded CAG repeat in the huntingtin gene ( HTT ) causes Huntington's disease (HD). Since the length of uninterrupted CAG repeat, not polyglutamine, determines the age-at-onset in HD, base editing strategies to convert CAG to CAA are anticipated to delay onset by shortening the uninterrupted CAG repeat. Here, we developed base editing strategies to convert CAG in the repeat to CAA and determined their molecular outcomes and effects on relevant disease phenotypes. Base editing strategies employing combinations of cytosine base editors and gRNAs efficiently converted CAG to CAA at various sites in the CAG repeat without generating significant indels, off-target edits, or transcriptome alterations, demonstrating their feasibility and specificity. Candidate BE strategies converted CAG to CAA on both expanded and non-expanded CAG repeats without altering HTT mRNA and protein levels. In addition, somatic CAG repeat expansion, which is the major disease driver in HD, was significantly decreased by a candidate BE strategy treatment in HD knock-in mice carrying canonical CAG repeats. Notably, CAG repeat expansion was abolished entirely in HD knock-in mice carrying CAA-interrupted repeats, supporting the therapeutic potential of CAG-to-CAA conversion base editing strategies in HD and potentially other repeat expansion disorders.
RESUMEN
Cell type-specific transcriptional differences between brain tissues from donors with Alzheimer's disease (AD) and unaffected controls have been well documented, but few studies have rigorously interrogated the regulatory mechanisms responsible for these alterations. We performed single nucleus multiomics (snRNA-seq plus snATAC-seq) on 105,332 nuclei isolated from cortical tissues from 7 AD and 8 unaffected donors to identify candidate cis-regulatory elements (CREs) involved in AD-associated transcriptional changes. We detected 319,861 significant correlations, or links, between gene expression and cell type-specific transposase accessible regions enriched for active CREs. Among these, 40,831 were unique to AD tissues. Validation experiments confirmed the activity of many regions, including several candidate regulators of APP expression. We identified ZEB1 and MAFB as candidate transcription factors playing important roles in AD-specific gene regulation in neurons and microglia, respectively. Microglia links were globally enriched for heritability of AD risk and previously identified active regulatory regions.
RESUMEN
Recent genome-wide association studies of age-at-onset in Huntington's disease (HD) point to distinct modes of potential disease modification: altering the rate of somatic expansion of the HTT CAG repeat or altering the resulting CAG threshold length-triggered toxicity process. Here, we evaluated the mouse orthologs of two HD age-at-onset modifier genes, FAN1 and RRM2B, for an influence on somatic instability of the expanded CAG repeat in Htt CAG knock-in mice. Fan1 knock-out increased somatic expansion of Htt CAG repeats, in the juvenile- and the adult-onset HD ranges, whereas knock-out of Rrm2b did not greatly alter somatic Htt CAG repeat instability. Simultaneous knock-out of Mlh1, the ortholog of a third HD age-at-onset modifier gene (MLH1), which suppresses somatic expansion of the Htt knock-in CAG repeat, blocked the Fan1 knock-out-induced acceleration of somatic CAG expansion. This genetic interaction indicates that functional MLH1 is required for the CAG repeat destabilizing effect of FAN1 loss. Thus, in HD, it is uncertain whether the RRM2B modifier effect on timing of onset may be due to a DNA instability mechanism. In contrast, the FAN1 modifier effects reveal that functional FAN1 acts to suppress somatic CAG repeat expansion, likely in genetic interaction with other DNA instability modifiers whose combined effects can hasten or delay onset and other CAG repeat length-driven phenotypes.
Asunto(s)
Proteínas de Ciclo Celular/genética , Endodesoxirribonucleasas/genética , Exodesoxirribonucleasas/genética , Proteína Huntingtina/genética , Enfermedad de Huntington/genética , Enzimas Multifuncionales/genética , Homólogo 1 de la Proteína MutL/genética , Ribonucleótido Reductasas/genética , Edad de Inicio , Animales , Modelos Animales de Enfermedad , Genes Modificadores/genética , Predisposición Genética a la Enfermedad , Estudio de Asociación del Genoma Completo , Humanos , Enfermedad de Huntington/patología , Ratones , Ratones Noqueados , Fenotipo , Expansión de Repetición de Trinucleótido/genéticaRESUMEN
The differentiation status of tumors is used as a prognostic indicator, with tumors comprised of less differentiated cells exhibiting higher levels of aggressiveness that correlate with a poor prognosis. Although oncogenes contribute to blocking differentiation, it is not clear how they globally alter miRNA expression during differentiation to achieve this result. The pediatric sarcoma Alveolar Rhabdomyosarcoma, which is primarily characterized by the expression of the PAX3-FOXO1 oncogenic fusion protein, consists of undifferentiated muscle cells. However, it is unclear what role PAX3-FOXO1 plays in promoting the undifferentiated state. We demonstrate that expression of PAX3-FOXO1 globally alters the expression of over 80 individual miRNA during early myogenic differentiation, resulting in three primary effects: 1) inhibition of the expression of 51 miRNA essential for promoting myogenesis, 2) promoting the aberrant expression of 43 miRNA not normally expressed during myogenesis, and 3) altering the expression pattern of 39 additional miRNA. Combined, these changes are predicted to have an overall negative effect on myogenic differentiation. This is one of the first studies describing how an oncogene globally alters miRNA expression to block differentiation and has clinical implications for the development of much needed multi-faceted tumor-specific therapeutic regimens.
RESUMEN
While many solid tumors are defined by the presence of a particular oncogene, the role that this oncogene plays in driving transformation through the acquisition of aneuploidy and overcoming growth arrest are often not known. Further, although aneuploidy is present in many solid tumors, it is not clear whether it is the cause or effect of malignant transformation. The childhood sarcoma, Alveolar Rhabdomyosarcoma (ARMS), is primarily defined by the t(2;13)(q35;q14) translocation, creating the PAX3-FOXO1 fusion protein. It is unclear what role PAX3-FOXO1 plays in the initial stages of tumor development through the acquisition and persistence of aneuploidy. In this study we demonstrate that PAX3-FOXO1 serves as a driver mutation to initiate a cascade of mRNA and miRNA changes that ultimately reprogram proliferating myoblasts to induce the formation of ARMS. We present evidence that cells containing PAX3-FOXO1 have changes in the expression of mRNA and miRNA essential for maintaining proper chromosome number and structure thereby promoting aneuploidy. Further, we demonstrate that the presence of PAX3-FOXO1 alters the expression of growth factor related mRNA and miRNA, thereby overriding aneuploid-dependent growth arrest. Finally, we present evidence that phosphorylation of PAX3-FOXO1 contributes to these changes. This is one of the first studies describing how an oncogene and post-translational modifications drive the development of a tumor through the acquisition and persistence of aneuploidy. This mechanism has implications for other solid tumors where large-scale genomics studies may elucidate how global alterations contribute to tumor phenotypes allowing the development of much needed multi-faceted tumor-specific therapeutic regimens.
Asunto(s)
Proteína Forkhead Box O1/metabolismo , Mutación , Proteínas de Fusión Oncogénica/genética , Factor de Transcripción PAX3/metabolismo , Rabdomiosarcoma Alveolar/genética , Aneuploidia , Animales , Ciclo Celular , Proliferación Celular , Transformación Celular Neoplásica/genética , Aberraciones Cromosómicas , Progresión de la Enfermedad , Proteína Forkhead Box O1/genética , Perfilación de la Expresión Génica , Humanos , Ratones , Ratones Endogámicos C57BL , MicroARNs/metabolismo , Mitosis , Desarrollo de Músculos , Mioblastos/metabolismo , Factor de Transcripción PAX3/genética , Fosforilación , Procesamiento Proteico-Postraduccional , ARN Mensajero/metabolismo , Rabdomiosarcoma Alveolar/metabolismo , Translocación GenéticaAsunto(s)
Melanoma/metabolismo , Factores de Transcripción Paired Box/metabolismo , Línea Celular Tumoral , Proliferación Celular , Transformación Celular Neoplásica , Humanos , Melanocitos/citología , Melanoma/patología , Invasividad Neoplásica , Metástasis de la Neoplasia , Factor de Transcripción PAX3 , Fenotipo , Fosforilación , Serina/químicaRESUMEN
The myogenic transcription factor Pax3, a member of the paired class homeodomain family of transcription factors, plays an essential role in early skeletal muscle development. We previously demonstrated that Pax3 is phosphorylated at three specific residues (Ser201, Ser205, and Ser209) and that the pattern of phosphorylation at these sites changes throughout early myogenesis. Further, we demonstrated that the protein kinase CK2 phosphorylates Pax3 at Ser205 and that this phosphorylation event is required for the subsequent phosphorylation of Ser201 by GSK3ß. However, the kinase that phosphorylates Pax3 at Ser209 has yet to be identified. In the present work we use standard purification methods and in vitro biochemical analyses to provide solid evidence identifying the protein kinase CK2 as phosphorylating Pax3 at Ser209. Further, we qualitatively demonstrate that the phosphorylation of Pax3 at Ser209 by CK2 is enhanced when Ser205 is previously phosphorylated. Taken together, our results allow us to propose a mechanism to describe the ordered phosphorylation of Pax3 throughout early myogenesis.
Asunto(s)
Quinasa de la Caseína II/metabolismo , Diferenciación Celular , Desarrollo de Músculos , Músculo Esquelético/crecimiento & desarrollo , Mioblastos Esqueléticos/citología , Factores de Transcripción Paired Box/metabolismo , Animales , Células Cultivadas , Ratones , Modelos Biológicos , Músculo Esquelético/metabolismo , Mutación , Mioblastos Esqueléticos/metabolismo , Factor de Transcripción PAX3 , Factores de Transcripción Paired Box/genética , Fosforilación , Serina/metabolismoRESUMEN
Pax3, a member of the paired class homeodomain family of transcription factors, is essential for early skeletal muscle development and is key in the development of the childhood solid muscle tumor alveolar rhabdomyosarcoma (ARMS). ARMS is primarily characterized by a t(2;13)(q35;q14) chromosomal translocation, which fuses the 5'-coding sequences of Pax3 with the 3'-coding sequence of the forkhead transcription factor FOXO1 generating the oncogenic fusion protein Pax3-FOXO1. We previously demonstrated that Pax3 and Pax3-FOXO1 are phosphorylated by the protein kinase CK2 at serine 205 in proliferating primary myoblasts and that this phosphorylation event is rapidly lost from Pax3, but not Pax3-FOXO1 upon the induction of differentiation. However, reports suggested that additional sites of phosphorylation might be present on Pax3. In this report we use in vitro and in vivo analyses to identify serines 201 and 209 as additional sites of phosphorylation and along with serine 205 are the only sites of phosphorylation on Pax3. We provide solid evidence supporting the role of the protein kinase GSK3ß as phosphorylating Pax3 at serine 201. Using phospho-specific antibodies we demonstrate a changing pattern of phosphorylation at serines 201, 205, and 209 throughout early myogenic differentiation and that this pattern of phosphorylation is different for Pax3-FOXO1 in primary myoblasts and in several ARMS cell lines. Taken together, our results allow us to propose a molecular model to describe the changing pattern of phosphorylation for Pax3 and the altered phosphorylation for Pax3-FOXO1 during early myogenic differentiation.