Molecular Phylogenetics of the Peromyscus boylii Species Group (Rodentia: Muridae) Based on Mitochondrial Cytochrome b Sequences

A maximum likelihood method for inferring evolutionary trees from DNA sequence data was developed by Felsenstein (1981). In evaluating the extent to which the maximum likelihood tree is a significantly better representation of the true tree, it is important to estimate the variance of the difference between log likelihood of different tree topologies. Bootstrap resampling can be used for this purpose (Hasegawa et al. 1988; Hasegawa and Kishino 1989), but it imposes a great computation burden. To overcome this difficulty, we developed a new method for estimating the variance by expressing it explicitly. The method was applied to DNA sequence data from primates in order to evaluate the maximum likelihood branching order among Hominoidea. It was shown that, although the orangutan is convincingly placed as an outgroup of a human and African apes clade, the branching order among human, chimpanzee, and gorilla cannot be determined confidently from the DNA sequence data presently available when the evolutionary rate constancy is not assumed.

A thermostable DNA polymerase was used in an in vitro DNA amplification procedure, the polymerase chain reaction. The enzyme, isolated from Thermus aquaticus, greatly simplifies the procedure and, by enabling the amplification reaction to be performed at higher temperatures, significantly improves the specificity, yield, sensitivity, and length of products that can be amplified. Single-copy genomic sequences were amplified by a factor of more than 10 million with very high specificity, and DNA segments up to 2000 base pairs were readily amplified. In addition, the method was used to amplify and detect a target DNA molecule present only once in a sample of 10(5) cells.

The complete sequence of the 16,295 bp mouse L cell mitochondrial DNA genome has been determined. Genes for the 12S and 16S ribosomal RNAs; 22 tRNAs; cytochrome c oxidase subunits I, II and III; ATPase subunit 6; cytochrome b; and eight unidentified proteins have been located. The genome displays exceptional economy of organization, with tRNA genes interspersed between rRNA and protein-coding genes with zero or few noncoding nucleotides between coding sequences. Only two significant portions of the genome, the 879 nucleotide displacement-loop region containing the origin of heavy-strand replication and the 32 nucleotide origin of light-strand replication, do not encode a functional RNA species. All of the remaining nucleotide sequence serves as a defined coding function, with the exception of 32 nucleotides, of which 18 occur at the 5' ends of open reading frames. Mouse mitochondrial DNA is unique in that the translational start codon is AUN, with any of the four nucleotides in the third position, whereas the only translational stop codon is the orthodox UAA. The mouse mitochondrial DNA genome is highly homologous in overall sequence and in gene organization to human mitochondrial DNA, with the descending order of conserved regions being tRNA genes; origin of light-strand replication; rRNA genes; known protein-coding genes; unidentified protein-coding genes; displacement-loop region.