DNA methylation signatures reveal that distinct combinations of transcription factors specify human immune cell epigenetic identity.

Epigenetic reprogramming underlies specification of immune cell lineages, but patterns that uniquely define immune cell types and the mechanisms by which they are established remain unclear. Here, we identified lineage-specific DNA methylation signatures of six immune cell types from human peripheral blood and determined their relationship to other epigenetic and transcriptomic patterns. Sites of lineage-specific hypomethylation were associated with distinct combinations of transcription factors in each cell type. By contrast, sites of lineage-specific hypermethylation were restricted mostly to adaptive immune cells. PU.1 binding sites were associated with lineage-specific hypo- and hypermethylation in different cell types, suggesting that it regulates DNA methylation in a context-dependent manner. These observations indicate that innate and adaptive immune lineages are specified by distinct epigenetic mechanisms via combinatorial and context-dependent use of key transcription factors. The cell-specific epigenomics and transcriptional patterns identified serve as a foundation for future studies on immune dysregulation in diseases and aging.

Increased PTCHD4 expression via m6A modification of PTCHD4 mRNA promotes senescent cell survival.

RNA modifications, including N6-methyladenosine (m6A), critically modulate protein expression programs in a range of cellular processes. Although the transcriptomes of cells undergoing senescence are strongly regulated, the landscape and impact of m6A modifications during senescence are poorly understood. Here, we report a robust m6A modification of PTCHD4 mRNA, encoding Patched Domain-Containing Protein 4, in senescent cells. The METTL3/METTL14 complex was found to incorporate the m6A modification on PTCHD4 mRNA; addition of m6A rendered PTCHD4 mRNA more stable and increased PTCHD4 production. MeRIP RT-qPCR and eCLIP analyses were used to map this m6A modification to the last exon of PTCHD4 mRNA. Further investigation identified IGF2BP1, but not other m6A readers, as responsible for the stabilization and increased abundance of m6A-modified PTCHD4 mRNA. Silencing PTCHD4, a transmembrane protein, enhanced growth arrest and DNA damage in pre-senescent cells and sensitized them to senolysis and apoptosis. Our results indicate that m6A modification of PTCHD4 mRNA increases the production of PTCHD4, a protein associated with senescent cell survival, supporting the notion that regulating m6A modification on specific mRNAs could be exploited to eliminate senescent cells for therapeutic benefit.

LncRNA OIP5-AS1-directed miR-7 degradation promotes MYMX production during human myogenesis.

Long noncoding RNAs (lncRNAs) and microRNAs (miRNAs) modulate gene expression programs in physiology and disease. Here, we report a noncoding RNA regulatory network that modulates myoblast fusion into multinucleated myotubes, a process that occurs during muscle development and muscle regeneration after injury. In early stages of human myogenesis, the levels of lncRNA OIP5-AS1 increased, while the levels of miR-7 decreased. Moreover, OIP5-AS1 bound and induced miR-7 decay via target RNA-directed miRNA decay; accordingly, loss of OIP5-AS1 attenuated, while antagonizing miR-7 accelerated, myotube formation. We found that the OIP5-AS1-mediated miR-7 degradation promoted myoblast fusion, as it derepressed the miR-7 target MYMX mRNA, which encodes the fusogenic protein myomixer (MYMX). Remarkably, an oligonucleotide site blocker interfered with the OIP5-AS1-directed miR-7 degradation, allowing miR-7 to accumulate, lowering MYMX production and suppressing myotube formation. These results highlight a mechanism whereby lncRNA OIP5-AS1-mediated miR-7 decay promotes myotube formation by stimulating a myogenic fusion program.

Enhanced myogenesis through lncFAM-mediated recruitment of HNRNPL to the MYBPC2 promoter.

The mammalian transcriptome comprises a vast family of long noncoding (lnc)RNAs implicated in physiologic processes such as myogenesis, through which muscle forms during embryonic development and regenerates in the adult. However, the specific molecular mechanisms by which lncRNAs regulate human myogenesis are poorly understood. Here, we identified a novel muscle-specific lncRNA, lncFAM71E1-2:2 (lncFAM), which increased robustly during early human myogenesis. Overexpression of lncFAM promoted differentiation of human myoblasts into myotubes, while silencing lncFAM suppressed this process. As lncFAM resides in the nucleus, chromatin isolation by RNA purification followed by mass spectrometry (ChIRP-MS) analysis was employed to identify the molecular mechanisms whereby it might promote myogenesis. Analysis of lncFAM-interacting proteins revealed that lncFAM recruited the RNA-binding protein HNRNPL to the promoter of MYBPC2, in turn increasing MYBPC2 mRNA transcription and enhancing production of the myogenic protein MYBPC2. These results highlight a mechanism whereby a novel ribonucleoprotein complex, lncFAM-HNRNPL, elevates MYBPC2 expression transcriptionally to promote myogenesis.

HuR and GRSF1 modulate the nuclear export and mitochondrial localization of the lncRNA RMRP.

Some mitochondrial long noncoding RNAs (lncRNAs) are encoded by nuclear DNA, but the mechanisms that mediate their transport to mitochondria are poorly characterized. Using affinity RNA pull-down followed by mass spectrometry analysis, we found two RNA-binding proteins (RBPs), HuR (human antigen R) and GRSF1 (G-rich RNA sequence-binding factor 1), that associated with the nuclear DNA-encoded lncRNA RMRP and mobilized it to mitochondria. In cultured human cells, HuR bound RMRP in the nucleus and mediated its CRM1 (chromosome region maintenance 1)-dependent export to the cytosol. After RMRP was imported into mitochondria, GRSF1 bound RMRP and increased its abundance in the matrix. Loss of GRSF1 lowered the mitochondrial levels of RMRP, in turn suppressing oxygen consumption rates and modestly reducing mitochondrial DNA replication priming. Our findings indicate that RBPs HuR and GRSF1 govern the cytoplasmic and mitochondrial localization of the lncRNA RMRP, which is encoded by nuclear DNA but has key functions in mitochondria.

Ribosome profiling analysis of human skeletal muscle identifies reduced translation of mitochondrial proteins with age.

With advancing age, human muscle loses strength and function, but the molecular causes of these losses are unknown. Skeletal muscle shows an age-dependent decline in the levels of different proteins, but whether such decline is associated with reduced translation has not been studied. To address this gap of knowledge, we used the technique of ribosome profiling to study translation in muscle from middle-aged and old individuals. Using ribosome occupancy as a measure of translation status, several mRNAs showed differential translation with age. Older age was associated with lower translation of myosin and titin isoforms and more broadly with the translation of proteins involved in oxidative phosphorylation encoded by the mitochondrial genome. Based on our findings, we propose that mitochondrial proteins are less translated in old skeletal muscle.

Transcriptome signature of cellular senescence.

Cellular senescence, an integral component of aging and cancer, arises in response to diverse triggers, including telomere attrition, macromolecular damage and signaling from activated oncogenes. At present, senescent cells are identified by the combined presence of multiple traits, such as senescence-associated protein expression and secretion, DNA damage and ß-galactosidase activity; unfortunately, these traits are neither exclusively nor universally present in senescent cells. To identify robust shared markers of senescence, we have performed RNA-sequencing analysis across eight diverse models of senescence triggered in human diploid fibroblasts (WI-38, IMR-90) and endothelial cells (HUVEC, HAEC) by replicative exhaustion, exposure to ionizing radiation or doxorubicin, and expression of the oncogene HRASG12V. The intersection of the altered transcriptomes revealed 50 RNAs consistently elevated and 18 RNAs consistently reduced across all senescence models, including many protein-coding mRNAs and some non-coding RNAs. We propose that these shared transcriptome profiles will enable the identification of senescent cells in vivo, the investigation of their roles in aging and malignancy and the development of strategies to target senescent cells therapeutically.

Variation in human chromosome 21 ribosomal RNA genes characterized by TAR cloning and long-read sequencing.

Despite the key role of the human ribosome in protein biosynthesis, little is known about the extent of sequence variation in ribosomal DNA (rDNA) or its pre-rRNA and rRNA products. We recovered ribosomal DNA segments from a single human chromosome 21 using transformation-associated recombination (TAR) cloning in yeast. Accurate long-read sequencing of 13 isolates covering ∼0.82 Mb of the chromosome 21 rDNA complement revealed substantial variation among tandem repeat rDNA copies, several palindromic structures and potential errors in the previous reference sequence. These clones revealed 101 variant positions in the 45S transcription unit and 235 in the intergenic spacer sequence. Approximately 60% of the 45S variants were confirmed in independent whole-genome or RNA-seq data, with 47 of these further observed in mature 18S/28S rRNA sequences. TAR cloning and long-read sequencing enabled the accurate reconstruction of multiple rDNA units and a new, high-quality 44 838 bp rDNA reference sequence, which we have annotated with variants detected from chromosome 21 of a single individual. The large number of variants observed reveal heterogeneity in human rDNA, opening up the possibility of corresponding variations in ribosome dynamics.

Molecular dynamics of Dkk4 modulates Wnt action and regulates meibomian gland development.

Secreted Dickkopf (Dkk) proteins are major Wnt pathway modulators during organ development. Dkk1 has been widely studied and acts as a general Wnt inhibitor. However, the molecular function of other Dkks remains largely unknown. Here, we show that Dkk4 selectively inhibits a subset of Wnts, but is further inactivated by proteolytic cleavage. Meibomian gland (MG) formation is employed as a model where Dkk4 and its Wnt targets are expressed. Skin-specific expression of Dkk4 arrests MG growth at early germ phase, which is similar to that observed in Eda-ablated Tabby mice. Consistent with transient Dkk4 action, intact Dkk4 inhibits MG extension but the cleaved form progressively increases during MG development with a concomitant upswing in Wnt activity. Furthermore, both Dkk4 and its receptor (and Wnt co-receptor) Lrp6 are direct Eda targets during MG induction. In cell and organotypic cultures, Dkk4 inhibition is eliminated by elevation of Lrp6. Also, Lrp6 upregulation restores MG formation in Tabby mice. Thus, the dynamic state of Dkk4 itself and its interaction with Lrp6 modulates Wnt function during MG development, with a novel limitation of Dkk4 action by proteolytic cleavage.

High-purity circular RNA isolation method (RPAD) reveals vast collection of intronic circRNAs.

High-throughput RNA sequencing methods coupled with specialized bioinformatic analyses have recently uncovered tens of thousands of unique circular (circ)RNAs, but their complete sequences, genes of origin and functions are largely unknown. Given that circRNAs lack free ends and are thus relatively stable, their association with microRNAs (miRNAs) and RNA-binding proteins (RBPs) can influence gene expression programs. While exoribonuclease treatment is widely used to degrade linear RNAs and enrich circRNAs in RNA samples, it does not efficiently eliminate all linear RNAs. Here, we describe a novel method for the isolation of highly pure circRNA populations involving RNase R treatment followed by Polyadenylation and poly(A)+ RNA Depletion (RPAD), which removes linear RNA to near completion. High-throughput sequencing of RNA prepared using RPAD from human cervical carcinoma HeLa cells and mouse C2C12 myoblasts led to two surprising discoveries: (i) many exonic circRNA (EcircRNA) isoforms share an identical backsplice sequence but have different body sizes and sequences, and (ii) thousands of novel intronic circular RNAs (IcircRNAs) are expressed in cells. In sum, isolating high-purity circRNAs using the RPAD method can enable quantitative and qualitative analyses of circRNA types and sequence composition, paving the way for the elucidation of circRNA functions.

Involvement of Wnt, Eda and Shh at defined stages of sweat gland development.

To maintain body temperature, sweat glands develop from embryonic ectoderm by a poorly defined mechanism. We demonstrate a temporal cascade of regulation during mouse sweat gland formation. Sweat gland induction failed completely when canonical Wnt signaling was blocked in skin epithelium, and was accompanied by sharp downregulation of downstream Wnt, Eda and Shh pathway genes. The Wnt antagonist Dkk4 appeared to inhibit this induction: Dkk4 was sharply downregulated in ß-catenin-ablated mice, indicating that it is induced by Wnt/ß-catenin; however, its overexpression repressed Wnt target genes and significantly reduced gland numbers. Eda signaling succeeded Wnt. Wnt signaling was still active and nascent sweat gland pre-germs were still seen in Eda-null mice, but the pre-germs failed to develop further and the downstream Shh pathway was not activated. When Wnt and Eda were intact but Shh was ablated, germ induction and subsequent duct formation occurred normally, but the final stage of secretory coil formation failed. Thus, sweat gland development shows a relay of regulatory steps initiated by Wnt/ß-catenin - itself modulated by Dkk4 - with subsequent participation of Eda and Shh pathways.

Forkhead transcription factor FoxA1 regulates sweat secretion through Bestrophin 2 anion channel and Na-K-Cl cotransporter 1.

Body temperature is maintained in a narrow range in mammals, primarily controlled by sweating. In humans, the dynamic thermoregulatory organ, comprised of 2-4 million sweat glands distributed over the body, can secrete up to 4 L of sweat per day, thereby making it possible to withstand high temperatures and endure prolonged physical stress (e.g., long-distance running). The genetic basis for sweat gland function, however, is largely unknown. We find that the forkhead transcription factor, FoxA1, is required to generate mouse sweating capacity. Despite continued sweat gland morphogenesis, ablation of FoxA1 in mice results in absolute anihidrosis (lack of sweating). This inability to sweat is accompanied by down-regulation of the Na-K-Cl cotransporter 1 (Nkcc1) and the Ca(2+)-activated anion channel Bestrophin 2 (Best2), as well as glycoprotein accumulation in gland lumens and ducts. Furthermore, Best2-deficient mice display comparable anhidrosis and glycoprotein accumulation. These findings link earlier observations that both sodium/potassium/chloride exchange and Ca(2+) are required for sweat production. FoxA1 is inferred to regulate two corresponding features of sweat secretion. One feature, via Best2, catalyzes a bicarbonate gradient that could help to drive calcium-associated ionic transport; the other, requiring Nkcc1, facilitates monovalent ion exchange into sweat. These mechanistic components can be pharmaceutical targets to defend against hyperthermia and alleviate defective thermoregulation in the elderly, and may provide a model relevant to more complex secretory processes.

Characterization of age-associated gene expression changes in mouse sweat glands.

Evaporation of sweat on the skin surface is the major mechanism for dissipating heat in humans. The secretory capacity of sweat glands (SWGs) declines during aging, leading to heat intolerance in the elderly, but the mechanisms responsible for this decline are poorly understood. We investigated the molecular changes accompanying SWG aging in mice, where sweat tests confirmed a significant reduction of active SWGs in old mice relative to young mice. We first identified SWG-enriched mRNAs by comparing the skin transcriptome of Eda mutant Tabby male mice, which lack SWGs, with that of wild-type control mice by RNA-sequencing analysis. This comparison revealed 171 mRNAs enriched in SWGs, including 47 mRNAs encoding 'core secretory' proteins such as transcription factors, ion channels, ion transporters, and trans-synaptic signaling proteins. Among these, 28 SWG-enriched mRNAs showed significantly altered abundance in the aged male footpad skin, and 11 of them, including Foxa1, Best2, Chrm3, and Foxc1 mRNAs, were found in the 'core secretory' category. Consistent with the changes in mRNA expression levels, immunohistology revealed that higher numbers of secretory cells from old SWGs express the transcription factor FOXC1, the protein product of Foxc1 mRNA. In sum, our study identified mRNAs enriched in SWGs, including those that encode core secretory proteins, and altered abundance of these mRNAs and proteins with aging in mouse SWGs.

Silencing or amplification of endocannabinoid signaling in blastocysts via CB1 compromises trophoblast cell migration.

Endocannabinoid signaling plays key roles in multiple female reproductive events. Previous studies have shown an interesting phenomenon, that mice with either silenced or elevated endocannabinoid signaling via Cnr1 encoding CB(1) show similar defects in several pregnancy events, including preimplantation embryo development. To unravel the downstream signaling of this phenomenon, microarray studies were performed using RNAs collected from WT, Cnr1(-/-), and Faah(-/-) mouse blastocysts on day 4 of pregnancy. The results indicate that about 100 genes show unidirectional changes under either silenced or elevated anandamide signaling via CB(1). Functional enrichment analysis of the microarray data predicted that multiple biological functions and pathways are affected under aberrant endocannabinoid signaling. Among them, genes enriched in cell migration are suppressed in Cnr1(-/-) or Faah(-/-) blastocysts. Cell migration assays validated the prediction of functional enrichment analysis that cell mobility and spreading of either Cnr1(-/-) or Faah(-/-) trophoblast stem cells are compromised. Either silenced or elevated endocannabinoid signaling via CB(1) causes similar changes in downstream targets in preimplantation embryos and trophoblast stem cells. This study provides evidence that a tightly regulated endocannabinoid signaling is critical to normal preimplantation embryo development and migration of trophoblast stem cells.

Single-cell transcriptomic analysis uncovers diverse and dynamic senescent cell populations.

Senescence is a state of enduring growth arrest triggered by sublethal cell damage. Given that senescent cells actively secrete proinflammatory and matrix-remodeling proteins, their accumulation in tissues of older persons has been linked to many diseases of aging. Despite intense interest in identifying robust markers of senescence, the highly heterogeneous and dynamic nature of the senescent phenotype has made this task difficult. Here, we set out to comprehensively analyze the senescent transcriptome of human diploid fibroblasts at the individual-cell scale by performing single-cell RNA-sequencing analysis through two approaches. First, we characterized the different cell states in cultures undergoing senescence triggered by different stresses, and found distinct cell subpopulations that expressed mRNAs encoding proteins with roles in growth arrest, survival, and the secretory phenotype. Second, we characterized the dynamic changes in the transcriptomes of cells as they developed etoposide-induced senescence; by tracking cell transitions across this process, we found two different senescence programs that developed divergently, one in which cells expressed traditional senescence markers such as p16 (CDKN2A) mRNA, and another in which cells expressed long noncoding RNAs and splicing was dysregulated. Finally, we obtained evidence that the proliferation status at the time of senescence initiation affected the path of senescence, as determined based on the expressed RNAs. We propose that a deeper understanding of the transcriptomes during the progression of different senescent cell phenotypes will help develop more effective interventions directed at this detrimental cell population.

Senescence suppresses the integrated stress response and activates a stress-enhanced secretory phenotype.

Senescence is a state of indefinite cell cycle arrest associated with aging, cancer, and age-related diseases. Here, using label-based mass spectrometry, ribosome profiling and nanopore direct RNA sequencing, we explore the coordinated interaction of translational and transcriptional programs of human cellular senescence. We find that translational deregulation and a corresponding maladaptive integrated stress response (ISR) is a hallmark of senescence that desensitizes senescent cells to stress. We present evidence that senescent cells maintain high levels of eIF2α phosphorylation, typical of ISR activation, but translationally repress production of the stress response transcription factor 4 (ATF4) by ineffective bypass of the inhibitory upstream open reading frames. Surprisingly, ATF4 translation remains inhibited even after acute proteotoxic and amino acid starvation stressors, resulting in a highly diminished stress response. Furthermore, absent a response, stress augments the senescence secretory phenotype, thus intensifying a proinflammatory state that exacerbates disease. Our results reveal a novel mechanism that senescent cells exploit to evade an adaptive stress response and remain viable.

The YAP-TEAD complex promotes senescent cell survival by lowering endoplasmic reticulum stress.

Sublethal cell damage can trigger senescence, a complex adaptive program characterized by growth arrest, resistance to apoptosis and a senescence-associated secretory phenotype (SASP). Here, a whole-genome CRISPR knockout screen revealed that proteins in the YAP-TEAD pathway influenced senescent cell viability. Accordingly, treating senescent cells with a drug that inhibited this pathway, verteporfin (VPF), selectively triggered apoptotic cell death largely by derepressing DDIT4, which in turn inhibited mTOR. Reducing mTOR function in senescent cells diminished endoplasmic reticulum (ER) biogenesis, triggering ER stress and apoptosis due to high demands on ER function by the SASP. Importantly, VPF treatment decreased the numbers of senescent cells in the organs of old mice and mice exhibiting doxorubicin-induced senescence. Moreover, VPF treatment reduced immune cell infiltration and pro-fibrotic transforming growth factor-ß signaling in aging mouse lungs, improving tissue homeostasis. We present an alternative senolytic strategy that eliminates senescent cells by hindering ER activity required for SASP production.

Pleiotropic effects of BAFF on the senescence-associated secretome and growth arrest.

Senescent cells release a variety of cytokines, proteases, and growth factors collectively known as the senescence-associated secretory phenotype (SASP). Sustained SASP contributes to a pattern of chronic inflammation associated with aging and implicated in many age-related diseases. Here, we investigated the expression and function of the immunomodulatory cytokine BAFF (B-cell activating factor; encoded by the TNFSF13B gene), a SASP protein, in multiple senescence models. We first characterized BAFF production across different senescence paradigms, including senescent human diploid fibroblasts (WI-38, IMR-90) and monocytic leukemia cells (THP-1), and tissues of mice induced to undergo senescence. We then identified IRF1 (interferon regulatory factor 1) as a transcription factor required for promoting TNFSF13B mRNA transcription in senescence. We discovered that suppressing BAFF production decreased the senescent phenotype of both fibroblasts and monocyte-like cells, reducing IL6 secretion and SA-ß-Gal staining. Importantly, however, the influence of BAFF on the senescence program was cell type-specific: in monocytes, BAFF promoted the early activation of NF-κB and general SASP secretion, while in fibroblasts, BAFF contributed to the production and function of TP53 (p53). We propose that BAFF is elevated across senescence models and is a potential target for senotherapy.