Lung adenocarcinoma is one of the deadliest tumors. Studies have shown that N6-methyladenosine RNA methylation regulators, as a dynamic chemical modification, affect the occurrence and development of lung adenocarcinoma. To investigate the relationship between mutations and expression levels of m6A regulators in lung adenocarcinoma, we investigated the mutations and expression levels of 38 m6A regulators. We found that mutations in m6A regulatory factors did not affect the changes in expression levels, and 19 differentially expressed genes were identified. All tumor samples were classified into two subtypes based on the expression levels of 19 differentially expressed m6A-regulated genes. Survival analysis showed significant differences in survival between the two subtypes. To explore the relationship between immune cell infiltration and survival in both subtypes, we calculated the infiltration of 23 immune cells in both subtypes, and we found that the subtype with high immune cell infiltration had better survival. We found that subtypes with low tumor purity and high stromal and immune scores had better survival. The m6A-related immune genes were identified by taking the intersection of differentially expressed genes and immune genes in the two isoforms and calculating the Pearson correlation coefficients between the intersecting immune genes and the differentially expressed m6A-regulated genes. Finally, a prognostic model associated with m6A and associated with immunity was developed using prognostic genes screened from m6A-associated immune genes. The predictive power of the model was evaluated and our model was able to achieve good prediction.

Abnormal histone modifications (HMs) can promote the occurrence of breast cancer. To elucidate the relationship between HMs and gene expression, we analyzed HM binding patterns and calculated their signal changes between breast tumor cells and normal cells. On this basis, the influences of HM signal changes on the expression changes of breast cancer-related genes were estimated by three different methods. The results showed that H3K79me2 and H3K36me3 may contribute more to gene expression changes. Subsequently, 2109 genes with differential H3K79me2 or H3K36me3 levels during cancerogenesis were identified by the Shannon entropy and submitted to perform functional enrichment analyses. Enrichment analyses displayed that these genes were involved in pathways in cancer, human papillomavirus infection, and viral carcinogenesis. Univariate Cox, LASSO, and multivariate Cox regression analyses were then adopted, and nine potential breast cancer-related driver genes were extracted from the genes with differential H3K79me2/H3K36me3 levels in the TCGA cohort. To facilitate the application, the expression levels of nine driver genes were transformed into a risk score model, and its robustness was tested via time-dependent receiver operating characteristic curves in the TCGA dataset and an independent GEO dataset. At last, the distribution levels of H3K79me2 and H3K36me3 in the nine driver genes were reanalyzed in the two cell lines and the regions with significant signal changes were located.

Hippo pathway can regulate cell division, differentiation and apoptosis, and control the shape and size of organs. To study the distribution patterns of histone modifications of Hippo pathway genes in embryonic stem cells is helpful to understand the molecular regulation mechanism of histone modification and Hippo pathway on stem cell self-renewal. In this study, 19 genes of Hippo pathway including YAP, TAZ, LATS1/2, MST1 and SAV1, and eight histone modifications in embryonic stem cells were chosen to study the spatial distribution patterns of histone modifications. It was found that there were obvious type specificity and the location preference of target regions in the distributions of histone modifications, and H3K4me3 and H3K36me3 played the most important regulatory roles. Through the correlation analysis of histone modifications, a histone modification functional cluster composed of H3K4ac, H3K4me3, H3K9ac and H3K27ac was detected in YAP. In addition, the spatial distribution patterns of histone modifications in Hippo pathway genes were obtained, which provided a new theoretical reference for elucidating the mechanism of histone modifications regulating the gene expression of Hippo pathway, and for revealing the molecular regulatory mechanism of histone modifications affecting the self-renewal of embryonic stem cells by regulating the Hippo pathway.

Enhancer–promoter (E–P) interaction is an essential component of cis-regulatory regulation for gene expression. However, to comprehensively study the gene expression with the regulation of long-range E–P interactions is a major challenge in the regulatory networks. As these types of gene expression are regulated by diverse genomic signatures, we presented a computational method to study the relationships between gene expression levels and diverse genomic signatures. In this paper, based on the datasets of long-range E–P interactions, we extracted feature parameters from multiple signatures (e.g., epigenetic marks, transcription factors) and used regression models to predict the gene expression levels. In our results, we found that the predicted expression values correlated well with the measured expression values in both the interacting and non-interacting sets, and the correlation values of the interacting set were higher than that of the corresponding non-interacting set in each cell line, which indicated that the distal enhancers would cooperate with diverse genomic signatures to facilitate the expression level of target genes. By comparing the important signature features for the gene expression levels between the interacting and non-interacting sets in the same cell line, we found that the important specific signatures affect the gene expression regulated by distal enhancers. Our research provided additional insights about the roles of diverse signatures in gene expression with the regulation of distal enhancers.