This webpage provides details on the following subjects.

Collinear sets

Hit popup

3-chr conserved

Collinear sets

Either using the Hit Filter, or right click in the hit region, set the highlighting for Collinear Sets. As shown below, the number of highlighted collinear hits and sets will be shown in the Information box (CoSets is the abbreviation for Collinear Sets).

In the image on the right, two different sets are shown where the pink highlighted hits is the 1st set, green is the 2nd; these are Color Wheel Highlight1 and Highlight2 colors. The sets were broken by an extra gene on Human Chr01. A collinear set either has all hits/genes to the same strand (+/+, -/-) or different (+/-, -/+). In the image on the right, all genes in both sets are on different strands. Note: The above is not guaranteed when using Algo1 (see Terminology), as the hit strands may be different than the strands assigned to the homologous genes. The collinear set algorithm does not consider the amount of overlap of the hit to a gene, or the similarity of the hit sequences. Sometimes a hit looks like it overlaps a gene, but the gene is actually in the gap of a clustered hit, so is not a gene hit.

The image below shows a collinear set of 4 genes and 6 brown hit-wires that are not part of the set. The following explains them, using the symbols g0 (hit 0 genes), g1 (hit 1 gene), g2 (hit paired genes).

1 g1 Ignored, not part of the set since it does not hit paired genes. 2 g0 Ignored. 3 g1 Ignored since the same gene has a hit to a paired gene. 4 g0 Ignored. 5 g1 This ends the set since it does not hit a paired gene. 6 g2 This could start a new set, but is not followed by another paired gene. Sometimes two genes look as one, so the Sequence Filter Gene Delimiter option is important to view these.

Hit popup

Go to top

Default view

Merge

Order

Default view

• Both lists are sorted by start coordinate.
• The # column in the two tables align to each, e.g. the rows #2 align to each other.
• The project Directory name that is alphabetically lower (e.g. Arab<Cabb) is the query and the other is the target. The query is numbered 1-N; the target is ordered to match the query.
• Two subhits may overlap (Gap<=0) on one chromosome but not the other.

Merge

In the Merged view, hits with Gap<=0 are merged with the hit they overlap with. The numbered hits are no longer 1-to-1, i.e. they do not necessarily align to each other like they do in the un-merged view. Totals are provided for the Len and Gap. The Merged view corresponds better to the visual view on the 2D track. The existence of the Merge button is a quick indication that there are overlapping hits.

Order

As stated above, each numbered subhit corresponds to the same number in the opposite list, e.g. the two #1 subhit aligns. An # on the far right implies that the order numbers are not sequential. The Order option produces another popup where the target list is sorted by the # column, and the coordinates are no longer sorted. The existence of the Order button is a quick indication that there are disordered subhits.

3-chr conserved

Chromosome Explorer

Queries

Chromosome Explorer

The Dot plot can show 3 chromosomes, as illustrated on the upper below (click image for a better view).

The 2D can show a 3-track view, as illustrated on the upper right and discussed below.

2D: The Sequence Filters for the reference track will have an extra section of filters as shown below.

g2x2 The genes conserved across the 3-chromosomes are highlighted along with the hit-wires. g2x1 The genes conserved on ONLY one chr-pair are highlighted along with the hit-wires. Show Only the g2x2 or g2x1 hit-wires will be shown.

The g2 implies the hit has a gene on both ends, and the x2 indicates it exists for both sets of chromosome pairs, the x1 indicates it exists for just one. When either g2x2 or g2x1 is selected, this algorithm is run to compute the values; it is only computed on the visible set.

There are two issues confusing the results:

• One side g2 filter:
  • The algorithm runs on the set of filtered hits. A g2 hit to a reference gene may be filtered in one chr-pair and not the other, making it appear to be a g1 hit.
  • In this case, the filtered hit will be forced displayed using the Highlight2 color (default pink).
  • For example: the lower middle image where the arrow is pointing.
  • If both sides are filtered, neither will be shown.
• g1 possible g2:
  • Both chr-pairs can have a hit to the same reference gene, where one chr-pair has a g2 hit-wire and the other has a g1 (but a possible g2 if it is missing the gene annotation).
  • In this case, the g1 hit-wire will be displayed in red.
  • For example: the lower right image where the arrow is pointing.
  • If the g2x1 hit is filtered, the adjacent g1 hit-wire will not be highlighted red.

The Sequence Filter of Show Gene# with Hit is used in the above display.

Caveats of display: Nothing is refreshed or recalculated!

• Suggestion: Set the filters from the "full sequence" view, then you can zoom in and back out and the highlight will remain.
• If the highlighting looks wrong after zooming, open Sequence Filter, click None and re-click g2x2 or g2x1.

Caveats of algorithm:

• For g2x2, say it shows genes A-B-C are conserved. The algorithm only checks the A-B hit and B-C hits exist, but does not check the A-C hit.
• There may be different numbers of highlighted genes on the different tracks or hit-wires between tracks. This is because multiple hits can align to a gene, and one gene can align to two different opposite genes.
• When there are overlapping genes, a hit is assigned to just one of them, which causes some conserved genes to have incorrect pairing.
• There are occasional tiny hits that just barely overlap genes; the current algorithm marks them as conserved genes. Also, the hits may align strictly in introns (especially for the long introns of mammalian genomes); the current algorithm marks these as conserved.

Queries: There is no query filters that give equivalent results.

Queries

→ For a query result, Report is only way to view more than two species in a query row.
→ As shown here, it is possible to bring up the 2D display with 3 chromosomes from the query table.

Select two rows with the following requirements: (1) Each row must have two genes. (2) There must be a shared gene. Different species Select View 2D with either the Region and Collinear options. The image on the right used the Collinear options. Same species Atha will be the reference since Gene# 5468 is listed in both rows, and Brap chromosomes 7 and 9 will be on its left and right (this 2D display is not shown).

Go to top

Email Comments To: cas1@arizona.edu