Fun-sized library prep: Miniaturization is a valid method for per-sample cost reduction in for targeted sequencing of antique and fresh and angiosperm DNA
Data files
Jul 16, 2026 version files 480.29 KB
-
BeadCleanUp1_JohnsonLab_10thVol.metadata
106 B
-
BeadCleanUp1_JohnsonLab_10thVol.protocol
133.83 KB
-
BeadCleanUp1_Reagents_JohnsonLab.dfx
1.36 KB
-
BeadCleanUp2_JohnsonLab_10thVol.metadata
106 B
-
BeadCleanUp2_JohnsonLab_10thVol.protocol
125.53 KB
-
BeadCleanUp2_Reagents_JohnsonLab.dfx
1.36 KB
-
EndPrep_AdaptorLigation_DF_noFS_JohnsonLab.dfx
2.38 KB
-
EndPrep_AdaptorLigation_DF_withFS_JohnsonLab.dfx
2.38 KB
-
FunSizedLibPrep_RCode.R
21.13 KB
-
FunSizedLibPrep_StatsFromSequencing_OverlapSet.csv
5.08 KB
-
FunSizedLibPrep_StatsFromSequencing.csv
163.38 KB
-
IndexPCR_Reagents_JohnsonLab.dfx
712 B
-
README.md
16.26 KB
-
Supp_TABLE_plastid_hybpiper_stats_combined.csv
6.70 KB
Abstract
Premise: Genomic analysis of population structure is important to the conservation of plant species of concern. A limitation of using genetic information in conservation is the cost of obtaining large datasets. Targeted sequencing and low-volume robotic liquid handlers can reduce library preparation reaction volumes and costs.
Methods: We used targeted sequencing via Angiosperms353 to obtain data for 768 samples, 18 of which were identical, at 0.5X and 0.1X reaction volumes. We calculated quality and quantity control statistics to compare the effects of tissue age and library reaction volume on sequencing results for on-target nuclear and off-target plastid genes.
Results: Library miniaturization to 0.1X reduces costs and generally performs comparably to 0.5X libraries. In the full dataset, Tenth-Volume Only libraries showed smaller insert sizes and fewer genes with mapped sequences, but no reduction in mapped reads. In the Overlap Set, 0.1X had equal or improved performance with no significant decrease in sequencing efficiency. Differences by tissue type likely reflected sampling variation.
Conclusions: Miniaturization to 0.1X substantially reduces per-sample costs while maintaining comparable sequencing quality across fresh and herbarium angiosperm DNA. Overall, miniaturization provides a reliable, cost-effective approach for targeted sequencing, increasing the feasibility of using herbarium collections and enabling broader access to population-level genomic studies.
Dataset DOI: 10.5061/dryad.4mw6m90qf
Description of the data and file structure
This dataset contains DNA sequencing data from 1,518 angiosperm specimens (243 species, 176 genera, 63 families), including fresh and herbarium tissue. DNA was extracted and dual-indexed libraries were prepared at 0.5X and 0.1X reaction volumes using automated liquid handling robots. Libraries were enriched with the Angiosperms353 probe set and sequenced on the Illumina NovaSeq 6000 platform. The dataset includes post-sequencing metrics such as insert fragment size, reads mapped, and gene recovery, along with library preparation cost and time data.
Files and Variables
File: FunSizedLibPrep_StatsFromSequencing.csv
Description: This file contains post-sequencing metrics generated from the HybPiper pipeline for angiosperm DNA libraries prepared at half-volume (0.5X) and tenth-volume (0.1X) reaction volumes. Metrics summarize library performance, including read counts, target capture efficiency, gene recovery, and quality control warnings.
Variables:
- Name: Unique sequence ID
- barcode: Unique herbarium-assigned barcode ID
- libVol: Library preparation volume category (
0.5Xor0.1X) - tissueType: Tissue origin (
freshorherbarium) - tissueVol: Combination of libVol and tissueType
- insertPeak: Average insert fragment size in base pairs
- NumReads: Total number of sequencing reads generated
- ReadsMapped: Number of reads mapped to targeted regions
- PctOnTarget: Percentage of reads mapping to target regions
- GenesMapped: Number of target genes with at least one read mapped
- GenesWithContigs: Number of genes assembled into contigs
- GenesWithSeqs: Number of genes with recovered sequences
- GenesAt25pct: Number of genes recovered at ≥25% of target length
- GenesAt50pct: Number of genes recovered at ≥50% of target length
- GenesAt75pct: Number of genes recovered at ≥75% of target length
- GenesAt150pct: Number of genes recovered at ≥150% of target length
- ParalogWarningsLong: Number of genes flagged for potential paralogy based on long sequences
- ParalogWarningsDepth: Number of genes flagged for potential paralogy based on read depth
- GenesWithoutStitchedContigs: Number of genes where contigs could not be stitched
- GenesWithStitchedContigs: Number of genes successfully stitched into contigs
- GenesWithStitchedContigsSkipped: Number of genes skipped during stitching
- GenesWithChimeraWarning: Number of genes flagged for potential chimeric assemblies
- Missing Values: N/A
File: FunSizedLibPrep_StatsFromSequencing_OverlapSet.csv
Description: This file contains post-sequencing metrics generated from the HybPiper pipeline for angiosperm DNA, specifically for a subset of samples prepared into libraries at both half-volume (0.5X) and tenth-volume (0.1X) reaction volumes. Metrics summarize library performance including read counts, target capture efficiency, gene recovery, and quality control warnings.
Variables:
- Name: Unique sequence ID
- barcode: Unique herbarium-assigned barcode ID
- libVol: Library preparation volume category (
0.5Xor0.1X) - tissueType: Tissue origin (
freshorherbarium) - tissueVol: Combination of libVol and tissueType
- insertPeak: Average insert fragment size in base pairs
- NumReads: Total number of sequencing reads generated
- ReadsMapped: Number of reads mapped to targeted regions
- PctOnTarget: Percentage of reads mapping to target regions
- GenesMapped: Number of target genes with at least one read mapped
- GenesWithContigs: Number of genes assembled into contigs
- GenesWithSeqs: Number of genes with recovered sequences
- GenesAt25pct: Number of genes recovered at ≥25% of target length
- GenesAt50pct: Number of genes recovered at ≥50% of target length
- GenesAt75pct: Number of genes recovered at ≥75% of target length
- GenesAt150pct: Number of genes recovered at ≥150% of target length
- ParalogWarningsLong: Number of genes flagged for potential paralogy based on long sequences
- ParalogWarningsDepth: Number of genes flagged for potential paralogy based on read depth
- GenesWithoutStitchedContigs: Number of genes where contigs could not be stitched
- GenesWithStitchedContigs: Number of genes successfully stitched into contigs
- GenesWithStitchedContigsSkipped: Number of genes skipped during stitching
- GenesWithChimeraWarning: Number of genes flagged for potential chimeric assemblies
- TotalBasesRecovered_Genes: Total number of base pairs recovered across all targeted coding regions for a sample, excluding flanking non-coding regions
- TotalBasesRecovered_SC: Total number of base pairs recovered across all target genes for a sample, including exons and flanking regions
- Missing Values: N/A
File: Supp_TABLE_plastid_hybpiper_stats_combined.csv
Description: This file contains post-sequencing metrics generated from the HybPiper pipeline for angiosperm DNA, specifically for a subset of samples prepared into libraries at both half-volume (0.5X) and tenth-volume (0.1X) reaction volumes. These metrics specifically pertain to plastid gene recovery from off-target reads, and summarize library performance including read counts, target capture efficiency, gene recovery, and quality control warnings.
Variables:
- Name: Unique sequence ID
- Volume: Library preparation reaction volume
- NumReads: Total number of sequencing reads generated
- ReadsMapped: Number of reads mapped to targeted regions
- PctOnTarget: Percentage of reads mapping to target regions
- GenesMapped: Number of target genes with at least one read mapped
- GenesWithContigs: Number of genes assembled into contigs
- GenesWithSeqs: Number of genes with recovered sequences
- GenesAt25pct: Number of genes recovered at ≥25% of target length
- GenesAt50pct: Number of genes recovered at ≥50% of target length
- GenesAt75pct: Number of genes recovered at ≥75% of target length
- GenesAt150pct: Number of genes recovered at ≥150% of target length
- ParalogWarningsLong: Number of genes flagged for potential paralogy based on long sequences
- ParalogWarningsDepth: Number of genes flagged for potential paralogy based on read depth
- GenesWithoutStitchedContigs: Number of genes where contigs could not be stitched
- GenesWithStitchedContigs: Number of genes successfully stitched into contigs
- GenesWithStitchedContigsSkipped: Number of genes skipped during stitching
- GenesWithChimeraWarning: Number of genes flagged for potential chimeric assemblies
- TotalBasesRecovered: Total number of base pairs recovered across all plastid regions for a sample, excluding flanking non-coding regions
- Missing Values: N/A
Code
File: FunSizedLibPrep_RCode.R
Description: This R code analyzes sequencing library performance by comparing half-volume (0.5X) and tenth-volume (0.1X) reactions, as well as tissue type (fresh vs. herbarium). It performs t-tests, one-way and two-way ANOVAs, and post-hoc Tukey HSD tests on metrics including average insert fragment size, total reads mapped, and number of genes with recovered sequences. The results are visualized using violin and box plots, with outliers highlighted, to assess differences across library volumes, tissue types, and their interactions.
All data analyses and visualizations were performed using R (R Project for Statistical Computing, version 4.3.1 or higher), a free and open-source statistical programming environment (https://www.r-project.org).
Required R packages:
ggplot2,gridExtra,ggpubr– for plots and figure layouttidyverse(includesdplyr,readr,tidyr) – for data manipulationbroom– to tidy statistical outputsAICcmodavg– for model selection
rstatix– for ANOVA and post-hoc tests
Software
Instrument Software Requirements
The automated liquid handling files included in this dataset are proprietary method files generated for SPT Labtech instruments.
.dfxfiles are method files for the SPT Labtech Dragonfly Discovery liquid handler and should be opened using Dragonfly control software..protocolfiles are method files for the SPT Labtech Mosquito HV liquid handler and should be opened using the Mosquito control software..metadatafiles accompany Mosquito protocol files and contain instrument configuration, consumable definitions, and plate layout information required by the Mosquito software. These files should remain paired with their associated.protocolfiles.
These files are provided to allow replication of the automated library preparation workflow described in the manuscript and are generally not intended to be viewed or edited using standard text editors (though .metadata files can be).
File: BeadCleanUp1_Reagents_JohnsonLab.dfx
Description: Dragonfly reagent dispensing method for the first SPRI bead cleanup following enzymatic library preparation. The file specifies the transfer of bead cleanup reagents (e.g., magnetic beads, 80% ethanol, and elution buffer) from source reservoirs to destination wells of a reagent plate according to the miniaturized protocol. The file defines dispense volumes, liquid classes, and plate layout required for this cleanup step.
Usage Notes: This file is intended for use with the SPT Labtech Dragonfly Discovery liquid handling system. The .dfx file is imported into the proprietary Dragonfly software to execute the reagent dispensing method. The file is not intended to be viewed or edited in a standard text editor.
File: BeadCleanUp2_Reagents_JohnsonLab.dfx
Description: Dragonfly reagent dispensing method for the second SPRI bead cleanup following PCR enrichment. The file specifies the transfer of bead cleanup reagents (e.g., magnetic beads, 80% ethanol, and elution buffer) from source reservoirs to destination wells of a reagent plate according to the miniaturized protocol. The file defines dispense volumes, liquid classes, and plate layout required for this cleanup step.
Usage Notes: This file is imported into the SPT Labtech Dragonfly Discovery software to execute the reagent dispensing workflow for the second SPRI bead cleanup. The .dfx file is imported into the proprietary Dragonfly software to execute the reagent dispensing method. The file is not intended to be viewed or edited in a standard text editor.
File: EndPrep_AdaptorLigation_DF_noFS_JohnsonLab.dfx
Description: Dragonfly dispensing method for the combined End Prep and Adapter Ligation workflow using the standard NEBNext Ultra II DNA Library Prep Kit (no-fragmentase version). The method dispenses the enzyme and buffer components required for end repair, dA-tailing, and adapter ligation into each sample well using miniaturized reaction volumes. This workflow is intended for DNA that has already been fragmented prior to library preparation.
Usage Notes: This file is opened within the SPT Labtech Dragonfly Discovery software and executed on a Dragonfly instrument. It contains the dispensing instructions for the miniaturized End Prep and Adapter Ligation workflow. The .dfx file is imported into the proprietary Dragonfly software to execute the reagent dispensing method. The file is not intended to be viewed or edited in a standard text editor.
File: EndPrep_AdaptorLigation_DF_withFS_JohnsonLab.dfx
Description: Dragonfly dispensing method for the combined End Prep and Adapter Ligation workflow using the NEBNext Ultra II FS DNA Library Prep Kit. In addition to end repair and adapter ligation reagents, this workflow accommodates the enzymatic fragmentation chemistry provided by the Ultra II FS kit. The file defines reagent transfers necessary for processing intact genomic DNA in a miniaturized reaction format.
Usage Notes: This file is imported into the SPT Labtech Dragonfly Discovery software for execution on a Dragonfly liquid handler. The .dfx file is imported into the proprietary Dragonfly software to execute the reagent dispensing method. The file is not intended to be viewed or edited in a standard text editor.
File: IndexPCR_Reagents_JohnsonLab.dfx
Description: Dragonfly method for dispensing reagents required for the library amplification (index PCR) step. The file distributes PCR master mix into sample wells prior to amplification while maintaining the reduced reaction volumes used in the miniaturized protocol. Sample-specific index primers are added separately via the Mosquito according to the experimental design.
Usage Notes: This file is executed using the SPT Labtech Dragonfly Discovery software on a Dragonfly liquid handler to dispense PCR reagents according to the miniaturized protocol. The .dfx file is imported into the proprietary Dragonfly software to execute the reagent dispensing method. The file is not intended to be viewed or edited in a standard text editor.
File: BeadCleanUp1_JohnsonLab_10thVol.protocol
Description: Mosquito liquid handling protocol implementing the first SPRI bead cleanup. The protocol contains the aspirate, dispense, and mixing commands executed by the Mosquito instrument to perform magnetic bead purification at one-tenth of the manufacturer's recommended reaction volume. The protocol is intended to be used following Dragonfly reagent dispensing.
Usage Notes: This file is used with the SPT Labtech Mosquito HV liquid handling system and is opened in the proprietary Mosquito software. The .protocol file contains the instrument commands executed during the run and is not intended to be viewed or edited in a standard text editor. It is intended to be used together with the associated .metadata file, which can be viewed in a standard text editor.
File: BeadCleanUp2_JohnsonLab_10thVol.protocol
Description: Mosquito liquid handling protocol implementing the second SPRI bead cleanup following index PCR amplification. The protocol automates bead transfers, washing, and elution steps using the miniaturized workflow and is designed to maximize library recovery while maintaining consistent handling across samples.
Usage Notes: This file is used with the SPT Labtech Mosquito HV liquid handling system and is opened in the proprietary Mosquito software. The .protocol file contains the instrument commands executed during the run and is not intended to be viewed or edited in a standard text editor. It is intended to be used together with the associated .metadata file, which can be viewed in a standard text editor.
File: BeadCleanUp1_JohnsonLab_10thVol.metadata
Description: Metadata file associated with the first Mosquito bead cleanup protocol. This file stores instrument configuration information, consumable definitions, plate layouts, and protocol settings required by the Mosquito software to execute the corresponding cleanup protocol correctly.
Usage Notes: This file accompanies the corresponding .protocol file and is automatically read by the SPT Labtech Mosquito software. It stores supporting configuration information required for successful execution of the protocol and is not intended for manual editing, but can be viewed in a standard text editor.
File: BeadCleanUp2_JohnsonLab_10thVol.metadata
Description: Metadata file associated with the second Mosquito bead cleanup protocol. The file contains the run configuration parameters and supporting information needed by the Mosquito software to accompany the second cleanup protocol, ensuring correct interpretation of labware, pipetting instructions, and workflow settings.
Usage Notes: This file accompanies the corresponding .protocol file and is automatically read by the SPT Labtech Mosquito software. It stores supporting configuration information required for successful execution of the protocol and is not intended for manual editing, but can be viewed in a standard text editor.
