Substitution rates of organelle and nuclear genes in sharks: implicating metabolic rate (again)

Rates of nucleotide substitution for nuclear genes are thought to be governed primarily by the number of germ line replication events (the so-called "generation time" hypothesis). In contrast, rates of mitochondrial DNA evolution appear to be set primarily by DNA damage pathways of mutation mediated by mutagenic by-products of oxidative phosphorylation (the so-called "metabolic-rate" hypothesis). Comparison of synonymous substitution rates estimated for the mitochondrial cytochrome b gene and nuclear-encoded dlx, hsp70, and RAG-1 genes in mammals and sharks shows that rates of molecular evolution for sharks are approximately an order of magnitude slower than those for mammals for both nuclear and mitochondrial genes. In addition, there is significant positive covariation of substitution rate for mitochondrial and nuclear genes within sharks. These results, interpreted in light of the pervasiveness of DNA damage by mutagenic by-products of oxygen metabolism to both nuclear and mitochondrial genes and coupled with increasing evidence for cross-genome activity of DNA repair enzymes, suggest that molecular clocks for mitochondrial and nuclear genes may be set primarily by common mutational mechanisms.      

Relative patterns and rates of evolution in heron nuclear and mitochondrial DNA

Mitochondrial cytochrome b sequence data from 15 species of herons (Aves: Ardeidae), representing 13 genera, were compared with DNA hybridization data of single-copy nuclear DNA (scnDNA) from the same species in a taxonomic congruence assessment of heron phylogeny. The two data sets produced a partially resolved, completely congruent estimate of phylogeny with the following basic structure: (Tigrisoma, Cochlearius, (((Zebrilus, (Ixobrychus, Botaurus)), (((Ardea, Casmerodius), Bubulcus), ((Egretta thula, Egretta caerulea, Egretta tricolor), Syrigma), Butorides, Nycticorax, Nyctanassa)))). Because congruence indicated similar phylogenetic information in the two data sets, we used the relatively unsaturated DNA hybridization distances as surrogates of time to examine graphically the patterns and rates of change in cytochrome b distances. Cytochrome b distances were computed either from whole sequences or from partitioned sequences consisting of transitions, transversions, specific codon site positions, or specific protein-coding regions. These graphical comparisons indicated that unpartitioned cytochrome b has evolved at 5-10 times the rate of scnDNA. Third-position transversions appeared to offer the most useful sequence partition for phylogenetic analysis because of their relatively fast rate of substitution (two times that of scnDNA) and negligible saturation. We also examined lineage-based rates of evolution by comparing branch length patterns between the nuclear and cytochrome b trees. The degree of correlation in corresponding branch lengths between cytochrome b and DNA hybridization trees depended on DNA sequence partitioning. When cytochrome b sequences were not partitioned, branch lengths in the cytochrome b and DNA hybridization trees were not correlated. However, when cytochrome b sequences were reduced to third-position transversions (i.e., unsaturated, relatively fast changing data), branch lengths were correlated. This finding suggests that lineage-based rates of DNA evolution in nuclear and mitochondrial genomes are influenced by common causes.     

Mitochondrial versus nuclear gene sequences in deep-level mammalian phylogeny reconstruction

Both mitochondrial and nuclear gene sequences have been employed in efforts to reconstruct deep-level phylogenetic relationships. A fundamental question in molecular systematics concerns the efficacy of different types of sequences in recovering clades at different taxonomic levels. We compared the performance of four mitochondrial data sets (cytochrome b, cytochrome oxidase II, NADH dehydrogenase subunit I, 12S rRNA-tRNA-16S rRNA) and eight nuclear data sets (exonic regions of alpha-2B adrenergic receptor, aquaporin, ss-casein, gamma-fibrinogen, interphotoreceptor retinoid binding protein, kappa-casein, protamine, von Willebrand Factor) in recovering deep-level mammalian clades. We employed parsimony and minimum-evolution with a variety of distance corrections for superimposed substitutions. In 32 different pairwise comparisons between these mitochondrial and nuclear data sets, we used the maximum set of overlapping taxa. In each case, the variable-length bootstrap was used to resample at the size of the smaller data set. The nuclear exons consistently performed better than mitochondrial protein and rRNA-tRNA coding genes on a per-residue basis in recovering benchmark clades. We also concatenated nuclear genes for overlapping taxa and made comparisons with concatenated mitochondrial protein-coding genes from complete mitochondrial genomes. The variable-length bootstrap was used to score the recovery of benchmark clades as a function of the number of resampled base pairs. In every case, the nuclear concatenations were more efficient than the mitochondrial concatenations in recovering benchmark clades. Among genes included in our study, the nuclear genes were much less affected by superimposed substitutions. Nuclear genes having appropriate rates of substitution should receive strong consideration in efforts to reconstruct deep-level phylogenetic relationships.     

Estimating hybridization in the presence of coalescence using phylogenetic intraspecific sampling

Background

Recent research has indicated a positive association between rates of molecular evolution and diversification in a number of taxa. However debate continues concerning the universality and cause of this relationship. Here, we present the first systematic investigation of this relationship within the mammals. We use phylogenetically independent sister-pair comparisons to test for a relationship between substitution rates and clade size at a number of taxonomic levels. Total, non-synonymous and synonymous substitution rates were estimated from mitochondrial and nuclear DNA sequences.

Results

We found no evidence for an association between clade size and substitution rates in mammals, for either the nuclear or the mitochondrial sequences. We found significant associations between body size and substitution rates, as previously reported.

Conclusions

Our results present a contrast to previous research, which has reported significant positive associations between substitution rates and diversification for birds, angiosperms and reptiles. There are three possible reasons for the differences between the observed results in mammals versus other clades. First, there may be no link between substitution rates and diversification in mammals. Second, this link may exist, but may be much weaker in mammals than in other clades. Third, the link between substitution rates and diversification may exist in mammals, but may be confounded by other variables.     

Phylogeny and comparative substitution rates of frogs inferred from sequences of three nuclear genes

Phylogenetic relationships among major clades of anuran amphibians were studied using partial sequences of three nuclear protein coding genes, Rag-1, Rag-2, and rhodopsin in 26 frog species from 18 families. The concatenated nuclear data set comprised 2,616 nucleotides and was complemented by sequences of the mitochondrial 12S and 16S rRNA genes for analyses of evolutionary rates. Separate and combined analyses of the nuclear markers supported the monophyly of modern frogs (Neobatrachia), whereas they did not provide support for the monophyly of archaic frog lineages (Archaeobatrachia), contrary to previous studies based on mitochondrial data. The Neobatrachia contain two well supported clades that correspond to the subfamilies Ranoidea (Hyperoliidae, Mantellidae, Microhylidae, Ranidae, and Rhacophoridae) and Hyloidea (Bufonidae, Hylidae, Leptodactylidae, and Pseudidae). Two other families (Heleophrynidae and Sooglossidae) occupied basal positions and probably represent ancient relicts within the Neobatrachia, which had been less clearly indicated by previous mitochondrial analyses. Branch lengths of archaeobatrachians were consistently shorter in all separate analyses, and nonparametric rate smoothing indicated accelerated substitution rates in neobatrachians. However, relative rate tests confirmed this tendency only for mitochondrial genes. In contrast, nuclear gene sequences from our study and from an additional GenBank survey showed no clear phylogenetic trends in terms of differences in rates of molecular evolution. Maximum likelihood trees based on Rag-1 and using only one neobatrachian and one archaeobatrachian sequence, respectively, even had longer archaeobatrachian branches averaged over all pairwise comparisons. More data are necessary to understand the significance of a possibly general assignation of short branches to basal and species-poor taxa by tree-reconstruction algorithms.     

Slow mitochondrial DNA sequence evolution in the Anthozoa (Cnidaria)

Mitochondrial genes have been used extensively in population genetic and phylogeographical analyses, in part due to a high rate of nucleotide substitution in animal mitochondrial DNA (mtDNA). Nucleotide sequences of anthozoan mitochondrial genes, however, are virtually invariant among conspecifics, even at third codon positions of protein-coding sequences. Hence, mtDNA markers are of limited use for population-level studies in these organisms. Mitochondrial gene sequence divergence among anthozoan species is also low relative to that exhibited in other animals, although higher level relationships can be resolved with these markers. Substitution rates in anthozoan nuclear genes are much higher than in mitochondrial genes, whereas nuclear genes in other metazoans usually evolve more slowly than, or similar to, mitochondrial genes. Although several mechanisms accounting for a slow rate of sequence evolution have been proposed, there is not yet a definitive explanation for this observation. Slow evolution and unique characteristics may be common in primitive metazoans, suggesting that patterns of mtDNA evolution in these organisms differ from that in other animal systems.     

Exploring patterns and extent of bias in estimating divergence time from mitochondrial DNA sequence data in a particular lineage: a case study of salamanders (order Caudata)

In the practice of molecular dating, substitution saturation will bias the results if not properly modeled. Date estimates based on commonly used mitochondrial DNA sequences likely suffer from this problem because of their high substitution rate. Nevertheless, the patterns and extent of such expected bias remain unknown for many major evolutionary lineages, which often differ in ages, available calibrations, and substitution rates of their mitochondrial genome. In this case study of salamanders, we used estimates based on multiple nuclear exons to assess the effects of saturation on dating divergences using mitochondrial genome sequences on a timescale of ~200-300 My. The results indicated that, due to saturation for older divergences and in the absence of younger effective calibration points, dates derived from the mitochondrial data were considerably overestimated and systematically biased toward the calibration point for the ingroup root. The overestimate might be as great as 3-10 times (about 20 My) older than actual divergence dates for recent splitting events and 40 My older for events that are more ancient. For deep divergences, dates estimated were strongly compressed together. Furthermore, excluding the third codon positions of protein-coding genes or only using the RNA genes or second codon positions did not considerably improve the performance. In the order Caudata, slowly evolving markers such as nuclear exons are preferred for dating a phylogeny covering a relatively wide time span. Dates estimated from these markers can be used as secondary calibrations for dating recent events based on rapidly evolving markers for which mitochondrial DNA sequences are attractive candidates due to their short coalescent time. In other groups, similar evaluation should be performed to facilitate the choice of markers for molecular dating and making inferences from the results.     

Evolutionary and biogeographic patterns of the Badidae (Teleostei: Perciformes) inferred from mitochondrial and nuclear DNA sequence data

We reconstructed phylogenetic relationships of the family Badidae using both mitochondrial and nuclear nucleotide sequence data to address badid systematics and to evaluate the role of vicariant speciation on their evolution and current distribution. Phy-logenetic hypotheses were derived from complete cytochrome b (1,140 base pairs) sequences of 33 individuals representing 13 badid species, and using three species of Nandidae as outgroups. Additionally, we sequenced the nuclear RAG1 (1,473 base pairs) and Tmo-4C4 (511 base pairs) genes from each of the badid species and one representative of the outgroup. Our molecular data provide the first phylogenetic hypothesis of badid intrarelationships. Analysis of the mitochondrial and nuclear nucleotide sequence data sets resulted in well-supported trees, indicating a basal split between the genera Dario and Badis, and further supporting the division of the genus Badis into five species groups as suggested by a previous taxonomic revision of the Badidae. Within the genus Badis, mitochondrial and nuclear phylogenies differed in the relative position of B. kyar. We also used our molecular phylogeny to test a vicariant speciation hypothesis derived from geological evidence of large-scale changes in drainage patterns in the Miocene affecting the Irrawaddy- and Tsangpo-Brahmaputra drainages, in the southeastern Himalaya. Within both genera, Badis and Dario, we observed a divergence into Irrawaddy- and Tsangpo-Brahmaputra clades. Using a cytb substitution rate of 8.2 x 10(-9) (substitutions x base pair(-1) x year(-1), we tentatively date this vicariant event at the Oligocene-Miocene boundary (19-24Myr). It is concordant with a hypothesized paleo connection of the Tsangpo river with the Irrawaddy drainage that was most likely interrupted during Miocene orogenic events through tectonic uplifts in eastern Tibet. Our data, therefore, indicate a substantial role of vicariant-based speciation shaping the current distribution patterns of badids.     

Primate numts and reticulate evolution of capped and golden leaf monkeys (Primates: Colobinae)

A recent phylogenetic study of langurs and leaf monkeys of South Asia suggested a reticulate evolution of capped and golden leaf monkeys through ancient hybridization between Semnopithecus and Trachypithecus. To test this hybridization scenario, I analysed nuclear copies of the mitochondrial cytochrome b gene (numts) from capped, golden and Phayre’s leaf monkeys. These numts were aligned with mitochondrial cytochrome b sequences of various species belonging to the genera Semnopithecus and Trachypithecus. In the phylogenetic tree derived from this alignment, the numts fell into three distinct clades (A, B and C) suggesting three independent integration events. Clade A was basal to Semnopithecus, and clades B and C were basal to Trachypithecus. Among the numts in clades A and C were sequences derived from species not represented in their respective sister mitochondrial groups. This unusual placement of certain numts is taken as additional support for the hybridization scenario. Based on the molecular dating of these integration events, hybridization is estimated to have occurred around 7.1 to 3.4 million years ago. Capped and golden leaf monkeys might have to be assigned to a new genus to reconcile their unique evolutionary history. Additionally, northeast India appears to be a ‘hot spot’ for lineages that might have evolved through reticulate evolution.     

Sequence evolution and phylogenetic signal in control-region and cytochrome b sequences of rainbow fishes (Melanotaeniidae)

The nucleotide sequences of segments of the cytochrome b gene (351 bp), the tRNA(Pro) gene (49 bp), and the control region (approximately 313 bp) of mitochondrial DNA were obtained from 26 fish representing different populations and species of Melanotaenia and one species of Glossolepis, freshwater rainbow fishes confined to Australia and New Guinea. The purpose was to investigate relative rates and patterns of sequence evolution. Overall levels of divergence were similar for the cytochrome b and tRNA control-region sequences, both ranging from < 1% within subspecies to 15%-19% between genera. However, the patterns of sequence evolution differed. For the cytochrome b gene, transitions consistently exceeded transversions, the bias ranging from 4.2:1 to 2:1, depending on the level of sequence divergence. However, in the control-region sequence, a bias toward transitions (2:1) was observed only in comparisons between very similar sequences, and transversions outnumbered transitions in comparisons of divergent sequences. Graphic comparisons suggested that the control region was saturated for transitions at relatively low levels of sequence divergence but accumulated transversions at a greater rate than did the cytochrome b sequence. These distinct patterns of base substitution are associated with differences in A+T content, which is 70% for the tRNA control- region segment versus 50% for cytochrome b. A test for skewness in the distribution of lengths of random trees indicated that both segments contained phylogenetic signal. Parsimony analyses of the data from the two regions, with or without weighting schemes appropriate to the respective patterns of sequence evolution, identified the same five groupings of sequences, but the relationships among the groups differed. However, in most cases the branches uniting different combinations of groups were poorly supported, and the differences among topologies were insignificant. Considering the observed patterns of base substitution and the results of the phylogenetic analyses, we deduce that both the control region and cytochrome b are appropriate for population genetic studies but that the control region is less effective than cytochrome b for resolving relationships among divergent lineages of rainbow fishes.      

Mitochondrial pseudogenes are pervasive and often insidious in the snapping shrimp genus Alpheus

Here we show that multiple DNA sequences, similar to the mitochondrial cytochrome oxidase I (COI) gene, occur within single individuals in at least 10 species of the snapping shrimp genus Alpheus. Cloning of amplified products revealed the presence of copies that differed in length and (more frequently) in base substitutions. Although multiple copies were amplified in individual shrimp from total genomic DNA (gDNA), only one sequence was amplified from cDNA. These results are best explained by the presence of nonfunctional duplications of a portion of the mtDNA, probably located in the nuclear genome, since transfer into the nuclear gene would render the COI gene nonfunctional due to differences in the nuclear and mitochondrial genetic codes. Analysis of codon variation suggests that there have been 21 independent transfer events in the 10 species examined. Within a single animal, differences between the sequences of these pseudogenes ranged from 0.2% to 20.6%, and those between the real mtDNA and pseudogene sequences ranged from 0.2% to 18.8% (uncorrected). The large number of integration events and the large range of divergences between pseudogenes and mtDNA sequences suggest that genetic material has been repeatedly transferred from the mtDNA to the nuclear genome of snapping shrimp. Unrecognized pseudogenes in phylogenetic or population studies may result in spurious results, although previous estimates of rates of molecular evolution based on Alpheus sister taxa separated by the Isthmus of Panama appear to remain valid. Especially worrisome for researchers are those pseudogenes that are not obviously recognizable as such. An effective solution may be to amplify transcribed copies of protein-coding mitochondrial genes from cDNA rather than using genomic DNA.     

Combined mitochondrial and nuclear sequences support the monophyly of forcipulatacean sea stars

Previous molecular phylogenetic analyses of forcipulatacean sea stars (Echinodermata: Asteroidea) have reconstructed a non-monophyletic order Forcipulatida, provided that two or more forcipulate families are included. This result could mean that one or more assumptions of the reconstruction method was violated, or else the traditional classification could be erroneous. The present molecular phylogenetic analysis included 12 non-forcipulatacean and 39 forcipulatacean sea stars, with multiple representatives of all but one of the forcipulate families and/or subfamilies. Bayesian analysis of approximately 4.2kb of sequence data representing seven partitions (nuclear 18S rRNA and 28S rRNA, mitochondrial 12S rRNA, 16S rRNA, 5 tRNAs and cytochrome oxidase I with first and second codon positions analyzed separately from third codon positions) recovered a consensus tree with three well-supported clades (78%-100% bootstrap support) that corresponded at least approximately to traditional taxonomic ranks: the superorder Forcipulatacea (Forcipulatida + Brisingida) + Pteraster, the Brisingida/Brisingidae and Asteriidae + Rathbunaster + Pycnopodia. When a molecular clock was enforced, the partitioned Bayesian analysis recovered the traditional Forcipulatacea. Five of six genera represented by two or more species were monophyletic with 100% bootstrap support. Most of the traditional subfamilial and familial groupings within the Forcipulatida were either unresolved or non-monophyletic. The separate partitions differed considerably in estimates of model parameters, mainly between nuclear sequences (with high GC content, low rates of sequence substitution and high transition/transversion rate ratios) and mitochondrial sequences.     

Decoupling of molecular and morphological evolution in deep lineages of a meiobenthic harpacticoid copepod

Molecular and biochemical genetic analyses have revealed that many marine invertebrate taxa, including some well-studied and presumably cosmopolitan species, are actually complexes of sibling species. When morphological differences are slight and estimated divergence times are old, data suggest either unusually high rates of sequence evolution or long-term morphological stasis. Here, five gene regions (mitochondrial cytochrome oxidase subunit I and large-subunit ribosomal 16S rDNA and nuclear ITS1, 5.8S rDNA, and ITS2) were analyzed in four geographic samples of the meiobenthic harpacticoid copepod Cletocamptus deitersi. Molecular sequences revealed four extremely differentiated molecular lineages with unalignable nuclear intergenic spacers and mitochondrial uncorrected divergences reaching 25% (cytochrome oxidase) and 36% (16S rDNA). These levels of divergence are greater than those reported previously for congeneric species in diverse invertebrate taxa, including crustaceans. The nominally intraspecific divergence matches or exceeds the corresponding divergence from a known congener (Cletocamptus helobius). A molecular clock applied to the cytochrome oxidase subunit I data suggests that these lineages split in the Miocene, consistent with the fossil record of a North American Cletocamptus from the same period. Morphological differences among the major lineages are subtle but congruent with the patterns of genetic differentiation. Our conclusion, based on concordant patterns of variation in two mitochondrial and three nuclear gene regions, as well as morphological observations, is that C. deitersi in North America is composed of at least four separate species by the genealogical concordance, phylogenetic, and morphological-species criteria. Alternative explanations for the deep phylogenetic nodes and apparent morphological stasis, including high rates of sequence evolution, balancing selection, and genetic signatures of historical events, are considered unlikely.     

The problems of molecular phylogenetics with the example of squamate reptiles: Mitochondrial DNA markers

The review considers the current problems of molecular phylogenetics based on mitochondrial and chromosomal DNA sequences. The emphasis is placed on mtDNA markers, which are widely employed in reconstructing molecular evolution, but often without a critical analysis of the physiological and biochemical features of mitochondria that affect the adequacy and reliability of the results. In addition to the factors that make mtDNA-based phylogenies difficult to interpret (unrecognized hybridization and introgression events, ancestral polymorphism, and nuclear paralogs of mtDNA sequences), attention is paid to the nonneutrality and unequal mutation rates of mtDNA genes and their fragments, violations of uniparental inheritance of mitochondria, recombination events, natural heteroplasmy, and mtDNA haplotypic diversity. These factors may influence the congruence of phylogenetic inferences and trees constructed for the same organisms with different mtDNA markers or with mitochondrial and nuclear markers. The review supports the viewpoint that mitochondrial genes and their fragments fail to provide reliable evolutionary markers when considered without a thorough study of the environmental conditions and life of the taxa. The influence of external conditions on the metabolism and physiology of mitochondria cannot be taken into account in full nor modeled well enough for phylogenetic applications. It is assumed that mtDNA is valuable as a phylogenetic marker primarily because its complete sequence may be analyzed to identify the apomorphic and synmorphic properties of a taxon and to search for informative nuclear paralogs of mtDNA for phylogeographical studies and estimations of relative evolution times.     

Patterns of nucleotide substitution among simultaneously duplicated gene pairs in Arabidopsis thaliana

We characterized rates and patterns of synonymous and nonsynonymous substitution in 242 duplicated gene pairs on chromosomes 2 and 4 of Arabidopsis thaliana. Based on their collinear order along the two chromosomes, the gene pairs were likely duplicated contemporaneously, and therefore comparison of genetic distances among gene pairs provides insights into the distribution of nucleotide substitution rates among plant nuclear genes. Rates of synonymous substitution varied 13.8-fold among the duplicated gene pairs, but 90% of gene pairs differed by less than 2.6-fold. Average nonsynonymous rates were approximately fivefold lower than average synonymous rates; this rate difference is lower than that of previously studied nonplant lineages. The coefficient of variation of rates among genes was 0.65 for nonsynonymous rates and 0.44 for synonymous rates, indicating that synonymous and nonsynonymous rates vary among genes to roughly the same extent. The causes underlying rate variation were explored. Our analyses tentatively suggest an effect of physical location on synonymous substitution rates but no similar effect on nonsynonymous rates. Nonsynonymous substitution rates were negatively correlated with GC content at synonymous third codon positions, and synonymous substitution rates were negatively correlated with codon bias, as observed in other systems. Finally, the 242 gene pairs permitted investigation of the processes underlying divergence between paralogs. We found no evidence of positive selection, little evidence that paralogs evolve at different rates, and no evidence of differential codon usage or third position GC content.     

Comparative evolution of the mitochondrial cytochrome b gene and nuclear beta-fibrinogen intron 7 in woodpeckers

Most molecular phylogenetic studies of vertebrates have been based on DNA sequences of mitochondrial-encoded genes. MtDNA evolves rapidly and is thus particularly useful for resolving relationships among recently evolved groups. However, it has the disadvantage that all of the mitochondrial genes are inherited as a single linkage group so that only one independent gene tree can be inferred regardless of the number of genes sequenced. Introns of nuclear genes are attractive candidates for independent sources of rapidly evolving DNA: they are pervasive, most of their nucleotides appear to be unconstrained by selection, and PCR primers can be designed for sequences in adjacent exons where nucleotide sequences are conserved. We sequenced intron 7 of the beta-fibrinogen gene (beta-fibint7) for a diversity of woodpeckers and compared the phylogenetic signal and nucleotide substitution properties of this DNA sequence with that of mitochondrial-encoded cytochrome b (cyt b) from a previous study. A few indels (insertions and deletions) were found in the beta-fibint7 sequences, but alignment was not difficult, and the indels were phylogentically informative. The beta-fibint7 and cyt b gene trees were nearly identical to each other but differed in significant ways from the traditional woodpecker classification. Cyt b evolves 2.8 times as fast as beta-fibint7 (14. 0 times as fast at third codon positions). Despite its relatively slow substitution rate, the phylogenetic signal in beta-fibint7 is comparable to that in cyt b for woodpeckers, because beta-fibint7 has less base composition bias and more uniform nucleotide substitution probabilities. As a consequence, compared with cyt b, beta-fibint7 nucleotide sites are expected to enter more distinct character states over the course of evolution and have fewer multiple substitutions and lower levels of homoplasy. Moreover, in contrast to cyt b, in which nearly two thirds of nucleotide sites rarely vary among closely related taxa, virtually all beta-fibint7 nucleotide sites appear free of selective constraints, which increases informative sites per unit sequenced. However, the estimated gamma distribution used to model rate variation among sites suggests constraints on some beta-fibint7 sites. This study suggests that introns will be useful for phylogenetic studies of recently evolved groups.     

Evolutionary rates of insertion and deletion in noncoding nucleotide sequences of primates

Insertions and deletions are responsible for gaps in aligned nucleotide sequences, but they have been usually ignored when the number of nucleotide substitutions was estimated. We compared six sets of nuclear and mitochondrial noncoding DNA sequences of primates and obtained the estimates of the evolutionary rate of insertion and deletion. The maximum-parsimony principle was applied to locate insertions and deletions on a given phylogenetic tree. Deletions were about twice as frequent as insertions for nuclear DNA, and single-nucleotide insertions and deletions were the most frequent in all events. The rate of insertion and deletion was found to be rather constant among branches of the phylogenetic tree, and the rate (approximately 2.0/kb/Myr) for mitochondrial DNA was found to be much higher than that (approximately 0.2/kb/Myr) for nuclear DNA. The rates of nucleotide substitution were about 10 times higher than the rate of insertion and deletion for both nuclear and mitochondrial DNA.      

Species status and phylogeography of two closely related Coptolabrus species (Coleoptera: Carabidae) in South Korea inferred from mitochondrial and nuclear gene sequences

We investigated the species status and intraspecific phylogeography in South Korea of two ground beetle species, Coptolabrus jankowskii and Coptolabrus smaragdinus (Coleoptera: Carabidae), using statistical parsimony networks and nested clade analyses based on sequences from the mitochondrial cytochrome oxidase subunit I (COI) and nuclear phosphoenolpyruvate carboxykinase (PepCK) and wingless (Wg) genes. Although traditional parsimony tree construction generally failed to resolve interspecific relationships and construct biologically meaningful genealogies, analysis using statistical parsimony networks yielded statistically significant inter- and intraspecific genealogical structures. We found that although these two species represent a notable case of trans-species polymorphisms in both mitochondrial and nuclear gene sequences, their status as separate species was evidenced by the nonrandom association between species and nested clades at various nesting levels. The exceptional occurrence of shared identical or very similar COI sequences was considered to be the result of introgressive hybridization. In addition, range expansion and fragmentation events across the Korean Peninsula and adjacent islands were inferred from nested clade phylogeographical analyses. The COI gene revealed the geographical divergence of major eastern and western clades and historical biogeographical events within each major clade, whereas the nuclear PepCK gene, which did not reveal corresponding east-west clades, indicated past fragmentation and range expansion across wide areas that may have been the result of older biogeographical events. Thus, phylogeographical inferences drawn from analyses of mitochondrial and nuclear genes can reveal different and potentially complementary information about phylogeographical processes.