- Methodology article
- Open Access

# A hidden Markov model approach for determining expression from genomic tiling micro arrays

- Kasper Munch†
^{1}Email author, - Paul P Gardner†
^{2}, - Peter Arctander
^{2}and - Anders Krogh
^{1}

**7**:239

https://doi.org/10.1186/1471-2105-7-239

© Munch et al; licensee BioMed Central Ltd. 2006

**Received:**11 November 2005**Accepted:**03 May 2006**Published:**03 May 2006

## Abstract

### Background

Genomic tiling micro arrays have great potential for identifying previously undiscovered coding as well as non-coding transcription. To-date, however, analyses of these data have been performed in an *ad hoc* fashion.

### Results

We present a probabilistic procedure, ExpressHMM, that adaptively models tiling data prior to predicting expression on genomic sequence. A hidden Markov model (HMM) is used to model the distributions of tiling array probe scores in expressed and non-expressed regions. The HMM is trained on sets of probes mapped to regions of annotated expression and non-expression. Subsequently, prediction of transcribed fragments is made on tiled genomic sequence. The prediction is accompanied by an expression probability curve for visual inspection of the supporting evidence. We test ExpressHMM on data from the Cheng *et al.* (2005) tiling array experiments on ten Human chromosomes [1]. Results can be downloaded and viewed from our web site [2].

### Conclusion

The value of adaptive modelling of fluorescence scores prior to categorisation into expressed and non-expressed probes is demonstrated. Our results indicate that our adaptive approach is superior to the previous analysis in terms of nucleotide sensitivity and transfrag specificity.

## Keywords

- Hide Markov Model
- Genomic Prediction
- Probe Score
- Tiling Array
- Class Hide Markov Model

## Background

Tiling micro arrays query genomic sequence in a manner not biased towards coding transcripts and has proven a valuable resource in the exploration of genomic expression. To date, the approach has been applied to Human, Arabidopsis, and Rice [3–11] and has been reviewed recently [12, 13]. The challenge of tiling micro array analysis is to categorise genomically consecutive probes as either expressed or non-expressed. The Affymetrix tiling array analyses [1, 3, 4] invoke extensive experience from gene arrays. Here, chip noise between replicates is dealt with using a quantile-normalisation procedure [14]. Background and probe sequence-specific signal is normalised using the MAS 5.0 method [15]. In order to convert expression values into transcribed fragments (transfrags) Affymetrix uses a smoothing window approach, whereby expression values within a genomic neighbourhood are averaged. A region is labelled as expressed if all probes exceed a threshold value. The threshold is determined from expressed RNA spiked into the micro array experiment. Neighbouring regions closer than a distance threshold are then joined, forming transfrag predictions. We will use the collective acronym ANTM for for Affymetrix's normalisation and transfrag methods. For a review of existing algorithms for tiling micro array analysis see Royce *et al.* [16].

Previous related work includes analysis of ChlP-chip experiments using Affymetrix tiling arrays to predict transcription factor binding sites [17]. These ChlP-enriched regions on the tiling arrays are predicted using log-odds scores for each probe generated by a two-state hidden Markov model. Tiling array intensities have been used together with a rudimentary HMM gene finder to establish likely exon-intron transcripts in windows of significant expression intensity [18]. This work, however, does not treat the probe intensities as observables of the HMM but as separate evidence supporting prediction of expression.

The approach to tiling array analysis presented here employs a hidden Markov model (HMM), trained directly on the correspondence between tiling array fluorescence scores and annotation of expression. The model is then used to categorise consecutive probes as expressed or non-expressed. Training the model using annotation of expression presents a problem since, in many cases, annotation does not represent actual expression. This may be due to a number of factors, such as non-expressed genes, alternative splicing, un-annotated exons, antisense transcription, and non-coding RNAs. This type of incomplete information presents a problem to machine learning approaches as missing information compromises the training of the model unless properly addressed. We describe a two-step modelling and training approach that allows the training algorithm to down-weight dubious annotation in the context of the relevant cell line. The adaptive nature of our method avoids the *ad hoc* thresholds characteristic of previous work. As a case study we analyse Human tiling array data of a single cell-line from Cheng *et al.* [1].

## Methods

### Data

The data we use is derived from the tiling array experiments on Human chromosomes 6, 7, 13, 14, 19, 20, 21, 22, X, and Y [1]. The tiling arrays use 25-mer probes querying genomic positions at five nucleotide intervals on the genomic sequence. Only non-repetitive sequence is tiled resulting in a coverage of approximately 30%. From the Cheng *et al.* raw array data we analyse the subset querying poly-adenylated RNA extracted from Human neuroblastoma cells (SK-N-AS cell line).

### Array and probe normalisation

The Affymetrix MAS 5.0 method addresses probe and array normalisation by, in addition to perfect match (PM) probes, also designing mismatch (MM) probes, where the complementary base at the 13th position of the 25 mer probe is used instead of the actual base. The mismatch probe is expected to capture much of the probe sequence and array-specific contributions to the signal as well as problematic probe eccentricities such as secondary-structure and cross-hybridisation effects while capturing little real signal. The difference in fluorescence score for the PM probe and an Idealised-MM probe score is then used to estimate expression signal from the queried genomic position.

We tested a number of strategies for array and probe normalisation. These included *PM* - *MM* with and without quantile normalisation and a novel approach inspired by Naef and Magnasco [19] that shares similarities with GC-RMA [20, 21].

For the normalisation inspired by Naef and Magnasco sequence and array specific signal is modelled in a way that allows this to be subtracted from the raw signal. For each chip a large sample of probe scores are selected and used to fit position-specific estimates of nucleotide contribution to the signal. Three free parameters for each position times the 25 parameters for each probe sequence position gives rise to 75 parameters. These are modelled by the following linear system:

where each nucleotide type (*A, C, G* or *T*) in position *k* in a 25-mer probe is modelled by a triple in the indicator function *a*_{
i
}(*a*_{3k-2}, *a*_{3k-1}, *a*_{3k}). I.e. an *A* at position *k* has corresponding coefficients (*a*_{3k-2}, *a*_{3k-1}, *a*_{3k}) = (1, 0, 0), a *C* has corresponding coefficients (0, 1, 0), a *G* has corresponding coefficients (0, 0, 1), and a *T* has corresponding coefficients (-1, -1, -1). A least-squares solution to this system of equations for *x*_{
i
}provides estimates of the expected contribution to the expression for any nucleotide at position *k*. That is, the estimated contribution to the expression level of an *A* at position *k* is given by *x*_{3k-2}, a *C* at position *k* is given by *x*_{3k-1}, a *G* at position *k* is given by *x*_{3k}, and a *T* at position *k* is given by -*x*_{3k-2}- *x*_{3k-1}- *x*_{3k}. To obtain a representative expression level over different replicates the mean log-likelihoods are computed:

*PM*-

*MM*approach rather than using an idealised

*MM*we use the raw values of both

*PM*and

*MM*. A representative expression level over different replicates using the

*MM*probe are computed as:

The performance of each normalisation strategy was evaluated (see discussion) and $\overline{D}$ was chosen as the more appropriate.

### Training data

### HMM Architecture and training

The correspondence between scores and annotated expression, presented by the training set, is captured by an HMM, see e.g. Durbin *et al.* [24]. An HMM is a probabilistic model that consists of a set of connected hidden states that each emit observables. In this case, states represent expression and non-expression and these states emit scores that constitute the sequence of probe scores over expressed or non-expressed regions. First order emission probabilities are used to model the dependency of scores for neighbouring-probes.

*E*

_{0}and

*N*

_{0}) model the bulk of probe scores in expressed and non-expressed regions, capturing the score distributions of the two cases. The loop states are connected by states that capture score gradients induced by genomically consecutive probes that overlap borders of expressed regions. The loop probabilities are tied during training of the HMM and thus estimated as one parameter. This is done in order to avoid bias of predictions towards typical exon length. The HMM used is discrete and expression scores are binned. Considering the distributions of scores for probes annotated as expressed and non-expressed, we choose 13 bins uniformly covering the scores from -50 to 50. Another two bins capture extreme positive and negative values, (see Figure 4).

The type of HMM used is a Class HMM [25]. The collection of states in a Class HMM is divided into classes, each designated by a label. The observables (i.e. probe scores) in the training set are each assigned a label corresponding to a class. An observable with a given label can only be emitted by a state belonging to the class designated by that label. Hence, the labelling of states and observables in the training set determines which parts of the training data are used to estimate which parts of the model. In our case, the labels are expression and non-expression and correspond to the red and blue colours in Figure 3. The advantage of the Class HMM is that it allows all model parameters to be estimated simultaneously from sequences of labelled probe scores. Alternatively, the individual parts would have to be estimated separately from expressed and non-expressed probe scores. The model is trained using the Baum-Welch algorithm [26].

### Self-supervised re-training

*X*

_{ i }, in the core model but does not contribute to the estimation of these parameters. The labelling of states in the shadow model differs from those of the core model in that the states ${{E}^{\prime}}_{i}$ model non-expressed probes labelled as expressed and the ${{N}^{\prime}}_{0}$ state models expressed probes labelled as non-expressed. Additional states impose a minimum on the number of consecutive probe scores that can be captured by the shadow model. Without this constraint the model is not able to distinguish natural variation and noise from dubious annotation. The shadow component models regions annotated as expressed using the score distribution initially learned from non-expressed regions. Analogously regions annotated as non-expressed are modelled using the score distribution initially learned from expressed regions. Re-training after adding the shadow model allows the Baum-Welch algorithm to weight the contribution of training information to parameter estimation in proportion to how characteristic this is: In the E-step of the Baum-Welch algorithm the posterior probabilities of each emission from state

*X*

_{ i }are summed over all observables in the training set. The ${{X}^{\prime}}_{i}$ state does not contribute to this sum. Hence, the posterior probability of using state

*X*

_{ i }relative to ${{X}^{\prime}}_{i}$ for a given observable effectively constitutes a weight reflecting how likely the annotation is at the given position in the training set given the probe score distributions learned in the initial training. As a result, each subsequence of labelled probe scores in the training set contributes to estimation of emission probabilities to an extent determined by its reliability. It should be noted that the maximum likelihood of the model is obtained when only transition probabilities in the slave model are positive – the situation where the model discards all training material. The pre-training, however, produces a starting point that allows the model to reach a desirable local optimum.

Once the parameters of the HMM have been estimated the parallel shadow architecture is removed. The sequences of probe scores corresponding to the tiling of each chromosome are then decoded using an N-best algorithm [25]. In this step the most likely sequences of contiguous expressed and non-expressed probes are established. In addition, the forward-backward algorithm [27] is used to assign a posterior probability of expression to each probe. The prediction for each sequence of probe scores is then mapped onto the nucleotide sequence of the corresponding chromosome.

The tiling of the chromosomes only covers non-repeat sequence. Some predictions may span untiled regions. To avoid this, gaps in the tiling larger than ten bases are removed from the predictions.

### Evaluation of performance

To establish the relative performance of ExpressHMM and Affymetrix's normalisation and transfrag method (ANTM) we test each method on the training set of probes covering RefSeq genes. To ensure that test and training material do not overlap, the performance of ExpressHMM is evaluated using six-fold cross-validation. The ANTM predictions are mapped onto the training set for comparison.

To acknowledge the most recent repeat predictions we remove, from all genomic predictions, all overlaps to the Human Repeat Masker annotation from the UCSC genome browser database (October 2005) [28]. Performance statistics are calculated using the Eval program [29]. Sensitivity is defined as true positives divided by the sum of true positives and false negatives but is calculated as the number of annotation objects (e.g. nucleotides or sequence blocks) overlapped by predictions, divided by the total number of annotation objects. Specificity is defined as true positives divided by the sum of true positives and false positives but is calculated as the number of prediction objects overlapping annotation objects divided by the total number of prediction objects. The statistics are calculated this way to ensure a sensible treatment of one-to-many overlaps. Since the motivation in tiling array analysis is to find un-annotated expression a high specificity based on known expression is not a goal in itself. For this reason we report 1 – specificity as "Novel Prediction".

To further evaluate the performance we compare the two sets of predictions to annotation of Human ESTs, mRNAs, and known genes (UCSC Oct. 2005). Only annotation of tiled genomic sequence is included. Before evaluation the expression annotation is pooled and all overlapping annotation is collapsed into a maximal set where each nucleotide is only represented once.

Within both sets of annotation only a subset are actually expressed in the SK-N-AS cell line, and a substantial fraction of expression remains un-annotated. Collapsing all overlapping annotation presents a similar problem especially affecting splice site statistics. Relying on this annotation will make estimates of sensitivity artificially low and the absolute values of this statistic should not be given too much weight. The relative values for ExpressHMM and ANTM, however, should accurately reflect relative performances.

## Results

Performance of ExpressHMM and ANTM on the training set of probes covering RefSeq genes. The statistics show to what extent the two approaches make predictions in accordance with the RefSeq annotation. The exon/transfrag statistics refer to the sets of contiguous probes representing exons/transfrags. The splice site sensitivities refer to the single probes constituting the borders of such sets.

ExpressHMM | ANTM | |
---|---|---|

Probe sensitivity | 65% | 44% |

Novel probes predicted | 65% | 49% |

Overlap exon Sensitivity | 58% | 58% |

Novel transfrags predicted | 27% | 58% |

5' Splice site sensitivity | 16% | 9% |

3' Splice site sensitivity | 16% | 13% |

Performance of ExpressHMM and ANTM on the tiled regions of the ten chromosomes. The statistics show to what extent the two approaches make predictions in accordance with the collection of Human EST, mRNA, and known gene annotation.

ExpressHMM | ANTM | |
---|---|---|

Nucleotide sensitivity | 35% | 20% |

Novel nucleotides | 55% | 46% |

Nucleotide Count | 43523440 | 21190998 |

Transfrag sensitivity | 27% | 22% |

Novel transfrags | 44% | 58% |

Transfrag Count | 88604 | 170788 |

Average Length | 491.21 | 124.08 |

Median Length | 325 | 79 |

In addition to the most likely categorisation of probes ExpressHMM also calculates a posterior probability of expression for each probe. When mapped to the nucleotide sequence this results in an expression probability curve parallel to the predictions allowing the user to evaluate the significance of predictions and to directly view the evidence of expression.

All predictions can be downloaded and viewed from our web site [2].

## Discussion

### Evaluation of normalisation strategies

In the process of developing our method we compared the performance of the *LL* and $\overline{D}$ normalisation strategies on the training set. Six-fold cross-validation was used to ensure no over-training. The nucleotide-level sensitivities of *LL* and $\overline{D}$ was approximately 38.8% and 44.6% respectively with corresponding specificities of 51.1% and 20.7%. Subsequent optimisation of the model resulted in further improvement of these values (see Table 2). The *LL* based normalisation is significantly more specific than the $\overline{D}$ approach but at a significant cost to sensitivity. We have used the $\overline{D}$ approach in this study because we view sensitivity as a more desirable trait than specificity for our predictor.

*LL*distribution for probes tiling coding regions. For the $\overline{D}$ method both the left and right tails of the distribution are heavy for probes tiling coding regions. The HMM only requires significantly different distributions between expressed and non-expressed regions in order to perform well. The two heavy tails generated by the $\overline{D}$ method seems to supply more discriminative information to the HMM than the single heavy tail generated by the

*LL*method. The effect of quantile normalisation was tested. This did not improve performance, and is not used for the results presented.

### Relative performances on the training set

As shown in Table 1, our probe sensitivity is higher than ANTM. Of the total number of probes annotated as expressed we predict 21% more than ANTM. Note that the upper limit to sensitivity is not 100% as only a fraction of the annotated sequence is actually expressed. For identification of exons ExpressHMM is equal to ANTM but better at predicting the borders of expressed regions correctly. The percentage of predicted probes that do not overlap RefSeq exons is 16% larger than ANTM. In contrast, the number of predicted transfrags not overlapping a RefSeq exon comprise only 27% of the ExpressHMM predictions whereas this number is 58% for ANTM. This means that even with the higher probe sensitivity ExpressHMM is more conservative than ANTM in predicting novel transfrags.

The HMM in ExpressHMM models the signal variability within expressed and non-expressed regions. As a result the algorithm does not predict very short expressed or non-expressed regions unless the relevant probe scores are highly characteristic. ANTM's threshold approach allows for little variability in score along a region. For these reasons ANTM tends to split exons into smaller predictions whilst ExpressHMM tends to join proximal exons into a single prediction. These differences are reflected in the length statistics given in Table 2. In effect, the two methods complement each other.

### Relative performances on the chromosome level

The comparison of genomic predictions to all Human EST, mRNA and known gene annotation within the tiled regions also reflects the good performance of ExpressHMM. The nucleotide sensitivity is 75% larger than ANTM's and our transfrag sensitivity is 23% larger. ExpressHMM identifies an 19% larger fraction of novel nucleotides. In contrast, however, the fraction of novel transfrags among ExpressHMM predictions is 24% smaller than among ANTM predictions.

ExpressHMM predicts roughly twice as many expressed nucleotides as ANTM. Still the total ExpressHMM predictions only amount to 1.5% of the genome. For comparison, recent cDNA based findings in Mouse indicates that 62% percent of the genome is transcribed and that the number of transcripts is at least an order of magnitude larger than the estimated number of genes [31].

As indicated by the nucleotide plot in Figure 7 the lower nucleotide specificity of ExpressHMM is not due to false positives uniformly distributed across the chromosomes. Rather, the distribution closely follows the density of both ANTM predictions and RefSeq exons. Note the excess of prediction in the sub-telomeric regions. This is most likely due to frequent duplication events in these regions as observed in the recent Chimpanzee-Human genome comparison [32]. These duplications are expected to cause extensive cross-hybridisation. The two large peaks around 4e+07 correspond to similar-sized peaks in the density of EST and mRNA annotation (data not shown).

The transfrag plot in Figure 7 shows the close correspondence between the distributions of ExpressHMM predictions and occurrence of RefSeq exons.

### Hidden Markov model

Individual values of probe scores form only a fraction of the signal characteristic of expressed and non-expressed regions. The distribution of scores within and between regions of expression constitutes an important source of information that is not utilised by simple threshold approaches. To capture as much relevant information as possible we use an adaptive approach. This learns directly from probe scores covering regions where expression is annotated. As a result the parameters of the model are optimal for the training data. These properties takes the place of the *ad hoc* cut-offs applied in previous approaches.

Before settling on the HMM described above, a variety of models were tested. These included different order Markov models with and without modelling of transfrag borders. None of these offered apparent advantages over the one presented here. The order of the HMM has a pronounced effect on the results. For zeroth and first order emissions the nucleotide statistics are very similar. On the transfrag level, however, the zeroth order model is more sensitive (72% identified RefSeq exons) but less specific (55% novel transfrags predicted). We also experimented with excluding score sequences from the training set that did not seem unambiguously expressed. This improved the performance of some models but not that presented here.

The training of the model parameters does not include any nucleotide sequence information. In addition, to avoid any discriminative length modelling characteristic of coding exons the trained intron and exon length distributions are geometric and estimated as one to make them identical. As a result, the training data does not bias ExpressHMM predictions towards the length distribution of coding exons.

The expression probability curve that accompanies the ExpressHMM prediction has a direct interpretation as the evidence of expression given the model. It is a valuable tool to visually assess the significance of each transfrag as well as to identify regions with evidence not sufficient to result in a transfrag. With slight modifications ExpressHMM can equally well be used for the analysis of tiling micro arrays with non-overlapping probes of various length.

## Conclusion

In addition to performing better than the approach presented by Cheng *et al.* [1], the adaptive approach used by ExpressHMM avoids *ad hoc* thresholds for the analysis of signal data. Decisions regarding prediction are learned from the data in an automated fashion. In addition to predicting the most likely categorisation into expression and non-expression each genomic position queried by a probe is assigned a probability of expression. The resulting graph has a clearer interpretation than a smoothed fluorescence score. We expect that tiling arrays will play an increasing role in the investigation of genomic output.

## Notes

## Declarations

### Acknowledgements

We thank Stefan Bekiranov, Tom Gingeras, Nils Hansen, and Alex Sczyrba for help and valuable advice.

## Authors’ Affiliations

## References

- Cheng J, Kapranov P, Drenkow J, Dike S, Brubaker S, Patel S, Long J, Stern D, Tammana H, Helt G, Sementchenko V, Piccolboni A, Bekiranov S, Bailey D, Ganesh M, Ghosh S, Bell I, Gerhard D, Gingeras T:
**Transcriptional maps of 10 human chromosomes at 5-nucleotide resolution.***Science*2005,**308**(5725):1149–1154. 10.1126/science.1108625View ArticlePubMedGoogle Scholar **ExpressHMM web site**[http://www.binf.ku.dk/~kasper/expresshmm]- Kapranov P, Cawley S, Drenkow J, Bekiranov S, Strausberg R, Fodor S, Gingeras T:
**Large-scale transcriptional activity in chromosomes 21 and 22.***Science*2002,**296**(5569):916–919. 10.1126/science.1068597View ArticlePubMedGoogle Scholar - Kampa D, Cheng J, Kapranov P, Yamanaka M, Brubaker S, Cawley S, Drenkow J, Piccolboni A, Bekiranov S, Helt G, Tammana H, Gingeras T:
**Novel RNAs identified from an in-depth analysis of the transcriptome of human chromosomes 21 and 22.***Genome Res*2004,**14**(3):331–342. 10.1101/gr.2094104PubMed CentralView ArticlePubMedGoogle Scholar - Selinger D, Cheung K, Mei R, Johansson E, Richmond C, Blattner F, Lockhart D, Church G:
**RNA expression analysis using a 30 base pair resolution Escherichia coli genome array.***Nat Biotechnol*2000,**18**(12):1262–1268. 10.1038/82367View ArticlePubMedGoogle Scholar - Shoemaker D, Schadt E, Armour C, He Y, Garrett-Engele P, McDonagh P, Loerch P, Leonardson A, Lum P, Cavet G, Wu L, Altschuler S, Edwards S, King J, Tsang J, Schimmack G, Schelter J, Koch J, Ziman M, Marton M, Li B, Cundiff P, Ward T, Castle J, Krolewski M, Meyer M, Mao M, Burchard J, Kidd M, Dai H, Phillips J, Linsley P, Stoughton R, Scherer S, Boguski M:
**Experimental annotation of the human genome using microarray technology.***Nature*2001,**409**(6822):922–927. 10.1038/35057141View ArticlePubMedGoogle Scholar - Rinn J, Euskirchen G, Bertone P, Martone R, Luscombe N, Hartman S, Harrison P, Nelson F, Miller P, Gerstein M, Weissman S, Snyder M:
**The transcriptional activity of human Chromosome 22.***Genes Dev*2003,**17**(4):529–540. 10.1101/gad.1055203PubMed CentralView ArticlePubMedGoogle Scholar - Yamada K, Lim J, Dale J,
*et al*.:**Empirical analysis of transcriptional activity in the Arabidopsis genome.***Science*2003,**302**(5646):842–846. 10.1126/science.1088305View ArticlePubMedGoogle Scholar - Stolc V, Gauhar Z, Mason C, Halasz G, van Batenburg M, Rifkin S, Hua S, Herreman T, Tongprasit W, Barbano P, Bussemaker H, White K:
**A gene expression map for the euchromatic genome of Drosophila melanogaster.***Science*2004,**306**(5696):655–660. 10.1126/science.1101312View ArticlePubMedGoogle Scholar - Ishkanian A, Malloff C, Watson S, DeLeeuw R, Chi B, Coe B, Snijders A, Albertson D, Pinkel D, Marra M, Ling V, MacAulay C, Lam W:
**A tiling resolution DNA microarray with complete coverage of the human genome.***Nat Genet*2004,**36**(3):299–303. 10.1038/ng1307View ArticlePubMedGoogle Scholar - Stolc V, Samanta M, Tongprasit W, Sethi H, Liang S, Nelson D, Hegeman A, Nelson C, Rancour D, Bednarek S, Ulrich E, Zhao Q, Wrobel R, Newman C, Fox B, Phillips G, Markley J, Sussman M:
**Identification of transcribed sequences in Arabidopsis thaliana by using high-resolution genome tiling arrays.***Proc Natl Acad Sci U S A*2005,**102**(12):4453–4458. 10.1073/pnas.0408203102PubMed CentralView ArticlePubMedGoogle Scholar - Johnson J, Edwards S, Shoemaker D, Schadt E:
**Dark matter in the genome: evidence of widespread transcription detected by microarray tiling experiments.***Trends Genet*2005,**21**(2):93–102. 10.1016/j.tig.2004.12.009View ArticlePubMedGoogle Scholar - Mockler TC, Chan S, Sundaresan A, Chen H, Jacobsen SE, Ecker JR:
**Applications of DNA tiling arrays for whole-genome analysis.***Genomics*2005,**85:**1–15. 10.1016/j.ygeno.2004.10.005View ArticlePubMedGoogle Scholar - Bolstad B, Irizarry R, Astrand M, Speed T:
**A comparison of normalization methods for high density oligonucleotide array data based on variance and bias.***Bioinformatics*2003,**19**(2):185–193. 10.1093/bioinformatics/19.2.185View ArticlePubMedGoogle Scholar - Affymetrix:
**Statistical algorithms description document.***Tech rep Affymetrix*2002. [http://www.affymetrix.com/support/technical/whitepapers/sadd_whitepaper.pdf]Google Scholar - Royce TE, Rozowsky JS, Bertone P, Samanta M, Stolc V, Weissman S, Snyder M, Gerstein M:
**Issues in the analysis of oligonucleotide tiling microarrays for transcript mapping.***Trends Genet*2005,**21**(8):466–475. 10.1016/j.tig.2005.06.007PubMed CentralView ArticlePubMedGoogle Scholar - Li W, Meyer CA, Liu XS:
**A hidden Markov model for analyzing ChlP-chip experiments on genome tiling arrays and its application to p53 binding sequences.***Bioinformatics*2005,**21**(Suppl 1):i274-i282. 10.1093/bioinformatics/bti1046View ArticlePubMedGoogle Scholar - Toyoda T, Shinozaki K:
**Tiling array-driven elucidation of transcriptional structures based on maximum-likelihood and Markov models.***Plant J*2005,**43**(4):611–621. 10.1111/j.1365-313X.2005.02470.xView ArticlePubMedGoogle Scholar - Naef F, Magnasco M:
**Solving the riddle of the bright mismatches: labeling and effective binding in oligonucleotide arrays.***Phys Rev E Stat Nonlin Soft Matter Phys*2003,**68**(1 Pt 1):011906–011906.View ArticlePubMedGoogle Scholar - Wu Z, RA I, Gentleman R, Martinez-Murillo F, Spencer F:
**A model-based background adjustment for oligonucleotide expression arrays.***Journal of the American Statistical Association*2004,**99**(468):909–917. 10.1198/016214504000000683View ArticleGoogle Scholar - Wu Z, Irizarry R:
**Stochastic models inspired by hybridization theory for short oligonucleotide arrays.***J Comput Biol*2005,**12**(6):882–893. 10.1089/cmb.2005.12.882View ArticlePubMedGoogle Scholar - Hubbard T, Andrews D, Caccamo M, Cameron G, Chen Y, Clamp M, Clarke L, Coates G, Cox T, Cunningham F, Curwen V, Cutts T, Down T, Durbin R, Fernandez-Suarez XM, Gilbert J, Hammond M, Herrero J, Hotz H, Howe K, Iyer V, Jekosch K, Kahari A, Kasprzyk A, Keefe D, Keenan S, Kokocinsci F, London D, Longden I, McVicker G, Melsopp C, Meidl P, Potter S, Proctor G, Rae M, Rios D, Schuster M, Searle S, Severin J, Slater G, Smedley D, Smith J, Spooner W, Stabenau A, Stalker J, Storey R, Trevanion S, Ureta-Vidal A, Vogel J, White S, Woodwark C, Birney E:
**Ensembl 2005.***Nucleic Acids Res*2005,**33**(Database issue):D447-D453. 10.1093/nar/gki138PubMed CentralView ArticlePubMedGoogle Scholar - Pruitt KD, Katz KS, Sicotte H, Maglott DR:
**Introducing RefSeq and LocusLink: curated human genome resources at the NCBI.***Trends Genet*2000,**16:**44–47. 10.1016/S0168-9525(99)01882-XView ArticlePubMedGoogle Scholar - Durbin R, Eddy S, Krogh A, Mitchison G:
*Biological Sequence Analysis: Probabilistic models of protein and nucleic acids*. Cambridge University Press; 1998.View ArticleGoogle Scholar - Krogh A:
**Two methods for improving performance of an HMM and their application for gene finding.***Proc Int Conf Intell Syst Mol Biol*1997,**5:**179–186.PubMedGoogle Scholar - Baum LE:
**An equality and associated maximization technique in statistical estimation for probabilistic functions of Markov processes.***Inequalities*1972,**3:**1–8.Google Scholar - Rabiner LR:
**A tutorial on hidden Markov models and selected applications in speech recognition.***Proc IEEE*1989,**77**(2):257–286. 10.1109/5.18626View ArticleGoogle Scholar - Karolchik D, Baertsch R, Diekhans M, Furey T, Hinrichs A, Lu Y, Roskin K, Schwartz M, Sugnet C, Thomas D, Weber R, Haussler D, Kent W:
**The UCSC Genome Browser Database.***Nucleic Acids Res*2003,**31:**51–54. 10.1093/nar/gkg129PubMed CentralView ArticlePubMedGoogle Scholar - Keibler E, Brent MR:
**Eval: a software package for analysis of genome annotations.***BMC Bioinformatics*2003,**4:**50. 10.1186/1471-2105-4-50PubMed CentralView ArticlePubMedGoogle Scholar - Siepel A, Bejerano G, Pedersen JS, Hinrichs AS, Hou M, Rosenbloom K, Clawson H, Spieth J, Hillier LW, Richards S, Weinstock GM, Wilson RK, Gibbs RA, Kent WJ, Miller W, Haussler D:
**Evolutionarily conserved elements in vertebrate, insect, worm, and yeast genomes.***Genome Res*2005,**15**(8):1034–1050. 10.1101/gr.3715005PubMed CentralView ArticlePubMedGoogle Scholar - Carninci P, Kasukawa T, Katayama S, Gough J, Frith MC, Maeda N, Oyama R, Ravasi T, Lenhard B, Wells C, Kodzius R, Shimokawa K, Bajic VB, Brenner SE, Batalov S, Forrest AR, Zavolan M, Davis MJ, Wilming LG, Aidinis V, Allen JE, Ambesi-Impiombato A, Apweiler R, Aturaliya RN, Bailey TL, Bansal M, Baxter L, Beisel KW, Bersano T, Bono H, Chalk AM, Chiu KP, Choudhary V, Christoffels A, Clutterbuck DR, Crowe ML, Dalla E, Dalrymple BP, de Bono B, Gatta GD, di Bernardo D, Down T, Engstrom P, Fagiolini M, Faulkner G, Fletcher CF, Fukushima T, Furuno M, Futaki S, Gariboldi M, Georgii-Hemming P, Gingeras TR, Gojobori T, Green RE, Gustincich S, Harbers M, Hayashi Y, Hensch TK, Hirokawa N, Hill D, Huminiecki L, Iacono M, Ikeo K, Iwama A, Ishikawa T, Jakt M, Kanapin A, Katoh M, Kawasawa Y, Kelso J, Kitamura H, Kitano H, Kollias G, Krishnan SP, Kruger A, Kummerfeld SK, Kurochkin IV, Lareau LF, Lazarevic D, Lipovich L, Liu J, Liuni S, McWilliam S, Babu MM, Madera M, Marchionni L, Matsuda H, Matsuzawa S, Miki H, Mignone F, Miyake S, Morris K, Mottagui-Tabar S, Mulder N, Nakano N, Nakauchi H, Ng P, Nilsson R, Nishiguchi S, Nishikawa S, Nori F, Ohara O, Okazaki Y, Orlando V, Pang KC, Pavan WJ, Pavesi G, Pesole G, Petrovsky N, Piazza S, Reed J, Reid JF, Ring BZ, Ringwald M, Rost B, Ruan Y, Salzberg SL, Sandelin A, Schneider C, Schonbach C, Sekiguchi K, Semple CA, Seno S, Sessa L, Sheng Y, Shibata Y, Shimada H, Shimada K, Silva D, Sinclair B, Sperling S, Stupka E, Sugiura K, Sultana R, Takenaka Y, Taki K, Tammoja K, Tan SL, Tang S, Taylor MS, Tegner J, Teichmann SA, Ueda HR, van Nimwegen E, Verardo R, Wei CL, Yagi K, Yamanishi H, Zabarovsky E, Zhu S, Zimmer A, Hide W, Bult C, Grimmond SM, Teasdale RD, Liu ET, Brusic V, Quackenbush J, Wahlestedt C, Mattick JS, Hume DA, Kai C, Sasaki D, Tomaru Y, Fukuda S, Kanamori-Katayama M, Suzuki M, Aoki J, Arakawa T, Iida J, Imamura K, Itoh M, Kato T, Kawaji H, Kawagashira N, Kawashima T, Kojima M, Kondo S, Konno H, Nakano K, Ninomiya N, Nishio T, Okada M, Plessy C, Shibata K, Shiraki T, Suzuki S, Tagami M, Waki K, Watahiki A, Okamura-Oho Y, Suzuki H, Kawai J, Hayashizaki Y:
**The transcriptional landscape of the mammalian genome.***Science*2005,**309**(5740):1559–1563. 10.1126/science.1112014View ArticlePubMedGoogle Scholar - Linardopoulou EV, Williams EM, Fan Y, Friedman C, Young JM, Trask BJ:
**Human subtelomeres are hot spots of interchromosomal recombination and segmental duplication.***Nature*2005,**437**(7055):94–100. 10.1038/nature04029PubMed CentralView ArticlePubMedGoogle Scholar - Griffiths-Jones S:
**The microRNA Registry.***Nucleic Acids Res*2004,**32**(Database issue):D109-D111. 10.1093/nar/gkh023PubMed CentralView ArticlePubMedGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.