Search

Coronaviruses constitute a global threat to human and animal health. It is essential to investigate the long-distance RNA-RNA interactions that approximate remote regulatory elements in strategies, including genome circularization, discontinuous transcription, and transcriptional enhancers, aimed at the rapid replication of their large genomes, pathogenicity, and immune evasion. Based on the primary sequences and modeled RNA-RNA interactions of two experimentally defined coronaviral enhancers, we detected via an in silico primary and secondary structural analysis potential enhancers in various coronaviruses, from the phylogenetically ancient avian infectious bronchitis virus (IBV) to the recently emerged SARS-CoV-2. These potential enhancers possess a core duplex-forming region that could transition between closed and open states, as molecular switches directed by viral or host factors. The duplex open state would pair with remote sequences in the viral genome and modulate the expression of downstream crucial genes involved in viral replication and host immune evasion. Consistently, variations in the predicted IBV enhancer region or its distant targets coincide with cases of viral attenuation, possibly driven by decreased open reading frame (ORF)3a immune evasion protein expression. If validated experimentally, the annotated enhancer sequences could inform structural prediction tools and antiviral interventions. [ABSTRACT FROM AUTHOR]

High-resolution melting (HRM) is a cost-efficient tool for targeted DNA methylation analysis. HRM yields the average methylation status across all CpGs in PCR products. Moreover, it provides information on the methylation pattern, e.g., the occurrence of monoallelic methylation. HRM assays have to be calibrated by analyzing DNA methylation standards of known methylation status and mixtures thereof. In general, DNA methylation levels determined by the classical calibration approach, including the whole temperature range in between normalization intervals, are in good agreement with the mean of the DNA methylation status of individual CpGs determined by pyrosequencing (PSQ), the gold standard of targeted DNA methylation analysis. However, the classical calibration approach leads to highly inaccurate results for samples with heterogeneous DNA methylation since they result in more complex melt curves, differing in their shape compared to those of DNA standards and mixtures thereof. Here, we present a novel calibration approach, i.e., temperature-wise calibration. By temperature-wise calibration, methylation profiles over temperature are obtained, which help in finding the optimal calibration range and thus increase the accuracy of HRM data, particularly for heterogeneous DNA methylation. For explaining the principle and demonstrating the potential of the novel calibration approach, we selected the promoter and two enhancers of MGMT, a gene encoding the repair protein MGMT. [ABSTRACT FROM AUTHOR]

Coronaviruses constitute a global threat to human and animal health. It is essential to investigate the long-distance RNA-RNA interactions that approximate remote regulatory elements in strategies, including genome circularization, discontinuous transcription, and transcriptional enhancers, aimed at the rapid replication of their large genomes, pathogenicity, and immune evasion. Based on the primary sequences and modeled RNA-RNA interactions of two experimentally defined coronaviral enhancers, we detected via an in silico primary and secondary structural analysis potential enhancers in various coronaviruses, from the phylogenetically ancient avian infectious bronchitis virus (IBV) to the recently emerged SARS-CoV-2. These potential enhancers possess a core duplex-forming region that could transition between closed and open states, as molecular switches directed by viral or host factors. The duplex open state would pair with remote sequences in the viral genome and modulate the expression of downstream crucial genes involved in viral replication and host immune evasion. Consistently, variations in the predicted IBV enhancer region or its distant targets coincide with cases of viral attenuation, possibly driven by decreased open reading frame (ORF)3a immune evasion protein expression. If validated experimentally, the annotated enhancer sequences could inform structural prediction tools and antiviral interventions.

Enhancers are located upstream or downstream of key deoxyribonucleic acid (DNA) sequences in genes and can adjust the transcription activity of neighboring genes. Identifying enhancers and determining their functions are important for understanding gene regulatory networks and expression regulatory mechanisms. However, traditional enhancer recognition relies on manual feature engineering, which is time-consuming and labor-intensive, making it difficult to perform large-scale recognition analysis. In addition, if the original dataset is too small, there is a risk of overfitting. In recent years, emerging methods, such as deep learning, have provided new insights for enhancing identification. However, these methods also present certain challenges. Deep learning models typically require a large amount of high-quality data, and data acquisition demands considerable time and resources. To address these challenges, in this paper, we propose a data-augmentation method based on generative adversarial networks to solve the problem of small datasets. Moreover, we used regularization methods such as weight decay to improve the generalizability of the model and alleviate overfitting. The Transformer encoder was used as the main component to capture the complex relationships and dependencies in enhancer sequences. The encoding layer was designed based on the principle of k-mers to preserve more information from the original DNA sequence. Compared with existing methods, the proposed approach made significant progress in enhancing the accuracy and strength of enhancer identification and prediction, demonstrating the effectiveness of the proposed method. This paper provides valuable insights for enhancer analysis and is of great significance for understanding gene regulatory mechanisms and studying disease correlations. [ABSTRACT FROM AUTHOR]

Ferroptosis is an iron-dependent form of cell death, which is reported to be associated with glioma progression and drug sensitivity. Targeting ferroptosis is a potential therapeutic approach for glioma. However, the molecular mechanism of glioma cell ferroptosis is not clear. In this study, we profile the change of 3D chromatin structure in glioblastoma ferroptosis by using HiChIP and study the 3D gene regulation network in glioblastoma ferroptosis. A combination of an analysis of HiChIP and RNA-seq data suggests that change of chromatin loops mediated by 3D chromatin structure regulates gene expressions in glioblastoma ferroptosis. Genes that are regulated by 3D chromatin structures include genes that were reported to function in ferroptosis, like HDM2 and TXNRD1. We propose a new regulatory mechanism governing glioblastoma cell ferroptosis by 3D chromatin structure. [ABSTRACT FROM AUTHOR]

The establishment and stable inheritance of individual patterns of gene expression in different cell types are required for the development of multicellular organisms. The important epigenetic regulators are the Polycomb group (PcG) and Trithorax group (TrxG) proteins, which control the silenced and active states of genes, respectively. In Drosophila, the PcG/TrxG group proteins are recruited to the DNA regulatory sequences termed the Polycomb response elements (PREs). The PREs are composed of the binding sites for different DNA-binding proteins, the so-called PcG recruiters. Currently, the role of the PcG recruiters in the targeting of the PcG proteins to PREs is well documented. However, there are examples where the PcG recruiters are also implicated in the active transcription and in the TrxG function. In addition, there is increasing evidence that the genome-wide PcG recruiters interact with the chromatin outside of the PREs and overlap with the proteins of differing regulatory classes. Recent studies of the interactomes of the PcG recruiters significantly expanded our understanding that they have numerous interactors besides the PcG proteins and that their functions extend beyond the regulation of the PRE repressive activity. Here, we summarize current data about the functions of the PcG recruiters. [ABSTRACT FROM AUTHOR]

Background: High expression of inhibitor of DNA binding 1 (ID1) correlates with poor prognosis in colorectal cancer (CRC). Aberrant enhancer activation in regulating ID1 transcription is limited. Methods: Immunohistochemistry (IHC), quantitative RT-PCR (RT-qPCR) and Western blotting (WB) were used to determine the expression of ID1. CRISPR-Cas9 was used to generate ID1 or enhancer E1 knockout cell lines. Dual-luciferase reporter assay, chromosome conformation capture assay and ChIP-qPCR were used to determine the active enhancers of ID1. Cell Counting Kit 8, colony-forming, transwell assays and tumorigenicity in nude mice were used to investigate the biological functions of ID1 and enhancer E1. Results: Human CRC tissues and cell lines expressed a higher level of ID1 than normal controls. ID1 promoted CRC cell proliferation and colony formation. Enhancer E1 actively regulated ID1 promoter activity. Signal transducer and activator of transcription 3 (STAT3) bound to ID1 promoter and enhancer E1 to regulate their activity. The inhibitor of STAT3 Stattic attenuated ID1 promoter and enhancer E1 activity and the expression of ID1. Enhancer E1 knockout down-regulated ID1 expression level and cell proliferation in vitro and in vivo. Conclusions: Enhancer E1 is positively regulated by STAT3 and contributes to the regulation of ID1 to promote CRC cell progression and might be a potential target for anti-CRC drug studies. [ABSTRACT FROM AUTHOR]

Keratins are a family of intermediate filament-forming proteins highly specific to epithelial cells. A combination of expressed keratin genes is a defining property of the epithelium belonging to a certain type, organ/tissue, cell differentiation potential, and at normal or pathological conditions. In a variety of processes such as differentiation and maturation, as well as during acute or chronic injury and malignant transformation, keratin expression undergoes switching: an initial keratin profile changes accordingly to changed cell functions and location within a tissue as well as other parameters of cellular phenotype and physiology. Tight control of keratin expression implies the presence of complex regulatory landscapes within the keratin gene loci. Here, we highlight patterns of keratin expression in different biological conditions and summarize disparate data on mechanisms controlling keratin expression at the level of genomic regulatory elements, transcription factors (TFs), and chromatin spatial structure. [ABSTRACT FROM AUTHOR]

Nuclear receptor subfamily 5 group A member 1 (NR5A1) is expressed in the pituitary gonadotrope and regulates their differentiation. Although several regulatory regions were implicated in Nr5a1 gene expression in the pituitary gland, none of these regions have been verified using mouse models. Furthermore, the molecular functions of NR5A1 in the pituitary gonadotrope have not been fully elucidated. In the present study, we generated mice lacking the pituitary enhancer located in the 6th intron of the Nr5a1 gene. These mice showed pituitary gland-specific disappearance of NR5A1, confirming the functional importance of the enhancer. Enhancer-deleted male mice demonstrated no defects at fetal stages. Meanwhile, androgen production decreased markedly in adult, and postnatal development of reproductive organs, such as the seminal vesicle, prostate, and penis was severely impaired. We further performed transcriptomic analyses of the whole pituitary gland of the enhancer-deleted mice and controls, as well as gonadotropes isolated from Ad4BP-BAC-EGFP mice. These analyses identified several genes showing gonadotrope-specific, NR5A1-dependent expressions, such as Spp1, Tgfbr3l, Grem1, and Nr0b2. These factors are thought to function downstream of NR5A1 and play important roles in reproductive organ development through regulation of pituitary gonadotrope functions. [ABSTRACT FROM AUTHOR]

Enhancers are located upstream or downstream of key deoxyribonucleic acid (DNA) sequences in genes and can adjust the transcription activity of neighboring genes. Identifying enhancers and determining their functions are important for understanding gene regulatory networks and expression regulatory mechanisms. However, traditional enhancer recognition relies on manual feature engineering, which is time-consuming and labor-intensive, making it difficult to perform large-scale recognition analysis. In addition, if the original dataset is too small, there is a risk of overfitting. In recent years, emerging methods, such as deep learning, have provided new insights for enhancing identification. However, these methods also present certain challenges. Deep learning models typically require a large amount of high-quality data, and data acquisition demands considerable time and resources. To address these challenges, in this paper, we propose a data-augmentation method based on generative adversarial networks to solve the problem of small datasets. Moreover, we used regularization methods such as weight decay to improve the generalizability of the model and alleviate overfitting. The Transformer encoder was used as the main component to capture the complex relationships and dependencies in enhancer sequences. The encoding layer was designed based on the principle of k-mers to preserve more information from the original DNA sequence. Compared with existing methods, the proposed approach made significant progress in enhancing the accuracy and strength of enhancer identification and prediction, demonstrating the effectiveness of the proposed method. This paper provides valuable insights for enhancer analysis and is of great significance for understanding gene regulatory mechanisms and studying disease correlations.

Ferroptosis is an iron-dependent form of cell death, which is reported to be associated with glioma progression and drug sensitivity. Targeting ferroptosis is a potential therapeutic approach for glioma. However, the molecular mechanism of glioma cell ferroptosis is not clear. In this study, we profile the change of 3D chromatin structure in glioblastoma ferroptosis by using HiChIP and study the 3D gene regulation network in glioblastoma ferroptosis. A combination of an analysis of HiChIP and RNA-seq data suggests that change of chromatin loops mediated by 3D chromatin structure regulates gene expressions in glioblastoma ferroptosis. Genes that are regulated by 3D chromatin structures include genes that were reported to function in ferroptosis, like HDM2 and TXNRD1. We propose a new regulatory mechanism governing glioblastoma cell ferroptosis by 3D chromatin structure.

A library of active genome regulatory elements (putative promoters and enhancers) from MIA PaCa-2 pancreatic adenocarcinoma cells was constructed using a specially designed lentiviral vector and a massive parallel reporter assay (ChIP-lentiMPRA). Chromatin immunoprecipitation of the cell genomic DNA by H3K27ac antibodies was used for primary enrichment of the library for regulatory elements. Totally, 11,264 unique genome regions, many of which are capable of enhancing the expression of the CopGFP reporter gene from the minimal CMV promoter, were identified. The regions tend to be located near promoters. Based on the proximity assay, we found an enrichment of highly expressed genes among those associated with three or more mapped distal regions (2 kb distant from the 5′-ends of genes). It was shown significant enrichment of genes related to carcinogenesis or Mia PaCa-2 cell identity genes in this group. In contrast, genes associated with 1–2 distal regions or only with proximal regions (within 2 kbp of the 5′-ends of genes) are more often related to housekeeping functions. Thus, ChIP-lentiMPRA is a useful strategy for creating libraries of regulatory elements for the study of tumor-specific gene transcription. [ABSTRACT FROM AUTHOR]

Enhancers in higher eukaryotes and upstream activating sequences (UASs) in yeast have been shown to recruit components of the RNA polymerase II (Pol II) transcription machinery. At least a fraction of Pol II recruited to enhancers in higher eukaryotes initiates transcription and generates enhancer RNA (eRNA). In contrast, UASs in yeast do not recruit transcription factor TFIIH, which is required for transcription initiation. For both yeast and mammalian systems, it was shown that Pol II is transferred from enhancers/UASs to promoters. We propose that there are two modes of Pol II recruitment to enhancers in higher eukaryotes. Pol II complexes that generate eRNAs are recruited via TFIID, similar to mechanisms operating at promoters. This may involve the binding of TFIID to acetylated nucleosomes flanking the enhancer. The resulting eRNA, together with enhancer-bound transcription factors and co-regulators, contributes to the second mode of Pol II recruitment through the formation of a transcription initiation domain. Transient contacts with target genes, governed by proteins and RNA, lead to the transfer of Pol II from enhancers to TFIID-bound promoters. [ABSTRACT FROM AUTHOR]

The dynamic balance of transcriptional and translational regulation together with degron-controlled proteolysis shapes the ever-changing cellular proteome. While a large variety of degradation signals has been characterized, our knowledge of cis-acting protein motifs that can in vivo stabilize otherwise short-lived proteins is very limited. We have identified and characterized a conserved 13-mer protein segment derived from the p54/Rpn10 ubiquitin receptor subunit of the Drosophila 26S proteasome, which fulfills all the characteristics of a protein stabilization motif (STABILON). Attachment of STABILON to various intracellular as well as medically relevant secreted model proteins resulted in a significant increase in their cellular or extracellular concentration in mammalian cells. We demonstrate that STABILON acts as a universal and dual function motif that, on the one hand, increases the concentration of the corresponding mRNAs and, on the other hand, prevents the degradation of short-lived fusion proteins. Therefore, STABILON may lead to a breakthrough in biomedical recombinant protein production. [ABSTRACT FROM AUTHOR]

Background: High expression of inhibitor of DNA binding 1 (ID1) correlates with poor prognosis in colorectal cancer (CRC). Aberrant enhancer activation in regulating ID1 transcription is limited. Methods: Immunohistochemistry (IHC), quantitative RT-PCR (RT-qPCR) and Western blotting (WB) were used to determine the expression of ID1. CRISPR-Cas9 was used to generate ID1 or enhancer E1 knockout cell lines. Dual-luciferase reporter assay, chromosome conformation capture assay and ChIP-qPCR were used to determine the active enhancers of ID1. Cell Counting Kit 8, colony-forming, transwell assays and tumorigenicity in nude mice were used to investigate the biological functions of ID1 and enhancer E1. Results: Human CRC tissues and cell lines expressed a higher level of ID1 than normal controls. ID1 promoted CRC cell proliferation and colony formation. Enhancer E1 actively regulated ID1 promoter activity. Signal transducer and activator of transcription 3 (STAT3) bound to ID1 promoter and enhancer E1 to regulate their activity. The inhibitor of STAT3 Stattic attenuated ID1 promoter and enhancer E1 activity and the expression of ID1. Enhancer E1 knockout down-regulated ID1 expression level and cell proliferation in vitro and in vivo. Conclusions: Enhancer E1 is positively regulated by STAT3 and contributes to the regulation of ID1 to promote CRC cell progression and might be a potential target for anti-CRC drug studies.

Myogenesis is a central step in prenatal myofiber formation, postnatal myofiber hypertrophy, and muscle damage repair in adulthood. RNA-Seq technology has greatly helped reveal the molecular mechanism of myogenesis, but batch effects in different experiments inevitably lead to misinterpretation of differentially expressed genes (DEGs). We previously applied the robust rank aggregation (RRA) method to effectively circumvent batch effects across multiple RNA-Seq datasets from 3T3-L1 cells. Here, we also used the RRA method to integrate nine RNA-Seq datasets from C2C12 cells and obtained 3140 robust DEGs between myoblasts and myotubes, which were then validated with array expression profiles and H3K27ac signals. The upregulated robust DEGs were highly enriched in gene ontology (GO) terms related to muscle cell differentiation and development. Considering that the cooperative binding of transcription factors (TFs) to enhancers to regulate downstream gene expression is a classical epigenetic mechanism, differentially expressed TFs (DETFs) were screened, and potential novel myogenic factors (MAF, BCL6, and ESR1) with high connection degree in protein–protein interaction (PPI) network were presented. Moreover, KLF5 cooperatively binds with the three key myogenic factors (MYOD, MYOG, and MEF2D) in C2C12 cells. Motif analysis speculates that the binding of MYOD and MYOG is KLF5-independent, while MEF2D is KLF5-dependent. It was revealed that KLF5-binding sites could be exploited to filter redundant MYOD-, MYOG-, and MEF2D-binding sites to focus on key enhancers for myogenesis. Further functional annotation of KLF5-binding sites suggested that KLF5 may regulate myogenesis through the PI3K-AKt signaling pathway, Rap1 signaling pathway, and the Hippo signaling pathway. In general, our study provides a wealth of untapped candidate targets for myogenesis and contributes new insights into the core regulatory mechanisms of myogenesis relying on KLF5-binding signal. [ABSTRACT FROM AUTHOR]

Genome-wide association studies for severe malaria (SM) have identified 30 genetic variants mostly located in non-coding regions. Here, we aimed to identify potential causal genetic variants located in these loci and demonstrate their functional activity. We systematically investigated the regulatory effect of the SNPs in linkage disequilibrium (LD) with the malaria-associated genetic variants. Annotating and prioritizing genetic variants led to the identification of a regulatory region containing five ATP2B4 SNPs in LD with rs10900585. We found significant associations between SM and rs10900585 and our candidate SNPs (rs11240734, rs1541252, rs1541253, rs1541254, and rs1541255) in a Senegalese population. Then, we demonstrated that both individual SNPs and the combination of SNPs had regulatory effects. Moreover, CRISPR/Cas9-mediated deletion of this region decreased ATP2B4 transcript and protein levels and increased Ca2+ intracellular concentration in the K562 cell line. Our data demonstrate that severe malaria-associated genetic variants alter the expression of ATP2B4 encoding a plasma membrane calcium-transporting ATPase 4 (PMCA4) expressed on red blood cells. Altering the activity of this regulatory element affects the risk of SM, likely through calcium concentration effect on parasitaemia. [ABSTRACT FROM AUTHOR]

Keratins are a family of intermediate filament-forming proteins highly specific to epithelial cells. A combination of expressed keratin genes is a defining property of the epithelium belonging to a certain type, organ/tissue, cell differentiation potential, and at normal or pathological conditions. In a variety of processes such as differentiation and maturation, as well as during acute or chronic injury and malignant transformation, keratin expression undergoes switching: an initial keratin profile changes accordingly to changed cell functions and location within a tissue as well as other parameters of cellular phenotype and physiology. Tight control of keratin expression implies the presence of complex regulatory landscapes within the keratin gene loci. Here, we highlight patterns of keratin expression in different biological conditions and summarize disparate data on mechanisms controlling keratin expression at the level of genomic regulatory elements, transcription factors (TFs), and chromatin spatial structure.

Nuclear receptor subfamily 5 group A member 1 (NR5A1) is expressed in the pituitary gonadotrope and regulates their differentiation. Although several regulatory regions were implicated in Nr5a1 gene expression in the pituitary gland, none of these regions have been verified using mouse models. Furthermore, the molecular functions of NR5A1 in the pituitary gonadotrope have not been fully elucidated. In the present study, we generated mice lacking the pituitary enhancer located in the 6th intron of the Nr5a1 gene. These mice showed pituitary gland-specific disappearance of NR5A1, confirming the functional importance of the enhancer. Enhancer-deleted male mice demonstrated no defects at fetal stages. Meanwhile, androgen production decreased markedly in adult, and postnatal development of reproductive organs, such as the seminal vesicle, prostate, and penis was severely impaired. We further performed transcriptomic analyses of the whole pituitary gland of the enhancer-deleted mice and controls, as well as gonadotropes isolated from Ad4BP-BAC-EGFP mice. These analyses identified several genes showing gonadotrope-specific, NR5A1-dependent expressions, such as Spp1, Tgfbr3l, Grem1, and Nr0b2. These factors are thought to function downstream of NR5A1 and play important roles in reproductive organ development through regulation of pituitary gonadotrope functions.

A library of active genome regulatory elements (putative promoters and enhancers) from MIA PaCa-2 pancreatic adenocarcinoma cells was constructed using a specially designed lentiviral vector and a massive parallel reporter assay (ChIP-lentiMPRA). Chromatin immunoprecipitation of the cell genomic DNA by H3K27ac antibodies was used for primary enrichment of the library for regulatory elements. Totally, 11,264 unique genome regions, many of which are capable of enhancing the expression of the CopGFP reporter gene from the minimal CMV promoter, were identified. The regions tend to be located near promoters. Based on the proximity assay, we found an enrichment of highly expressed genes among those associated with three or more mapped distal regions (2 kb distant from the 5′-ends of genes). It was shown significant enrichment of genes related to carcinogenesis or Mia PaCa-2 cell identity genes in this group. In contrast, genes associated with 1–2 distal regions or only with proximal regions (within 2 kbp of the 5′-ends of genes) are more often related to housekeeping functions. Thus, ChIP-lentiMPRA is a useful strategy for creating libraries of regulatory elements for the study of tumor-specific gene transcription.

Enhancers in higher eukaryotes and upstream activating sequences (UASs) in yeast have been shown to recruit components of the RNA polymerase II (Pol II) transcription machinery. At least a fraction of Pol II recruited to enhancers in higher eukaryotes initiates transcription and generates enhancer RNA (eRNA). In contrast, UASs in yeast do not recruit transcription factor TFIIH, which is required for transcription initiation. For both yeast and mammalian systems, it was shown that Pol II is transferred from enhancers/UASs to promoters. We propose that there are two modes of Pol II recruitment to enhancers in higher eukaryotes. Pol II complexes that generate eRNAs are recruited via TFIID, similar to mechanisms operating at promoters. This may involve the binding of TFIID to acetylated nucleosomes flanking the enhancer. The resulting eRNA, together with enhancer-bound transcription factors and co-regulators, contributes to the second mode of Pol II recruitment through the formation of a transcription initiation domain. Transient contacts with target genes, governed by proteins and RNA, lead to the transfer of Pol II from enhancers to TFIID-bound promoters.

The dynamic balance of transcriptional and translational regulation together with degron-controlled proteolysis shapes the ever-changing cellular proteome. While a large variety of degradation signals has been characterized, our knowledge of cis-acting protein motifs that can in vivo stabilize otherwise short-lived proteins is very limited. We have identified and characterized a conserved 13-mer protein segment derived from the p54/Rpn10 ubiquitin receptor subunit of the Drosophila 26S proteasome, which fulfills all the characteristics of a protein stabilization motif (STABILON). Attachment of STABILON to various intracellular as well as medically relevant secreted model proteins resulted in a significant increase in their cellular or extracellular concentration in mammalian cells. We demonstrate that STABILON acts as a universal and dual function motif that, on the one hand, increases the concentration of the corresponding mRNAs and, on the other hand, prevents the degradation of short-lived fusion proteins. Therefore, STABILON may lead to a breakthrough in biomedical recombinant protein production.

Myogenesis is a central step in prenatal myofiber formation, postnatal myofiber hypertrophy, and muscle damage repair in adulthood. RNA-Seq technology has greatly helped reveal the molecular mechanism of myogenesis, but batch effects in different experiments inevitably lead to misinterpretation of differentially expressed genes (DEGs). We previously applied the robust rank aggregation (RRA) method to effectively circumvent batch effects across multiple RNA-Seq datasets from 3T3-L1 cells. Here, we also used the RRA method to integrate nine RNA-Seq datasets from C2C12 cells and obtained 3140 robust DEGs between myoblasts and myotubes, which were then validated with array expression profiles and H3K27ac signals. The upregulated robust DEGs were highly enriched in gene ontology (GO) terms related to muscle cell differentiation and development. Considering that the cooperative binding of transcription factors (TFs) to enhancers to regulate downstream gene expression is a classical epigenetic mechanism, differentially expressed TFs (DETFs) were screened, and potential novel myogenic factors (MAF, BCL6, and ESR1) with high connection degree in protein–protein interaction (PPI) network were presented. Moreover, KLF5 cooperatively binds with the three key myogenic factors (MYOD, MYOG, and MEF2D) in C2C12 cells. Motif analysis speculates that the binding of MYOD and MYOG is KLF5-independent, while MEF2D is KLF5-dependent. It was revealed that KLF5-binding sites could be exploited to filter redundant MYOD-, MYOG-, and MEF2D-binding sites to focus on key enhancers for myogenesis. Further functional annotation of KLF5-binding sites suggested that KLF5 may regulate myogenesis through the PI3K-AKt signaling pathway, Rap1 signaling pathway, and the Hippo signaling pathway. In general, our study provides a wealth of untapped candidate targets for myogenesis and contributes new insights into the core regulatory mechanisms of myogenesis relying on KLF5-binding signal.

Genome-wide association studies for severe malaria (SM) have identified 30 genetic variants mostly located in non-coding regions. Here, we aimed to identify potential causal genetic variants located in these loci and demonstrate their functional activity. We systematically investigated the regulatory effect of the SNPs in linkage disequilibrium (LD) with the malaria-associated genetic variants. Annotating and prioritizing genetic variants led to the identification of a regulatory region containing five ATP2B4 SNPs in LD with rs10900585. We found significant associations between SM and rs10900585 and our candidate SNPs (rs11240734, rs1541252, rs1541253, rs1541254, and rs1541255) in a Senegalese population. Then, we demonstrated that both individual SNPs and the combination of SNPs had regulatory effects. Moreover, CRISPR/Cas9-mediated deletion of this region decreased ATP2B4 transcript and protein levels and increased Ca2+ intracellular concentration in the K562 cell line. Our data demonstrate that severe malaria-associated genetic variants alter the expression of ATP2B4 encoding a plasma membrane calcium-transporting ATPase 4 (PMCA4) expressed on red blood cells. Altering the activity of this regulatory element affects the risk of SM, likely through calcium concentration effect on parasitaemia.