DUX4 is a transcription factor whose misexpression in skeletal muscle causes facioscapulohumeral muscular dystrophy (FSHD). While DUX4's transcriptional activity has been extensively characterized, the DUX4-induced proteome remains undescribed. Here, we report concurrent measurement of RNA and protein levels in DUX4-expressing cells via RNA-seq and quantitative mass spectrometry. DUX4 transcriptional targets were robustly translated, confirming the likely clinical relevance of proposed FSHD biomarkers. However, a multitude of mRNAs and proteins exhibited discordant expression changes upon DUX4 expression. Our dataset revealed unexpected proteomic, but not transcriptomic, dysregulation of diverse molecular pathways, including Golgi apparatus fragmentation, as well as extensive post-transcriptional buffering of stress response genes. Key components of RNA degradation machineries, including UPF1, UPF3B, and XRN1, exhibited suppressed protein, but not mRNA, levels, explaining the build-up of aberrant RNAs that characterizes DUX4-expressing cells. Our results provide a resource for the FSHD community and illustrate the importance of post-transcriptional processes to DUX4-induced pathology.
Mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository (Vizcaino et al., 2013) with the dataset identifier PXD010221. To enable easy access to processed peptide-spectrum match data, we have also deposited peptide-level data to Dryad (doi:10.5061/dryad.ck06k75). The previously published RNA-seq data are available through the NCBI SRA database under accession number GSE85461
Proteomics analysis of DUX4-expressing myoblastsPRIDE PRoteomics IDEntifications database, PXD010221.
Data from: Quantitative proteomics reveals key roles for post-transcriptional gene regulation in the molecular pathology of FSHD.Dryad Digital Repository, doi:10.5061/dryad.ck06k75.
Model systems of DUX4 expression recapitulate the transcriptional profile of FSHD cellsNCBI Gene Expression Omnibus, GSE85461.
- Stephen J Tapscott
- Robert K Bradley
- Sujatha Jagannathan
- Robert K Bradley
- Yuko Ogata
- Philip R Gafken
- Yuko Ogata
- Philip R Gafken
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
- Nahum Sonenberg, McGill University, Canada
© 2019, Jagannathan et al.
This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.
Alternative polyadenylation yields many mRNA isoforms whose 3' termini occur disproportionately in clusters within 3' UTRs. Previously, we showed that profiles of poly(A) site usage are regulated by the rate of transcriptional elongation by RNA polymerase (Pol) II (Geisberg et., 2020). Pol II derivatives with slow elongation rates confer an upstream-shifted poly(A) profile, whereas fast Pol II strains confer a downstream-shifted poly(A) profile. Within yeast isoform clusters, these shifts occur steadily from one isoform to the next across nucleotide distances. In contrast, the shift between clusters from the last isoform of one cluster to the first isoform of the next - is much less pronounced, even over large distances. GC content in a region 13-30 nt downstream from isoform clusters correlates with their sensitivity to Pol II elongation rate. In human cells, the upstream shift caused by a slow Pol II mutant also occurs continuously at the nucleotide level within clusters, but not between them. Pol II occupancy increases just downstream of the most speed-sensitive poly(A) sites, suggesting a linkage between reduced elongation rate and cluster formation. These observations suggest that 1) Pol II elongation speed affects the nucleotide-level dwell time allowing polyadenylation to occur, 2) poly(A) site clusters are linked to the local elongation rate and hence do not arise simply by intrinsically imprecise cleavage and polyadenylation of the RNA substrate, 3) DNA sequence elements can affect Pol II elongation and poly(A) profiles, and 4) the cleavage/polyadenylation and Pol II elongation complexes are spatially, and perhaps physically, coupled so that polyadenylation occurs rapidly upon emergence of the nascent RNA from the Pol II elongation complex.
South Asian women are at increased risk of developing gestational diabetes mellitus (GDM). Few studies have investigated the genetic contributions to GDM risk. We investigated the association of a type 2 diabetes (T2D) polygenic risk score (PRS), on its own, and with GDM risk factors, on GDM-related traits using data from two birth cohorts in which South Asian women were enrolled during pregnancy. 837 and 4372 pregnant South Asian women from the SouTh Asian BiRth CohorT (START) and Born in Bradford (BiB) cohort studies underwent a 75-g glucose tolerance test. PRSs were derived using genome-wide association study results from an independent multi-ethnic study (~18% South Asians). Associations with fasting plasma glucose (FPG); 2 hr post-load glucose (2hG); area under the curve glucose; and GDM were tested using linear and logistic regressions. The population attributable fraction (PAF) of the PRS was calculated. Every 1 SD increase in the PRS was associated with a 0.085 mmol/L increase in FPG ([95% confidence interval, CI=0.07–0.10], p=2.85×10−20); 0.21 mmol/L increase in 2hG ([95% CI=0.16–0.26], p=5.49×10−16); and a 45% increase in the risk of GDM ([95% CI=32–60%], p=2.27×10−14), independent of parental history of diabetes and other GDM risk factors. PRS tertile 3 accounted for 12.5% of the population’s GDM alone, and 21.7% when combined with family history. A few weak PRS and GDM risk factors interactions modulating FPG and GDM were observed. Taken together, these results show that a T2D PRS and family history of diabetes are strongly and independently associated with multiple GDM-related traits in women of South Asian descent, an effect that could be modulated by other environmental factors.