How the Easter Egg Weevils Got Their Spots: Phylogenomics Reveals Müllerian Mimicry in Pachyrhynchus (Coleoptera, Curculionidae).

The evolutionary origins of mimicry in the Easter egg weevil, Pachyrhynchus, have fascinated researchers since first noted more than a century ago by Alfred Russel Wallace. Müllerian mimicry, or mimicry in which 2 or more distasteful species look similar, is widespread throughout the animal kingdom. Given the varied but discrete color patterns in Pachyrhynchus, this genus presents one of the best opportunities to study the evolution of both perfect and imperfect mimicry. We analyzed more than 10,000 UCE loci using a novel partitioning strategy to resolve the relationships of closely related species in the genus. Our results indicate that many of the mimetic color patterns observed in sympatric species are due to convergent evolution. We suggest that this convergence is driven by positive frequency-dependent selection. [Biogeography, discrete traits, frequency-dependent selection, mimicry, partitioning, Philippines, polymorphic, UCE.].

The Easter Egg Weevil (Pachyrhynchus) genome reveals syntenic patterns in Coleoptera across 200 million years of evolution.

Patterns of genomic architecture across insects remain largely undocumented or decoupled from a broader phylogenetic context. For instance, it is unknown whether translocation rates differ between insect orders. We address broad scale patterns of genome architecture across Insecta by examining synteny in a phylogenetic framework from open-source insect genomes. To accomplish this, we add a chromosome level genome to a crucial lineage, Coleoptera. Our assembly of the Pachyrhynchus sulphureomaculatus genome is the first chromosome scale genome for the hyperdiverse Phytophaga lineage and currently the largest insect genome assembled to this scale. The genome is significantly larger than those of other weevils, and this increase in size is caused by repetitive elements. Our results also indicate that, among beetles, there are instances of long-lasting (>200 Ma) localization of genes to a particular chromosome with few translocation events. While some chromosomes have a paucity of translocations, intra-chromosomal synteny was almost absent, with gene order thoroughly shuffled along a chromosome. This large amount of reshuffling within chromosomes with few inter-chromosomal events contrasts with patterns seen in mammals in which the chromosomes tend to exchange larger blocks of material more readily. To place our findings in an evolutionary context, we compared syntenic patterns across Insecta in a phylogenetic framework. For the first time, we find that synteny decays at an exponential rate relative to phylogenetic distance. Additionally, there are significant differences in decay rates between insect orders, this pattern was not driven by Lepidoptera alone which has a substantially different rate.

High-Quality, Chromosome-Level Reference Genomes of the Viviparous Caribbean Skinks Spondylurus nitidus and S. culebrae.

New World mabuyine skinks are a diverse radiation of morphologically cryptic lizards with unique reproductive biologies. Recent studies examining population-level data (morphological, ecological, and genomic) have uncovered novel biodiversity and phenotypes, including the description of dozens of new species and insights into the evolution of their highly complex placental structures. Beyond the potential for this diverse group to serve as a model for the evolution of viviparity in lizards, much of the taxonomic diversity is concentrated in regions experiencing increasing environmental instability from climate and anthropogenic change. Consequently, a better understanding of genome structure and diversity will be an important tool in the adaptive management and conservation of this group. Skinks endemic to Caribbean islands are particularly vulnerable to global change with several species already considered likely extinct and several remaining species either endangered or threatened. Combining PacBio long-read sequencing, Hi-C, and RNAseq data, here we present the first genomic resources for this group by describing new chromosome-level reference genomes for the Puerto Rican Skink Spondylurus nitidus and the Culebra Skink S. culebrae. Results indicate two high quality genomes, both ∼1.4 Gb, assembled nearly telomere to telomere with complete mitochondrion assembly and annotation.

Genome Assembly of Pyrocephalus nanus: A Step Toward the Genetic Conservation of the Endangered Little Vermilion Flycatcher of the Galapagos Islands.

Incredibly powerful whole genome studies of conservation genetics, evolution, and biogeography become possible for non-model organisms when reference genomes are available. Here, we report the sequence and assembly of the whole genome of the little vermilion flycatcher (Pyrocephalus nanus; family Tyrannidae), which is an endemic, endangered, and declining species of the Galapagos Islands. Using PacBio HiFi reads to assemble long contigs and Hi-C reads for scaffolding, we assembled a genome of 1.07 Gb comprising 267 contigs in 152 scaffolds, scaffold N50 74 M, contig N50 17.8 M, with 98.9% assigned to candidate chromosomal sequences and 99.72% of the BUSCO passeriformes 10,844 single-copy orthologs present. In addition, we used the novel HiFiMiTie pipeline to fully assemble and verify all portions of the mitochondrial genome from HiFi reads, obtaining a mitogenome of 17,151 bases, containing 13 protein-coding genes, 22 tRNAs, 2 rRNAs, two control regions, and a unique structure of control region duplication and repeats. These genomes will be a critical tool for much-needed studies of phylogenetics, population genetics, biogeography, and conservation genetics of Pyrocephalus and related genera. This genome and other studies that use it will be able to provide recommendations for conservation management, taxonomic improvement, and to understand the evolution and diversification of this genus within the Galapagos Islands.

Biogeography confounds the signal of cospeciation in Batesian mimicry.

Since the inception of the field of evolution, mimicry has yielded insights into foundational evolutionary processes, including adaptive peak shifts, speciation, and the emergence and maintenance of phenotypic polymorphisms.1,2,3 In recent years, the coevolutionary processes generating mimicry have gained increasing attention from researchers. Despite significant advances in understanding Batesian and Müllerian mimicry in Lepidopteran systems, few other mimetic systems have received similar detailed research. Here, we present a Batesian mimicry complex involving flightless, armored Pachyrhynchus weevils and their winged Doliops longhorn beetle mimics and examine their coevolutionary patterns within the Philippine archipelagos. Pachyrhynchus weevils are primarily found in the Philippines, where distinct species radiations have occurred on different islands, each with unique color patterns serving as a warning to predators. This defensive trait and mimicry between unrelated species were first described by Wallace in 1889. Notably, the distantly related longhorn beetle Doliops, despite being soft-bodied and ostensibly palatable, mimics the heavily armored, flightless Pachyrhynchus. To address mimicry in this system, we reconstructed the phylogeny of Doliops using a probe set consisting of 38,000 ultraconserved elements. Our study examines the following questions central to understanding the Pachyrhynchus-Doliops mimicry system: (1) to what extent are coevolutionary interactions conserved (i.e., lineage-constrained) and (2) are the codiversification patterns primarily driven by biotic or abiotic factors?4,5,6 To assess color mimicry and cospeciation, we examined the evolution of nanostructure-based warning colors and the effect of island biogeography on cospeciation. Our findings demonstrate the beetle's ability to repeatedly evolve multiple solutions to similar evolutionary challenges, evolving similar color patterns using different types of photonic crystals with varying degrees of order. We revealed that the observed pattern of cospeciation is driven mainly by abiotic factors from their biogeographic history. Unlike the patterns of coevolution seen between angiosperms and insect lineages,7 most ecological interactions do not persist longer than a few million years, leading to patterns of modularity rather than ecological nestedness.4,6,7.

Phylogenomics of Tetraopes longhorn beetles unravels their evolutionary history and biogeographic origins.

Tetraopes longhorn beetles are known for their resistance to milkweed plant toxins and their coevolutionary dynamics with milkweed plants (Asclepias). This association is considered a textbook example of coevolution, in which each species of Tetraopes is specialized to feed on one or a few species of Asclepias. A major challenge to investigating coevolutionary hypotheses and conducting molecular ecology studies lies in the limited understanding of the evolutionary history and biogeographical patterns of Tetraopes. By integrating genomic, morphological, paleontological, and geographical data, we present a robust phylogeny of Tetraopes and their relatives, using three inference methods with varying subsets of data, encompassing 2-12 thousand UCE loci. We elucidate the diversification patterns of Tetraopes species across major biogeographical regions and their colonization of the American continent. Our findings suggest that the genus originated in Central America approximately 21 million years ago during the Miocene and diversified from the Mid-Miocene to the Pleistocene. These events coincided with intense geological activity in Central America. Additionally, independent colonization events in North America occurred from the Late Miocene to the early Pleistocene, potentially contributing to the early diversification of the group. Our data suggest that a common ancestor of Tetraopini migrated into North America, likely facilitated by North Atlantic land bridges, while closely related tribes diverged in Asia and Europe during the Paleocene. Establishing a robust and densely sampled phylogeny of Tetraopes beetles provides a foundation for investigating micro- and macroevolutionary phenomena, including clinal variation, coevolution, and detoxification mechanisms in this ecologically important group.

A range-wide genetic bottleneck overwhelms contemporary landscape factors and local abundance in shaping genetic patterns of an alpine butterfly (Lepidoptera: Pieridae: Colias behrii).

Spatial and environmental heterogeneity are major factors in structuring species distributions in alpine landscapes. These landscapes have also been affected by glacial advances and retreats, causing alpine taxa to undergo range shifts and demographic changes. These nonequilibrium population dynamics have the potential to obscure the effects of environmental factors on the distribution of genetic variation. Here, we investigate how demographic change and environmental factors influence genetic variation in the alpine butterfly Colias behrii. Data from 14 microsatellite loci provide evidence of bottlenecks in all population samples. We test several alternative models of demography using approximate Bayesian computation (ABC), with the results favouring a model in which a recent bottleneck precedes rapid population growth. Applying independent calibrations to microsatellite loci and a nuclear gene, we estimate that this bottleneck affected both northern and southern populations 531-281 years ago, coinciding with a period of global cooling. Using regression approaches, we attempt to separate the effects of population structure, geographical distance and landscape on patterns of population genetic differentiation. Only 40% of the variation in F(ST) is explained by these models, with geographical distance and least-cost distance among meadow patches selected as the best predictors. Various measures of genetic diversity within populations are also decoupled from estimates of local abundance and habitat patch characteristics. Our results demonstrate that demographic change can have a disproportionate influence on genetic diversity in alpine species, contrasting with other studies that suggest landscape features control contemporary demographic processes in high-elevation environments.

Metagenomic clustering reveals microbial contamination as an essential consideration in ultraconserved element design for phylogenomics with insect museum specimens.

Phylogenomics via ultraconserved elements (UCEs) has led to improved phylogenetic reconstructions across the tree of life. However, inadvertently incorporating non-targeted DNA into the UCE marker design will lead to misinformation being incorporated into subsequent analyses. To date, the effectiveness of basic metagenomic filtering strategies has not been assessed in arthropods. Designing markers from museum specimens requires careful consideration of methods due to the high levels of microbial contamination typically found in such specimens. We investigate if contaminant sequences are carried forward into a UCE marker set we developed from insect museum specimens using a standard bioinformatics pipeline. We find that the methods currently employed by most researchers do not exclude contamination from the final set of targets. Lastly, we highlight several paths forward for reducing contamination in UCE marker design.

An Annotated Chromosome-Level Reference Genome of the Red-Eared Slider Turtle (Trachemys scripta elegans).

Among vertebrates, turtles have many unique characteristics providing biologists with opportunities to study novel evolutionary innovations and processes. We present here a high-quality, partially phased, and chromosome-level Red-Eared Slider (Trachemys scripta elegans, TSE) genome as a reference for future research on turtle and tetrapod evolution. This TSE assembly is 2.269 Gb in length, has one of the highest scaffold N50 and N90 values of any published turtle genome to date (N50 = 129.68 Mb and N90 = 19 Mb), and has a total of 28,415 annotated genes. We introduce synteny analyses using BUSCO single-copy orthologs, which reveal two chromosome fusion events accounting for differences in chromosome counts between emydids and other cryptodire turtles and reveal many fission/fusion events for birds, crocodiles, and snakes relative to TSE. This annotated chromosome-level genome will provide an important reference genome for future studies on turtle, vertebrate, and chromosome evolution.

Correction: Ultraconserved elements (UCEs) resolve the phylogeny of Australasian smurf-weevils.

[This corrects the article DOI: 10.1371/journal.pone.0188044.].

Ultraconserved elements (UCEs) resolve the phylogeny of Australasian smurf-weevils.

Weevils (Curculionoidea) comprise one of the most diverse groups of organisms on earth. There is hardly a vascular plant or plant part without its own species of weevil feeding on it and weevil species diversity is greater than the number of fishes, birds, reptiles, amphibians and mammals combined. Here, we employ ultraconserved elements (UCEs) designed for beetles and a novel partitioning strategy of loci to help resolve phylogenetic relationships within the radiation of Australasian smurf-weevils (Eupholini). Despite being emblematic of the New Guinea fauna, no previous phylogenetic studies have been conducted on the Eupholini. In addition to a comprehensive collection of fresh specimens, we supplement our taxon sampling with museum specimens, and this study is the first target enrichment phylogenomic dataset incorporating beetle specimens from museum collections. We use both concatenated and species tree analyses to examine the relationships and taxonomy of this group. For species tree analyses we present a novel partitioning strategy to better model the molecular evolutionary process in UCEs. We found that the current taxonomy is problematic, largely grouping species on the basis of similar color patterns. Finally, our results show that most loci required multiple partitions for nucleotide rate substitution, suggesting that single partitions may not be the optimal partitioning strategy to accommodate rate heterogeneity for UCE loci.

Permanent Genetic Resources added to Molecular Ecology Resources Database 1 October 2009-30 November 2009.

This article documents the addition of 411 microsatellite marker loci and 15 pairs of Single Nucleotide Polymorphism (SNP) sequencing primers to the Molecular Ecology Resources Database. Loci were developed for the following species: Acanthopagrus schlegeli, Anopheles lesteri, Aspergillus clavatus, Aspergillus flavus, Aspergillus fumigatus, Aspergillus oryzae, Aspergillus terreus, Branchiostoma japonicum, Branchiostoma belcheri, Colias behrii, Coryphopterus personatus, Cynogolssus semilaevis, Cynoglossus semilaevis, Dendrobium officinale, Dendrobium officinale, Dysoxylum malabaricum, Metrioptera roeselii, Myrmeciza exsul, Ochotona thibetana, Neosartorya fischeri, Nothofagus pumilio, Onychodactylus fischeri, Phoenicopterus roseus, Salvia officinalis L., Scylla paramamosain, Silene latifo, Sula sula, and Vulpes vulpes. These loci were cross-tested on the following species: Aspergillus giganteus, Colias pelidne, Colias interior, Colias meadii, Colias eurytheme, Coryphopterus lipernes, Coryphopterus glaucofrenum, Coryphopterus eidolon, Gnatholepis thompsoni, Elacatinus evelynae, Dendrobium loddigesii Dendrobium devonianum, Dysoxylum binectariferum, Nothofagus antarctica, Nothofagus dombeyii, Nothofagus nervosa, Nothofagus obliqua, Sula nebouxii, and Sula variegata. This article also documents the addition of 39 sequencing primer pairs and 15 allele specific primers or probes for Paralithodes camtschaticus.