Defective Integrator activity shapes the transcriptome of patients with multiple sclerosis.

HP1α/CBX5 is an epigenetic regulator with a suspected role in multiple sclerosis (MS). Here, using high-depth RNA sequencing on monocytes, we identified a subset of MS patients with reduced CBX5 expression, correlating with progressive stages of the disease and extensive transcriptomic alterations. Examination of rare non-coding RNA species in these patients revealed impaired maturation/degradation of U snRNAs and enhancer RNAs, indicative of reduced activity of the Integrator, a complex with suspected links to increased MS risk. At protein-coding genes, compromised Integrator activity manifested in reduced pre-mRNA splicing efficiency and altered expression of genes regulated by RNA polymerase II pause-release. Inactivation of Cbx5 in the mouse mirrored most of these transcriptional defects and resulted in hypersensitivity to experimental autoimmune encephalomyelitis. Collectively, our observations suggested a major contribution of the Integrator complex in safeguarding against transcriptional anomalies characteristic of MS, with HP1α/CBX5 emerging as an unexpected regulator of this complex's activity. These findings bring novel insights into the transcriptional aspects of MS and provide potential new criteria for patient stratification.

Design and evaluation of antisense sequence length for modified mouse U7 small nuclear RNA to induce efficient pre-messenger RNA splicing modulation in vitro.

Pre-messenger RNA (pre-mRNA) splicing modulation is an attractive approach for investigating the mechanisms of genetic disorders caused by mis-splicing. Previous reports have indicated that a modified U7 small nuclear RNA (U7 snRNA) is a prospective tool for modulating splicing both in vitro and in vivo. To date, very few studies have investigated the role of antisense sequence length in modified U7 snRNA. In this study, we designed a series of antisense sequences with various lengths and evaluated their efficiency in inducing splicing modulation. To express modified U7 snRNAs, we constructed a series of plasmid DNA sequences which codes cytomegalovirus (CMV) enhancer, human U1 promoter, and modified mouse U7 snRNAs with antisense sequences of different lengths. We evaluated in vitro splicing modulation efficiency using a luciferase reporter system for simple and precise evaluation as well as reverse transcription-polymerase chain reaction to monitor splicing patterns. Our in vitro assay findings suggest that antisense sequences of modified mouse U7 snRNAs have an optimal length for efficient splicing modulation, which depends on the target exon. In addition, antisense sequences that were either too long or too short decreased splicing modulation efficiency. To confirm reproducibility, we performed an in vitro assay using two target genes, mouse Fas and mouse Dmd. Together, our data suggests that the antisense sequence length should be optimized for modified mouse U7 snRNAs to induce efficient splicing modulation.

Therapy-induced secretion of spliceosomal components mediates pro-survival crosstalk between ovarian cancer cells.

Ovarian cancer often develops resistance to conventional therapies, hampering their effectiveness. Here, using ex vivo paired ovarian cancer ascites obtained before and after chemotherapy and in vitro therapy-induced secretomes, we show that molecules secreted by ovarian cancer cells upon therapy promote cisplatin resistance and enhance DNA damage repair in recipient cancer cells. Even a short-term incubation of chemonaive ovarian cancer cells with therapy-induced secretomes induces changes resembling those that are observed in chemoresistant patient-derived tumor cells after long-term therapy. Using integrative omics techniques, we find that both ex vivo and in vitro therapy-induced secretomes are enriched with spliceosomal components, which relocalize from the nucleus to the cytoplasm and subsequently into the extracellular vesicles upon treatment. We demonstrate that these molecules substantially contribute to the phenotypic effects of therapy-induced secretomes. Thus, SNU13 and SYNCRIP spliceosomal proteins promote therapy resistance, while the exogenous U12 and U6atac snRNAs stimulate tumor growth. These findings demonstrate the significance of spliceosomal network perturbation during therapy and further highlight that extracellular signaling might be a key factor contributing to the emergence of ovarian cancer therapy resistance.

Functional analysis of the zinc finger modules of the Saccharomyces cerevisiae splicing factor Luc7.

Identification of splice sites is a critical step in pre-messenger RNA (pre-mRNA) splicing because the definition of the exon/intron boundaries controls what nucleotides are incorporated into mature mRNAs. The intron boundary with the upstream exon is initially identified through interactions with the U1 small nuclear ribonucleoprotein (snRNP). This involves both base-pairing between the U1 snRNA and the pre-mRNA as well as snRNP proteins interacting with the 5' splice site (5'ss)/snRNA duplex. In yeast, this duplex is buttressed by two conserved protein factors, Yhc1 and Luc7. Luc7 has three human paralogs (LUC7L, LUC7L2, and LUC7L3), which play roles in alternative splicing. What domains of these paralogs promote splicing at particular sites is not yet clear. Here, we humanized the zinc finger (ZnF) domains of the yeast Luc7 protein in order to understand their roles in splice site selection using reporter assays, transcriptome analysis, and genetic interactions. Although we were unable to determine a function for the first ZnF domain, humanization of the second ZnF domain to mirror that found in LUC7L or LUC7L2 resulted in altered usage of nonconsensus 5'ss. In contrast, the corresponding ZnF domain of LUC7L3 could not support yeast viability. Further, humanization of Luc7 can suppress mutation of the ATPase Prp28, which is involved in U1 release and exchange for U6 at the 5'ss. Our work reveals a role for the second ZnF of Luc7 in splice site selection and suggests that different ZnF domains may have different ATPase requirements for release by Prp28.

A role for the C. elegans Argonaute protein CSR-1 in small nuclear RNA 3' processing.

The Integrator is a multi-subunit protein complex that catalyzes the maturation of snRNA transcripts via 3' cleavage, a step required for snRNA incorporation with snRNP for spliceosome biogenesis. Here we developed a GFP based in vivo snRNA misprocessing reporter as a readout of Integrator function and performed a genome-wide RNAi screen for Integrator regulators. We found that loss of the Argonaute encoding csr-1 gene resulted in widespread 3' misprocessing of snRNA transcripts that is accompanied by a significant increase in alternative splicing. Loss of the csr-1 gene down-regulates the germline expression of Integrator subunits 4 and 6 and is accompanied by a reduced protein translation efficiency of multiple Integrator catalytic and non-catalytic subunits. Through isoform and motif mutant analysis, we determined that CSR-1's effect on snRNA processing is dependent on its catalytic slicer activity but does not involve the CSR-1a isoform. Moreover, mRNA-sequencing revealed high similarity in the transcriptome profile between csr-1 and Integrator subunit knockdown via RNAi. Together, our findings reveal CSR-1 as a new regulator of the Integrator complex and implicate a novel role of this Argonaute protein in snRNA 3' processing.

Single-cell and single-nucleus RNA-sequencing from paired normal-adenocarcinoma lung samples provide both common and discordant biological insights.

Whether single-cell RNA-sequencing (scRNA-seq) captures the same biological information as single-nucleus RNA-sequencing (snRNA-seq) remains uncertain and likely to be context-dependent. Herein, a head-to-head comparison was performed in matched normal-adenocarcinoma human lung samples to assess biological insights derived from scRNA-seq versus snRNA-seq and better understand the cellular transition that occurs from normal to tumoral tissue. Here, the transcriptome of 160,621 cells/nuclei was obtained. In non-tumor lung, cell type proportions varied widely between scRNA-seq and snRNA-seq with a predominance of immune cells in the former (81.5%) and epithelial cells (69.9%) in the later. Similar results were observed in adenocarcinomas, in addition to an overall increase in cell type heterogeneity and a greater prevalence of copy number variants in cells of epithelial origin, which suggests malignant assignment. The cell type transition that occurs from normal lung tissue to adenocarcinoma was not always concordant whether cells or nuclei were examined. As expected, large differential expression of the whole-cell and nuclear transcriptome was observed, but cell-type specific changes of paired normal and tumor lung samples revealed a set of common genes in the cells and nuclei involved in cancer-related pathways. In addition, we showed that the ligand-receptor interactome landscape of lung adenocarcinoma was largely different whether cells or nuclei were evaluated. Immune cell depletion in fresh specimens partly mitigated the difference in cell type composition observed between cells and nuclei. However, the extra manipulations affected cell viability and amplified the transcriptional signatures associated with stress responses. In conclusion, research applications focussing on mapping the immune landscape of lung adenocarcinoma benefit from scRNA-seq in fresh samples, whereas snRNA-seq of frozen samples provide a low-cost alternative to profile more epithelial and cancer cells, and yield cell type proportions that more closely match tissue content.

Spliceosomal snRNAs, the Essential Players in pre-mRNA Processing in Eukaryotic Nucleus: From Biogenesis to Functions and Spatiotemporal Characteristics.

Spliceosomal small nuclear RNAs (snRNAs) are a fundamental class of non-coding small RNAs abundant in the nucleoplasm of eukaryotic cells, playing a crucial role in splicing precursor messenger RNAs (pre-mRNAs). They are transcribed by DNA-dependent RNA polymerase II (Pol II) or III (Pol III), and undergo subsequent processing and 3' end cleavage to become mature snRNAs. Numerous protein factors are involved in the transcription initiation, elongation, termination, splicing, cellular localization, and terminal modification processes of snRNAs. The transcription and processing of snRNAs are regulated spatiotemporally by various mechanisms, and the homeostatic balance of snRNAs within cells is of great significance for the growth and development of organisms. snRNAs assemble with specific accessory proteins to form small nuclear ribonucleoprotein particles (snRNPs) that are the basal components of spliceosomes responsible for pre-mRNA maturation. This article provides an overview of the biological functions, biosynthesis, terminal structure, and tissue-specific regulation of snRNAs.

PRPF8-mediated dysregulation of hBrr2 helicase disrupts human spliceosome kinetics and 5´-splice-site selection causing tissue-specific defects.

The carboxy-terminus of the spliceosomal protein PRPF8, which regulates the RNA helicase Brr2, is a hotspot for mutations causing retinitis pigmentosa-type 13, with unclear role in human splicing and tissue-specificity mechanism. We used patient induced pluripotent stem cells-derived cells, carrying the heterozygous PRPF8 c.6926 A > C (p.H2309P) mutation to demonstrate retinal-specific endophenotypes comprising photoreceptor loss, apical-basal polarity and ciliary defects. Comprehensive molecular, transcriptomic, and proteomic analyses revealed a role of the PRPF8/Brr2 regulation in 5'-splice site (5'SS) selection by spliceosomes, for which disruption impaired alternative splicing and weak/suboptimal 5'SS selection, and enhanced cryptic splicing, predominantly in ciliary and retinal-specific transcripts. Altered splicing efficiency, nuclear speckles organisation, and PRPF8 interaction with U6 snRNA, caused accumulation of active spliceosomes and poly(A)+ mRNAs in unique splicing clusters located at the nuclear periphery of photoreceptors. Collectively these elucidate the role of PRPF8/Brr2 regulatory mechanisms in splicing and the molecular basis of retinal disease, informing therapeutic approaches.

Structural insights into human exon-defined spliceosome prior to activation.

Spliceosome is often assembled across an exon and undergoes rearrangement to span a neighboring intron. Most states of the intron-defined spliceosome have been structurally characterized. However, the structure of a fully assembled exon-defined spliceosome remains at large. During spliceosome assembly, the pre-catalytic state (B complex) is converted from its precursor (pre-B complex). Here we report atomic structures of the exon-defined human spliceosome in four sequential states: mature pre-B, late pre-B, early B, and mature B. In the previously unknown late pre-B state, U1 snRNP is already released but the remaining proteins are still in the pre-B state; unexpectedly, the RNAs are in the B state, with U6 snRNA forming a duplex with 5'-splice site and U5 snRNA recognizing the 3'-end of the exon. In the early and mature B complexes, the B-specific factors are stepwise recruited and specifically recognize the exon 3'-region. Our study reveals key insights into the assembly of the exon-defined spliceosomes and identifies mechanistic steps of the pre-B-to-B transition.

Locus-specific detection of pseudouridine with CRISPR-Cas13a.

We combined the CRISPR-Cas13a system with CMC chemical labeling, developing an approach that enables precise identification of pseudouridine (Ψ) sites at specific loci within ribosomal RNA (rRNA), messenger RNA (mRNA) and small nuclear RNAs (snRNA). This method, with good efficiency and simplicity, detects Ψ sites through fluorescence measurement, providing a straightforward and fast validation for targeted Ψ sites of interest.

A role for SNU66 in maintaining 5' splice site identity during spliceosome assembly.

In spliceosome assembly, the 5' splice site is initially recognized by U1 snRNA. U1 leaves the spliceosome during the assembly process, therefore other factors contribute to the maintenance of 5' splice site identity as it is loaded into the catalytic site. Recent structural data suggest that human tri-snRNP 27K (SNRP27) M141 and SNU66 H734 interact to stabilize the U4/U6 quasi-pseudo knot at the base of the U6 snRNA ACAGAGA box in pre-B complex. Previously, we found that mutations in Caenorhabditis elegans at SNRP-27 M141 promote changes in alternative 5'ss usage. We tested whether the potential interaction between SNRP-27 M141 and SNU-66 H765 (the C. elegans equivalent position to human SNU66 H734) contributes to maintaining 5' splice site identity during spliceosome assembly. We find that SNU-66 H765 mutants promote alternative 5' splice site usage. Many of the alternative 5' splicing events affected by SNU-66(H765G) overlap with those affected SNRP-27(M141T). Double mutants of snrp-27(M141T) and snu-66(H765G) are homozygous lethal. We hypothesize that mutations at either SNRP-27 M141 or SNU-66 H765 allow the spliceosome to load alternative 5' splice sites into the active site. Tests with mutant U1 snRNA and swapped 5' splice sites indicate that the ability of SNRP-27 M141 and SNU-66 H765 mutants to affect a particular 5' splice alternative splicing event is dependent on both the presence of a weaker consensus 5'ss nearby and potentially nearby splicing factor binding sites. Our findings confirm a new role for the C terminus of SNU-66 in maintenance of 5' splice site identity during spliceosome assembly.

Splicing analysis of STAT3 tandem donor suggests non-canonical binding registers for U1 and U6 snRNAs.

Tandem donor splice sites (5'ss) are unique regions with at least two GU dinucleotides serving as splicing cleavage sites. The Δ3 tandem 5'ss are a specific subclass of 5'ss separated by 3 nucleotides which can affect protein function by inserting/deleting a single amino acid. One 5'ss is typically preferred, yet factors governing particular 5'ss choice are not fully understood. A highly conserved exon 21 of the STAT3 gene was chosen as a model to study Δ3 tandem 5'ss splicing mechanisms. Based on multiple lines of experimental evidence, endogenous U1 snRNA most likely binds only to the upstream 5'ss. However, the downstream 5'ss is used preferentially, and the splice site choice is not dependent on the exact U1 snRNA binding position. Downstream 5'ss usage was sensitive to exact nucleotide composition and dependent on the presence of downstream regulatory region. The downstream 5'ss usage could be best explained by two novel interactions with endogenous U6 snRNA. U6 snRNA enables the downstream 5'ss usage in STAT3 exon 21 by two mechanisms: (i) binding in a novel non-canonical register and (ii) establishing extended Watson-Crick base pairing with the downstream regulatory region. This study suggests that U6:5'ss interaction is more flexible than previously thought.

Unconventional features in the transcription and processing of spliceosomal small nuclear RNAs in the protozoan parasite Trichomonas vaginalis.

Trichomonas vaginalis is a medically important protozoan parasite, and a deep-branching, evolutionarily divergent unicellular eukaryote that has conserved several key features of eukaryotic gene expression. Trichomonas vaginalis possesses a metazoan/plant-like capping apparatus, mRNAs with a cap 1 structure and spliceosomes containing the five small nuclear RNAs (snRNAs). However, in contrast to metazoan and plant snRNAs, the structurally conserved T. vaginalis snRNAs were initially identified as lacking the canonical guanosine cap nucleotide. To explain this unusual condition, we sought to investigate transcriptional and processing features of the spliceosomal snRNAs in this protist. Here, we show that T. vaginalis spliceosomal snRNA genes mostly lack typical eukaryotic promoters. In contrast to other eukaryotes, the putative TATA box in the T. vaginalis U6 snRNA gene was found to be dispensable for transcription or RNA polymerase selectivity. Moreover, U6 transcription in T. vaginalis was virtually insensitive to tagetitoxin compared with other cellular transcripts produced by the same RNA polymerase III. Most important and unexpected, snRNA transcription in T. vaginalis appears to bypass capping as we show that these transcripts retain their original 5'-triphosphate groups. In conclusion, transcription and processing of spliceosomal snRNAs in T. vaginalis deviate considerably from the conventional rules of other eukaryotes.

Distinct functions for the paralogous RBM41 and U11/U12-65K proteins in the minor spliceosome.

Here, we identify RBM41 as a novel unique protein component of the minor spliceosome. RBM41 has no previously recognized cellular function but has been identified as a paralog of U11/U12-65K, a known unique component of the U11/U12 di-snRNP. Both proteins use their highly similar C-terminal RRMs to bind to 3'-terminal stem-loops in U12 and U6atac snRNAs with comparable affinity. Our BioID data indicate that the unique N-terminal domain of RBM41 is necessary for its association with complexes containing DHX8, an RNA helicase, which in the major spliceosome drives the release of mature mRNA from the spliceosome. Consistently, we show that RBM41 associates with excised U12-type intron lariats, is present in the U12 mono-snRNP, and is enriched in Cajal bodies, together suggesting that RBM41 functions in the post-splicing steps of the minor spliceosome assembly/disassembly cycle. This contrasts with U11/U12-65K, which uses its N-terminal region to interact with U11 snRNP during intron recognition. Finally, while RBM41 knockout cells are viable, they show alterations in U12-type 3' splice site usage. Together, our results highlight the role of the 3'-terminal stem-loop of U12 snRNA as a dynamic binding platform for the U11/U12-65K and RBM41 proteins, which function at distinct stages of the assembly/disassembly cycle.

scaRNA20 promotes pseudouridylatory modification of small nuclear snRNA U12 and improves cardiomyogenesis.

Non-coding RNAs, particularly small Cajal-body associated RNAs (scaRNAs), play a significant role in spliceosomal RNA modifications. While their involvement in ischemic myocardium regeneration is known, their role in cardiac development is unexplored. We investigated scaRNA20's role in iPSC differentiation into cardiomyocytes (iCMCs) via overexpression and knockdown assays. We measured scaRNA20-OE-iCMCs and scaRNA20-KD-iCMCs contractility using Particle Image Velocimetry (PIV), comparing them to control iCMCs. We explored scaRNA20's impact on alternative splicing via pseudouridylation (Ψ) of snRNA U12, analyzing its functional consequences in cardiac differentiation. scaRNA20-OE-iPSC differentiation increased beating colonies, upregulated cardiac-specific genes, activated TP53 and STAT3, and preserved contractility under hypoxia. Conversely, scaRNA20-KD-iCMCs exhibited poor differentiation and contractility. STAT3 inhibition in scaRNA20-OE-iPSCs hindered cardiac differentiation. RNA immunoprecipitation revealed increased Ψ at the 28th uridine of U12 RNA in scaRNA20-OE iCMCs. U12-KD iCMCs had reduced cardiac differentiation, which improved upon U12 RNA introduction. In summary, scaRNA20-OE in iPSCs enhances cardiomyogenesis, preserves iCMC function under hypoxia, and may have implications for ischemic myocardium regeneration.

Sm complex assembly and 5' cap trimethylation promote selective processing of snRNAs by the 3' exonuclease TOE1.

Competing exonucleases that promote 3' end maturation or degradation direct quality control of small non-coding RNAs, but how these enzymes distinguish normal from aberrant RNAs is poorly understood. The Pontocerebellar Hypoplasia 7 (PCH7)-associated 3' exonuclease TOE1 promotes maturation of canonical small nuclear RNAs (snRNAs). Here, we demonstrate that TOE1 achieves specificity toward canonical snRNAs through their Sm complex assembly and cap trimethylation, two features that distinguish snRNAs undergoing correct biogenesis from other small non-coding RNAs. Indeed, disruption of Sm complex assembly via snRNA mutations or protein depletions obstructs snRNA processing by TOE1, and in vitro snRNA processing by TOE1 is stimulated by a trimethylated cap. An unstable snRNA variant that normally fails to undergo maturation becomes fully processed by TOE1 when its degenerate Sm binding motif is converted into a canonical one. Our findings uncover the molecular basis for how TOE1 distinguishes snRNAs from other small non-coding RNAs and explain how TOE1 promotes maturation specifically of canonical snRNAs undergoing proper processing.

Examining the capacity of human U1 snRNA variants to facilitate pre-mRNA splicing.

The human U1 snRNA is encoded by a multigene family consisting of transcribed variants and defective pseudogenes. Many variant U1 (vU1) snRNAs have been demonstrated to not only be transcribed but also processed by the addition of a trimethylated guanosine cap, packaged into snRNPs, and assembled into spliceosomes; however, their capacity to facilitate pre-mRNA splicing has, so far, not been tested. A recent systematic analysis of the human snRNA genes identified 178 U1 snRNA genes that are present in the genome as either tandem arrays or single genes on multiple chromosomes. Of these, 15 were found to be expressed in human tissues and cell lines, although at significantly low levels from their endogenous loci, <0.001% of the canonical U1 snRNA. In this study, we found that placing the variants in the context of the regulatory elements of the RNU1-1 gene improves the expression of many variants to levels comparable to the canonical U1 snRNA. Application of a previously established HeLa cell-based minigene reporter assay to examine the capacity of the vU1 snRNAs to support pre-mRNA splicing revealed that even though the exogenously expressed variant snRNAs were enriched in the nucleus, only a few had a measurable effect on splicing.

Visualizing a two-state conformational ensemble in stem-loop 3 of the transcriptional regulator 7SK RNA.

Structural plasticity is integral to RNA function; however, there are currently few methods to quantitatively resolve RNAs that have multiple structural states. NMR spectroscopy is a powerful approach for resolving conformational ensembles but is size-limited. Chemical probing is well-suited for large RNAs but provides limited structural and kinetics information. Here, we integrate the two approaches to visualize a two-state conformational ensemble for the central stem-loop 3 (SL3) of 7SK RNA, a critical element for 7SK RNA function in transcription regulation. We find that the SL3 distal end exchanges between two equally populated yet structurally distinct states in both isolated SL3 constructs and full-length 7SK RNA. We rationally designed constructs that lock SL3 into a single state and demonstrate that both chemical probing and NMR data fit to a linear combination of the two states. Comparison of vertebrate 7SK RNA sequences shows either or both states are highly conserved. These results provide new insights into 7SK RNA structural dynamics and demonstrate the utility of integrating chemical probing with NMR spectroscopy to gain quantitative insights into RNA conformational ensembles.

Global analysis of binding sites of U2AF1 and ZRSR2 reveals RNA elements required for mutually exclusive splicing by the U2- and U12-type spliceosome.

Recurring mutations in genes encoding 3' splice-site recognition proteins, U2AF1 and ZRSR2 are associated with human cancers. Here, we determined binding sites of the proteins to reveal that U2-type and U12-type splice sites are recognized by U2AF1 and ZRSR2, respectively. However, some sites are spliced by both the U2-type and U12-type spliceosomes, indicating that well-conserved consensus motifs in some U12-type introns could be recognized by the U2-type spliceosome. Nucleotides flanking splice sites of U12-type introns are different from those flanking U2-type introns. Remarkably, the AG dinucleotide at the positions -1 and -2 of 5' splice sites of U12-type introns with GT-AG termini is not present. AG next to 5' splice site introduced by a single nucleotide substitution at the -2 position could convert a U12-type splice site to a U2-type site. The class switch of introns by a single mutation and the bias against G at the -1 position of U12-type 5' splice site support the notion that the identities of nucleotides in exonic regions adjacent to splice sites are fine-tuned to avoid recognition by the U2-type spliceosome. These findings may shed light on the mechanism of selectivity in U12-type intron splicing and the mutations that affect splicing.

snRNA-seq analysis in multinucleated myogenic FSHD cells identifies heterogeneous FSHD transcriptome signatures associated with embryonic-like program activation and oxidative stress-induced apoptosis.

The sporadic nature of DUX4 expression in FSHD muscle challenges comparative transcriptome analyses between FSHD and control samples. A variety of DUX4 and FSHD-associated transcriptional changes have been identified, but bulk RNA-seq strategies prohibit comprehensive analysis of their spatiotemporal relation, interdependence and role in the disease process. In this study, we used single-nucleus RNA-sequencing of nuclei isolated from patient- and control-derived multinucleated primary myotubes to investigate the cellular heterogeneity in FSHD. Taking advantage of the increased resolution in snRNA-sequencing of fully differentiated myotubes, two distinct populations of DUX4-affected nuclei could be defined by their transcriptional profiles. Our data provides insights into the differences between these two populations and suggests heterogeneity in two well-known FSHD-associated transcriptional aberrations: increased oxidative stress and inhibition of myogenic differentiation. Additionally, we provide evidence that DUX4-affected nuclei share transcriptome features with early embryonic cells beyond the well-described cleavage stage, progressing into the 8-cell and blastocyst stages. Altogether, our data suggests that the FSHD transcriptional profile is defined by a mixture of individual and sometimes mutually exclusive DUX4-induced responses and cellular state-dependent downstream effects.