m6A-Related lncRNAs Predict Overall Survival of Patients and Regulate the Tumor Immune Microenvironment in Osteosarcoma.

Background: m6A-related lncRNAs have demonstrated great potential tumor diagnostic and therapeutic targets. The goal of this work was to find m6A-regulated lncRNAs in osteosarcoma patients. Method: The Cancer Genome Atlas (TCGA) database was used to retrieve RNA sequencing and medical information from osteosarcoma sufferers. The Pearson's correlation test was used to identify the m6A-related lncRNAs. A risk model was built using univariate and multivariable Cox regression analysis. Kaplan-Meier survival analysis and receiver functional requirements were used to assess the risk model's performance (ROC). By using the CIBERSORT method, the associations between the relative risks and different immune cell infiltration were investigated. Lastly, the bioactivities of high-risk and low-risk subgroups were investigated using Gene Set Enrichment Analysis (GSEA). Result: A total of 531 m6A-related lncRNAs were obtained from TCGA. Seven lncRNAs have demonstrated prognostic values. A total of 88 OS patients were separated into cluster 1, cluster 2, and cluster 3. The overall survival rate of OS patients in cluster 3 was more favorable than that of those in cluster 1 and cluster 2. The average Stromal score was much higher in cluster 1 than in cluster 2 and cluster 3 (P < 0.05). The expression levels of lncRNAs used in the construction of the risk prediction model in the high-risk group were generally lower than those in the low-risk group. Analysis of patient survival indicated that the survival of the low-risk group was higher than that of the high-risk group (P < 0.0001) and the area under the curve (AUC) of the ROC curve was 0.719. Using the CIBERSORT algorithm, the results revealed that Macrophages M0, Macrophages M2, and T cells CD4 memory resting accounted for a large proportion of immune cell infiltration. By GSEA analysis, our results implied that the high-risk group was mainly involved in unfolded protein response, DNA repair signaling, and epithelial-mesenchymal transition signaling pathway and glycolysis pathway; meanwhile, the low-risk group was mainly involved in estrogen response early and KRAS signaling pathway. Conclusion: Our investigation showed that m6A-related lncRNAs remained tightly connected to the immunological microenvironment of osteosarcoma tumors, potentially influencing carcinogenesis and development. The immune microenvironment and immune-related biochemical pathways can be changed by regulating the transcription of M6A modulators or lncRNAs. In addition, we looked for risk-related signaling of m6A-related lncRNAs in osteosarcomas and built and validated the risk prediction system. The findings of our current analysis will facilitate the assessment of outcomes and the development of immunotherapies for sufferers of osteosarcomas.

Single-cell analyses identify dysfunctional CD16+ CD8 T cells in smokers.

Tobacco smoke exposure contributes to the global burden of communicable and chronic diseases. To identify immune cells affected by smoking, we use single-cell RNA sequencing on peripheral blood from smokers and nonsmokers. Transcriptomes reveal a subpopulation of FCGR3A (CD16)-expressing Natural Killer (NK)-like CD8 T lymphocytes that increase in smokers. Mass cytometry confirms elevated CD16+ CD8 T cells in smokers. Inferred as highly differentiated by pseudotime analysis, NK-like CD8 T cells express markers characteristic of effector memory re-expressing CD45RA T (TEMRA) cells. Indicative of immune aging, smokers' CD8 T cells are biased toward differentiated cells and smokers have fewer naïve cells than nonsmokers. DNA methylation-based models show that smoking dose is associated with accelerated aging and decreased telomere length, a biomarker of T cell senescence. Immune aging accompanies T cell senescence, which can ultimately lead to impaired immune function. This suggests a role for smoking-induced, senescence-associated immune dysregulation in smoking-mediated pathologies.

Identification of Smoking-Associated Differentially Methylated Regions Using Reduced Representation Bisulfite Sequencing and Cell type-Specific Enhancer Activation and Gene Expression.

BACKGROUND: Cigarette smoke is a causal factor in cancers and cardiovascular disease. Smoking-associated differentially methylated regions (SM-DMRs) have been observed in disease studies, but the causal link between altered DNA methylation and transcriptional change is obscure. OBJECTIVE: Our objectives were to finely resolve SM-DMRs and to interrogate the mechanistic link between SM-DMRs and altered transcription of enhancer noncoding RNA (eRNA) and mRNA in human circulating monocytes. METHOD: We integrated SM-DMRs identified by reduced representation bisulfite sequencing (RRBS) of circulating CD14+ monocyte DNA collected from two independent human studies [n=38 from Clinical Research Unit (CRU) and n=55 from the Multi-Ethnic Study of Atherosclerosis (MESA), about half of whom were active smokers] with gene expression for protein-coding genes and noncoding RNAs measured by RT-PCR or RNA sequencing. Candidate SM-DMRs were compared with RRBS of purified CD4+ T cells, CD8+ T cells, CD15+ granulocytes, CD19+ B cells, and CD56+ NK cells (n=19 females, CRU). DMRs were validated using pyrosequencing or bisulfite amplicon sequencing in up to 85 CRU volunteers, who also provided saliva DNA. RESULTS: RRBS identified monocyte SM-DMRs frequently located in putative gene regulatory regions. The most significant monocyte DMR occurred at a poised enhancer in the aryl-hydrocarbon receptor repressor gene (AHRR) and it was also detected in both granulocytes and saliva DNA. To our knowledge, we identify for the first time that SM-DMRs in or near AHRR, C5orf55-EXOC-AS, and SASH1 were associated with increased noncoding eRNA as well as mRNA in monocytes. Functionally, the AHRR SM-DMR appeared to up-regulate AHRR mRNA through activating the AHRR enhancer, as suggested by increased eRNA in the monocytes, but not granulocytes, from smokers compared with nonsmokers. CONCLUSIONS: Our findings suggest that AHRR SM-DMR up-regulates AHRR mRNA in a monocyte-specific manner by activating the AHRR enhancer. Cell type-specific activation of enhancers at SM-DMRs may represent a mechanism driving smoking-related disease. https://doi.org/10.1289/EHP2395.

Tobacco exposure-related alterations in DNA methylation and gene expression in human monocytes: the Multi-Ethnic Study of Atherosclerosis (MESA).

Alterations in DNA methylation and gene expression in blood leukocytes are potential biomarkers of harm and mediators of the deleterious effects of tobacco exposure. However, methodological issues, including the use of self-reported smoking status and mixed cell types have made previously identified alterations in DNA methylation and gene expression difficult to interpret. In this study, we examined associations of tobacco exposure with DNA methylation and gene expression, utilizing a biomarker of tobacco exposure (urine cotinine) and CD14+ purified monocyte samples from 934 participants of the community-based Multi-Ethnic Study of Atherosclerosis (MESA). Urine cotinine levels were measured using an immunoassay. DNA methylation and gene expression were measured with microarrays. Multivariate linear regression was used to test for associations adjusting for age, sex, race/ethnicity, education, and study site. Urine cotinine levels were associated with methylation of 176 CpGs [false discovery rate (FDR)<0.01]. Four CpGs not previously identified by studies of non-purified blood samples nominally replicated (P value<0.05) with plasma cotinine-associated methylation in 128 independent monocyte samples. Urine cotinine levels associated with expression of 12 genes (FDR<0.01), including increased expression of P2RY6 (Beta ± standard error = 0.078 ± 0.008, P = 1.99 × 10-22), a gene previously identified to be involved in the release of pro-inflammatory cytokines. No cotinine-associated (FDR<0.01) methylation profiles significantly (FDR<0.01) correlated with cotinine-associated (FDR<0.01) gene expression profiles. In conclusion, our findings i) identify potential monocyte-specific smoking-associated methylation patterns and ii) suggest that alterations in methylation may not be a main mechanism regulating gene expression in monocytes in response to cigarette smoking.

Correction: Distinct Epigenetic Effects of Tobacco Smoking in Whole Blood and among Leukocyte Subtypes.

[This corrects the article DOI: 10.1371/journal.pone.0166486.].

Distinct Epigenetic Effects of Tobacco Smoking in Whole Blood and among Leukocyte Subtypes.

Tobacco smoke exposure dramatically alters DNA methylation in blood cells and may mediate smoking-associated complex diseases through effects on immune cell function. However, knowledge of smoking effects in specific leukocyte subtypes is limited. To better characterize smoking-associated methylation changes in whole blood and leukocyte subtypes, we used Illumina 450K arrays and Reduced Representation Bisulfite Sequencing (RRBS) to assess genome-wide DNA methylation. Differential methylation analysis in whole blood DNA from 172 smokers and 81 nonsmokers revealed 738 CpGs, including 616 previously unreported CpGs, genome-wide significantly associated with current smoking (p <1.2x10-7, Bonferroni correction). Several CpGs (MTSS1, NKX6-2, BTG2) were associated with smoking duration among heavy smokers (>22 cigarettes/day, n = 86) which might relate to long-term heavy-smoking pathology. In purified leukocyte subtypes from an independent group of 20 smokers and 14 nonsmokers we further examined methylation and gene expression for selected genes among CD14+ monocytes, CD15+ granulocytes, CD19+ B cells, and CD2+ T cells. In 10 smokers and 10 nonsmokers we used RRBS to fine map differential methylation in CD4+ T cells, CD8+ T cells, CD14+, CD15+, CD19+, and CD56+ natural killer cells. Distinct cell-type differences in smoking-associated methylation and gene expression were identified. AHRR (cg05575921), ALPPL2 (cg21566642), GFI1 (cg09935388), IER3 (cg06126421) and F2RL3 (cg03636183) showed a distinct pattern of significant smoking-associated methylation differences across cell types: granulocytes> monocytes>> B cells. In contrast GPR15 (cg19859270) was highly significant in T and B cells and ITGAL (cg09099830) significant only in T cells. Numerous other CpGs displayed distinctive cell-type responses to tobacco smoke exposure that were not apparent in whole blood DNA. Assessing the overlap between these CpG sites and differential methylated regions (DMRs) with RRBS in 6 cell types, we confirmed cell-type specificity in the context of DMRs. We identified new CpGs associated with current smoking, pack-years, duration, and revealed unique profiles of smoking-associated DNA methylation and gene expression among immune cell types, providing potential clues to hematopoietic lineage-specific effects in disease etiology.

A Polymorphic Antioxidant Response Element Links NRF2/sMAF Binding to Enhanced MAPT Expression and Reduced Risk of Parkinsonian Disorders.

The NRF2/sMAF protein complex regulates the oxidative stress response by occupying cis-acting enhancers containing an antioxidant response element (ARE). Integrating genome-wide maps of NRF2/sMAF occupancy with disease-susceptibility loci, we discovered eight polymorphic AREs linked to 14 highly ranked disease-risk SNPs in individuals of European ancestry. Among these SNPs was rs242561, located within a regulatory region of the MAPT gene (encoding microtubule-associated protein Tau). It was consistently occupied by NRF2/sMAF in multiple experiments and its strong-binding allele associated with higher mRNA levels in cell lines and human brain tissue. Induction of MAPT transcription by NRF2 was confirmed using a human neuroblastoma cell line and a Nrf2-deficient mouse model. Most importantly, rs242561 displayed complete linkage disequilibrium with a highly protective allele identified in multiple GWASs of progressive supranuclear palsy, Parkinson's disease, and corticobasal degeneration. These observations suggest a potential role for NRF2/sMAF in tauopathies and a possible role for NRF2 pathway activators in disease prevention.

DNA Methylation of the Aryl Hydrocarbon Receptor Repressor Associations With Cigarette Smoking and Subclinical Atherosclerosis.

BACKGROUND: Tobacco smoke contains numerous agonists of the aryl hydrocarbon receptor (AhR) pathway, and activation of the AhR pathway was shown to promote atherosclerosis in mice. Intriguingly, cigarette smoking is most strongly and robustly associated with DNA modifications to an AhR pathway gene, the AhR repressor (AHRR). We hypothesized that altered AHRR methylation in monocytes, a cell type sensitive to cigarette smoking and involved in atherogenesis, may be a part of the biological link between cigarette smoking and atherosclerosis. METHODS AND RESULTS: DNA methylation profiles of AHRR in monocytes (542 CpG sites ± 150 kb of AHRR, using Illumina 450K array) were integrated with smoking habits and ultrasound-measured carotid plaque scores from 1256 participants of the Multi-Ethnic Study of Atherosclerosis (MESA). Methylation of cg05575921 significantly associated (P=6.1 × 10(-134)) with smoking status (current versus never). Novel associations between cg05575921 methylation and carotid plaque scores (P=3.1 × 10(-10)) were identified, which remained significant in current and former smokers even after adjusting for self-reported smoking habits, urinary cotinine, and well-known cardiovascular disease risk factors. This association replicated in an independent cohort using hepatic DNA (n=141). Functionally, cg05575921 was located in a predicted gene expression regulatory element (enhancer) and had methylation correlated with AHRR mRNA profiles (P=1.4 × 10(-17)) obtained from RNA sequencing conducted on a subset (n=373) of the samples. CONCLUSIONS: These findings suggest that AHRR methylation may be functionally related to AHRR expression in monocytes and represents a potential biomarker of subclinical atherosclerosis in smokers.

Ten-eleven translocation 1 (Tet1) is regulated by O-linked N-acetylglucosamine transferase (Ogt) for target gene repression in mouse embryonic stem cells.

As a member of the Tet (Ten-eleven translocation) family proteins that can convert 5-methylcytosine (5mC) to 5-hydroxylmethylcytosine (5hmC), Tet1 has been implicated in regulating global DNA demethylation and gene expression. Tet1 is highly expressed in embryonic stem (ES) cells and appears primarily to repress developmental genes for maintaining pluripotency. To understand how Tet1 may regulate gene expression, we conducted large scale immunoprecipitation followed by mass spectrometry of endogenous Tet1 in mouse ES cells. We found that Tet1 could interact with multiple chromatin regulators, including Sin3A and NuRD complexes. In addition, we showed that Tet1 could also interact with the O-GlcNAc transferase (Ogt) and be O-GlcNAcylated. Depletion of Ogt led to reduced Tet1 and 5hmC levels on Tet1-target genes, whereas ectopic expression of wild-type but not enzymatically inactive Ogt increased Tet1 levels. Mutation of the putative O-GlcNAcylation site on Tet1 led to decreased O-GlcNAcylation and level of the Tet1 protein. Our results suggest that O-GlcNAcylation can positively regulate Tet1 protein concentration and indicate that Tet1-mediated 5hmC modification and target repression is controlled by Ogt.

Kdm2b maintains murine embryonic stem cell status by recruiting PRC1 complex to CpG islands of developmental genes.

Polycomb group (PcG) proteins play important roles in repressing lineage-specific genes and maintaining the undifferentiated state of mouse embryonic stem cells (mESCs). However, how PcG proteins are recruited to their target genes is largely unknown. Here, we show that the H3K36-specific histone demethylase Kdm2b is highly expressed in mESCs and regulated by the pluripotent factors Oct4 and Sox2 directly. Depletion of Kdm2b in mESCs causes de-repression of lineage-specific genes and induces early differentiation. The function of Kdm2b depends on its CxxC-ZF domain, which mediates its genome-wide binding to CpG islands (CGIs). Kdm2b interacts with the core components of polycomb repressive complex 1 (PRC1) and recruits the complex to the CGIs of early lineage-specific genes. Thus, our study not only reveals an Oct4-Sox2-Kdm2b-PRC1-CGI regulatory axis and its function in maintaining the undifferentiated state of mESCs, but also demonstrates a critical function of Kdm2b in recruiting PRC1 to the CGIs of lineage-specific genes to repress their expression.

Telomeres-structure, function, and regulation.

In mammals, maintenance of the linear chromosome ends (or telomeres) involves faithful replication of genetic materials and protection against DNA damage signals, to ensure genome stability and integrity. These tasks are carried out by the telomerase holoenzyme and a unique nucleoprotein structure in which an array of telomere-associated proteins bind to telomeric DNA to form special protein/DNA complexes. The telomerase complex, which is comprised of telomeric reverse transcriptase (TERT), telomeric RNA component (TERC), and other assistant factors, is responsible for adding telomeric repeats to the ends of chromosomes. Without proper telomere maintenance, telomere length will shorten with successive round of DNA replication due to the so-called end replication problem. Aberrant regulation of telomeric proteins and/or telomerase may lead to abnormalities that can result in diseases such as dyskeratosis congenita (DC) and cancers. Understanding the mechanisms that regulate telomere homeostasis and the factors that contribute to telomere dysfunction should aid us in developing diagnostic and therapeutic tools for these diseases.

The trithorax group protein Ash2l is essential for pluripotency and maintaining open chromatin in embryonic stem cells.

Embryonic stem (ES) cells exhibit general characteristics of open chromatin, a state that may be necessary for ES cells to efficiently self-renew while remaining poised for differentiation. Histone H3K4 and H3K9 trimethylation associate as a general rule, with open and silenced chromatin, respectively, for ES cell pluripotency maintenance. However, how histone modifications are regulated to maintain open chromatin in ES cells remains largely unknown. Here, we demonstrate that trithorax protein Ash2l, homologue of the Drosophila Ash2 (absent, small, homeotic-2) protein, is a key regulator of open chromatin in ES cells. Consistent with Ash2l being a core subunit of mixed lineage leukemia methyltransferase complex, RNAi knockdown of Ash2l was sufficient to reduce H3K4 methylation levels and drive ES cells to a silenced chromatin state with high H3K9 trimethylation. Genome-wide ChIP-seq analysis indicated that Ash2l is recruited to target loci through two distinct modes and enriched at a family of genes implicated in open chromatin regulation, including chromatin remodeler Cdh7, transcription factor c-Myc, and H3K9 demethylase Kdm4c. Our results underscore the importance of Ash2l in open chromatin regulation and provide insight into how the open chromatin landscape is maintained in ES cells.

SGK3 is associated with estrogen receptor expression in breast cancer.

While breast cancer mortality rate has seen a steady decline in the last few decades, advances in better treatment and diagnostic tools remain important as we come into the age of personalized therapy. In this report, we describe our studies of SGK3's role in breast cancer. SGK3 (also known as CISK) is a member of the AGC family of kinases. Our previous work indicates that SGK3 functions downstream of the PI 3-kinase cascade and shares molecular and biochemical similarities with Akt. Here, we show that SGK3 expression is linked to estrogen receptor (ER) both in breast caner cell lines and in primary tumor samples. Our analysis also indicated a positive correlation between SGK3 expression and tumor prognosis. Importantly, our immunochemistry analysis of human tumor samples established a clinical link between SGK3 expression and ER+ tumors. These findings implicate SGK3 as an additional component to a complex and heterogeneous disease, and point to the potential benefits of incorporating SGK3 into the process of breast cancer diagnosis and treatment.

Convergent transcription through a long CAG tract destabilizes repeats and induces apoptosis.

Short repetitive sequences are common in the human genome, and many fall within transcription units. We have previously shown that transcription through CAG repeat tracts destabilizes them in a way that depends on transcription-coupled nucleotide excision repair and mismatch repair. Recent observations that antisense transcription accompanies sense transcription in many human genes led us to test the effects of antisense transcription on triplet repeat instability in human cells. Here, we report that simultaneous sense and antisense transcription (convergent transcription) initiated from two inducible promoters flanking a CAG95 tract in a nonessential gene enhances repeat instability synergistically, arrests the cell cycle, and causes massive cell death via apoptosis. Using chemical inhibitors and small interfering RNA (siRNA) knockdowns, we identified the ATR (ataxia-telangiectasia mutated [ATM] and Rad3 related) signaling pathway as a key mediator of this cellular response. RNA polymerase II, replication protein A (RPA), and components of the ATR signaling pathway accumulate at convergently transcribed repeat tracts, accompanied by phosphorylation of ATR, CHK1, and p53. Cell death depends on simultaneous sense and antisense transcription and is proportional to their relative levels, it requires the presence of the repeat tract, and it occurs in both proliferating and nonproliferating cells. Convergent transcription through a CAG repeat represents a novel mechanism for triggering a cellular stress response, one that is initiated by events at a single locus in the genome and resembles the response to DNA damage.

OB fold-containing protein 1 (OBFC1), a human homolog of yeast Stn1, associates with TPP1 and is implicated in telomere length regulation.

The telosome/shelterin, a six-protein complex formed by TRF1, TRF2, RAP1, TIN2, POT1, and TPP1, functions as the core of the telomere interactome, acting as the molecular platform for the assembly of higher order complexes and coordinating cross-talks between various protein subcomplexes. Within the telosome, there are two oligonucleotide- or oligosaccharide-binding (OB) fold-containing proteins, TPP1 and POT1. They can form heterodimers that bind to the telomeric single-stranded DNA, an activity that is central for telomere end capping and telomerase recruitment. Through proteomic analyses, we found that in addition to POT1, TPP1 can associate with another OB fold-containing protein, OBFC1/AAF44. The yeast homolog of OBFC1 is Stn1, which plays a critical role in telomere regulation. We show here that OBFC1/AAF44 can localize to telomeres in human cells and bind to telomeric single-stranded DNA in vitro. Furthermore, overexpression of an OBFC1 mutant resulted in elongated telomeres in human cells, implicating OBFC1/AAF4 in telomere length regulation. Taken together, our studies suggest that OBFC1/AAF44 represents a new player in the telomere interactome for telomere maintenance.

TPP1 is a homologue of ciliate TEBP-beta and interacts with POT1 to recruit telomerase.

Telomere dysfunction may result in chromosomal abnormalities, DNA damage responses, and even cancer. Early studies in lower organisms have helped to establish the crucial role of telomerase and telomeric proteins in maintaining telomere length and protecting telomere ends. In Oxytricha nova, telomere G-overhangs are protected by the TEBP-alpha/beta heterodimer. Human telomeres contain duplex telomeric repeats with 3' single-stranded G-overhangs, and may fold into a t-loop structure that helps to shield them from being recognized as DNA breaks. Additionally, the TEBP-alpha homologue, POT1, which binds telomeric single-stranded DNA (ssDNA), associates with multiple telomeric proteins (for example, TPP1, TIN2, TRF1, TRF2 and RAP1) to form the six-protein telosome/shelterin and other subcomplexes. These telomeric protein complexes in turn interact with diverse pathways to form the telomere interactome for telomere maintenance. However, the mechanisms by which the POT1-containing telosome communicates with telomerase to regulate telomeres remain to be elucidated. Here we demonstrate that TPP1 is a putative mammalian homologue of TEBP-beta and contains a predicted amino-terminal oligonucleotide/oligosaccharide binding (OB) fold. TPP1-POT1 association enhanced POT1 affinity for telomeric ssDNA. In addition, the TPP1 OB fold, as well as POT1-TPP1 binding, seemed critical for POT1-mediated telomere-length control and telomere-end protection in human cells. Disruption of POT1-TPP1 interaction by dominant negative TPP1 expression or RNA interference (RNAi) resulted in telomere-length alteration and DNA damage responses. Furthermore, we offer evidence that TPP1 associates with the telomerase in a TPP1-OB-fold-dependent manner, providing a physical link between telomerase and the telosome/shelterin complex. Our findings highlight the critical role of TPP1 in telomere maintenance, and support a yin-yang model in which TPP1 and POT1 function as a unit to protect human telomeres, by both positively and negatively regulating telomerase access to telomere DNA.