FIRMA: a method for detection of alternative splicing from exon array data.

MOTIVATION: Analyses of EST data show that alternative splicing is much more widespread than once thought. The advent of exon and tiling microarrays means that researchers now have the capacity to experimentally measure alternative splicing on a genome wide level. New methods are needed to analyze the data from these arrays. RESULTS: We present a method, finding isoforms using robust multichip analysis (FIRMA), for detecting differential alternative splicing in exon array data. FIRMA has been developed for Affymetrix exon arrays, but could in principle be extended to other exon arrays, tiling arrays or splice junction arrays. We have evaluated the method using simulated data, and have also applied it to two datasets: a panel of 11 human tissues and a set of 10 pairs of matched normal and tumor colon tissue. FIRMA is able to detect exons in several genes confirmed by reverse transcriptase PCR. AVAILABILITY: R code implementing our methods is contributed to the package aroma.affymetrix.

A human endogenous retrovirus-like (HERV) LTR formed more than 10 million years ago due to an insertion of HERV-H LTR into the 5' LTR of HERV-K is situated on human chromosomes 10, 19 and Y.

A chimeric long terminal repeat (LTR) containing the whole LTR of a human endogenous retrovirus-like element of the H family (HERV-H) inserted downstream of the core enhancer region of the 5' LTR of a HERV-K retroelement was detected and sequenced in the human 19p12 locus, known to be enriched with genes encoding zinc finger proteins. Similar chimeras were also detected in human chromosomes 10 and Y in human-hamster hybrid cells containing individual human chromosomes. This finding was interpreted as evidence of transpositions of the chimera in the genome. PCR analyses detected the chimera in the genomes of chimpanzee and gorilla, but not in that of orangutan. These data demonstrate that the chimera appeared in the primate germ cells more than 10 million years ago, before divergence of the human/chimpanzee and the gorilla lineages. The combination of the two LTRs forms a new regulatory system that can be involved in nearby gene expression.

Subfamilies and nearest-neighbour dendrogram for the LTRs of human endogenous retroviruses HERV-K mapped on human chromosome 19: physical neighbourhood does not correlate with identity level.

Sequences of 45 long terminal repeats (LTRs) of the human endogenous retroviruses HERV-K family, precisely mapped by us earlier on human chromosome 19, were determined and a nearest-neighbour dendrogram was constructed. No correlation was observed between the degree of identity of the LTR pairs and their relative positions on the chromosome. Thus, sequences of distantly located LTRs, even positioned on different chromosome arms, could be highly similar to each other, whereas those of closely located LTRs could differ significantly. We conclude that the LTRs have randomly transposed across the chromosome in the course of evolution. The alignment of the LTR sequences allowed us to assign most of the LTRs to two major subfamilies. The LTRs belonging to the first subfamily (LTR-I) are characterised by higher intrasubfamily sequence divergence than those of the second subfamily (LTR-II). The two subfamilies are easily distinguished by the presence of characteristic deletions/insertions in the LTR sequences. The higher divergence of the first subfamily members suggests that their propagation started at earlier stages of evolution, probably soon after the insertion of their ancestral sequence into the primate genome. In turn, each of the subfamilies includes several distinct branches with various degrees of intragroup divergence and with characteristic diagnostic features, suggesting that the members of the branches represent amplified copies of particular master genes which had appeared at different periods of evolution. The sequences of the LTRs demonstrate a characteristic distribution of conservative and variable regions, indicating that the LTRs might have some sequence-dependent functions in the primate genome.