Identifying cis Elements for Spatiotemporal Control of Mammalian DNA Replication.

The temporal order of DNA replication (replication timing [RT]) is highly coupled with genome architecture, but cis-elements regulating either remain elusive. We created a series of CRISPR-mediated deletions and inversions of a pluripotency-associated topologically associating domain (TAD) in mouse ESCs. CTCF-associated domain boundaries were dispensable for RT. CTCF protein depletion weakened most TAD boundaries but had no effect on RT or A/B compartmentalization genome-wide. By contrast, deletion of three intra-TAD CTCF-independent 3D contact sites caused a domain-wide early-to-late RT shift, an A-to-B compartment switch, weakening of TAD architecture, and loss of transcription. The dispensability of TAD boundaries and the necessity of these "early replication control elements" (ERCEs) was validated by deletions and inversions at additional domains. Our results demonstrate that discrete cis-regulatory elements orchestrate domain-wide RT, A/B compartmentalization, TAD architecture, and transcription, revealing fundamental principles linking genome structure and function.

miRNA analysis reveals novel dysregulated pathways in amyotrophic lateral sclerosis.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease. Its complex pathogenesis and phenotypic heterogeneity hinder therapeutic development and early diagnosis. Altered RNA metabolism is a recurrent pathophysiologic theme, including distinct microRNA (miRNA) profiles in ALS tissues. We profiled miRNAs in accessible biosamples, including skin fibroblasts and whole blood and compared them in age- and sex-matched healthy controls versus ALS participants with and without repeat expansions to chromosome 9 open reading frame 72 (C9orf72; C9-ALS and nonC9-ALS), the most frequent ALS mutation. We identified unique and shared profiles of differential miRNA (DmiRNA) levels in each C9-ALS and nonC9-ALS tissues versus controls. Fibroblast DmiRNAs were validated by quantitative real-time PCR and their target mRNAs by 5-bromouridine and 5-bromouridine-chase sequencing. We also performed pathway analysis to infer biological meaning, revealing anticipated, tissue-specific pathways and pathways previously linked to ALS, as well as novel pathways that could inform future research directions. Overall, we report a comprehensive study of a miRNA profile dataset from C9-ALS and nonC9-ALS participants across two accessible biosamples, providing evidence of dysregulated miRNAs in ALS and possible targets of interest. Distinct miRNA patterns in accessible tissues may also be leveraged to distinguish ALS participants from healthy controls for earlier diagnosis. Future directions may look at potential correlations of miRNA profiles with clinical parameters.

Type I interferons and MAVS signaling are necessary for tissue resident memory CD8+ T cell responses to RSV infection.

Respiratory syncytial virus (RSV) can cause bronchiolitis and viral pneumonia in young children and the elderly. Lack of vaccines and recurrence of RSV infection indicate the difficulty in eliciting protective memory immune responses. Tissue resident memory T cells (TRM) can confer protection from pathogen re-infection and, in human experimental RSV infection, the presence of lung CD8+ TRM cells correlates with a better outcome. However, the requirements for generating and maintaining lung TRM cells during RSV infection are not fully understood. Here, we use mouse models to assess the impact of innate immune response determinants in the generation and subsequent expansion of the TRM cell pool during RSV infection. We show that CD8+ TRM cells expand independently from systemic CD8+ T cells after RSV re-infection. Re-infected MAVS and MyD88/TRIF deficient mice, lacking key components involved in innate immune recognition of RSV and induction of type I interferons (IFN-α/ß), display impaired expansion of CD8+ TRM cells and reduction in antigen specific production of granzyme B and IFN-γ. IFN-α treatment of MAVS deficient mice during primary RSV infection restored TRM cell expansion upon re-challenge but failed to recover TRM cell functionality. Our data reveal how innate immunity, including the axis controlling type I IFN induction, instructs and regulates CD8+ TRM cell responses to RSV infection, suggesting possible mechanisms for therapeutic intervention.

KMT2D links TGF-ß signaling to noncanonical activin pathway and regulates pancreatic cancer cell plasticity.

Although KMT2D, also known as MLL2, is known to play an essential role in development, differentiation, and tumor suppression, its role in pancreatic cancer development is not well understood. Here, we discovered a novel signaling axis mediated by KMT2D, which links TGF-ß to the activin A pathway. We found that TGF-ß upregulates a microRNA, miR-147b, which in turn leads to post-transcriptional silencing of KMT2D. Loss of KMT2D induces the expression and secretion of activin A, which activates a noncanonical p38 MAPK-mediated pathway to modulate cancer cell plasticity, promote a mesenchymal phenotype, and enhance tumor invasion and metastasis in mice. We observed a decreased KMT2D expression in human primary and metastatic pancreatic cancer. Furthermore, inhibition or knockdown of activin A reversed the protumoral role of KMT2D loss. These findings support a tumor-suppressive role of KMT2D in pancreatic cancer and identify miR-147b and activin A as novel therapeutic targets.

Co-transcriptional splicing efficiencies differ within genes and between cell types.

Pre-mRNA splicing is carried out by the spliceosome and involves splice site recognition, removal of introns, and ligation of exons. Components of the spliceosome have been shown to interact with the elongating RNA polymerase II (RNAPII) which is thought to allow splicing to occur concurrently with transcription. However, little is known about the regulation and efficiency of co-transcriptional splicing in human cells. In this study, we used Bru-seq and BruChase-seq to determine the co-transcriptional splicing efficiencies of 17,000 introns expressed across 6 human cell lines. We found that less than half of all introns across these 6 cell lines were co-transcriptionally spliced. Splicing efficiencies for individual introns showed variations across cell lines, suggesting that splicing may be regulated in a cell-type specific manner. Moreover, the splicing efficiency of introns varied within genes. The efficiency of co-transcriptional splicing did not correlate with gene length, intron position, splice site strengths, or the intron/neighboring exons GC content. However, we identified binding signals from multiple RNA binding proteins (RBPs) that correlated with splicing efficiency, including core spliceosomal machinery components-such as SF3B4, U2AF1 and U2AF2 showing higher binding signals in poorly spliced introns. In addition, multiple RBPs, such as BUD13, PUM1 and SND1, showed preferential binding in exons that flank introns with high splicing efficiencies. The nascent RNA splicing patterns presented here across multiple cell types add to our understanding of the complexity in RNA splicing, wherein RNA-binding proteins may play important roles in determining splicing outcomes in a cell type- and intron-specific manner.

GRHL2-controlled gene expression networks in luminal breast cancer.

Grainyhead like 2 (GRHL2) is an essential transcription factor for development and function of epithelial tissues. It has dual roles in cancer by supporting tumor growth while suppressing epithelial to mesenchymal transitions (EMT). GRHL2 cooperates with androgen and estrogen receptors (ER) to regulate gene expression. We explore genome wide GRHL2 binding sites conserved in three ER⍺/GRHL2 positive luminal breast cancer cell lines by ChIP-Seq. Interaction with the ER⍺/FOXA1/GATA3 complex is observed, however, only for a minor fraction of conserved GRHL2 peaks. We determine genome wide transcriptional dynamics in response to loss of GRHL2 by nascent RNA Bru-seq using an MCF7 conditional knockout model. Integration of ChIP- and Bru-seq pinpoints candidate direct GRHL2 target genes in luminal breast cancer. Multiple connections between GRHL2 and proliferation are uncovered, including transcriptional activation of ETS and E2F transcription factors. Among EMT-related genes, direct regulation of CLDN4 is corroborated but several targets identified in other cells (including CDH1 and ZEB1) are ruled out by both ChIP- and Bru-seq as being directly controlled by GRHL2 in luminal breast cancer cells. Gene clusters correlating positively (including known GRHL2 targets such as ErbB3, CLDN4/7) or negatively (including TGFB1 and TGFBR2) with GRHL2 in the MCF7 knockout model, display similar correlation with GRHL2 in ER positive as well as ER negative breast cancer patients. Altogether, this study uncovers gene sets regulated directly or indirectly by GRHL2 in luminal breast cancer, identifies novel GRHL2-regulated genes, and points to distinct GRHL2 regulation of EMT in luminal breast cancer cells. Video Abstract.

Principles of mRNA control by human PUM proteins elucidated from multimodal experiments and integrative data analysis.

The human PUF-family proteins, PUM1 and PUM2, posttranscriptionally regulate gene expression by binding to a PUM recognition element (PRE) in the 3'-UTR of target mRNAs. Hundreds of PUM1/2 targets have been identified from changes in steady-state RNA levels; however, prior studies could not differentiate between the contributions of changes in transcription and RNA decay rates. We applied metabolic labeling to measure changes in RNA turnover in response to depletion of PUM1/2, showing that human PUM proteins regulate expression almost exclusively by changing RNA stability. We also applied an in vitro selection workflow to precisely identify the binding preferences of PUM1 and PUM2. By integrating our results with prior knowledge, we developed a "rulebook" of key contextual features that differentiate functional versus nonfunctional PREs, allowing us to train machine learning models that accurately predict the functional regulation of RNA targets by the human PUM proteins.

Lack of IL-1 Receptor Signaling Reduces Spontaneous Airway Eosinophilia in Juvenile Mice with Muco-Obstructive Lung Disease.

Previous studies demonstrated spontaneous type 2 airway inflammation with eosinophilia in juvenile Scnn1b (sodium channel, non-voltage-gated 1, ß-subunit)-transgenic (Scnn1b-Tg) mice with muco-obstructive lung disease. IL-1 receptor (IL-1R) signaling has been implicated in allergen-driven airway disease; however, its role in eosinophilic inflammation in muco-obstructive lung disease remains unknown. In this study, we examined the role of IL-1R signaling in the development of airway eosinophilia and type 2 inflammation in juvenile Scnn1b-Tg mice. We determined effects of genetic deletion of Il1r1 (IL-1 receptor type I) on eosinophil counts, transcript levels of key type 2 cytokines, markers of eosinophil activation and apoptosis, and tissue morphology in lungs of Scnn1b-Tg mice at different time points during neonatal development. Furthermore, we measured endothelial surface expression of intercellular adhesion molecule 1 (ICAM-1), an integrin involved in eosinophil transendothelial migration, and determined effects of eosinophil depletion using an anti-IL-5 antibody on lung morphology. Lack of IL-1R reduced airway eosinophilia and structural lung damage, but it did not reduce concentrations of type 2 cytokines and associated eosinophil activation in Scnn1b-Tg mice. Structural lung damage in Scnn1b-Tg mice was also reduced by eosinophil depletion. Lack of IL-1R was associated with reduced expression of ICAM-1 on lung endothelial cells and reduced eosinophil counts in lungs from Scnn1b-Tg mice. We conclude that IL-1R signaling is implicated in airway eosinophilia independent of type 2 cytokines in juvenile Scnn1b-Tg mice. Our data suggest that IL-1R signaling may be relevant in the pathogenesis of eosinophilic airway inflammation in muco-obstructive lung diseases, which may be mediated in part by ICAM-1-dependent transmigration of eosinophils into the lungs.

Large transcription units unify copy number variants and common fragile sites arising under replication stress.

Copy number variants (CNVs) resulting from genomic deletions and duplications and common fragile sites (CFSs) seen as breaks on metaphase chromosomes are distinct forms of structural chromosome instability precipitated by replication inhibition. Although they share a common induction mechanism, it is not known how CNVs and CFSs are related or why some genomic loci are much more prone to their occurrence. Here we compare large sets of de novo CNVs and CFSs in several experimental cell systems to each other and to overlapping genomic features. We first show that CNV hotpots and CFSs occurred at the same human loci within a given cultured cell line. Bru-seq nascent RNA sequencing further demonstrated that although genomic regions with low CNV frequencies were enriched in transcribed genes, the CNV hotpots that matched CFSs specifically corresponded to the largest active transcription units in both human and mouse cells. Consistently, active transcription units >1 Mb were robust cell-type-specific predictors of induced CNV hotspots and CFS loci. Unlike most transcribed genes, these very large transcription units replicated late and organized deletion and duplication CNVs into their transcribed and flanking regions, respectively, supporting a role for transcription in replication-dependent lesion formation. These results indicate that active large transcription units drive extreme locus- and cell-type-specific genomic instability under replication stress, resulting in both CNVs and CFSs as different manifestations of perturbed replication dynamics.

Oncogenic fusion protein EWS-FLI1 is a network hub that regulates alternative splicing.

The synthesis and processing of mRNA, from transcription to translation initiation, often requires splicing of intragenic material. The final mRNA composition varies based on proteins that modulate splice site selection. EWS-FLI1 is an Ewing sarcoma (ES) oncoprotein with an interactome that we demonstrate to have multiple partners in spliceosomal complexes. We evaluate the effect of EWS-FLI1 on posttranscriptional gene regulation using both exon array and RNA-seq. Genes that potentially regulate oncogenesis, including CLK1, CASP3, PPFIBP1, and TERT, validate as alternatively spliced by EWS-FLI1. In a CLIP-seq experiment, we find that EWS-FLI1 RNA-binding motifs most frequently occur adjacent to intron-exon boundaries. EWS-FLI1 also alters splicing by directly binding to known splicing factors including DDX5, hnRNP K, and PRPF6. Reduction of EWS-FLI1 produces an isoform of γ-TERT that has increased telomerase activity compared with wild-type (WT) TERT. The small molecule YK-4-279 is an inhibitor of EWS-FLI1 oncogenic function that disrupts specific protein interactions, including helicases DDX5 and RNA helicase A (RHA) that alters RNA-splicing ratios. As such, YK-4-279 validates the splicing mechanism of EWS-FLI1, showing alternatively spliced gene patterns that significantly overlap with EWS-FLI1 reduction and WT human mesenchymal stem cells (hMSC). Exon array analysis of 75 ES patient samples shows similar isoform expression patterns to cell line models expressing EWS-FLI1, supporting the clinical relevance of our findings. These experiments establish systemic alternative splicing as an oncogenic process modulated by EWS-FLI1. EWS-FLI1 modulation of mRNA splicing may provide insight into the contribution of splicing toward oncogenesis, and, reciprocally, EWS-FLI1 interactions with splicing proteins may inform the splicing code.

Rate of elongation by RNA polymerase II is associated with specific gene features and epigenetic modifications.

The rate of transcription elongation plays an important role in the timing of expression of full-length transcripts as well as in the regulation of alternative splicing. In this study, we coupled Bru-seq technology with 5,6-dichlorobenzimidazole 1-ß-D-ribofuranoside (DRB) to estimate the elongation rates of over 2000 individual genes in human cells. This technique, BruDRB-seq, revealed gene-specific differences in elongation rates with a median rate of around 1.5 kb/min. We found that genes with rapid elongation rates showed higher densities of H3K79me2 and H4K20me1 histone marks compared to slower elongating genes. Furthermore, high elongation rates had a positive correlation with gene length, low complexity DNA sequence, and distance from the nearest active transcription unit. Features that negatively correlated with elongation rate included the density of exons, long terminal repeats, GC content of the gene, and DNA methylation density in the bodies of genes. Our results suggest that some static gene features influence transcription elongation rates and that cells may alter elongation rates by epigenetic regulation. The BruDRB-seq technique offers new opportunities to interrogate mechanisms of regulation of transcription elongation.

DNA repair and recovery of RNA synthesis following exposure to ultraviolet light are delayed in long genes.

The kinetics of DNA repair and RNA synthesis recovery in human cells following UV-irradiation were assessed using nascent RNA Bru-seq and quantitative long PCR. It was found that UV light inhibited transcription elongation and that recovery of RNA synthesis occurred as a wave in the 5'-3' direction with slow recovery and TC-NER at the 3' end of long genes. RNA synthesis resumed fully at the 3'-end of genes after a 24 h recovery in wild-type fibroblasts, but not in cells deficient in transcription-coupled nucleotide excision repair (TC-NER) or global genomic NER (GG-NER). Different transcription recovery profiles were found for individual genes but these differences did not fully correlate to differences in DNA repair of these genes. Our study gives the first genome-wide view of how UV-induced lesions affect transcription and how the recovery of RNA synthesis of large genes are particularly delayed by the apparent lack of resumption of transcription by arrested polymerases.

Zinc finger protein 407 (ZFP407) regulates insulin-stimulated glucose uptake and glucose transporter 4 (Glut4) mRNA.

The glucose transporter GLUT4 facilitates insulin-stimulated glucose uptake in peripheral tissues including adipose, muscle, and heart. GLUT4 function is impaired in obesity and type 2 diabetes leading to hyperglycemia and an increased risk of cardiovascular disease and neuropathy. To better understand the regulation of GLUT4 function, a targeted siRNA screen was performed and led to the discovery that ZFP407 regulates insulin-stimulated glucose uptake in adipocytes. The decrease in insulin-stimulated glucose uptake due to ZFP407 deficiency was attributed to a reduction in GLUT4 mRNA and protein levels. The decrease in GLUT4 was due to both decreased transcription of Glut4 mRNA and decreased efficiency of Glut4 pre-mRNA splicing. Interestingly, ZFP407 coordinately regulated this decrease in transcription with an increase in the stability of Glut4 mRNA, resulting in opposing effects on steady-state Glut4 mRNA levels. More broadly, transcriptome analysis revealed that ZFP407 regulates many peroxisome proliferator-activated receptor (PPAR) γ target genes beyond Glut4. ZFP407 was required for the PPARγ agonist rosiglitazone to increase Glut4 expression, but was not sufficient to increase expression of a PPARγ target gene reporter construct. However, ZFP407 and PPARγ co-overexpression synergistically activated a PPARγ reporter construct beyond the level of PPARγ alone. Thus, ZFP407 may represent a new modulator of the PPARγ signaling pathway.

Coordinated regulation of synthesis and stability of RNA during the acute TNF-induced proinflammatory response.

Steady-state gene expression is a coordination of synthesis and decay of RNA through epigenetic regulation, transcription factors, micro RNAs (miRNAs), and RNA-binding proteins. Here, we present bromouride labeling and sequencing (Bru-Seq) and bromouridine pulse-chase and sequencing (BruChase-Seq) to assess genome-wide changes to RNA synthesis and stability in human fibroblasts at homeostasis and after exposure to the proinflammatory tumor necrosis factor (TNF). The inflammatory response in human cells involves rapid and dramatic changes in gene expression, and the Bru-Seq and BruChase-Seq techniques revealed a coordinated and complex regulation of gene expression both at the transcriptional and posttranscriptional levels. The combinatory analysis of both RNA synthesis and stability using Bru-Seq and BruChase-Seq allows for a much deeper understanding of mechanisms of gene regulation than afforded by the analysis of steady-state total RNA and should be useful in many biological settings.

Use of Bru-Seq and BruChase-Seq for genome-wide assessment of the synthesis and stability of RNA.

Gene expression studies commonly examine total cellular RNA, which only provides information about its steady-state pool of RNA. It remains unclear whether differences in the steady-state reflects variable rates of transcription or RNA degradation. To specifically monitor RNA synthesis and degradation genome-wide, we developed Bru-Seq and BruChase-Seq. These assays are based on metabolic pulse-chase labeling of RNA using bromouridine (Bru). In Bru-Seq, recently labeled RNAs are sequenced to reveal spans of nascent transcription in the genome. In BruChase-Seq, cells are chased in uridine for different periods of time following Bru-labeling, allowing for the isolation of RNA populations of specific ages. Here we describe these methodologies in detail and highlight their usefulness in assessing RNA synthesis and stability as well as splicing kinetics with examples of specific genes from different human cell lines.

Isoform and pathway-specific regulation of post-transcriptional RNA processing in human cells.

Steady-state levels of RNA transcripts are controlled by their rates of synthesis and degradation. Here we used nascent RNA Bru-seq and BruChase-seq to profile RNA dynamics across 16 human cell lines as part of ENCODE4 Deeply Profiled Cell Lines collection. We show that RNA turnover dynamics differ widely between transcripts of different genes and between different classes of RNA. Gene set enrichment analysis (GSEA) revealed that transcripts encoding proteins belonging to the same pathway often show similar turnover dynamics. Furthermore, transcript isoforms show distinct dynamics suggesting that RNA turnover is important in regulating mRNA isoform choice. Finally, splicing across newly made transcripts appears to be cooperative with either all or none type splicing. These data sets generated as part of ENCODE4 illustrate the intricate and coordinated regulation of RNA dynamics in controlling gene expression to allow for the precise coordination of cellular functions.

Characterizing nascent transcription patterns of PROMPTs, eRNAs, and readthrough transcripts in the ENCODE4 deeply profiled cell lines.

Arising as co-products of canonical gene expression, transcription-associated lincRNAs, such as promoter upstream transcripts (PROMPTs), enhancer RNAs (eRNAs), and readthrough (RT) transcripts, are often regarded as byproducts of transcription, although they may be important for the expression of nearby genes. We identified regions of nascent expression of these lincRNA in 16 human cell lines using Bru-seq techniques, and found distinctly regulated patterns of PROMPT, eRNA, and RT transcription using the diverse biochemical approaches in the ENCODE4 deeply profiled cell lines collection. Transcription of these lincRNAs was influenced by sequence-specific features and the local or 3D chromatin landscape. However, these sequence and chromatin features do not describe the full spectrum of lincRNA expression variability we identify, highlighting the complexity of their regulation. This may suggest that transcription-associated lincRNAs are not merely byproducts, but rather that the transcript itself, or the act of its transcription, is important for genomic function.

Implementation of a Multicomponent Otago-Based Virtual Fall Reduction, Education, and Exercise Program (MOVing FREEly) in Older Veterans.

PURPOSE: The COVID-19 pandemic limited access to community fall prevention programs, thus establishing the need for virtual interventions. Herein, we describe the feasibility, effectiveness, and acceptability of a virtual, multicomponent fall prevention program (MOVing FREEly). METHODS: A team of clinical falls prevention experts developed a six-week multicomponent fall prevention exercise and education class for older community-dwelling adults at risk of falling. Feasibility was measured through class attendance; effectiveness was measured through changes in performance measures, self-report of falling risk, and concern about falling; acceptability was assessed through questionnaires completed immediately upon program completion and at a three-month follow up. RESULTS: A total of 32 patients participated in the MOVing FREEly program. Attendance for education and exercise classes on average was greater than 80% with little attrition. Patient reported reduced concern of falling, improvement in the falls efficacy scale-international (FES-I) short form, and had statistically significant improvement in 30 s sit-to-stand and single-leg balance tests. The program was well received by participants, saving them significant time and costs of travel. CONCLUSIONS: A virtual, multicomponent fall prevention program is feasible and acceptable and effective as reducing falling risk. Future studies can explore the ability of this program to reduce falling incident and injury.

CDK12 regulates co-transcriptional splicing and RNA turnover in human cells.

The cyclin-dependent kinase CDK12 has garnered interest as a cancer therapeutic target as DNA damage response genes are particularly suppressed by loss of CDK12 activity. In this study, we assessed the acute effects of CDK12 inhibition on transcription and RNA processing using nascent RNA Bru-seq and BruChase-seq. Acute transcriptional changes were overall small after CDK12 inhibition but over 600 genes showed intragenic premature termination, including DNA repair and cell cycle genes. Furthermore, many genes showed reduced transcriptional readthrough past the end of genes in the absence of CDK12 activity. RNA turnover was dramatically affected by CDK12 inhibition and importantly, caused increased degradation of many transcripts from DNA damage response genes. We also show that co-transcriptional splicing was suppressed by CDK12 inhibition. Taken together, these studies reveal the roles of CDK12 in regulating transcription elongation, transcription termination, co-transcriptional splicing, and RNA turnover.

Induction of Genes Implicated in Stress Response and Autophagy by a Novel Quinolin-8-yl-nicotinamide QN523 in Pancreatic Cancer.

Using a cytotoxicity-based phenotypic screen of a highly diverse library of 20,000 small-molecule compounds, we identified a quinolin-8-yl-nicotinamide, QN519, as a promising lead. QN519 represents a novel scaffold with drug-like properties, showing potent in vitro cytotoxicity in a panel of 12 cancer cell lines. Subsequently, lead optimization campaign generated compounds with IC50 values < 1 µM. An optimized compound, QN523, shows significant in vivo efficacy in a pancreatic cancer xenograft model. QN523 treatment significantly increased the expression of HSPA5, DDIT3, TRIB3, and ATF3 genes, suggesting activation of the stress response pathway. We also observed a significant increase in the expression of WIPI1, HERPUD1, GABARAPL1, and MAP1LC3B, implicating autophagy as a major mechanism of action. Due to the lack of effective treatments for pancreatic cancer, discovery of novel agents such as the QN series of compounds with unique mechanism of action has the potential to fulfill a clear unmet medical need.