Inferences from a rapidly moving hybrid zone

Anartia fatima and Anartia amathea (Lepidoptera: Nymphalidae) are sister taxa whose ranges abut in a narrow hybrid zone in eastern Panama. At the center of the zone, hybrids are abundant, although deviations from Hardy-Weinberg and linkage disequilibria are strong, due in part to assortative mating. We measured differences across the zone in four wing color-pattern characters, three allozyme loci, and mitochondrial haplotype. Wing pattern, allozyme, and mitochondrial clines were coincident (i.e., had the same positions) and concordant (i.e.. all markers had similar cline shapes, about 28 km wide). Repeated samples demonstrated that the hybrid zone has been moving eastwards at an average rate of 2.5 km/year over the past 20 years, accompanied by an equivalent movement of the mtDNA cline. No introgression of mtDNA haplotypes were found in the 'wake" of the moving cline, as might be expected for a neutral marker. The concordance of morphological and mtDNA clines between 1994 and 2000, in spite of hybrid zone movement, suggests strong epistasis between the mitochondrial genome and nuclear loci. Cline movement is achieved mainly by pure fatima immigrating into amathea populations; hybrids had little effect, and were presumably outcompeted by fitter pure fatima genotypes. This movement can be explained if random dispersal of 7-19 km x gen(1/2) is coupled with a competitive advantage to A. fatima genomes of 2-5%. Hybrid zone motion is equivalent to Phase III of Wright's shifting balance. Hybrid zone movement has rarely been considered likely in the past, but our results show that it may be more important in biogeography and evolution than generally realized.     

The historical biogeography of two Caribbean butterflies (Lepidoptera: Heliconiidae) as inferred from genetic variation at multiple loci

Mitochondrial DNA and allozyme variation was examined in populations of two Neotropical butterflies, Heliconius charithonia and Dryas iulia. On the mainland, both species showed evidence of considerable gene flow over huge distances. The island populations, however, revealed significant genetic divergence across some, but not all, ocean passages. Despite the phylogenetic relatedness and broadly similar ecologies of these two butterflies, their intraspecific biogeography clearly differed. Phylogenetic analyses of mitochondrial DNA sequences revealed that populations of D. iulia north of St. Vincent are monophyletic and were probably derived from South America. By contrast, the Jamaican subspecies of H. charithonia rendered West Indian H. charithonia polyphyletic with respect to the mainland populations; thus, H. charithonia seems to have colonized the Greater Antilles on at least two separate occasions from Central America. Colonization velocity does not correlate with subsequent levels of gene flow in either species. Even where range expansion seems to have been instantaneous on a geological timescale, significant allele frequency differences at allozyme loci demonstrate that gene flow is severely curtailed across narrow ocean passages. Stochastic extinction, rapid (re)colonization, but low gene flow probably explain why, in the same species, some islands support genetically distinct and nonexpanding populations, while nearby a single lineage is distributed across several islands. Despite the differences, some common biogeographic patterns were evident between these butterflies and other West Indian taxa; such congruence suggests that intraspecific evolution in the West Indies has been somewhat constrained by earth history events, such as changes in sea level.     

Molecular phylogeography reveals island colonization history and diversification of western Indian Ocean sunbirds (Nectarinia: Nectariniidae)

We constructed a phylogenetic hypothesis for western Indian Ocean sunbirds (Nectarinia) and used this to investigate the geographic pattern of their diversification among the islands of the Indian Ocean. A total of 1309 bp of mitochondrial sequence data was collected from the island sunbird taxa of the western Indian Ocean region, combined with sequence data from a selection of continental (African and Asian) sunbirds. Topological and branch length information combined with estimated divergence times are used to present hypotheses for the direction and sequence of colonization events in relation to the geological history of the Indian Ocean region. Indian Ocean sunbirds fall into two well-supported clades, consistent with two independent colonizations from Africa within the last 3.9 million years. The first clade contains island populations representing the species Nectarinia notata, while the second includes Nectarinia souimanga, Nectarinia humbloti, Nectarinia dussumieri, and Nectarinia coquereli. With respect to the latter clade, application of Bremer's [Syst. Biol. 41 (1992) 436] ancestral areas method permits us to posit the Comoros archipelago as the point of initial colonization in the Indian Ocean. The subsequent expansion of the souimanga clade across its Indian Ocean range occurred rapidly, with descendants of this early expansion remaining on the Comoros and granitic Seychelles. The data suggest that a more recent expansion from Anjouan in the Comoros group led to the colonization of Madagascar by sunbirds representing the souimanga clade. In concordance with the very young geological age of the Aldabra group, the sunbirds of this archipelago have diverged little from the Madagascar population; this is attributed to colonization of the Aldabra archipelago in recent times, in one or possibly two or more waves originating from Madagascar. The overall pattern of sunbird radiation across Indian Ocean islands indicates that these birds disperse across ocean barriers with relative ease, but that their subsequent evolutionary success probably depends on a variety of factors including prior island occupation by competing species.     

A genetic linkage map of the mimetic butterfly Heliconius melpomene

Heliconius melpomene is a mimetic butterfly that exhibits great geographic variation in color pattern. We present here a genetic linkage map based on analysis of genetic markers in 73 individuals from a single F(2) family, offspring of a cross between H. m. cythera from western Ecuador and H. m. melpomene from French Guiana. A novel "three-step method" is described for the analysis of dominant markers in an F(2) cross, using outbred parental strains and taking advantage of the lack of crossing over in female Lepidoptera. This method is likely to prove useful for future mapping studies in outbred species with crossing over restricted to one sex, such as the Lepidoptera and Drosophila. The resulting linkage map has 21 linkage groups corresponding to the 21 chromosomes of H. melpomene and includes 219 AFLP markers, 23 microsatellites, 19 single-copy nuclear genes, and the color pattern switch genes Yb and Sb. The marker density is high, averaging >1/7 cM. The total map length is 1616 cM and the average chromosome length is 77 cM. The genome size of H. melpomene was estimated to be 292 Mb, giving a relationship of physical-to-map distance of 180 kb/cM. This map forms the basis for future comparative linkage analysis of color pattern evolution in Heliconius.     

Spatial genetic structure of Simarouba amara Aubl. (Simaroubaceae), a dioecious, animal-dispersed Neotropical tree, on Barro Colorado Island, Panama

Simarouba amara (Simaroubaceae) is a vertebrate-dispersed, insect-pollinated Neotropical tree found in lowland moist forest from upper Mesoamerica to the Amazon basin. We assessed the spatial genetic structure of S. amara within the 50-ha Forest Dynamics Plot on Barro Colorado Island in the Republic of Panama. A total of 300 individuals were genotyped using five microsatellite loci, representing 100 individuals with a dbh>or=10 cm, 100 individuals of 1-10 cm dbh, and 100 individuals of <1 cm dbh. The 200 individuals in the two larger size classes were also genotyped with 155 AFLP loci. Spatial autocorrelation analysis using Moran's Index detected significant genotypic association at the smallest distance classes for 1-10 cm dbh (0-20 m) and >10 cm dbh (0-40 m) size categories. Significant spatial autocorrelations were detected over larger scales (0-140 m) in <1 cm dbh individuals. The relatively weak genetic structure of S. amara, in comparison to other recent studies, may be explained by pollen and seed dispersal over the 50 ha plot, overlapping seed shadows, and postrecruitment mortality.     

Long-distance gene flow and cross-Andean dispersal of lowland rainforest bees (Apidae: Euglossini) revealed by comparative mitochondrial DNA phylogeography

Euglossine bees (Apidae; Euglossini) exclusively pollinate hundreds of orchid species and comprise up to 25% of bee species richness in neotropical rainforests. As one of the first studies of comparative phylogeography in a neotropical insect group, we performed a mitochondrial DNA (mtDNA)-based analysis of 14 euglossine species represented by populations sampled across the Andes and/or across the Amazon basin. The mtDNA divergences within species were consistently low; across the 12 monophyletic species the mean intraspecific divergence among haplotypes was 0.9% (range of means, 0-1.9%). The cytochrome oxidase 1 (CO1) divergence among populations separated by the Andes (N = 11 species) averaged 1.1% (range 0.0-2.0%). The mtDNA CO1 data set displayed homogeneous rates of nucleotide substitution, permitting us to infer dispersal across the cordillera long after the final Andean uplift based on arthropod molecular clocks of 1.2-1.5% divergence per million years. Gene flow across the 3000-km breadth of the Amazon basin was inferred from identical cross-Amazon haplotypes found in five species. Although mtDNA haplotypes for 12 of the 14 euglossine species were monophyletic, a reticulate CO1 phylogeny was recovered in Euglossa cognata and E. mixta, suggesting large ancestral populations and recent speciation. Reference to closely related outgroups suggested recent speciation for the majority of species. Phylogeographical structure across a broad spatial scale is weaker in euglossine bees than in any neotropical group previously examined, and may derive from a combination of Quaternary speciation, population expansion and/or long-distance gene flow.     

Population structure and biogeography of migratory freshwater fishes (Prochilodus: Characiformes) in major South American rivers

Mitochondrial DNA (mtDNA) sequences were used to infer the phylogenetic relationships of Prochilodus species in the Paraná, Amazonas, Orinoco, and Magdalena basins. Sequences of ATPase subunits 6 and 8 (total 840 bp) were obtained for 21 Prochilodus specimens from the four river systems. Using Semaprochilodus as an outgroup, phylogenetic analyses showed that: (i) each river basin contains a monophyletic group of mtDNA lineages; and (ii) the branching order places Magdalena in a basal position with subsequent branching of Orinoco, Amazon and Paraná. The mitochondrial control region was sequenced for 26 P. lineatus (from the Paraná basin) and six other Prochilodus specimens from the Magdalena, Orinoco and Amazon. All 26 control region haplotypes were unique with sequence divergence ranging from 0.3 to 3.6%. The control region phylogeny is well resolved but phylogenetic structure is not associated with geography. For example, mtDNA haplotypes from the upper Paraná (Mogui Guassú) and the upper Bermejo, separated by at least 2600 km, have close genealogical ties. Phylogeographic analyses, including nested clade analysis, suggest high levels of gene flow within this basin.     

c-mos variation in songbirds: molecular evolution, phylogenetic implications, and comparisons with mitochondrial differentiation

Nucleotide sequences from the c-mos proto-oncogene have previously been used to reconstruct the phylogenetic relationships between distantly related vertebrate taxa. To explore c-mos variation at shallower levels of avian divergence, we compared c-mos sequences from representative passerine taxa that span a range of evolutionary differentiation, from basal passerine lineages to closely allied genera. Phylogenetic reconstructions based on these c-mos sequences recovered topologies congruent with previous DNA-DNA hybridization-based reconstructions, with many nodes receiving high support, as indicated by bootstrap and reliability values. One exception was the relationship of Acanthisitta to the remaining passerines, where the c-mos-based searches indicated a three-way polytomy involving the Acanthisitta lineage and the suboscine and oscine passerine clades. We also compared levels of c-mos and mitochondrial differentiation across eight oscine passerine taxa and found that c-mos nucleotide substitutions accumulate at a rate similar to that of transversion substitutions in mitochondrial protein-coding genes. These comparisons suggest that nuclear-encoded loci such as c-mos provide a temporal window of phylogenetic resolution that overlaps the temporal range where mitochondrial protein-coding sequences have their greatest utility and that c-mos substitutions and mtDNA transversions can serve as complementary, informative, and independent phylogenetic markers for the study of avian relationships.     

Tropical Understory Piper Shrubs Maintain High Levels of Genotypic Diversity Despite Frequent Asexual Recruitment

Many plant species have the capacity to regenerate asexually by resprouting from stem and leaf fragments. In the pan‐tropical shrub genus Piper, this tendency is thought to be higher in shade‐tolerant than light‐demanding species, and to represent a trade‐off with annual seed production. Here we use molecular markers to identify clones in five Piper species varying in light requirements. We test predictions that (i) asexual recruitment success is highest in shade‐tolerant species, and (ii) that consequently, shade‐tolerant species are characterized by lower genotypic diversity than light‐demanding Piper. We found that two shade‐tolerant Piper species recruited asexually more frequently (36–42% of sampled shoots were of asexual origin) than, two light‐demanding and one shade‐tolerant species (0–26%). Furthermore, as predicted, genotypic diversity was negatively correlated with the frequency of asexual recruitment in the population. Nonetheless, genotypic diversity of Piper was high compared with other clonal plants. The proportion of unique genotypes found per population ranged from 0.58 to 1.0 and the genotypic Simpson's diversity ranged from 0.93 to 1.0 for all five species. Our results suggest that even though asexual reproduction plays an important role in maintaining local populations of Piper in the understory, it does not seem to reduce genotypic diversity to levels that will threaten these species ability to respond to environmental change. Abstract in Spanish is available in the online version of this article.