DNAharvester: A Nextflow Pipeline for Analysing Highly Degraded DNA from Ancient and Historical Specimens

DNAharvester is a modular, Nextflow-based pipeline designed to address the bioinformatic challenges of analyzing highly degraded ancient DNA by integrating advanced filtering, flexible mapping strategies, and comprehensive downstream workflows to maximize authentic data recovery while minimizing contamination and reference bias.

The Gist

Imagine you are a detective trying to solve a mystery from thousands of years ago. Your only clues are tiny, crumbled scraps of paper (DNA) found in an old, dusty attic (ancient bones or artifacts). These scraps are in terrible shape: they are torn into tiny pieces, covered in mold, mixed with trash from other people who visited the room, and some of the ink has faded or changed color over time.

If you try to read these clues using a standard method, you might accidentally piece together the wrong story because the clues are so messy. This is the biggest headache for scientists studying ancient DNA.

Enter "DNAharvester," the new super-tool described in this paper.

Think of DNAharvester as a high-tech, automated detective squad built specifically to handle these messy, ancient clues. Here is how it works, using simple analogies:

1. The "Gold Panning" Filter

Before the detective even looks at the clues, they need to separate the gold from the dirt. Ancient samples are full of "dirt" (bacteria, fungi, and modern human contamination). DNAharvester has a special metagenomic filter that acts like a fine-mesh sieve. It washes away the junk and keeps only the precious ancient DNA, ensuring the scientists aren't wasting time on fake leads.

2. The "Chameleon" Map Reader

Once the clues are clean, they need to be matched to a master map (a reference genome). But ancient DNA is so broken that a standard map reader might force a torn piece into the wrong spot just to make it fit.
DNAharvester is smart enough to change its strategy. It can switch between different "map readers" (like BWA-aln, BWA-mem, or Bowtie2) depending on how broken the clues are. It's like having a team of translators who can speak different dialects; if one translator can't understand the broken text, another one steps in to make sure the clues are placed in the right spot, not just any spot.

3. The "Truth Detector"

Sometimes, the map itself is biased—it might make the ancient DNA look more like the modern map than it really is. DNAharvester has a built-in lie detector. It constantly checks its own work to make sure it isn't accidentally "forcing" the ancient clues to look like the modern ones, ensuring the final story is authentic.

4. The "Swiss Army Knife" of Analysis

Once the clues are sorted and mapped, DNAharvester doesn't stop. It has a whole toolbox of sub-tasks ready to go:

• Rebuilding the Family Tree: It can stitch together the tiny scraps to rebuild the full mitochondrial DNA (the family history passed down from mothers).
• Identifying the Invaders: It looks at the "trash" that was filtered out to see if ancient pathogens (like plague or tuberculosis) were hiding there.
• Guessing the Gender: It can tell if the ancient person was male or female just by looking at specific chromosome clues.
• Reading the Faded Ink: Since the DNA is so damaged, it uses different math tricks to guess what the original genetic code said, even if the letters are half-erased.

5. The "Plug-and-Play" Factory

The best part? This whole operation runs on Nextflow, which is like a universal power adapter for computers.

• Scalable: It can run on a small laptop for a single bone or on a massive supercomputer for a thousand bones.
• Portable: It's "containerized," meaning all the tools are packed in a box that works on any computer, no matter what software is installed.
• Reproducible: If another scientist uses the same settings, they get the exact same result. No more "it worked on my computer" excuses.

The Bottom Line

Before DNAharvester, analyzing ancient DNA was like trying to fix a shattered vase with duct tape and guesswork. It was messy, prone to errors, and hard to repeat.

DNAharvester turns this into a streamlined, automated assembly line. It takes the most broken, contaminated, and difficult ancient samples and turns them into clear, reliable genetic stories. It makes high-level ancient DNA research accessible to everyone, not just the bioinformatics wizards, ensuring that the voices of our ancestors are heard clearly, even after thousands of years of silence.

Technical Summary

1. Problem Statement

The field of ancient DNA (aDNA) research faces significant bioinformatic hurdles when processing highly degraded DNA from palaeontological and archaeological specimens. Unlike modern DNA, aDNA is characterized by:

• Short fragment lengths and low endogenous content (the proportion of DNA belonging to the target organism).
• Post-mortem damage (chemical modifications that mimic mutations).
• High levels of cross-species contamination (environmental microbes, modern human handling).

These factors lead to spurious mapping (reads mapping incorrectly to the reference genome) and reference bias (systematic preference for reference alleles), which severely compromise the accuracy of downstream population genetic inferences. Existing tools often lack the flexibility to handle this specific combination of challenges simultaneously.

2. Methodology

The authors developed DNAharvester, a modular, reproducible pipeline built on the Nextflow workflow management system. The pipeline employs a multi-stage strategy designed to maximize authentic data recovery:

• Pre-mapping Processing:

  • Metagenomic Filtering: Reads are screened against a database of common contaminants before mapping to remove non-target sequences early.
  • Competitive Mapping: Reads are mapped against multiple references to distinguish between the target genome and potential contaminants.
  • Adaptive Aligner Selection: The pipeline dynamically selects the most appropriate alignment algorithm based on the data characteristics, supporting BWA-aln, BWA-mem, and Bowtie2.

• Mapping and Bias Mitigation:

  • Implements systematic evaluation and correction for reference bias and spurious mapping.
  • Utilizes flexible mapping and filtering strategies that can be tuned based on the specific preservation state of the sample.

• Specialized Subworkflows:

  • Iterative Assembly: For reconstructing mitogenomes.
  • Feature Identification: Detects genomic repeats and CpG sites.
  • Taxonomic Classification: Classifies unmapped reads to identify microbial or pathogen presence.
  • Genetic Sex Determination: Analyzes coverage ratios on sex chromosomes.
  • Variant Calling: Offers three distinct strategies to accommodate varying sequencing depths:
    1. Diploid variant calling.
    2. Genotype likelihood estimation.
    3. Pseudo-haploid random allele calling (crucial for low-coverage data).

• Technical Implementation:

  • The entire framework is containerized (likely using Docker/Singularity) to ensure portability and environment consistency.
  • Configuration is managed through simple files, allowing researchers to adapt the pipeline without deep coding knowledge.

3. Key Contributions

• Specialized Architecture: Unlike general-purpose pipelines, DNAharvester is explicitly designed for the unique constraints of highly degraded DNA, prioritizing the maximization of authentic signal recovery.
• Modularity and Flexibility: The pipeline allows researchers to swap aligners, adjust filtering thresholds, and choose variant calling methods based on sample quality and sequencing depth.
• Comprehensive Analysis Suite: It consolidates diverse analytical tasks (from pathogen screening to mitogenome assembly) into a single, cohesive workflow, reducing the need for multiple disparate tools.
• Reproducibility and Scalability: By leveraging Nextflow and containerization, the pipeline ensures that analyses are reproducible across different computing environments and can scale from small datasets to large population studies.

4. Results

The authors validated DNAharvester through:

• Simulated Scenarios: Testing across a gradient of degradation levels and contamination rates.
• Empirical Datasets: Application to real-world ancient and historical specimens.

Key Findings:

• The pipeline successfully mitigated complex background contamination, effectively distinguishing authentic endogenous DNA from environmental noise.
• It demonstrated the ability to preserve authentic genomic signals even in the most challenging, highly degraded samples where other methods might fail.
• The flexible configuration allowed for optimal performance across varying sequencing depths, ensuring robust variant calling regardless of data coverage.

5. Significance

DNAharvester represents a major step forward in standardizing aDNA analysis. By streamlining complex bioinformatic tasks into a user-friendly, configuration-driven framework, it:

• Democratizes Access: Makes rigorous genomic analysis of ancient remains accessible to a broader research community, including those without extensive bioinformatics expertise.
• Enhances Data Quality: Provides a standardized approach to minimize reference bias and contamination errors, leading to more reliable population genetic inferences.
• Future-Proofs Research: Its modular design ensures it can adapt to new sequencing technologies and analytical algorithms as the field of paleogenomics continues to evolve.