Whole-slide laser microdissection for tumour enrichment.

The practical application of genome-scale technologies to precision oncology research requires flexible tissue processing strategies that can be used to differentially select both tumour and normal cell populations from formalin-fixed, paraffin-embedded tissues. As tumour sequencing scales towards clinical implementation, practical difficulties in scheduling and obtaining fresh tissue biopsies at scale, including blood samples as surrogates for matched 'normal' DNA, have focused attention on the use of formalin-preserved clinical samples collected routinely for diagnostic purposes. In practice, such samples often contain both tumour and normal cells which, if correctly partitioned, could be used to profile both tumour and normal genomes, thus identifying somatic alterations. Here we report a semi-automated method for laser microdissecting entire slide-mounted tissue sections to enrich for cells of interest with sufficient yield for whole genome and transcriptome sequencing. Using this method, we demonstrated enrichment of tumour material from mixed tumour-normal samples by up to 67%. Leveraging new methods that allow for the extraction of high-quality nucleic acids from small amounts of formalin-fixed tissues, we further showed that the method was successful in yielding sequence data of sufficient quality for use in BC Cancer's Personalized OncoGenomics (POG) program. © 2020 The Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Sources of erroneous sequences and artifact chimeric reads in next generation sequencing of genomic DNA from formalin-fixed paraffin-embedded samples.

Tissues used in pathology laboratories are typically stored in the form of formalin-fixed, paraffin-embedded (FFPE) samples. One important consideration in repurposing FFPE material for next generation sequencing (NGS) analysis is the sequencing artifacts that can arise from the significant damage to nucleic acids due to treatment with formalin, storage at room temperature and extraction. One such class of artifacts consists of chimeric reads that appear to be derived from non-contiguous portions of the genome. Here, we show that a major proportion of such chimeric reads align to both the 'Watson' and 'Crick' strands of the reference genome. We refer to these as strand-split artifact reads (SSARs). This study provides a conceptual framework for the mechanistic basis of the genesis of SSARs and other chimeric artifacts along with supporting experimental evidence, which have led to approaches to reduce the levels of such artifacts. We demonstrate that one of these approaches, involving S1 nuclease-mediated removal of single-stranded fragments and overhangs, also reduces sequence bias, base error rates, and false positive detection of copy number and single nucleotide variants. Finally, we describe an analytical approach for quantifying SSARs from NGS data.

Increasing quality, throughput and speed of sample preparation for strand-specific messenger RNA sequencing.

BACKGROUND: RNA-Sequencing (RNA-seq) is now commonly used to reveal quantitative spatiotemporal snapshots of the transcriptome, the structures of transcripts (splice variants and fusions) and landscapes of expressed mutations. However, standard approaches for library construction typically require relatively high amounts of input RNA, are labor intensive, and are time consuming. METHODS: Here, we report the outcome of a systematic effort to optimize and streamline steps in strand-specific RNA-seq library construction. RESULTS: This work has resulted in the identification of an optimized messenger RNA isolation protocol, a potent reverse transcriptase for cDNA synthesis, and an efficient chemistry and a simplified formulation of library construction reagents. We also present an optimization of bead-based purification and size selection designed to maximize the recovery of cDNA fragments. CONCLUSIONS: These developments have allowed us to assemble a rapid high throughput pipeline that produces high quality data from amounts of total RNA as low as 25 ng. While the focus of this study is on RNA-seq sample preparation, some of these developments are also relevant to other next-generation sequencing library types.

Assembling the 20 Gb white spruce (Picea glauca) genome from whole-genome shotgun sequencing data.

UNLABELLED: White spruce (Picea glauca) is a dominant conifer of the boreal forests of North America, and providing genomics resources for this commercially valuable tree will help improve forest management and conservation efforts. Sequencing and assembling the large and highly repetitive spruce genome though pushes the boundaries of the current technology. Here, we describe a whole-genome shotgun sequencing strategy using two Illumina sequencing platforms and an assembly approach using the ABySS software. We report a 20.8 giga base pairs draft genome in 4.9 million scaffolds, with a scaffold N50 of 20,356 bp. We demonstrate how recent improvements in the sequencing technology, especially increasing read lengths and paired end reads from longer fragments have a major impact on the assembly contiguity. We also note that scalable bioinformatics tools are instrumental in providing rapid draft assemblies. AVAILABILITY: The Picea glauca genome sequencing and assembly data are available through NCBI (Accession#: ALWZ0100000000 PID: PRJNA83435). http://www.ncbi.nlm.nih.gov/bioproject/83435.

Herbivory and floral signaling: phenotypic plasticity and tradeoffs between reproduction and indirect defense.

Plant defense against herbivores may compromise attraction of mutualists, yet information remains limited about the mechanisms underlying such signaling tradeoffs. Here, we investigated the effects of foliar herbivory by two herbivore species on defense compounds, floral signaling, pollinator and parasitoid attraction, and seed production. Herbivory generally reduced the quantity of many floral volatile organic compounds VOCs) in Brassica rapa. By contrast, floral color, flower diameter, and plant height remained unaffected. The decreased amounts of floral volatiles led to reduced attractiveness of flowers to pollinators, but increased the attractiveness of herbivore-infested plants to parasitoids. Plants infested with the native butterfly Pieris brassicae produced more flowers during early flowering, effectively compensating for the lower olfactory attractiveness. Herbivory by the invasive Spodoptera littoralis increased the amounts of glucobrassicanapin, and led to delayed flowering. These plants tended to attract fewer pollinators and to produce fewer seeds. Our study indicates a tradeoff between pollinator attraction and indirect defense (parasitoid attraction), which can be mitigated by reduced floral VOC emission and production of more early flowers. We suggest that this compensatory mechanism is specific to plant-herbivore associations with a coevolutionary history.

Worldwide population genetic structure of the oriental fruit moth (Grapholita molesta), a globally invasive pest.

BACKGROUND: Invasive pest species have large impacts on agricultural crop yields, and understanding their population dynamics is important for ensuring food security. The oriental fruit moth Grapholita molesta is a cosmopolitan pest of stone and pome fruit species including peach and apple, and historical records indicate that it has invaded North and South America, Europe, Australia and Africa from its putative native range in Asia over the past century. RESULTS: We used 13 microsatellite loci, including nine newly developed markers, to characterize global population structure of G. molesta. Approximately 15 individuals from each of 26 globally distributed populations were genotyped. A weak but significant global pattern of isolation-by-distance was found, and G. molesta populations were geographically structured on a continental scale. Evidence does not support that G. molesta was introduced to North America from Japan as previously proposed. However, G. molesta was probably introduced from North America to The Azores, South Africa, and Brazil, and from East Asia to Australia. Shared ancestry was inferred between populations from Western Europe and from Brazil, although it remains unresolved whether an introduction occurred from Europe to Brazil, or vice versa. Both genetic diversity and levels of inbreeding were surprisingly high across the range of G. molesta and were not higher or lower overall in introduced areas compared to native areas. There is little evidence for multiple introductions to each continent (except in the case of South America), or for admixture between populations from different origins. CONCLUSIONS: Cross-continental introductions of G. molesta appear to be infrequent, which is surprising given its rapid worldwide expansion over the past century. We suggest that area-wide spread via transport of fruits and other plant materials is a major mechanism of ongoing invasion, and management efforts should therefore target local and regional farming communities and distribution networks.

The genome sequence of the Loggerhead sea turtle, Caretta caretta Linnaeus 1758.

We present a genome assembly of Caretta caretta (the Loggerhead sea turtle; Chordata, Testudines, Cheloniidae), generated from genomic data from two unrelated females. The genome sequence is 2.13 gigabases in size. The assembly has a busco completion score of 96.1% and N50 of 130.95 Mb. The majority of the assembly is scaffolded into 28 chromosomal representations with a remaining 2% of the assembly being excluded from these.

Assembly and annotation of the black spruce genome provide insights on spruce phylogeny and evolution of stress response.

Black spruce (Picea mariana [Mill.] B.S.P.) is a dominant conifer species in the North American boreal forest that plays important ecological and economic roles. Here, we present the first genome assembly of P. mariana with a reconstructed genome size of 18.3 Gbp and NG50 scaffold length of 36.0 kbp. A total of 66,332 protein-coding sequences were predicted in silico and annotated based on sequence homology. We analyzed the evolutionary relationships between P. mariana and 5 other spruces for which complete nuclear and organelle genome sequences were available. The phylogenetic tree estimated from mitochondrial genome sequences agrees with biogeography; specifically, P. mariana was strongly supported as a sister lineage to P. glauca and 3 other taxa found in western North America, followed by the European Picea abies. We obtained mixed topologies with weaker statistical support in phylogenetic trees estimated from nuclear and chloroplast genome sequences, indicative of ancient reticulate evolution affecting these 2 genomes. Clustering of protein-coding sequences from the 6 Picea taxa and 2 Pinus species resulted in 34,776 orthogroups, 560 of which appeared to be specific to P. mariana. Analysis of these specific orthogroups and dN/dS analysis of positive selection signatures for 497 single-copy orthogroups identified gene functions mostly related to plant development and stress response. The P. mariana genome assembly and annotation provides a valuable resource for forest genetics research and applications in this broadly distributed species, especially in relation to climate adaptation.

Can plant resistance to specialist herbivores be explained by plant chemistry or resource use strategy?

At both a macro- and micro-evolutionary level, selection of and performance on host plants by specialist herbivores are thought to be governed partially by host plant chemistry. Thus far, there is little evidence to suggest that specialists can detect small structural differences in secondary metabolites of their hosts, or that such differences affect host choice or performance of specialists. We tested whether phytochemical differences between closely related plant species are correlated with specialist host choice. We conducted no-choice feeding trials using 17 plant species of three genera of tribe Senecioneae (Jacobaea, Packera, and Senecio; Asteraceae) and a more distantly related species (Cynoglossum officinale; Boraginaceae) containing pyrrolizidine alkaloids (PAs), and four PA-sequestering specialist herbivores of the genus Longitarsus (Chrysomelidae). We also assessed whether variation in feeding by specialist herbivores is attributable to different resource use strategies of the tested plant species. Plant resource use strategy was quantified by measuring leaf dry matter content, which is related to both plant nutritive value and to plant investment in quantitative defences. We found no evidence that intra-generic differences in PA profiles affect feeding by specialist herbivores. Instead, our results indicate that decisions to begin feeding are related to plant resource use strategy, while decisions to continue feeding are not based on any plant characteristics measured in this study. These findings imply that PA composition does not significantly affect host choice by these specialist herbivores. Leaf dry matter content is somewhat phylogenetically conserved, indicating that plants may have difficulty altering resource use strategy in response to selection pressure by herbivores and other environmental factors on an evolutionary time scale.

Pyrrolizidine alkaloid variation in shoots and roots of segregating hybrids between Jacobaea vulgaris and Jacobaea aquatica.

Hybridization can lead to novel qualitative or quantitative variation of secondary metabolite (SM) expression that can have ecological and evolutionary consequences. We measured pyrrolizidine alkaloid (PA) expression in the shoots and roots of a family including one Jacobaea vulgaris genotype and one Jacobaea aquatica genotype (parental genotypes), two F(1) hybrid genotypes, and 102 F(2) hybrid genotypes using liquid chromatography-tandem mass spectrometry (LC-MS/MS). We detected 37 PAs in the roots and shoots of J. vulgaris, J. aquatica and the hybrids. PA concentrations and compositions differed between genotypes, and between roots and shoots. Three otosenine-like PAs that only occurred in the shoots of parental genotypes were present in the roots of F(2) hybrids; PA compositions were sometimes novel in F(2) hybrids compared with parental genotypes, and in some cases transgressive PA expression occurred. We also found that PAs from within structural groups covaried both in the roots and in the shoots, and that PA expression was correlated between shoots and roots. Considerable and novel variation present among F(2) hybrids indicates that hybridization has a potential role in the evolution of PA diversity in the genus Jacobaea, and this hybrid system is useful for studying the genetic control of PA expression.

Long-distance dispersal and high genetic diversity are implicated in the invasive spread of the common reed, Phragmites australis (Poaceae), in northeastern North America.

PREMISE OF THE STUDY: The Eurasian subspecies of the common reed (Phragmites australis subsp. australis, hereafter abbreviated as P. a. australis) was introduced to North America in the late 18(th) century and rapidly expanded its range, posing an ecological threat to wetlands. In this study, we aimed to determine whether admixture among multiple lineages, dispersal mechanisms, and high genetic diversity have contributed to the invasion of P. a. australis in the northeastern part of its range. Understanding mechanisms of the P. a. australis invasion will 1) contribute to a broader understanding of the factors that facilitate plant invasion, and 2) help us to develop effective management strategies for wetlands threatened by P. a. australis invasion. METHODS: We used a population genetics approach incorporating nine microsatellite loci to study genetic diversity and population structure in relation to biogeography of introduced North American Phragmites a. australis stands in the northeastern continental region. KEY RESULTS: Phragmites a. australis is genetically diverse in the region studied here. Significant population structure exists, and population structure is likely influenced by both long-distance dispersal via major waterways, and short-distance dispersal overland. Different lineages sometimes colonize geographically proximate locations leading to opportunities for admixture. Clonal reproduction likely exaggerates geographical structure among some stands, although high genetic and clonal diversity within some stands implies that sexual reproduction occurs frequently in P. a. australis. CONCLUSIONS: A variety of factors, including admixture among multiple lineages, multiple modes of dispersal, and plasticity in reproductive strategy promote the invasion success of Phragmites a. australis. Wetland managers in the St. Lawrence River/Great Lakes region should focus monitoring efforts on the shores of conservation lands to prevent the establishment of propagules from novel lineages.

The relationship between structurally different pyrrolizidine alkaloids and western flower thrips resistance in F(2) hybrids of Jacobaea vulgaris and Jacobaea aquatica.

Segregating plant hybrids often have more ecological and molecular variability compared to parental species, and are therefore useful for studying relationships between different traits, and the adaptive significance of trait variation. Hybrid systems have been used to study the relationship between the expression of plant defense compounds and herbivore susceptibility. We conducted a western flower thrips (WFT) bioassay using a hybrid family and investigated the relationship between WFT resistance and pyrrolizidine alkaloid (PA) variation. The hybrid family consisted of two parental (Jacobaea vulgaris and Jacobaea aquatica) genotypes, two F(1) genotypes, and 94 F(2) hybrid lines. The J. aquatica genotype was more susceptible to thrips attack than the J. vulgaris genotype, the two F(1) hybrids were as susceptible as J. aquatica, and susceptibility to WFT differed among F(2) hybrid lines: 69 F(2) lines were equally susceptible compared to J. aquatica, 10 F(2) lines were more susceptible than J. aquatica and 15 F(2) lines were as resistant as J. vulgaris or were intermediate to the two parental genotypes. Among 37 individual PAs that were derived from four structural groups (senecionine-, jacobine-, erucifoline- and otosenine-like PAs), the N-oxides of jacobine, jaconine, and jacoline were negatively correlated with feeding damage caused by WFT, and the tertiary amines of jacobine, jaconine, jacoline, and other PAs did not relate to feeding damage. Total PA concentration was negatively correlated with feeding damage. Among the four PA groups, only the total concentration of the jacobine-like PAs was negatively correlated with feeding damage. Multiple regression tests suggested that jacobine-like PAs play a greater role in WFT resistance than PAs from other structural groups. We found no evidence for synergistic effects of different PAs on WFT resistance. The relationship between PA variation and WFT feeding damage in the Jacobaea hybrids suggests a role for PAs in resistance to generalist insects.

Applications and implications of neutral versus non-neutral markers in molecular ecology.

The field of molecular ecology has expanded enormously in the past two decades, largely because of the growing ease with which neutral molecular genetic data can be obtained from virtually any taxonomic group. However, there is also a growing awareness that neutral molecular data can provide only partial insight into parameters such as genetic diversity, local adaptation, evolutionary potential, effective population size, and taxonomic designations. Here we review some of the applications of neutral versus adaptive markers in molecular ecology, discuss some of the advantages that can be obtained by supplementing studies of molecular ecology with data from non-neutral molecular markers, and summarize new methods that are enabling researchers to generate data from genes that are under selection.

Species by environment interactions affect pyrrolizidine alkaloid expression in Senecio jacobaea, Senecio aquaticus, and their hybrids.

We examined the effects of water and nutrient availability on the expression of the defense pyrrolizidine alkaloids (PAs) in Senecio jacobaea and S. aquaticus. Senecio jacobaea, and S. aquaticus are adapted to different natural habitats, characterized by differing abiotic conditions and different selection pressures from natural enemies. We tested if PA concentration and diversity are plastic over a range of water and nutrient treatments, and also whether such plasticity is dependent on plant species. We also tested the hypothesis that hybridization may contribute to PA diversity within plants, by comparing PA expression in parental species to that in artificially generated F(1) hybrids, and also in later generation natural hybrids between S. jacobaea and S. aquaticus. We showed that total PA concentration in roots and shoots is not dependent on species, but that species determines the pattern of PA diversification. Pyrrolizidine alkaloid diversity and concentration are both dependent on environmental factors. Hybrids produce a putatively novel PA, and this PA is conserved in natural hybrids, that are backcrossed to S. jacobaea. Natural hybrids that are backcrossed several times to S. jacobaea are with regard to PA diversity significantly different from S. jacobaea but not from S. aquaticus, while F(1) hybrids are in all cases more similar to S. jacobaea. These results collectively suggest that PA diversity is under the influence of natural selection.

Complete Mitochondrial Genome of a Gymnosperm, Sitka Spruce (Picea sitchensis), Indicates a Complex Physical Structure.

Plant mitochondrial genomes vary widely in size. Although many plant mitochondrial genomes have been sequenced and assembled, the vast majority are of angiosperms, and few are of gymnosperms. Most plant mitochondrial genomes are smaller than a megabase, with a few notable exceptions. We have sequenced and assembled the complete 5.5-Mb mitochondrial genome of Sitka spruce (Picea sitchensis), to date, one of the largest mitochondrial genomes of a gymnosperm. We sequenced the whole genome using Oxford Nanopore MinION, and then identified contigs of mitochondrial origin assembled from these long reads based on sequence homology to the white spruce mitochondrial genome. The assembly graph shows a multipartite genome structure, composed of one smaller 168-kb circular segment of DNA, and a larger 5.4-Mb single component with a branching structure. The assembly graph gives insight into a putative complex physical genome structure, and its branching points may represent active sites of recombination.

Complete Chloroplast Genome Sequence of a Black Spruce (Picea mariana) from Eastern Canada.

Here, we present the chloroplast genome sequence of black spruce (Picea mariana), a conifer widely distributed throughout North American boreal forests. This complete and annotated chloroplast sequence is 123,961 bp long and will contribute to future studies on the genetic basis of evolutionary change in spruce and adaptation in conifers.

A novel operative technique to manage a symptomatic synovial cyst associated with an os odontoideum.

Synovial cysts in association with an os odontoideum are rare. All previous cases have been treated via a posterior approach with good results. The authors present a patient with a synovial cyst in association with an os odontoideum and discuss a planned two-stage procedure to excise the cyst via a neuronavigation-assisted transoral-transpharyngeal approach, followed by a posterior C1/2 fixation.

De Novo Sequencing, Assembly, and Annotation of Four Threespine Stickleback Genomes Based on Microfluidic Partitioned DNA Libraries.

: The threespine stickleback is a geographically widespread and ecologically highly diverse fish that has emerged as a powerful model system for evolutionary genomics and developmental biology. Investigations in this species currently rely on a single high-quality reference genome, but would benefit from the availability of additional, independently sequenced and assembled genomes. We present here the assembly of four new stickleback genomes, based on the sequencing of microfluidic partitioned DNA libraries. The base pair lengths of the four genomes reach 92-101% of the standard reference genome length. Together with their de novo gene annotation, these assemblies offer a resource enhancing genomic investigations in stickleback. The genomes and their annotations are available from the Dryad Digital Repository (https://doi.org/10.5061/dryad.113j3h7).

Evaluation of protocols for rRNA depletion-based RNA sequencing of nanogram inputs of mammalian total RNA.

Next generation RNA-sequencing (RNA-seq) is a flexible approach that can be applied to a range of applications including global quantification of transcript expression, the characterization of RNA structure such as splicing patterns and profiling of expressed mutations. Many RNA-seq protocols require up to microgram levels of total RNA input amounts to generate high quality data, and thus remain impractical for the limited starting material amounts typically obtained from rare cell populations, such as those from early developmental stages or from laser micro-dissected clinical samples. Here, we present an assessment of the contemporary ribosomal RNA depletion-based protocols, and identify those that are suitable for inputs as low as 1-10 ng of intact total RNA and 100-500 ng of partially degraded RNA from formalin-fixed paraffin-embedded tissues.

A high-throughput protocol for isolating cell-free circulating tumor DNA from peripheral blood.

The analysis of cell-free circulating tumor DNA (ctDNA) is potentially a less invasive, more dynamic assessment of cancer progression and treatment response than characterizing solid tumor biopsies. Standard isolation methods require separation of plasma by centrifugation, a time-consuming step that complicates automation. To address these limitations, we present an automatable magnetic bead-based ctDNA isolation method that eliminates centrifugation to purify ctDNA directly from peripheral blood (PB). To develop and test our method, ctDNA from cancer patients was purified from PB and plasma. We found that allelic fractions of somatic single-nucleotide variants from target gene capture libraries were comparable, indicating that the PB ctDNA purification method may be a suitable replacement for the plasma-based protocols currently in use.