Historical and contemporary factors generate unique butterfly communities on islands

The Intriguing Biogeographic Pattern of the Italian Wall Lizard Podarcis siculus (Squamata: Lacertidae) in the Tuscan Archipelago Reveals the Existence of a New Ancient Insular Clade

The Tuscan Archipelago is one of the most ancient and ecologically heterogeneous island systems in the Mediterranean. The biodiversity of these islands was strongly shaped by the Pliocene and Pleistocene sea regressions and transgression, resulting in different waves of colonization and isolation of species coming from the mainland. The Italian wall lizard, Podarcis siculus, is present on the following islands of the Tuscan Archipelago: Elba, Giglio, Giannutri, Capraia, Montecristo and Cerboli. The species in the area displays a relatively high morphological variability that in the past led to the description of several subspecies. In this study, both the genetic and morphological diversity of P. siculus of the Tuscan Archipelago were investigated. Specifically, the meristic characters and the dorsal pattern were analyzed, while the genetic relationships among these populations were explored with mtDNA and microsatellite nuclear markers to reconstruct the colonization history of the Archipelago. Our results converge in the identification of at least two different waves of colonization in the Archipelago: Elba, and the populations of Cerboli and Montecristo probably originate from historical introductions from mainland Tuscany, while those of Giglio and Capraia are surviving populations of an ancient lineage which colonized the Tuscan Archipelago during the Pliocene and which shares a common ancestry with the P. siculus populations of south-eastern Italy. Giannutri perhaps represents an interesting case of hybridization between the populations from mainland Tuscany and the Giglio-Capraia clade. Based on the high phenotypic and molecular distinctiveness of this ancient clade, these populations should be treated as distinct units deserving conservation and management efforts as well as further investigation to assess their taxonomic status.

Population genetics of a recent range expansion and subsequent loss of migration in monarch butterflies

Range expansions—whether permanent or transient—strongly influence the distribution of genetic variation in space. Monarch butterflies are best known for long‐distance seasonal migration within North America but are also established as nonmigratory populations around the world, including on Pacific Islands. Previous research has highlighted stepwise expansion across the Pacific, though questions remain about expansion timing and the population genetic consequences of migration loss. Here, we present reduced‐representation sequencing data for 275 monarchs from North America (n = 85), 12 Pacific Islands (n = 136) and three locations in Australia (n = 54), with the goal of understanding (i) how the monarch's Pacific expansion has shaped patterns of population genetic variation and (ii) how loss of migration has influenced spatial patterns of differentiation. We find support for previously described stepwise dispersal across the Pacific and document an additional expansion from Hawaii into the Mariana Islands. Nonmigratory monarchs within the Mariana Islands show strong patterns of differentiation, despite their proximity; by contrast, migratory North American samples form a single genetically panmictic population across the continent. Estimates of Pacific establishment timing are highly uncertain (~100–1,000,000 years ago) but overlap with historical records that indicate a recent expansion. Our data support (i) a recent expansion across the Pacific whose timing overlaps with available historical records of establishment and (ii) a strong role for seasonal migration in determining patterns of spatial genetic variation. Our results are noteworthy because they demonstrate how the evolution of partial migration can drive population differentiation over contemporary timescales.

Local adaptation and migratory habits balance spatial-genetic structure between continental and insular chestnut tiger butterflies in East Asia

Geographic and climatic differences between islands and continents may affect the evolution of their biota, and promote divergent selection in species distributed in both landscapes. To assess spatial-genetic structure, we genotyped 18 expressed sequence tag-simple sequence repeat (EST-SSR) loci and sequenced two mtDNA markers (ND5 and COI) and one nuclear marker (EF1α) in two subspecies of the butterfly Parantica sita. Compared with nuclear markers, mtDNA had a stronger signal of population structure. Approximate Bayesian computation (ABC) suggested that a continuous-gene-flow model best described the data. According to this model, the two subspecies diverged approximately 23.1 kya, with ten times more introgression from the continental (ssp. sita) to the insular subspecies (ssp. niphonica) than vice versa. Ecological niche modeling was performed to predict the paleo- and current potential distributions and elucidate the geohistorical process, which revealed a northeastern, insular origin. Winter precipitation and annual temperature range were the main determinants of the subspecies distributions. Maximum-likelihood population-effects models showed that the population differentiation of the insular and continental subspecies was primarily affected by environmental resistance and local climate. Sex-biased migration capacity and long-term precipitation-driven divergence between the continental and insular lineages shaped the current genetic structure of P. sita. Evidence from the nuclear markers confirmed inter-subspecific gene flow despite adaptive divergence between the subspecies. These results imply that the continental subspecies is still capable of returning to the island and introgressing with the insular subspecies.

A DNA barcode library for the butterflies of North America

Although the butterflies of North America have received considerable taxonomic attention, overlooked species and instances of hybridization continue to be revealed. The present study assembles a DNA barcode reference library for this fauna to identify groups whose patterns of sequence variation suggest the need for further taxonomic study. Based on 14,626 records from 814 species, DNA barcodes were obtained for 96% of the fauna. The maximum intraspecific distance averaged 1/4 the minimum distance to the nearest neighbor, producing a barcode gap in 76% of the species. Most species (80%) were monophyletic, the others were para- or polyphyletic. Although 15% of currently recognized species shared barcodes, the incidence of such taxa was far higher in regions exposed to Pleistocene glaciations than in those that were ice-free. Nearly 10% of species displayed high intraspecific variation (>2.5%), suggesting the need for further investigation to assess potential cryptic diversity. Aside from aiding the identification of all life stages of North American butterflies, the reference library has provided new perspectives on the incidence of both cryptic and potentially over-split species, setting the stage for future studies that can further explore the evolutionary dynamics of this group.

Wolbachia affects mitochondrial population structure in two systems of closely related Palaearctic blue butterflies

The bacterium Wolbachia infects many insect species and spreads by diverse vertical and horizontal means. As co-inherited organisms, these bacteria often cause problems in mitochondrial phylogeny inference. The phylogenetic relationships of many closely related Palaearctic blue butterflies (Lepidoptera: Lycaenidae: Polyommatinae) are ambiguous. We considered the patterns of Wolbachia infection and mitochondrial diversity in two systems: Aricia agestis/Aricia artaxerxes and the Pseudophilotes baton species complex. We sampled butterflies across their distribution ranges and sequenced one butterfly mitochondrial gene and two Wolbachia genes. Both butterfly systems had uninfected and infected populations, and harboured several Wolbachia strains. Wolbachia was highly prevalent in A. artaxerxes and the host's mitochondrial structure was shallow, in contrast to A. agestis. Similar bacterial alleles infected both Aricia species from nearby sites, pointing to a possible horizontal transfer. Mitochondrial history of the P. baton species complex mirrored its Wolbachia infection and not the taxonomical division. Pseudophilotes baton and P. vicrama formed a hybrid zone in Europe. Wolbachia could obscure mitochondrial history, but knowledge on the infection helps us to understand the observed patterns. Testing for Wolbachia should be routine in mitochondrial DNA studies.

How long is 3 km for a butterfly? Ecological constraints and functional traits explain high mitochondrial genetic diversity between Sicily and the Italian Peninsula

Populations inhabiting Mediterranean islands often show contrasting genetic lineages, even on islands that were connected to the mainland during glacial maxima. This pattern is generated by forces acting in historical and contemporary times. Understanding these phenomena requires comparative studies integrating genetic structure, functional traits and dispersal constraints. Using as a model the butterfly species living across the Messina strait (3 km wide) separating Sicily from the Italian Peninsula, we aimed to unravel the mechanisms limiting the dispersal of matrilines and generating genetic differentiation across a narrow sea strait. We analysed the mitochondrial COI gene of 84 butterfly species out of 90 documented in Sicily and compared them with populations from the neighbouring southern Italian Peninsula (1,398 sequences) and from the entire Palaearctic region (8,093 sequences). For each species, we regressed 13 functional traits and 2 ecological constraints to dispersal (winds experienced at the strait and climatic suitability) against genetic differentiation between Sicily and Italian Peninsula to understand the factors limiting dispersal. More than a third of the species showed different haplogroups across the strait and most of them also represented endemic haplogroups for this island. One fifth of Sicilian populations (and 32.3% of endemic lineages) had their closest relatives in distant areas, instead of the neighbouring Italian Peninsula, which suggests high relictuality. Haplotype diversity was significantly explained by the length of the flight period, an intrinsic phenology trait, while genetic differentiation was explained by both intrinsic traits (wingspan and degree of generalism) and contemporary local constraints (winds experienced at the strait and climatic suitability). A relatively narrow sea strait can produce considerable differentiation among butterfly matrilines and this phenomenon showed a largely deterministic fingerprint. Because of unfavourable winds, populations of the less dispersive Sicilian butterflies tended to differentiate into endemic variants or to maintain relict populations. Understanding these phenomena required the integration of DNA sequences, species traits and physical constraints for a large taxon at continental scale. Future studies may reveal if the patterns here shown for mitochondrial DNA are also reflected in the nuclear genome or, alternatively, are the product of limited female dispersal.

Global invasion history of the agricultural pest butterfly Pieris rapae revealed with genomics and citizen science

The small cabbage white butterfly, Pieris rapae, is a major agricultural pest of cruciferous crops and has been introduced to every continent except South America and Antarctica as a result of human activities. In an effort to reconstruct the near-global invasion history of P. rapae, we developed a citizen science project, the "Pieris Project," and successfully amassed thousands of specimens from 32 countries worldwide. We then generated and analyzed nuclear (double-digest restriction site-associated DNA fragment procedure [ddRAD]) and mitochondrial DNA sequence data for these samples to reconstruct and compare different global invasion history scenarios. Our results bolster historical accounts of the global spread and timing of P. rapae introductions. We provide molecular evidence supporting the hypothesis that the ongoing divergence of the European and Asian subspecies of P. rapae (∼1,200 y B.P.) coincides with the diversification of brassicaceous crops and the development of human trade routes such as the Silk Route (Silk Road). The further spread of P. rapae over the last ∼160 y was facilitated by human movement and trade, resulting in an almost linear series of at least 4 founding events, with each introduced population going through a severe bottleneck and serving as the source for the next introduction. Management efforts of this agricultural pest may need to consider the current existence of multiple genetically distinct populations. Finally, the international success of the Pieris Project demonstrates the power of the public to aid scientists in collections-based research addressing important questions in invasion biology, and in ecology and evolutionary biology more broadly.

The conundrum of species delimitation: a genomic perspective on a mitogenetically super-variable butterfly

The Palaearctic butterfly Melitaea didyma stands out as one of the most striking cases of intraspecific genetic differentiation detected in Lepidoptera: 11 partially sympatric mitochondrial lineages have been reported, displaying levels of divergence of up to 7.4%. To better understand the evolutionary processes underlying the diversity observed in mtDNA, we compared mtDNA and genome-wide SNP data using double-digest restriction site-associated DNA sequencing (ddRADseq) results from 93 specimens of M. didyma ranging from Morocco to eastern Kazakhstan. We found that, between ddRADseq and mtDNA results, there is a match only in populations that probably remained allopatric for long periods of time. Other mtDNA lineages may have resulted from introgression events and were probably affected by Wolbachia infection. The five main ddRADseq clades supported by STRUCTURE were parapatric or allopatric and showed high pairwise F_ST values, but some were also estimated to display various levels of gene flow. Melitaea didyma represents one of the first cases of deep mtDNA splits among European butterflies assessed by a genome-wide DNA analysis and reveals that the interpretation of patterns remains challenging even when a high amount of genomic data is available. These findings actualize the ongoing debate of species delimitation in allopatry, an issue probably of relevance to a significant proportion of global biodiversity.

Use of genetic, climatic, and microbiological data to inform reintroduction of a regionally extinct butterfly

Species reintroductions are increasingly used as means of mitigating biodiversity loss. Besides habitat quality at the site targeted for reintroduction, the choice of source population can be critical for success. The butterfly Melanargia russiae (Esper´s marbled white) was extirpated from Hungary over 100 years ago, and a reintroduction program has recently been approved. We used museum specimens of this butterfly, mitochondrial DNA data (mtDNA), endosymbiont screening, and climatic-similarity analyses to determine which extant populations should be used for its reintroduction. The species displayed 2 main mtDNA lineages across its range: 1 restricted to Iberia and southern France (Iberian lineage) and another found throughout the rest of its range (Eurasian lineage). These 2 lineages possessed highly divergent wsp alleles of the bacterial endosymbiont Wolbachia. The century-old Hungarian specimens represented an endemic haplotype belonging to the Eurasian lineage, differing by one mutation from the Balkan and eastern European populations. The Hungarian populations of M. russiae occurred in areas with a colder and drier climate relative to most sites with extant known populations. Our results suggest the populations used for reintroduction to Hungary should belong to the Eurasian lineage, preferably from eastern Ukraine (genetically close and living in areas with the highest climatic similarity). Materials stored in museum collections can provide unique opportunities to document historical genetic diversity and help direct conservation.

DNA barcoding unravels contrasting evolutionary history of two widespread Asian tiger moth species during the Late Pleistocene

Populations of widespread pest insects in tropical areas are characterized by a complex evolutionary history, with overlapping natural and human-mediated dispersal events, sudden expansions, and bottlenecks. Here, we provide biogeographic reconstructions for two widespread pest species in the tiger moth genus Creatonotos (Lepidoptera: Erebidae: Arctiinae) based on the mitochondrial cytochrome c oxidase subunit I (COI) gene. The Asian Creatonotos transiens reveals shallow genetic divergence between distant populations that does not support its current intraspecific systematics with several local subspecies. In contrast, the more widespread Creatonotos gangis comprises at least three divergent subclades corresponding to certain geographic areas, i.e. Australia, Arabia + South Asia and Southeast Asia. With respect to our approximate Bayesian computation (ABC) model, the expansion of Creatonotos gangis into Australia is placed in the Late Pleistocene (~65–63 ka). This dating coincide with an approximate time of the earliest human migration into the continent (~65–54 ka) and the period of intervisibility between Timor and Australia (~65–62 ka). Our findings highlight that the drying Sunda and Sahul shelf areas likely support successful migrations of Asian taxa into Australia during the Pleistocene. The phylogeographic patterns discovered in this study can be used to improve the effectiveness of integrated pest control programs that is a task of substantial practical importance to a broad range of agricultural stakeholders.

Genomics of extreme ecological specialists: multiple convergent evolution but no genetic divergence between ecotypes of Maculinea alcon butterflies

Biotic interactions are often acknowledged as catalysers of genetic divergence and eventual explanation of processes driving species richness. We address the question, whether extreme ecological specialization is always associated with lineage sorting, by analysing polymorphisms in morphologically similar ecotypes of the myrmecophilous butterfly Maculinea alcon. The ecotypes occur in either hygric or xeric habitats, use different larval host plants and ant species, but no significant distinctive molecular traits have been revealed so far. We apply genome-wide RAD-sequencing to specimens originating from both habitats across Europe in order to get a view of the potential evolutionary processes at work. Our results confirm that genetic variation is mainly structured geographically but not ecologically - specimens from close localities are more related to each other than populations of each ecotype from distant localities. However, we found two loci for which the association with xeric versus hygric habitats is supported by segregating alleles, suggesting convergent evolution of habitat preference. Thus, ecological divergence between the forms probably does not represent an early stage of speciation, but may result from independent recurring adaptations involving few genes. We discuss the implications of these results for conservation and suggest preserving biotic interactions and main genetic clusters.

DNA Barcoding of an Assembly of Montane Andean Butterflies (Satyrinae): Geographical Scale and Identification Performance

DNA barcoding is a technique used primarily for the documentation and identification of biological diversity based on mitochondrial DNA sequences. Butterflies have received particular attention in DNA barcoding studies, although varied performance may be obtained due to different scales of geographic sampling and speciation processes in various groups. The montane Andean Satyrinae constitutes a challenging study group for taxonomy. The group displays high richness, with more of 550 species, and remarkable morphological similarity among taxa, which renders their identification difficult. In the present study, we evaluated the effectiveness of DNA barcodes in the identification of montane Andean satyrines and the effect of increased geographical scale of sampling on identification performance. Mitochondrial sequences were obtained from 104 specimens of 39 species and 16 genera, collected in a forest remnant in the northwest Andes. DNA barcoding has proved to be a useful tool for the identification of the specimens, with a well-defined gap and producing clusters with unambiguous identifications for all the morphospecies in the study area. The expansion of the geographical scale with published data increased genetic distances within species and reduced those among species, but did not generally reduce the success of specimen identification. Only in Forsterinaria rustica (Butler, 1868), a taxon with high intraspecific variation, the barcode gap was lost and low support for monophyly was obtained. Likewise, expanded sampling resulted in a substantial increase in the intraspecific distance in Morpho sulkowskyi (Kollar, 1850); Panyapedaliodes drymaea (Hewitson, 1858); Lymanopoda obsoleta (Westwood, 1851); and Lymanopoda labda Hewitson, 1861; but for these species, the barcode gap was maintained. These divergent lineages are nonetheless worth a detailed study of external and genitalic morphology variation, as well as ecological features, in order to determine the potential existence of cryptic species. Even including these cases, DNA barcoding performance in specimen identification was 100% successful based on monophyly, an unexpected result in such a taxonomically complicated group.

Disentangling the influence of environmental and anthropogenic factors on the distribution of endemic vascular plants in Sardinia

Due to the impelling urgency of plant conservation and the increasing availability of high resolution spatially interpolated (e.g. climate variables) and categorical data (e.g. land cover and vegetation type), many recent studies have examined relationships among plant species distributions and a diversified set of explanatory factors; nevertheless, global and regional patterns of endemic plant richness remain in many cases unexplained. One such pattern is the 294 endemic vascular plant taxa recorded on a 1 km resolution grid on the environmentally heterogeneous island of Sardinia. Sixteen predictors, including topographic, geological, climatic and anthropogenic factors, were used to model local (number of taxa inside each 1 km grid cell) Endemic Vascular Plant Richness (EVPR). Generalized Linear Models were used to evaluate how each factor affected the distribution of local EVPR. Significant relationships with local EVPR and topographic, geological, climatic and anthropogenic factors were found. In particular, elevation explained the larger fraction of variation in endemic richness but other environmental factors (e.g. precipitation seasonality and slope) and human-related factors (e.g. the Human Influence Index (HII) and the proportion of anthropogenic land uses) were, respectively, positively and negatively correlated with local EVPR. Regional EVPR (number of endemic taxa inside each 100 m elevation interval) was also measured to compare local and regional EVPR patterns along the elevation gradient. In contrast to local, regional EVPR tended to decrease with altitude partly due to the decreasing area covered along altitude. The contrasting results between local and regional patterns suggest that local richness increases as a result of increased interspecific aggregation along altitude, whereas regional richness may depend on the interaction between area and altitude. This suggests that the shape and magnitude of the species-area relationship might vary with elevation. This work provides-for the first time in Sardinia-a comprehensive analysis of the influence of environmental factors on the pattern of EVPR in the entire territory, from sea level to the highest peaks. Elevation, as well as other environmental and human-related variables, were confirmed to be influencing factors. In addition, variations of EVPR patterns at regional-to-local spatial scales inspire next investigations on the possible interaction between elevation and area in explaining patterns of plant species richness.