Skip to main content

​
Overview

PROTÉGÉ PD implements the PhyloTag approach described by Caro-Quintero et al. (2015), a systematic method for designing degenerate primers targeting conserved protein-coding genes for phylogenetic analysis and taxonomic identification.
The PhyloTag method leverages the dual nature of protein-coding genes: they are conserved enough at the amino acid level to allow universal primer design, yet variable enough at the nucleotide level to provide phylogenetic resolution.

​
Why Protein-Coding Genes?

Protein-coding genes serve as ideal phylogenetic markers for several key reasons:

​
Conservation at the Amino Acid Level

Functional proteins face selective pressure to maintain their amino acid sequences, especially in essential metabolic genes. This conservation allows primers to bind across diverse taxonomic groups.

​
Variation at the Nucleotide Level

Due to the degeneracy of the genetic code (multiple codons encode the same amino acid), synonymous substitutions accumulate over evolutionary time. These silent mutations provide the sequence variation needed for phylogenetic analysis.

​
Structured Evolution

Protein-coding genes evolve in a structured manner:
  • Third codon positions accumulate changes faster (wobble position)
  • First and second positions are more conserved
  • This predictable pattern aids alignment and phylogenetic inference

​
How PROTÉGÉ PD Implements PhyloTag

PROTÉGÉ PD follows a systematic workflow that translates the PhyloTag methodology into practical primer design:
The tool reads nucleotide sequences from a FASTA file and translates them to protein sequences:
# From protege.py:132-145
for seq_record in SeqIO.parse(genFile, 'fasta'):
    nuc = seq_record.seq._data.decode("utf-8")
    amino = nuc.translate()
This translation ensures that alignment occurs at the more conserved protein level.
Protein sequences are aligned using MUSCLE (Multiple Sequence Comparison by Log-Expectation):
muscle_lin -in translated_seqs.fas -out aligned_muscle.fas
Aligning at the protein level accounts for codon structure and avoids misalignments caused by synonymous substitutions.
The amino acid alignment is back-translated to nucleotide codons (protege.py:221-237), preserving alignment gaps while restoring the original DNA sequence:
Protein: M  A  T  -  K
Codons:  ATG GCT ACC --- AAA
This creates a codon-aware nucleotide alignment.
A consensus sequence is built position-by-position (protege.py:246-300), incorporating:
  • Nucleotide frequencies at each alignment position
  • Degeneracy codes when multiple nucleotides exceed the consensus threshold
  • Gap handling based on the --nogapconsensus flag
  • Consensus threshold (default 90%) for determining conserved positions
See Primer Degeneracy for details on how degeneracies are calculated.
The consensus sequence is scanned with a sliding window (default 21 nucleotides = 7 codons) to identify potential primer regions (protege.py:305-323).Each primer candidate is evaluated for:
  • Number of degeneracies
  • Absence of gaps
  • Melting temperature characteristics
See Melting Temperature for Tm calculation methods.

​
Key Parameters

ParameterFlagDefaultDescription
Consensus percentage-c / --consensus90%Minimum frequency for consensus calling
Gap handling-g / --nogapconsensusTrueWhether to allow gaps in consensus
Primer length-d / --codon7 codonsLength of primer in codons (21 nt)
Consensus threshold trade-off: Higher thresholds (over 90%) produce fewer degenerate primers but may miss conserved regions. Lower thresholds (under 80%) increase degeneracy but improve taxonomic coverage.

​
Scientific Foundation

The PhyloTag approach is grounded in the following publication: Caro-Quintero A, Konstantinidis KT. (2015) Improving the Taxonomic Resolution of Prokaryotic 16S rRNA Gene Databases by Pairwise Alignment. Genome Biology and Evolution, 7(12):3416-3424. https://academic.oup.com/gbe/article/7/12/3416/2467318 This methodology has been successfully applied to design primers for:
  • Bacterial housekeeping genes (gyrB, rpoB, recA)
  • Fungal phylogenetic markers
  • Environmental DNA metabarcoding

​
Advantages Over Traditional Methods

Compared to manual primer design or 16S rRNA-based approaches: ✓ Systematic and reproducible - algorithmic workflow eliminates subjective choices ✓ Higher phylogenetic resolution - protein-coding genes evolve faster than rRNA ✓ Broad taxonomic coverage - degenerate primers accommodate sequence variation ✓ Codon-aware alignment - respects genetic code structure ✓ Visual interface - interactive visualization helps select optimal primer pairs

​
Next Steps

To understand how PROTÉGÉ PD implements specific aspects of the PhyloTag approach:

Build docs developers (and LLMs) love

Get started for freeTalk to us