Comprehensive phylogeny of Pieridae butterflies reveals strong correlation between diversification and temperature

Temperature is thought to be a key factor influencing global species richness patterns. We investigate the link between temperature and diversification in the butterfly family Pieridae by combining next generation DNA sequences and published molecular data with fine-grained distribution data. We sampled nearly 600 pierid butterfly species to infer the most comprehensive molecular phylogeny of the family and curated a distribution dataset of more than 800,000 occurrences. We found strong evidence that species in environments with more stable daily temperatures or cooler maximum temperatures in the warm seasons have higher speciation rates. Furthermore, speciation and extinction rates decreased in tandem with global temperatures through geological time, resulting in a constant net diversification.

Large-scale genomic data reveal the phylogeny and evolution of owlet moths (Noctuoidea)

The owlet moths (Noctuoidea; ~43-45K described species) are one of the most ecologically diverse and speciose superfamilies of animals. Moreover, they comprise some of the world's most notorious pests of agriculture and forestry. Despite their contributions to terrestrial biodiversity and impacts on ecosystems and economies, the evolutionary history of Noctuoidea remains unclear because the superfamily lacks a statistically robust phylogenetic and temporal framework. We reconstructed the phylogeny of Noctuoidea using data from 1234 genes (946.4 kb nucleotides) obtained from the genome and transcriptome sequences of 76 species. The relationships among the six families of Noctuoidea were well resolved and consistently recovered based on both concatenation and gene coalescence approaches, supporting the following relationships: Oenosandridae + (Notodontidae + (Erebidae + (Nolidae + (Euteliidae + Noctuidae)))). A Yule tree prior with three unlinked molecular clocks was identified as the preferred BEAST analysis using marginal-likelihood estimations. The crown age of Noctuoidea was estimated at 74.5 Ma, with most families originating before the end of the Paleogene (23 Ma). Our study provides the first statistically robust phylogenetic and temporal framework for Noctuoidea, including all families of owlet moths, based on large-scale genomic data.

HerediGene Population Study IT infrastructure: A model to support genomic research recruitment and precision public health

The HerediGene Population Study is a large research study focused on identifying new genetic biomarkers for disease prevention, diagnosis, prognosis, and development of new therapeutics. A substantial IT infrastructure evolved to reach enrollment targets and return results to participants. More than 170,000 participants have been enrolled in the study to date, with 5.87% of those whole genome sequenced and 0.46% of those genotyped harboring pathogenic variants. Among other purposes, this infrastructure supports: (1) identifying candidates from clinical criteria, (2) monitoring for qualifying clinical events (e.g., blood draw), (3) contacting candidates, (4) obtaining consent electronically, (5) initiating lab orders, (6) integrating consent and lab orders into clinical workflow, (7) de-identifying samples and clinical data, (8) shipping/transmitting samples and clinical data, (9) genotyping/sequencing samples, (10) and re-identifying and returning results for participants where applicable. This study may serve as a model for similar genomic research and precision public health initiatives.

A global phylogeny of butterflies reveals their evolutionary history, ancestral hosts and biogeographic origins

Butterflies are a diverse and charismatic insect group that are thought to have evolved with plants and dispersed throughout the world in response to key geological events. However, these hypotheses have not been extensively tested because a comprehensive phylogenetic framework and datasets for butterfly larval hosts and global distributions are lacking. We sequenced 391 genes from nearly 2,300 butterfly species, sampled from 90 countries and 28 specimen collections, to reconstruct a new phylogenomic tree of butterflies representing 92% of all genera. Our phylogeny has strong support for nearly all nodes and demonstrates that at least 36 butterfly tribes require reclassification. Divergence time analyses imply an origin ~100 million years ago for butterflies and indicate that all but one family were present before the K/Pg extinction event. We aggregated larval host datasets and global distribution records and found that butterflies are likely to have first fed on Fabaceae and originated in what is now the Americas. Soon after the Cretaceous Thermal Maximum, butterflies crossed Beringia and diversified in the Palaeotropics. Our results also reveal that most butterfly species are specialists that feed on only one larval host plant family. However, generalist butterflies that consume two or more plant families usually feed on closely related plants.

The Development of an Infrastructure to Facilitate the Use of Whole Genome Sequencing for Population Health

The clinical use of genomic analysis has expanded rapidly resulting in an increased availability and utility of genomic information in clinical care. We have developed an infrastructure utilizing informatics tools and clinical processes to facilitate the use of whole genome sequencing data for population health management across the healthcare system. Our resulting framework scaled well to multiple clinical domains in both pediatric and adult care, although there were domain specific challenges that arose. Our infrastructure was complementary to existing clinical processes and well-received by care providers and patients. Informatics solutions were critical to the successful deployment and scaling of this program. Implementation of genomics at the scale of population health utilizes complicated technologies and processes that for many health systems are not supported by current information systems or in existing clinical workflows. To scale such a system requires a substantial clinical framework backed by informatics tools to facilitate the flow and management of data. Our work represents an early model that has been successful in scaling to 29 different genes with associated genetic conditions in four clinical domains. Work is ongoing to optimize informatics tools; and to identify best practices for translation to smaller healthcare systems.

The western redcedar genome reveals low genetic diversity in a self-compatible conifer

We assembled the 9.8 Gbp genome of western redcedar (WRC, Thuja plicata), an ecologically and economically important conifer species of the Cupressaceae. The genome assembly, derived from a uniquely inbred tree produced through five generations of self-fertilization (selfing), was determined to be 86% complete by BUSCO analysis - one of the most complete genome assemblies for a conifer. Population genomic analysis revealed WRC to be one of the most genetically depauperate wild plant species, with an effective population size of approximately 300 and no significant genetic differentiation across its geographic range. Nucleotide diversity, π, is low for a continuous tree species, with many loci exhibiting zero diversity, and the ratio of π at zero- to four-fold degenerate sites is relatively high (~ 0.33), suggestive of weak purifying selection. Using an array of genetic lines derived from up to five generations of selfing, we explored the relationship between genetic diversity and mating system. While overall heterozygosity was found to decline faster than expected during selfing, heterozygosity persisted at many loci, and nearly 100 loci were found to deviate from expectations of genetic drift, suggestive of associative overdominance. Nonreference alleles at such loci often harbor deleterious mutations and are rare in natural populations, implying that balanced polymorphisms are maintained by linkage to dominant beneficial alleles. This may account for how WRC remains responsive to natural and artificial selection, despite low genetic diversity.

Genomic selection reveals hidden relatedness and increased breeding efficiency in western redcedar polycross breeding

Western redcedar (WRC) is an ecologically and economically important forest tree species characterized by low genetic diversity with high self-compatibility and high heartwood durability. Using sequence capture genotyping of target genic and non-genic regions, we genotyped 44 parent trees and 1520 offspring trees representing 26 polycross (PX) families collected from three progeny test sites using 45,378 SNPs. Trees were phenotyped for eight traits related to growth, heartwood and foliar chemistry associated with wood durability and deer browse resistance. We used the genomic realized relationship matrix for paternity assignment, maternal pedigree correction, and to estimate genetic parameters. We compared genomics-based (GBLUP) and two pedigree-based (ABLUP: polycross and reconstructed full-sib [FS] pedigrees) models. Models were extended to estimate dominance genetic effects. Pedigree reconstruction revealed significant unequal male contribution and separated the 26 PX families into 438 FS families. Traditional maternal PX pedigree analysis resulted in up to 51% overestimation in genetic gain and 44% in diversity. Genomic analysis resulted in up to 22% improvement in offspring breeding value (BV) theoretical accuracy, 35% increase in expected genetic gain for forward selection, and doubled selection intensity for backward selection. Overall, all traits showed low to moderate heritability (0.09-0.28), moderate genotype by environment interaction (type-B genetic correlation: 0.51-0.80), low to high expected genetic gain (6.01%-55%), and no significant negative genetic correlation reflecting no large trade-offs for multi-trait selection. Only three traits showed a significant dominance effect. GBLUP resulted in smaller but more accurate heritability estimates for five traits, but larger estimates for the wood traits. Comparison between all, genic-coding, genic-non-coding and intergenic SNPs showed little difference in genetic estimates. In summary, we show that GBLUP overcomes the PX limitations, successfully captures expected historical and hidden relatedness as well as linkage disequilibrium (LD), and results in increased breeding efficiency in WRC.

Phylogeny and classification of Odonata using targeted genomics

Dragonflies and damselflies are a charismatic, medium-sized insect order (~6300 species) with a unique potential to approach comparative research questions. Their taxonomy and many ecological traits for a large fraction of extant species are relatively well understood. However, until now, the lack of a large-scale phylogeny based on high throughput data with the potential to connect both perspectives has precluded comparative evolutionary questions for these insects. Here, we provide an ordinal hypothesis of classification based on anchored hybrid enrichment using a total of 136 species representing 46 of the 48 families or incertae sedis, and a total of 478 target loci. Our analyses recovered the monophyly for all three suborders: Anisoptera, Anisozygoptera and Zygoptera. Although the backbone of the topology was reinforced and showed the highest support values to date, our genomic data was unable to stronglyresolve portions of the topology. In addition, a quartet sampling approach highlights the potential evolutionary scenarios that may have shaped evolutionary phylogeny (e.g., incomplete lineage sorting and introgression) of this taxon. Finally, in light of our phylogenomic reconstruction and previous morphological and molecular information we proposed an updated odonate classification and define five new families (Amanipodagrionidae fam. nov., Mesagrionidae fam. nov., Mesopodagrionidae fam. nov., Priscagrionidae fam. nov., Protolestidae fam. nov.) and reinstate another two (Rhipidolestidae stat. res., Tatocnemididae stat. res.). Additionally, we feature the problematic taxonomic groupings for examination in future studies to improve our current phylogenetic hypothesis.

A target enrichment probe set for resolving the flagellate land plant tree of life

PREMISE

New sequencing technologies facilitate the generation of large-scale molecular data sets for constructing the plant tree of life. We describe a new probe set for target enrichment sequencing to generate nuclear sequence data to build phylogenetic trees with any flagellate land plants, including hornworts, liverworts, mosses, lycophytes, ferns, and all gymnosperms.

METHODS

We leveraged existing transcriptome and genome sequence data to design the GoFlag 451 probes, a set of 56,989 probes for target enrichment sequencing of 451 exons that are found in 248 single-copy or low-copy nuclear genes across flagellate plant lineages.

RESULTS

Our results indicate that target enrichment using the GoFlag451 probe set can provide large nuclear data sets that can be used to resolve relationships among both distantly and closely related taxa across the flagellate land plants. We also describe the GoFlag 408 probes, an optimized probe set covering 408 of the 451 exons from the GoFlag 451 probe set that is commercialized by RAPiD Genomics.

CONCLUSIONS

A target enrichment approach using the new probe set provides a relatively low-cost solution to obtain large-scale nuclear sequence data for inferring phylogenetic relationships across flagellate land plants.

A target enrichment probe set for resolving the flagellate land plant tree of life

PREMISE

New sequencing technologies facilitate the generation of large-scale molecular data sets for constructing the plant tree of life. We describe a new probe set for target enrichment sequencing to generate nuclear sequence data to build phylogenetic trees with any flagellate land plants, including hornworts, liverworts, mosses, lycophytes, ferns, and all gymnosperms.

METHODS

We leveraged existing transcriptome and genome sequence data to design the GoFlag 451 probes, a set of 56,989 probes for target enrichment sequencing of 451 exons that are found in 248 single-copy or low-copy nuclear genes across flagellate plant lineages.

RESULTS

Our results indicate that target enrichment using the GoFlag451 probe set can provide large nuclear data sets that can be used to resolve relationships among both distantly and closely related taxa across the flagellate land plants. We also describe the GoFlag 408 probes, an optimized probe set covering 408 of the 451 exons from the GoFlag 451 probe set that is commercialized by RAPiD Genomics.

CONCLUSIONS

A target enrichment approach using the new probe set provides a relatively low-cost solution to obtain large-scale nuclear sequence data for inferring phylogenetic relationships across flagellate land plants.

A phylogenomic framework, evolutionary timeline and genomic resources for comparative studies of decapod crustaceans

Comprising over 15 000 living species, decapods (crabs, shrimp and lobsters) are the most instantly recognizable crustaceans, representing a considerable global food source. Although decapod systematics have received much study, limitations of morphological and Sanger sequence data have yet to produce a consensus for higher-level relationships. Here, we introduce a new anchored hybrid enrichment kit for decapod phylogenetics designed from genomic and transcriptomic sequences that we used to capture new high-throughput sequence data from 94 species, including 58 of 179 extant decapod families, and 11 of 12 major lineages. The enrichment kit yields 410 loci (greater than 86 000 bp) conserved across all lineages of Decapoda, more clade-specific molecular data than any prior study. Phylogenomic analyses recover a robust decapod tree of life strongly supporting the monophyly of all infraorders, and monophyly of each of the reptant, 'lobster' and 'crab' groups, with some results supporting pleocyemate monophyly. We show that crown decapods diverged in the Late Ordovician and most crown lineages diverged in the Triassic-Jurassic, highlighting a cryptic Palaeozoic history, and post-extinction diversification. New insights into decapod relationships provide a phylogenomic window into morphology and behaviour, and a basis to rapidly and cheaply expand sampling in this economically and ecologically significant invertebrate clade.

Phylogenomics resolves major relationships and reveals significant diversification rate shifts in the evolution of silk moths and relatives

BACKGROUND

Silkmoths and their relatives constitute the ecologically and taxonomically diverse superfamily Bombycoidea, which includes some of the most charismatic species of Lepidoptera. Despite displaying spectacular forms and diverse ecological traits, relatively little attention has been given to understanding their evolution and drivers of their diversity. To begin to address this problem, we created a new Bombycoidea-specific Anchored Hybrid Enrichment (AHE) probe set and sampled up to 571 loci for 117 taxa across all major lineages of the Bombycoidea, with a newly developed DNA extraction protocol that allows Lepidoptera specimens to be readily sequenced from pinned natural history collections.

RESULTS

The well-supported tree was overall consistent with prior morphological and molecular studies, although some taxa were misplaced. The bombycid Arotros Schaus was formally transferred to Apatelodidae. We identified important evolutionary patterns (e.g., morphology, biogeography, and differences in speciation and extinction), and our analysis of diversification rates highlights the stark increases that exist within the Sphingidae (hawkmoths) and Saturniidae (wild silkmoths).

CONCLUSIONS

Our study establishes a backbone for future evolutionary, comparative, and taxonomic studies of Bombycoidea. We postulate that the rate shifts identified are due to the well-documented bat-moth "arms race". Our research highlights the flexibility of AHE to generate genomic data from a wide range of museum specimens, both age and preservation method, and will allow researchers to tap into the wealth of biological data residing in natural history collections around the globe.

Origin and macroevolution of micro-moths on sunken Hawaiian Islands

The origins and evolution of Hawaiian biodiversity are a matter of controversy, and the mechanisms of lineage diversification for many organisms on this remote archipelago remain unclear. Here we focus on the poorly known endemic leaf-mining moth genus Philodoria (Lepidoptera, Gracillariidae), whose species feed on a diversity of Hawaiian plant lineages, many of which are critically endangered. We use anchored hybrid enrichment to assemble the first phylogenomic dataset (507 loci) for any Hawaiian animal taxon. To uncover the timing and pattern of diversification of these moths, we apply two frequently used dating calibration strategies, biogeographic calibrations and secondary calibrations. Island calibrations on their own resulted in much younger and unrealistic dates compared to strategies that relied on secondary calibrations. Philodoria probably originated on the now partially sunken islands of Laysan or Lisianski, approximately 21 Ma, and were associated with host plants in the families Ebenaceae, Malvaceae or Primulaceae. Major feeding groups associated with specific host-plant families originated soon after the plants colonized the islands. Allopatric isolation and host shifts, in concert and independently, probably play major roles in the diversification of Philodoria Our dating results indicate that Philodoria is among the oldest known Hawaiian arthropod lineages, and that island calibrations alone can lead to unrealistically young dates.

Anchored phylogenomics illuminates the skipper butterfly tree of life

BACKGROUND

Butterflies (Papilionoidea) are perhaps the most charismatic insect lineage, yet phylogenetic relationships among them remain incompletely studied and controversial. This is especially true for skippers (Hesperiidae), one of the most species-rich and poorly studied butterfly families.

METHODS

To infer a robust phylogenomic hypothesis for Hesperiidae, we sequenced nearly 400 loci using Anchored Hybrid Enrichment and sampled all tribes and more than 120 genera of skippers. Molecular datasets were analyzed using maximum-likelihood, parsimony and coalescent multi-species phylogenetic methods.

RESULTS

All analyses converged on a novel, robust phylogenetic hypothesis for skippers. Different optimality criteria and methodologies recovered almost identical phylogenetic trees with strong nodal support at nearly all nodes and all taxonomic levels. Our results support Coeliadinae as the sister group to the remaining skippers, the monotypic Euschemoninae as the sister group to all other subfamilies but Coeliadinae, and the monophyly of Eudaminae plus Pyrginae. Within Pyrginae, Celaenorrhinini and Tagiadini are sister groups, the Neotropical firetips, Pyrrhopygini, are sister to all other tribes but Celaenorrhinini and Tagiadini. Achlyodini is recovered as the sister group to Carcharodini, and Erynnini as sister group to Pyrgini. Within the grass skippers (Hesperiinae), there is strong support for the monophyly of Aeromachini plus remaining Hesperiinae. The giant skippers (Agathymus and Megathymus) once classified as a subfamily, are recovered as monophyletic with strong support, but are deeply nested within Hesperiinae.

CONCLUSIONS

Anchored Hybrid Enrichment sequencing resulted in a large amount of data that built the foundation for a new, robust evolutionary tree of skippers. The newly inferred phylogenetic tree resolves long-standing systematic issues and changes our understanding of the skipper tree of life. These resultsenhance understanding of the evolution of one of the most species-rich butterfly families.

Phylogenetics of moth-like butterflies (Papilionoidea: Hedylidae) based on a new 13-locus target capture probe set

The Neotropical moth-like butterflies (Hedylidae) are perhaps the most unusual butterfly family. In addition to being species-poor, this family is predominantly nocturnal and has anti-bat ultrasound hearing organs. Evolutionary relationships among the 36 described species are largely unexplored. A new, target capture, anchored hybrid enrichment probe set ('BUTTERFLY2.0') was developed to infer relationships of hedylids and some of their butterfly relatives. The probe set includes 13 genes that have historically been used in butterfly phylogenetics. Our dataset comprised of up to 10,898 aligned base pairs from 22 hedylid species and 19 outgroups. Eleven of the thirteen loci were successfully captured from all samples, and the remaining loci were captured from ≥94% of samples. The inferred phylogeny was consistent with recent molecular studies by placing Hedylidae sister to Hesperiidae, and the tree had robust support for 80% of nodes. Our results are also consistent with morphological studies, with Macrosoma tipulata as the sister species to all remaining hedylids, followed by M. semiermis sister to the remaining species in the genus. We tested the hypothesis that nocturnality evolved once from diurnality in Hedylidae, and demonstrate that the ancestral condition was likely diurnal, with a shift to nocturnality early in the diversification of this family. The BUTTERFLY2.0 probe set includes standard butterfly phylogenetics markers, captures sequences from decades-old museum specimens, and is a cost-effective technique to infer phylogenetic relationships of the butterfly tree of life.

Resolving Relationships among the Megadiverse Butterflies and Moths with a Novel Pipeline for Anchored Phylogenomics

The advent of next-generation sequencing technology has allowed for thecollection of large portions of the genome for phylogenetic analysis. Hybrid enrichment and transcriptomics are two techniques that leverage next-generation sequencing and have shown much promise. However, methods for processing hybrid enrichment data are still limited. We developed a pipeline for anchored hybrid enrichment (AHE) read assembly, orthology determination, contamination screening, and data processing for sequences flanking the target "probe" region. We apply this approach to study the phylogeny of butterflies and moths (Lepidoptera), a megadiverse group of more than 157,000 described species with poorly understood deep-level phylogenetic relationships. We introduce a new, 855 locus AHE kit for Lepidoptera phylogenetics and compare resulting trees to those from transcriptomes. The enrichment kit was designed from existing genomes, transcriptomes, and expressed sequence tags and was used to capture sequence data from 54 species from 23 lepidopteran families. Phylogenies estimated from AHE data were largely congruent with trees generated from transcriptomes, with strong support for relationships at all but the deepest taxonomic levels. We combine AHE and transcriptomic data to generate a new Lepidoptera phylogeny, representing 76 exemplar species in 42 families. The tree provides robust support for many relationships, including those among the seven butterfly families. The addition of AHE data to an existing transcriptomic dataset lowers node support along the Lepidoptera backbone, but firmly places taxa with AHE data on the phylogeny. Combining taxa sequenced for AHE with existing transcriptomes and genomes resulted in a tree with strong support for (Calliduloidea $+$ Gelechioidea $+$ Thyridoidea) $+$ (Papilionoidea $+$ Pyraloidea $+$ Macroheterocera). To examine the efficacy of AHE at a shallow taxonomic level, phylogenetic analyses were also conducted on a sister group representing a more recent divergence, the Saturniidae and Sphingidae. These analyses utilized sequences from the probe region and data flanking it, nearly doubled the size of the dataset; resulting trees supported new phylogenetics relationships, especially within the Saturniidae and Sphingidae (e.g., Hemarina derived in the latter). We hope that our data processing pipeline, hybrid enrichment gene set, and approach of combining AHE data with transcriptomes will be useful for the broader systematics community.

Resolving Relationships among the Megadiverse Butterflies and Moths with a Novel Pipeline for Anchored Phylogenomics

The advent of next-generation sequencing technology has allowed for thecollection of large portions of the genome for phylogenetic analysis. Hybrid enrichment and transcriptomics are two techniques that leverage next-generation sequencing and have shown much promise. However, methods for processing hybrid enrichment data are still limited. We developed a pipeline for anchored hybrid enrichment (AHE) read assembly, orthology determination, contamination screening, and data processing for sequences flanking the target "probe" region. We apply this approach to study the phylogeny of butterflies and moths (Lepidoptera), a megadiverse group of more than 157,000 described species with poorly understood deep-level phylogenetic relationships. We introduce a new, 855 locus AHE kit for Lepidoptera phylogenetics and compare resulting trees to those from transcriptomes. The enrichment kit was designed from existing genomes, transcriptomes, and expressed sequence tags and was used to capture sequence data from 54 species from 23 lepidopteran families. Phylogenies estimated from AHE data were largely congruent with trees generated from transcriptomes, with strong support for relationships at all but the deepest taxonomic levels. We combine AHE and transcriptomic data to generate a new Lepidoptera phylogeny, representing 76 exemplar species in 42 families. The tree provides robust support for many relationships, including those among the seven butterfly families. The addition of AHE data to an existing transcriptomic dataset lowers node support along the Lepidoptera backbone, but firmly places taxa with AHE data on the phylogeny. Combining taxa sequenced for AHE with existing transcriptomes and genomes resulted in a tree with strong support for (Calliduloidea $+$ Gelechioidea $+$ Thyridoidea) $+$ (Papilionoidea $+$ Pyraloidea $+$ Macroheterocera). To examine the efficacy of AHE at a shallow taxonomic level, phylogenetic analyses were also conducted on a sister group representing a more recent divergence, the Saturniidae and Sphingidae. These analyses utilized sequences from the probe region and data flanking it, nearly doubled the size of the dataset; resulting trees supported new phylogenetics relationships, especially within the Saturniidae and Sphingidae (e.g., Hemarina derived in the latter). We hope that our data processing pipeline, hybrid enrichment gene set, and approach of combining AHE data with transcriptomes will be useful for the broader systematics community.

A phylogenomic analysis of lichen-feeding tiger moths uncovers evolutionary origins of host chemical sequestration

Host species utilize a variety of defenses to deter feeding, including secondary chemicals. Some phytophagous insects have evolved tolerance to these chemical defenses, and can sequester secondary defense compounds for use against their own predators and parasitoids. While numerous studies have examined plant-insect interactions, little is known about lichen-insect interactions. Our study focused on reconstructing the evolution of lichen phenolic sequestration in the tiger moth tribe Lithosiini (Lepidoptera: Erebidae: Arctiinae), the most diverse lineage of lichen-feeding moths, with 3000 described species. We built an RNA-Seq dataset and examined the adult metabolome for the presence of lichen-derived phenolics. Using the transcriptomic dataset, we recover a well-resolved phylogeny of the Lithosiini, and determine that the metabolomes within species are more similar than those among species. Results from an initial ancestral state reconstruction suggest that the ability to sequester phenolics produced by a single chemical pathway preceded generalist sequestration of phenolics produced by multiple chemical pathways. We conclude that phenolics are consistently and selectively sequestered within Lithosiini. Furthermore, sequestration of compounds from a single chemical pathway may represent a synapomorphy of the tribe, and the ability to sequester phenolics produced by multiple pathways arose later. These findings expand on our understanding of the interactions between Lepidoptera and their lichen hosts.

Evidence for common horizontal transmission of Wolbachia among butterflies and moths

Background

Wolbachia is one of the most widespread bacteria on Earth. Previous research on Wolbachia-host interactions indicates that the bacterium is typically transferred vertically, from mother to offspring, through the egg cytoplasm. Although horizontal transmission of Wolbachia from one species to another is reported to be common in arthropods, limited direct ecological evidence is available. In this study, we examine horizontal transmission of Wolbachia using a multilocus sequence typing (MLST) strains dataset and used Wolbachia and Lepidoptera genomes to search for evidence for lateral gene transfer (LGT) in Lepidoptera, one of the most diverse cosmopolitan insect orders. We constructed a phylogeny of arthropod-associated MLST Wolbachia strains and calibrated the age of Wolbachia strains associated with lepidopteran species.

Results

Our results reveal inter-specific, inter-generic, inter-familial, and inter-ordinal horizontal transmission of Wolbachia strains, without discernible geographic patterns. We found at least seven probable cases of horizontal transmission among 31 species within Lepidoptera and between Lepidoptera and other arthropod hosts. The divergence time analysis revealed that Wolbachia is recently (22.6–4.7 mya, 95 % HPD) introduced in Lepidoptera. Analysis of nine Lepidoptera genomes (Bombyx mori, Danaus plexippus, Heliconius melpomene, Manduca sexta, Melitaea cinxia, Papilio glaucus, P. polytes, P. xuthus and Plutella xylostella) yielded one possible instance of Wolbachia LGT.

Conclusions

Our results provide evidence of high incidence of identical and multiple strains of Wolbachia among butterflies and moths, adding Lepidoptera to the growing body of evidence for common horizontal transmission of Wolbachia. This study demonstrates interesting dynamics of this remarkable and influential microorganism.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0660-x) contains supplementary material, which is available to authorized users.

Wallacellus is Euwallacea: molecular phylogenetics settles generic relationships (Coleoptera: Curculionidae: Scolytinae: Xyleborini)

Euwallacea Hopkins and Wallacellus Hulcr & Cognato are ambrosia beetle genera within the tribe Xyleborini (Coleoptera: Curculionidae: Scolytinae). Several species have recently received attention due to their establishment in non-native regions with serious ecological and economic consequences. To clarify generic placement of these species, we tested reciprocal monophyly of the two genera and the placement of several species using molecular phylogenetics. We sequenced, or re-used published sequences of, three markers (COI mtDNA, 28S nuclear rDNA and ArgK single-copy nuclear) from representatives of Euwallacea, Wallacellus, the Ambrosiodmus clade, and the clade containing Xyleborus s. str., and inferred their relationships with a Bayesian approach. We also tested explicit alternative topologies, and examined taxonomic utility of characters used for the delimitation of the genera.        All species of Euwallacea, Wallacellus, and two species of Xyleborus were monophyletic with high phylogenetic support. Based on the analysis and shared morphological characters, we transferred the following species to Euwallacea: Xyleborus declivispinatus (Schedl), Wallacellus piceus (Motschulsky), Xyleborus posticus (Eichhoff), Wallacellus similis (Ferrari), and Wallacellus striatulus (Browne). The genus Wallacellus was made a junior synonym of Euwallacea and morphological diagnosis of Euwallacea was updated. The results demonstrated that Euwallacea has a pantropical distribution.

Genetic Record for a Recent Invasion of Phenacoccus solenopsis (Hemiptera: Pseudococcidae) in Asia

The cotton mealybug Phenacoccus solenopsis Tinsley is an emerging invasive insect pest. Since its first report as a pest in the United States in 1991, it has invaded and colonized more than 23 countries over the past century. It was first recorded from Pakistan in 2006 and from China in 2008. In this study, we performed field surveys from 2010 to 2012 and obtained mtCOI sequences from specimens across China and Pakistan, then compared them with already available mtCOI sequences from additional Asian and North American countries. Our genetic analysis provides evidence that P. solenopsis should be classified into two groups, one of which is found only in the United States, and the other found only in Asia. The Asian group contains nine unique haplotypes, two of which have invaded and spread across China, Pakistan, India, and Vietnam over the last 4-6 yr. Our genetic analysis also indicates that P. solenopsis has a close relationship with the parasitoid wasp Aenasius bambawalei Hayat, providing preliminary evidence of a congruent spread of this mealybug and its parasitoids across China.

Moth tails divert bat attack: evolution of acoustic deflection

Adaptations to divert the attacks of visually guided predators have evolved repeatedly in animals. Using high-speed infrared videography, we show that luna moths (Actias luna) generate an acoustic diversion with spinning hindwing tails to deflect echolocating bat attacks away from their body and toward these nonessential appendages. We pit luna moths against big brown bats (Eptesicus fuscus) and demonstrate a survival advantage of ∼ 47% for moths with tails versus those that had their tails removed. The benefit of hindwing tails is equivalent to the advantage conferred to moths by bat-detecting ears. Moth tails lured bat attacks to these wing regions during 55% of interactions between bats and intact luna moths. We analyzed flight kinematics of moths with and without hindwing tails and suggest that tails have a minimal role in flight performance. Using a robust phylogeny, we find that long spatulate tails have independently evolved four times in saturniid moths, further supporting the selective advantage of this anti-bat strategy. Diversionary tactics are perhaps more common than appreciated in predator-prey interactions. Our finding suggests that focusing on the sensory ecologies of key predators will reveal such countermeasures in prey.

Underground evolution: new roots for the old tree of lumbricid earthworms

Earthworms belonging to the family Lumbricidae are extremely abundant in terrestrial temperate regions. They affect soil properties and nutrient cycling, thus shaping plant community composition and aboveground food webs. Some lumbricids are also model organisms in ecology and toxicology. Despite the intense research efforts dedicated to lumbricids over the last 130years, the evolutionary relationships and taxonomic classification of these organisms are still subject to great debate. Resolution of their systematics is hampered by the structural simplicity of the earthworm body plan and the existence of cryptic species. We sampled 160 earthworm specimens belonging to 84 lumbricid species (28 genera) and 22 Lumbricoidea outgroups, sequenced two nuclear genes, four mitochondrial genes and seven mitochondrial tRNAs and examined 22 morphological characters. We then applied a combination of phylogenetic methods to generate the first robust genus-level phylogeny of the Lumbricidae. Our results show that the current Lumbricidae classification and the underlying hypotheses of character evolution must be revised. Our chronogram suggests that lumbricids emerged in the Lower Cretaceous in the holarctic region and that their diversification has been driven by tectonic processes (e.g. Laurasia split) and geographical isolation. Our chronogram and character reconstruction analysis reveal that spermathecae number does not follow a gradual pattern of reduction and that parthenogenesis arose from sexual relatives multiple times in the group; the same analysis also indicates that both epigeic and anecic earthworms evolved from endogeic ancestors. These findings emphasize the strong and multiple changes to which morphological and ecological characters are subjected, challenging the hypothesis of character stasis in Lumbricidae.

Using phylogenetically-informed annotation (PIA) to search for light-interacting genes in transcriptomes from non-model organisms

Background

Tools for high throughput sequencing and de novo assembly make the analysis of transcriptomes (i.e. the suite of genes expressed in a tissue) feasible for almost any organism. Yet a challenge for biologists is that it can be difficult to assign identities to gene sequences, especially from non-model organisms. Phylogenetic analyses are one useful method for assigning identities to these sequences, but such methods tend to be time-consuming because of the need to re-calculate trees for every gene of interest and each time a new data set is analyzed. In response, we employed existing tools for phylogenetic analysis to produce a computationally efficient, tree-based approach for annotating transcriptomes or new genomes that we term Phylogenetically-Informed Annotation (PIA), which places uncharacterized genes into pre-calculated phylogenies of gene families.

Results

We generated maximum likelihood trees for 109 genes from a Light Interaction Toolkit (LIT), a collection of genes that underlie the function or development of light-interacting structures in metazoans. To do so, we searched protein sequences predicted from 29 fully-sequenced genomes and built trees using tools for phylogenetic analysis in the Osiris package of Galaxy (an open-source workflow management system). Next, to rapidly annotate transcriptomes from organisms that lack sequenced genomes, we repurposed a maximum likelihood-based Evolutionary Placement Algorithm (implemented in RAxML) to place sequences of potential LIT genes on to our pre-calculated gene trees. Finally, we implemented PIA in Galaxy and used it to search for LIT genes in 28 newly-sequenced transcriptomes from the light-interacting tissues of a range of cephalopod mollusks, arthropods, and cubozoan cnidarians. Our new trees for LIT genes are available on the Bitbucket public repository (http://bitbucket.org/osiris_phylogenetics/pia/) and we demonstrate PIA on a publicly-accessible web server (http://galaxy-dev.cnsi.ucsb.edu/pia/).

Conclusions

Our new trees for LIT genes will be a valuable resource for researchers studying the evolution of eyes or other light-interacting structures. We also introduce PIA, a high throughput method for using phylogenetic relationships to identify LIT genes in transcriptomes from non-model organisms. With simple modifications, our methods may be used to search for different sets of genes or to annotate data sets from taxa outside of Metazoa.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0350-x) contains supplementary material, which is available to authorized users.

Phylotranscriptomics: saturated third codon positions radically influence the estimation of trees based on next-gen data

Recent advancements in molecular sequencing techniques have led to a surge in the number of phylogenetic studies that incorporate large amounts of genetic data. We test the assumption that analyzing large number of genes will lead to improvements in tree resolution and branch support using moths in the superfamily Bombycoidea, a group with some interfamilial relationships that have been difficult to resolve. Specifically, we use a next-gen data set that included 19 taxa and 938 genes (∼1.2M bp) to examine how codon position and saturation might influence resolution and node support among three key families. Maximum likelihood, parsimony, and species tree analysis using gene tree parsimony, on different nucleotide and amino acid data sets, resulted in largely congruent topologies with high bootstrap support compared with prior studies that included fewer loci. However, for a few shallow nodes, nucleotide and amino acid data provided high support for conflicting relationships. The third codon position was saturated and phylogenetic analysis of this position alone supported a completely different, potentially misleading sister group relationship. We used the program RADICAL to assess the number of genes needed to fix some of these difficult nodes. One such node originally needed a total of 850 genes but only required 250 when synonymous signal was removed. Our study shows that, in order to effectively use next-gen data to correctly resolve difficult phylogenetic relationships, it is necessary to assess the effects of synonymous substitutions and third codon positions.

Phylogenomics provides strong evidence for relationships of butterflies and moths

Butterflies and moths constitute some of the most popular and charismatic insects. Lepidoptera include approximately 160 000 described species, many of which are important model organisms. Previous studies on the evolution of Lepidoptera did not confidently place butterflies, and many relationships among superfamilies in the megadiverse clade Ditrysia remain largely uncertain. We generated a molecular dataset with 46 taxa, combining 33 new transcriptomes with 13 available genomes, transcriptomes and expressed sequence tags (ESTs). Using HaMStR with a Lepidoptera-specific core-orthologue set of single copy loci, we identified 2696 genes for inclusion into the phylogenomic analysis. Nucleotides and amino acids of the all-gene, all-taxon dataset yielded nearly identical, well-supported trees. Monophyly of butterflies (Papilionoidea) was strongly supported, and the group included skippers (Hesperiidae) and the enigmatic butterfly-moths (Hedylidae). Butterflies were placed sister to the remaining obtectomeran Lepidoptera, and the latter was grouped with greater than or equal to 87% bootstrap support. Establishing confident relationships among the four most diverse macroheteroceran superfamilies was previously challenging, but we recovered 100% bootstrap support for the following relationships: ((Geometroidea, Noctuoidea), (Bombycoidea, Lasiocampoidea)). We present the first robust, transcriptome-based tree of Lepidoptera that strongly contradicts historical placement of butterflies, and provide an evolutionary framework for genomic, developmental and ecological studies on this diverse insect order.

The emergence of lobsters: phylogenetic relationships, morphological evolution and divergence time comparisons of an ancient group (decapoda: achelata, astacidea, glypheidea, polychelida)

Lobsters are a ubiquitous and economically important group of decapod crustaceans that include the infraorders Polychelida, Glypheidea, Astacidea and Achelata. They include familiar forms such as the spiny, slipper, clawed lobsters and crayfish and unfamiliar forms such as the deep-sea and "living fossil" species. The high degree of morphological diversity among these infraorders has led to a dynamic classification and conflicting hypotheses of evolutionary relationships. In this study, we estimated phylogenetic relationships among the major groups of all lobster families and 94% of the genera using six genes (mitochondrial and nuclear) and 195 morphological characters across 173 species of lobsters for the most comprehensive sampling to date. Lobsters were recovered as a non-monophyletic assemblage in the combined (molecular + morphology) analysis. All families were monophyletic, with the exception of Cambaridae, and 7 of 79 genera were recovered as poly- or paraphyletic. A rich fossil history coupled with dense taxon coverage allowed us to estimate and compare divergence times and origins of major lineages using two drastically different approaches. Age priors were constructed and/or included based on fossil age information or fossil discovery, age, and extant species count data. Results from the two approaches were largely congruent across deep to shallow taxonomic divergences across major lineages. The origin of the first lobster-like decapod (Polychelida) was estimated in the Devonian (∼409-372 Ma) with all infraorders present in the Carboniferous (∼353-318 Ma). Fossil calibration subsampling studies examined the influence of sampling density (number of fossils) and placement (deep, middle, and shallow) on divergence time estimates. Results from our study suggest including at least 1 fossil per 10 operational taxonomic units (OTUs) in divergence dating analyses. [Dating; decapods; divergence; lobsters; molecular; morphology; phylogenetics.].

Evolution of Manduca sexta hornworms and relatives: biogeographical analysis reveals an ancestral diversification in Central America

The hawkmoth genus Manduca is a diverse group of very large, conspicuous moths that has served as an important model across many biological disciplines. Two species in particular, the tobacco hornworm (Manduca sexta) and the tomato hornworm (Manduca quinquemaculatus) have been researched extensively. Studies across biological fields have referred to these two species as being closely related or even sister species, but the extent to which these two model organisms are related remains largely unknown. We conducted a comprehensive multi-gene phylogenetic analysis of Manduca, based on both an ML and Bayesian framework, which resulted in a monophyletic Manduca but only when two other genera, Dolba and Euryglottis are included. We tentatively conclude that the sister group to Manduca sexta comprises the Caribbean M. afflicta and M. johanni, and the sister lineage to this clade includes M. quinquemaculatus and the Hawaiian M. blackburni. Thus, M. sexta and M. quinquemaculatus are closely related, but are not sister species. Biogeographical analyses reveal an ancestral center of diversification in Central America, and Manduca appears to have subsequently colonized North and South America. Our phylogeny provides an important foundation for comparative studies of two model organisms and their relatives.

Phylogeny and evolutionary patterns in the Dwarf crayfish subfamily (Decapoda: Cambarellinae)

The Dwarf crayfish or Cambarellinae, is a morphologically singular subfamily of decapod crustaceans that contains only one genus, Cambarellus. Its intriguing distribution, along the river basins of the Gulf Coast of United States (Gulf Group) and into Central México (Mexican Group), has until now lacked of satisfactory explanation. This study provides a comprehensive sampling of most of the extant species of Cambarellus and sheds light on its evolutionary history, systematics and biogeography. We tested the impact of Gulf Group versus Mexican Group geography on rates of cladogenesis using a maximum likelihood framework, testing different models of birth/extinction of lineages. We propose a comprehensive phylogenetic hypothesis for the subfamily based on mitochondrial and nuclear loci (3,833 bp) using Bayesian and Maximum Likelihood methods. The phylogenetic structure found two phylogenetic groups associated to the two main geographic components (Gulf Group and Mexican Group) and is partially consistent with the historical structure of river basins. The previous hypothesis, which divided the genus into three subgenera based on genitalia morphology was only partially supported (P = 0.047), resulting in a paraphyletic subgenus Pandicambarus. We found at least two cases in which phylogenetic structure failed to recover monophyly of recognized species while detecting several cases of cryptic diversity, corresponding to lineages not assigned to any described species. Cladogenetic patterns in the entire subfamily are better explained by an allopatric model of speciation. Diversification analyses showed similar cladogenesis patterns between both groups and did not significantly differ from the constant rate models. While cladogenesis in the Gulf Group is coincident in time with changes in the sea levels, in the Mexican Group, cladogenesis is congruent with the formation of the Trans-Mexican Volcanic Belt. Our results show how similar allopatric divergence in freshwater organisms can be promoted through diverse vicariant factors.

Testing phylogenetic hypotheses of the subgenera of the freshwater crayfish genus Cambarus (Decapoda: Cambaridae)

Background

The genus Cambarus is one of three most species rich crayfish genera in the Northern Hemisphere. The genus has its center of diversity in the Southern Appalachians of the United States and has been divided into 12 subgenera. Using Cambarus we test the correspondence of subgeneric designations based on morphology used in traditional crayfish taxonomy to the underlying evolutionary history for these crayfish. We further test for significant correlation and explanatory power of geographic distance, taxonomic model, and a habitat model to estimated phylogenetic distance with multiple variable regression.

Methodology/Principal Findings

We use three mitochondrial and one nuclear gene regions to estimate the phylogenetic relationships for species within the genus Cambarus and test evolutionary hypotheses of relationships and associated morphological and biogeographical hypotheses. Our resulting phylogeny indicates that the genus Cambarus is polyphyletic, however we fail to reject the monophyly of Cambarus with a topology test. The majority of the Cambarus subgenera are rejected as monophyletic, suggesting the morphological characters used to define those taxa are subject to convergent evolution. While we found incongruence between taxonomy and estimated phylogenetic relationships, a multiple model regression analysis indicates that taxonomy had more explanatory power of genetic relationships than either habitat or geographic distance.

Conclusions

We find convergent evolution has impacted the morphological features used to delimit Cambarus subgenera. Studies of the crayfish genus Orconectes have shown gonopod morphology used to delimit subgenera is also affected by convergent evolution. This suggests that morphological diagnoses based on traditional crayfish taxonomy might be confounded by convergent evolution across the cambarids and has little utility in diagnosing relationships or defining natural groups. We further suggest that convergent morphological evolution appears to be a common occurrence in invertebrates suggesting the need for careful phylogenetically based interpretations of morphological evolution in invertebrate systematics.

An earthworm riddle: systematics and phylogeography of the Spanish lumbricid Postandrilus

Background

As currently defined, the genus Postandrilus Qui and Bouché, 1998, (Lumbricidae) includes six earthworm species, five occurring in Majorca (Baleares Islands, western Mediterranean) and another in Galicia (NW Spain). This disjunct and restricted distribution raises some interesting phylogeographic questions: (1) Is Postandrilus distribution the result of the separation of the Baleares-Kabylies (BK) microplate from the proto-Iberian Peninsula in the Late Oligocene (30–28 Mya) – vicariant hypothesis? (2) Did Postandrilus diversify in Spain and then colonize the Baleares during the Messinian salinity crisis (MSC) 5.96–5.33 Mya – dispersal hypothesis? (3) Is the distribution the result of a two-step process – vicariance with subsequent dispersal?

Methodology/Principal Findings

To answer these questions and assess Postandrilus evolutionary relationships and systematics, we collected all of the six Postandrilus species (46 specimens – 16 locations) and used Aporrectodea morenoe and three Prosellodrilus and two Cataladrilus species as the outgroup. Regions of the nuclear 28S rDNA and mitochondrial 16S rDNA, 12S rDNA, ND1, COII and tRNA genes (4,666 bp) were sequenced and analyzed using maximum likelihood and Bayesian methods of phylogenetic and divergence time estimation. The resulting trees revealed six new Postandrilus species in Majorca that clustered with the other five species already described. This Majorcan clade was sister to an Iberian clade including A. morenoe (outgroup) and Postandrilus bertae. Our phylogeny and divergence time estimates indicated that the split between the Iberian and Majorcan Postandrilus clades took place 30.1 Mya, in concordance with the break of the BK microplate from the proto-Iberian Peninsula, and that the present Majorcan clade diversified 5.7 Mya, during the MSC.

Conclusions

Postandrilus is highly diverse including multiple cryptic species in Majorca. The genus is not monophyletic and invalid as currently defined. Postandrilus is of vicariant origin and its radiation began in the Late Oligocene.

Population genetic structure of an endangered Utah endemic, Astragalus ampullarioides (Fabaceae)

The endangered Shivwits milkvetch, Astragalus ampullarioides, is a perennial, herbaceous plant. This Utah endemic was federally listed as endangered in 2001 because of its high habitat specificity and low numbers of individuals and populations. All habitat currently occupied by A. ampullarioides was designated as critical by the U.S. Fish and Wildlife Service in 2006 as a result of conservation litigation. We used AFLP markers to assess genetic differentiation among the seven extant populations and quantified genetic diversity in each. Six different AFLP markers resulted in 217 unambiguous polymorphic loci. We used multiple methods to examine any changes in population genetic structure in this species over time. Results indicate that A. ampullarioides had much higher gene flow among populations in the past, but has since fragmented into regional genetic units. These regions further fragmented genetically, and extant populations have differentiated through genetic drift. Populations had low levels of gene flow, even between geographically close populations. Rapid urban development reduces gene flow among regions and encroaches on populations of A. ampullarioides and remaining patches of unoccupied habitat. The genetic makeup of each of the extant populations should be carefully considered in management decisions such as population establishment or augmentation.