Download Detecting and Modeling Long Range Correlation in Genomic

Survey
yes no Was this document useful for you?
   Thank you for your participation!

* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project

Document related concepts

Transcriptional regulation wikipedia , lookup

DNA repair wikipedia , lookup

Promoter (genetics) wikipedia , lookup

Mutation wikipedia , lookup

Agarose gel electrophoresis wikipedia , lookup

DNA sequencing wikipedia , lookup

DNA barcoding wikipedia , lookup

Maurice Wilkins wikipedia , lookup

Whole genome sequencing wikipedia , lookup

Replisome wikipedia , lookup

Comparative genomic hybridization wikipedia , lookup

Gel electrophoresis of nucleic acids wikipedia , lookup

Bisulfite sequencing wikipedia , lookup

Nucleic acid analogue wikipedia , lookup

Point mutation wikipedia , lookup

Molecular cloning wikipedia , lookup

Community fingerprinting wikipedia , lookup

DNA supercoil wikipedia , lookup

Endogenous retrovirus wikipedia , lookup

Genome evolution wikipedia , lookup

Artificial gene synthesis wikipedia , lookup

Genomic library wikipedia , lookup

Cre-Lox recombination wikipedia , lookup

Deoxyribozyme wikipedia , lookup

Non-coding DNA wikipedia , lookup

Molecular evolution wikipedia , lookup

Transcript
Detecting and modeling long range correlation in
genomic sequences
Yi Zhou*, Archisman Rudra*, Salvatore Paxia* and Bud Mishra$*
$ CSHL
* NYU Courant Bioinformatics Group
Abstract
A genome encodes information that is needed to create complex machineries combining DNA, RNA and
proteins. However, this structure has evolved by certain basic biological processes that modify the genome in
a specific but stochastic manner, and has been shaped by selection pressure.
With complete sequences of many genomes available, it is now possible to question whether all such genome
evolution processes are adequately understood. In particular, we measure the long-range correlation (LRC) of
DNA sequences in the hope of distinguishing between different models of DNA evolution.
In order to study DNA sequence LRC, we view the DNA sequences as being generated from a random walk
model. We map a whole genomic sequence using a purine-pyrimidine binary rule. This creates a `DNA walk'
along the genome. The degree of LRC in the sequence is characterized by the Hurst exponent (H), which can
be estimated using various methods. (For infinite length: H = 0.5, no LRC; H > 0.5, positive LRC; H < 0.5,
negative LRC.)
We have analyzed various genomes using VALIS: bacteria, invertebrate and vertebrate. We observe a
consistently higher H value in the non-coding regions compared to the coding regions. Thus, the DNA walks
down the non-coding region sequences possess stronger positive LRC than those in the coding regions. In
addition, the H values in different regions increase with the evolutionary positions of the corresponding
organisms. This suggests that some cellular events tend to make DNA sequences more positively correlated
as evolution proceeds.
Based on our observations, we hypothesize that the differences in the strengths of LRC in DNA sequences
are caused by a spectrum of events affecting DNA evolution. Those include DNA polymerase stuttering,
transposons and recombination, which tend to add deletions and insertions, and natural selection and DNA
repair mechanisms, which try to eliminate the changes in the sequences. The differences in the distribution of
such spectrum in coding and non-coding regions and in different organisms cause the differences in the
degree of LRC in DNA sequences.
The hypothesis can be tested 'in silico' using Polya's Urn model. In our model, each basic DNA sequence
change is modeled using several probability distribution functions. The functions can decide the
insertion/deletion positions of the DNA fragments, the copy number of the inserted fragments and the
sequence of the inserted/deleted pieces. Moreover, those functions can be interdependent. The combination
of those basic DNA sequence changes can represent most of the natural DNA evolution events: deletion,
insertion, point mutation, tandem repeats, transposition, etc.
Our analysis and simulation were carried out on two novel computational tools: VALIS, a bioinformatics
environment for genome analysis and `Genome Grammar', a system for simulating genome evolution. Our
`Genome Grammar' can handle stochastic grammars and primitives for many kinds of mathematical
probability distributions. It allows one to apply hypothesized processes on sequences from different sources.
In particular, it enables us to conduct our experiments on DNA evolution based on Polya's Urn model.
Finally, the 'in silico' experimental results can be verified 'in vivo' using microbial mutants in the
corresponding cellular processes.