<p>I'm not sure how to generate the additional mutations, but I would say that HGMD is not the way to find all the pathogenic variants. I would probably filter <a href="http://www.ensembl.org/Homo_sapiens/Gene/Variation_Gene/Table?db=core;g=ENSG00000108055;r=10:110567695-110606048" rel="nofollow noreferrer">this table</a> by either Clinical significance or Evidence-&gt;Phenotype association.</p>

CC BY-SA 4.0
2021-01-12T13:31:39.163Z
2021-01-12T13:31:39.163Z
Emily_Ensembl
<p>I would like to create a list (probably .vcf format would be good) of all known pathogenic missense mutations in a human gene of interest and then add other variants that could lead to the same pathogenic amino acid substitutions. I intend to search targeted sequencing data of my gene of interest for variants that would cause these particular substitions.</p>
<p>I'm very new to bioinformatics so I'm not so familiar with the appropriate databases and tools to use, but this is a rough idea of how it could be done:</p>
<ol>
<li>extract all mutations listed in Ensembl within my gene of interest that are labelled with <code>HGMD_MUTATION</code>  (is there an option to download a vcf file with all HGMD mutations within my gene? I've only managed to do so for individual mutations.)</li>
<li>annotate list with reference and mutant amino acids for every mutation in the list</li>
<li>add additional mutations (which may not be present in HGMD) which could lead to the same substitution. E.g. I'm interested in the substitution Proline (reference codon: CCG) to Leucine (TTR or CTN) at AA position 102, so I would like to list all hypothetical mutations that would create a leucine codon in this position, not only the one given in Ensembl (which is CCG -&gt; CTG).</li>
</ol>
<p>Could you point me towards the best tools to use for these purposes? Unfortunately I haven't got a HGMD subscription. Many thanks!</p>

CC BY-SA 4.0
2021-06-30T06:02:02.673Z
Why HGMD specifically (not clear to me)? You can simply use the ClinVar database to find all pathogenic mutation within regions of interest/gene (you can check the genomic coordinates of any gene using either UCSC or ensemble).
Lot_to_learn
<vcf><genomics><ensembl><mutations>
How do I create a VCF file of all known pathogenic mutations in a gene of interest?
# How do I create a VCF file of all known pathogenic mutations in a gene of interest?

(Asked by: username_0 on 2021-01-10)

I would like to create a list (probably .vcf format would be good) of all known pathogenic missense mutations in a human gene of interest and then add other variants that could lead to the same pathogenic amino acid substitutions. I intend to search targeted sequencing data of my gene of interest for variants that would cause these particular substitions.

I'm very new to bioinformatics so I'm not so familiar with the appropriate databases and tools to use, but this is a rough idea of how it could be done:

1. extract all mutations listed in Ensembl within my gene of interest that are labelled with `HGMD_MUTATION` (is there an option to download a vcf file with all HGMD mutations within my gene? I've only managed to do so for individual mutations.)
2. annotate list with reference and mutant amino acids for every mutation in the list
3. add additional mutations (which may not be present in HGMD) which could lead to the same substitution. E.g. I'm interested in the substitution Proline (reference codon: CCG) to Leucine (TTR or CTN) at AA position 102, so I would like to list all hypothetical mutations that would create a leucine codon in this position, not only the one given in Ensembl (which is CCG -> CTG).

Could you point me towards the best tools to use for these purposes? Unfortunately I haven't got a HGMD subscription. Many thanks!

**Tags:** vcf, genomics, ensembl, mutations
> **Comments:**
> **username_1**: Why HGMD specifically (not clear to me)? You can simply use the ClinVar database to find all pathogenic mutation within regions of interest/gene (you can check the genomic coordinates of any gene using either UCSC or ensemble).  

---

## 1 answers

### Answer by username_2 on 2021-01-12. Score: 1

I'm not sure how to generate the additional mutations, but I would say that HGMD is not the way to find all the pathogenic variants. I would probably filter [this table](http://www.ensembl.org/Homo_sapiens/Gene/Variation_Gene/Table?db=core;g=ENSG00000108055;r=10:110567695-110606048) by either Clinical significance or Evidence->Phenotype association.


---

