Genotype imputation without a reference panel
1
0
Entering edit mode
12 days ago
Mohadese • 0

I have genotyping array data, for some of the individuals we have missing data even on the genotyping array. We want to impute those missing data based on the rest of the dataset. What tools can I use?

genotyping-array-data reference-panel genotype-imputation • 350 views
ADD COMMENT
0
Entering edit mode

Be aware the individual may have good reason to not get a call at that site. Maybe the gene is deleted or duplicated and has no clear diploid genotype. Imputation will put in a likely result, but not the sure result.

ADD REPLY
0
Entering edit mode
12 days ago
Pine ▴ 20

I might be wrong but I am not sure if using the same genotyping data internally would work in a classic imputation setting. My impression is that you need a separate reference panel for imputation beside the genotype data you want to analyze. If lack of reference panel is an issue, some online servers provide preloaded reference panels, for instance:

Michigan imputation server: https://imputationserver.sph.umich.edu/index.html#!

Guidelines here: https://imputationserver.readthedocs.io/en/latest/

I suspect you might be able to fill in the missing dosage values with population average for instance (how the missing values are sometimes treated for non-genetic variables) using packages like SNPStat in R which give you full access to genotype matrices. But I am not personally familiar with studies that have used this approach for genotyping data.

ADD COMMENT
0
Entering edit mode

Unfortunately we are not able to use online servers like Michigan or TOPMed.

ADD REPLY
0
Entering edit mode

Many genotype reference panels and imputation packages are publicly available [including the ones used by Michigan imputation server] and can be downloaded for a local imputation. My understanding is that local imputation is also a common practice, but unfortunately I do not have personal experience with it.

ADD REPLY

Login before adding your answer.

Traffic: 1179 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6