Dataset Information

Analysis of mRNA Expression and DNA Methylation Datasets According to the Genomic Distribution of CpG Sites in Osteoarthritis.

ABSTRACT: Objectives Transcriptional changes in cartilage can impact function by causing degradation such as that which occurs during the development of osteoarthritis (OA). Epigenetic regulation may be a key factor leading to transcriptional changes in OA. In this study, we performed a combined analysis of DNA methylation and gene expression microarray datasets and identified key transcription factors (TFs) central to the regulation of gene expression in OA. Methods A DNA methylation profile dataset (GSE63106) and a gene expression profiling dataset (GSE114007) were extracted from the Gene Expression Omnibus (GEO). We used ChAMP methylation analysis and the Limma package to identify differentially methylation genes (DMGs) and differentially expressed genes (DEGs) from normal and human knee cartilage samples in OA. Function enrichment analysis of DMGs was conducted using the DAVID database. A combined analysis of DEGs and DMGs was conducted to identify key TFs in OA. We then validated the mRNA expression of selected TFs in normal and OA cartilage by RT-qPCR. Primary chondrocytes were cultured and treated with the DNA methylation inhibitor 5-Aza-2-deoxycytidine (5-Aza) for functional validation. Results We identified 2,170 differential methylation sites (DMS) containing 1005 genes and 1986 DEGs between normal human and OA cartilage. Functional analysis of DMGs revealed that focal adhesion, extracellular matrix (ECM)-receptor interactions, the PI3K-Akt signaling pathway, and the FoxO signaling pathway were involved in OA. Integrated analysis showed a subset of 17 TFs. Four TFs (ELF3, SOX11, RARA, and FOXD2) were validated. RT-qPCR results showed the mRNA expression of SOX11, RARA, and FOXD2 were consistent with the results from the mRNA expression data. However, the expression of ELF3 could not be validated. Upon 5-Aza-2'-deoxycytidine (5-Aza) treatment, the mRNA levels of ELF3 and SOX11 were down-regulated, whilst RARA was up-regulated, and FOXD2 showed no significant change in expression level. Conclusions the effect of DNA methylation on the transcriptional regulation is related to the distribution of methylated sites across the genome. Epigenetic studies on the positions of DMS in transcriptional units can inform a better understanding of the function of DNA methylation and its transcription regulation.

SUBMITTER: Yi P

PROVIDER: S-EPMC8082497 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Genetics and epigenetics are important subjects in the field of osteoarthritis (OA) research. DNA methylation may affect gene transcription, but the specific mechanisms have remained to be fully elucidated. In the present study, the ChAMP methylation analysis package was used to identify differentially methylated genes (DMGs) from the dataset GSE63695 from the Gene Expression Omnibus (GEO) database. The distribution of differentially methylated sites (DMS) and the total array sites across the genome were analyzed by enrichment analysis. Subsequently, two mRNA expression profiling datasets, GSE114007 and GSE113825, were obtained from the GEO database and common differentially expressed genes (DEGs) were identified using the Limma package. Key genes were screened by analyzing the distribution of DMS across the genome consisting of DEGs and DMGs. A total of 1,662 and 1,986 DEGs were identified between OA and normal human cartilage from the GSE113825 and GSE114007 dataset, respectively. A further screening revealed 292 genes with common differences between the two datasets. A total of 574 DMS containing 394 DMGs were observed between OA and normal cartilage. Integrative analysis revealed a corresponding subset of 15 genes. Of these, 6 genes were verified by reverse transcription-quantitative PCR, confirming that the mRNA expression of 5 genes (MAP1B, FNDC1, ANLN, SCNN1A and STC2) in OA cartilage was consistent with the mRNA expression from the analysis of the datasets. Upon treatment with the DNA methylation inhibitor 5-aza-2'-deoxycytidine, the mRNA levels of FNDC1 and SCNN1A were decreased, and no significant alteration in the mRNA levels of MAP1B, ANLN, KCNN4 and STC2 was observed. The incidence of differential methylation varied in subregions of the genome and the effects on transcription were associated with the distribution of DEGs across the genome. The regulation of this appears more complex than initially postulated. Combining the data on epigenetic differences of OA with the genome or transcriptome data for analysis may improve the understanding of the pathophysiological processes of OA. FNDC1 and SCNN1A may potentially be valuable biomarkers for OA.

Project description:BACKGROUND: The Human Papillomavirus (HPV) genome is divided into early and late coding sequences, including 8 open reading frames (ORFs) and a regulatory region (LCR). Viral gene expression may be regulated through epigenetic mechanisms, including cytosine methylation at CpG dinucleotides. We have analyzed the distribution of CpG sites and CpG islands/clusters (CGI) among 92 different HPV genomes grouped in function of their preferential tropism: cutaneous or mucosal. We calculated the proportion of CpG sites (PCS) for each ORF and calculated the expected CpG values for each viral type. RESULTS: CpGs are underrepresented in viral genomes. We found a positive correlation between CpG observed and expected values, with mucosal high-risk (HR) virus types showing the smallest O/E ratios. The ranges of the PCS were similar for most genomic regions except E4, where the majority of CpGs are found within islands/clusters. At least one CGI belongs to each E2/E4 region. We found positive correlations between PCS for each viral ORF when compared with the others, except for the LCR against four ORFs and E6 against three other ORFs. The distribution of CpG islands/clusters among HPV groups is heterogeneous and mucosal HR-HPV types exhibit both lower number and shorter island sizes compared to cutaneous and mucosal Low-risk (LR) HPVs (all of them significantly different). CONCLUSIONS: There is a difference between viral and cellular CpG underrepresentation. There are significant correlations between complete genome PCS and a lack of correlations between several genomic region pairs, especially those involving LCR and E6. L2 and L1 ORF behavior is opposite to that of oncogenes E6 and E7. The first pair possesses relatively low numbers of CpG sites clustered in CGIs while the oncogenes possess a relatively high number of CpG sites not associated to CGIs. In all HPVs, E2/E4 is the only region with at least one CGI and shows a higher content of CpG sites in every HPV type with an identified E4. The mucosal HR-HPVs show either the shortest CGI size, followed by the mucosal LR-HPVs and lastly by the cutaneous viral subgroup, and a trend to the lowest CGI number, followed by the cutaneous viral subgroup and lastly by the mucosal LR-HPVs.

Dataset Information

Analysis of mRNA Expression and DNA Methylation Datasets According to the Genomic Distribution of CpG Sites in Osteoarthritis.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets