The number of nucleotides required to determine the branching order of three species,with special reference to the human-chimpanzee-gorilla divergence

Summary A mathematical theory for computing the probabilities of various nucleotide configurations among related species is developed, and the probability of obtaining the correct tree (topology) from nucleotide sequence data is evaluated using models of evolutionary trees that are close to the tree of mitochondrial DNAs from human, chimpanzee, gorilla, orangutan, and gibbon. Special attention is given to the number of nucleotides required to resolve the branching order among the three most closely related organisms (human, chimpanzee, and gorilla). If the extent of DNA divergence is close to that obtained by Brown et al. for mitochondrial DNA and if sequence data are available only for the three most closely related organisms, the number of nucleotides (m^*) required to obtain the correct tree with a probability of 95% is about 4700. If sequence data for two outgroup species (orangutan and gibbon) are available, m^* becomes about 2600–2700 when the transformed distance, distance-Wagner, maximum parsimony, or compatibility method is used. In the unweighted pair-group method, m^* is not affected by the availability of data from outgroup species. When these five different tree-making methods, as well as Fitch and Margoliash's method, are applied to the mitochondrial DNA data (1834 bp) obtained by Brown et al. and by Hixson and Brown, they all give the same phylogenetic tree, in which human and chimpanzee are most closely related. However, the trees considered here are gene trees, and to obtain the correct species tree, sequence data for several independent loci must be used.     

Mitochondrial DNA sequences of primates: Tempo and mode of evolution

Summary We cloned and sequenced a segment of mitochondrial DNA from human, chimpanzee, gorilla, orangutan, and gibbon. This segment is 896 bp in length, contains the genes for three transfer RNAs and parts of two proteins, and is homologous in all 5 primates. The 5 sequences differ from one another by base substitutions at 283 positions and by a deletion of one base pair. The sequence differences range from 9 to 19% among species, in agreement with estimates from cleavage map comparisons, thus confirming that the rate of mtDNA evolution in primates is 5 to 10 times higher than in nuclear DNA. The most striking new finding to emerge from these comparisons is that transitions greatly outnumber transversions. Ninety-two percent of the differences among the most closely related species (human, chimpanzee, and gorilla) are transitions. For pairs of species with longer divergence times, the observed percentage of transitions falls until, in the case of comparisons between primates and non-primates, it reaches a value of 45. The time dependence is probably due to obliteration of the record of transitions by multiple substitutions at the same nucleotide site. This finding illustrates the importance of choosing closely related species for analysis of the evolutionary process. The remarkable bias toward transitions in mtDNA evolution necessitates the revision of equations that correct for multiple substitutions at the same site. With revised equations, we calculated the incidence of silent and replacement substitutions in the two protein-coding genes. The silent substitution rate is 4 to 6 times higher than the replacement rate, indicating strong functional constraints at replacement sites. Moreover, the silent rate for these two genes is about 10% per million years, a value 10 times higher than the silent rate for the nuclear genes studied so far. In addition, the mean substitution rate in the three mitochondrial tRNA genes is at least 100 times higher than in nuclear tRNA genes. Finally, genealogical analysis of the sequence differences supports the view that the human lineage branched off only slightly before the gorilla and chimpanzee lineages diverged and strengthens the hypothesis that humans are more related to gorillas and chimpanzees than is the orangutan.Abbreviations mtDNA mitochondrial DNA - bp base pair - URF unidentified reading frame     

Neighboring base composition is strongly correlated with base substitution bias in a region of the chloroplast genome

Nucleotide sequence from a region of the chloroplast genome is presented for 12 species spanning four subfamilies of the grass family. The region contains the coding sequence for the rbcL gene and the intergenic spacer between the gene coding the large subunit of ribulose-1,5-bisphosphate carboxylase/oxygenase rbcL and the photosystem I gene psal. This intergenic spacer contains a pseudogene for rpl23 as well as two noncoding segments with different A+T contents. Using the sequence of rbcL a chloroplast phylogeny of this family was constructed by parsimony. Variable sites of the two noncoding segments were traced onto the phylogeny to study the dynamics of base substitution. This was also performed for the fourfold-degenerate sites of the rbcL gene. A wide variation in transversion/transition is observed between the two noncoding segments and between the noncoding DNA and the fourfold-degenerate sites of rbcL This variation is correlated with regional A+T content. As regional A+T content decreases, the ratio of transversions to transitions also decreases. Substitutions were then scored in relation to neighboring base composition. The composition of the two bases immediately flanking each substitution is highly correlated with the transversion/transition bias. When both the 5 and 3 flanking bases are an A or a T, transversions are observed 2.2 times as frequently as transitions. When either or both neighbors are a C or a G, the opposite trend is found; transitions are observed 1.5 times more frequently than transversions. Correspondence to: Brian R. Morton     

Current understanding of UV-induced base pair substitution mutation in E. coli with particular reference to the DNA polymerase III complex

UV mutagenesis in E. coli is believed to occur in two discrete steps. The second step involves continued DNA synthesis beyond a blocking lesion in the template strand. This bypass step requires induced levels of umuD and umuC gene products and activated recA protein. DNA polymerase III may be involved since a dnaE mutator strain (believed to have defective base selection) is associated with enhanced UV mutagenesis in conjunction with a genetic background permitting the bypass step. In non-UV-mutable umu and lexA strains, UV mutagenesis can be demonstrated if delayed photorevesal is given. This is interpreted as indicating that an earlier misincorporation step can occur in such strains but the resulting mutations do not survive because the bypass step is blocked. The misincorporation step does not require any induced SOS gene products and can occur either at the replication fork or during repair replication following excision of a DNA lesion. Neither a dnaE mutator gene (leading to a defective subunit of DNA polymerase III holoenzyme) nor a mutD5 mutator gene (leading to a defective ε proofreading subunit) had any effect on he misincorporation step. Although this is consistent with DNA polymerase III holoenzyme not being involved in the misincorporation step, other interpretations involving the inhibition of ε proofreading activity by recA protein are possible.

In vitro studies are reported in which sites of termination of synthesis by DNA polymerase III holoenzyme on UV-irradiated M13 mp8 DNA were examined in the presence of inhibitors of the 3′–5′ proofreading exonuclease (including recA protein). No evidence was found for incorporation of bases opposite photoproducts suggesting that either inhibition is more complete in the cell and/or that other factors are involved in the misincorporation step.     

Modelling mitochondrial site polymorphisms to infer the number of segregating units and mutation rate

We present a mathematical model of mitochondrial inheritance evolving under neutral evolution to interpret the heteroplasmies observed at some sites. A comparison of the levels of heteroplasmies transmitted from mother to her offspring allows us to estimate the number N_x of inherited mitochondrial genomes (segregating units). The model demonstrates the necessity of accounting for both the multiplicity of an unknown number N_x, and the threshold θ, below which heteroplasmy cannot be detected reliably, in order to estimate the mitochondrial mutation rate μ_m in the maternal line of descent. Our model is applicable to pedigree studies of any eukaryotic species where site heteroplasmies are observed in regions of the mitochondria, provided neutrality can be assumed. The model is illustrated with an analysis of site heteroplasmies in the first hypervariable region of mitochondrial sequence data sampled from Adélie penguin families, providing an estimate N_x and μ_m. This estimate of μ_m was found to be consistent with earlier estimates from ancient DNA analysis.     

The application of molecular markers to the study and conservation of fish populations, with special reference to Salmo

The main molecular techniques which can be used to generate genetic markers, and the applications of these markers to studies of fish populations are outlined. Published and ongoing studies, in the authors' laboratories, on brown trout and Atlantic salmon are used to compare the resolution and applicability of allozyme, mitochondrial DNA and minisatellite (variable number of tandem repeats) markers for studies on population structuring, genetic variation within populations, and the impact of the accidental and deliberate introduction of non-native salmonids on the genetic make-up of natural populations.     

Restriction map of Chinese hamster mitochondrial DNA containing replication coordinates: comparison with Syrian hamster mitochondrial genome

A precise physical map, containing the structurally and operationally defined D-loop origin, terminal region, and direction of heavy-strand replication, has been constructed for mitochondrial DNA (mtDNA) from ovary (CHO-KI) and lung cells of Chinese hamster (Cricetulus griseus 2 N = 22), and compared with our previously established genome coordinates for mtDNA from Syrian hamster ( Mesocricetus auratus 2 N = 44). All four HpaI sites in Cricetulus are conserved in Mesocricetus (8 sites). Extensive variation exists for hexanucleotides cleaved by EcoRI HindIII PstI. KpnI and BamHI. Sequence divergence between Chinese and Syrian hamster mtDNAs, as reflected from analysis of the mapped recognition sites for these six endonucleases, is estimated as 5–9% base substitutions. mtDNAs from both hamster and several other mammalian species contain a commonly conserved HpaI site in the region of light strand initiation.     

Drosophila mitochondrial DNA: Conserved sequences in the A+T-rich region and supporting evidence for a secondary structure model of the small ribosomal RNA

Summary The sequence of a segment of theDrosophila virilis mitochondrial DNA (mtDNA) molecule that contains the A+T-rich region, the small rRNA gene, the tRNA^f-met, tRNA^gln, and tRNA^ile genes, and portions of the ND2 and tRNA^val genes is presented and compared with the corresponding segment of theD. yakuba mtDNA molecule. The A+T-rich regions ofD. virilis andD. yakuba contain two correspondingly located sequences of 49 and 276/274 nucleotides that appear to have been conserved during evolution. In each species the replication origin of the mtDNA molecule is calculated to lie within a region that overlaps the larger conserved sequence, and within this overlap is found a potential hairpin structure. Substitutions between the larger conserved sequences of the A+T-rich regions, the small mt-rRNA genes, and the ND2 genes are biased in favor of transversions, 71–97% of which are AT changes. There is a 13.8 times higher frequency of nucleotide differences between the 5 halves than between the 3 halves of theD. virilis andD. yakuba small mt-rRNA genes. Considerations of the effects of observed substitutions and deletion/insertions on possible nucleotide pairing within the small mt-rRNA genes ofD. virilis andD. yakuba strongly support the secondary structure model for theDrosophila small mt-rRNA that we previously proposed.     

A molecular perspective on the phylogeny of placental mammals based on mitochondrial 12S rDNA sequences,with special reference to the problem of the Paenungulata

Part of the 12S rDNA gene was amplified and sequenced for 11 placental mammals, 3 marsupials, and 2 monotremes. Multiple alignments for these sequences and nine additional placental sequences taken from GenBank were obtained using CLUSTAL. Phylogenetic analyses were performed using standard parsimony, transversion parsimony, and Lake's method of invariants. All of our analyses uniteLoxodontia withDugong. Procavia, in turn, is a sister group to these taxa, thus supporting the monophyly of the Paenungulata. Perissodactyls are a sister group to paenungulates when transitions and transversions are both included but not when transitions are omitted. Likewise, cetaceans are a sister group to artiodactyls on minimum length trees under standard parsimony but not under transversion parsimony. Rodent monophyly and bat monophyly also receive mixed support, as does a putative alliance between primates and lagomorphs. Interestingly, the percentage divergence between the echidna and the platypus is less than for the rat and mouse.     

The Dictyostelium discoideum mitochondrial genome: A primordial system using the universal code and encoding hydrophilic proteins atypical of metazoan mitochondrial DNA

A 3,345-bp fragment of Dictyostelium discoideum mitochondrial DNA (mtDNA) has been sequenced. This fragment contained the 80-kDa subunit of complex I (NADH:ubiquinone oxidoreductase), encoding a predicted amino acid sequence of 688 residues and a molecular mass of 79,805 daltons which is nuclear encoded in other metazoa. The C-terminus of the D. discoideum complex I gene shared a 10-bp overlap with NADH:ubiquinone oxidoreductase chain 5 (ND5), while 21 by 5 were three tRNA genes (two isoleucine and a histidine) and a further 25 by 5 of these genes is the partial sequence (104 residues) of an unidentified open reading frame (ORF104). Both the 80-kDa subunit and the ORF104 were hydrophilic and highly charged, suggesting they are not membrane associated, unlike most mitochondrially encoded proteins in the metazoa. Sequence analysis of the 80-kDa subunit, its adjacent ND5 gene, and ORF104 indicates the universal stop codon TGA, which codes for tryptophan in nearly all nonplant mtDNA, is either unassigned or coding for a stop codon in D. discoideum. The large size of the mitochondrial genome (54 kb), the lack of intergenic sequence, and the apparent use of the universal code suggest D. discoideum mtDNA may encode many primitive genes that are nuclear encoded in higher organisms.Correspondence to: K.L. WilliamsData deposition: GeneBank     

Estimating the distribution of selection coefficients from phylogenetic data with applications to mitochondrial and viral DNA

The distribution of selection coefficients of new mutations is of key interest in population genetics. In this paper we explore how codon-based likelihood models can be used to estimate the distribution of selection coefficients of new amino acid replacement mutations from phylogenetic data. To obtain such estimates we assume that all mutations at the same site have the same selection coefficient. We first estimate the distribution of selection coefficients from two large viral data sets under the assumption that the viral population size is the same along all lineages of the phylogeny and that the selection coefficients vary among sites. We then implement several new models in which the lineages of the phylogeny may have different population sizes. We apply the new models to a data set consisting of the coding regions from eight primate mitochondrial genomes. The results suggest that there might be little power to determine the exact shape of the distribution of selection coefficient but that the normal and gamma distributions fit the data significantly better than the exponential distribution.     

Current Intraspecific Dynamics of Sequence Evolution Differs from Long-Term Trends and Can Account for the AT-Richness of Honeybee Mitochondrial DNA

The very high AT content of hymenopteran mtDNA has warranted speculation about nucleotide substitution processes in this group. Here we investigate the pattern of honeybee, Apis mellifera, mtDNA nucleotide polymorphisms inferred from phylogeny in terms of differences between the ATPase6, COI, COII, COIII, cytochrome b, and ND2 genes and strand asymmetry in mutation rates. The observed transition/transversion ratios and the distribution of nonsynonymous substitutions between regions differed significantly. The pattern of differences between genes leading to these heterogeneities (the ATPase6 and COIII genes group apart from the rest) differed markedly from that predicted on the basis of long-term evolutionary change and may indicate differences between current and long-term dynamics of sequence evolution. Also, there is strong strand asymmetry in substitutions, which probably results in a mutability of G and C sufficiently high to account for the AT-richness of honeybee mtDNA. Received: 21 October 1998 / Accepted: 27 January 1999     

The complete mitochondrial genome of the leafminer Liriomyza sativae (Diptera: Agromyzidae): great difference in the A+T-rich region compared to Liriomyza trifolii

The complete mitochondrial genome sequence of Liriomyza sativae Blanchard (15,551 bp) was determined and analyzed in this study. The circular genome contained 37 genes including 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes and an A + T-rich region. The initiation codons of COI and ND1 were ‘ATCA’ and ‘GTG’, respectively. ND2 gene used the truncated termination codon ‘T’. All the tRNA genes had the typical cloverleaf secondary structures except for tRNA^Ser(AGN) gene, which was found with the absence of a DHU arm. In addition, a tRNA-like secondary structure (tRNA^Met) was found in the A + T-rich region. The great difference was that the length of L. sativae A + T-rich region was 597 bp shorter than that of Liriomyza trifolii (Burgess). Meanwhile, some minor differences such as ‘TATA’ block were also observed in L. sativae in contrast to ‘TACA’ block in L. trifolii. There were also some essential structure elements such as ‘TATA’ block, ‘G(A)_nT’ block, poly-T stretch and stem-and-loop structure in the A + T-rich region of L. sativae mitochondrial genome.     

Vector competence of South African Culicoides species for bluetongue virus serotype 1 (BTV-1) with special reference to the effect of temperature on the rate of virus replication in C. imicola and C. bolitinos

Abstract. The oral susceptibility of 22 South African livestock associated Culicoides species to infection with bluetongue virus serotype 1 (BTV‐1) and its replication rate in C. imicola Kieffer and C. bolitinos Meiswinkel (Diptera: Ceratopogonidae) over a range of different incubation periods and temperatures are reported. Field‐collected Culicoides were fed on sheep blood containing 7.5 log₁₀TCID₅₀/mL of BTV‐1, and then held at constant different temperatures. Virus replication was measured over time by assaying individual flies in BHK‐21 cells using a microtitration procedure. Regardless of the incubation temperatures (10, 15, 18, 23.5 and 30°C) the mean virus titre/midge, infection rates (IR) and the proportion of infected females with transmission potential (TP = virus titre/midge ≥ 3 log₁₀ TCID₅₀) were found to be significantly higher in C. bolitinos than in C. imicola. Results from days 4–10 post‐infection (dpi), at 15–30°C, shows that the mean IR and TP values in C. bolitinos ranged from 36.7 to 87.8%, and from 8.4 to 87.7%, respectively; in C. imicola the respective values were 11.0–13.7% and 0–46.8%. In both species the highest IR was recorded at 25°C and the highest TP at 30°C. The time required for the development of TP in C. bolitinos ranged from 2 dpi at 25°C to 8 dpi at 15°C. In C. imicola it ranged from 4 dpi at 30°C to 10 dpi at 23.5°C; no individuals with TP were detected at 15°C. There was no evidence of virus replication in flies held at 10°C. When, at various points of incubation, individual flies were transferred from 10°C to 23.5°C and then assayed 4–10 days later, virus was recovered from both species. The mean virus titres/midge, and proportion of individuals with TP and IR, were again significantly higher in C. bolitinos than in C. imicola. Also the infection prevalence in C. magnus Colaço was higher than in C. imicola. Low infection prevalences were found in C. bedfordi Ingram & Macfie, C. leucostictus Kieffer, C. pycnostictus Ingram & Macfie, C. gulbenkiani Caeiro and C. milnei Austen. BTV‐1 was not detected in 14 other Culicoides species tested; however, some of these were tested in limited numbers. The present study indicates a multivector potential for BTV transmission in South Africa. In C. imicola and C. bolitinos the replication rates are distinct and are significantly influenced by temperature. These findings are discussed in relation to the epidemiology of bluetongue in South Africa.     

Characteristics of Nucleotide Substitution in the Hepatitis C Virus Genome: Constraints on Sequence Change in Coding Regions at Both Ends of the Genome

Comparison of complete genome sequences for different variants of hepatitis C virus (HCV) reveals several different constraints on sequence change. Synonymous changes are suppressed in coding regions at both 5′ and 3′ ends of the genome. No evidence was found for the existence of alternative reading frames or for a lower mutation frequency in these regions. Instead, suppression may be due to constraints imposed by RNA secondary structures identified within the core and NS5b genes. Nonsynonymous substitutions are less frequent than synonymous ones except in the hypervariable region of E2 and, to a lesser extent, in E1, NS2, and NS5b. Transitions are more frequent than transversions, particularly at the third position of codons where the bias is 16:1. In addition, nucleotide substitutions may not occur symmetrically since there is a bias toward G or C at the third position of codons, while T ↔ C transitions were twice as frequent as A ↔ G transitions. These different biases do not affect the phylogenetic analysis of HCV variants but need to be taken into account in interpreting sequence change in longitudinal studies. Received: 9 September 1996 / Accepted: 20 April 1997     

Complete mitochondrial genomes from four subspecies of common chaffinch (Fringilla coelebs): New inferences about mitochondrial rate heterogeneity,neutral theory,and phylogenetic relationships within the order Passeriformes

We describe whole mitochondrial genome sequences from four subspecies of the common chaffinch (Fringilla coelebs), and compare them to 31 publicly available mitochondrial genome sequences from other Passeriformes. Rates and patterns of mitochondrial gene evolution are analyzed at different taxonomic levels within this avian order, and evidence is adduced for and against the nearly neutral theory of molecular evolution and the role of positive selection in shaping genetic variation of this small but critical genome. We find evidence of mitochondrial rate heterogeneity in birds as in other vertebrates, likely due to differences in mutational pressure across the genome. Unlike in gadine fish and some of the human mitochondrial work we do not observe strong support for the nearly neutral theory of molecular evolution; instead evidence from molecular clocks, distribution of dN/dS ratios at different levels of the taxonomic hierarchy and in different lineages, McDonald–Kreitman tests within Fringillidae, and site-specific tests of selection within Passeriformes, all point to a role for positive selection, especially for the complex I NADH dehydrogenase genes. The protein-coding mitogenome phylogeny of the order Passeriformes is broadly consistent with previously-reported molecular findings, but provides support for a sister relationship between the superfamilies Muscicapoidea and Passeroidea on a short basal internode of the Passerida where relationships have been difficult to resolve. An unexpected placement of the Paridae (represented by Hume's groundpecker) within the Muscicapoidea was observed. Consistent with other molecular studies the mtDNA phylogeny reveals paraphyly within the Muscicapoidea and a sister relationship of Fringilla with Carduelis rather than Emberiza.