AutEx: Regions of autozygosity

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