Introgression of wing pattern alleles and speciation via homoploid hybridization in Heliconius butterflies: a review of evidence from the genome

Molecular phylogeny of Hesperiidae (Lepidoptera) with an emphasis on Asian and African genera

Despite considerable research efforts in recent years, the deeper phylogenetic relationships among skipper butterflies (Hesperiidae) remain unresolved. This is primarily because of limited sampling, especially within Asian and African lineages. In this study, we consolidated previous data and extensively sampled Asian and African taxa to elucidate the phylogenetic relationships within Hesperiidae. The molecular dataset comprised sequences from two mitochondrial and two nuclear gene regions from 563 species that represented 353 genera. Our analyses revealed seven subfamilies within Hesperiidae: Coeliadinae, Euschemoninae, Eudaminae, Pyrginae, Heteropterinae, Trapezitinae, and Hesperiinae. The systematics of most tribes and genera aligned with those of prior studies. However, notable differences were observed in several tribes and genera. Overall, the position of taxa assigned to incertae sedis in Hesperiinae is largely clarified in this study. Our results strongly support the monophyly of the tribe Tagiadini (Pyrginae), and the systematics of some genera are clarified with comprehensive discussion. We recognize 15 tribes within the subfamily Hesperiinae. Of these, nine tribes are discussed in detail: Aeromachini, Astictopterini, Erionotini, Unkanini (new status), Ancistroidini, Ismini (confirmed status), Plastingini (new status), Gretnini (confirmed status), and Eetionini (confirmed status). We propose four subtribes within Astictopterini: Hypoleucina subtrib.n., Aclerosina, Cupithina, and Astictopterina. Furthermore, we describe a new genus (Hyarotoidesgen.n.) and reinstate two genera (Zeareinst.stat. and Separeinst.stat.) as valid. Additionally, we propose several new combinations: Zea mythecacomb.n.,Sepa bononiacomb.n. & reinst.stat., and Sepa umbrosacomb.n. This study, with extensive sampling of Asian and African taxa, greatly enhances the understanding of the knowledge of the skipper tree of life.

Comparative transcriptomics of albino and warningly-coloured caterpillars

Coloration is perhaps one of the most prominent adaptations for survival and reproduction of many taxa. Coloration is of particular importance for aposematic species, which rely on their coloring and patterning acting as a warning signal to deter predators. Most research has focused on the evolution of warning coloration by natural selection. However, little information is available for color mutants of aposematic species, particularly at the genomic level. Here, I compare the transcriptomes of albino mutant caterpillars of the aposematic wood tiger moth (Arctia plantaginis) to those of their full sibs having their distinctive orange-black warning coloration. The results showed >290 differentially expressed genes genome-wide. Genes involved in the immune system, structural constituents of cuticular, and immunity were mostly downregulated in the albino caterpillars. Surprisingly, higher expression was observed in core melanin genes from albino caterpillars, suggesting that melanin synthesis may be disrupted in terminal ends of the pathway during its final conversion. Taken together, these results suggest that caterpillar albinism may not be due to a depletion of melanin precursor genes. In contrast, the albino condition may result from the combination of faulty melanin conversion late in its synthesis and structural deficiencies in the cuticular preventing its deposition. The results are discussed in the context of how albinism may impact individuals of aposematic species in the wild.

Genomics of Evolutionary Novelty in Hybrids and Polyploids

It has long been recognized that hybridization and polyploidy are prominent processes in plant evolution. Although classically recognized as significant in speciation and adaptation, recognition of the importance of interspecific gene flow has dramatically increased during the genomics era, concomitant with an unending flood of empirical examples, with or without genome doubling. Interspecific gene flow is thus increasingly thought to lead to evolutionary innovation and diversification, via adaptive introgression, homoploid hybrid speciation and allopolyploid speciation. Less well understood, however, are the suite of genetic and genomic mechanisms set in motion by the merger of differentiated genomes, and the temporal scale over which recombinational complexity mediated by gene flow might be expressed and exposed to natural selection. We focus on these issues here, considering the types of molecular genetic and genomic processes that might be set in motion by the saltational event of genome merger between two diverged species, either with or without genome doubling, and how these various processes can contribute to novel phenotypes. Genetic mechanisms include the infusion of new alleles and the genesis of novel structural variation including translocations and inversions, homoeologous exchanges, transposable element mobilization and novel insertional effects, presence-absence variation and copy number variation. Polyploidy generates massive transcriptomic and regulatory alteration, presumably set in motion by disrupted stoichiometries of regulatory factors, small RNAs and other genome interactions that cascade from single-gene expression change up through entire networks of transformed regulatory modules. We highlight both these novel combinatorial possibilities and the range of temporal scales over which such complexity might be generated, and thus exposed to natural selection and drift.

Mechanisms for Color Convergence in a Mimetic Radiation of Poison Frogs

In animals, bright colors often evolve to mimic other species when a resemblance is selectively favored. Understanding the proximate mechanisms underlying such color mimicry can give insights into how mimicry evolves-for example, whether color convergence evolves from a shared set of mechanisms or through the evolution of novel color production mechanisms. We studied color production mechanisms in poison frogs (Dendrobatidae), focusing on the mimicry complex of Ranitomeya imitator. Using reflectance spectrometry, skin pigment analysis, electron microscopy, and color modeling, we found that the bright colors of these frogs, both within and outside the mimicry complex, are largely structural and produced by iridophores but that color production depends crucially on interactions with pigments. Color variation and mimicry are regulated predominantly by iridophore platelet thickness and, to a lesser extent, concentration of the red pteridine pigment drosopterin. Compared with each of the four morphs of model species that it resembles, R. imitator displays greater variation in both structural and pigmentary mechanisms, which may have facilitated phenotypic divergence in this species. Analyses of nonmimetic dendrobatids in other genera demonstrate that these mechanisms are widespread within the family and that poison frogs share a complex physiological "color palette" that can produce diverse and highly reflective colors.

Missing data, clade support and "reticulation": the molecular systematics of Heliconius and related genera (Lepidoptera: Nymphalidae) re-examined

Kozak et al. (2015, Syst. Biol., 64: 505) portrayed the inference of evolutionary history among Heliconius and allied butterfly genera as a particularly difficult problem for systematics due to prevalent gene conflict caused by interspecific reticulation. To control for this, Kozak et al. conducted a series of multispecies coalescent phylogenetic analyses that they claimed revealed pervasive conflict among markers, but ultimately chose as their preferred hypothesis a phylogenetic tree generated by the traditional supermatrix approach. Intrigued by this seemingly contradictory set of conclusions, we conducted further analyses focusing on two prevalent aspects of the data set: missing data and the uneven contribution of phylogenetic signal among markers. Here, we demonstrate that Kozak et al. overstated their findings of reticulation and that evidence of gene-tree conflict is largely lacking. The distribution of intrinsic homoplasy and incongruence homoplasy in their data set does not follow the pattern expected if phylogenetic history had been obscured by pervasive horizontal gene flow; in fact, noise within individual gene partitions is ten times higher than the incongruence among gene partitions. We show that the patterns explained by Kozak et al. as a result of reticulation can be accounted for by missing data and homoplasy. We also find that although the preferred topology is resilient to missing data, measures of support are sensitive to, and strongly eroded by too many empty cells in the data matrix. Perhaps more importantly, we show that when some taxa are missing almost all characters, adding more genes to the data set provides little or no increase in support for the tree.

Quantified reproductive isolation in Heliconius butterflies: Implications for introgression and hybrid speciation

Heliconius butterflies have become a model for the study of speciation with gene flow. For adaptive introgression to take place, there must be incomplete barriers to gene exchange that allow interspecific hybridization and multiple generations of backcrossing. The recent publication of estimates of individual components of reproductive isolation between several species of butterflies in the Heliconius melpomene-H. cydno clade allowed us to calculate total reproductive isolation estimates for these species. According to these estimates, the butterflies are not as promiscuous as has been implied. Differences between species are maintained by intrinsic mechanisms, while reproductive isolation of geographical races within species is mainly due to allopatry. We discuss the implications of this strong isolation for basic aspects of the hybrid speciation with introgression hypothesis.

Living Organisms Author Their Read-Write Genomes in Evolution

Evolutionary variations generating phenotypic adaptations and novel taxa resulted from complex cellular activities altering genome content and expression: (i) Symbiogenetic cell mergers producing the mitochondrion-bearing ancestor of eukaryotes and chloroplast-bearing ancestors of photosynthetic eukaryotes; (ii) interspecific hybridizations and genome doublings generating new species and adaptive radiations of higher plants and animals; and, (iii) interspecific horizontal DNA transfer encoding virtually all of the cellular functions between organisms and their viruses in all domains of life. Consequently, assuming that evolutionary processes occur in isolated genomes of individual species has become an unrealistic abstraction. Adaptive variations also involved natural genetic engineering of mobile DNA elements to rewire regulatory networks. In the most highly evolved organisms, biological complexity scales with "non-coding" DNA content more closely than with protein-coding capacity. Coincidentally, we have learned how so-called "non-coding" RNAs that are rich in repetitive mobile DNA sequences are key regulators of complex phenotypes. Both biotic and abiotic ecological challenges serve as triggers for episodes of elevated genome change. The intersections of cell activities, biosphere interactions, horizontal DNA transfers, and non-random Read-Write genome modifications by natural genetic engineering provide a rich molecular and biological foundation for understanding how ecological disruptions can stimulate productive, often abrupt, evolutionary transformations.

Hybrids and horizontal transfer: introgression allows adaptive allele discovery

Evolution has devised countless remarkable solutions to diverse challenges. Understanding the mechanistic basis of these solutions provides insights into how biological systems can be subtly tweaked without maladaptive consequences. The knowledge gained from illuminating these mechanisms is equally important to our understanding of fundamental evolutionary mechanisms as it is to our hopes of developing truly rational plant breeding and synthetic biology. In particular, modern population genomic approaches are proving very powerful in the detection of candidate alleles for mediating consequential adaptations that can be tested functionally. Especially striking are signals gained from contexts involving genetic transfers between populations, closely related species, or indeed between kingdoms. Here we discuss two major classes of these scenarios, adaptive introgression and horizontal gene flow, illustrating discoveries made across kingdoms.

Phytogeographic patterns and cryptic diversity in an aposematic toad from NW Argentina

The Yungas Redbelly Toad, Melanophryniscus rubriventris, is patchily distributed in Argentina, confined to the upland portion (1000-2000m above sea level) of the montane forests of northern and central regions of Salta, and in central-eastern and south-eastern Jujuy. This species is known for its striking aposematic color variation across its geographic distribution, and was once treated as a complex of three subspecies based on distinctive color patterns. Here we assess the geographical genetic variation within M. rubriventris and quantify divergence in color and pattern among individuals sampled from Northwestern Argentina. We compare multi-gene phylogeography of M. rubriventris to patterns of dorsal and ventral coloration to test whether evolutionary affinities predict variation in warning color. Our results reveal two well-supported species lineages: one confined to the extreme northern portion of our sampling area, and the other extending over most of the Argentine portion of the species' range, within which there are two populations. However, these well-supported evolutionary relationships do not mirror the marked variation in warning coloration. This discordance between DNA genealogy and warning color variation may reflect selection brought about by differences in local predation pressures, potentially coupled with effects of sexual selection and thermoregulation.

Phylogenetic information content of Copepoda ribosomal DNA repeat units: ITS1 and ITS2 impact

The utility of various regions of the ribosomal repeat unit for phylogenetic analysis was examined in 16 species representing four families, nine genera, and two orders of the subclass Copepoda (Crustacea). Fragments approximately 2000 bp in length containing the ribosomal DNA (rDNA) 18S and 28S gene fragments, the 5.8S gene, and the internal transcribed spacer regions I and II (ITS1 and ITS2) were amplified and analyzed. The DAMBE (Data Analysis in Molecular Biology and Evolution) software was used to analyze the saturation of nucleotide substitutions; this test revealed the suitability of both the 28S gene fragment and the ITS1/ITS2 rDNA regions for the reconstruction of phylogenetic trees. Distance (minimum evolution) and probabilistic (maximum likelihood, Bayesian) analyses of the data revealed that the 28S rDNA and the ITS1 and ITS2 regions are informative markers for inferring phylogenetic relationships among families of copepods and within the Cyclopidae family and associated genera. Split-graph analysis of concatenated ITS1/ITS2 rDNA regions of cyclopoid copepods suggested that the Mesocyclops, Thermocyclops, and Macrocyclops genera share complex evolutionary relationships. This study revealed that the ITS1 and ITS2 regions potentially represent different phylogenetic signals.

Phylogeography of Heliconius cydno and its closest relatives: disentangling their origin and diversification

The origins of the extraordinary diversity within the Neotropics have long fascinated biologists and naturalists. Yet, the underlying factors that have given rise to this diversity remain controversial. To test the relative importance of Quaternary climatic change and Neogene tectonic and paleogeographic reorganizations in the generation of biodiversity, we examine intraspecific variation across the Heliconius cydno radiation and compare this variation to that within the closely related Heliconius melpomene and Heliconius timareta radiations. Our data, which consist of both mtDNA and genome-scan data from nearly 2250 amplified fragment length polymorphism (AFLP) loci, reveal a complex history of differentiation and admixture at different geographic scales. Both mtDNA and AFLP phylogenies suggest that H. timareta and H. cydno are probably geographic extremes of the same radiation that probably diverged from H. melpomene prior to the Pliocene-Pleistocene boundary, consistent with hypotheses of diversification that rely on geological events in the Pliocene. The mtDNA suggests that this radiation originated in Central America or the northwestern region of South America, with a subsequent colonization of the eastern and western slopes of the Andes. Our genome-scan data indicate significant admixture among sympatric H. cydno/H. timareta and H. melpomene populations across the extensive geographic ranges of the two radiations. Within H. cydno, both mtDNA and AFLP data indicate significant population structure at local scales, with strong genetic differences even among adjacent H. cydno colour pattern races. These genetic patterns highlight the importance of past geoclimatic events, intraspecific gene flow, and local population differentiation in the origin and establishment of new adaptive forms.

Do Heliconius butterfly species exchange mimicry alleles?

Hybridization has the potential to transfer beneficial alleles across species boundaries, and there are a growing number of examples in which this has apparently occurred. Recent studies suggest that Heliconius butterflies have transferred wing pattern mimicry alleles between species via hybridization, but ancestral polymorphism could also produce a signature of shared ancestry around mimicry genes. To distinguish between these alternative hypotheses, we measured DNA sequence divergence around putatively introgressed mimicry loci and compared this with the rest of the genome. Our results reveal that putatively introgressed regions show strongly reduced sequence divergence between co-mimetic species, suggesting that their divergence times are younger than the rest of the genome. This is consistent with introgression and not ancestral variation. We further show that this signature of introgression occurs at sites throughout the genome, not just around mimicry genes.