Large scale genome skimming from herbarium material for accurate plant identification and phylogenomics

Genome-wide RAD sequencing resolves the evolutionary history of serrate leaf Juniperus and reveals discordance with chloroplast phylogeny

Juniper (Juniperus) is an ecologically important conifer genus of the Northern Hemisphere, the members of which are often foundational tree species of arid regions. The serrate leaf margin clade is native to topologically variable regions in North America, where hybridization has likely played a prominent role in their diversification. Here we use a reduced-representation sequencing approach (ddRADseq) to generate a phylogenomic data set for 68 accessions representing all 22 species in the serrate leaf margin clade, as well as a number of close and distant relatives, to improve understanding of diversification in this group. Phylogenetic analyses using three methods (SVDquartets, maximum likelihood, and Bayesian) yielded highly congruent and well-resolved topologies. These phylogenies provided improved resolution relative to past analyses based on Sanger sequencing of nuclear and chloroplast DNA, and were largely consistent with taxonomic expectations based on geography and morphology. Calibration of a Bayesian phylogeny with fossil evidence produced divergence time estimates for the clade consistent with a late Oligocene origin in North America, followed by a period of elevated diversification between 12 and 5 Mya. Comparison of the ddRADseq phylogenies with a phylogeny based on Sanger-sequenced chloroplast DNA revealed five instances of pronounced discordance, illustrating the potential for chloroplast introgression, chloroplast transfer, or incomplete lineage sorting to influence organellar phylogeny. Our results improve understanding of the pattern and tempo of diversification in Juniperus, and highlight the utility of reduced-representation sequencing for resolving phylogenetic relationships in non-model organisms with reticulation and recent divergence.

Molecular identification and geographic origin of a post-Medieval elephant finding from southwestern Portugal using high-throughput sequencing

Molecular species identification plays a crucial role in archaeology and palaeontology, especially when diagnostic morphological characters are unavailable. Molecular markers have been used in forensic science to trace the geographic origin of wildlife products, such as ivory. So far, only a few studies have applied genetic methods to both identify the species and circumscribe the provenance of historic wildlife trade material. Here, by combining ancient DNA methods and genome skimming on a historical elephantid tooth found in southwestern Portugal, we aimed to identify its species, infer its placement in the elephantid phylogenetic tree, and triangulate its geographic origin. According to our results the specimen dates back to the eighteenth century CE and belongs to a female African forest elephant (non-hybrid Loxodonta cyclotis individual) geographically originated from west—west-central Africa, from areas where one of the four major mitochondrial clades of L. cyclotis is distributed. Historical evidence supports our inference, pointing out that the tooth should be considered as post-Medieval raw ivory trade material between West Africa and Portugal. Our study provides a comprehensive approach to study historical products and artefacts using archaeogenetics and contributes towards enlightening cultural and biological historical aspects of ivory trade in western Europe.

Genome skimming and NMR chemical fingerprinting provide quality assurance biotechnology to validate Sarsaparilla identity and purity

Sarsaparilla is a popular natural health product (NHP) that has been reported to be one of the most adulterated botanicals in the marketplace. Several plausible explanations are documented including economically motivated product substitution, unintentional errors due to ambiguous trade name associated with several different taxa, and wild harvesting of incorrect non-commercial plants. Unfortunately, this includes the case of an adulterant species Decalepis hamiltonii, a Red listed medicinal plant species by the International Union for Conservation of Nature (IUCN) and declared as a species with high conservation concern by the National Biodiversity Authority of India (NBA). This study provides validated genomic (genome skimming & DNA probes) and metabolomic (NMR chemical fingerprints) biotechnology solutions to prevent adulteration on both raw materials and finished products. This is also the first use of Oxford Nanopore on herbal products enabling the use of genome skimming as a tool for quality assurance within the supply chain of botanical ingredients. The validation of both genomics and metabolomics approach provided quality assurance perspective for both product identity and purity. This research enables manufactures and retailers to verify their supply chain is authentic and that consumers can enjoy safe, healthy products.

Mitochondrial Genome of Fagus sylvatica L. as a Source for Taxonomic Marker Development in the Fagales

European beech, Fagus sylvatica L., is one of the most important and widespread deciduous tree species in Central Europe and is widely managed for its hard wood. The complete DNA sequence of the mitochondrial genome of Fagus sylvatica L. was assembled and annotated based on Illumina MiSeq reads and validated using long reads from nanopore MinION sequencing. The genome assembled into a single DNA sequence of 504,715 bp in length containing 58 genes with predicted function, including 35 protein-coding, 20 tRNA and three rRNA genes. Additionally, 23 putative protein-coding genes were predicted supported by RNA-Seq data. Aiming at the development of taxon-specific mitochondrial genetic markers, the tool SNPtax was developed and applied to select genic SNPs potentially specific for different taxa within the Fagales. Further validation of a small SNP set resulted in the development of four CAPS markers specific for Fagus, Fagaceae, or Fagales, respectively, when considering over 100 individuals from a total of 69 species of deciduous trees and conifers from up to 15 families included in the marker validation. The CAPS marker set is suitable to identify the genus Fagus in DNA samples from tree tissues or wood products, including wood composite products.

Integrated Puncture Score: force-displacement weighted rind penetration tests improve stalk lodging resistance estimations in maize

BACKGROUND

Stalk lodging (breaking of agricultural plant stalks prior to harvest) is a multi-billion dollar a year problem. Rind penetration resistance tests have been used by plant scientists and breeders to estimate the stalk lodging resistance of maize for nearly a hundred years. However, the rind puncture method has two key limitations: (1) the predictive power of the test decreases significantly when measuring elite or pre-commercial hybrids, and (2) using rind penetration measurements as a breeding metric does not necessarily create stronger stalks. In this study, we present a new rind penetration method called the Integrated Puncture Score, which uses a modified rind penetration testing protocol and a physics-based model to provide a robust measure of stalk lodging resistance.

RESULTS

Two datasets, one with a diverse array of maize hybrids and one with only elite hybrids, were evaluated by comparing traditional rind penetration testing and the Integrated Puncture Score method to measurements of stalk bending strength. When evaluating the diverse set of hybrids, both methods were good predictors of stalk bending strength (R2 values of 0.67). However, when evaluating elite hybrids, the Integrated Puncture Score had an R2 value of 0.74 whereas the traditional method had an R2 value of 0.48. Additionally, the Integrated Puncture Score was able to differentiate between the strongest and weakest hybrids in the elite hybrid data set whereas the traditional rind penetration method was not. Additional experiments revealed strong evidence in favor of the data aggregation steps utilized to compute the Integrated Puncture Score.

CONCLUSIONS

This study presents a new method for evaluating rind penetration resistance that highly correlates with stalk bending strength and can possibly be used as a breeding index for assessing stalk lodging resistance. This research lays the foundation required to develop a field-based high-throughput phenotyping device for stalk lodging resistance.

Beyond DNA barcoding: The unrealized potential of genome skim data in sample identification

Genetic tools are increasingly used to identify and discriminate between species. One key transition in this process was the recognition of the potential of the ca 658bp fragment of the organelle cytochrome c oxidase I (COI) as a barcode region, which revolutionized animal bioidentification and lead, among others, to the instigation of the Barcode of Life Database (BOLD), containing currently barcodes from >7.9 million specimens. Following this discovery, suggestions for other organellar regions and markers, and the primers with which to amplify them, have been continuously proposed. Most recently, the field has taken the leap from PCR-based generation of DNA references into shotgun sequencing-based "genome skimming" alternatives, with the ultimate goal of assembling organellar reference genomes. Unfortunately, in genome skimming approaches, much of the nuclear genome (as much as 99% of the sequence data) is discarded, which is not only wasteful, but can also limit the power of discrimination at, or below, the species level. Here, we advocate that the full shotgun sequence data can be used to assign an identity (that we term for convenience its "DNA-mark") for both voucher and query samples, without requiring any computationally intensive pretreatment (e.g. assembly) of reads. We argue that if reference databases are populated with such "DNA-marks," it will enable future DNA-based taxonomic identification to complement, or even replace PCR of barcodes with genome skimming, and we discuss how such methodology ultimately could enable identification to population, or even individual, level.

The Treasure Vault Can be Opened: Large-Scale Genome Skimming Works Well Using Herbarium and Silica Gel Dried Material

Genome skimming has the potential for generating large data sets for DNA barcoding and wider biodiversity genomic studies, particularly via the assembly and annotation of full chloroplast (cpDNA) and nuclear ribosomal DNA (nrDNA) sequences. We compare the success of genome skims of 2051 herbarium specimens from Norway/Polar regions with 4604 freshly collected, silica gel dried specimens mainly from the European Alps and the Carpathians. Overall, we were able to assemble the full chloroplast genome for 67% of the samples and the full nrDNA cluster for 86%. Average insert length, cover and full cpDNA and rDNA assembly were considerably higher for silica gel dried than herbarium-preserved material. However, complete plastid genomes were still assembled for 54% of herbarium samples compared to 70% of silica dried samples. Moreover, there was comparable recovery of coding genes from both tissue sources (121 for silica gel dried and 118 for herbarium material) and only minor differences in assembly success of standard barcodes between silica dried (89% ITS2, 96% matK and rbcL) and herbarium material (87% ITS2, 98% matK and rbcL). The success rate was > 90% for all three markers in 1034 of 1036 genera in 160 families, and only Boraginaceae worked poorly, with 7 genera failing. Our study shows that large-scale genome skims are feasible and work well across most of the land plant families and genera we tested, independently of material type. It is therefore an efficient method for increasing the availability of plant biodiversity genomic data to support a multitude of downstream applications.

Testing multiple substrates for terrestrial biodiversity monitoring using environmental DNA metabarcoding

Biological surveys based on visual identification of the biota are challenging, expensive and time consuming, yet crucial for effective biomonitoring. DNA metabarcoding is a rapidly developing technology that can also facilitate biological surveys. This method involves the use of next generation sequencing technology to determine the community composition of a sample. However, it is uncertain as to what biological substrate should be the primary focus of metabarcoding surveys. This study aims to test multiple sample substrates (soil, scat, plant material and bulk arthropods) to determine what organisms can be detected from each and where they overlap. Samples (n = 200) were collected in the Pilbara (hot desert climate) and Swan Coastal Plain (hot Mediterranean climate) regions of Western Australia. Soil samples yielded little plant or animal DNA, especially in the Pilbara, probably due to conditions not conducive to long-term preservation. In contrast, scat samples contained the highest overall diversity with 131 plant, vertebrate and invertebrate families detected. Invertebrate and plant sequences were detected in the plant (86 families), pitfall (127 families) and vane trap (126 families) samples. In total, 278 families were recovered from the survey, 217 in the Swan Coastal Plain and 156 in the Pilbara. Aside from soil, 22%-43% of the families detected were unique to the particular substrate, and community composition varied significantly between substrates. These results demonstrate the importance of selecting appropriate metabarcoding substrates when undertaking terrestrial surveys. If the aim is to broadly capture all biota then multiple substrates will be required.

Museomics Unveil the Phylogeny and Biogeography of the Neglected Juan Fernandez Archipelago Megalachne and Podophorus Endemic Grasses and Their Connection With Relict Pampean-Ventanian Fescues

Oceanic islands constitute natural laboratories to study plant speciation and biogeographic patterns of island endemics. Juan Fernandez is a southern Pacific archipelago consisting of three small oceanic islands located 600-700 km west of the Chilean coastline. Exposed to current cold seasonal oceanic climate, these 5.8-1 Ma old islands harbor a remarkable endemic flora. All known Fernandezian endemic grass species belong to two genera, Megalachne and Podophorus, of uncertain taxonomic adscription. Classical and modern classifications have placed them either in Bromeae (Bromus), Duthieinae, Aveneae/Poeae, or Loliinae (fine-leaved Festuca); however, none of them have clarified their evolutionary relationships with respect to their closest Festuca relatives. Megalachne includes four species, which are endemic to Masatierra (Robinson Crusoe island) (M. berteroniana and M. robinsoniana) and to Masafuera (Alejandro Selkirk island) (M. masafuerana and M. dantonii). The monotypic Podophorus bromoides is a rare endemic species to Masatierra which is only known from its type locality and is currently considered extinct. We have used museomic approaches to uncover the challenging evolutionary history of these endemic grasses and to infer the divergence and dispersal patterns from their ancestors. Genome skimming data were produced from herbarium samples of M. berteroniana and M. masafuerana, and the 164 years old type specimen of P. bromoides, as well as for a collection of 33 species representing the main broad- and fine-leaved Loliinae lineages. Paired-end reads were successfully mapped to plastomes and nuclear ribosomal cistrons of reference Festuca species and used to reconstruct phylogenetic trees. Filtered ITS and trnTLF sequences from these genomes were further combined with our large Loliinae data sets for accurate biogeographic reconstruction. Nuclear and plastome data recovered a strongly supported fine-leaved Fernandezian clade where Podophorus was resolved as sister to Megalachne. Bayesian divergence dating and dispersal-extinction-cladogenesis range evolution analyses estimated the split of the Fernandezian clade from its ancestral southern American Pampas-Ventanian Loliinae lineage in the Miocene-Pliocene transition, following a long distance dispersal from the continent to the uplifted volcanic palaeo-island of Santa Clara-Masatierra. Consecutive Pliocene-Pleistocene splits and a Masatierra-to-Masafuera dispersal paved the way for in situ speciation of Podophorus and Megalachne taxa.

All-in-one sequencing: an improved library preparation method for cost-effective and high-throughput next-generation sequencing

BACKGROUND

Next generation sequencing (NGS) has been widely used in biological research, due to its rapid decrease in cost and increasing ability to generate data. However, while the sequence generation step has seen many improvements over time, the library preparation step has not, resulting in low-efficiency library preparation methods, especially for the most time-consuming and labor-intensive steps: size-selection and quantification. Consequently, there can be bottlenecks in projects with large sample cohorts.

RESULTS

We have described the all-in-one sequencing (AIO-seq) method, where instead of performing size-selection and quantification for samples individually, one sample one tube, up to 116 samples are pooled and analyzed in a single tube, 'All-In-One'. The AIO-seq method pools libraries based on the samples' expected data yields and the calculated concentrations of the size selected regions (target region), which can easily be obtained with the Agilent 2100 Bioanalyzer and Qubit Fluorometer. AIO-seq was applied to whole genome sequencing and RNA-seq libraries successfully, and it is envisaged that it could be applied to any type of NGS library, such as chromatin immunoprecipitation coupled with massively parallel sequencing, assays for transposase-accessible chromatin with high-throughput sequencing, and high-throughput chromosome conformation capture. We also demonstrated that for genetic population samples with low coverage sequences, like recombinant inbred lines (RIL), AIO-seq could be further simplified, by mixing the libraries immediately after PCR, without calculating the target region concentrations.

CONCLUSIONS

The AIO-seq method is thus labor saving and cost effective, and suitable for projects with large sample cohorts, like those used in plant breeding or population genetics research.