Thank you, but I have more intergenic regions in my starting list.
• 0 views
•
link
written
by
arsala521 •
Getting human-specific regions of the human genome
Hi everyone,
I have a list of human genomic regions (bed file enlisting genomic coordinates). From that list I want to extract those regions which are human-specific. I want to ask what can be appropriate pipeline for that.
I found one study using this pipeline: Using liftOver utility for converting the coordinates to mouse genome (mm9), then those not converted to the marmoset genome (CalJac3), and again those not converted to the chimp genome (PanTro2) to finally get human-specific regions.
Should I follow this, or could there be some alternate or more efficient way?
Thanks in advance
• 1,975 views
•
link
updated
by
Alex Reynolds
•
written
by
arsala521 •
•
written
by
arsala521 •
1 answer
Starting with orthologous genes ("orthologs"), you might take the subset of human genes which do not overlap or contain those orthologs, by some threshold of base or percentage overlap.
• 0 views
•
link
written
by
Alex Reynolds
Perhaps perform a BLAST nucleotide search between genomes. High-scoring regions would presumably share similarity and may be orthologous. You could also incorporate phyloP or phastcons conservation scores over your intergenic regions, looking for lower-than-background scoring regions as further evidence of speciation among primate alignment.
• 0 views
•
link
written
by
Alex Reynolds
Okay, thank you for the guide. I will look into these details.
• 0 views
•
link
written
by
arsala521 •
Log in to answer this question.
More posts like this
-
Subtracting bed file from whole genome
written by arsala521 •
-
How to obtain genomic coordinates of novel transcripts?
written by nattzy94 •
-
Convert variant coordinates from one genome build to the other
written by Gene_MMP8
-
Methylation Status of CpG islands for a specific gene
written by kspata
-
human disease specific gene
written by rajesh
-
SNP coordinate systems in UCSC LiftOver (hg19), NCBI Remap (GRCh37(hg19)) and dbSNP (GRCh38(hg38))
written by curious •
-
How does liftover browser rank hits of coordinates?
written by jrxu.bioinf •
-
Small genomic region proximity enrichment
written by Nicolas Rosewick
-
liftover a fasta file
written by diviya.smith •
-
Any Method Of Converting Bigwig File Format Into Bed Format?
written by samwise.techno •
I think you are looking at mappability to identify unique regions. You can take a look at this recent paper.
Thank you. Just looked into this paper. What I understood is, it is about identifying unique regions within a single genome.
This is a little tricky because of the "orthology" That Alex mentions below. There will always be some sequence similarity between humans and our close relatives (monkeys, mice etc) so there may be no true
human-specificregions that are coding.Something like 96% of human and chimpanzee sequence is identical, by one measure. This may be a very difficult problem, without more detail in the question.
Thank you for the guide. I will look into more details.