Skip to content

Modified methylation_sites resource to point to a resource with merged methylation data with autosomes and sex chromosomes and updated transform_grch38_methylation accordingly.#801

Merged
jkgoodrich merged 1 commit intomainfrom
jg/create_sex_chr_methy_resources
Sep 11, 2025
Merged

Modified methylation_sites resource to point to a resource with merged methylation data with autosomes and sex chromosomes and updated transform_grch38_methylation accordingly.#801
jkgoodrich merged 1 commit intomainfrom
jg/create_sex_chr_methy_resources

Conversation

@jkgoodrich
Copy link
Contributor

This is to add the methylation data for sex chromosomes that Ruchit gave us.

The current file is located here: gs://gnomad/v4.1/constraint/resources/methylation_all.ht

If it looks good I will move it to the staging location.

The file was created like this:

from gnomad.resources.grch38.reference_data import methylation_sites

methyl_autosomes = methylation_sites.ht()
methyl_chrx_par = hl.read_table("gs://gnomad/v4.1/constraint/resources/methylation_chrX_par.ht")
methyl_chrx_nonpar = hl.read_table("gs://gnomad/v4.1/constraint/resources/methylation_chrX.ht")
methyl_chry_nonpar = hl.read_table("gs://gnomad/v4.1/constraint/resources/methylation_chrY.ht")

methyl_chrx_par = methyl_chrx_par.rename({"methy_level": "methylation_level"}).select("methylation_level")
methyl_chrx_nonpar = methyl_chrx_nonpar.rename({"methy_level": "methylation_level"}).select("methylation_level")
methyl_chry_nonpar = methyl_chry_nonpar.rename({"methy_level": "methylation_level"}).select("methylation_level")

methyl_ht = methyl_autosomes.union(methyl_chrx_par, methyl_chrx_nonpar, methyl_chry_nonpar)
methyl_ht = methyl_ht.naive_coalesce(1000).checkpoint("gs://gnomad/v4.1/constraint/resources/methylation_all.ht", overwrite=True)

…ged methylation data with autosomes and sex chromosomes and updated `transform_grch38_methylation` accordingly.
Copy link
Contributor

@ch-kr ch-kr left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

one minor comment

@jkgoodrich jkgoodrich merged commit 7a749c6 into main Sep 11, 2025
6 checks passed
@jkgoodrich jkgoodrich deleted the jg/create_sex_chr_methy_resources branch September 11, 2025 20:15
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants