VTBuilder

VTBuilder assembles non-chimeric transcripts from complex sequencing read datasets to reconstruct isoform diversity and preserve relationships between co-evolving sites.

Key Features:

• Non-Chimeric Transcript Construction: Constructs accurate non-chimeric transcripts from read data while maintaining isoform integrity.
• Handling Isoform Complexity: Manages datasets with high isoform variation, applicable to multi-isoform transcriptomes and rapidly evolving viral populations such as HIV-1.
• Preservation of Co-evolving Sites: Maintains relationships between co-evolving sites within reconstructed transcripts.
• Validation against Sanger SSTs and Comparison to Trinity: Validated on simulated reads from 54 Sanger sequenced transcripts (SSTs), achieving >99% sequence similarity to 48 SSTs and outperforming Trinity in the same comparison.
• Length-Distribution Fidelity: Produces transcript length distributions closely matching SSTs, in contrast to shorter outputs observed from Trinity.
• Scalability and Sequencing Input Support: Applied to assemble 5 million Illumina MiSeq paired-end reads into 1481 transcripts and recovered major toxin genes from a venom gland transcriptome.

Scientific Applications:

• Transcriptomics: Reconstructs complex sequence and isoform diversity in transcriptomes, including snake venom gland datasets.
• Phylogenetics: Preserves co-evolving site relationships to support evolutionary and phylogenetic analyses.
• Viral Population Studies: Applicable to assembly and variant reconstruction in rapidly evolving viral populations such as HIV-1.
• Drug-Resistance and Parasite Studies: Supports analysis of genetic variation in drug-resistant parasite populations.
• Toxin Gene Recovery: Enables reconstruction and recovery of major toxin genes from venom gland sequencing data (e.g., Bitis arietans, Echis ocellatus).

Methodology:

Validated using simulated reads from 54 Sanger sequenced transcripts (SSTs) expressed in the venom gland of Echis ocellatus (constructed 55 transcripts with >99% sequence similarity to 48 SSTs and compared to Trinity which constructed 14 matching transcripts), evaluated transcript length distributions, and applied assembly to 5 million Illumina MiSeq paired-end reads from the venom gland of Bitis arietans yielding 1481 transcripts and recovery of major toxin genes.