You can use the Identify Rare Disease Causing Mutations in a Family of Four (WGS) ready-to-use workflow to identify de novo and compound heterozygous variants from an extended family of four, where the fourth individual is not affected.
To run the Identify Rare Disease Causing Mutations in a Family of Four (WGS) workflow, go to:
Toolbox | Ready-to-Use Workflows | Whole Genome Sequencing () | Hereditary Disease () | Identify Rare Disease Causing Mutations in a Family of Four (WGS ()
- Double-click on the Identify Rare Disease Causing Mutations in a Family of Four (WGS) tool to start the analysis. If you are connected to a server, you will first be asked where you would like to run the analysis.
- The trimmed sequencing reads from the different family members are specified one at a time in successive dialogs (figure 20.41).
- Select which reference data set should be used to identify causal inherited variants (figure 20.42).
- Specify the Hapmap populations that should be used for filtering out variants found in Hapmap for each family member successively, with de novo and recessive being the proband variants (figure 20.43).
- Similarly, specify the parameters for the Fixed Ploidy Variant Detection tool for each family member successively (figure 20.44).
The parameters used by the Fixed Ploidy Variant Detection tool can be adjusted. We have optimized the parameters to the individual analyses, but you may want to tweak some of the parameters to fit your particular sequencing data. A good starting point could be to run an analysis with the default settings.
The parameters that can be set are:
- Required variant probability is the minimum probability value of the 'variant site' required for the variant to be called. Note that it is not the minimum value of the probability of the individual variant. For the Fixed Ploidy Variant detector, if a variant site - and not the variant itself - passes the variant probability threshold, then the variant with the highest probability at that site will be reported even if the probability of that particular variant might be less than the threshold. For example if the required variant probability is set to 0.9 then the individual probability of the variant called might be less than 0.9 as long as the probability of the entire variant site is greater than 0.9.
- Ignore broken pairs: When ticked, reads from broken pairs are ignored. Broken pairs may arise for a number of reasons, one being erroneous mapping of the reads. In general, variants based on broken pair reads are likely to be less reliable, so ignoring them may reduce the number of spurious variants called. However, broken pairs may also arise for biological reasons (e.g. due to structural variants) and if they are ignored some true variants may go undetected. Please note that ignored broken pair reads will not be considered for any non-specific match filters.
- Minimum coverage: Only variants in regions covered by at least this many reads are called.
- Minimum count: Only variants that are present in at least this many reads are called.
- Minimum frequency: Only variants that are present at least at the specified frequency (calculated as 'count'/'coverage') are called.
- Specify the affected sibling's gender (figure 20.45)
- In the last wizard step you can check the selected settings by clicking on the button labeled Preview All Parameters.
In the Preview All Parameters wizard you can only check the settings, and if you wish to make changes you have to use the Previous button from the wizard to edit parameters in the relevant windows.
- Choose to Save your results and click on the button labeled Finish.
The following outputs are generated:
- Read Mapping One for each family member. The reads mapped to the reference sequence.
- Filtered Variant Track One for each family member. The variants identified in each of the family members. The variant track can be opened in table view to see all information about the variants.
- Read Mapping Report One for each family member. The report consists of a number of tables and graphs that in different ways provide information about the mapped reads from each sample.
- De novo variants Variant track showing de novo variants in the proband. The variant track can be opened in table view to see all information about the variants.
- Recessive variants Variant track showing recessive variants in the proband. The variant track can be opened in table view to see all information about the variants.
- Identified Compound Heterozygous Genes Proband Gene track with the identified putative compound heterozygous Variants in the proband. The gene track can be opened in table view to see the gene names.
- Gene List with de novo Variants Gene track with the identified putative compound heterozygous Variants in the proband. The gene track can be opened in table view to see the gene names.
- Gene List with recessive Variants Gene track with the identified recessive variants in the proband. The gene track can be opened in table view to see the gene names.
- Track List This is a collection of tracks shown together in a view that makes it easy to compare information from the individual tracks, such as compare the identified variants with the read mappings and information from databases.
- De novo Mutations Amino Acid Track
- Recessive Variants Amino Acid Track