AutEx: Regions of autozygosity
Description
Prerequisites
Dialog
Practical tips
Algorithm
Description
Identification of autozygous regions in an individual with consanguineous parents.
Prerequisites
The AutEx algorithm can be applied to any variant file loadable in Filtus, as long as the genotype format is correctly given in the input settings.
However, the following points should be kept in mind to maximize the power of AutEx:
- The more variants, the better! A typical single-sample file with 15-20 000 exonic variants is ok, but will leave some parts of the genome poorly covered.
If available, use files where intronic/intergenic variants are included.
- If the variant file has VCF-like genotype format, and contains 0/0 genotypes, remember to check the keep 0/0 box.
The 0/0 genotypes are just as informative for autozygosity mapping, as the 1/1 variants!
- Correct allele frequencies are important for this algorithm to perform optimally. Without a frequency column in your variant file,
you can still run AutEx (specifying a default frequency value) but the output will be less accurate.
Dialog
To open the AutEx dialog, first select a single sample in the Loaded samples window (at the top left of the main view of Filtus) and then choose AutEx: Regions of autozygosity in the Analysis menu. Any
filters (e.g. PASS) should be applied before opening the dialog.
- Model parameters
-
Use the parental relation parameter to indicate the degree of inbreeding of the sample under analysis.
The algorithm uses this as prior information on the number and length of autozygous regions it can expect.
In our experience, the results are not sensitive to this choice, in particular when the number of variants is large. This is not surprising:
Any true autozygous segment covered by the sequencing is likely to contain many variants. The collective evidence of these variants quickly drowns the impact of the prior.
- Allele frequencies
-
Indicate a column containing frequencies for
the ALT alleles. The Missing
entry value is substituted whenever the column does not
contain a number. If your variant file does not have frequency data,
the Missing entry value
will be used for all variants.
- Output
-
TODO
- Summary
-
A summary of the findings is printed here. The identified autozygous regions are shown in the main Filtus window; to inspect them you must close the AutEx dialog.
TODO: Explain inbreeding coefficient estimate
- Plot
-
This functionality plots a single chromosome, showing all the variants, the posterior probability curve, and the identified autozygous regions.
Practical tips
Tip 1: Before running AutEx it is recommended to apply filters removing variants of low quality (base/genotype/mapping).
Tip 2: The algorithm only works on the autosomes. You don't have to filter out X and Y yourself, this is taken care of internally.
Tip 3: To save the regions to a file, first close the AutEx dialog and then select Save main window content in the File menu.
Tip 4: After saving the regions to a file (previous tip) this file can be directly used as a Region filter.
This two-step approach is the common way to perform homozygosity/autozygosity mapping in Filtus.
Algorithm
todo