Molecular systematics of the butterfly genus Ithomia (Lepidoptera: Ithomiinae): a composite phylogenetic hypothesis based on seven genes

Contrasting patterns of Andean diversification among three diverse clades of Neotropical clearwing butterflies

The Neotropical region is the most biodiverse on Earth, in a large part due to the highly diverse tropical Andean biota. The Andes are a potentially important driver of diversification within the mountains and for neighboring regions. We compared the role of the Andes in diversification among three subtribes of Ithomiini butterflies endemic to the Neotropics, Dircennina, Oleriina, and Godyridina. The diversification patterns of Godyridina have been studied previously. Here, we generate the first time‐calibrated phylogeny for the largest ithomiine subtribe, Dircennina, and we reanalyze a published phylogeny of Oleriina to test different biogeographic scenarios involving the Andes within an identical framework. We found common diversification patterns across the three subtribes, as well as major differences. In Dircennina and Oleriina, our results reveal a congruent pattern of diversification related to the Andes with an Andean origin, which contrasts with the Amazonian origin and multiple Andean colonizations of Godyridina. In each of the three subtribes, a clade diversified in the Northern Andes at a faster rate. Diversification within Amazonia occurred in Oleriina and Godyridina, while virtually no speciation occurred in Dircennina in this region. Dircennina was therefore characterized by higher diversification rates within the Andes compared to non‐Andean regions, while in Oleriina and Godyridina, we found no difference between these regions. Our results and discussion highlight the importance of comparative approaches in biogeographic studies.

Choosing and Using Introns in Molecular Phylogenetics

Introns are now commonly used in molecular phylogenetics in an attempt to recover gene trees that are concordant with species trees, but there are a range of genomic, logistical and analytical considerations that are infrequently discussed in empirical studies that utilize intron data. This review outlines expedient approaches for locus selection, overcoming paralogy problems, recombination detection methods and the identification and incorporation of LVHs in molecular systematics. A range of parsimony and Bayesian analytical approaches are also described in order to highlight the methods that can currently be employed to align sequences and treat indels in subsequent analyses. By covering the main points associated with the generation and analysis of intron data, this review aims to provide a comprehensive introduction to using introns (or any non-coding nuclear data partition) in contemporary phylogenetics.

North Andean origin and diversification of the largest ithomiine butterfly genus

The Neotropics harbour the most diverse flora and fauna on Earth. The Andes are a major centre of diversification and source of diversity for adjacent areas in plants and vertebrates, but studies on insects remain scarce, even though they constitute the largest fraction of terrestrial biodiversity. Here, we combine molecular and morphological characters to generate a dated phylogeny of the butterfly genus Pteronymia (Nymphalidae: Danainae), which we use to infer spatial, elevational and temporal diversification patterns. We first propose six taxonomic changes that raise the generic species total to 53, making Pteronymia the most diverse genus of the tribe Ithomiini. Our biogeographic reconstruction shows that Pteronymia originated in the Northern Andes, where it diversified extensively. Some lineages colonized lowlands and adjacent montane areas, but diversification in those areas remained scarce. The recent colonization of lowland areas was reflected by an increase in the rate of evolution of species' elevational ranges towards present. By contrast, speciation rate decelerated with time, with no extinction. The geological history of the Andes and adjacent regions have likely contributed to Pteronymia diversification by providing compartmentalized habitats and an array of biotic and abiotic conditions, and by limiting dispersal between some areas while promoting interchange across others.

Identification and comparison of gonadal transcripts of testis and ovary of adult common carp Cyprinus carpio using suppression subtractive hybridization

The limited number of gonad-specific and gonad-related genes that have been identified in fish represents a major obstacle in the study of fish gonad development and sex differentiation. In common carp Cyprinus carpio from China's Yellow River, the ovary and testis differ in volume and weight in adult fish of the same age. Comparing sperm, egg, and somatic cell transcripts in this carp may provide insight into the mechanisms of its gonad development and sex differentiation. In the present work, gene expression patterns in the carp ovary and testis were compared using suppression subtractive hybridization. Two bidirectional subtracted complementary DNA (cDNA) libraries were analyzed in parallel using testis or ovary as testers. Eighteen nonredundant clones were identified in the male library, including 15 known cDNAs. The expression patterns of selected genes in testis and ovary were analyzed using reverse transcriptase polymerase chain reaction. Tektin-1, GAPDS, FGFIBP, IGFBP-5, and an unknown gene from the Ccmg4 clone were observed to be expressed only in testis. GSDF, BMI1b, Wt1a, and an unknown gene from the Ccme2 clone were expressed at higher levels in testis than in ovary at sexual maturity. Thirty functional expressed sequence tags (ESTs) were identified in 43 sequenced clones in the female library, including 28 known cDNAs, one uncharacterized cDNA (EST clone), and one novel sequence. Eight identified ESTs showed significant differences in expression between the testis and the ovary. ZP3C and Psmb2 were expressed exclusively in ovary, whereas the expression levels of IFIPGL-1, Setd6, ATP-6, CDC45, AIF-1, and an unknown gene from the Ccfh2 clone were more strongly expressed in ovary than in testis. In addition, the expression of ZP3C, Wt1a, and Setd6 was analyzed in male and female gonads, heart, liver, kidney, and brain. ZP3C was expressed only in ovary. Setd6 expression was significantly stronger in female tissues than that in the male, except in the liver, and Wt1a expression showed sexual dimorphism in the kidney and liver. Results suggest that these genes could play key roles during carp growth, both in the gonad and other tissues. The results provide a resource for further investigation of molecular mechanisms responsible for gonad development and sex differentiation in Yellow River common carp.

Molecular systematics of the Simulium jenningsi species group (Diptera: Simuliidae), with three new fast-evolving nuclear genes for phylogenetic inference

A molecular phylogeny was inferred for the 22 nominal species of black flies in the Simulium jenningsi species group, which includes major pests of humans and livestock in North America. Females are structurally monomorphic, presenting a problem for identification of the pests. For each species, we sequenced approximately two kilobases from the mitochondrial genome (ND2, Cox I, proximal one-half of Cox II) and about six kilobases from the nuclear genome (ca. 2 kilobases each from 3 rapidly evolving nuclear genes: big zinc finger [BZF], "5-intron gene" [5intG], and elongation complex protein 1 [ECP1]) and analyzed them phylogenetically using maximum likelihood and Bayesian methods. The three nuclear loci have not previously been used in phylogenetic studies. The mitochondrial region recovered 6 group members as monophyletic. BZF, 5intG, and ECP1 sequences each permitted identification of 13 species and recovered the S. fibrinflatum and S. taxodium subgroups. Simulium aranti Stone and Snoddy and S. luggeri Nicholson and Mickel were consistently recovered at the base of the group. Simulium ozarkense Moulton and Adler, S. dixiense Stone and Snoddy, S. krebsorum Moulton and Adler, and S. haysi Stone and Snoddy branched off before two well-supported sister groups of the remaining species. This remainder consisted of species occupying slow, sandy lowland streams-S. definitum Moulton and Adler, S. jonesi Stone and Snoddy, and the S. taxodium subgroup (S. taxodium Snoddy and Beshear, S. chlorum Moulton and Adler, S. confusum Moulton and Adler, and S. lakei Snoddy)-as sister to two clades of species inhabiting swift, rocky upland streams-the S. fibrinflatum subgroup (S. fibrinflatum Twinn, S. notiale Stone and Snoddy, and S. snowi Stone and Snoddy) and a clade comprised of S. anchistinum Moulton and Adler, S. jenningsi Malloch, and S. nyssa Stone and Snoddy, plus species having cocoons without anterolateral apertures (S. infenestrum Moulton and Adler, S. podostemi Snoddy, S. penobscotense Snoddy and Bauer, and S. remissum Moulton and Adler). Simulium snowi Stone and Snoddy is here considered a synonym of S. notiale Stone and Snoddy. Trees inferred from BZF and 5intG were largely concordant with those from ECP1, but slightly less resolved. Combining mitochondrial and nuclear data sets did not greatly improve the performance of the ECP1 data set alone. We, therefore, propose ECP1 as the gold standard for identification of members of the S. jenningsi group. Maximum likelihood analysis of combined sequences from all three nuclear genes, with three morphological constraints imposed, yielded a tree proposed as the best hypothesis of relationships among group members, based on all available data.

Next-generation sequencing for high-throughput molecular ecology: a step-by-step protocol for targeted multilocus genotyping by pyrosequencing

Next-generation sequencing technology can now provide population biologists and phylogeographers with information at the genomic scale; however, many pertinent questions in population genetics and phylogeography can be answered effectively with modest levels of genomic information. For the past two decades, most population-level studies have lacked nuclear DNA (nDNA) sequence data due to the complications and cost of amplifying and sequencing diploid loci. However, pyrosequencing of emulsion PCR reactions, amplifying from only one molecule at a time, can generate megabases of clonally amplified loci at high coverage, thereby greatly simplifying allelic sequence determination. Here, we present a step-by-step methodology for utilizing the 454 GS FLX Titanium pyrosequencing platform to simultaneously sequence 16 populations (at 20 individuals per population) at 10 different nDNA loci (3,200 loci in total) in one plate of sequencing for less than the cost of traditional Sanger sequencing.

Next-generation phylogeography: a targeted approach for multilocus sequencing of non-model organisms

The field of phylogeography has long since realized the need and utility of incorporating nuclear DNA (nDNA) sequences into analyses. However, the use of nDNA sequence data, at the population level, has been hindered by technical laboratory difficulty, sequencing costs, and problematic analytical methods dealing with genotypic sequence data, especially in non-model organisms. Here, we present a method utilizing the 454 GS-FLX Titanium pyrosequencing platform with the capacity to simultaneously sequence two species of sea star (Meridiastra calcar and Parvulastra exigua) at five different nDNA loci across 16 different populations of 20 individuals each per species. We compare results from 3 populations with traditional Sanger sequencing based methods, and demonstrate that this next-generation sequencing platform is more time and cost effective and more sensitive to rare variants than Sanger based sequencing. A crucial advantage is that the high coverage of clonally amplified sequences simplifies haplotype determination, even in highly polymorphic species. This targeted next-generation approach can greatly increase the use of nDNA sequence loci in phylogeographic and population genetic studies by mitigating many of the time, cost, and analytical issues associated with highly polymorphic, diploid sequence markers.

Phylogenetics and biogeography of a spectacular Old World radiation of butterflies: the subtribe Mycalesina (Lepidoptera: Nymphalidae: Satyrini)

Background

Butterflies of the subtribe Mycalesina (Nymphalidae: Satyrinae) are important model organisms in ecology and evolution. This group has radiated spectacularly in the Old World tropics and presents an exciting opportunity to better understand processes of invertebrate rapid radiations. However, the generic-level taxonomy of the subtribe has been in a constant state of flux, and relationships among genera are unknown. There are six currently recognized genera in the group. Mycalesis, Lohora and Nirvanopsis are found in the Oriental region, the first of which is the most speciose genus among mycalesines, and extends into the Australasian region. Hallelesis and Bicyclus are found in mainland Africa, while Heteropsis is primarily Madagascan, with a few species in Africa. We infer the phylogeny of the group with data from three genes (total of 3139 bp) and use these data to reconstruct events in the biogeographic history of the group.

Results

The results indicate that the group Mycalesina radiated rapidly around the Oligocene-Miocene boundary. Basal relationships are unresolved, but we recover six well-supported clades. Some species of Mycalesis are nested within a primarily Madagascan clade of Heteropsis, while Nirvanopsis is nested within Lohora. The phylogeny suggests that the group had its origin either in Asia or Africa, and diversified through dispersals between the two regions, during the late Oligocene and early Miocene. The current dataset tentatively suggests that the Madagascan fauna comprises two independent radiations. The Australasian radiation shares a common ancestor derived from Asia. We discuss factors that are likely to have played a key role in the diversification of the group.

Conclusions

We propose a significantly revised classification scheme for Mycalesina. We conclude that the group originated and radiated from an ancestor that was found either in Asia or Africa, with dispersals between the two regions and to Australasia. Our phylogeny paves the way for further comparative studies on this group that will help us understand the processes underlying diversification in rapid radiations of invertebrates.

Molecular phylogenetics of the neotropical butterfly subtribe Oleriina (Nymphalidae: Danainae: Ithomiini)

The Oleriina is one of the most speciose subtribes of the neotropical nymphalid butterfly tribe Ithomiini. They are widely distributed across the Andes and Amazonian lowlands and like other ithomiines they are involved in complex mimicry rings. This subtribe is of particular interest because it contains the most diverse ithomiine genus, Oleria, as well as two genera, Megoleria and Hyposcada, that feed on hostplants not utilized elsewhere in the tribe. Here we present the first comprehensive species-level phylogeny for the Oleriina, representing 83% of recognised species in the group, and based on 6698bp from eight mitochondrial (mt) and nuclear (nc) genes. Topologies are largely congruent for ncDNA and the concatenated dataset and the genera Oleria, Hyposcada and Megoleria are recovered and well-supported, although strongly discordant genealogy between mtDNA and ncDNA suggest possible introgression among Hyposcada and Megoleria. A fourth clade containing the type species of Ollantaya is consistently recovered, and this recently synonymized name is resurrected. Clear subdivisions within Oleria separate the genus into four species groups, onega, amalda, makrena and aegle, which also correspond to differing biogeographic and elevation range characteristics. Unlike other ithomiine genera, the Oleriina show homogeneity in mimetic wing pattern, in sharp contrast to the emerging paradigm that mimetic shifts have enhanced diversification in the tribe. Our results show a potentially more important role for geographic isolation in the diversification of the Oleriina compared to other Ithomiini studied to date and provide a framework for more detailed biogeographical studies, in addition to a rare opportunity for comparative analyses with other neotropical groups.

Mitochondrial DNA barcoding detects some species that are real, and some that are not

Mimicry and extensive geographical subspecies polymorphism combine to make species in the ithomiine butterfly genus Mechanitis (Lepidoptera; Nymphalidae) difficult to determine. We use mitochondrial DNA (mtDNA) barcoding, nuclear sequences and amplified fragment length polymorphism (AFLP) genotyping to investigate species limits in this genus. Although earlier biosystematic studies based on morphology described only four species, mtDNA barcoding revealed eight well-differentiated haplogroups, suggesting the presence of four new putative 'cryptic species'. However, AFLP markers supported only one of these four new 'cryptic species' as biologically meaningful. We demonstrate that in this genus, deep genetic divisions expected on the basis of mtDNA barcoding are not always reflected in the nuclear genome, and advocate the use of AFLP markers as a check when mtDNA barcoding gives unexpected results.

Natural selection in action during speciation

The role of natural selection in speciation, first described by Darwin, has finally been widely accepted. Yet, the nature and time course of the genetic changes that result in speciation remain mysterious. To date, genetic analyses of speciation have focused almost exclusively on retrospective analyses of reproductive isolation between species or subspecies and on hybrid sterility or inviability rather than on ecologically based barriers to gene flow. However, if we are to fully understand the origin of species, we must analyze the process from additional vantage points. By studying the genetic causes of partial reproductive isolation between specialized ecological races, early barriers to gene flow can be identified before they become confounded with other species differences. This population-level approach can reveal patterns that become invisible over time, such as the mosaic nature of the genome early in speciation. Under divergent selection in sympatry, the genomes of incipient species become temporary genetic mosaics in which ecologically important genomic regions resist gene exchange, even as gene flow continues over most of the genome. Analysis of such mosaic genomes suggests that surprisingly large genomic regions around divergently selected quantitative trait loci can be protected from interrace recombination by "divergence hitchhiking." Here, I describe the formation of the genetic mosaic during early ecological speciation, consider the establishment, effects, and transitory nature of divergence hitchhiking around key ecologically important genes, and describe a 2-stage model for genetic divergence during ecological speciation with gene flow.

Pinniped phylogenetic relationships inferred using AFLP markers

Amplified fragment length polymorphisms (AFLPs) are widely used for phylogenetic reconstruction in plants but their use in animal taxa has been limited, and in mammals rare. In addition, their use has been largely limited to shallow relationships amongst species or subspecies. Here, we genotype 23 pinniped species for 310 AFLP markers and find a strong phylogenetic signal, with individuals coclustering within species, and overall a good agreement between our phylogeny and those constructed using mitochondrial DNA and nuclear sequences even at nodes approximately 15 million years old. Although supporting the existing ideas about pinniped relationships, our data shed light on relationships within the hitherto relatively unresolved Phocine species group, and provide further supporting evidence for raising two subspecies of Zalophus californianus, Z. c. californianus and Z. c. wollebaeki, to species level. Plotting AFLP divergence time estimates against those based on both mtDNA and nuclear sequences we find strong linear relationships, suggesting that the different markers are evolving in a clocklike fashion. These data further emphasize the utility of AFLP markers as general tools for phylogenetic reconstruction.

Genomic outposts serve the phylogenomic pioneers: designing novel nuclear markers for genomic DNA extractions of lepidoptera

Increasing the number of characters used in phylogenetic studies is the next crucial step towards generating robust and stable phylogenetic hypotheses - i.e., strongly supported and consistent across reconstruction method. Here we describe a genomic approach to finding new protein-coding genes for systematics in nonmodel taxa, which can be PCR amplified from standard, slightly degraded genomic DNA extracts. We test this approach on Lepidoptera, searching the draft genomic sequence of the silk moth Bombyx mori, for exons > 500 bp in length, removing annotated gene families, and compared remaining exons with butterfly EST databases to identify conserved regions for primer design. These primers were tested on a set of 65 taxa primarily in the butterfly family Nymphalidae. We were able to identify and amplify six previously unused gene regions (Arginine Kinase, GAPDH, IDH, MDH, RpS2, and RpS5) and two rarely used gene regions (CAD and DDC) that when added to the three traditional gene regions (COI, EF-1alpha and wingless) gave a data set of 8114 bp. Phylogenetic robustness and stability increased with increasing numbers of genes. Smaller taxanomic subsets were also robust when using the full gene data set. The full 11-gene data set was robust and stable across reconstruction methods, recovering the major lineages and strongly supporting relationships within them. Our methods and insights should be applicable to taxonomic groups having a single genomic reference species and several EST databases from taxa that diverged less than 100 million years ago.

The genetic mosaic suggests a new role for hitchhiking in ecological speciation

Early in ecological speciation, the genomically localized effects of divergent selection cause heterogeneity among loci in divergence between incipient species. We call this pattern of genomic variability in divergence the 'genetic mosaic of speciation'. Previous studies have used F(ST) outliers as a way to identify divergently selected genomic regions, but the nature of the relationship between outlier loci and quantitative trait loci (QTL) involved in reproductive isolation has not yet been quantified. Here, we show that F(ST) outliers between a pair of incipient species are significantly clustered around QTL for traits that cause ecologically based reproductive isolation. Around these key QTL, extensive 'divergence hitchhiking' occurs because reduced inter-race mating and negative selection decrease the opportunity for recombination between chromosomes bearing different locally adapted QTL alleles. Divergence hitchhiking is likely to greatly increase the opportunity for speciation in populations that are sympatric, regardless of whether initial divergence was sympatric or allopatric. Early in ecological speciation, analyses of population structure, gene flow or phylogeography based on different random or arbitrarily chosen neutral markers should be expected to conflict--only markers in divergently selected genomic regions will reveal the evolutionary history of adaptive divergence and ecologically based reproductive isolation. Species retain mosaic genomes for a very long time, and gene exchange in hybrid zones can vary dramatically among loci. However, in hybridizing species, the genomic regions that affect ecologically based reproductive isolation are difficult to distinguish from regions that have diverged for other reasons.

The tektin family of microtubule-stabilizing proteins

Tektins are insoluble a-helical proteins essential for the construction of cilia and flagella and are found throughout the eukaryotes apart from higher plants.

Forest-obligate Sabethes mosquitoes suggest palaeoecological perturbations

The origin of tropical forest diversity has been hotly debated for decades. Although specific mechanisms vary, many such explanations propose some vicariance in the distribution of species during glacial cycles and several have been supported by genetic evidence in Neotropical taxa. However, no consensus exists with regard to the extent or time frame of the vicariance events. Here, we analyse the cytochrome oxidase II mitochondrial gene of 250 Sabethes albiprivus B mosquitoes sampled from western Sao Paulo in Brazil. There was very low population structuring among collection sites (Phi(ST)=0.03, P=0.04). Historic demographic analyses and the contemporary geographic distribution of genetic diversity suggest that the populations sampled are not at demographic equilibrium. Three distinct mitochondrial clades were observed in the samples, one of which differed significantly in its geographic distribution relative to the other two within a small sampling area (approximately 70 x 35 km). This fact, supported by the inability of maximum likelihood analyses to achieve adequate fits to simple models for the population demography of the species, suggests a more complex history, possibly involving disjunct forest refugia. This hypothesis is supported by a genetic signal of recent population growth, which is expected if population sizes of this forest-obligate insect increased during the forest expansions that followed glacial periods. Although a time frame cannot be reliably inferred for the vicariance event leading to the three genetic clades, molecular clock estimates place this at approximately 1 Myr before present.

Limited performance of DNA barcoding in a diverse community of tropical butterflies

DNA 'barcoding' relies on a short fragment of mitochondrial DNA to infer identification of specimens. The method depends on genetic diversity being markedly lower within than between species. Closely related species are most likely to share genetic variation in communities where speciation rates are rapid and effective population sizes are large, such that coalescence times are long. We assessed the applicability of DNA barcoding (here the 5' half of the cytochrome c oxidase I) to a diverse community of butterflies from the upper Amazon, using a group with a well-established morphological taxonomy to serve as a reference. Only 77% of species could be accurately identified using the barcode data, a figure that dropped to 68% in species represented in the analyses by more than one geographical race and at least one congener. The use of additional mitochondrial sequence data hardly improved species identification, while a fragment of a nuclear gene resolved issues in some of the problematic species. We acknowledge the utility of barcodes when morphological characters are ambiguous or unknown, but we also recommend the addition of nuclear sequence data, and caution that species-level identification rates might be lower in the most diverse habitats of our planet.

The phylogenetic pattern of speciation and wing pattern change in neotropical Ithomia butterflies (Lepidoptera: nymphalidae)

Species level phylogenetic hypotheses can be used to explore patterns of divergence and speciation. In the tropics, speciation is commonly attributed to either vicariance, perhaps within climate-induced forest refugia, or ecological speciation caused by niche adaptation. Mimetic butterflies have been used to identify forest refugia as well as in studies of ecological speciation, so they are ideal for discriminating between these two models. The genus Ithomia contains 24 species of warningly colored mimetic butterflies found in South and Central America, and here we use a phylogenetic hypothesis based on seven genes for 23 species to investigate speciation in this group. The history of wing color pattern evolution in the genus was reconstructed using both parsimony and likelihood. The ancestral pattern for the group was almost certainly a transparent butterfly, and there is strong evidence for convergent evolution due to mimicry. A punctuationist model of pattern evolution was a significantly better fit to the data than a gradualist model, demonstrating that pattern changes above the species level were associated with cladogenesis and supporting a model of ecological speciation driven by mimicry adaptation. However, there was only one case of sister species unambiguously differing in pattern, suggesting that some recent speciation events have occurred without pattern shifts. The pattern of geographic overlap between clades over time shows that closely related species are mostly sympatric or, in one case, parapatric. This is consistent with modes of speciation with ongoing gene flow, although rapid range changes following allopatric speciation could give a similar pattern. Patterns of lineage accumulation through time differed significantly from that expected at random, and show that most of the extant species were present by the beginning of the Pleistocene at the latest. Hence Pleistocene refugia are unlikely to have played a major role in Ithomia diversification.

Tektin interactions and a model for molecular functions

Tektins from echinoderm flagella were analyzed for microheterogeneity, self-associations and association with tubulin, resulting in a general model of tektin filament structure and function applicable to most eukaryotic cilia and flagella. Using a new antibody to tektin consensus peptide RPNVELCRD, well-characterized chain-specific antibodies and quantitative gel densitometry, tektins A, B and C were found to be present in equimolar amounts in Sarkosyl-urea-stable filaments. In addition, two isoforms of tektin A are present in half-molar ratios to tektins B and C. Cross-linking of AB filaments indicates in situ nearest neighbor associations of tektin A1B and A2B heterodimers, -trimers, -tetramers and higher oligomers. Soluble purified tektin C is cross-linked as homodimers, trimers and tetramers, but not higher oligomers. Tektin filaments associate with both loosely bound and tightly bound tubulin, and with the latter in a 1:1 molar ratio, implying a specific, periodic association of tightly bound tubulin along the tektin axis. Similarly, in tektin-containing Sarkosyl-stable protofilament ribbons, two polypeptides ( approximately 67/73 kDa, homologues of rib72, efhc1 and efhc2) are present in equimolar ratios to each other and to individual tektins, co-fractionating with loosely bound tubulin. These results suggest a super-coiled arrangement of tektin filaments, the organization of which has important implications for the evolution, assembly and functions of cilia and flagella.

Strikingly variable divergence times inferred across an Amazonian butterfly 'suture zone'

'Suture zones' are areas where hybrid and contact zones of multiple taxa are clustered. Such zones have been regarded as strong evidence for allopatric divergence by proponents of the Pleistocene forest refugia theory, a vicariance hypothesis frequently used to explain diversification in the Amazon basin. A central prediction of the refugia and other vicariance theories is that the taxa should have a common history so that divergence times should be coincident among taxa. A suture zone for Ithomiinae butterflies near Tarapoto, NE Peru, was therefore studied to examine divergence times of taxa in contact across the zone. We sequenced 1619bp of the mitochondrial COI/COII region in 172 individuals of 31 species from across the suture zone. Inferred divergence times differed remarkably, with divergence between some pairs of widespread species (each of which may have two or more subspecies interacting in the zone, as in the genus Melinaea) being considerably less than that between hybridizing subspecies in other genera (for instance in Oleria). Our data therefore strongly refute a simple hypothesis of simultaneous vicariance and suggest that ongoing parapatric or other modes of differentiation in continuous forest may be important in driving diversification in Amazonia.

Linking patterns and processes of species diversification in the cone flies Strobilomyia (Diptera: Anthomyiidae)

Phytophagous insects provide useful models for the study of ecological speciation. Much attention has been paid to host shifts, whereas situations where closely related lineages of insects use the same plant during different time periods have been relatively neglected in previous studies of insect diversification. Flies of the genus Strobilomyia are major pests of conifers in Eurasia and North America. They are specialized feeders in cones and seeds of Abies (fir), Larix (larch) ,and Picea (spruce). This close association is accompanied by a large number of sympatric Strobilomyia species coexisting within each tree genus. We constructed a molecular phylogeny with a 1320 base-pair fragment of mitochondrial DNA that demonstrated contrasting patterns of speciation in larch cone flies, as opposed to spruce and fir cone flies; this despite their comparable geographic distributions and similar resource quality of the host. Species diversity is the highest on larch, and speciation is primarily driven by within-host phenological shifts, followed by allopatric speciation during geographical expansion. By contrast, fewer species exploit spruce and fir, and within-host phenological shifts did not occur. This study illustrates within-host adaptive radiation through phenological shifts, a neglected mode of sympatric speciation.