Whole-genome amplification in double-digest RADseq results in adequate libraries but fewer sequenced loci

Population genomics and connectivity of Vazella pourtalesii sponge grounds of the northwest Atlantic with conservation implications of deep sea vulnerable marine ecosystems

Sponges are key ecosystem engineers that shape, structure and enhance the biodiversity of marine benthic communities globally. Sponge aggregations and reefs are recognized as vulnerable marine ecosystems (or VMEs) due to their susceptibility to damage from bottom-contact fishing gears. Ensuring their long-term sustainability, preservation, and ecosystem functions requires the implementation of sound scientific conservation tools. Here, the genetic diversity, structure, and connectivity of the deep-sea glass sponge, Vazella pourtalesii (Schmidt, 1870), was investigated using 1,102 neutral SNPs obtained in RADseq. This species is distributed across the northwest Atlantic from Florida, USA to Nova Scotia, Canada and we sequenced samples covering this full distribution and provided evidence of strong genetic structure with two distinct clusters: Florida together with the Carolina Shelves and the Scotian Shelf. We estimated moderate levels of diversity with low migration across large distances (> 1000 kms) and high connectivity at smaller scales (< 300 kms). Further, fishing pressure on genetic diversity was evaluated, within two Sponge Conservation Areas (SCAs) on the Scotian Shelf. Those areas have different disturbance histories, and cumulative fishing pressure. Slightly lower levels of genetic diversity were found inside the SCAs, and yet they encompassed a high proportion of the diversity observed within the Scotian Shelf. We provide baseline data for future monitoring of the SCAs, discussing our findings in the light of existing area-based management tools.

Generation of genome-wide SNP markers from minimally invasive sampling in endangered animals and applications in species ecology and conservation

High-density genotyping methods have revolutionized the field of population and conservation genetics in the past decade. To exploit the technological and analytical advances in the field, access to high-quality genetic material is a key component. However, access to such samples in endangered and rare animals is often challenging or even impossible. Here, we used a minimally invasive sampling method (MIS) in the endangered cave salamander Proteus anguinus, the olm, to generate thousands of genetic markers using ddRADseq for population and conservation genomic analyses. Using tail clips and MIS skin swabs taken from the same individual, we investigated genotyping data properties of the two different sampling types. We found that sufficient DNA can be extracted from swab samples to generate up to 200,000 polymorphic SNPs in divergent Proteus lineages. Swab and tissue samples were highly reproducible exhibiting low SNP genotyping error rates. We found that SNPs were most frequently (~50%) located within genic regions, while the rest mapped to mostly flanking regions of repetitive DNA. The vast majority of DNA recovered from swabbing was host DNA. However, a fraction of DNA recovered from swabs contained additional ecological information on the species, including eDNA from the surrounding environment and bacterial skin fauna. Most exogenous DNA recovered from swabs were bacteria (~80%), followed by vertebrates (~20%). Our results demonstrate that MIS can be used to (i) generate tens of thousands of ddRADseq markers for conservation and population genomic analyses and (ii) inform on the species health status and ecology from exogenous DNA.

Using genome-wide data to ascertain taxonomic status and assess population genetic structure for Houston toads (Bufo [= Anaxyrus] houstonensis)

The Houston toad (Bufo [= Anaxyrus] houstonensis) is an endangered amphibian with a small geographic range. Land-use changes have primarily driven decline in B. houstonensis with population supplementation predominant among efforts to reduce its current extinction risk. However, there has been historic uncertainty regarding the evolutionary and conservation significance of B. houstonensis. To this end, we used 1170 genome-wide nuclear DNA markers to examine phylogenetic relationships between our focal taxon, representatives of the Nearctic B. americanus group, and B. nebulifer, a sympatric Middle American species. Phylogenetic analyses indicate B. houstonensis is a taxon that is distinct from B. americanus. We corroborated such genetic distinctiveness with an admixture analysis that provided support for recent reproductive isolation between B. americanus and B. houstonensis. However, ABBA-BABA tests for ancient admixture indicated historic gene flow between Nearctic species while no signal of historic gene flow was detected between Nearctic and Middle-American species. We used an admixture analysis to recognize four Management Units (MU) based on observed genetic differentiation within B. houstonensis and recommend captive propagation, population supplementation, and habitat restoration efforts specific to each MU. Our results re-affirm the evolutionary novelty of an endangered relict.

Continent-island boundary and environment-shaped evolution in the marine amphipod Ampithoe marcuzzii complex (Crustacea: Eumalacostraca: Ampithoidae)

Marine amphipods are crustaceans that lack a larval phase and consequently have low dispersion rates. Despite that, these crustaceans present a remarkable ability to be transported by rafting on natural floating substrata, especially macroalgae, where they find shelter, food and a mating ground. The species Ampithoe marcuzzii is widely distributed throughout the western Atlantic Ocean. Here, it was used as a model to study seascape genomics and phylogeography in invertebrates with low dispersion capacities. We anticipated that the lineages would present isolation-by-distance patterns. However, surface currents and other abiotic variables could facilitate connectivity among distant sites. Based on mitochondrial and nuclear genes, SNPs, and environmental associations, we observed the presence of a species complex within A. marcuzzii, separating mainland and insular populations. Each species showed an independent evolutionary history, with a strong latitudinal population structure and evidence of isolation-by-distance and isolation-by-environment, characterizing the 'continent' species. Historical expansion and environmental variables were observed associated with the southeastern population, and ecological niche modeling corroborated the region as a paleorefuge. Conversely, populations from 'islands' presented complicated evolutionary histories, with closer localities genetically isolated and distant localities connected. These findings indicate that insular populations with low dispersion capacity might be more susceptible to spatial connectivity by floating substrata and to changes in surface currents. In contrast, mainland populations might be more vulnerable to local climate changes due to lack of gene flow.

Genomic evidence for homoploid hybrid speciation in a marine mammal apex predator

Hybridization is widespread and constitutes an important source of genetic variability and evolution. In animals, its role in generating novel and independent lineages (hybrid speciation) has been strongly debated, with only a few cases supported by genomic data. The South American fur seal (SAfs) Arctocephalus australis is a marine apex predator of Pacific and Atlantic waters, with a disjunct set of populations in Peru and Northern Chile [Peruvian fur seal (Pfs)] with controversial taxonomic status. We demonstrate, using complete genome and reduced representation sequencing, that the Pfs is a genetically distinct species with an admixed genome that originated from hybridization between the SAfs and the Galapagos fur seal (Arctocephalus galapagoensis) ~400,000 years ago. Our results strongly support the origin of Pfs by homoploid hybrid speciation over alternative introgression scenarios. This study highlights the role of hybridization in promoting species-level biodiversity in large vertebrates.

Testing assertions of widespread introgressive hybridization in a clade of neotropical toads with low mate selectivity (Rhinella granulosa species group)

Discordance between different genomic regions, often identified through multilocus sequencing of selected markers, presents particular difficulties in identifying historical processes which drive species diversity and boundaries. Mechanisms causing discordance, such as incomplete lineage sorting or introgression due to interspecific hybridization, are better identified based on population-level genomic datasets. In the toads of the Rhinella granulosa species group, patterns of mito-nuclear discordance and potential hybridization have been reported by several studies. However, these patterns were proposed based on few loci, such that alternative mechanisms behind gene-tree heterogeneity cannot be ruled out. Using genome-wide ddRADseq loci from a subset of species within this clade, we found only partial concordance between currently recognized species-level taxon boundaries and patterns of genetic structure. While most taxa within the R. granulosa group correspond to clades, genetic clustering analyses sometimes grouped distinct taxonomic units into a single cluster. Moreover, levels of admixture between inferred clusters were limited and restricted to a single taxon pair which is best explained by incomplete lineage sorting as opposed to introgressive hybridization, according to D-statistics results. These findings contradict previous assertions of widespread cryptic diversity and gene flow within the R. granulosa clade. Lastly, our analyses suggest that diversification events within the Rhinella granulosa group mostly dated back to the early Pliocene, being generally younger than species divergences in other closely related clades that present high levels of cross-species gene flow. This finding uniquely contradicts common assertions that this young clade of toads exhibits interspecific hybridization.

Cost-saving population genomic investigation of Daphnia longispina complex resting eggs using whole-genome amplification and pre-sequencing screening

Resting stages of aquatic organisms that accumulate in the sediment over time are an exceptional resource that allows direct insights into past populations and addressing evolutionary questions. This is of particular interest in taxa that face relatively new environmental challenges, e.g., climate change and eutrophication, such as the Daphnia longispina species complex, a keystone zooplankton group in European freshwater ecosystems. However, genomic analysis might be challenging as DNA yield from many of these resting stages can be low and the material degraded. To reliably allow the resequencing of single Daphnia resting eggs from different sediment layers and characterize genomic changes through time, we performed whole-genome amplification to obtain DNA amounts suitable for genome resequencing and tested multiple protocols involving egg isolation, whole-genome amplification kits, and library preparation. A pre-sequencing contamination screening was developed, consisting of amplifying mitochondrial Daphnia and bacterial markers, to quickly assess and exclude possibly contaminated samples. In total, we successfully amplified and sequenced nine genomes from Daphnia resting eggs that could be identified as Daphnia longispina species. We analyzed the genome coverage and heterozygosity of these samples to optimize this method for future projects involving population genomic investigation of the resting egg bank.

Genome-wide footprints in the carob tree (Ceratonia siliqua) unveil a new domestication pattern of a fruit tree in the Mediterranean

Intense research efforts over the last two decades have renewed our understanding of plant phylogeography and domestication in the Mediterranean basin. Here we aim to investigate the evolutionary history and the origin of domestication of the carob tree (Ceratonia siliqua), which has been cultivated for millennia for food and fodder. We used >1000 microsatellite genotypes to delimit seven carob evolutionary units (CEUs). We investigated genome‐wide diversity and evolutionary patterns of the CEUs with 3557 single nucleotide polymorphisms generated by restriction‐site associated DNA sequencing (RADseq). To address the complex wild vs. cultivated status of sampled trees, we classified 56 sampled populations across the Mediterranean basin as wild, seminatural or cultivated. Nuclear and cytoplasmic loci were identified from RADseq data and separated for analyses. Phylogenetic analyses of these genomic‐wide data allowed us to resolve west‐to‐east expansions from a single long‐term refugium probably located in the foothills of the High Atlas Mountains near the Atlantic coast. Our findings support multiple origins of domestication with a low impact on the genetic diversity at range‐wide level. The carob was mostly domesticated from locally selected wild genotypes and scattered long‐distance westward dispersals of domesticated varieties by humans, concomitant with major historical migrations by Romans, Greeks and Arabs. Ex situ efforts to preserve carob genetic resources should prioritize accessions from both western and eastern populations, with emphasis on the most differentiated CEUs situated in southwest Morocco, south Spain and eastern Mediterranean. Our study highlights the relevance of wild and seminatural habitats in the conservation of genetic resources for cultivated trees.

Fine-scale genome-wide signature of Pleistocene glaciation in Thitarodes moths (Lepidoptera: Hepialidae), host of Ophiocordyceps fungus in the Hengduan Mountains

The Hengduan Mountains region is a biodiversity hotspot known for its topologically complex, deep valleys and high mountains. While landscape and glacial refugia have been evoked to explain patterns of inter-species divergence, the accumulation of intra-species (i.e. population level) genetic divergence across the mountain-valley landscape in this region has received less attention. We used genome-wide restriction site-associated DNA sequencing (RADseq) to reveal signatures of Pleistocene glaciation in populations of Thitarodes shambalaensis (Lepidoptera: Hepialidae), the host moth of parasitic Ophiocordyceps sinensis (Hypocreales: Ophiocordycipitaceae) or "caterpillar fungus" endemic to the glacier of eastern Mt. Gongga. We used moraine history along the glacier valleys to model the distribution and environmental barriers to gene flow across populations of T. shambalaensis. We found that moth populations separated by less than 10 km exhibited valley-based population genetic clustering and isolation-by-distance (IBD), while gene flow among populations was best explained by models using information about their distributions at the local last glacial maximum (LGM_L , 58 kya), not their contemporary distribution. Maximum likelihood lineage history among populations, and among subpopulations as little as 500 meters apart, recapitulated glaciation history across the landscape. We also found signals of isolated population expansion following the retreat of LGM_L glaciers. These results reveal the fine-scale, long-term historical influence of landscape and glaciation on the genetic structuring of populations of an endangered and economically important insect species. Similar mechanisms, given enough time and continued isolation, could explain the contribution of glacier refugia to the generation of species diversity among the Hengduan Mountains.

Phylogenomics, introgression, and demographic history of South American true toads (Rhinella)

The effects of genetic introgression on species boundaries and how they affect species' integrity and persistence over evolutionary time have received increased attention. The increasing availability of genomic data has revealed contrasting patterns of gene flow across genomic regions, which impose challenges to inferences of evolutionary relationships and of patterns of genetic admixture across lineages. By characterizing patterns of variation across thousands of genomic loci in a widespread complex of true toads (Rhinella), we assess the true extent of genetic introgression across species thought to hybridize to extreme degrees based on natural history observations and multi-locus analyses. Comprehensive geographic sampling of five large-ranged Neotropical taxa revealed multiple distinct evolutionary lineages that span large geographic areas and, at times, distinct biomes. The inferred major clades and genetic clusters largely correspond to currently recognized taxa; however, we also found evidence of cryptic diversity within taxa. While previous phylogenetic studies revealed extensive mito-nuclear discordance, our genetic clustering analyses uncovered several admixed individuals within major genetic groups. Accordingly, historical demographic analyses supported that the evolutionary history of these toads involved cross-taxon gene flow both at ancient and recent times. Lastly, ABBA-BABA tests revealed widespread allele sharing across species boundaries, a pattern that can be confidently attributed to genetic introgression as opposed to incomplete lineage sorting. These results confirm previous assertions that the evolutionary history of Rhinella was characterized by various levels of hybridization even across environmentally heterogeneous regions, posing exciting questions about what factors prevent complete fusion of diverging yet highly interdependent evolutionary lineages.

A role of asynchrony of seasons in explaining genetic differentiation in a Neotropical toad

The process of diversification can be studied at the phylogeographic level by attempting to identify the environmental features that promote and maintain population divergence. Here we investigate diversification in Rhinella granulosa, a Neotropical toad from northeastern Brazil, by testing a range of hypotheses that encompass different putative mechanisms reducing gene flow among populations. We sequenced single nucleotide polymorphisms and examined individual predictions related to the role of geographic barriers (rivers), ecological gradients, historical habitat stability, and spatial variation in climate seasonality, also known as the asynchrony of seasons hypothesis. This hypothesis postulates that temporal asynchrony of wet and dry seasons over short distances causes parapatric populations to become isolated by time. After determining genetic structure, inferring past distributions, ranking demographic models, and estimating the power of monthly climatic variables, our results identified two populations that are not associated with geographic barriers, biome gradients, or historical refugia. Instead, they are predicted by spatial variation in monthly rainfall and minimum temperature, consistent with the asynchrony of seasons hypothesis, supported also by our comparative framework using multiple matrix regression and linear mixed effects modeling. Due to the toad's life history, climate likely mediates gene flow directly, with genetic differentiation being provoked by neutral mechanisms related to climate driven population isolation, and/or by natural selection against migrants from populations with different breeding times. The asynchrony of seasons hypothesis is seldom considered in phylogeographic studies, but our results indicate that it should be tested in systems where breeding is tightly coupled with climate.

Facilitating population genomics of non-model organisms through optimized experimental design for reduced representation sequencing

BACKGROUND

Genome-wide data are invaluable to characterize differentiation and adaptation of natural populations. Reduced representation sequencing (RRS) subsamples a genome repeatedly across many individuals. However, RRS requires careful optimization and fine-tuning to deliver high marker density while being cost-efficient. The number of genomic fragments created through restriction enzyme digestion and the sequencing library setup must match to achieve sufficient sequencing coverage per locus. Here, we present a workflow based on published information and computational and experimental procedures to investigate and streamline the applicability of RRS.

RESULTS

In an iterative process genome size estimates, restriction enzymes and size selection windows were tested and scaled in six classes of Antarctic animals (Ostracoda, Malacostraca, Bivalvia, Asteroidea, Actinopterygii, Aves). Achieving high marker density would be expensive in amphipods, the malacostracan target taxon, due to the large genome size. We propose alternative approaches such as mitogenome or target capture sequencing for this group. Pilot libraries were sequenced for all other target taxa. Ostracods, bivalves, sea stars, and fish showed overall good coverage and marker numbers for downstream population genomic analyses. In contrast, the bird test library produced low coverage and few polymorphic loci, likely due to degraded DNA.

CONCLUSIONS

Prior testing and optimization are important to identify which groups are amenable for RRS and where alternative methods may currently offer better cost-benefit ratios. The steps outlined here are easy to follow for other non-model taxa with little genomic resources, thus stimulating efficient resource use for the many pressing research questions in molecular ecology.

Genome-wide assessment elucidates connectivity and the evolutionary history of the highly dispersive marine invertebrate Littoraria flava (Littorinidae: Gastropoda)

An important goal of marine population genetics is to understand how spatial connectivity patterns are influenced by historical and evolutionary factors. In this study, we evaluate the demographic history and population structure of Littoraria flava, a highly dispersive marine gastropod in the Brazilian intertidal zone. To test the hypotheses that the species has (1) historically high levels of gene flow on a macrogeographical spatial scale and (2) a distribution in rocky shores that consists of subpopulations, we collected specimens along the Brazilian coastline and combined different sets of genetic markers (mitochondrial DNA, ITS-2 and single nucleotide polymorphisms) with niche-based modelling to predict its palaeodistribution. Low genetic structure was observed, as well as high gene flow over long distances. The demographic analyses suggest that L. flava has had periods of population bottlenecks followed by expansion. According to both palaeodistribution and coalescent simulations, these expansion events occurred during the Pleistocene interglacial cycles (21 kya) and the associated climatic changes were the probable drivers of the distribution of the species. This is the first phylogeographical study of a marine gastropod on the South American coast based on genomic markers associated with niche modelling.

Genomic library preparation and hybridization capture of formalin-fixed tissues and allozyme supernatant for population genomics and considerations for combining capture- and RADseq-based single nucleotide polymorphism data sets

Until recently many historical museum specimens were largely inaccessible to genomic inquiry, but high-throughput sequencing (HTS) approaches have allowed researchers to successfully sequence genomic DNA from dried and fluid-preserved museum specimens. In addition to preserved specimens, many museums contain large series of allozyme supernatant samples, but the amenability of these samples to HTS has not yet been assessed. Here, we compared the performance of a target-capture approach using alternative sources of genomic DNA from 10 specimens of spring salamanders (Plethodontidae: Gyrinophilus porphyriticus) collected between 1985 and 1990: allozyme supernatants, allozyme homogenate pellets and formalin-fixed tissues. We designed capture probes based on double-digest restriction-site associated sequencing (RADseq) derived loci from frozen blood samples available for seven of the specimens and assessed the success and consistency of capture and RADseq approaches. This study design enabled direct comparisons of data quality and potential biases among the different data sets for phylogenomic and population genomic analyses. We found that in phylogenetic analyses, all enrichment types for a given specimen clustered together. In principal component space all capture-based samples clustered together, but RADseq samples did not cluster with corresponding capture-based samples. Single nucleotide polymorphism calls were on average 18.3% different between enrichment types for a given individual, but these discrepancies were primarily due to differences in heterozygous/homozygous single nucleotide polymorphism calls. We demonstrate that both allozyme supernatant and formalin-fixed samples can be successfully used for population genomic analyses and we discuss ways to identify and reduce biases associated with combining capture and RADseq data.

Evolutionary drivers of sexual signal variation in Amazon Slender Anoles

Phenotypic variation among populations, as seen in the signaling traits of many species, provides an opportunity to test whether similar factors generate repeated phenotypic patterns in different parts of a species' range. We investigated whether genetic divergence, abiotic gradients, and sympatry with closely related species explain variation in the dewlap colors of Amazon Slender Anoles, Anolis fuscoauratus. To this aim, we characterized dewlap diversity in the field with respect to population genetic structure and evolutionary relationships, assessed whether dewlap phenotypes are associated with climate or landscape variables, and tested for nonrandom associations in the distributions of A. fuscoauratus phenotypes and sympatric Anolis species. We found that dewlap colors vary among but not within sites in A. fuscoauratus. Regional genetic clusters included multiple phenotypes, while populations with similar dewlaps were often distantly related. Phenotypes did not segregate in environmental space, providing no support for optimized signal transmission at a local scale. Instead, we found a negative association between certain phenotypes and sympatric Anolis species with similar dewlap color attributes, suggesting that interactions with closely related species promoted dewlap divergence among A. fuscoauratus populations. Amazon Slender Anoles emerge as a promising system to address questions about parallel trait evolution and the contribution of signaling traits to speciation.

Incomplete lineage sorting and ancient admixture, and speciation without morphological change in ghost-worm cryptic species

Morphologically similar species, that is cryptic species, may be similar or quasi-similar owing to the deceleration of morphological evolution and stasis. While the factors underlying the deceleration of morphological evolution or stasis in cryptic species remain unknown, decades of research in the field of paleontology on punctuated equilibrium have originated clear hypotheses. Species are expected to remain morphologically identical in scenarios of shared genetic variation, such as hybridization and incomplete lineage sorting, or in scenarios where bottlenecks reduce genetic variation and constrain the evolution of morphology. Here, focusing on three morphologically similar Stygocapitella species, we employ a whole-genome amplification method (WGA) coupled with double-digestion restriction-site associated DNA sequencing (ddRAD) to reconstruct the evolutionary history of the species complex. We explore population structure, use population-level statistics to determine the degree of connectivity between populations and species, and determine the most likely demographic scenarios which generally reject for recent hybridization. We find that the combination of WGA and ddRAD allowed us to obtain genomic-level data from microscopic eukaryotes (∼1 millimetre) opening up opportunities for those working with population genomics and phylogenomics in such taxa. The three species share genetic variance, likely from incomplete lineage sorting and ancient admixture. We speculate that the degree of shared variation might underlie morphological similarity in the Atlantic species complex.

Species delimitation integrative approach reveals three new species in the Nemertopsis bivittata complex

The presence of cryptic species is fairly frequent in many invertebrate groups and even more so among invertebrates with simple morphology, such as nemerteans. Consequently, the use of molecular methods for species delimitation has become a needed tool to complement morphological analyses to better recognise such species. Nemertopsis bivittata is one example of species with subtle morphological variation, but ample geographic distribution, being a good candidate for a species complex study. Here we applied two mitochondrial genes, and 2903 single nucleotide polymorphism (SNP) variants in addition to morphological characters to investigate the presence of cryptic species among specimens previously identified as N. bivittata along the Brazilian Coast. To do so, specimens were collected at 15 different sites in the north-east, south-east and southern regions. Three new species of Nemertopsis are described based on morphological and molecular analyses: Nemertopsis caete sp. nov., Nemertopsis pamelaroeae sp. nov. and Nemertopsis berthalutzae sp. nov. The species N. pamelaroeae and N. berthalutzae present broad distributions from north-east to south-east; N. caete, however, is restricted to the north-east coast. This is the first study to use this combined approach in nemerteans and shows the advantages of integrating genomic markers with classical taxonomy, and applying objective approaches to delimiting species as independently evolving entities.

Evaluating insect-host interactions as a driver of species divergence in palm flower weevils

Plants and their specialized flower visitors provide valuable insights into the evolutionary consequences of species interactions. In particular, antagonistic interactions between insects and plants have often been invoked as a major driver of diversification. Here we use a tropical community of palms and their specialized insect flower visitors to test whether antagonisms lead to higher population divergence. Interactions between palms and the insects visiting their flowers range from brood pollination to florivory and commensalism, with the latter being species that feed on decaying-and presumably undefended-plant tissues. We test the role of insect-host interactions in the early stages of diversification of nine species of beetles sharing host plants and geographical ranges by first delimiting cryptic species and then using models of genetic isolation by environment. The degree to which insect populations are structured by the genetic divergence of plant populations varies. A hierarchical model reveals that this variation is largely uncorrelated with the kind of interaction, showing that antagonistic interactions are not associated with higher genetic differentiation. Other aspects of host use that affect plant-associated insects regardless of the outcomes of their interactions, such as sensory biases, are likely more general drivers of insect population divergence.

Disjunction between canola distribution and the genetic structure of its recently described pest, the canola flower midge (Contarinia brassicola)

Population genomics is a useful tool to support integrated pest management as it can elucidate population dynamics, demography, and histories of invasion. Here, we use a restriction site-associated DNA sequencing approach combined with whole-genome amplification (WGA) to assess genomic population structure of a newly described pest of canola, the diminutive canola flower midge, Contarinia brassicola. Clustering analyses recovered little geographic structure across the main canola production region but differentiated several geographically disparate populations at edges of the agricultural zone. Given a lack of alternative hypotheses for this pattern, we suggest these data support alternative hosts for this species and thus our canola-centric view of this midge as a pest has limited our understanding of its biology. These results speak to the need for increased surveying efforts across multiple habitats and other potential hosts within Brassicaceae to improve both our ecological and evolutionary knowledge of this species and contribute to effective management strategies. We additionally found that use of WGA prior to library preparation was an effective method for increasing DNA quantity of these small insects prior to restriction site-associated DNA sequencing and had no discernible impact on genotyping consistency for population genetic analysis; WGA is therefore likely to be tractable for other similar studies that seek to randomly sample markers across the genome in small organisms.

Impacts of the Toba eruption and montane forest expansion on diversification in Sumatran parachuting frogs (Rhacophorus)

Catastrophic events, such as volcanic eruptions, can have profound impacts on the demographic histories of resident taxa. Due to its presumed effect on biodiversity, the Pleistocene eruption of super-volcano Toba has received abundant attention. We test the effects of the Toba eruption on the diversification, genetic diversity, and demography of three co-distributed species of parachuting frogs (Genus Rhacophorus) on Sumatra. We generate target-capture data (~950 loci and ~440,000 bp) for three species of parachuting frogs and use these data paired with previously generated double digest restriction-site associated DNA (ddRADseq) data to estimate population structure and genetic diversity, to test for population size changes using demographic modelling, and to estimate the temporal clustering of size change events using a full-likelihood Bayesian method. We find that populations around Toba exhibit reduced genetic diversity compared with southern populations, and that northern populations exhibit a shift in effective population size around the time of the eruption (~80 kya). However, we infer a stronger signal of expansion in southern populations around ~400 kya, and at least two of the northern populations may have also expanded at this time. Taken together, these findings suggest that the Toba eruption precipitated population declines in northern populations, but that the demographic history of these three species was also strongly impacted by mid-Pleistocene forest expansion during glacial periods. We propose local rather than regional effects of the Toba eruption, and emphasize the dynamic nature of diversification on the Sunda Shelf.

Thitarodes shambalaensis sp. nov. (Lepidoptera, Hepialidae): a new host of the caterpillar fungus Ophiocordyceps sinensis supported by genome-wide SNP data

A new species of ghost moth, Thitarodesshambalaensissp. nov., is described from Yanzigou glacier, Mt. Gongga, Sichuan, China. The species is a host of the economically important caterpillar fungus Ophiocordycepssinensis. Establishment of this new species is supported by morphology and genetic differentiation measured in a CO1 phylogeny and in genome-wide SNP coverage. A summary tree from 538 sequences of different genetic markers from Thitarodes (including sequences extracted from caterpillar fungus sclerotium samples) support the genus Thitarodes as a monophyletic group, and indicate that Thitarodes is the host genus for O.sinensis. Sampling efforts so far have centered on half of the known phylogenetic diversity of Thitarodes, with some species-level clusters (separated by < 2.5% genetic distance) containing 17 described species. Fifteen clusters are known from either a single “orphan taxon” or a single sequence from a caterpillar fungus sclerotium sample. We provide suggestions for building a more robust phylogeny of the genus Thitarodes and highlight some of the conservation threats that species from this genus face due to unprecedented habitat exploitation.

Morphological and Molecular Perspectives on the Phylogeny, Evolution, and Classification of Weevils (Coleoptera: Curculionoidea): Proceedings from the 2016 International Weevil Meeting

The 2016 International Weevil Meeting was held immediately after the International Congress of Entomology (ICE). It built on the topics and content of the 2016 ICE weevil symposium Phylogeny and Evolution of Weevils (Coleoptera: Curculionoidea): A Symposium in Honor of Dr. Guillermo "Willy” Kuschel. Beyond catalyzing research and collaboration, the meeting was intended to serve as a forum for identifying priorities and goals for those who study weevils. The meeting consisted of 46 invited and contributed lectures, discussion sessions and introductory remarks presented by 23 speakers along with eight contributed research posters. These were organized into three convened sessions, each lasting one day: (1) weevil morphology; (2) weevil fossils, biogeography and host/habitat associations; and (3) molecular phylogenetics and classification of weevils. Some of the topics covered included the 1K Weevils Project, major morphological character systems of adult and larval weevils, weevil morphological terminology, prospects for future morphological character discovery, phylogenetic analysis of morphological character data, the current status of weevil molecular phylogenetics and evolution, resources available for phylogenetic and comparative genomic studies of weevils, the weevil fossil record, weevil biogeography and evolution, weevil host plants, evolutionary development of the weevil rostrum, resources available for weevil identification and the current status of and challenges in weevil classification.