Museum Genomics Confirms that the Lord Howe Island Stick Insect Survived Extinction

Genome skimming for improved phylogenetics of Taiwanese phasmids (Insecta: Phasmatodea)

Taiwan, a relatively young continental island, harbors a high proportion of endemic phasmids, reflecting its unique evolutionary history. However, a comprehensive phylogenetic framework to clarify these phasmids is still lacking. In this study, we sequenced ten of eleven valid genera and two undescribed species of Taiwanese phasmids (total 16 species) using the genome-skimming approach. We also integrated these sequences with public databases to create two aligned datasets: one comprising 92 taxa (mitogenomes) and the other 606 taxa (seven nuclear and mitochondrial genes), enabling us to examine their phylogenetic relationships using longer sequences and more samples. Our analyses show that Taiwanese phasmids should be categorized into six families, with a revised number of genera to 13. Furthermore, four species require taxonomic treatments: namely Micadina honei (Günther, 1940) comb. nov., Micadina truncatum (Shiraki, 1935) comb. nov., Otraleus okunii (Shiraki, 1935) comb. nov., and Ramulus granulatus (Shiraki, 1935) syn. nov. now recognized as Ramulus artemis (Westwood, 1859). While some Taiwanese genera exhibit polyphyletic relationships, our findings highlight the importance of taxon sampling, particularly for type species in resolving these systematic issues. The genome-skimming approach has proven to be an excellent method for producing comparable sequence datasets, facilitating the investigation of highly diverse insects, even when samples are old, small, or have highly fragmented DNAs.

Purging of Highly Deleterious Mutations Through an Extreme Bottleneck

Transitions to captivity often produce population bottlenecks. On the one hand, bottlenecks increase inbreeding and decrease effective population size, thus increasing extinction risk. On the other hand, elevated homozygosity associated with inbreeding may purge deleterious mutations. Previous empirical studies of purging in captive breeding programs have focused on phenotypic measurements. We test natural selection's ability to purge deleterious mutations following an extreme population bottleneck by analyzing patterns of genetic diversity in wild and captive-bred individuals of the Lord Howe Island stick insect, Dryococelus australis. Dryococelus australis has been bred in captivity for two decades, having passed through an extreme bottleneck-only two mating pairs with few new additions since then. The magnitude of the bottleneck together with high female fecundity but low offspring recruitment set up nearly ideal conditions for the purging of deleterious mutations. As expected, captive-bred individuals had a greater number of long runs of homozygosity compared with wild individuals, implying strong inbreeding in captivity which would facilitate purging in homozygous regions. Stop-codon mutations were preferentially depleted in captivity compared with other mutations in coding and noncoding regions. The more deleterious a mutation was predicted to be, the more likely it was found outside of runs of homozygosity, implying that inbreeding facilitates the expression and thus removal of deleterious mutations, even after such an extreme bottleneck and under the benign conditions of captivity. These data implicate inbreeding and recessive deleterious mutation load in fitness variation among captive and wild D. australis.

Developing the Protocol Infrastructure for DNA Sequencing Natural History Collections

Intentionally preserved biological material in natural history collections represents a vast repository of biodiversity. Advances in laboratory and sequencing technologies have made these specimens increasingly accessible for genomic analyses, offering a window into the genetic past of species and often permitting access to information that can no longer be sampled in the wild. Due to their age, preparation and storage conditions, DNA retrieved from museum and herbarium specimens is often poor in yield, heavily fragmented and biochemically modified. This not only poses methodological challenges in recovering nucleotide sequences, but also makes such investigations susceptible to environmental and laboratory contamination. In this paper, we review the practical challenges associated with making the recovery of DNA sequence data from museum collections more routine. We first review key operational principles and issues to address, to guide the decision-making process and dialogue between researchers and curators about when and how to sample museum specimens for genomic analyses. We then outline the range of steps that can be taken to reduce the likelihood of contamination including laboratory set-ups, workflows and working practices. We finish by presenting a series of case studies, each focusing on protocol practicalities for the application of different mainstream methodologies to museum specimens including: (i) shotgun sequencing of insect mitogenomes, (ii) whole genome sequencing of insects, (iii) genome skimming to recover plant plastid genomes from herbarium specimens, (iv) target capture of multi-locus nuclear sequences from herbarium specimens, (v) RAD-sequencing of bird specimens and (vi) shotgun sequencing of ancient bovid bone samples.

Genomes from historical Drosophila melanogaster specimens illuminate adaptive and demographic changes across more than 200 years of evolution

The ability to perform genomic sequencing on long-dead organisms is opening new frontiers in evolutionary research. These opportunities are especially notable in the case of museum collections, from which countless documented specimens may now be suitable for genomic analysis-if data of sufficient quality can be obtained. Here, we report 25 newly sequenced genomes from museum specimens of the model organism Drosophila melanogaster, including the oldest extant specimens of this species. By comparing historical samples ranging from the early 1800s to 1933 against modern-day genomes, we document evolution across thousands of generations, including time periods that encompass the species' initial occupation of northern Europe and an era of rapidly increasing human activity. We also find that the Lund, Sweden population underwent local genetic differentiation during the early 1800s to 1933 interval (potentially due to drift in a small population) but then became more similar to other European populations thereafter (potentially due to increased migration). Within each century-scale time period, our temporal sampling allows us to document compelling candidates for recent natural selection. In some cases, we gain insights regarding previously implicated selection candidates, such as ChKov1, for which our inferred timing of selection favors the hypothesis of antiviral resistance over insecticide resistance. Other candidates are novel, such as the circadian-related gene Ahcy, which yields a selection signal that rivals that of the DDT resistance gene Cyp6g1. These insights deepen our understanding of recent evolution in a model system, and highlight the potential of future museomic studies.

The oral repellent - science fiction or common sense? Insects, vector-borne diseases, failing strategies, and a bold proposition

Over the last decades, unimaginable amounts of money have gone into research and development of vector control measures, repellents, treatment, and vaccines for vector borne diseases. Technological progress and scientific breakthroughs allowed for ever more sophisticated and futuristic strategies. Yet, each year, millions of people still die or suffer from potentially serious consequences of malaria or dengue to more recent infections, such as zika or chikungunya, or of debilitating consequences of neglected tropical diseases. This does not seem value for money. In addition, all current vector control strategies and personal protection methods have shortcomings, some serious, that are either destructive to non-target species or unsatisfactory in their effectiveness. On the other hand, the rapid decline in insect populations and their predators reflects decades-long aggressive and indiscriminate vector control. This major disruption of biodiversity has an impact on human life not anticipated by the well-meaning killing of invertebrates. The objective of this paper is to re-examine current control methods, their effectiveness, their impact on biodiversity, human and animal health, and to call for scientific courage in the pursuit of fresh ideas. This paper brings together topics that are usually presented in isolation, thereby missing important links that offer potential solutions to long-standing problems in global health. First, it serves as a reminder of the importance of insects to human life and discusses the few that play a role in transmitting disease. Next, it examines critically the many currently employed vector control strategies and personal protection methods. Finally, based on new insights into insect chemo-sensation and attractants, this perspective makes a case for revisiting a previously abandoned idea, the oral repellent, and its use via currently successful methods of mass-application. The call is out for focused research to provide a powerful tool for public health, tropical medicine, and travel medicine.

Genomics, Population Divergence, and Historical Demography of the World's Largest and Endangered Butterfly, The Queen Alexandra's Birdwing

The world's largest butterfly is the microendemic Papua New Guinean Ornithoptera alexandrae. Despite years of conservation efforts to protect its habitat and breed this up-to-28-cm butterfly, this species still figures as endangered in the IUCN Red List and is only known from two allopatric populations occupying a total of only ∼140 km². Here we aim at assembling reference genomes for this species to investigate its genomic diversity, historical demography and determine whether the population is structured, which could provide guidance for conservation programs attempting to (inter)breed the two populations. Using a combination of long and short DNA reads and RNA sequencing, we assembled six reference genomes of the tribe Troidini, with four annotated genomes of O. alexandrae and two genomes of related species Ornithoptera priamus and Troides oblongomaculatus. We estimated the genomic diversity of the three species, and we proposed scenarios for the historical population demography using two polymorphism-based methods taking into account the characteristics of low-polymorphic invertebrates. Indeed, chromosome-scale assemblies reveal very low levels of nuclear heterozygosity across Troidini, which appears to be exceptionally low for O. alexandrae (lower than 0.01%). Demographic analyses demonstrate low and steadily declining Ne throughout O. alexandrae history, with a divergence into two distinct populations about 10,000 years ago. These results suggest that O. alexandrae distribution has been microendemic for a long time. It should also make local conservation programs aware of the genomic divergence of the two populations, which should not be ignored if any attempt is made to cross the two populations.

Maximizing Molecular Data From Low-Quality Fluid-Preserved Specimens in Natural History Collections

Over the past decade, museum genomics studies have focused on obtaining DNA of sufficient quality and quantity for sequencing from fluid-preserved natural history specimens, primarily to be used in systematic studies. While these studies have opened windows to evolutionary and biodiversity knowledge of many species worldwide, published works often focus on the success of these DNA sequencing efforts, which is undoubtedly less common than obtaining minimal or sometimes no DNA or unusable sequence data from specimens in natural history collections. Here, we attempt to obtain and sequence DNA extracts from 115 fresh and 41 degraded samples of homalopsid snakes, as well as from two degraded samples of a poorly known snake, Hydrablabes periops. Hydrablabes has been suggested to belong to at least two different families (Natricidae and Homalopsidae) and with no fresh tissues known to be available, intractable museum specimens currently provide the only opportunity to determine this snake’s taxonomic affinity. Although our aim was to generate a target-capture dataset for these samples, to be included in a broader phylogenetic study, results were less than ideal due to large amounts of missing data, especially using the same downstream methods as with standard, high-quality samples. However, rather than discount results entirely, we used mapping methods with references and pseudoreferences, along with phylogenetic analyses, to maximize any usable molecular data from our sequencing efforts, identify the taxonomic affinity of H. periops, and compare sequencing success between fresh and degraded tissue samples. This resulted in largely complete mitochondrial genomes for five specimens and hundreds to thousands of nuclear loci (ultra-conserved loci, anchored-hybrid enrichment loci, and a variety of loci frequently used in squamate phylogenetic studies) from fluid-preserved snakes, including a specimen of H. periops from the Field Museum of Natural History collection. We combined our H. periops data with previously published genomic and Sanger-sequenced datasets to confirm the familial designation of this taxon, reject previous taxonomic hypotheses, and make biogeographic inferences for Hydrablabes. A second H. periops specimen, despite being seemingly similar for initial raw sequencing results and after being put through the same protocols, resulted in little usable molecular data. We discuss the successes and failures of using different pipelines and methods to maximize the products from these data and provide expectations for others who are looking to use DNA sequencing efforts on specimens that likely have degraded DNA.

Life Science Identifier (Hydrablabes periops)

urn:lsid:zoobank.org:pub:F2AA44 E2-D2EF-4747-972A-652C34C2C09D.

Assessing species coverage and assembly quality of rapidly accumulating sequenced genomes

BACKGROUND

Ambitious initiatives to coordinate genome sequencing of Earth's biodiversity mean that the accumulation of genomic data is growing rapidly. In addition to cataloguing biodiversity, these data provide the basis for understanding biological function and evolution. Accurate and complete genome assemblies offer a comprehensive and reliable foundation upon which to advance our understanding of organismal biology at genetic, species, and ecosystem levels. However, ever-changing sequencing technologies and analysis methods mean that available data are often heterogeneous in quality. To guide forthcoming genome generation efforts and promote efficient prioritization of resources, it is thus essential to define and monitor taxonomic coverage and quality of the data.

FINDINGS

Here we present an automated analysis workflow that surveys genome assemblies from the United States NCBI, assesses their completeness using the relevant BUSCO datasets, and collates the results into an interactively browsable resource. We apply our workflow to produce a community resource of available assemblies from the phylum Arthropoda, the Arthropoda Assembly Assessment Catalogue. Using this resource, we survey current taxonomic coverage and assembly quality at the NCBI, examine how key assembly metrics relate to gene content completeness, and compare results from using different BUSCO lineage datasets.

CONCLUSIONS

These results demonstrate how the workflow can be used to build a community resource that enables large-scale assessments to survey species coverage and data quality of available genome assemblies, and to guide prioritizations for ongoing and future sampling, sequencing, and genome generation initiatives.

The phylogenic position of aschiphasmatidae in euphasmatodea based on mitochondrial genomic evidence

The mitochondrial genome (mitogenome) has been regarded as significant source of data to better understand the phylogenetic relationships within the Euphasmatodea, but no mitogenome in Aschiphasmatoidea has been sequenced to date. In this study, two mitogenomes of Orthomeria smaragdinum and Nanhuaphasma hamicercum of Aschiphasmatidae were sequenced and annotated for the first time. The same mitochondrial gene rearrangement structure was present in the two mitogenomes sequenced, showing as the translocation of tRNA-Arg and tRNA-Asn, which conformed to the tandem duplication-random loss and could be used as a possible synapomorphy for Aschiphasmatidae. The phylogenetic results based on the maximum likelihood (ML) and bayesian inference (BI) methods both showed that Aschiphasmatidae and Neophasmatodea in Euphasmatodea are sister taxa. Although the monophyly of Oriophasmata, Occidophasmata, Diapheromeridae, Phasmatidae, Lonchodidae and Bacilloidea has not been solved, the monophyly of Neophasmatodea and Phyllioidea was well supported.

Genotyping-in-Thousands by sequencing of archival fish scales reveals maintenance of genetic variation following a severe demographic contraction in kokanee salmon

Historical DNA analysis of archival samples has added new dimensions to population genetic studies, enabling spatiotemporal approaches for reconstructing population history and informing conservation management. Here we tested the efficacy of Genotyping-in-Thousands by sequencing (GT-seq) for collecting targeted single nucleotide polymorphism genotypic data from archival scale samples, and applied this approach to a study of kokanee salmon (Oncorhynchus nerka) in Kluane National Park and Reserve (KNPR; Yukon, Canada) that underwent a severe 12-year population decline followed by a rapid rebound. We genotyped archival scales sampled pre-crash and contemporary fin clips collected post-crash, revealing high coverage (> 90% average genotyping across all individuals) and low genotyping error (< 0.01% within-libraries, 0.60% among-libraries) despite the relatively poor quality of recovered DNA. We observed slight decreases in expected heterozygosity, allelic diversity, and effective population size post-crash, but none were significant, suggesting genetic diversity was retained despite the severe demographic contraction. Genotypic data also revealed the genetic distinctiveness of a now extirpated population just outside of KNPR, revealing biodiversity loss at the northern edge of the species distribution. More broadly, we demonstrated GT-seq as a valuable tool for collecting genome-wide data from archival samples to address basic questions in ecology and evolution, and inform applied research in wildlife conservation and fisheries management.

The mitochondrial genome of Eurycantha calcarata Lucas, 1869 (Phasmatodea: Lonchodinae) and its phylogeny

The Lonchodinae (Phasmatodea: Phasmatidae) is rich in insect species with more than 330 species of 40 genera. The phylogenetic relationships within Lonchodinae have been under debate. We successfully sequenced the complete mitogenome of Eurycantha calcarata Lucas, 1869 (Phasmatodea: Lonchodinae) with a length of 16,280 bp, which had the same genes and gene arrangements as those of various published papers on stick insects. The whole mitogenome and control region of E. calcarata had a high AT content of 78.2 and 85.9%, respectively. All PCGs used ATN as the start codon, and most PCGs used TAA/TAG as the stop codons excluding COX2 (T), COX3 (TA), and ND5 (TA). To discuss the phylogeny of Lonchodinae, we reconstructed the phylogenetic relationships of 27 species of Phasmatodea including E. calcarata and two species of Embioptera used as outgroups. In BI and ML trees, the monophyly of Lonchodinae and Necrosciinae was well supported, whereas the monophyly of Clitumninae was not recovered. These results indicated that Lonchodinae was a sister clade to Phylliinae and E. calcarata was a sister clade to Phraortes genus.

Three Complete Mitochondrial Genomes of Orestes guangxiensis, Peruphasma schultei, and Phryganistria guangxiensis (Insecta: Phasmatodea) and Their Phylogeny

Insects of the order Phasmatodea are mainly distributed in the tropics and subtropics and are best known for their remarkable camouflage as plants. In this study, we sequenced three complete mitochondrial genomes from three different families: Orestes guangxiensis, Peruphasma schultei, and Phryganistria guangxiensis. The lengths of the three mitochondrial genomes were 15,896 bp, 16,869 bp, and 17,005 bp, respectively, and the gene composition and structure of the three stick insects were identical to those of the most recent common ancestor of insects. The phylogenetic relationships among stick insects have been chaotic for a long time. In order to discuss the intra- and inter-ordinal relationship of Phasmatodea, we used the 13 protein-coding genes (PCGs) of 85 species for maximum likelihood (ML) and Bayesian inference (BI) analyses. Results showed that the internal topological structure of Phasmatodea had a few differences in both ML and BI trees and long-branch attraction (LBA) appeared between Embioptera and Zoraptera, which led to a non-monophyletic Phasmatodea. Consequently, after removal of the Embioptera and Zoraptera species, we re-performed ML and BI analyses with the remaining 81 species, which showed identical topology except for the position of Tectarchus ovobessus (Phasmatodea). We recovered the monophyly of Phasmatodea and the sister-group relationship between Phasmatodea and Mantophasmatodea. Our analyses also recovered the monophyly of Heteropterygidae and the paraphyly of Diapheromeridae, Phasmatidae, Lonchodidae, Lonchodinae, and Clitumninae. In this study, Peruphasma schultei (Pseudophasmatidae), Phraortes sp. YW-2014 (Lonchodidae), and species of Diapheromeridae clustered into the clade of Phasmatidae. Within Heteropterygidae, O. guangxiensis was the sister clade to O. mouhotii belonging to Dataminae, and the relationship of (Heteropteryginae + (Dataminae + Obriminae)) was recovered.

Museomics dissects the genetic basis for adaptive seasonal colouration in the least weasel

Dissecting the link between genetic variation and adaptive phenotypes provides outstanding opportunities to understand fundamental evolutionary processes. Here, we use a museomics approach to investigate the genetic basis and evolution of winter coat colouration morphs in least weasels (Mustela nivalis), a repeated adaptation for camouflage in mammals with seasonal pelage colour moults across regions with varying winter snow. Whole-genome sequence data was obtained from biological collections and mapped onto a newly assembled reference genome for the species. Sampling represented two replicate transition zones between nivalis and vulgaris colouration morphs in Europe, which typically develop white or brown winter coats, respectively. Population analyses showed that the morph distribution across transition zones is not a by-product of historical structure. Association scans linked a 200 kb genomic region to colouration morph, which was validated by genotyping museum specimens from inter-morph experimental crosses. Genotyping the wild populations narrowed down the association to pigmentation gene MC1R and pinpointed a candidate amino acid change co-segregating with colouration morph. This polymorphism replaces an ancestral leucine residue by lysine at the start of the first extracellular loop of the protein in the vulgaris morph. A selective sweep signature overlapped the association region in vulgaris, suggesting that past adaptation favoured winter-brown morphs and can anchor future adaptive responses to decreasing winter snow. Using biological collections as valuable resources to study natural adaptations, our study showed a new evolutionary route generating winter colour variation in mammals and that seasonal camouflage can be modulated by changes at single key genes.

The mitochondrial genome of a leaf insect Phyllium westwoodii (Phasmatodea: Phylliidae) in Southeast Asia

The nearly complete mitochondrial genome (mitogenome) of Phyllium westwoodii, a typical leaf mimic insect in Phasmatodea, was obtained in this study. This mitogenome is 17,222 bp in length and contains 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and almost complete control regions. All PCGs initiate with 'ATN' except for NAD4L that uses 'TTG' as the start codon, and terminate with 'TAA' except for COX2 that uses a single 'T' residue as the stop codon. The phylogenetic analysis based on the concatenated sequences of 13 PCGs and two rRNAs shows that P. westwoodii is closer to Phyllium tibetense than Phyllium giganteum.

A SMRT approach for targeted amplicon sequencing of museum specimens (Lepidoptera)-patterns of nucleotide misincorporation

Natural history collections are a valuable resource for molecular taxonomic studies and for examining patterns of evolutionary diversification, particularly in the case of rare or extinct species. However, the recovery of sequence information is often complicated by DNA degradation. This article describes use of the Sequel platform (Pacific Biosciences) to recover the 658 bp barcode region of the mitochondrial cytochrome c oxidase I (COI) gene from 380 butterflies with an average age of 50 years. Nested multiplex PCR was employed for library preparation to facilitate sequence recovery from extracts with low concentrations of highly degraded DNA. By employing circular consensus sequencing (CCS) of short amplicons (circa 150 bp), full-length barcodes could be assembled without a reference sequence, an important advance from earlier protocols which required reference sequences to guide contig assembly. The Sequel protocol recovered COI sequences (499 bp on average) from 318 of 380 specimens (84%), much higher than for Sanger sequencing (26%). Because each read derives from a single molecule, it was also possible to quantify the incidence of substitutions arising from DNA damage. In agreement with past work on sequence changes induced by DNA degradation, the transition C/G → T/A was the most prevalent category of change, but its rate of occurrence (4.58E−4) was so low that it did not impede the recovery of reliable sequences. Because the current protocol recovers COI sequence from most museum specimens, and because sequence fidelity is unaffected by nucleotide misincorporations, large-scale sequence characterization of museum specimens is feasible.

Positively selected variants in functionally important regions of TLR7 in Alouatta guariba clamitans with yellow fever virus exposure in Northern Argentina

BACKGROUND

In 2007-2009, a major yellow fever virus (YFV) outbreak in Northern Argentina decimated the local howler monkey (Alouatta) population.

AIMS

To evaluate whether the surviving howler monkeys possess advantageous genetic variants inherited from monkeys alive prior to the YFV outbreak, we explored the relationship between Toll-like receptor (TLR) 7 and TLR8 gene variation and YFV susceptibility.

METHODS

We used samples from Alouatta individuals in Misiones, Argentina alive before the YFV outbreak, individuals that died during the outbreak, and individuals that survived the outbreak and are alive today. We measured genetic divergence between Alouatta YFV exposure groups and evaluated Alouatta-specific substitutions for functional consequences.

RESULTS

We did not find different allele frequencies in the post-YFV exposure Alouatta group compared to the pre-exposure group. We identified three nonsynonymous variants in TLR7 in Alouatta guariba clamitans. Two of these substitutions are under positive selection in functionally important regions of the gene.

DISCUSSION AND CONCLUSIONS

Our results did not indicate that surviving howler monkey spossess advantageous genetic variants at greater frequency than those alive before the YFV outbreak. However, the positively selected unique coding differences in A. guariba clamitans are in the region important in pathogen detection which may affect YFV resistance. Morework is necessary to fully explore this hypothesis.

Effects of sample age on data quality from targeted sequencing of museum specimens: what are we capturing in time?

Background

Next generation sequencing (NGS) can recover DNA data from valuable extant and extinct museum specimens. However, archived or preserved DNA is difficult to sequence because of its fragmented, damaged nature, such that the most successful NGS methods for preserved specimens remain sub-optimal. Improving wet-lab protocols and comprehensively determining the effects of sample age on NGS library quality are therefore of vital importance. Here, I examine the relationship between sample age and several indicators of library quality following targeted NGS sequencing of ~ 1300 loci using 271 samples of pinned moth specimens (Helicoverpa armigera) ranging in age from 5 to 117 years.

Results

I find that older samples have lower DNA concentrations following extraction and thus require a higher number of indexing PCR cycles during library preparation. When sequenced reads are aligned to a reference genome or to only the targeted region, older samples have a lower number of sequenced and mapped reads, lower mean coverage, and lower estimated library sizes, while the percentage of adapters in sequenced reads increases significantly as samples become older. Older samples also show the poorest capture success, with lower enrichment and a higher improved coverage anticipated from further sequencing.

Conclusions

Sample age has significant, measurable impacts on the quality of NGS data following targeted enrichment. However, incorporating a uracil-removing enzyme into the blunt end-repair step during library preparation could help to repair DNA damage, and using a method that prevents adapter-dimer formation may result in improved data yields.

Utilizing field collected insects for next generation sequencing: Effects of sampling, storage, and DNA extraction methods

DNA sequencing technologies continue to advance the biological sciences, expanding opportunities for genomic studies of non-model organisms for basic and applied questions. Despite these opportunities, many next generation sequencing protocols have been developed assuming a substantial quantity of high molecular weight DNA (>100 ng), which can be difficult to obtain for many study systems. In particular, the ability to sequence field-collected specimens that exhibit varying levels of DNA degradation remains largely unexplored. In this study we investigate the influence of five traditional insect capture and curation methods on Double-Digest Restriction Enzyme Associated DNA (ddRAD) sequencing success for three wild bee species. We sequenced a total of 105 specimens (between 7-13 specimens per species and treatment). We additionally investigated how different DNA quality metrics (including pre-sequence concentration and contamination) predicted downstream sequencing success, and also compared two DNA extraction methods. We report successful library preparation for all specimens, with all treatments and extraction methods producing enough highly reliable loci for population genetic analyses. Although results varied between species, we found that specimens collected by net sampling directly into 100% EtOH, or by passive trapping followed by 100% EtOH storage before pinning tended to produce higher quality ddRAD assemblies, likely as a result of rapid specimen desiccation. Surprisingly, we found that specimens preserved in propylene glycol during field sampling exhibited lower-quality assemblies. We provide recommendations for each treatment, extraction method, and DNA quality assessment, and further encourage researchers to consider utilizing a wider variety of specimens for genomic analyses.

Genomic impact of severe population decline in a nomadic songbird

Uncovering the population genetic histories of non-model organisms is increasingly possible through advances in next generation sequencing and DNA sampling of museum specimens. This new information can inform conservation of threatened species, particularly those for which historical and contemporary population data are unavailable or challenging to obtain. The critically endangered, nomadic regent honeyeater Anthochaera phrygia was abundant and widespread throughout south-eastern Australia prior to a rapid population decline and range contraction since the 1970s. A current estimated population of 250-400 individuals is distributed sparsely across 600,000 km2 from northern Victoria to southern Queensland. Using hybridization RAD (hyRAD) techniques, we obtained a SNP dataset from 64 museum specimens (date 1879-1960), 102 'recent' (1989-2012) and 52 'current' (2015-2016) wild birds sampled throughout the historical and contemporary range. We aimed to estimate population genetic structure, genetic diversity and population size of the regent honeyeater prior to its rapid decline. We then assessed the impact of the decline on recent and current population size, structure and genetic diversity. Museum sampling showed population structure in regent honeyeaters was historically low, which remains the case despite a severe fragmentation of the breeding range. Population decline has led to minimal loss of genetic diversity since the 1980's. Capacity to quantify the overall magnitude of both genetic diversity loss and population decline was limited by the poorer quality of genomic data derived from museum specimens. A rapid population decline, coupled with the regent honeyeater's high mobility, means a detectable genomic impact of this decline has not yet manifested. Extinction may occur in this nomadic species before a detectable genomic impact of small population size is realised. We discuss the implications for genetic management of endangered mobile species and enhancing the value of museum specimens in population genomic studies.

Complete mitochondrial genomes from transcriptomes: assessing pros and cons of data mining for assembling new mitogenomes

Thousands of eukaryotes transcriptomes have been generated, mainly to investigate nuclear genes expression, and the amount of available data is constantly increasing. A neglected but promising use of this large amount of data is to assemble organelle genomes. To assess the reliability of this approach, we attempted to reconstruct complete mitochondrial genomes from RNA-Seq experiments of Reticulitermes termite species, for which transcriptomes and conspecific mitogenomes are available. We successfully assembled complete molecules, although a few gaps corresponding to tRNAs had to be filled manually. We also reconstructed, for the first time, the mitogenome of Reticulitermes banyulensis. The accuracy and completeness of mitogenomes reconstruction appeared independent from transcriptome size, read length and sequencing design (single/paired end), and using reference genomes from congeneric or intra-familial taxa did not significantly affect the assembly. Transcriptome-derived mitogenomes were found highly similar to the conspecific ones obtained from genome sequencing (nucleotide divergence ranging from 0% to 3.5%) and yielded a congruent phylogenetic tree. Reads from contaminants and nuclear transcripts, although slowing down the process, did not result in chimeric sequence reconstruction. We suggest that the described approach has the potential to increase the number of available mitogenomes by exploiting the rapidly increasing number of transcriptomes.

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.

Mitochondrial genome of Phalantus geniculatus (Hemiptera: Reduviidae): trnT duplication and phylogenetic implications

Reduviidae is the second largest family of Heteroptera and most of them are important natural enemies of agricultural and forest pests. Most of the sequenced mitochondrial (mt) genomes in this family have the typical gene arrangement of insects and encode 37 coding genes (13 protein-coding genes, 22 tRNA genes and two rRNA genes). In the present study, we sequenced the mt genome of Phalantus geniculatus from the subfamily Peiratinae through high-throughput sequencing and encountered the duplication of tRNA genes for the first time in this subfamily. We identified 23 tRNA genes, including 22 tRNAs commonly found in insect mt genomes and an extra trnT (trnT2), which has high sequence similarity (96.9%) to trnT1. The presence of a "pseudo-trnP" in the non-coding region between trnT1 and trnT2 supports the hypothesis that the presence of an extra trnT can be explained by the tandem duplication-random loss (TDRL) model. Phylogenetic results inferred from mt genome sequences supported a sister relationship between Phymatinae and the remaining sampled subfamilies, as well as a paraphyletic Reduviinae. The present study highlights the utility of mt genomes in the phylogenetic study of Reduviidae based on the large scale taxon sampling in the future.

The Draft Genome of the Invasive Walking Stick, Medauroidea extradendata, Reveals Extensive Lineage-Specific Gene Family Expansions of Cell Wall Degrading Enzymes in Phasmatodea

Plant cell wall components are the most abundant macromolecules on Earth. The study of the breakdown of these molecules is thus a central question in biology. Surprisingly, plant cell wall breakdown by herbivores is relatively poorly understood, as nearly all early work focused on the mechanisms used by symbiotic microbes to breakdown plant cell walls in insects such as termites. Recently, however, it has been shown that many organisms make endogenous cellulases. Insects, and other arthropods, in particular have been shown to express a variety of plant cell wall degrading enzymes in many gene families with the ability to break down all the major components of the plant cell wall. Here we report the genome of a walking stick, Medauroidea extradentata, an obligate herbivore that makes uses of endogenously produced plant cell wall degrading enzymes. We present a draft of the 3.3Gbp genome along with an official gene set that contains a diversity of plant cell wall degrading enzymes. We show that at least one of the major families of plant cell wall degrading enzymes, the pectinases, have undergone a striking lineage-specific gene family expansion in the Phasmatodea. This genome will be a useful resource for comparative evolutionary studies with herbivores in many other clades and will help elucidate the mechanisms by which metazoans breakdown plant cell wall components.

Assembling large genomes: analysis of the stick insect (Clitarchus hookeri) genome reveals a high repeat content and sex-biased genes associated with reproduction

BACKGROUND

Stick insects (Phasmatodea) have a high incidence of parthenogenesis and other alternative reproductive strategies, yet the genetic basis of reproduction is poorly understood. Phasmatodea includes nearly 3000 species, yet only the genome of Timema cristinae has been published to date. Clitarchus hookeri is a geographical parthenogenetic stick insect distributed across New Zealand. Sexual reproduction dominates in northern habitats but is replaced by parthenogenesis in the south. Here, we present a de novo genome assembly of a female C. hookeri and use it to detect candidate genes associated with gamete production and development in females and males. We also explore the factors underlying large genome size in stick insects.

RESULTS

The C. hookeri genome assembly was 4.2 Gb, similar to the flow cytometry estimate, making it the second largest insect genome sequenced and assembled to date. Like the large genome of Locusta migratoria, the genome of C. hookeri is also highly repetitive and the predicted gene models are much longer than those from most other sequenced insect genomes, largely due to longer introns. Miniature inverted repeat transposable elements (MITEs), absent in the much smaller T. cristinae genome, is the most abundant repeat type in the C. hookeri genome assembly. Mapping RNA-Seq reads from female and male gonadal transcriptomes onto the genome assembly resulted in the identification of 39,940 gene loci, 15.8% and 37.6% of which showed female-biased and male-biased expression, respectively. The genes that were over-expressed in females were mostly associated with molecular transportation, developmental process, oocyte growth and reproductive process; whereas, the male-biased genes were enriched in rhythmic process, molecular transducer activity and synapse. Several genes involved in the juvenile hormone synthesis pathway were also identified.

CONCLUSIONS

The evolution of large insect genomes such as L. migratoria and C. hookeri genomes is most likely due to the accumulation of repetitive regions and intron elongation. MITEs contributed significantly to the growth of C. hookeri genome size yet are surprisingly absent from the T. cristinae genome. Sex-biased genes identified from gonadal tissues, including genes involved in juvenile hormone synthesis, provide interesting candidates for the further study of flexible reproduction in stick insects.

Conservation Biology: A Walking Stick's Redux on Lord Howe Island

Arguably the world's rarest insect - the Lord Howe Island 'tree lobster' - is being brought back from the brink. A recent study has confirmed the identity of this species using genomic data, which backstops its reintroduction to this World Heritage listed oceanic island.