Ecological specialization is associated with genetic structure in the ant-associated butterfly family Lycaenidae

To and fro in the archipelago: Repeated inter-island dispersal and New Guinea's orogeny affect diversification of Delias, the world's largest butterfly genus

The world's largest butterfly genus Delias, commonly known as Jezebels, comprises ca. 251 species found throughout Asia, Australia, and Melanesia. Most species are endemic to islands in the Indo-Australian Archipelago or to New Guinea and nearby islands in Melanesia, and many species are restricted to montane habitats over 1200 m. We inferred an extensively sampled and well-supported molecular phylogeny of the group to better understand the spatial and temporal dimensions of its diversification. The remarkable diversity of Delias evolved in just ca. 15-16 Myr (crown age). The most recent common ancestor of a clade with most of the species dispersed out of New Guinea ca. 14 Mya, but at least six subsequently diverging lineages dispersed back to the island. Diversification was associated with frequent dispersal of lineages among the islands of the Indo-Australian Archipelago, and the divergence of sister taxa on a single landmass was rare and occurred only on the largest islands, most notably on New Guinea. We conclude that frequent inter-island dispersal during the Neogene-likely facilitated by frequent sea level change-sparked much diversification during that period. Many extant New Guinea lineages started diversifying 5 Mya, suggesting that orogeny facilitated their diversification. Our results largely agree with the most recently proposed species group classification system, and we use our large taxon sample to extend this system to all described species. Finally, we summarize recent insights to speculate how wing pattern evolution, mimicry, and sexual selection might also contribute to these butterflies' rapid speciation and diversification.

Rapid radiation of ant parasitic butterflies during the Miocene aridification of Africa

Africa has undergone a progressive aridification during the last 20 My that presumably impacted organisms and fostered the evolution of life history adaptations. We test the hypothesis that shift to living in ant nests and feeding on ant brood by larvae of phyto-predaceous Lepidochrysops butterflies was an adaptive response to the aridification of Africa that facilitated the subsequent radiation of butterflies in this genus. Using anchored hybrid enrichment we constructed a time-calibrated phylogeny for Lepidochrysops and its closest, non-parasitic relatives in the Euchrysops section (Poloyommatini). We estimated ancestral areas across the phylogeny with process-based biogeographical models and diversification rates relying on time-variable and clade-heterogeneous birth-death models. The Euchrysops section originated with the emerging Miombo woodlands about 22 million years ago (Mya) and spread to drier biomes as they became available in the late Miocene. The diversification of the non-parasitic lineages decreased as aridification intensified around 10 Mya, culminating in diversity decline. In contrast, the diversification of the phyto-predaceous Lepidochrysops lineage proceeded rapidly from about 6.5 Mya when this unusual life history likely first evolved. The Miombo woodlands were the cradle for diversification of the Euchrysops section, and our findings are consistent with the hypothesis that aridification during the Miocene selected for a phyto-predaceous life history in species of Lepidochrysops, with ant nests likely providing caterpillars a safe refuge from fire and a source of food when vegetation was scarce.

Comparative mitogenomic and evolutionary analysis of Lycaenidae (Insecta: Lepidoptera): Potential association with high-altitude adaptation

Harsh environments (e.g., hypoxia and cold temperatures) of the Qinghai-Tibetan Plateau have a substantial influence on adaptive evolution in various species. Some species in Lycaenidae, a large and widely distributed family of butterflies, are adapted to the Qinghai-Tibetan Plateau. Here, we sequenced four mitogenomes of two lycaenid species in the Qinghai-Tibetan Plateau and performed a detailed comparative mitogenomic analysis including nine other lycaenid mitogenomes (nine species) to explore the molecular basis of high-altitude adaptation. Based on mitogenomic data, Bayesian inference, and maximum likelihood methods, we recovered a lycaenid phylogeny of [Curetinae + (Aphnaeinae + (Lycaeninae + (Theclinae + Polyommatinae)))]. The gene content, gene arrangement, base composition, codon usage, and transfer RNA genes (sequence and structure) were highly conserved within Lycaenidae. TrnS1 not only lacked the dihydrouridine arm but also showed anticodon and copy number diversity. The ratios of non-synonymous substitutions to synonymous substitutions of 13 protein-coding genes (PCGs) were less than 1.0, indicating that all PCGs evolved under purifying selection. However, signals of positive selection were detected in cox1 in the two Qinghai-Tibetan Plateau lycaenid species, indicating that this gene may be associated with high-altitude adaptation. Three large non-coding regions, i.e., rrnS-trnM (control region), trnQ-nad2, and trnS2-nad1, were found in the mitogenomes of all lycaenid species. Conserved motifs in three non-coding regions (trnE-trnF, trnS1-trnE, and trnP-nad6) and long sequences in two non-coding regions (nad6-cob and cob-trnS2) were detected in the Qinghai-Tibetan Plateau lycaenid species, suggesting that these non-coding regions were involved in high-altitude adaptation. In addition to the characterization of Lycaenidae mitogenomes, this study highlights the importance of both PCGs and non-coding regions in high-altitude adaptation.

Phylogeographic structure of the dunes sagebrush lizard, an endemic habitat specialist

Phylogeographic divergence and population genetic diversity within species reflect the impacts of habitat connectivity, demographics, and landscape level processes in both the recent and distant past. Characterizing patterns of differentiation across the geographic range of a species provides insight on the roles of organismal and environmental traits in evolutionary divergence and future population persistence. This is particularly true of habitat specialists where habitat availability and resource dependence may result in pronounced genetic structure as well as increased population vulnerability. We use DNA sequence data as well as microsatellite genotypes to estimate range-wide phylogeographic divergence, historical population connectivity, and historical demographics in an endemic habitat specialist, the dunes sagebrush lizard (Sceloporus arenicolus). This species is found exclusively in dune blowouts and patches of open sand within the shinnery oak-sand dune ecosystem of southeastern New Mexico and adjacent Texas. We find evidence of phylogeographic structure consistent with breaks and constrictions in suitable habitat at the range-wide scale. In addition, we find support for a dynamic and variable evolutionary history across the range of S. arenicolus. Populations in the Monahans Sandhills have deeply divergent lineages consistent with long-term demographic stability. In contrast, populations in the Mescalero Sands are not highly differentiated, though we do find evidence of demographic expansion in some regions and relative demographic stability in others. Phylogeographic history and population genetic differentiation in this species has been shaped by the configuration of habitat patches within a geologically complex and historically dynamic landscape. Our findings identify regions as genetically distinctive conservation units as well as underscore the genetic and demographic history of different lineages of S. arenicolus.

Combining stable isotope analysis with DNA metabarcoding improves inferences of trophic ecology

Knowing what animals eat is fundamental to our ability to understand and manage biodiversity and ecosystems, but researchers often must rely on indirect methods to infer trophic position and food intake. Using an approach that combines evidence from stable isotope analysis and DNA metabarcoding, we assessed the diet and trophic position of Anthene usamba butterflies, for which there are no known direct observations of larval feeding. An earlier study that analyzed adults rather than caterpillars of A. usamba inferred that this butterfly was aphytophagous, but we found that the larval guts of A. usamba and two known herbivorous lycaenid species contain chloroplast 16S sequences. Moreover, chloroplast barcoding revealed high sequence similarity between chloroplasts found in A. usamba guts and the chloroplasts of the Vachellia drepanolobium trees on which the caterpillars live. Stable isotope analysis provided further evidence that A. usamba caterpillars feed on V. drepanolobium, and the possibilities of strict herbivory versus limited omnivory in this species are discussed. These results highlight the importance of combining multiple approaches and considering ontogeny when using stable isotopes to infer trophic ecology where direct observations are difficult or impossible.