How should gaps be treated in parsimony? A comparison of approaches using simulation

Simulation with indels was used to produce alignments where true site homologies in DNA sequences were known; the gaps from these datasets were removed and the sequences were then aligned to produce hypothesized alignments. Both alignments were then analyzed under three widely used methods of treating gaps during tree reconstruction under the maximum parsimony principle. With the true alignments, for many cases (82%), there was no difference in topological accuracy for the different methods of gap coding. However, in cases where a difference was present, coding gaps as a fifth state character or as separate presence/absence characters outperformed treating gaps as unknown/missing data nearly 90% of the time. For the hypothesized alignments, on average, all gap treatment approaches performed equally well. Data sets with higher sequence divergence and more pectinate tree shapes with variable branch lengths are more affected by gap coding than datasets associated with shallower non-pectinate tree shapes.     

Islands within an island: phylogeography and conservation genetics of the endangered Hawaiian tree snail Achatinella mustelina

Mitochondrial DNA (mtDNA) sequences were used to evaluate phylogeographic structure within and among populations of three endangered Hawaiian tree snail species (n = 86). The primary focus of this investigation was on setting conservation priorities for Achatinella mustelina. Limited data sets for two additional endangered Hawaiian tree snails, A. livida and A. sowerbyana, were also developed for comparative purposes. Pairwise genetic distance matrices and phylogenetic trees were generated, and an analysis of molecular variance was performed on 675-base pair cytochrome oxidase I gene sequences from multiple populations of Hawaiian tree snails. Sequence data were analysed under distance-based maximum-likelihood, and maximum-parsimony optimality criteria. Within the focal species, A. mustelina, numbers of variable and parsimony informative sites were 90 and 69, respectively. Pairwise intraspecific mtDNA sequence divergence ranged from 0 to 5.3% in A. mustelina, from 0 to 1.0% in A. livida and from 0 to 1.9% in A. sowerbyana. For A. mustelina, population genetic structure and mountain topography were strongly correlated. Maximum genetic distances were observed across deep, largely deforested valleys, and steep mountain peaks, independent of geographical distance. However, in certain areas where forest cover is presently fragmented, little mtDNA sequence divergence exists despite large geographical scales (8 km). Genetic data were used to define evolutionarily significant units for conservation purposes including decisions regarding placement of predator exclusion fences, captive propagation, re-introduction and translocation.     

Mitogenomic analysis for coelacanths (Latimeria chalumnae) caught in Tanzania

In recent years, a large number of individuals of the species Latimeria chalumnae, one of the living fossil coelacanths, have been landed off the coast of Tanzania. Although L. chalumnae specimens have also been landed at other localities in the western Indian Ocean, so far, viable populations of this species have been identified only at two localities, Comoros and South Africa. Therefore, the recent active catch off Tanzania suggests a new habitat for L. chalumnae. To examine the genetic background of the Tanzanian fish, we analyzed complete mtDNA sequences of two Tanzanian individuals (Kigombe-9 and Songo Mnara-1) collected from the north and south coasts of Tanzania. Using the recently reported criteria for six haplotypes established in a population genetic study for coelacanths living in the western Indian Ocean [Schartl, M., Hornung, U., Hissman, K., Schauer, J., Fricke, H., 2005. Relatedness among east African coelacanths. Nature 435, 901.], we characterized Songo Mnara-1 as haplotype 1 and Kigombe-9 as haplotype 5. We suggest that the Songo Mnara specimen is a member of the Comoran group, but was swept away by the South Equatorial current. The individual from Kigombe may be a member of an undiscovered population that exists near the boundary between Tanzania and Kenya. Further analysis using more than 19 individuals recently captured off the north coast of Tanzania will reveal whether a new population exists there. Our sequence data suggest additional variable sites in the mtDNA sequence that may define the population structure of coelacanths in the western Indian Ocean and also raise the possibility that the previously published Comoran coelacanth mtDNA sequence contains several critical errors including base changes and indels.     

Molecular evolution of insertions and deletion in the chloroplast genome of silene

Insertions, deletions, and inversions in the chloroplast genome of higher plants have been shown to be extremely useful for resolving phylogenetic relationships both between closely related taxa and among more basal lineages. Introns and intergenic spacers from the chloroplast genome are now increasingly used for phylogenetic and population genetic studies of populations from a single species, and it is therefore interesting to know whether indels can provide useful data and hence increase the power of intraspecific studies. Here, we show that indels in three cpDNA intergenic spacers and one cpDNA intron for two species of Silene evolve at slightly higher rates than base pair substitutions. Repeat indels appear to have the highest rate of evolution and are thus more prone to homoplasy. We show that coded indel data have high information content for phylogenetic analysis, and indels thus provide useful information to infer phylogenetic relationships at the intraspecific level.     

Evolution and phylogenetic utility of alignment gaps within intron sequences of three nuclear genes in bumble bees (Bombus)

To test whether gaps resulting from sequence alignment contain phylogenetic signal concordant with those of base substitutions, we analyzed the occurrence of indel mutations upon a well-resolved, substitution-based tree for three nuclear genes in bumble bees (Bombus, Apidae: Bombini). The regions analyzed were exon and intron sequences of long-wavelength rhodopsin (LW Rh), arginine kinase (ArgK), and elongation factor-1alpha (EF-1alpha) F2 copy genes. LW Rh intron had only a few uninformative gaps, ArgK intron had relatively long gaps that were easily aligned, and EF-1alpha intron had many short gaps, resulting in multiple optimal alignments. The unambiguously aligned gaps within ArgK intron sequences showed no homoplasy upon the substitution-based tree, and phylogenetic signals within ambiguously aligned regions of EF-1alpha intron were highly congruent with those of base substitutions. We further analyzed the contribution of gap characters to phylogenetic reconstruction by incorporating them in parsimony analysis. Inclusion of gap characters consistently improved support for nodes recovered by substitutions, and inclusion of ambiguously aligned regions of EF-1alpha intron resolved several additional nodes, most of which were apical on the phylogeny. We conclude that gaps are an exceptionally reliable source of phylogenetic information that can be used to corroborate and refine phylogenies hypothesized by base substitutions, at least at lower taxonomic levels. At present, full use of gaps in phylogenetic reconstruction is best achieved in parsimony analysis, pending development of well-justified and generally applicable methods for incorporating indels in explicitly model-based methods.     

Genetic structure in Iberian and Moroccan populations of the globally threatened great bustard Otis tarda: a microsatellite perspective

Patterns of genetic structure and gene flow among populations help us understand population dynamics and properly manage species of concern. Matrilineal mtDNA sequence data have been instrumental in revealing genetic structure at the intraspecific level, but bi‐parentally inherited markers are needed to confirm patterns at the genome level and to assess the potential role of sex‐biased dispersal on gene flow, particularly in species where males are known to be the main dispersing sex. Here we use microsatellite loci to examine patterns of genetic structure across the range of the great bustard in Iberia and Morocco, an area representing 70% of the world population of this globally threatened species. We used population differentiation statistics and Bayesian analysis of population structure to analyse data from 14 microsatellite loci. These data provide greater resolution than mtDNA sequences, and results reveal the existence of three main genetic units corresponding to Morocco, the northeastern part of Spain, and the rest of the Iberian Peninsula. Our results, together with previous mtDNA data, confirm the genetic differentiation of the northern Africa population and the importance of the Strait of Gibraltar as a barrier to gene flow for both males and females, rendering the Moroccan population a separate management unit of high conservation concern.     

Genetic homogeneity among colonies of the white-nest swiftlet (Aerodramus fuciphagus) in Thailand

The white-nest swiftlet, Aerodramus fuciphagus, originally lived in large colonies in natural caves, but now it also occurs in man-made buildings. We investigated the patterns of genetic differentiation in two mitochondrial DNA genes (cyt-b and ND2) and eight microsatellite loci among and within colonies of A. fuciphagus from across recently established man-made colonies in Thailand. Ten white-nest swiftlet colonies were sampled along the coast of the Gulf of Thailand and the Andaman Sea in Thailand during 2003-2006. The genetic diversity of mtDNA was very low, and few significant PhiST values were found between pairs of colonies. Analyses of haplotype relationships did not show genetic structure across the sampled distribution. The level of genetic diversity for microsatellite loci was high, but FST values were not significant. However, due to small sample sizes for some colonies that could limit conclusions on genetic differentiation from PhiST and FST, we also analyzed the microsatellite data using STRUCTURE and found that number of subpopulations of white-nest swiftlets in sampled colonies was one. The lack of genetic differentiation among swiftlet house colonies could be a result of high gene flow between colonies and large population sizes. Our results suggest that A. fuciphagus living in recently established man-made colonies in Thailand should be considered members of a single panmictic population. Future work will be necessary to determine whether this panmixia is stable or a temporary result of the recent explosive expansion of the number of colonies, and comparisons to natural colonies may provide an understanding of mechanisms producing the lack of genetic structure in swiftlet house colonies.     

Highly conserved nuclear copies of the mitochondrial control region in the desert locust Schistocerca gregaria: some implications for population studies

Animal mitochondrial DNA has proved a valuable marker in intraspecific systematic studies. However, if nucleotide sequence heterogeneity exists at the individual level, its usefulness will be much reduced. This study demonstrates that the presence of highly conserved non-coding mitochondrial sequences in the nuclear genome of Schistocerca gregaria greatly impairs the use of mtDNA in population genetic studies. Caution is called for in other organisms; and it seems necessary to check for conserved nuclear copies of mitochondrial sequences before launching into a large scale analysis of populations using mtDNA as a genetic marker. Experimental procedures are suggested for this purpose.     

Mitochondrial DNA phylogeography of Acrossocheilus paradoxus (Cyprinidae) in Taiwan

Nucleotide sequences of 3' end of the cytochrome b gene, tRNA genes, D-loop control region, and the 5' end of the 12S rRNA of mitochondrial DNA (mtDNA) were used to assess the genetic and phylogeographic structure of Acrossocheilus paradoxus populations, a Cyprinidae fish of Taiwan. A hierarchical examination of populations in 12 major streams from three geographical regions using an analysis of molecular variance (AMOVA) indicates high genetic differentiation both among populations (PhiST = 0.511, P < 0.001) and among regions (PhiCT = 0.368, P < 0.001). Limited migration largely contributed to the genetic differentiation. High nucleotide diversity (1.13%) and haplotype diversity (0.80%) were detected among populations. The degree of genetic differentiation was correlated with geographical distance between populations, a result consistent with the one-dimensional stepping stone models. A neighbour-joining tree recovered by (DAMBE) supports the pattern of isolation by distance and reveals a closer relationship between populations of the central and southern regions. A minimum spanning network based on nucleotide substitutions reflected migration routes from populations of the central region to the northern and southern regions, respectively. Postglacial colonization and expansion can explain the phylogeographical pattern. Single and ancient migration events may have allowed the northern region to attain the monophyly of mtDNA alleles. In contrast, most populations within geographical regions are either paraphyletic or polyphyletic due to the relatively shorter time period for coalescence. Both low haplotype number and genetic variability suggest a bottleneck event in the Chingmei population of northern Taiwan. Based on coalescence theory, the monophyly of the Tungkang population of the southern region may be associated with a founder event.     

Phylogenetic Relationships Among Fragmented Asian Black Bear (Ursus Thibetanus) Populations in Western Japan

The number of Asian black bears (Ursus thibetanus) in Japan has been reduced and their habitats fragmented and isolated because of human activities. Our previous study examining microsatellite DNA loci revealed significant genetic differentiation among four local populations in the western part of Honshu. Here, an approximate 700-bp nucleotide sequence of mitochondrial DNA (mtDNA) control region was analysed in 119 bears to infer the evolutionary history of these populations. Thirteen variable sites and variation in the number of Ts at a T-repeat site were observed among the analysed sequences, which defined 20 mtDNA haplotypes with the average sequence divergence of 0.0051 (SD = 0.00001). The observed haplotype frequencies differed significantly among the four populations. Phylogeographic analysis of the haplotypes suggested that black bears in this region have gone through two different colonisation histories, since the observed haplotypes belonged to two major monophyletic lineages and the lineages were distributed with an apparent border. The spatial genetic structure revealed by using mtDNA was different from that observed using microsatellite DNA markers, probably due to female philopatry and male-biased dispersal. Since nuclear genetic diversity will be lost in the three western populations because of the small population size and genetic isolation, their habitats need to be preserved, and these four populations should be linked to each other by corridors to promote gene flow from the easternmost population with higher nuclear genetic diversity.     

mtDNA diversity in rhesus monkeys reveals overestimates of divergence time and paraphyly with neighboring species

Reconstructions of the human-African great ape phylogeny by using mitochondrial DNA (mtDNA) have been subject to considerable debate. One confounding factor may be the lack of data on intraspecific variation. To test this hypothesis, we examined the effect of intraspecific mtDNA diversity on the phylogenetic reconstruction of another Plio- Pleistocene radiation of higher primates, the fascicularis group of macaque (Macaca) monkey species. Fifteen endonucleases were used to identify 10 haplotypes of 40-47 restriction sites in M. mulatta, which were compared with similar data for the other members of this species group. Interpopulational, intraspecific mtDNA diversity was large (0.5%- 4.5%), and estimates of divergence time and branching order incorporating this variation were substantially different from those based on single representatives of each species. We conclude that intraspecific mtDNA diversity is substantial in at least some primate species. Consequently, without prior information on the extent of genetic diversity within a particular species, intraspecific variation must be assessed and accounted for when reconstructing primate phylogenies. Further, we question the reliability of hominoid mtDNA phylogenies, based as they are on one or a few representatives of each species, in an already depauperate superfamily of primates.      

Phylogenetic network for European mtDNA

The sequence in the first hypervariable segment (HVS-I) of the control region has been used as a source of evolutionary information in most phylogenetic analyses of mtDNA. Population genetic inference would benefit from a better understanding of the variation in the mtDNA coding region, but, thus far, complete mtDNA sequences have been rare. We determined the nucleotide sequence in the coding region of mtDNA from 121 Finns, by conformation-sensitive gel electrophoresis and subsequent sequencing and by direct sequencing of the D loop. Furthermore, 71 sequences from our previous reports were included, so that the samples represented all the mtDNA haplogroups present in the Finnish population. We found a total of 297 variable sites in the coding region, which allowed the compilation of unambiguous phylogenetic networks. The D loop harbored 104 variable sites, and, in most cases, these could be localized within the coding-region networks, without discrepancies. Interestingly, many homoplasies were detected in the coding region. Nucleotide variation in the rRNA and tRNA genes was 6%, and that in the third nucleotide positions of structural genes amounted to 22% of that in the HVS-I. The complete networks enabled the relationships between the mtDNA haplogroups to be analyzed. Phylogenetic networks based on the entire coding-region sequence in mtDNA provide a rich source for further population genetic studies, and complete sequences make it easier to differentiate between disease-causing mutations and rare polymorphisms.     

Closing the gap: Avian lineage splits at a young,narrow seaway imply a protracted history of mixed population response

The evolutionary significance of spatial habitat gaps has been well recognized since Alfred Russel Wallace compared the faunas of Bali and Lombok. Gaps between islands influence population structuring of some species, and flightless birds are expected to show strong partitioning even where habitat gaps are narrow. We examined the population structure of the most numerous living flightless land bird in New Zealand, Weka (Gallirallus australis). We surveyed Weka and their feather lice in native and introduced populations using genetic data gathered from DNA sequences of mitochondrial genes and nuclear β‐fibrinogen and five microsatellite loci. We found low genetic diversity among extant Weka population samples. Two genetic clusters were evident in the mtDNA from Weka and their lice, but partitioning at nuclear loci was less abrupt. Many formerly recognized subspecies/species were not supported; instead, we infer one subspecies for each of the two main New Zealand islands. Although currently range restricted, North Island Weka have higher mtDNA diversity than the more wide‐ranging southern Weka. Mismatch and neutrality statistics indicate North Island Weka experienced rapid and recent population reduction, while South Island Weka display the signature of recent expansion. Similar haplotype data from a widespread flying relative of Weka and other New Zealand birds revealed instances of North Island—South Island partitioning associated with a narrow habitat gap (Cook Strait). However, contrasting patterns indicate priority effects and other ecological factors have a strong influence on spatial exchange at this scale.     

Effects of Asymmetric Nuclear Introgression,Introgressive Mitochondrial Sweep,and Purifying Selection on Phylogenetic Reconstruction and Divergence Estimates in the Pacific Clade of Locustella Warblers

When isolated but reproductively compatible populations expand geographically and meet, simulations predict asymmetric introgression of neutral loci from a local to invading taxon. Genetic introgression may affect phylogenetic reconstruction by obscuring topology and divergence estimates. We combined phylogenetic analysis of sequences from one mtDNA and 12 nuDNA loci with analysis of gene flow among 5 species of Pacific Locustella warblers to test for presence of genetic introgression and its effects on tree topology and divergence estimates. Our data showed that nuDNA introgression was substantial and asymmetrical among all members of superspecies groups whereas mtDNA showed no introgression except a single species pair where the invader''s mtDNA was swept by mtDNA of the local species. This introgressive sweep of mtDNA had the opposite direction of the nuDNA introgression and resulted in the paraphyly of the local species'' mtDNA haplotypes with respect to those of the invader. Тhe multilocus nuDNA species tree resolved all inter- and intraspecific relationships despite substantial introgression. However, the node ages on the species tree may be underestimated as suggested by the differences in node age estimates based on non-introgressing mtDNA and introgressing nuDNA. In turn, the introgressive sweep and strong purifying selection appear to elongate internal branches in the mtDNA gene tree.     

Toward a more accurate time scale for the human mitochondrial DNA tree

Several estimates of the time of occurrence of the most recent common mitochondrial DNA (mtDNA) ancestor of modern humans have been made. Estimates derived from noncoding regions based on a model that classifies sites into two categories (variable and invariable) have been consistently older than those derived from the third positions of codons. This discrepancy can be attributed to a violation of the assumption of rate homogeneity among variable sites when analyzing the noncoding regions. Additional data from the partial control region sequences allow us to take into account some of this further heterogeneity. By assigning the sites to three classes (highly variable, moderately variable, and invariable) and by assuming that the last common mtDNA ancestor of humans and chimpanzees lived 4 million years ago, the most recent common mtDNA ancestor of humans is estimated to have occurred 211,000 ±111,000 years ago (±1 SE), consistent with the estimate, 101,000 ± 52,000 years, made from third positions of codons and also with those proposed previously. We used the same technique to estimate when a putative expansion of modern humans out of Africa took place and estimated a time of 89,000 ± 69,000 years ago. Even though the standard errors of these estimates are large, they allow us to reject the multiregional hypothesis of modern human origin.Deceased July 21, 1991 Correspondence to: M. Hasegawa     

Comparative sequencing of a microsatellite locus reveals size homoplasy within and between european oak species (<Emphasis Type="Italic">Quercus</Emphasis> spp.)

Microsatellites (simple sequence repeats [SSRs]) are highly variable molecular markers that are a rich and readily assayed source of variation for population genetic studies. Cross-amplification between closely related species is possible when there are no (or few) sequence differences in the primer binding sites. The occurrence of nonhomologous fragments of the same size (size homoplasy) is a contraint of microsatellites. Size homoplasy can be caused by insertions/deletions (indels) in SSR flanking regions. We found that size variation in locus ssrQZAG9 is due to different repeat numbers of the SSR motifs but also to indels in SSR flanking regions. Indels were found within species belonging to sectionsRobur andCerris of genusQuercus and also between species of the 2 sections. In sectionRobur (Quercis robur L.,Quercus petraea [Matt.] Liebl.,Quercus pubescens Willd.), we detected rare alleles with an indel of 57 bp or 62 bp followed by a smaller indel of 12 bp in the SSR flanking regions. These alleles show a size range overlapping with that of alleles amplified inQuercus cerris L. (sectionCerris). Multiple alignments with sequences of sectionRobur revealed the same SSR repeat motif but multiple indels in SSR flanking regions inQ. cerris. We discuss the effects of size homoplasy of SSR loci for the study of interspecific gene flow and on estimates of population differentiation.     

Size, frequency, and phylogenetic signal of multiple-residue indels in sequence alignment of introns

Indels in DNA sequences frequently affect more than a single nucleotide, creating problems for alignment, character coding and phylogenetic analysis. However, the size and frequency of multiple‐residue indels is not usually tested, and with popular alignment packages their reconstruction is indirectly acheived by reducing the affine (gap extension) cost. We explored the length distribution of indels in intron sequences of the gene Mp20 by modifying the gap opening and gap extension costs. Given a “known” tree for the study group, global homology levels were greatest under low gap cost, with gap extension costs of roughly 0.4‐fold the opening cost. Different approaches to gap coding and weighting suggested that taxonomic congruence was correlated with high frequencies of multiple‐position indels, with a maximum indel length of 2–5 bp and few indels above 15 bp, but also including a proportion of indels > 100 bp. Only a small minority of indels could be reconstructed as single‐position indels. Consequently, tree topologies improved when homologous multinucleotide indels were recoded as binary characters which are otherwise highly homoplastic and weighted characters in single‐position coding. In tree‐generating alignment procedures as implemented in POY, where gap penalty determines the character weight during tree search, the problem of assigning inappropriately high weight to multiple‐residue indels could partly be overcome by setting the extension costs to about 0.4‐fold lower than gap opening costs. We conclude that multiple consecutive gap positions are not independent characters and hence methods for parsimony reconstruction of long indels are required. Finally, we also observed a general lack of correlation between taxonomic and character congruence, demonstrating the difficulties of applying congruence criteria to decide among competing alignments. This highlights the value of recent model‐based alignment procedures which can implement the statistical distributions of indel size classes, and do not rely on potentially circular strategies for optimizing overall congruence. © The Willi Hennig Society 2006.     

Complete rRNA gene sequences reveal that the microsporidium Nosema bombi infects diverse bumblebee (Bombus spp.) hosts and contains multiple polymorphic sites

Characterisation of microsporidian species and differentiation among genetic variants of the same species has typically relied on ribosomal RNA (rRNA) gene sequences. We characterised the entire rRNA gene of a microsporidium from 11 isolates representing eight different European bumblebee (Bombus) species. We demonstrate that the microsporidium Nosema bombi infected all hosts that originated from a wide geographic area. A total of 16 variable sites (all single nucleotid polymorphisms (SNPs)) was detected in the small subunit (SSU) rRNA gene and 42 (39 SNPs and 3 indels) in the large subunit (LSU) rRNA sequence. Direct sequencing of PCR-amplified DNA products of the internal transcribed spacer (ITS) region revealed identical sequences in all isolates. In contrast, ITS fragment length determined by PAGE and sequencing of cloned amplicons gave better resolution of sequences and revealed multiple SNPs across isolates and two fragment sizes in each isolate (six short and seven long amplicon variants). Genetic variants were not unique to individual host species. Moreover, two or more sequence variants were obtained from individual bumblebee hosts, suggesting the existence of multiple, variable copies of rRNA in the same microsporidium, and contrary to that expected for a class of multi-gene family under concerted evolution theory. Our data on within-genome rRNA variability call into question the usefulness of rRNA sequences to characterise intraspecific genetic variants in the Microsporidia and other groups of unicellular organisms.     

Morphological variation versus genetic divergence: a taxonomic implication for Mogannia species (Cicadidae: Cicadinae)

Species showing intraspecific morphological variation tend to be very difficult to identify using morphological characters. One such example is the cicada genus Mogannia where some species show considerable intraspecific variation mainly exhibited by wing pattern and body colouration. Thirty-one variants covering different putative species of Mogannia were recognized and illustrated in the present paper. Molecular data of mitochondrial COI and Cytb sequences were employed to test the level of variation and phylogeny of them. The existence of a ‘barcoding gap’ between intraspecific and interspecific genetic divergences and the reciprocally monophyletic clades indicate that all the closely related variants represent a single species, and that all these variants correspond to six species, respectively. However, the evolutionary relationships of intraspecific variants are not resolved possibly due to insufficient genetic variation among them. Our results indicated that some morphological characters, especially the wing pattern and body colouration, and even the number of apical processes of the aedeagus in a couple of related species, must be used with great caution in delimiting Mogannia species and their relatives. The factors responsible for intraspecific morphological variation and phylogeny of Mogannia spp. are preliminarily discussed.     

A trans-Amazonian screening of mtDNA reveals deep intraspecific divergence in forest birds and suggests a vast underestimation of species diversity

The Amazonian avifauna remains severely understudied relative to that of the temperate zone, and its species richness is thought to be underestimated by current taxonomy. Recent molecular systematic studies using mtDNA sequence reveal that traditionally accepted species-level taxa often conceal genetically divergent subspecific lineages found to represent new species upon close taxonomic scrutiny, suggesting that intraspecific mtDNA variation could be useful in species discovery. Surveys of mtDNA variation in Holarctic species have revealed patterns of variation that are largely congruent with species boundaries. However, little information exists on intraspecific divergence in most Amazonian species. Here we screen intraspecific mtDNA genetic variation in 41 Amazonian forest understory species belonging to 36 genera and 17 families in 6 orders, using 758 individual samples from Ecuador and French Guiana. For 13 of these species, we also analyzed trans-Andean populations from the Ecuadorian Chocó. A consistent pattern of deep intraspecific divergence among trans-Amazonian haplogroups was found for 33 of the 41 taxa, and genetic differentiation and genetic diversity among them was highly variable, suggesting a complex range of evolutionary histories. Mean sequence divergence within families was the same as that found in North American birds (13%), yet mean intraspecific divergence in Neotropical species was an order of magnitude larger (2.13% vs. 0.23%), with mean distance between intraspecific lineages reaching 3.56%. We found no clear relationship between genetic distances and differentiation in plumage color. Our results identify numerous genetically and phenotypically divergent lineages which may result in new species-level designations upon closer taxonomic scrutiny and thorough sampling, although lineages in the tropical region could be older than those in the temperate zone without necessarily representing separate species. In-depth phylogeographic surveys are urgently needed to avoid underestimating tropical diversity, and the use of mtDNA markers can be instrumental in identifying and prioritizing taxa for species discovery.