HomozygosityMapper

Documentation

HomozygosityMapper stores all data in a database. To upload data, you first have to create a profile (i.e. to register).
It is not possible to choose a user name that is already in use, HomozygosityMapper will raise an error in such cases. You should always enter a valid e-mail address because we will inform you in advance when major changes to the website are planned or when we have to archive old projects due to disk space limitations.
If you forgot your password, please send us an e-mail.

After you have created a password, you will be automatically logged in with these data. If you return to HomozygosityMapper, please click on the login link on the homepage and enter your credentials.

Genotypes from the following chips can be imported:

• Affymetrix 10K Xba142
• Affymetrix 50K Hind240
• Affymetrix 50K Xba240
• Affymetrix 250K_Nsp
• Affymetrix 250K_Sty
• Affymetrix GenomeWideSNP5
• Affymetrix GenomeWideSNP6
• Illumina (all SNP chips)

To import the genotypes, specify a new project name under new project name or select one of your existing projects (add to existing project) if you want to add new genotypes. Adding genotypes to an existing project is also the solution for projects run on chip sets consisting of more than one chip (see 'mixing' below). You can also decide whether or not you want to make the project public (checkbox access restriction) - this setting can be changed at any time.

After choosing the project, you have to select the respective genotyping chip (drop-down menu chip), point to the genotype file (genotype file) and eventually submit the data.

Genotype files must be tabular with the samples as columns and the SNPs as rows, they can (or rather should) also be zipped or gzipped. Lines starting with the number sign (#) will be ignored. Sample files to demonstrate the allowed formats can be found here. The genotypes must be written in standard Affymetrix or Illumina notation, i.e. (NoCall, AA, AB, BB or -1, 0, 1, 2 for Affymetrix and --, AA, AB, BB for Illumina, respectively).

The data import may take seconds or even half an hour, depending on the size of the file, the upload speed and the working load of our server.

When the upload has been finished, you can proceed to the next step.

Mixing of arrays with different markers per sample is possible. However, HomozygosityMapper only considers markers that were genotypes for a sample. That means, that if you mix different arrays you will automatically obtain lower homozygosity scores for samples with fewer genotypes.
We believe that any other solution (using block length in megabases instead of SNP numbers, treating non-typed markers as unknown genotypes, or use only markers available in all samples) would merely complicate the application and make it less performant in more realistic scenarios. If, however, you are in such a situation, please drop us an e-mail. At attempt to solve this might be worthwile, challenging and exciting... :-)

After the genotypes were imported, you can start the homozygosity mapping. First, you must select the project you want to analyse (drop-down menu project) and provide a new for the analysis (analysis name). You can specify further details under analysis description which will be displayed together with the analyses' names to indicate specific models etc. which might be too complex to include them in the name.

Under cases and controls, you can specify the samples that should be regarded as cases or controls, respectively. Cases are, of course, obligatory. Please note that sample IDs must be written exactly as in the genotype file; HomozygosityMapper will, however, complain, if it cannot find a sample.

If your cases are from a single family or if you are absolutely sure that there is no genetic heterogeneity among the samples you plan to analyse as cases, you can tick the require genetic homogeneity box. In this case, homozygosity in all samples is required for a region to be included.

If you are interested in the excess of homozygous over heterozygous genotypes, you can change the allele frequencies drop-down menu to compare homozygosity in affected indivuals with homozygosity in your controls or one of the four HapMap populations. These frequencies are only plotted in the graphical output and will not affect the analysis in any way, though.

The last parameter to be set is the maximum block length (limit block length). If this field is left empty, a default setting for the set of markers will be used. If you don't want to limit the block length, specify a negative value. While the default settings are suitable for most analyses, after a visual inspection of the results you may decide to re-analyse your data according to the degree of homozygous blocks in your data. A small value might be appropriate if you want to find a founder mutation in a very large sample consisting of individuals from non-related families. Default values are:

number of SNPs    block length limit 
    > 800000         1000
    > 400000          500
    > 200000          250
    >  45000           80
    >      0           15  	

The analysis might take several minutes to complete.

This interface allows you to select one of your own analyses (coloured blue) to inspect homozygosity scores and genotypes. In addition, analyses of other researchers made public (green) or accessible for you (black) will be listed.

Simply click on the right analysis to continue.

Here, the homozygosity scores are plotted against the physical position. The interface can either display the whole genome, a single chromosome or a selected region on a single chromosome. Interesting regions are indicated by red bars.
Above the bar chart, small red and blue bars are shown. These represent the excess (red) or the shortage (blue) of the found homozygosity against the expectation calculated from controls or population values as selected when the data was analysed. This view becomes more interesting when
Below the bar chart, all scores above 80% of the maximum are listed (sorted by their score) and direct links to these regions are provided. On top of this table, 'broad' regions are displayed - in these, smaller decreases of the score within a homozygous region are neglected. At the bottom of the table, 'narrow' regions with sharp limits follow. If you expect some degree of heterozygosity, you should consider the 'broad' regions.

genome-wide homozygosity

In the whole genome view, clicking on a chromosome will zoom on this chromosome.

chromosomal homozygosity

After clicking on a chromosome, the view zooms into this chromosome. Now only the interesting regions located within the selected region (here: the chromosome) will be displayes as links.

chomosomal homozygosity - zoom in

After clicking on the edges of the region the user want to study in details, a pop-up will appear and offer links to further zoom in or to see the underlying genotypes.

The genotype view depicts the single genotypes of all samples marked as cases or controls. The genotypes are indicated as coloured boxes; where blue codes for heterozygous, grey for unknown and red for homozygous genotypes. Different shades of red are used, longer homozygous stretches will be drawn in a 'deeper' red than single homozygous markers. Genotypes homozygous for the minor allele are marked with a black diagonal bar.

The homozygous region found by HomozygosityMapper is shown as a rectangle, it's limits can be shifted to other markers by clicking on these.

After the region has been defined, you can go ahead by clicking on the GeneDistiller button to find candidate genes (see below) within the region.

HomozygosityMapper integrates with GeneDistiller. When clicking the GeneDistiller button, GeneDistiller's query interface will open with the correct region filled in. Clicking on submit in this interface will list all genes within the region; for a description how GeneDistiller can be used to focus on genes of interest take a look at it's manual:
http://www.genedistiller.org/manual.html

When GeneDistiller is queried from HomozygosityMapper, each gene will carry another link, 'HomozygosityMapper', which will open the genotypes view and indicate the position of the selected genes within the genotypes. Bookmarking GeneDistiller's settings (on the bottom of it's results list) will maintain the hyperlinks to the respective homozygosity mapping.

You can use this interface to fine-tune access to your data. It will display a table with all of your projects as rows and all user logins as columns. For each user and project, you can select whether he or she shall be allowed to query your data. You can further grant the right to perform new analyses on your data to collaborators (the project and the analyses will remain under your control though). At the end of each row, you have the option to grant (or revoke) query permissions to the public.

Click on SetPermissions when you're done.

Use this interface to delete single analyses within a project or complete projects.

All of your projects and analyses will be listed. Checking a project/analysis will mark it for deletion; if you delete a project, all analyses performed with that dataset will be deleted as well.

Please note that deleted data will be lost forever.

HomozygosityMapper is based on a PostgreSQL database, all interfaces were written in Perl. The database schema and further details can be found in the technical documentation.

The application was developed using various flavours of Mozilla Firefox. It was further tested with Microsoft Internet Explorer 7 and should work

In addition to the web-based version, a stand-alone application exists. This program does not make use of a database and interfaces were created with Perl/Tk. While it offers speed gains for small chips (up to 50,000 SNPs), the performance decreases with the number of markers being studied. This version is not supported anymore but we will gladly share the source code. Please send us an e-mail if you are interested.

If you are unsure which medium best suits your contact needs: E-Mails are preferred.

Why can't I further zoom in on a single chromosome? Why can't I change the limits of a homozygous region?
Sounds as if JavaScript was not enabled. With Mozilla Firefox, select Tools -> Options -> Content and check the 'Enable JavaScript' button.

I forgot my login and/or my password. What can I do?
Please send an e-mail to dominik.seelow(at)charite.de.

I am using a chip that is not listed, how can I upload my genotypes?
So far, HomozygosityMapper only supports the most common chips manufactured by Affymetrix and Illumina. If you want to import genotypes derived by other chips (or means), you can either bring your genotype file in the format as used by Affymetrix or Illumina and select a chip containing all your markers or write us an e-mail. The integration of completely new chips and file formats might take some time, though.

I made a terrible mistake and deleted my project. Is there any way to restore it?
Well, yes, at least there might be one. If you are lucky, the data is included in one of our backups. Please write an e-mail to us and specify the exact name of the project. We will then see what we can do.

Is the data within HomozygosityMapper stored forever?
Hm. Certainly not forever but at least we do not plan to include an automatic time-out after which your data will be deleted. We might, for disk space reasons, decide to archive 'old projects'. However, we will inform the respective project's owners in advance.

I have created an account but I cannot log in. Why?
Please ensure that you allowed this site to set cookies in your browser. This is the default setting in all current browsers - if you deactivated cookies you could change the setting to 'ask' instead of 'block' or add an exception for this website.