Recent developments in the statistical analysis of DNA sequences are reviewed. The pace with which sequence data are being generated and analysed has increased with the growth of the human genome project. Two areas of activity are emphasized: attention to error rates in recorded sequences, and heterogeneity in structure of sequences. There is now empirical evidence suggesting error rates in the range 0.1% ∼ 1%, and such rates will affect evolutionary studies since these are about the rates at which DNA sequences from different individuals are expected to differ. Heterogeneity for such quantities as base composition, or lengths between successive subsequences of specified types, may be sufficient to account for observed long-range correlations between bases. The need for statistical models and analyses of DNA sequence data will continue, and will offer interesting challenges.
Get full access to this article
View all access options for this article.
References
1.
Cavalli-Sforza LL, Bodmer WFThe genetics of human populations. WH Freeman:San Francisco, 1971.
2.
National Center for Biotechnology Information.Distribution CD-ROM release notes. Bethesda, MD: National Center for Biotechnology Information, 1993.
3.
US Department of Health and Human Services and US Department of Energy.Understanding our genetic inheritance. The US human genome project: the first five years FY 1991-1995. National Technical Information Service Report DOE/ER-0452P, 1990.
4.
Oliver SG et al. The complete DNA sequence of yeast chromosome III. Nature1992; 357: 38-46.
5.
Weir BSStatistical analysis of DNA sequences. Journal of the National Cancer Institute1988; 80: 395-406.
6.
US Congress, Office of Technology Assessment.Mapping our genes - genome projecis: how big, how fast? Washington, DC: US Government Printing Office, 1988: OTA-BA-373.
7.
National Research Council.Mapping and sequencing the human genome. Washington, DC: National Academy Press, 1988.
8.
Roberts L.NIH takes new tack on gene mapping. Science1992; 258: 1573.
9.
NIH/CEPH Collaborative Mapping Group.A comprehensive genetic linkage map of the human genome. Science1992; 258: 67-86.
10.
Chumakov I. et al. Continuum of overlapping clones spanning the entire human chromosome 21. Nature1992 ; 359: 380-87.
11.
Vollrath D. , Foote S., Hilton A.The human Y chromosome: a 43-interval map based on naturally occurring deletions. Science1992; 258: 52-59.
12.
Maxam A., Gilbert W.A new method for sequencing DNA. Proceedings of the National Academy of Sciences USA1977; 74: 560-64.
13.
Sanger F., Nicklen S., Coulson ARDNA sequencing with chain-terminating inhibitors. Proceedings of the National Academy of Sciences USA1977; 74: 5463-67.
14.
Endo I., Soeda E., Murakami Y., Nishi K.Human genome analysis system. Nature1991; 352: 89-90.
15.
Kieleczawa J., Dunn JJ, Studier FWDNA sequencing by primer walking with strings of contiguous hexamers. Science1992; 258: 1787-91.
16.
Edgington SMBreaking open the bottlenecks in genomic DNA sequencing. Bio/Techniques1993; 11: 39-42.
17.
Hunkapiller T., Kaiser RJ, Koop BF, Hood L.Large-scale and automated DNA sequence determination . Science1991; 254: 59-67.
18.
Little P.Complementary questions. Nature1991;352:20-21.
19.
Guyer RL, Koshland DE Jr. The molecule of the year. Science1989; 246: 1543-44.
20.
Applied Biosystems.373 DNA Sequencer. Foster City, CA: Applied Biosystems, 1992.
21.
Li H., Gyllensten UB, Ciu X., Saiki RK, Erlich HA, Arnheim N. Amplification and analysis of DNA sequences in single human sperm and diploid cells. Nature1988; 335: 414-17.
22.
Helmuth R., Fildes N., Blake E.HLA-DQα allele and genotype frequencies in various human populations, determined by using enzymatic amplification and oligonucleotide probes. American Journal of Human Genetics1990; 47: 515-23.
23.
DeSalle R., Gatesy J., Wheeler W., Grimaldi D.DNA sequences from a fossil termite in oligomiocene amber and their phylogenetic implications. Science1992; 257: 1933-36.
24.
Golenberg EM , Giannasi DE, Clegg MTChloroplast DNA sequence from a Miocene Magnolia species. Nature1990; 344: 656-58.
25.
Burks C., Cinkosky MJ, Gilna P. et al. GenBank: current status and future directions. Methods in Enzymology1991; 183: 3-22.
26.
Kahn P., Cameron G.EMBL data library. Methods in Enzymology1991 ; 183: 29-31.
27.
National Center for Biotechnology Information.NCBI News. Bethesda, MD: National Center for Biotechnology Information, 1992.
Anonymous.DNA-base reading software and DNA-analysis software. Journal of NIH Research1991; 3: 93-100.
30.
Cuticchia AJ , Arnold J., Brody H., Timberlake WECMAP: contig mapping and analysis package, a relational database for chromosome reconstruction. CABIOS1992; 8: 467-74.
31.
Shin D-G., Lee C., Zhang J., Rudd KE, Berg CMRedesigning, implementing and integrating Escherichia coli genome software tools with an object-oriented database system. CABIOS1992; 8: 227-38.
32.
Pearson PL, Maidak B., Chipperfield M., Robbins R.The human genome initiative - do databases reflect current progress? Science1991; 254: 214-15.
33.
Savakis C., Doelz R.Contamination of cDNA sequences in databases. Science1993; 259: 1677-78.
34.
von Heijne G.Sequence analysis in molecular biology: treasure trove or trivial pursuit ? New York: Academic Press, 1987.
35.
von Heijne G.Getting sense out of sequence data. Nature1988; 333: 605-607.
36.
Weir BSStatistical analysis of molecular genetic data. IMA Journal of Mathematics Applied to Medical Biology1985; 2: 1-39.
37.
Burge C., Campbell AM, Karlin S.Over- and under-representation of short oligonucleotides in DNA sequences . Proceedings of the National Academy of Sciences USA1992; 89: 1358-62.
38.
Davies K.Too many trinucleotide repeats. Nature1992 ; 358: 91.
39.
The Huntington's Disease Collaborative Research Group.A novel gene containing a trinucleotide repeat that is expanded and unstable on Huntington's disease chromosomes. Cell1993; 72: 971-83.
40.
Weir BSPopulation genetics in the forensic DNA debate. Proceedings of the National Academy of Sciences USA 1992; 89: 11654-59.
41.
Arnold J., Cuticchia AJ, Newsome DA, Jennings WW III, Ivarie R.Mono- through hexanucleotide composition of the sense strand of yeast DNA: a Markov chain analysis. Nucleic Acids Research1988; 16: 7145-58.
42.
Phillips GJ , Arnold J., Ivarie R.Mono-through hexanucleotide composition of the Escherichia coli genome: a Markov chain analysis. Nucleic Acids Research1987; 15: 2611-26.
43.
Phillips GJ , Arnold J., Ivarie R.The effect of codon usage on the oligonucleotide composition of the E. coli genome and identification of over- and underrepresented sequences by Markov chain analysis. Nucleic Acids Research1987; 15: 2627-38.
44.
Cuticchia AJ , Ivarie R., Arnold J.The application of Markov chain analysis to oligonucleotide frequency prediction and physical mapping of Drosophila melanogaster . Nucleic Acids Research1992; 20: 3651-57.
45.
Churchill GAA stochastic model for heterogeneous DNA sequences. Bulletin of MathematicalBiology1989; 51: 79-94.
46.
Tavaré S., Giddings BWSome statistical aspects of the primary structure of nucleotide sequences . In: Waterman MS ed. Mathematical methods for DNA sequences. Boca Raton: CRC Press, 1989: 117-32.
47.
Cowan R.Expected frequencies of DNA patterns using Whittle's formula. Journal of Applied Probability. 1991; 28: 886-92.
48.
Karlin S., Brendel V. Patchiness and correlations inDNA sequences. Science1993; 259: 677-80.
49.
Bernadi G.The isochore organization of the human genome. Annual Review of Genetics1989 ; 23: 637-61.
Nee S.Uncorrelated DNA walks. Nature1992; 357: 450.
55.
Chatzidimitrou-Dreismann CA, Larhammar D.Long-range correlations in DNA. Nature1993; 361: 212-13.
56.
Baer R. et al. DNA sequence and expression of the B95-8 Epstein-Barr virus genome. Nature1984; 310: 207-11.b
57.
Daniels DL, Plunkett G III, Burland V., Blattner FRAnalysis of the Escherichia coli genome: DNA sequence of the region from 84.5 to 86.5 minutes. Science1992; 257: 771-77.
58.
Posfai J., Roberts RJFinding errors in DNA sequences. Proceedings of the National Academy of Sciences USA1992; 89: 4698-702.
59.
Churchill GA , Waterman MSThe accuracy of DNA sequences: estimating sequence quality. Genomics1992; 14: 89-98.
60.
Howell N., McCullough DA, Kubacka I., Halvorson S., Mackey D.The sequence of human mtDNA: the question of errors versus polymorphisms. American Journal of Human Genetics1992; 50: 1333-36.
61.
Marzuki S., Lertrit P., Noer AS et al. Reply to Howell et al.: The need for a joint effort in the construction of a reference data base for normal sequence variants of human mtDNA. American Journal of Human Genetics1992; 50: 1337-40.
62.
Lincoln SE, Lander ESSystematic detection of errors in genetic linkage data. Genomics1992; 14: 604-10.
63.
Clark AG, Whittam TSSequencing errors and molecular evolutionary analysis. Molecular BiologyEvolution1992; 94: 744-52.
64.
Morton NE, Collins A.Standard maps of chromosome 10. Annals of Human Genetics1990; 54: 235-51.
65.
Beutow KHInfluence of aberrant observations on high-resolution linkage analysis outcomes . American Journal of Human Genetics1991; 49: 985-94.
66.
Searls DBThe linguistics of DNA. American Scientist1992; 80: 579-91.
67.
Karlin S., Macken C.Some statistical problems in the assessment of inhomogeneities of DNA sequence data. Journal of the American Statistical Association1991; 86: 27-35.
68.
Karlin S., Brendel V.Chance and statistical significance in protein and DNA sequence analysis. Science1992; 257: 39-49.
69.
Karlin S., Altschul SFMethods for assessing the statistical significance of molecular sequence features by using general scoring schemes. Proceedings of the National Academy of Sciences USA 1990; 87: 2264-68.
70.
Hodgman TCThe elucidation of protein function by sequence motif analysis. CABIOS1989; 5: 1-13.
Gonnet GH, Benner SAA word in your protein. Nature1993; 361: 121.
73.
Queen CL, Korn LJComputer analysis of nucleic acids and proteins. Methods in Enzymology1980; 65: 595-609.
74.
Needleman SB , Wunsch CDA general method applicable to the search for similarities in the amino acid sequence of two proteins. Journal of MolecularBiology1970; 48: 443-53.
75.
Smith TF, Waterman MSIdentification of common molecular subsequences. Journal of MolecularBiology1981; 147: 195-97.
76.
Waterman MS ed. Mathematical methods for DNA sequences. CRC Press: Boca Raton, 1989: 79-81.
77.
Brooks LD, Weir BS, Schaffer HEThe probabilities of similarities in DNA sequence comparisons. Genomics1988; 3: 207-16.
78.
Mott R.Maximum-likelihood estimation of the statistical distribution of Smith-Waterman local sequence similarity scores. Bulletin of MathematicalBiology1992; 54: 59-75.
79.
Gonnet GH, Cohen MA, Benner SAExhaustive matching of the entire protein sequence database. Science1992; 256: 1443-45.
80.
Barton GJComputer speed and sequence comparison. Science1992; 257: 1609.
81.
Benner SA, Cohen MA, Gonnet GHResponse to Barton. Science1992; 257: 1609-10.
82.
Taylor WRSequencing and computer time. Science1992 ; 258: 1421.
83.
Jones R.Sequence pattern matching on a massively parallel computer. CABIOS1992; 8: 377-83.
84.
Fitch WM, Smith TFOptimal sequence alignments. Proceedings of the National Academy of Sciences USA1983; 80: 1382-86.
85.
Gotoh O.Optimal sequence alignment allowing for long gaps. Bulletin of MathematicalBiology1990; 52: 359-73.
86.
Waterman MS , Eggert M., Lander E.Parametric sequence comparisons. Proceedings of the National Academy of Sciences USA1992; 89: 6090-93.
87.
Taylor WRA flexible method to align large numbers of biological sequences. Journal of MolecularEvolution1988; 2: 161-69.
88.
Carillo H., Lipman D.The multiple sequence alignment problem in biology. SIAM Journal of Applied Mathematics1988; 48: 1073-82.
89.
Altschul SF , Lipman DJTrees, stars, and multiple biological sequence alignment. SIAM Journal of Applied Mathematics1989; 49: 197-207.
90.
Altschul SFGap costs for multiple sequence alignment. Journal of TheoreticalBiology1989; 138: 297-309.
91.
Lander ESFinding similarities and differences among genomes. Nature Genetics1993; 4: 5-6.
92.
Nelson SF, McCusker JH, Sander MA, Kee Y., Modrich P., Brown POGenomic mismatch scanning: a new approach to genetic linkage mapping. Nature Genetics1993; 4: 11-17.
93.
Lisitsyn N. , Lisitsyn N., Wigler M.Cloning the differences between two complex genomes. Science1993; 259: 946-53.
94.
Miyamoto MM, Cracraft J eds. Phylogenetic analysis of DNA sequences. New York: Oxford University Press, 1991.
95.
Swofford DL, Olsen GJPhylogeny reconstruction. In: Hillis DM, Moritz C eds. Molecular systematics. Sunderland, MA: Sinauer, 1990: 411-501.
Vigilant L. , Stoneking M., Harpending H., Hawkes K., Wilson ACAfrican populations and the evolution of human mitochondrial DNA. Science1991; 253: 1503-07.
98.
Templeton ARHuman origins and analysis of mitochondrial DNA sequences. Science1992; 255: 737.
99.
Hedges SB, Kumar S., Tamura K., Stoneking M.Human origins and analysis of mitochondrial DNA sequences. Science1992; 255: 737-39.
100.
Goldman N.Maximum likelihood inference of phylogenetic trees, with special reference to a Poisson process model of DNA substitution and to parsimony analysis. Systematic Zoology1990; 39: 345-61.
101.
Weir BS, Basten CJSampling strategies for distances between DNA sequences. Biometrics1990; 46: 551-82.
102.
Nei M.Molecular evolutionary genetics. New York: Columbia University Press, 1987.
103.
Weir BSGenetic data analysis. Sunderland, MA: Sinauer, 1990.
104.
Felsenstein J.Evolutionary trees from DNA sequences: a maximum likelihood approach. Journal of MolecularEvolution1981; 17: 368-76.
105.
Muse S., Weir BSTesting for equality of evolutionary rates. Genetics1992; 132: 269-76.
106.
Wu C-I. , Li W-H.Evidence for higher rates of nucleotide substitution in rodents than in man. Proceedings of the National Academy of Sciences USA1985; 82: 1741-45.
107.
Felsenstein J.Confidence limits on phylogenies: an approach using the bootstrap. Evolution1985; 39: 783-91.
108.
Navidi WC, Churchill GA, von Haesler A.Methods for inferring phylogenies from nucleic acid sequence data by using maximum likelihood and linear invariants. Molecular Biology Evolution1991; 8: 128-43.
109.
Hudson RRGene genealogies and the coalescent process. In: Futuyma D, Antonovics J eds. Oxford surveys in evolutionary biology, Volume 7. Oxford: Oxford University Press , 1990: 1-44.
110.
Lundstrom R., Tavaré S., Ward RHEstimating substitution rates from molecular data using the coalescent. Proceedings of the National Academy of Sciences USA 1992; 89: 5961-65.
111.
Felsenstein J.Estimating effective population size from samples of sequences: inefficiency of pairwise and segregating sites as compared to phylogenetic estimates. Genetical Research1992; 59: 139-47.
112.
Felsenstein J.Estimating effective population size from samples of sequences: a bootstrap Monte Carlo integration method. Genetical Research1992; 60: 209-20.
113.
Mazzarella R., Montanaro V., Kere J. et al. Conserved sequence-tagged sites: a phylogenetic approach to genome mapping. Proceedings of the National Academy of Sciences USA1992; 89: 3681-85.
114.
Colombo P. , Yon J., Garson K., Fried M.Conservation of the organization of five tightly clustered genes over 600 million years of divergent evolution. Proceedings of the National Academy of Sciences USA1992; 89: 6358-62.
115.
Davisson MT , Lalley PA, Peters J., Doolittle DP, Hillyard AL, Searle AGReport of the comparative committee for human, mouse and other rodents. Cytogenetic Cell Genetics1991; 58: 1152-89.
116.
Green P., Lipman D., Hillier L., Waterston R., States D., Claverie J-M.Ancient conserved regions in new gene sequences and protein databases. Science1993; 259: 1711-16.
117.
Thorne JL, Kishino H., Felsenstein J.An evolutionary model for the maximum likelihood alignment of DNA sequences . Journal of Molecular Evolution1991; 33: 114-24.
118.
Thorne JL, Kishino H., Felsenstein J.Inching towards reality: an improved likelihood model of sequence evolution. Journal of MolecularEvolution1992; 34: 3-16.
119.
Ou C-Y. et al. Molecular epidemiology of HIV transmission in a dental practice. Science1992; 256: 1165-71.
120.
Smith TF, Waterman MSThe continuing case of the Florida dentist. Science1992; 256: 1155-71.
121.
Pesole G., Attimonelli M., Preparata G., Saccone C.A statistical method for detecting regions with different evolutionary dynamics in multialigned sequences. Molecular Phylogenetics andEvolution1992; 1: 91-96.
122.
Bodmer WFThe human genome sequence and the analysis of multifactorial traits. Molecular approaches to human polygenic disease. Chichester: Wiley, 1987: 215-28.
