Targeted exon sequencing by in-solution hybrid selection

Overview of Admixture Mapping

Admixture mapping is a powerful method of gene mapping for diseases or traits that show differential risk by ancestry. Admixture mapping has been applied most often to Americans who trace ancestry to various combinations of Native Americans, Europeans, and West Africans. Recent developments in admixture mapping include improvements in the methods and reference data needed to make inferences about ancestry as well as extensions of the mapping approach in the framework of linear mixed models. This overview outlines the key concepts of admixture mapping. It describes approaches for inferring local ancestry and provides strategies for performing admixture mapping depending on the study design. Finally, linkage analysis, association analysis, and admixture mapping are compared, with an emphasis on integrating admixture mapping and association testing. Published 2023. This article is a U.S. Government work and is in the public domain in the USA.

Phylogenomics and historical biogeography of the cleptoparasitic bee genus Nomada (Hymenoptera: Apidae) using ultraconserved elements

The genus Nomada Scopoli (Hymenoptera: Apidae) is the largest genus of brood parasitic bees with nearly 800 species found across the globe and in nearly all biogeographic realms except Antarctica. There is no previous molecular phylogeny focused on Nomada despite their high species abundance nor is there an existing comprehensive biogeography for the genus. Using ultraconserved element (UCE) phylogenomic data, we constructed the first molecular phylogeny for the genus Nomada and tested the monophyly of 16 morphologically established species groups. We also estimated divergence dates using fossil calibration points and inferred the origin and of this genus around the globe. Our phylogeny recovered 14 of the 16 previously established species groups as monophyletic. The superba and ruficornis groups, however, were recovered as non-monophyletic and need to be re-evaluated using morphology. Divergence dating and historic biogeographic analyses performed on the phylogenetic reconstruction indicates that Nomada most likely originated in the Holarctic ∼65 Mya. Geodispersal into the southern hemisphere occurred three times; once during the Eocene into the Afrotropics, once during the Oligocene into the Neotropics, and once during the Miocene into Australasia. Geodispersal across the Holarctic was most frequent and occurred repeatedly throughout the Cenozoic era, using the De Geer, Thulean, and the Bering Land Bridges. This is the first instance of a bee using both the Thulean and De Geer land bridges and has implications of how early bee species dispersed throughout the Palearctic in the late Cretaceous and early Paleogene.

Phylogenetic relationships and the evolution of host preferences in the largest clade of brood parasitic bees (Apidae: Nomadinae)

Brood parasites (also known as cleptoparasites) represent a substantial fraction of global bee diversity. Rather than constructing their own nests, these species instead invade those of host bees to lay their eggs. Larvae then hatch and consume the food provisions intended for the host's offspring. While this life history strategy has evolved numerous times across the phylogeny of bees, the oldest and most speciose parasitic clade is the subfamily Nomadinae (Apidae). However, the phylogenetic relationships among brood parasitic apids both within and outside the Nomadinae have not been fully resolved. Here, we present new findings on the phylogeny of this diverse group of brood parasites based on ultraconserved element (UCE) sequence data and extensive taxon sampling with 114 nomadine species representing all tribes. We suggest a broader definition of the subfamily Nomadinae to describe a clade that includes almost all parasitic members of the family Apidae. The tribe Melectini forms the sister group to all other Nomadinae, while the remainder of the subfamily is composed of two sister clades: a "nomadine line" representing the former Nomadinae sensu stricto, and an "ericrocidine line" that unites several mostly Neotropical lineages. We find the tribe Osirini Handlirsch to be polyphyletic, and divide it into three lineages, including the newly described Parepeolini trib. nov. In addition to our taxonomic findings, we use our phylogeny to explore the evolution of different modes of parasitism, detecting two independent transitions from closed-cell to open-cell parasitism. Finally, we examine how nomadine host-parasite associations have evolved over time. In support of Emery's rule, which suggests close relationships between hosts and parasites, we confirm that the earliest nomadines were parasites of their close free-living relatives within the family Apidae, but that over time their host range broadened to include more distantly related hosts spanning the diversity of bees. This expanded breadth of host taxa may also be associated with the transition to open-cell parasitism.

Museomics and phylogenomics with protein-encoding ultraconserved elements illuminate the evolution of life history and phallic morphology of flesh flies (Diptera: Sarcophagidae)

BACKGROUND

The common name of the Flesh flies (Sarcophagidae) usually relates them with organisms feeding on decomposing organic matter, although the biology of one of the largest radiations among insects also includes predation, coprophagy, and even kleptoparasitism. The question of whether the ancestor of all sarcophagids was a predator or a decomposer, or in association to which host have sarcophagids evolved, has thus always piqued the curiosity of flesh fly specialists. Such curiosity has often been hindered by both the impossibility of having a well-supported phylogeny of Sarcophagidae and its sister group to trace live habits and the scarcity of information on the biology of the group. Using a phylogenomic dataset of protein-encoding ultraconserved elements from representatives of all three subfamilies of Sarcophagidae as ingroup and a large Calyptratae outgroup, a robust phylogenetic framework and timescale are generated to understand flesh fly systematics and the evolution of their life histories.

RESULTS

The evolutionary history for Sarcophagidae reconstructed here differs considerably from previous hypotheses. Within subfamily Sarcophaginae, a group of predatory flies, including genera Lepidodexia and Boettcheria, emerged as sister-group to the rest of Sarcophaginae. The genera Oxysarcodexia, Ravinia, and Tricharaea, long considered archaic and early-branching coprophagous and sarcosaprophagous lineages, were found nested well within the Sarcophaginae as sister-group to the sarcosaprophagous Microcerella. Predation on invertebrates is suggested as the ancestral and dominant strategy throughout the early evolution of flesh flies. Several transitions from predation to sarcosaprophagy and coprophagy occur across the sarcophagid phylogenetic tree, in contrast with almost no transitions from sarcosaprophagy or coprophagy to predatory habits. Regarding the morphological evolution of flesh flies, there might be a concerted evolution of male genitalia traits, such as the phallotrema position and the juxta, or the vesica and the folding of the phallotrema. One diversification rate shift was inferred in the evolution of sarcophagids, which is related to the origin of genus Sarcophaga.

CONCLUSIONS

This study has a significant impact on understanding sarcophagid evolution and highlights the importance of having a robust phylogenetic framework to reconstruct the ancestral character state of biological and morphological characters. I discuss the evolution of life histories of the family in relation to their hosts or substrates and outline how sarcosaprophagy, coprophagy, and kleptoparasitism behavior on various hosts may have evolved from predation on invertebrates. This study provides a phylogenetic framework for further physiological and comparative genomic work between predatory, sarcosaprophagous, coprophagous, and kleptoparasitic lineages, which could also have significant implications for the evolution of diverse life histories in other Diptera.

Comprehensive phylogenomic analyses re-write the evolution of parasitism within cynipoid wasps

BACKGROUND

Parasitoidism, a specialized life strategy in which a parasite eventually kills its host, is frequently found within the insect order Hymenoptera (wasps, ants and bees). A parasitoid lifestyle is one of two dominant life strategies within the hymenopteran superfamily Cynipoidea, with the other being an unusual plant-feeding behavior known as galling. Less commonly, cynipoid wasps exhibit inquilinism, a strategy where some species have adapted to usurp other species' galls instead of inducing their own. Using a phylogenomic data set of ultraconserved elements from nearly all lineages of Cynipoidea, we here generate a robust phylogenetic framework and timescale to understand cynipoid systematics and the evolution of these life histories.

RESULTS

Our reconstructed evolutionary history for Cynipoidea differs considerably from previous hypotheses. Rooting our analyses with non-cynipoid outgroups, the Paraulacini, a group of inquilines, emerged as sister-group to the rest of Cynipoidea, rendering the gall wasp family Cynipidae paraphyletic. The families Ibaliidae and Liopteridae, long considered archaic and early-branching parasitoid lineages, were found nested well within the Cynipoidea as sister-group to the parasitoid Figitidae. Cynipoidea originated in the early Jurassic around 190 Ma. Either inquilinism or parasitoidism is suggested as the ancestral and dominant strategy throughout the early evolution of cynipoids, depending on whether a simple (three states: parasitoidism, inquilinism and galling) or more complex (seven states: parasitoidism, inquilinism and galling split by host use) model is employed.

CONCLUSIONS

Our study has significant impact on understanding cynipoid evolution and highlights the importance of adequate outgroup sampling. We discuss the evolutionary timescale of the superfamily in relation to their insect hosts and host plants, and outline how phytophagous galling behavior may have evolved from entomophagous, parasitoid cynipoids. Our study has established the framework for further physiological and comparative genomic work between gall-making, inquiline and parasitoid lineages, which could also have significant implications for the evolution of diverse life histories in other Hymenoptera.

UCE data reveal multiple origins of rose gallers in North America: Global phylogeny of Diplolepis Geoffroy (Hymenoptera: Cynipidae)

Gall wasps in the genus Diplolepis Geoffroy are specialized herbivores that induce galls exclusively on roses. Despite their wide distribution across the Holarctic, little is known about their evolutionary history. Here we present the first phylogenomic tree of global Diplolepis reconstructed using Ultraconserved Elements (UCEs), resulting in a robust phylogeny based on 757 genes. Results support the existence of two principal clades: a Nearctic stem-galler clade, and a Holarctic leaf-galler clade that further splits into two Palearctic groups and one Nearctic group. This topology is congruent with a previous study based on the mitochondrial gene COI, an unexpected result given the common occurrence of mitonuclear discordance in closely related oak gall wasp lineages. Most Diplolepis species were recovered as reciprocally monophyletic, with some notable exceptions such as the D. polita and the D. ignota complex, for which species boundaries remain unresolved. Historical biogeographic reconstruction was unable to pinpoint the origin of Diplolepis, but confirms two independent incursions into the Nearctic. Ancestral state reconstruction analysis highlights the conservatism of gall location on the host plants, as shifts to different host organs are relatively rare. We suggest that Diplolepis were originally leaf gallers, with a Nearctic stem-galler clade undergoing a major plant organ switch onto rose stems. Host organ switch or reversal is uncommon, which suggests a level of conservatism. Our study showcases the resolving power of UCEs at the species level while also suggesting improvements to advance future Cynipoidea phylogenomics. Our results also highlight the additional sampling needed to clarify taxonomic relationships in the Nearctic and eastern Palearctic regions.

Phylogenomic species delimitation in microendemic frogs of the Brazilian Atlantic Forest

The advent of next-generation sequencing allows researchers to use large-scale datasets for species delimitation analyses, yet one can envision an inflection point where the added accuracy of including more loci does not offset the increased computational burden. One alternative to including all loci could be to prioritize the analysis of loci for which there is an expectation of high informativeness. Here, we explore the issue of species delimitation and locus selection with montane species from two anuran genera that have been isolated in sky islands across the southern Brazilian Atlantic Forest: Melanophryniscus (Bufonidae) and Brachycephalus (Brachycephalidae). To delimit species, we obtained genetic data using target enrichment of ultraconserved elements from 32 populations (13 for Melanophryniscus and 19 for Brachycephalus), and we were able to create datasets that included over 800 loci with no missing data. We ranked loci according to their number of parsimony-informative sites, and we performed species delimitation analyses using BPP with the most informative 10, 20, 40, 80, 160, 320, and 640 loci. We identified three types of phylogenetic node: nodes with either consistently high or low support regardless of the number of loci or their informativeness and nodes that were initially poorly supported where support became stronger as we included more data. When viewed across all sensitivity analyses, our results suggest that the current species richness in both genera is likely underestimated. In addition, our results show the effects of different sampling strategies on species delimitation using phylogenomic datasets.

Multilocus molecular systematics of the circumtropical reef-fish genus Abudefduf (Pomacentridae): history, geography and ecology of speciation

We investigated a pantropical sub-family and genus of damselfishes, the sergeant-majors (Pomacentridae: Abudefdufinae: Abudefduf), to identify the tempo and mechanisms of speciation in the lineage. We examined sequence capture data from 500 loci and 20 species, with multiple individuals sampled from across the geographic ranges of widespread species. Utilizing a maximum likelihood framework, as well as a time-calibrated Bayesian phylogeny, the following key questions are addressed: What is the historical tempo of speciation? What are the relative contributions of vicariant, peripatric and parapatric speciation to sergeant-major diversity? How is speciation related to major variation in trophic ecology? The approximately 20 species of sergeant-majors fall into three main lineages. The ancestral condition appears to be benthivory, which is predominant in two lineages comprising six species. The remaining species of sergeant-majors, of which there are at least 15, fall within a clade composed entirely of planktivores. This clade is sister to a benthivore clade that included one species, Abudefduf notatus, in transition to planktivory. Most speciation of sergeant-majors, which appeared ∼24 million years ago, occurred in the last 10 million years. Present distributional patterns indicate vicariant speciation precipitated by the closure of land barriers between both sides of the Atlantic and the Pacific, and the emergence of land between the Indian and Pacific Oceans. Within this backdrop, frequent oscillations in sea level over the last 10 million years also appear to have generated conditions suitable for both peripatric and vicariant speciation, and most speciation within the genus appears linked to these changes in sea level. Diversification within the genus has been concentrated in planktivorous seargeant-majors rather than benthivores. The root cause is unclear, but does not appear to be related to differences in dispersal potential, which is greater in the planktivorous species, due to the ability of their post-larval juveniles to raft with floating debris. This elevated speciation rate in planktivores and their propensity to form local endemics may reflect relaxation of selective pressures (e.g., on crypticity) that limit speciation in benthivorous sergeant-majors. Finally, our data allow us to clarify relationships of geminate sergeant-major species, indicating that there are subdivisions within the Atlantic for both benthivore and planktivore geminate pairs that may have misled previous studies.

Explosive diversification of marine fishes at the Cretaceous-Palaeogene boundary

The Cretaceous-Palaeogene (K-Pg) mass extinction is linked to the rapid emergence of ecologically divergent higher taxa (for example, families and orders) across terrestrial vertebrates, but its impact on the diversification of marine vertebrates is less clear. Spiny-rayed fishes (Acanthomorpha) provide an ideal system for exploring the effects of the K-Pg on fish diversification, yet despite decades of morphological and molecular phylogenetic efforts, resolution of both early diverging lineages and enormously diverse subclades remains problematic. Recent multilocus studies have provided the first resolved phylogenetic backbone for acanthomorphs and suggested novel relationships among major lineages. However, these new relationships and associated timescales have not been interrogated using phylogenomic approaches. Here, we use targeted enrichment of >1,000 ultraconserved elements in conjunction with a divergence time analysis to resolve relationships among 120 major acanthomorph lineages and provide a new timescale for acanthomorph radiation. Our results include a well-supported topology that strongly resolves relationships along the acanthomorph backbone and the recovery of several new relationships within six major percomorph subclades. Divergence time analyses also reveal that crown ages for five of these subclades, and for the bulk of the species diversity in the sixth, coincide with the K-Pg boundary, with divergences between anatomically and ecologically distinctive suprafamilial clades concentrated in the first 10 million years of the Cenozoic.

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

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

Capturing Darwin's dream

Evolutionary biologists from Darwin forward have dreamed of having data that would elucidate our understanding of evolutionary history and the diversity of life. Sequence capture is a relatively old DNA technology, but its use is growing rapidly due to advances in (i) massively parallel DNA sequencing approaches and instruments, (ii) massively parallel bait construction, (iii) methods to identify target regions and (iv) sample preparation. We give a little historical context to these developments, summarize some of the important advances reported in this special issue and point to further advances that can be made to help fulfill Darwin's dream.

Sequence Capture and Phylogenetic Utility of Genomic Ultraconserved Elements Obtained from Pinned Insect Specimens

Obtaining sequence data from historical museum specimens has been a growing research interest, invigorated by next-generation sequencing methods that allow inputs of highly degraded DNA. We applied a target enrichment and next-generation sequencing protocol to generate ultraconserved elements (UCEs) from 51 large carpenter bee specimens (genus Xylocopa), representing 25 species with specimen ages ranging from 2-121 years. We measured the correlation between specimen age and DNA yield (pre- and post-library preparation DNA concentration) and several UCE sequence capture statistics (raw read count, UCE reads on target, UCE mean contig length and UCE locus count) with linear regression models. We performed piecewise regression to test for specific breakpoints in the relationship of specimen age and DNA yield and sequence capture variables. Additionally, we compared UCE data from newer and older specimens of the same species and reconstructed their phylogeny in order to confirm the validity of our data. We recovered 6-972 UCE loci from samples with pre-library DNA concentrations ranging from 0.06-9.8 ng/μL. All investigated DNA yield and sequence capture variables were significantly but only moderately negatively correlated with specimen age. Specimens of age 20 years or less had significantly higher pre- and post-library concentrations, UCE contig lengths, and locus counts compared to specimens older than 20 years. We found breakpoints in our data indicating a decrease of the initial detrimental effect of specimen age on pre- and post-library DNA concentration and UCE contig length starting around 21-39 years after preservation. Our phylogenetic results confirmed the integrity of our data, giving preliminary insights into relationships within Xylocopa. We consider the effect of additional factors not measured in this study on our age-related sequence capture results, such as DNA fragmentation and preservation method, and discuss the promise of the UCE approach for large-scale projects in insect phylogenomics using museum specimens.

Discovery of an expanded set of avian leukosis subgroup E proviruses in chickens using Vermillion, a novel sequence capture and analysis pipeline [corrected]

Transposable elements (TEs), such as endogenous retroviruses (ERVs), are common in the genomes of vertebrates. ERVs result from retroviral infections of germ-line cells, and once integrated into host DNA they become part of the host's heritable genetic material. ERVs have been ascribed positive effects on host physiology such as the generation of novel, adaptive genetic variation and resistance to infection, as well as negative effects as agents of tumorigenesis and disease. The avian leukosis virus subgroup E family (ALVE) of endogenous viruses of chickens has been used as a model system for studying the effects of ERVs on host physiology, and approximately 30 distinct ALVE proviruses have been described in the Gallus gallus genome. In this report we describe the development of a software tool, which we call Vermillion, and the use of this tool in combination with targeted next-generation sequencing (NGS) to increase the number of known proviruses belonging to the ALVE family of ERVs in the chicken genome by 4-fold, including expanding the number of known ALVE elements on chromosome 1 (Gga1) from the current 9 to a total of 40. Although we focused on the discovery of ALVE elements in chickens, with appropriate selection of target sequences Vermillion can be used to develop profiles of other families of ERVs and TEs in chickens as well as in species other than the chicken.

BaitFisher: A Software Package for Multispecies Target DNA Enrichment Probe Design

Target DNA enrichment combined with high-throughput sequencing technologies is a powerful approach to probing a large number of loci in genomes of interest. However, software algorithms that explicitly consider nucleotide sequence information of target loci in multiple reference species for optimizing design of target enrichment baits to be applicable across a wide range of species have not been developed. Here we present an algorithm that infers target DNA enrichment baits from multiple nucleotide sequence alignments. By applying clustering methods and the combinatorial 1-center sequence optimization to bait design, we are able to minimize the total number of baits required to efficiently probe target loci in multiple species. Consequently, more loci can be probed across species with a given number of baits. Using transcript sequences of 24 apoid wasps (Hymenoptera: Crabronidae, Sphecidae) from the 1KITE project and the gene models of Nasonia vitripennis, we inferred 57,650, 120-bp-long baits for capturing 378 coding sequence sections of 282 genes in apoid wasps. Illumina reduced-representation library sequencing confirmed successful enrichment of the target DNA when applying these baits to DNA of various apoid wasps. The designed baits furthermore enriched a major fraction of the target DNA in distantly related Hymenoptera, such as Formicidae and Chalcidoidea, highlighting the baits' broad taxonomic applicability. The availability of baits with broad taxonomic applicability is of major interest in numerous disciplines, ranging from phylogenetics to biodiversity monitoring. We implemented our new approach in a software package, called BaitFisher, which is open source and freely available at https://github.com/cmayer/BaitFisher-package.git.

Phylogenomic methods outperform traditional multi-locus approaches in resolving deep evolutionary history: a case study of formicine ants

BACKGROUND

Ultraconserved elements (UCEs) have been successfully used in phylogenomics for a variety of taxa, but their power in phylogenetic inference has yet to be extensively compared with that of traditional Sanger sequencing data sets. Moreover, UCE data on invertebrates, including insects, are sparse. We compared the phylogenetic informativeness of 959 UCE loci with a multi-locus data set of ten nuclear markers obtained via Sanger sequencing, testing the ability of these two types of data to resolve and date the evolutionary history of the second most species-rich subfamily of ants in the world, the Formicinae.

RESULTS

Phylogenetic analyses show that UCEs are superior in resolving ancient and shallow relationships in formicine ants, demonstrated by increased node support and a more resolved phylogeny. Phylogenetic informativeness metrics indicate a twofold improvement relative to the 10-gene data matrix generated from the identical set of taxa. We were able to significantly improve formicine classification based on our comprehensive UCE phylogeny. Our divergence age estimations, using both UCE and Sanger data, indicate that crown-group Formicinae are older (104-117 Ma) than previously suggested. Biogeographic analyses infer that the diversification of the subfamily has occurred on all continents with no particular hub of cladogenesis.

CONCLUSIONS

We found UCEs to be far superior to the multi-locus data set in estimating formicine relationships. The early history of the clade remains uncertain due to ancient rapid divergence events that are unresolvable even with our genomic-scale data, although this might be largely an effect of several problematic taxa subtended by long branches. Our comparison of divergence ages from both Sanger and UCE data demonstrates the effectiveness of UCEs for dating analyses. This comparative study highlights both the promise and limitations of UCEs for insect phylogenomics, and will prove useful to the growing number of evolutionary biologists considering the transition from Sanger to next-generation sequencing approaches.

Tumor neoantigens: building a framework for personalized cancer immunotherapy

It is now well established that the immune system can recognize developing cancers and that therapeutic manipulation of immunity can induce tumor regression. The capacity to manifest remarkably durable responses in some patients has been ascribed in part to T cells that can (a) kill tumor cells directly, (b) orchestrate diverse antitumor immune responses, (c) manifest long-lasting memory, and (d) display remarkable specificity for tumor-derived proteins. This specificity stems from fundamental differences between cancer cells and their normal counterparts in that the former develop protein-altering mutations and undergo epigenetic and genetic alterations, resulting in aberrant protein expression. These events can result in formation of tumor antigens. The identification of mutated and aberrantly expressed self-tumor antigens has historically been time consuming and laborious. While mutant antigens are usually expressed in a tumor-specific manner, aberrantly expressed antigens are often shared between cancers and, therefore, in the past, have been the major focus of therapeutic cancer vaccines. However, advances in next-generation sequencing and epitope prediction now permit the rapid identification of mutant tumor neoantigens. This review focuses on a discussion of mutant tumor neoantigens and their use in personalizing cancer immunotherapies.

Target enrichment of ultraconserved elements from arthropods provides a genomic perspective on relationships among Hymenoptera

Gaining a genomic perspective on phylogeny requires the collection of data from many putatively independent loci across the genome. Among insects, an increasingly common approach to collecting this class of data involves transcriptome sequencing, because few insects have high-quality genome sequences available; assembling new genomes remains a limiting factor; the transcribed portion of the genome is a reasonable, reduced subset of the genome to target; and the data collected from transcribed portions of the genome are similar in composition to the types of data with which biologists have traditionally worked (e.g. exons). However, molecular techniques requiring RNA as a template, including transcriptome sequencing, are limited to using very high-quality source materials, which are often unavailable from a large proportion of biologically important insect samples. Recent research suggests that DNA-based target enrichment of conserved genomic elements offers another path to collecting phylogenomic data across insect taxa, provided that conserved elements are present in and can be collected from insect genomes. Here, we identify a large set (n = 1510) of ultraconserved elements (UCEs) shared among the insect order Hymenoptera. We used in silico analyses to show that these loci accurately reconstruct relationships among genome-enabled hymenoptera, and we designed a set of RNA baits (n = 2749) for enriching these loci that researchers can use with DNA templates extracted from a variety of sources. We used our UCE bait set to enrich an average of 721 UCE loci from 30 hymenopteran taxa, and we used these UCE loci to reconstruct phylogenetic relationships spanning very old (≥220 Ma) to very young (≤1 Ma) divergences among hymenopteran lineages. In contrast to a recent study addressing hymenopteran phylogeny using transcriptome data, we found ants to be sister to all remaining aculeate lineages with complete support, although this result could be explained by factors such as taxon sampling. We discuss this approach and our results in the context of elucidating the evolutionary history of one of the most diverse and speciose animal orders.

Genetic and clonal dissection of murine small cell lung carcinoma progression by genome sequencing

Small cell lung carcinoma (SCLC) is a highly lethal, smoking-associated cancer with few known targetable genetic alterations. Using genome sequencing, we characterized the somatic evolution of a genetically engineered mouse model (GEMM) of SCLC initiated by loss of Trp53 and Rb1. We identified alterations in DNA copy number and complex genomic rearrangements and demonstrated a low somatic point mutation frequency in the absence of tobacco mutagens. Alterations targeting the tumor suppressor Pten occurred in the majority of murine SCLC studied, and engineered Pten deletion accelerated murine SCLC and abrogated loss of Chr19 in Trp53; Rb1; Pten compound mutant tumors. Finally, we found evidence for polyclonal and sequential metastatic spread of murine SCLC by comparative sequencing of families of related primary tumors and metastases. We propose a temporal model of SCLC tumorigenesis with implications for human SCLC therapeutics and the nature of cancer-genome evolution in GEMMs.

Creation and application of immortalized bait libraries for targeted enrichment and next-generation sequencing

Since the introduction of next-generation sequencing, several techniques have been developed to selectively enrich and sequence specific parts of the genome at high coverage. These techniques include enzymatic methods employing molecular inversion probes, PCR based approaches, hybrid capture, and in-solution capture. In-solution capture employs RNA probes transcribed from a pool of DNA template oligos designed to match regions of interest to specifically bind and enrich genomic DNA fragments. This method is highly efficient, especially if genomic target regions are large in size or quantity. Diverse in-solution capture kits are available, but are costly when large sample numbers need to be analyzed. Here we present a cost-effective strategy for the design of custom DNA libraries, their transcription into RNA libraries, and application for in-solution capture. We show the efficacy by comparing the method to a commercial kit and further demonstrate that emulsion PCR can be used for bias free amplification and virtual immortalization of DNA template libraries.

Not all sequence tags are created equal: designing and validating sequence identification tags robust to indels

Ligating adapters with unique synthetic oligonucleotide sequences (sequence tags) onto individual DNA samples before massively parallel sequencing is a popular and efficient way to obtain sequence data from many individual samples. Tag sequences should be numerous and sufficiently different to ensure sequencing, replication, and oligonucleotide synthesis errors do not cause tags to be unrecoverable or confused. However, many design approaches only protect against substitution errors during sequencing and extant tag sets contain too few tag sequences. We developed an open-source software package to validate sequence tags for conformance to two distance metrics and design sequence tags robust to indel and substitution errors. We use this software package to evaluate several commercial and non-commercial sequence tag sets, design several large sets (max(count) = 7,198) of edit metric sequence tags having different lengths and degrees of error correction, and integrate a subset of these edit metric tags to polymerase chain reaction (PCR) primers and sequencing adapters. We validate a subset of these edit metric tagged PCR primers and sequencing adapters by sequencing on several platforms and subsequent comparison to commercially available alternatives. We find that several commonly used sets of sequence tags or design methodologies used to produce sequence tags do not meet the minimum expectations of their underlying distance metric, and we find that PCR primers and sequencing adapters incorporating edit metric sequence tags designed by our software package perform as well as their commercial counterparts. We suggest that researchers evaluate sequence tags prior to use or evaluate tags that they have been using. The sequence tag sets we design improve on extant sets because they are large, valid across the set, and robust to the suite of substitution, insertion, and deletion errors affecting massively parallel sequencing workflows on all currently used platforms.

Rare target enrichment for ultrasensitive PCR detection using cot-rehybridization and duplex-specific nuclease

Nucleic acid detection by polymerase chain reaction (PCR) is invaluable for the detection of dilute and rare sequences, including pathogens and infrequent species in complex clinical and environmental backgrounds. The presence of excess complex background nucleic acid can reduce sensitivity and specificity. This is because mispriming can cause failure of the amplification reaction. Here we describe a new approach to ultrasensitive PCR detection, using enrichment of rare target nucleic acid from abundant background by combining the classic technique of cot-rehybridization to convert the abundant background to double-stranded form, with the use of a newly described, highly processive duplex-specific crab nuclease. We show that trace sequences in a vast excess of background DNA can be undetectable by PCR, independent of the amount of the mixture added to the PCR, and that these sequences can be made detectable by background suppression using this method.

High-throughput detection of actionable genomic alterations in clinical tumor samples by targeted, massively parallel sequencing

Knowledge of "actionable" somatic genomic alterations present in each tumor (e.g., point mutations, small insertions/deletions, and copy-number alterations that direct therapeutic options) should facilitate individualized approaches to cancer treatment. However, clinical implementation of systematic genomic profiling has rarely been achieved beyond limited numbers of oncogene point mutations. To address this challenge, we utilized a targeted, massively parallel sequencing approach to detect tumor genomic alterations in formalin-fixed, paraffin-embedded (FFPE) tumor samples. Nearly 400-fold mean sequence coverage was achieved, and single-nucleotide sequence variants, small insertions/deletions, and chromosomal copynumber alterations were detected simultaneously with high accuracy compared with other methods in clinical use. Putatively actionable genomic alterations, including those that predict sensitivity or resistance to established and experimental therapies, were detected in each tumor sample tested. Thus, targeted deep sequencing of clinical tumor material may enable mutation-driven clinical trials and, ultimately, "personalized" cancer treatment.

SIGNIFICANCE

Despite the rapid proliferation of targeted therapeutic agents, systematic methods to profile clinically relevant tumor genomic alterations remain underdeveloped. We describe a sequencingbased approach to identifying genomic alterations in FFPE tumor samples. These studies affirm the feasibility and clinical utility of targeted sequencing in the oncology arena and provide a foundation for genomics-based stratification of cancer patients.