This is a test version of Biostars. For the public version, visit https://www.biostars.org.
get a set of reliable SNPs when knowing parent and child relationship

I have a set of SNP data in vcf format and have information regarding parent and offspring relationship. I want to use this kinship relationship to subset SNPs to a smaller number to ensure those SNPs are really reliable. Could someone tell me if there is some existing tools to do so?

snp identity-by-descent
• 2,748 views • link
updated by Ram Ram • written by lilepisorus • lilepisorus

Hi, Pierre

Could you follow up with this question? How to subtract variants in violations.vcf from all variants? Any existing tools? I checked the -selectVariant in GATK, it seems that it can go for the common or shared ones, not the unique ones.

Thanks
Li

• 0 views • link
updated by Ram Ram • written by lilepisorus • lilepisorus

2 answers

https://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_gatk_tools_walkers_variantutils_SelectVariants.php "Generating a VCF of all the variants that are mendelian violations:"

java -Xmx2g -jar GenomeAnalysisTK.jar \
   -R ref.fasta \
   -T SelectVariants \
   --variant input.vcf \
   -ped family.ped \
   -mvq 50 \
   -o violations.vcf

I wrote a tool https://github.com/lindenb/jvarkit/wiki/VCFTrio to annotate a VCF with mendelian incompatibilities in trios.

• 0 views • link
updated by Ram Ram • written by Pierre Lindenbaum Pierre Lindenbaum

Hi, Pierre

I tried both ways you suggested.

The SelectVariants in GATK gave me the full set of SNPs in my input, suggesting all SNPs violated the Mendel heritability rule.

I also tried VCFTrio, in the result file, I saw comments of tag "MENDEL" in the header, but for each SNP, I did not see "MENDEL" in the INFO field. So I cannot filter sites which violate Mendelian rules based on the INFO field.

I am wondering if it is true that all of sites violate or I did anything wrong in my command. My SNP set is from SNP calling for the whole genome in GATK. So, I doubted none of the SNPs obey the Mendelian Rules based on the provided pedigree. Here is my command for both analyses:

java -Xmx240g -jar GenomeAnalysisTK.jar \
  --pedigreeValidationType SILENT \
  -R Glycine_max.V2.fasta \
  -T SelectVariants \
  --variant:VCF combined.REDUCED.1.realign.123013.uniq.sorted.SNPsOnly.vcf \
  -ped soybean.ped \
  -mvq 50 \
  -o violations.vcf

java -jar ./jvarkit/dist-1.128/vcftrio.jar \
  -p soybean.ped.vcftrio.txt \
  ../vcfFiles/combined.REDUCED.1.realign.123013.uniq.sorted.SNPsOnly.vcf > vcftrio.result
• 0 views • link
updated by Ram Ram • written by lilepisorus • lilepisorus

It's possible that the pedigree tools are making assumptions that the reference is a human genome, in order to understand the inheritance patterns? (E.g. in the RTG solution, the reference.sdf should contain a configuration file that specifies the autosomes and the sex chromosome inheritance patterns for the reference genome of interest, so that these constraints can be applied appropriately to the individual samples during variant calling and mendelianness checking).

I'm not familiar with the soybean genome -- is it diploid (which may be an assumption of the mendelianness tools) or polyploid?

• 1 views • link
written by Len Trigg Len Trigg

Here's another way, which outputs just the consistent variants, using Real Time Genomics:

rtg mendelian -t reference.sdf --pedigree family.ped -i input.vcf \
    --output-consistent goodvariants.vcf.gz

There is also --output-inconsistent if you want a file with just the violations, or regular --output if you want the mendelian consistency status just added as annotations. It may help to use --lenient if the variants were not jointly called.

• 0 views • link
updated by Ram Ram • written by Len Trigg Len Trigg

Log in to answer this question.

More posts like this