DNA Methylome and RNA-seq intersect
Entering edit mode
2.9 years ago
Morris_Chair ▴ 360

Hello everyone, I have a methylome dataset and a RNA-seq dataset and I want intersect each other in order to define the methylation status of differentially expressed genes of my interest.

To define the promoter region of all genes in the genome I used this code:

allproms <- promoters(TxDb.Hsapiens.UCSC.hg38.knownGene, downstream = 0, upstream = 2584)

but I’m not sure about the range to give to extract the informations,

My objective is to find out which promoters have methylated bases within their range. I am attempting to do this with the following lines of code:

metproms <- join_overlap_inner(allproms, metdat)

Where allproms is the object containing the coordinates and ranges of the promoters, and metdat the object with the ranges for the methylated bases.

Is there a better way to do that?

Thank you

RNA-seq Methylome DNA • 678 views
Entering edit mode
2.9 years ago
Tian ▴ 50

Based on my shallow experience, you may try:

  • FEM, which statistically integrates: RNA-seq, Methylation, and PPI network.
  • MOFA: Multi-Omics Factor Analysis. I used it for multi-omic clustering only.

I will encounter this problem like you soon in next month, so I am also looking for better tools or pure algorithms.


Login before adding your answer.

Traffic: 2306 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6