Does effective population size affect rates of molecular evolution: Mitochondrial data for host/parasite species pairs in bees suggests not

Adaptive evolutionary theory argues that organisms with larger effective population size (N e) should have higher rates of adaptive evolution and therefore greater capacity to win evolutionary arm races. However, in some certain cases, species with much smaller N e may be able to survive besides their opponents for an extensive evolutionary time. Neutral theory predicts that accelerated rates of molecular evolution in organisms with exceedingly small N e are due to the effects of genetic drift and fixation of slightly deleterious mutations. We test this prediction in two obligate social parasite species and their respective host species from the bee tribe Allodapini. The parasites (genus Inquilina) have been locked into tight coevolutionary arm races with their exclusive hosts (genus Exoneura) for ~15 million years, even though Inquilina exhibit N e that are an order of magnitude smaller than their host. In this study, we compared rates of molecular evolution between host and parasite using nonsynonymous to synonymous substitution rate ratios (dN/dS) of eleven mitochondrial protein-coding genes sequenced from transcriptomes. Tests of selection on mitochondrial genes indicated no significant differences between host and parasite dN/dS, with evidence for purifying selection acting on all mitochondrial genes of host and parasite species. Several potential factors which could weaken the inverse relationship between N e and rate of molecular evolution are discussed.

Inferring Evolutionary Timescales without Independent Timing Information: An Assessment of "Universal" Insect Rates to Calibrate a Collembola (Hexapoda) Molecular Clock

Previous estimates of nucleotide substitution rates are routinely applied as secondary or “universal” molecular clock calibrations for estimating evolutionary timescales in groups that lack independent timing information. A major limitation of this approach is that rates can vary considerably among taxonomic groups, but the assumption of rate constancy is rarely evaluated prior to using secondary rate calibrations. Here I evaluate whether an insect mitochondrial DNA clock is appropriate for estimating timescales in Collembola—a group of insect-like arthropods characterized by high levels of cryptic diversity. Relative rates of substitution in cytochrome oxidase subunit 1 (COI) were inferred via Bayesian analysis across a topologically constrained Hexapod phylogeny using a relaxed molecular clock model. Rates for Collembola did not differ significantly from the average rate or from the rates estimated for most other groups (25 of 30), suggesting that (1) their apparent cryptic diversity cannot be explained by accelerated rates of molecular evolution and (2) clocks calibrated using “universal” insect rates may be appropriate for estimating evolutionary timescales in this group. However, of the 31 groups investigated, 10 had rates that deviated significantly from the average (6 higher, 4 lower), underscoring the need for caution and careful consideration when applying secondary insect rate calibrations. Lastly, this study exemplifies a relatively simple approach for evaluating rate constancy within a taxonomic group to determine whether the use of secondary rates are appropriate for molecular clock calibrations.

Conserved numts mask a highly divergent mitochondrial-COI gene in a species complex of Australian stingless bees Tetragonula (Hymenoptera: Apidae)

Tetragonula carbonaria, Tetragonula davenporti, Tetragonula hockingsi and Tetragonula mellipes comprise a species complex of Australian stingless bee species known as the 'Carbonaria' group. The species are difficult to distinguish morphologically and the major species-defining characters relate to comb architecture and nest entrance ornamentation. The taxonomy of the group is further complicated by likely nuclear mitochondrial pseudogenes (numts) and inter-specific hybrids. Here we demonstrate the existence of COI numts and isolate and characterize the 'true' mt-COI gene in T. carbonaria and T. hockingsi. Numts were isolated from enriched-nuclear DNA extraction followed by PCR amplification and Sanger sequencing, and were recognized by the presence of deletions and/or premature stop codons in the translated sequences. The mt-COI sequences were obtained from NGS sequencing using purified mtDNA. In T. carbonaria, two numts (numt1 and numt2) were identified and a third (numt3) was identified in T. hockingsi. Numt2 and numt3 are similar (1.2% sequence divergence), indicating a recent common origin. The genetic distance between the mt-COI of the two Tetragonula species was higher than might be expected for closely related species, 16.5%, corroborating previous studies in which T. carbonaria and T. hockingsi were regarded as separate species. The three numts are more similar to the COI of other stingless bee species, including Australian Austroplebia australis and South American Melipona bicolor (81.7-83.9%) than to the mt-COI of their own species (70-71.4%). This is because the mt-COI of T. carbonaria and T. hockingsi differ greatly from other Meliponinae. Our findings explain some formerly puzzling aspects of Carbonaria biogeography, and misinterpreted amplifications.

Founder effects initiated rapid species radiation in Hawaiian cave planthoppers

The Hawaiian Islands provide the venue of one of nature's grand experiments in evolution. Here, we present morphological, behavioral, genetic, and geologic data from a young subterranean insect lineage in lava tube caves on Hawai'i Island. The Oliarus polyphemus species complex has the potential to become a model for studying rapid speciation by stochastic events. All species in this lineage live in extremely similar environments but show strong differentiation in behavioral and morphometric characters, which are random with respect to cave age and geographic distribution. Our observation that phenotypic variability within populations decreases with increasing cave age challenges traditional views on founder effects. Furthermore, these cave populations are natural replicates that can be used to test the contradictory hypotheses. Moreover, Hawaiian cave planthoppers exhibit one of the highest speciation rates among animals and, thus, radically shift our perception on the evolutionary potential of obligate cavernicoles.

Examining the role of effective population size on mitochondrial and multilocus divergence time discordance in a songbird

Estimates of speciation times are subject to a number of potential errors. One source of bias is that effective population size (Ne) has been shown to influence substitution rates. This issue is of particular interest for phylogeographic studies because population sizes can vary dramatically among genetically structured populations across species' ranges. In this study, we used multilocus data to examine temporal phylogeographic patterns in a widespread North American songbird, the Northern Cardinal (Cardinalis cardinalis). Species tree estimation indicated that the phylogeographic structure of C. cardinalis was comprised of four well-supported mainland lineages with large population sizes (large Ne) and two island lineages comprised of much smaller populations (small Ne). We inferred speciation times from mtDNA and multilocus data and found there was discordance between events that represented island-mainland divergences, whereas both estimates were similar for divergences among mainland lineages. We performed coalescent simulations and found that the difference in speciation times could be attributed to stochasticity for a recently diverged island lineage. However, the magnitude of the change between speciation times estimated from mtDNA and multilocus data of an older island lineage was substantially greater than predicted by coalescent simulations. For this divergence, we found the discordance in time estimates was due to a substantial increase in the mtDNA substitution rate in the small island population. These findings indicate that in phylogeographic studies the relative tempo of evolution between mtDNA and nuclear DNA can become highly discordant in small populations.

Testing the effect of metabolic rate on DNA variability at the intra-specific level

We tested the metabolic rate hypothesis (whereby rates of mtDNA evolution are postulated to be mediated primarily by mutagenic by-products of respiration) by examining whether mass-specific metabolic rate was correlated with root-to-tip distance on a set of mtDNA trees for the springtail Cryptopygus antarcticus travei from sub-Antarctic Marion Island.Using Bayesian analyses and a novel application of the comparative phylogenetic method, we did not find significant evidence that contemporary metabolic rates directly correlate with mutation rate (i.e., root-to-tip distance) once the underlying phylogeny is taken into account. However, we did find significant evidence that metabolic rate is dependent on the underlying mtDNA tree, or in other words, lineages with related mtDNA also have similar metabolic rates.We anticipate that future analyses which apply this methodology to datasets with longer sequences, more taxa, or greater variability will have more power to detect a significant direct correlation between metabolic rate and mutation rate. We conclude with suggestions for future analyses that would extend the preliminary approach applied here, in particular highlighting ways to tease apart oxidative stress effects from the effects of population size and/or selection coefficients operating on the molecular evolutionary rate.

Near neutrality: leading edge of the neutral theory of molecular evolution

The nearly neutral theory represents a development of Kimura's neutral theory of molecular evolution that makes testable predictions that go beyond a mere null model. Recent evidence has strongly supported several of these predictions, including the prediction that slightly deleterious variants will accumulate in a species that has undergone a severe bottleneck or in cases where recombination is reduced or absent. Because bottlenecks often occur in speciation and slightly deleterious mutations in coding regions will usually be nonsynonymous, we should expect that the ratio of nonsynonymous to synonymous fixed differences between species should often exceed the ratio of nonsynonymous to synonymous polymorphisms within species. Many data support this prediction, although they have often been wrongly interpreted as evidence for positive Darwinian selection. The use of conceptually flawed tests for positive selection has become widespread in recent years, seriously harming the quest for an understanding of genome evolution. When properly analyzed, many (probably most) claimed cases of positive selection will turn out to involve the fixation of slightly deleterious mutations by genetic drift in bottlenecked populations. Slightly deleterious variants are a transient feature of evolution in the long term, but they have substantially affected contemporary species, including our own.

Molecular and quantitative trait variation across the native range of the invasive speciesHypericum canariense: evidence for ancient patterns of colonization via pre-adaptation?

To understand the success of invasive species, it is important to know whether colonization events are facilitated by adaptive evolution or are limited to sites where a species is pre-adapted to thrive. Studies of the ancient colonization patterns of an invader in its native range provide an opportunity to examine its natural history of adaptation and colonization. This study uses molecular (internal transcribed spacer sequence and amplified fragment length polymorphism) and common garden approaches to assess the ancient patterns of establishment and quantitative trait evolution in the invasive shrub Hypericum canariense. This species has an unusually small and discrete native range in the Canary Islands. Our data reveal two genetic varieties with divergent life histories and different colonization patterns across the islands. Although molecular divergence within each variety is large (pairwise FST from 0.18 to 0.32 between islands) and nearly as great as divergence between them, life-history traits show striking uniformity within varieties. The discrepancy between molecular and life-history trait divergence points to the action of stabilizing selection within varieties and the influence of pre-adaptation on patterns of colonization. The colonization history of H. canariense reflects how the relationship between selective environments in founding and source populations can dictate establishment by particular lineages and their subsequent evolutionary stasis or change.

Looking for Darwin in all the wrong places: the misguided quest for positive selection at the nucleotide sequence level

Recent years have seen an explosion of interest in evidence for positive Darwinian selection at the molecular level. This quest has been hampered by the use of statistical methods that fail adequately to rule out alternative hypotheses, particularly the relaxation of purifying selection and the effects of population bottlenecks, during which the effectiveness of purifying selection is reduced. A further problem has been the assumption that positive selection will generally involve repeated amino-acid changes to a single protein. This model was derived from the case of the vertebrate major histocompatibility complex (MHC), but the MHC proteins are unusual in being involved in protein-protein recognition and in a co-evolutionary process of pathogens. There is no reason to suppose that repeated amino-acid changes to a single protein are involved in selectively advantageous phenotypes in general. Rather adaptive phenotypes are much more likely to result from other causes, including single amino-acid changes; deletion or silencing of genes or changes in the pattern of gene expression.

Comparative phylogeography of coastal limpets across a marine disjunction in New Zealand

Cook Strait, which separates the North and South Island of New Zealand, has been a transient, but re-occurring feature of the New Zealand land mass throughout the Pleistocene, maintaining its current width and depth for the past 5000 years. Historic land fragmentation coupled with the complex hydrography of the Greater Cook Strait region has created both biogeographic and phylogeographic disjunctions between the North and South Island in several marine species. Here we use mitochondrial cytochrome b DNA sequences of three endemic intertidal limpets, Cellana ornata, Cellana radians and Cellana flava to assess intraspecific phylogeographic patterns across Cook Strait and to look for interspecific concordance of ecological and evolutionary processes among closely related taxa. We sequenced 328-359 bp in 85-321 individuals from 8-31 populations spanning the biogeographic range of the three species. Intraspecific phylogeographic analyses show moderate to strong genetic discontinuity among North and South Island populations due to allopatric fragmentation. This pattern was broadly concordant across the three species and the observed divergence among this group of intertidal limpets (0.3-2.0%) is similar to that of previously studied subtidal organisms. For each species, divergence time calculations suggest contemporary North and South Island lineages diverged from their respective most recent common ancestor approximately 200 000 to 300 000 years before present (bp), significantly earlier than previous estimates in other coastal marine taxa that arose from a miscalculation of divergence time.

Higher-order organization and compartmentalization of satellite DNA PIM357 in species of the coleopteran genus Pimelia

The PIM357 satellite DNA family is present in 26 Pimelia taxa (Tenebrionidae, Coleoptera) with endemic congeneric species from the Canary Islands showing higher interrepeat variability than continental ones. In this paper, we compare the repetitive DNA sequences of a Canarian species that has distinct subfamilies of repeat units, P. radula ascendens, with another without such subfamilies, P. sparsa sparsa. The chromosomal localization of the repeat units and the comparison of the variability of randomly cloned monomers to the one estimated by comparing repeat units from dimers and trimers suggest the absence of satellite subfamilies in P. sparsa sparsa. Hence, the repeat units of this species seem to be uniformly and randomly distributed throughout all chromosomes out of one chromosomal pair. On the contrary, P. radula ascendens shows four divergent subfamilies of repeat units supported by several diagnostic nucleotide substitutions. These subfamilies seem to form four distinct repeat units: monomer subfamily 1, monomer subfamily 4 and two higher-order units (dimer linking subfamily 1 and 4, and dimer linking subfamily 2 and 3). Moreover, monomers of subfamily 1 are present in three chromosomal pairs only. We discuss the effect of different potential factors acting in the concerted evolution and the genomic organization of stDNA sequences in these taxa.

Rates and patterns of molecular evolution in inbred and outbred Arabidopsis

The evolution of self-fertilization is associated with a large reduction in the effective rate of recombination and a corresponding decline in effective population size. If many spontaneous mutations are slightly deleterious, this shift in the breeding system is expected to lead to a reduced efficacy of natural selection and genome-wide changes in the rates of molecular evolution. Here, we investigate the effects of the breeding system on molecular evolution in the highly self-fertilizing plant Arabidopsis thaliana by comparing its coding and noncoding genomic regions with those of its close outcrossing relative, the self-incompatible A. lyrata. More distantly related species in the Brassicaceae are used as outgroups to polarize the substitutions along each lineage. In contrast to expectations, no significant difference in the rates of protein evolution is observed between selfing and outcrossing Arabidopsis species. Similarly, no consistent overall difference in codon bias is observed between the species, although for low-biased genes A. lyrata shows significantly higher major codon usage. There is also evidence of intron size evolution in A. thaliana, which has consistently smaller introns than its outcrossing congener, potentially reflecting directional selection on intron size. The results are discussed in the context of heterogeneity in selection coefficients across loci and the effects of life history and population structure on rates of molecular evolution. Using estimates of substitution rates in coding regions and approximate estimates of divergence and generation times, the genomic deleterious mutation rate (U) for amino acid substitutions in Arabidopsis is estimated to be approximately 0.2-0.6 per generation.

Changes in the recombinational environment affect divergence in the yellow gene of Drosophila

The complete coding region of the yellow (y) gene was sequenced in different Drosophila species. In the species of the melanogaster subgroup (D. melanogaster, D. simulans, D. mauritiana, D. yakuba, and D. erecta), this gene is located at the tip of the X chromosome in a region with a strong reduction in recombination rate. In contrast, in D. ananassae (included in the ananassae subgroup of the melanogaster group) and in the obscura group species (D. subobscura, D. madeirensis, D. guanche, and D. pseudoobscura), the y gene is located in regions with normal recombination rates. As predicted by the hitchhiking and background selection models, this change in the recombinational environment affected synonymous divergence in the y-gene-coding region. Estimates of the number of synonymous substitutions per site were much lower between the obscura group species and D. ananassae than between the species of the obscura group and the melanogaster subgroup. In fact, a highly significant increase in the rate of synonymous substitution was detected in all lineages leading to the species of the melanogaster subgroup relative to the D. ananassae lineage. This increase can be explained by a higher fixation rate of mutations from preferred to unpreferred codons (slightly deleterious mutations). The lower codon bias detected in all species of the melanogaster subgroup relative to D. ananassae (or to the obscura group species) would be consistent with this proposal. Therefore, at least in Drosophila, changes in the recombination rate in different lineages might cause deviations of the molecular-clock hypothesis and contribute to the overdispersion of the rate of synonymous substitution. In contrast, the change in the recombinational environment of the y gene has no detectable effect on the rate of amino acid replacement in the Yellow protein.

Elevated rates of nonsynonymous substitution in island birds

Slightly deleterious mutations are expected to fix at relatively higher rates in small populations than in large populations. Support for this prediction of the nearly-neutral theory of molecular evolution comes from many cases in which lineages inferred to differ in long-term average population size have different rates of nonsynonymous substitution. However, in most of these cases, the lineages differ in many other ways as well, leaving open the possibility that some factor other than population size might have caused the difference in substitution rates. We compared synonymous and nonsynonymous substitutions in the mitochondrial cyt b and ND2 genes of nine closely related island and mainland lineages of ducks and doves. We assumed that island taxa had smaller average population sizes than those of their mainland sister taxa for most of the time since they were established. In all nine cases, more nonsynonymous substitutions occurred on the island branch, but synonymous substitutions showed no significant bias. As in previous comparisons of this kind, the lineages with smaller populations might differ in other respects that tend to increase rates of nonsynonymous substitution, but here such differences are expected to be slight owing to the relatively recent origins of the island taxa. An examination of changes to apparently "preferred" and "unpreferred" synonymous codons revealed no consistent difference between island and mainland lineages.

Evolution on a volcanic conveyor belt: using phylogeographic reconstructions and K-Ar-based ages of the Hawaiian Islands to estimate molecular evolutionary rates

The Hawaiian Islands form as the Pacific Plate moves over a 'hot spot' in the earth's mantle where magma extrudes through the crust to build huge shield volcanos. The islands subside and erode as the plate carries them to the north-west, eventually to become coral atolls and seamounts. Thus islands are ordered linearly by age, with the oldest islands in the north-west (e.g. Kauai at 5.1 Ma) and the youngest in the south-east (e.g. Hawaii at 0.43 Ma). K-Ar estimates of the date of an island's formation provide a maximum age for the taxa inhabiting the island. These ages can be used to calibrate rates of molecular change under the following assumptions: (i) K-Ar dates are accurate; (ii) tree topologies show that derivation of taxa parallels the timing of island formation; (iii) populations do not colonize long after island emergence; (iv) the coalescent point for sister taxa does not greatly predate the formation of the colonized younger island; (v) saturation effects and (vi) among-lineage rate variation are minimal or correctable; and (vii) unbiased standard errors of distances and regressions can be estimated from multiple pairwise comparisons. We use the approach to obtain overall corrected rate calibrations for: (i) part of the mitochondrial cytochrome b gene in Hawaiian drepanidines (0.016 sequence divergence/Myr); (ii) the Yp1 gene in Hawaiian Drosophila (0.019/Myr Kambysellis et al. 1995); and (iii) parts of the mitochondrial 12S and 16S rRNA and tRNAval in Laupala crickets (0.024-0.102/Myr, Shaw 1996). We discuss the reliability of the estimates given the assumptions (i-vii) above and contrast the results with previous calibrations of Adh in Hawaiian Drosophila and chloroplast DNA in lobeliods.

Departure from neutrality at the mitochondrial NADH dehydrogenase subunit 2 gene in humans, but not in chimpanzees

To test whether patterns of mitochondrial DNA (mtDNA) variation are consistent with a neutral model of molecular evolution, nucleotide sequences were determined for the 1041 bp of the NADH dehydrogenase subunit 2 (ND2) gene in 20 geographically diverse humans and 20 common chimpanzees. Contingency tests of neutrality were performed using four mutational categories for the ND2 molecule: synonymous and nonsynonymous mutations in the transmembrane regions, and synonymous and nonsynonymous mutations in the surface regions. The following three topological mutational categories were also used: intraspecific tips, intraspecific interiors, and interspecific fixed differences. The analyses reveal a significantly greater number of nonsynonymous polymorphisms within human transmembrane regions than expected based on interspecific comparisons, and they are inconsistent with a neutral equilibrium model. This pattern of excess nonsynonymous polymorphism is not seen within chimpanzees. Statistical tests of neutrality, such as TAJIMA's D test, and the D and F tests proposed by FU and LI, indicate an excess of low frequency polymorphisms in the human data, but not in the chimpanzee data. This is consistent with recent directional selection, a population bottleneck or background selection of slightly deleterious mutations in human mtDNA samples. The analyses further support the idea that mitochondrial genome evolution is governed by selective forces that have the potential to affect its use as a "neutral" marker in evolutionary and population genetic studies.

Increased genetic diversity in mitochondrial genes is correlated with the evolution of parasitism in the Hymenoptera

A higher AT content and rate of mtDNA sequence divergence was found in parasitic wasps (Apocrita) compared with nonparasitic wasps (Symphyta). The compositional bias was reflected in extreme codon bias for a cytochrome oxidase I protein coding gene fragment as well as in the types of amino acid substitutions that have occurred during the evolution of this gene fragment. In some instances, compositional bias influenced the definition of a conservative amino acid change. The increased rate of mtDNA sequence evolution probably arose during the early Jurassic, coincident with the first appearance of parasitic wasps in the fossil record. Our results suggest a causal link between the rate of sequence divergence and the parasitic lifestyle.

Mitochondrial DNA phylogeny and sequential colonization of Canary Islands by darkling beetles of the genus Pimelia (Tenebrionidae)

Fifteen species of the darkling beetle genus Pimelia (Tenebrionidae: Coleoptera) have been sequenced for a 365 b.p. portion of the mitochondrial Cytochrome Oxidase I gene. Thirteen of these are endemic species inhibiting the Canarian archipelago in the Atlantic and the other two are continental conspecifics. This data set has been phylogenetically analysed by maximum parsimony and distance approaches and the resulting trees used to deduce sequential interisland colonization. This suggests patterns compatible with the geological dating of the islands, but with increasing uncertainty when older events are considered. A colonization sequence from Fuerteventura to Tenerife followed by Tenerife to Gran Canarian and La Gomera, and from the latter to La Palma and then to El Hierro is proposed for the genus. A relatively recent secondary colonization from Gran Canaria to Gomera is deduced.

ADH evolution and the phylogenetic footprint

The evolution of any given protein reflects the interplay between proximal selective forces involving the conservation of protein structure and function and more general populational factors that shape the action and efficiency of natural selection. In an attempt to address that interplay, we have analyzed patterns of amino acid replacement within a well-conserved molecule, alcohol dehydrogenase (ADH), in the Drosophilidae. A sliding window, moved along the protein sequence in order to quantify the extent of change at each amino acid position, reveals heterogeneous amounts of replacement across the molecule when all ADH sequences are analyzed simultaneously. Surprisingly, the replacement profile for ADH differs significantly in the melanogaster, mulleri, and Hawaiian subgroups, reflecting the imprint of the differing evolutionary histories of each of these assemblages on the evolution of this conservative molecule.

Organ and cell allometry in Hawaiian Drosophila: how to make a big fly

The importance of body size in predicting many aspects of an animal's biology has become well established in recent years. However, little is known about how body size evolves at the cellular level. Some published data suggest that it is cell number and not cell size that accompanies changes in organ and body size across taxa. We examined organ and cell allometry in the wing, eye and basitarsus of adult Hawaiian Drosophila, ranging in body length from 0.2 mm to 0.8 mm. Linear measurements of all three structures exhibit a positive allometry with body length. Exponents of the allometric equation were 0.96, 0.55 and 1.50 for wing, eye and basitarsus, respectively. Surface markers were used to quantify cell size of each organ. The allometric exponents for cell size as a function of organ size were 0.53, 0.68 and 0.33 for wing, eye and basitarsus, respectively. In contrast to reports in the literature on other systems, our results for Hawaiian Drosophila indicate that cell size may contribute between one third and two thirds to evolutionary changes in organ and body size.

Evolutionary analysis of cytochrome b sequences in some Perciformes: evidence for a slower rate of evolution than in mammals

To obtain information relative to the phylogenesis and microevolutionary rate of fish mitochondrial DNA, the nucleotide sequence of cytochrome b gene in seven fish species belonging to the order of Perciformes was determined. Sequence analysis showed that fish mitochondrial DNA has a nucleotide compositional bias similar to that of sharks but lower compared to mammals and birds. Quantitative evolutionary analysis, carried out by using a markovian stochastic model, clarifies some phylogenetic relationships within the Perciformes order, particularly in the Scombridae family, and between Perciformes, Gadiformes, Cypriniformes, and Acipenseriformes. The molecular clock of mitochondrial DNA was calibrated with the nucleotide substitution rate of cytochrome b gene in five shark species having divergence times inferred from paleontological estimates. The results of such analysis showed that Acipenseriformes diverged from Perciformes by about 200 MY, that the Perciformes common ancestor dates back to 150 MY, and that fish mitochondrial DNA has a nucleotide substitution rate three to five times lower than that of mammals.

Naja siamensis, a cryptic species of venomous snake revealed by mtDNA sequencing

Because of possible variation in venom composition, an understanding of venomous snake systematics is of great importance for the optimization of antivenom treatment of snakebite patients. Intraspecific variation in the morphology of many venomous snakes complicates the definition and identification of some species when allopatric populations are involved. Selectively neutral or near-neutral mtDNA sequences can reveal evolutionary relationships obscured by ecogenetically-caused morphological variation. We use comparative sequencing of the cytochrome oxidase subunit 1 gene to reveal the existence of a widespread, cryptic species of spitting cobra from southeast Asia. This species, Naja siamensis, is widely sympatric with other Asiatic cobra species. This may be of considerable medical significance, and calls for further research into venom composition in Asiatic cobras.

Complementary DNA-DNA hybridization in Drosophila

We have performed DNA-DNA hybridization experiments among several species of Drosophila using the evolutionarily conserved portion of the genome representing sequences coding for amino acids of proteins. This was done by using as tracer, radioactively labeled complementary DNA that was reverse transcribed from adult mRNA. We show that this procedure extends phylogenetically the distance over which the technique can be applied to fast-evolving groups such as Drosophila. The major phylogenetic conclusions are (1) the subgenus Sophophora is a monophyletic lineage; (2) within Sophophora the melanogaster subgroup is closer to the obscura group than either group is to the willistoni group; (3) the subgenus Drosophila is complex with most major lineages originating deep in the phylogeny; the subgenus may not be monophyletic; (4) as with most groups classically placed in Drosophila, the Hawaiian Drosophila originate early, supporting the notion that this lineage is older than the extant islands; and (5) the virilis/repleta lineage is monophyletic within Drosophila.

Evolution of the autosomal chorion cluster in Drosophila. IV. The Hawaiian Drosophila: rapid protein evolution and constancy in the rate of DNA divergence

Autosomal chorion genes s18, s15, and s19 are shown to diverge at extremely rapid rates in closely related taxa of Hawaiian Drosophila. Their nucleotide divergence rates are at least as fast as those of intergenic regions that are known to evolve more extensively between distantly related species. Their amino acid divergence rates are the fastest known to date. There are two nucleotide replacement substitutions for every synonymous one. The molecular basis for observed length and substitution mutations is analyzed. Length mutations are strongly associated with direct repeats in general, and with tandem repeats in particular, whereas the rate for an average transition is twice that for an average transversion. The DNA sequence of the cluster was used to construct a phylogenetic tree for five taxa of the Hawaiian picture-winged species group of Drosophila. Assignment of observed base substitutions occurring in various branches of the tree reveals an excess of would-be homoplasies in a centrally localized 1.8-kb segment containing the s15 gene. This observation may be a reflection of ancestral excess polymorphisms in the segment. The chorion cluster appears to evolve at a constant rate regardless of whether the central 1.8-kb segment is included or not in the analysis. Assuming that the time of divergence of Drosophila grimshawi and the planitibia subgroup coincides with the emergence of the island of Kauai, the overall rate of base substitution in the cluster is estimated to be 0.8% million years, whereas synonymous sites are substituted at a rate of 1.2% million years.

Evolution of human mitochondrial DNA: evidence for departure from a pure neutral model of populations at equilibrium

Human mitochondrial DNA (mt-DNA) data from 18 populations have been carefully reexamined. A phylogeny of 77 mtDNA types found among the 1389 individuals analyzed for restriction fragment length polymorphisms (RFLPs) was established using the parsimony principle and compared to a UPGMA tree of the 18 populations. Both analyses agreed in separating African samples from the other populations, though the mtDNA type phylogeny suggested close relations between Africans and other continental groups. Conformity of observed mtDNA type frequency distributions with the "infinite allele" model was studied for 31 human populations. Several Oriental and Caucasoid populations were found to be overly homogeneous, generally due to an elevated frequency of one particular type. Contrastingly, all African samples conformed to the neutral model of populations at equilibrium and presented more diversified distributions. This suggested that part of the apparent African divergence was due to heterogeneous evolutionary processes and confirmed that some diversity reducing factors were at work in Caucasoids and Orientals. Several nonexclusive hypotheses accounting for the rejection of the neutrality tests were discussed. Alternative hypotheses concerning modern human emergence were also reviewed in the light of present results.