Histone lysine demethylase KDM5B facilitates proliferation and suppresses apoptosis in human acute myeloid leukemia cells through the miR-140-3p/BCL2 axis.

The histone lysine demethylase KDM5B is frequently up-regulated in various human cancer cells. However, its expression and functional role in human acute myeloid leukemia (AML) cells remain unclear. Here, we found that the expression level of KDM5B is high in primary human AML cells. We have demonstrated that knocking down KDM5B leads to apoptosis and impairs proliferation in primary human AML and some human AML cell lines. We further identified miR-140-3p as a downstream target gene of KDM5B. KDM5B expression was inversely correlated with the miR-140-3p level in primary human AML cells. Molecular studies showed that silencing KDM5B enhanced H3K4 trimethylation (H3K4me3) at the promoter of miR-140-3p, leading to high expression of miR-140-3p, which in turn inhibited B-cell CLL/lymphoma 2 (BCL2) expression. Finally, we demonstrate that the defective proliferation induced by KDM5B knockdown (KD) can be rescued with the miR-140-3p inhibitor or enhanced by combining KDM5B KD with a BCL2 inhibitor. Altogether, our data support the conclusion that KDM5B promotes tumorigenesis in human AML cells through the miR-140-3p/BCL2 axis. Targeting the KDM5B/miR-140-3p/BCL2 pathway may hold therapeutic promise for treating human AML.

Directed Differentiation of Human Induced Pluripotent Stem Cells to Heart Valve Cells.

BACKGROUND: A main obstacle in current valvular heart disease research is the lack of high-quality homogeneous functional heart valve cells. Human induced pluripotent stem cells (hiPSCs)-derived heart valve cells may help with this dilemma. However, there are no well-established protocols to induce hiPSCs to differentiate into functional heart valve cells, and the networks that mediate the differentiation have not been fully elucidated. METHODS: To generate heart valve cells from hiPSCs, we sequentially activated the Wnt, BMP4, VEGF (vascular endothelial growth factor), and NFATc1 signaling pathways using CHIR-99021, BMP4, VEGF-165, and forskolin, respectively. The transcriptional and functional similarity of hiPSC-derived heart valve cells compared with primary heart valve cells were characterized. Longitudinal single-cell RNA sequencing was used to uncover the trajectory, switch genes, pathways, and transcription factors of the differentiation. RESULTS: An efficient protocol was developed to induce hiPSCs to differentiate into functional hiPSC-derived valve endothelial-like cells and hiPSC-derived valve interstitial-like cells. After 6-day differentiation and CD144 magnetic bead sorting, ≈70% CD144+ cells and 30% CD144- cells were obtained. On the basis of single-cell RNA sequencing data, the CD144+ cells and CD144- cells were found to be highly similar to primary heart valve endothelial cells and primary heart valve interstitial cells in gene expression profile. Furthermore, CD144+ cells had the typical function of primary heart valve endothelial cells, including tube formation, uptake of low-density lipoprotein, generation of endothelial nitric oxide synthase, and response to shear stress. Meanwhile, CD144- cells could secret collagen and matrix metalloproteinases, and differentiate into osteogenic or adipogenic lineages like primary heart valve interstitial cells. Therefore, we identified CD144+ cells and CD144- cells as hiPSC-derived valve endothelial-like cells and hiPSC-derived valve interstitial-like cells, respectively. Using single-cell RNA sequencing analysis, we demonstrated that the trajectory of heart valve cell differentiation was consistent with embryonic valve development. We identified the main switch genes (NOTCH1, HEY1, and MEF2C), signaling pathways (TGF-ß, Wnt, and NOTCH), and transcription factors (MSX1, SP5, and MECOM) that mediated the differentiation. Finally, we found that hiPSC-derived valve interstitial-like cells might derive from hiPSC-derived valve endothelial-like cells undergoing endocardial-mesenchymal transition. CONCLUSIONS: In summary, this is the first study to report an efficient strategy to generate functional hiPSC-derived valve endothelial-like cells and hiPSC-derived valve interstitial-like cells from hiPSCs, as well as to elucidate the differentiation trajectory and transcriptional dynamics of hiPSCs differentiated into heart valve cells.

Ppan is essential for preimplantation development in mice†.

PETER PAN (PPAN), located to nucleoli and mitochondria, is a member of the Brix domain protein family, involved in rRNA processing through its rRNA binding motif and mitochondrial apoptosis by protecting mitochondria structure and suppressing basal autophagic flux. Ppan is important for cell proliferation and viability, and mutation of Ppan in Drosophila caused larval lethality and oogenesis failure. Yet, its role in mammalian reproduction remains unclear. In this study, we explored the function of Ppan in oocyte maturation and early embryogenesis using conditional knockout mouse model. Deficiency of maternal Ppan significantly downregulated the expression level of 5.8S rRNA, 18S rRNA, and 28S rRNA, though it had no effect on oocyte maturation or preimplantation embryo development. However, depletion of both maternal and zygotic Ppan blocked embryonic development at morula stage. Similar phenotype was obtained when only zygotic Ppan was depleted. We further identified no DNA binding activity of PPAN in mouse embryonic stem cells, and depletion of Ppan had minimum impact on transcriptome but decreased expression of 5.8S rRNA, 18S rRNA, and 28S rRNA nevertheless. Our findings demonstrate that Ppan is indispensable for early embryogenesis in mice.

Regulatory roles of alternative splicing at Ezh2 gene in mouse oocytes.

BACKGROUND: Enhancer of zeste homologue 2 (EZH2), the core member of polycomb repressive complex 2 (PRC2), has multiple splicing modes and performs various physiological functions. However, function and mechanism of alternative splicing at Ezh2 exon 3 in reproduction are unknown. METHODS: We generated Ezh2Long and Ezh2Short mouse models with different point mutations at the Ezh2 exon 3 alternative splicing site, and each mutant mouse model expressed either the long or the short isoform of Ezh2. We examined mutant mouse fertility and oocyte development to assess the function of Ezh2 alternative splicing at exon 3 in the reproductive system. RESULTS: We found that Ezh2Long female mice had normal fertility. However, Ezh2Short female mice had significantly decreased fertility and obstructed oogenesis, with compromised mitochondrial function in Ezh2Short oocytes. Interestingly, increased EZH2 protein abundance and accumulated H3K27me3 were observed in Ezh2Short oocytes. CONCLUSIONS: Our results demonstrate that correct Ezh2 alternative splicing at exon 3 is important for mouse oogenesis.

Exposure of mouse oocytes to N,N-dimethylformamide impairs mitochondrial functions and reduces oocyte quality.

N,N-dimethylformamide (DMF) is a widely-used solvent for the synthesis of synthetic fibers such as polyacrylonitrile fiber, and can also be used to make medicine. Although this organic solvent has multipurpose applications, its biological toxicity cannot be ignored and its impact on mammalian reproduction remains largely unexplored. Our study found that DMF exposure inhibited oocyte maturation and fertilization ability. Transcriptomic analysis indicated that DMF exposure changed the expression of genes and transposable elements in oocytes. Subcellular structure examination found that DMF exposure caused mitochondrial dysfunction, abnormal aggregation of mitochondria and decreased mitochondrial membrane potential in mouse oocytes. Its exposure also caused abnormal distribution of Golgi apparatus and endoplasmic reticulum which formed large number of clusters. In addition, oxidative stress occurs in oocytes exposed to DMF, which was manifested by an increase in the level of reactive oxygen species. We found that DMF exposure induced disordered spindle and chromosomes abnormality. Meanwhile, we examined various histone modification levels in oocytes exposed to DMF and found that DMF exposure reduced H3K9me3, H3K9ac, H3K27ac, and H4K16ac levels in mouse oocytes. Moreover, DMF-treated oocytes failed to form pronuclei after fusion with normal sperm. Collectively, DMF exposure caused mitochondrial damage, oxidative stress, spindle assembly and chromosome arrangement disorder, leading to oocyte maturation arrest and fertilization failure.

Cleavage-embryo genes and transposable elements are regulated by histone variant H2A.X.

During mammalian preimplantation development, stimulation of zygotic genome activation (ZGA) and transposable elements (TEs) shapes totipotency profiling. A rare mouse embryonic stem cells (mESCs) subpopulation is capable of transiently entering a state resembling 2-cell stage embryos, with subtypes of TEs expressed and ZGA genes transiently activated. In this study, we found that deletion of H2A.X in mESCs led to a significant upregulation of ZGA genes and misregulated TEs. ChIP-seq analysis indicated a direct association of H2A.X at the Dux locus for silencing the Dux gene and its downstream ZGA genes in mESCs. We also demonstrated that histone variant H2A.X is highly enriched in human cleavage embryos when ZGA genes and TEs are active. Therefore, we propose that H2A.X plays an important role in regulating ZGA genes and TEs to establish totipotency.

Bookmarking by specific and nonspecific binding of FoxA1 pioneer factor to mitotic chromosomes.

While most transcription factors exit the chromatin during mitosis and the genome becomes silent, a subset of factors remains and "bookmarks" genes for rapid reactivation as cells progress through the cell cycle. However, it is unknown whether such bookmarking factors bind to chromatin similarly in mitosis and how different binding capacities among them relate to function. We compared a diverse set of transcription factors involved in liver differentiation and found markedly different extents of mitotic chromosome binding. Among them, the pioneer factor FoxA1 exhibits the greatest extent of mitotic chromosome binding. Genomically, ~15% of the FoxA1 interphase target sites are bound in mitosis, including at genes that are important for liver differentiation. Biophysical, genome mapping, and mutagenesis studies of FoxA1 reveals two different modes of binding to mitotic chromatin. Specific binding in mitosis occurs at sites that continue to be bound from interphase. Nonspecific binding in mitosis occurs across the chromosome due to the intrinsic chromatin affinity of FoxA1. Both specific and nonspecific binding contribute to timely reactivation of target genes post-mitosis. These studies reveal an unexpected diversity in the mechanisms by which transcription factors help retain cell identity during mitosis.

Stable intronic sequence RNAs (sisRNAs) are selected regions in introns with distinct properties.

BACKGROUND: Stable introns and intronic fragments make up the largest population of RNA in the oocyte nucleus of the frog Xenopus tropicalis. These stable intronic sequence RNAs (sisRNAs) persist through the onset of zygotic transcription when synchronous cell division has ended, and the developing embryo consists of approximately 8000 cells. Despite their abundance, the sequence properties and biological function of sisRNAs are just beginning to be understood. RESULTS: To characterize this population of non-coding RNA, we identified all of the sisRNAs in the X. tropicalis oocyte nucleus using published high-throughput RNA sequencing data. Our analysis revealed that sisRNAs, have an average length of ~ 360 nt, are widely expressed from genes with multiple introns, and are derived from specific regions of introns that are GC and TG rich, while CpG poor. They are enriched in introns at both ends of transcripts but preferentially at the 3' end. The consensus binding sites of specific transcription factors such as Stat3 are enriched in sisRNAs, suggesting an association between sisRNAs and transcription factors involved in early development. Evolutionary conservation analysis of sisRNA sequences in seven vertebrate genomes indicates that sisRNAs are as conserved as other parts of introns, but much less conserved than exons. CONCLUSION: In total, our results indicate sisRNAs are selected intron regions with distinct properties and may play a role in gene expression regulation.

The Epstein-Barr Virus B-ZIP Protein Zta Recognizes Specific DNA Sequences Containing 5-Methylcytosine and 5-Hydroxymethylcytosine.

The Epstein-Barr virus (EBV) B-ZIP transcription factor Zta binds to many DNA sequences containing methylated CG dinucleotides. Using protein binding microarrays (PBMs), we analyzed the sequence specific DNA binding of Zta to four kinds of double-stranded DNA (dsDNA): (1) DNA containing cytosine in both strands, (2) DNA with 5-methylcytosine (5mC) in one strand and cytosine in the second strand, (3) DNA with 5-hydroxymethylcytosine (5hmC) in one strand and cytosine in the second strand, and (4) DNA in which both cytosines in all CG dinucleotides contain 5mC. We compared these data to PBM data for three additional B-ZIP proteins (CREB1 and CEBPB homodimers and cJun|cFos heterodimers). With cytosine, Zta binds the TRE motif TGAC/GTCA as previously reported. With CG dinucleotides containing 5mC on both strands, many TRE motif variants containing a methylated CG dinucleotide at two positions in the motif, such as MGAGTCA and TGAGMGA (where M = 5mC), were preferentially bound. 5mC inhibits binding of Zta to both TRE motif half-sites GTCA and CTCA. Like the CREB1 homodimer, the Zta homodimer and the cJun|cFos heterodimer more strongly bind the C/EBP half-site tetranucleotide GCAA when it contains 5mC. Zta also binds dsDNA sequences containing 5hmC in one strand, although the effect is less dramatic than that observed for 5mC. Our results identify new DNA sequences that are well-bound by the viral B-ZIP protein Zta only when they contain 5mC or 5hmC, uncovering the potential for discovery of new viral and host regulatory programs controlled by EBV.

High levels of serum ß2-microglobulin predict severity of coronary artery disease.

BACKGROUND: The identification of new risk factors for coronary artery disease (CAD) is increasingly sought in an effort to tackle this threatening disease. ß2-microglobulin (B2M) is reported to associate with peripheral arterial disease and adverse cardiovascular outcomes. However, the association between B2M and cardiovascular disease remains under-researched. This study evaluated the effects of B2M on CAD without renal dysfunction. METHODS: One thousand seven hundred sixty-two subjects (403 non-CAD subjects and 1,359 CAD subjects) were investigated. Fasting samples were collected to determine B2M level. The Gensini and SYNTAX scores were used to assess the severity of CAD. RESULTS: CAD subjects were significantly higher in serum B2M level comparing with non-CAD subjects (1.25 ± 0.46 vs 1.14 ± 0.28 mg/L, p < 0.001). Serum B2M level was a risk factor of CAD after adjusting potential confounders (Odds Ratio (OR) = 2.363, 95% confidence interval (CI): 1.467-3.906, p = 0.001). Receiver operating characteristics (ROC) showed B2M level moderately predicted diagnosis of CAD (the area under the ROC curve (AUC) = 0.608, 95% CI: 0.577-0.639, p < 0.001). Furthermore, serum B2M level was positively associated with Gensini score system, SYNTAX score system and the number of disease vessels (NDV ≥ 2). CONCLUSIONS: The significant association between serum B2M and CAD suggests that B2M could be a biomarker for CAD.

5-Methylcytosine (5mC) and 5-Hydroxymethylcytosine (5hmC) Enhance the DNA Binding of CREB1 to the C/EBP Half-Site Tetranucleotide GCAA.

In human and mouse stem cells and brain, 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) can occur outside of CG dinucleotides. Using protein binding microarrays (PBMs) containing 60-mer DNA probes, we evaluated the effect of 5mC and 5hmC on one DNA strand on the double-stranded DNA binding of the mouse B-ZIP transcription factors (TFs) CREB1, ATF1, and JUND. 5mC inhibited binding of CREB1 to the canonical CRE half-site |GTCA but enhanced binding to the C/EBP half-site |GCAA. 5hmC inhibited binding of CREB1 to all 8-mers except TGAT|GCAA, where binding is enhanced. We observed similar DNA binding patterns with ATF1, a closely related B-ZIP domain. In contrast, both 5mC and 5hmC inhibited binding of JUND. These results identify new DNA sequences that are well-bound by CREB1 and ATF1 only when they contain 5mC or 5hmC. Analysis of two X-ray structures examines the consequences of 5mC and 5hmC on DNA binding by CREB and FOS|JUN.

CG methylated microarrays identify a novel methylated sequence bound by the CEBPB|ATF4 heterodimer that is active in vivo.

To evaluate the effect of CG methylation on DNA binding of sequence-specific B-ZIP transcription factors (TFs) in a high-throughput manner, we enzymatically methylated the cytosine in the CG dinucleotide on protein binding microarrays. Two Agilent DNA array designs were used. One contained 40,000 features using de Bruijn sequences where each 8-mer occurs 32 times in various positions in the DNA sequence. The second contained 180,000 features with each CG containing 8-mer occurring three times. The first design was better for identification of binding motifs, while the second was better for quantification. Using this novel technology, we show that CG methylation enhanced binding for CEBPA and CEBPB and inhibited binding for CREB, ATF4, JUN, JUND, CEBPD, and CEBPG. The CEBPB|ATF4 heterodimer bound a novel motif CGAT|GCAA 10-fold better when methylated. The electrophoretic mobility shift assay (EMSA) confirmed these results. CEBPB ChIP-seq data using primary female mouse dermal fibroblasts with 50× methylome coverage for each strand indicate that the methylated sequences well-bound on the arrays are also bound in vivo. CEBPB bound 39% of the methylated canonical 10-mers ATTGC|GCAAT in the mouse genome. After ATF4 protein induction by thapsigargin which results in ER stress, CEBPB binds methylated CGAT|GCAA in vivo, recapitulating what was observed on the arrays. This methodology can be used to identify new methylated DNA sequences preferentially bound by TFs, which may be functional in vivo.

Investigating the Potential Shared Molecular Mechanisms between COVID-19 and Alzheimer's Disease via Transcriptomic Analysis.

SARS-CoV-2 caused the COVID-19 pandemic. COVID-19 may elevate the risk of cognitive impairment and even cause dementia in infected individuals; it may accelerate cognitive decline in elderly patients with dementia, possibly in Alzheimer's disease (AD) patients. However, the mechanisms underlying the interplay between AD and COVID-19 are still unclear. To investigate the underlying mechanisms and associations between AD progression and SARS-CoV-2 infection, we conducted a series of bioinformatics research into SARS-CoV-2-infected cells, COVID-19 patients, AD patients, and SARS-CoV-2-infected AD patients. We identified the common differentially expressed genes (DEGs) in COVID-19 patients, AD patients, and SARS-CoV-2-infected cells, and these DEGs are enriched in certain pathways, such as immune responses and cytokine storms. We constructed the gene interaction network with the signaling transduction module in the center and identified IRF7, STAT1, STAT2, and OAS1 as the hub genes. We also checked the correlations between several key transcription factors and the SARS-CoV-2 and COVID-19 pathway-related genes. We observed that ACE2 expression is positively correlated with IRF7 expression in AD and coronavirus infections, and interestingly, IRF7 is significantly upregulated in response to different RNA virus infections. Further snRNA-seq analysis indicates that NRGN neurons or endothelial cells may be responsible for the increase in ACE2 and IRF7 expression after SARS-CoV-2 infection. The positive correlation between ACE2 and IRF7 expressions is confirmed in the hippocampal formation (HF) of SARS-CoV-2-infected AD patients. Our findings could contribute to the investigation of the molecular mechanisms underlying the interplay between AD and COVID-19 and to the development of effective therapeutic strategies for AD patients with COVID-19.

PARP1 interacts with WDR5 to enhance target gene recognition and facilitate tumorigenesis.

Poly (ADP-ribose) polymerase-1 (PARP1) is a nuclear protein that attaches negatively charged poly (ADP-ribose) (PAR) to itself and other target proteins. While its function in DNA damage repair is well established, its role in target chromatin recognition and regulation of gene expression remains to be better understood. This study showed that PARP1 interacts with SET1/MLL complexes by binding directly to WDR5. Notably, although PARP1 does not modulate WDR5 PARylation or the global level of H3K4 methylation, it exerts locus-specific effects on WDR5 binding and H3K4 methylation. Interestingly, PARP1 and WDR5 show extensive co-localization on chromatin, with WDR5 facilitating the recognition and expression of target genes regulated by PARP1. Furthermore, we demonstrated that inhibition of the WDR5 Win site impedes the interaction between PARP1 and WDR5, thereby inhibiting PARP1 from binding to target genes. Finally, the combined inhibition of the WDR5 Win site and PARP shows a profound inhibitory effect on the proliferation of cancer cells. These findings illuminate intricate mechanisms underlying chromatin recognition, gene transcription, and tumorigenesis, shedding light on previously unrecognized roles of PARP1 and WDR5 in these processes.

Single cell transcriptome profiling of infrapatellar fat pad highlights the role of interstitial inflammatory fibroblasts in osteoarthritis.

OBJECTIVES: Osteoarthritis (OA) is a whole-joint disease in which the role of the infrapatellar fat pad (IFP) in its pathogenesis is unclear. Our study explored the cellular heterogeneity of IFP to understand OA and identify therapeutic targets. METHODS: Single-cell and single-nuclei RNA sequencing were used to analyze 10 IFP samples, comprising 5 from OA patients and 5 from healthy controls. Analyses included differential gene expression, enrichment, pseudotime trajectory, and cellular communication, along with comparative studies with visceral and subcutaneous fats. Key subcluster and pathways were validated using multiplex immunohistochemistry. RESULTS: The scRNA-seq performed on the IFPs of the OA and control group profiled the gene expressions of over 49,674 cells belonging to 11 major cell types. We discovered that adipose stem and progenitor cells (ASPCs), contributing to the formation of both adipocytes and synovial-lining fibroblasts (SLF). Interstitial inflammatory fibroblasts (iiFBs) were a subcluster of ASPCs that exhibit notable pro-inflammatory and proliferative characteristics. We identified four adipocyte subtypes, with one subtype showing a reduced lipid synthesis ability. Furthermore, iiFBs modulated the activities of macrophages and T cells in the IFP. Compared to subcutaneous and visceral adipose tissues, iiFBs represented a distinctive subpopulation of ASPCs in IFP that regulated cartilage proliferation through the MK pathway. CONCLUSION: This study presents a comprehensive single-cell transcriptomic atlas of IFP, uncovering its complex cellular landscape and potential impact on OA progression. Our findings highlight the role of iiFBs in OA, especially through MK pathway, opening new avenues for understanding OA pathogenesis and developing novel targeted therapies.

Chromatin Modifier EP400 Regulates Oocyte Quality and Zygotic Genome Activation in Mice.

Epigenetic modifiers that accumulate in oocytes, play a crucial role in steering the developmental program of cleavage embryos and initiating life. However, the identification of key maternal epigenetic regulators remains elusive. In the findings, the essential role of maternal Ep400, a chaperone for H3.3, in oocyte quality and early embryo development in mice is highlighted. Depletion of Ep400 in oocytes resulted in a decline in oocyte quality and abnormalities in fertilization. Preimplantation embryos lacking maternal Ep400 exhibited reduced major zygotic genome activation (ZGA) and experienced developmental arrest at the 2-to-4-cell stage. The study shows that EP400 forms protein complex with NFYA, occupies promoters of major ZGA genes, modulates H3.3 distribution between euchromatin and heterochromatin, promotes transcription elongation, activates the expression of genes regulating mitochondrial functions, and facilitates the expression of rate-limiting enzymes of the TCA cycle. This intricate process driven by Ep400 ensures the proper execution of the developmental program, emphasizing its critical role in maternal-to-embryonic transition.

Contribution of nucleosome binding preferences and co-occurring DNA sequences to transcription factor binding.

BACKGROUND: Chromatin plays a critical role in regulating transcription factors (TFs) binding to their canonical transcription factor binding sites (TFBS). Recent studies in vertebrates show that many TFs preferentially bind to genomic regions that are well bound by nucleosomes in vitro. Co-occurring secondary motifs sometimes correlated with functional TFBS. RESULTS: We used a logistic regression to evaluate how well the propensity for nucleosome binding and co-occurrence of a secondary motif identify which canonical motifs are bound in vivo. We used ChIP-seq data for three transcription factors binding to their canonical motifs: c-Jun binding the AP-1 motif (TGA(C)/(G)TCA), GR (glucocorticoid receptor) binding the GR motif (G-ACA---(T)/(C)GT-C), and Hoxa2 (homeobox a2) binding the Pbx (Pre-B-cell leukemia homeobox) motif (TGATTGAT). For all canonical TFBS in the mouse genome, we calculated intrinsic nucleosome occupancy scores (INOS) for its surrounding 150-bps DNA and examined the relationship with in vivo TF binding. In mouse mammary 3134 cells, c-Jun and GR proteins preferentially bound regions calculated to be well-bound by nucleosomes in vitro with the canonical AP-1 and GR motifs themselves contributing to the high INOS. Functional GR motifs are enriched for AP-1 motifs if they are within a nucleosome-sized 150-bps region. GR and Hoxa2 also bind motifs with low INOS, perhaps indicating a different mechanism of action. CONCLUSION: Our analysis quantified the contribution of INOS and co-occurring sequence to the identification of functional canonical motifs in the genome. This analysis revealed an inherent competition between some TFs and nucleosomes for binding canonical TFBS. GR and c-Jun cooperate if they are within 150-bps. Binding of Hoxa2 and a fraction of GR to motifs with low INOS values suggesting they are not in competition with nucleosomes and may function using different mechanisms.

Evolution of an antifreeze protein by neofunctionalization under escape from adaptive conflict.

The evolutionary model escape from adaptive conflict (EAC) posits that adaptive conflict between the old and an emerging new function within a single gene could drive the fixation of gene duplication, where each duplicate can freely optimize one of the functions. Although EAC has been suggested as a common process in functional evolution, definitive cases of neofunctionalization under EAC are lacking, and the molecular mechanisms leading to functional innovation are not well-understood. We report here clear experimental evidence for EAC-driven evolution of type III antifreeze protein gene from an old sialic acid synthase (SAS) gene in an Antarctic zoarcid fish. We found that an SAS gene, having both sialic acid synthase and rudimentary ice-binding activities, became duplicated. In one duplicate, the N-terminal SAS domain was deleted and replaced with a nascent signal peptide, removing pleiotropic structural conflict between SAS and ice-binding functions and allowing rapid optimization of the C-terminal domain to become a secreted protein capable of noncolligative freezing-point depression. This study reveals how minor functionalities in an old gene can be transformed into a distinct survival protein and provides insights into how gene duplicates facing presumed identical selection and mutation pressures at birth could take divergent evolutionary paths.

Corrigendum: A genome-wide investigation of effects of aberrant DNA methylation on the usage of alternative promoters in hepatocellular carcinoma.

[This corrects the article DOI: 10.3389/fonc.2021.780266.].

Construction and Verification of an RNA-Binding Protein-Associated Prognostic Model for Gliomas.

OBJECTIVE: To construct and verificate an RNA-binding protein (RBP)-associated prognostic model for gliomas using integrated bioinformatics analysis. METHODS: RNA-sequencing and clinic pathological data of glioma patients from The Cancer Genome Atlas (TCGA) database and the Chinese Glioma Genome Atlas database (CGGA) were downloaded. The aberrantly expressed RBPs were investigated between gliomas and normal samples in TCGA database. We then identified prognosis related hub genes and constructed a prognostic model. This model was further validated in the CGGA-693 and CGGA-325 cohorts. RESULTS: Totally 174 differently expressed genes-encoded RBPs were identified, containing 85 down-regulated and 89 up-regulated genes. We identified five genes-encoded RBPs (ERI1, RPS2, BRCA1, NXT1, and TRIM21) as prognosis related key genes and constructed a prognostic model. Overall survival (OS) analysis revealed that the patients in the high-risk subgroup based on the model were worse than those in the low-risk subgroup. The area under the receiver operator characteristic curve (AUC) of the prognostic model was 0.836 in the TCGA dataset and 0.708 in the CGGA-693 dataset, demonstrating a favorable prognostic model. Survival analyses of the five RBPs in the CGGA-325 cohort validated the findings. A nomogram was constructed based on the five genes and validated in the TCGA cohort, confirming a promising discriminating ability for gliomas. CONCLUSION: The prognostic model of the five RBPs might serve as an independent prognostic algorithm for gliomas.