A Data-Analysis Pipeline for High-Throughput Systematic Evolution of Ligands by Exponential Enrichment (HT-SELEX) in the Characterization of Telomeric Proteins

Key findings from this study

• The study validated the analysis pipeline by applying it to the DNA-binding domain of human POT1, yielding results consistent with previous SELEX investigations.
• The authors demonstrate that all four C. elegans telomeric proteins examined—POT-1, POT-2, POT-3, and MRT-1—show a binding preference for G-enriched DNA sequences.
• The study found that C. elegans POT-1 additionally binds secondary structural elements in DNA, distinguishing it from the other three telomeric proteins analyzed.

Overview

This study develops and validates a comprehensive bioinformatics pipeline for analyzing high-throughput systematic evolution of ligands by exponential enrichment (HT-SELEX) data generated through next-generation sequencing. The pipeline was designed to characterize the nucleic acid-binding properties of telomeric proteins, specifically addressing the need for deeper analysis of SELEX datasets beyond what earlier methods provided. The research builds upon previous SELEX investigations of POT1, a component of the mammalian shelterin complex that protects single-stranded telomeric DNA and prevents inappropriate ATR kinase-mediated DNA damage responses at chromosome ends.

The investigation focuses on the DNA-binding domain of human POT1 (hPOT1) and extends to four Caenorhabditis elegans proteins considered homologs of hPOT1: POT-1, POT-2, POT-3, and MRT-1. Previous work had identified a binding pocket in POT1 termed the POT-hole, which recognizes the 5' phosphorylated dC at the telomeric double-stranded to single-stranded junction. The conservation of POT-hole residues across eukaryotes, including C. elegans, provided rationale for examining these homologous proteins using the newly developed analytical framework. The pipeline was constructed to be accessible and broadly applicable to diverse proteins that bind nucleic acids.

Methods and approach

The authors developed an extensive bioinformatics analysis pipeline specifically designed for processing SELEX datasets generated by next-generation sequencing platforms. The pipeline was validated through application to the DNA-binding domain of human POT1, serving as a benchmark against results from a previous SELEX study to ensure methodological consistency and accuracy. This validation step established the reliability of the computational approach before extending analysis to less-characterized proteins.

The pipeline was subsequently applied to characterize the DNA-binding activity of four C. elegans telomeric proteins: POT-1, POT-2, POT-3, and MRT-1. These proteins were selected based on their designation as homologs of human POT1, allowing comparative analysis of DNA-binding preferences across evolutionarily divergent organisms. The analytical framework was designed to enable deeper sequence analysis than previous methods, extracting more comprehensive information about binding specificity and preferences from the enriched sequence pools generated through SELEX rounds.

Results

Application of the pipeline to the DNA-binding domain of human POT1 yielded results consistent with previous SELEX investigations, validating the accuracy and reliability of the new analytical approach. This concordance with established findings demonstrated that the pipeline could reproduce known binding characteristics while potentially revealing additional insights from the same dataset through more extensive analysis.

Analysis of the four C. elegans telomeric proteins revealed distinct binding preferences. All four proteins—POT-1, POT-2, POT-3, and MRT-1—exhibited a binding preference for G-enriched DNA sequences, suggesting a common recognition motif among these telomeric proteins. POT-1 from C. elegans displayed an additional binding capacity not observed in the other three proteins: recognition of secondary structural elements in DNA. This distinction indicated that while these proteins share homology and certain binding characteristics, they have evolved specialized recognition properties that may reflect distinct functional roles at telomeres.

Implications

The development of this bioinformatics pipeline addresses a methodological need in the characterization of protein-nucleic acid interactions using HT-SELEX data. The validation with human POT1 and successful application to C. elegans telomeric proteins demonstrates that the pipeline can be deployed across different protein systems and organisms, providing a standardized approach for extracting comprehensive binding information from SELEX experiments. The accessibility of the pipeline, as emphasized by the authors, positions it as a tool that can be adopted by laboratories investigating diverse nucleic acid-binding proteins beyond telomeric factors.

The findings regarding C. elegans telomeric proteins contribute to understanding telomere protection mechanisms in non-mammalian eukaryotes. The shared preference for G-enriched sequences among POT-1, POT-2, POT-3, and MRT-1 aligns with the guanine-rich composition of telomeric DNA, while POT-1's additional capacity to recognize secondary structures suggests functional specialization. These results provide a foundation for investigating how these proteins coordinate to protect C. elegans telomeres and whether the mechanistic principles differ from the mammalian shelterin complex. The conservation of certain binding characteristics alongside divergent features suggests both conserved and lineage-specific aspects of telomere protection across eukaryotic evolution.

Scope and limitations

This summary is based on the study abstract and available metadata. It does not include a full analysis of the complete paper, supplementary materials, or underlying datasets unless explicitly stated. Findings should be interpreted in the context of the original publication.