Chemical signals act as the main reproductive barrier between sister and mimetic Heliconius butterflies.

Colour pattern is the main trait that drives mate recognition between Heliconius species that are phylogenetically close. However, when this cue is compromised such as in cases of mimetic, sympatric and closely related species, alternative mating signals must evolve to ensure reproductive isolation and species integrity. The closely related species Heliconius melpomene malleti and H. timareta florencia occur in the same geographical region, and despite being co-mimics, they display strong reproductive isolation. In order to test which cues differ between species, and potentially contribute to reproductive isolation, we quantified differences in the wing phenotype and the male chemical profile. As expected, the wing colour pattern was indistinguishable between the two species, while the chemical profile of the androconial and genital males' extracts showed marked differences. We then conducted behavioural experiments to study the importance of these signals in mate recognition by females. In agreement with our previous results, we found that chemical blends and not wing colour pattern drive the preference of females for conspecific males. Also, experiments with hybrid males and females suggested an important genetic component for both chemical production and preference. Altogether, these results suggest that chemicals are the major reproductive barrier opposing gene flow between these two sister and co-mimic species.

Temporal changes in hamlet communities (Hypoplectrus spp., Serranidae) over 17 years.

Transect surveys of hamlet communities (Hypoplectrus spp., Serranidae) covering 14 000 m2 across 16 reefs off La Parguera, Puerto Rico, are presented and compared with a previous survey conducted in the year 2000. The hamlet community has noticeably changed over 17 years, with a > 30% increase in relative abundance of the yellowtail hamlet Hypoplectrus chlorurus on the inner reefs at the expense of the other hamlet species. The data also suggest that the density of H. chlorurus has declined and that its distribution has shifted towards shallower depths. Considering that H. chlorurus has been previously identified as one of the few fish showing a positive association with seawater turbidity on the inner reefs of La Parguera and that sedimentation of terrestrial origin has increased over recent decades on these reefs, it is proposed that turbidity may constitute an important but so far overlooked ecological driver of hamlet communities.

Recombination in the eggs and sperm in a simultaneously hermaphroditic vertebrate.

When there is no recombination (achiasmy) in one sex, it is in the heterogametic one. This observation is so consistent that it constitutes one of the few patterns in biology that may be regarded as a 'rule' and Haldane (Haldane 1922 J. Genet. 12, 101-109. (doi:10.1007/BF02983075)) proposed that it might be driven by selection against recombination in the sex chromosomes. Yet differences in recombination rates between the sexes (heterochiasmy) have also been reported in hermaphroditic species that lack sex chromosomes. In plants-the vast majority of which are hermaphroditic-selection at the haploid stage has been proposed to drive heterochiasmy. Yet few data are available for hermaphroditic animals, and barely any for hermaphroditic vertebrates. Here, we leverage reciprocal crosses between two black hamlets (Hypoplectrus nigricans, Serranidae), simultaneously hermaphroditic reef fishes from the wider Caribbean, to generate high-density egg- and sperm-specific linkage maps for each parent. We find globally higher recombination rates in the eggs, with dramatically pronounced heterochiasmy at the chromosome peripheries. We suggest that this pattern may be due to female meiotic drive, and that this process may be an important source of heterochiasmy in animals. We also identify a large non-recombining region that may play a role in speciation and local adaptation in Hypoplectrus.

Conservatism and novelty in the genetic architecture of adaptation in Heliconius butterflies.

Understanding the genetic architecture of adaptive traits has been at the centre of modern evolutionary biology since Fisher; however, evaluating how the genetic architecture of ecologically important traits influences their diversification has been hampered by the scarcity of empirical data. Now, high-throughput genomics facilitates the detailed exploration of variation in the genome-to-phenotype map among closely related taxa. Here, we investigate the evolution of wing pattern diversity in Heliconius, a clade of neotropical butterflies that have undergone an adaptive radiation for wing-pattern mimicry and are influenced by distinct selection regimes. Using crosses between natural wing-pattern variants, we used genome-wide restriction site-associated DNA (RAD) genotyping, traditional linkage mapping and multivariate image analysis to study the evolution of the architecture of adaptive variation in two closely related species: Heliconius hecale and H. ismenius. We implemented a new morphometric procedure for the analysis of whole-wing pattern variation, which allows visualising spatial heatmaps of genotype-to-phenotype association for each quantitative trait locus separately. We used the H. melpomene reference genome to fine-map variation for each major wing-patterning region uncovered, evaluated the role of candidate genes and compared genetic architectures across the genus. Our results show that, although the loci responding to mimicry selection are highly conserved between species, their effect size and phenotypic action vary throughout the clade. Multilocus architecture is ancestral and maintained across species under directional selection, whereas the single-locus (supergene) inheritance controlling polymorphism in H. numata appears to have evolved only once. Nevertheless, the conservatism in the wing-patterning toolkit found throughout the genus does not appear to constrain phenotypic evolution towards local adaptive optima.

Genomic atolls of differentiation in coral reef fishes (Hypoplectrus spp., Serranidae).

Because the vast majority of species are well diverged, relatively little is known about the genomic architecture of speciation during the early stages of divergence. Species within recent evolutionary radiations are often minimally diverged from a genomic perspective, and therefore provide rare opportunities to address this question. Here, we leverage the hamlet radiation (Hypoplectrus spp., brightly coloured reef fishes from the tropical western Atlantic) to characterize genomic divergence during the early stages of speciation. Transect surveys and spawning observations in Belize, Honduras and Panama confirm that sympatric barred (H. puella), black (H. nigricans) and butter (H. unicolor) hamlets are phenotypically distinct and reproductively isolated, although hybrid spawnings and individuals with intermediate phenotypes are seen on rare occasions. A survey of approximately 100 000 restriction site-associated SNPs in 126 samples from the three species across the three replicate populations reveals extremely slight genomewide divergence among species (FST  = 0.0038), indicating that ecomorphological differences and functional reproductive isolation are maintained in sympatry in a backdrop of extraordinary genomic similarity. Nonetheless, a very small proportion of SNPs (0.05% on average) are identified as FST outliers among sympatric species. Remarkably, a single SNP is identified as an outlier in repeated populations for the same species pair. A minicontig assembled de novo around this SNP falls into the genomic region containing the HoxCa10 and HoxCa11 genes in 10 teleost species, suggesting an important role for Hox gene evolution in this radiation. This finding, if confirmed, would provide a better understanding of the links between micro- and macroevolutionary processes.

On the spatial scale of dispersal in coral reef fishes.

Marine biologists have gone through a paradigm shift, from the assumption that marine populations are largely 'open' owing to extensive larval dispersal to the realization that marine dispersal is 'more restricted than previously thought'. Yet, population genetic studies often reveal low levels of genetic structure across large geographic areas. On the other side, more direct approaches such as mark-recapture provide evidence of localized dispersal. To what extent can direct and indirect studies of marine dispersal be reconciled? One approach consists in applying genetic methods that have been validated with direct estimates of dispersal. Here, we use such an approach-genetic isolation by distance between individuals in continuous populations-to estimate the spatial scale of dispersal in five species of coral reef fish presenting low levels of genetic structure across the Caribbean. Individuals were sampled continuously along a 220-km transect following the Mesoamerican Barrier Reef, population densities were estimated from surveys covering 17 200 m(2) of reef, and samples were genotyped at a total of 58 microsatellite loci. A small but positive isolation-by-distance slope was observed in the five species, providing mean parent-offspring dispersal estimates ranging between 7 and 42 km (CI 1-113 km) and suggesting that there might be a correlation between minimum/maximum pelagic larval duration and dispersal in coral reef fishes. Coalescent-based simulations indicate that these results are robust to a variety of dispersal distributions and sampling designs. We conclude that low levels of genetic structure across large geographic areas are not necessarily indicative of extensive dispersal at ecological timescales.

Sharp genetic discontinuity across a unimodal Heliconius hybrid zone.

Hybrid zones are powerful natural systems to study evolutionary processes to gain an understanding of adaptation and speciation. In the Cauca Valley (Colombia), two butterfly races, Heliconius cydno cydnides and Heliconius cydno weymeri, meet and hybridize. We characterized this hybrid zone using a combination of mitochondrial DNA (mtDNA) sequences, amplified fragment length polymorphisms (AFLPs), microsatellites and sequences for nuclear loci within and outside of the genomic regions that cause differences in wing colour pattern. The hybrid zone is largely composed of individuals of mixed ancestry. However, there is strong genetic discontinuity between the hybridizing races in mtDNA and, to a lesser extent, in all nuclear markers surveyed. The mtDNA clustering of H. c. cydnides with the H. cydno race from the Magdalena Valley and H. c. weymeri with the H. cydno race from the pacific coast suggests that H. c. cydnides colonized the Cauca Valley from the north, whereas H. c. weymeri did so by crossing the Andes in the southern part, implying a secondary contact origin. Colonization of the valley by H. cydno was accompanied by mimicry shift. Strong ecological isolation, driven by locally adaptive differences in mimetic wing patterns, is playing an important role in maintaining the hybrid zone. However, selection on wing pattern alone is not sufficient to explain the genetic discontinuity observed. There is evidence for differences in male mating preference, but the contribution of additional barriers needs further investigation. Overall, our results support the idea that speciation is a cumulative process, where the combination of multiple isolation barriers, combined with major phenotypic differences, facilitates population divergence in face of gene flow.

Butterfly genomics eclosing.

Technological and conceptual advances of the last decade have led to an explosion of genomic data and the emergence of new research avenues. Evolutionary and ecological functional genomics, with its focus on the genes that affect ecological success and adaptation in natural populations, benefits immensely from a phylogenetically widespread sampling of biological patterns and processes. Among those organisms outside established model systems, butterflies offer exceptional opportunities for multidisciplinary research on the processes generating and maintaining variation in ecologically relevant traits. Here we highlight research on wing color pattern variation in two groups of Nymphalid butterflies, the African species Bicyclus anynana (subfamily Satyrinae) and species of the South American genus Heliconius (subfamily Heliconiinae), which are emerging as important systems for studying the nature and origins of functional diversity. Growing genomic resources including genomic and cDNA libraries, dense genetic maps, high-density gene arrays, and genetic transformation techniques are extending current gene mapping and expression profiling analysis and enabling the next generation of research questions linking genes, development, form, and fitness. Efforts to develop such resources in Bicyclus and Heliconius underscore the general challenges facing the larger research community and highlight the need for a community-wide effort to extend ongoing functional genomic research on butterflies.

Multilocus molecular systematics of the circumtropical reef-fish genus Abudefduf (Pomacentridae): history, geography and ecology of speciation.

We investigated a pantropical sub-family and genus of damselfishes, the sergeant-majors (Pomacentridae: Abudefdufinae: Abudefduf), to identify the tempo and mechanisms of speciation in the lineage. We examined sequence capture data from 500 loci and 20 species, with multiple individuals sampled from across the geographic ranges of widespread species. Utilizing a maximum likelihood framework, as well as a time-calibrated Bayesian phylogeny, the following key questions are addressed: What is the historical tempo of speciation? What are the relative contributions of vicariant, peripatric and parapatric speciation to sergeant-major diversity? How is speciation related to major variation in trophic ecology? The approximately 20 species of sergeant-majors fall into three main lineages. The ancestral condition appears to be benthivory, which is predominant in two lineages comprising six species. The remaining species of sergeant-majors, of which there are at least 15, fall within a clade composed entirely of planktivores. This clade is sister to a benthivore clade that included one species, Abudefduf notatus, in transition to planktivory. Most speciation of sergeant-majors, which appeared ∼24 million years ago, occurred in the last 10 million years. Present distributional patterns indicate vicariant speciation precipitated by the closure of land barriers between both sides of the Atlantic and the Pacific, and the emergence of land between the Indian and Pacific Oceans. Within this backdrop, frequent oscillations in sea level over the last 10 million years also appear to have generated conditions suitable for both peripatric and vicariant speciation, and most speciation within the genus appears linked to these changes in sea level. Diversification within the genus has been concentrated in planktivorous seargeant-majors rather than benthivores. The root cause is unclear, but does not appear to be related to differences in dispersal potential, which is greater in the planktivorous species, due to the ability of their post-larval juveniles to raft with floating debris. This elevated speciation rate in planktivores and their propensity to form local endemics may reflect relaxation of selective pressures (e.g., on crypticity) that limit speciation in benthivorous sergeant-majors. Finally, our data allow us to clarify relationships of geminate sergeant-major species, indicating that there are subdivisions within the Atlantic for both benthivore and planktivore geminate pairs that may have misled previous studies.

Mimicry meets the mitochondrion. Evolution.

A recent molecular study of the evolution of mimicry in tropical butterflies of the genus Heliconius proves that the mimics adapted to previously diverged 'model' species, but does not clearly distinguish between opposing views of how the model species diverged.

Concordant evolutionary patterns among Indo-West Pacific butterflyfishes.

Genetic differences within a 495 base pair section of the mitochondrial cytochrome b gene reveal a striking concordance among species in two monophyletic groups of Indo-west Pacific butterflyfishes. In both species groups, an approximately 2.0% genetic break clearly partitions individuals between the Indian Ocean and Pacific Ocean. However, levels of intra-Pacific mtDNA variation are low, on average less than 1.0%, and fail to cluster by species boundaries defined by colour pattern. Individuals from different species, separated by thousands of kilometers, often possess identical cytochrome b sequences, whereas conspecifics from the same reefs can show up to 1.5% difference. The discrepancy between the mtDNA gene tree and species boundaries may reflect retained ancestral variation or may be the result of hybridization. The strong temporal and phylogenetic concordance between these two independent species groups suggests that genetic differentiation was influenced by common environmental factors. Low levels of within- and between-species genetic differences imply a recent divergence time and suggest a link between speciation within each group and Pleistocene climatic fluctuations. These results paint a turbulent picture of the recent evolutionary history of the Indo-West Pacific.

Development and characterization of 11 microsatellite loci for the Mona Island iguana (Cyclura cornuta stejnegeri).

We isolated and characterized 11 microsatellite loci in the Mona Island iguana (Cyclura cornuta stejnegeri). Eleven loci exhibit moderate to high allelic diversity (two to 12 alleles, mean = 4.5) and polymorphism (mean observed heterozygosity, 0.56; range, 0.26 to 0.78) in 41 adults. This marker set has low probability of identity and high parentage exclusion power and will be suitable for studies of paternity, social organization and relatedness in this species.

Heliconius wing patterns: an evo-devo model for understanding phenotypic diversity.

Evolutionary Developmental Biology aims for a mechanistic understanding of phenotypic diversity, and present knowledge is largely based on gene expression and interaction patterns from a small number of well-known model organisms. However, our understanding of biological diversification depends on our ability to pinpoint the causes of natural variation at a micro-evolutionary level, and therefore requires the isolation of genetic and developmental variation in a controlled genetic background. The colour patterns of Heliconius butterflies (Nymphalidae: Heliconiinae) provide a rich suite of naturally occurring variants with striking phenotypic diversity and multiple taxonomic levels of variation. Diversification in the genus is well known for its dramatic colour-pattern divergence between races or closely related species, and for Müllerian mimicry convergence between distantly related species, providing a unique system to study the development basis of colour-pattern evolution. A long history of genetic studies has showed that pattern variation is based on allelic combinations at a surprisingly small number of loci, and recent developmental evidence suggests that pattern development in Heliconius is different from the eyespot determination of other butterflies. Fine-scale genetic mapping studies have shown that a shared toolkit of genes is used to produce both convergent and divergent phenotypes. These exciting results and the development of new genomic resources make Heliconius a very promising evo-devo model for the study of adaptive change.

First-generation linkage map of the warningly colored butterfly Heliconius erato.

We report the first genetic linkage map of Heliconius erato, a species that shows remarkable variation in its warningly colored wing patterns. We use crosses between H. erato and its sister species, H. himera, to place two major color pattern genes, D and Cr, on a linkage map containing AFLP, allozyme, microsatellite and single-copy nuclear loci. We identified all 21 linkage groups in an initial genetic screen of 22 progeny from an F1 female x male H. himera family. Of the 229 markers, 87 used to identify linkage groups were also informative in 35 progeny from a sibling backcross (H. himera female x F1 male). With these, and an additional 33 markers informative in the second family, we constructed recombinational maps for 19 of the 21 linkage groups. These maps varied in length from 18.1 to 431.1 centimorgans (cM) and yielded an estimated total length of 2400 cM. The average distance between markers was 23 cM, and eight of the 19 linkage groups, including the sex chromosome (Z) and the chromosome containing the Cr locus, contained two or more codominant anchor loci. Of the three potential candidate genes mapped here, Cubitus interruptus (Ci), Decapentaplegic (Dpp) and Wingless (Wg), only Ci was linked, although loosely, to a known Heliconius color pattern locus. This work is an important first step for constructing a denser genetic map of the H. erato color pattern radiation and for a comparative genomic study of the architecture of mimicry in Heliconius butterflies.

Rapid rate of control-region evolution in Pacific butterflyfishes (Chaetodontidae).

Sequence differences in the tRNA-proline (tRNApro) end of the mitochondrial control-region of three species of Pacific butterflyfishes accumulated 33-43 times more rapidly than did changes within the mitochondrial cytochrome b gene (cytb). Rapid evolution in this region was accompanied by strong transition/transversion bias and large variation in the probability of a DNA substitution among sites. These substitution constraints placed an absolute ceiling on the magnitude of sequence divergence that could be detected between individuals. This divergence "ceiling" was reached rapidly and led to a decay in the relative rate of control-region/cytb b evolution. A high rate of evolution in this section of the control-region of butterflyfishes stands in marked contrast to the patterns reported in some other fish lineages. Although the mechanism underlying rate variation remains unclear, all taxa with rapid evolution in the 5'-end of the control-region showed extreme transition biases. By contrast, in taxa with slower control-region evolution, transitions accumulated at nearly the same rate as transversions. More information is needed to understand the relationship between nucleotide bias and the rate of evolution in the 5'-end of the control-region. Despite strong constraints on sequence change, phylogenetic information was preserved in the group of recently differentiated species and supported the clustering of sequences into three major mtDNA groupings. Within these groups, very similar control-region sequences were widely distributed across the Pacific Ocean and were shared between recognized species, indicating a lack of mitochondrial sequence monophyly among species.

The phylogeographic pattern of mitochondrial DNA variation in the Dall's porpoise Phocoenoides dalli.

We used 11 restriction endonucleases to study mtDNA variation in 101 Dall's porpoises Phocoenoides dalli from the Bering Sea and western North Pacific. There was little phylogeographic patterning among the 34 mtDNA haplotypes identified in this analysis, suggesting a strong historical connection among populations across this region. Nonetheless, mtDNA variation does not appear to be randomly distributed in this species. Both GST and AMOVA uncovered significant differences in the distribution of mtDNA variation between the Bering Sea and western North Pacific populations. These mtDNA results, coupled with differences in allozyme variation and parasite infestation, support the demographic distinctiveness of Bering Sea and western North Pacific stocks of Dall's porpoise. The lack of a strong phylogeographic orientation of mtDNA haplotypes within the Dall's porpoise is similar to the pattern reported in other vertebrates such as coyotes, blackbirds, chickadees, marine catfish, and catadromous eels. Like Dall's porpoise, these species are broadly distributed, and have large populations linked by moderate to high levels of gene flow. However, the more complex, deeply branched phylogenetic network of mtDNA haplotypes within Dall's porpoise, relative to these other vertebrates, suggests important differences between these species in the forces shaping mtDNA variation. One such force is the effective size of female populations, which appears to have been comparatively large and stable in Dall's porpoise.

Patterns of genetic diversity and biogeographical history of the tropical wetland tree, Pterocarpus officinalis (Jacq.), in the Caribbean basin.

Studies examining intraspecific variation in plant species with widespread distributions and disjunct populations have mainly concentrated on temperate species. Here, we determined the genetic structure of a broadly distributed wetland tropical tree, Pterocarpus officinalis (Jacq.), from eight Neotropical populations using amplified length fragment polymorphisms (AFLP). AFLPs proved highly variable with almost half (48%) of the genetic variation at these loci occurring among individuals within populations. Nonetheless, there was a strong geographical pattern in the distribution of AFLP variation within P. officinalis. Caribbean and continental populations fell into two well-defined genetic clusters supported by the presence of a number of unique AFLP bands. Within these two regions, there were also strong genetic differences among populations, caused mainly by frequency differences in AFLP bands, making it difficult to determine the evolutionary relationships among populations. In addition, our analysis of P. officinalis revealed striking differences in the levels of AFLP variation among the eight populations sampled. In general, Caribbean populations had lower genetic diversity than continental populations. Moreover, there was a clear loss in AFLP diversity with distance from the continent among Caribbean populations. The overall genetic pattern within P. officinalis suggests that past colonization history, coupled with genetic drift within local populations, rather than contemporary gene flow are the major forces shaping variation within this species.

What initiates speciation in passion-vine butterflies?

Studies of the continuum between geographic races and species provide the clearest insights into the causes of speciation. Here we report on mate choice and hybrid viability experiments in a pair of warningly colored butterflies, Heliconius erato and Heliconius himera, that maintain their genetic integrity in the face of hybridization. Hybrid sterility and inviability have been unimportant in the early stages of speciation of these two Heliconius. We find no evidence of reduced fecundity, egg hatch, or larval survival nor increases in developmental time in three generations of hybrid crosses. Instead, speciation in this pair appears to have been catalyzed by the association of strong mating preferences with divergence in warning coloration and ecology. In mate choice experiments, matings between the two species are a tenth as likely as matings within species. F1 hybrids of both sexes mate frequently with both pure forms. However, male F1 progeny from crosses between H. himera mothers and H. erato fathers have somewhat reduced mating success. The strong barrier to gene flow provided by divergence in mate preference is probably enhanced by frequency-dependent predation against hybrids similar to the type known to occur across interracial hybrid zones of H. erato. In addition, the transition between this pair falls at the boundary between wet and dry forest, and rare hybrids may also be selected against because they are poorly adapted to either biotope. These results add to a growing body of evidence that challenge the importance of genomic incompatibilities in the earliest stages of speciation.