Subtelomere

A subtelomere is the dynamic, structurally polymorphic transition domain located on linear eukaryote chromosomes immediately proximal to the terminal (TTAGGG)_n hexameric telomere repeats and distal to chromosome-specific single-copy euchromatin. Spanning approximately 100 to 500 kilobases (kb) at each chromosome arm terminus in the human genome (such as Chr1p36.33, Chr7p22.1, Chr11p15.5, Chr16p13.3, and Chr19p13.3), subtelomeres are mosaic patchworks of GC-rich segmental duplications, variable number tandem repeats (VNTRs), degenerate telomere-like interstitial tracts (ITS), and rapidly evolving multigene families.

Unlike the conserved single-copy core of mammalian chromosomes, subtelomeric regions experience elevated rates of non-allelic homologous recombination (NAHR), sister-chromatid exchange, and double-strand break repair. Across human and microbial eukaryotes, subtelomeres are enriched in non-canonical secondary structures—specifically parallel G-quadruplex (G4) clusters and alternating purine–pyrimidine Z-DNA tracts—that partition replication timing, nucleate telomeric repeat-containing RNA (TERRA) transcription, and organize three-dimensional chromatin loops. In parasitic protozoa and pathogenic fungi, subtelomeric instability serves as an evolutionary engine for antigenic variation and host immune evasion, whereas in humans, subtelomeric rearrangements cause clinically significant neurodevelopmental microdeletion syndromes and facioscapulohumeral muscular dystrophy (FSHD).

Distal and Proximal Subtelomeric Subdomains

High-resolution long-read assemblies, culminating in the Telomere-to-Telomere (T2T-CHM13) and GRCh38 reference genomes, divide human subtelomeres into two structurally distinct compartments: a distal subtelomeric subdomain (extending ~10 to 100 kb inward from the terminal (TTAGGG)_n array) and a proximal subtelomeric subdomain (extending an additional ~100 to 400 kb to the boundary of unique chromosome-specific sequence). The distal zone contains short monomer repeats (Srpt), CpG islands, and subtelomeric promoters equipped with 29-bp, 37-bp, and 61-bp repeat elements that drive transcription toward the chromosome tip to produce the long non-coding RNA TERRA.

Both distal and proximal zones contain blocks of subtelomeric segmental duplications (paralogous blocks of 10–50 kb sharing 90–99.8% sequence identity across multiple non-homologous chromosome ends, such as 1p, 5q, 6p, 16p, and 19p). Interspersed between these paralogous blocks are interstitial telomeric sequences (ITS)—head-to-tail or head-to-head relics of ancient telomere-to-telomere fusion events and break-induced replication repairs—alongside dense GC-rich macrosatellites such as the 3.3-kb D4Z4 repeat array on Chr4q35 and Chr10q26.

Dual Enrichment of G-Quadruplex and Z-DNA Motifs

Comparative multi-chromosome window scans across the distal 500 kb of human autosomes (GRCh38 Chr1p36.33, Chr7p22.1, Chr11p15.5, Chr16p13.3, and Chr19p13.3) demonstrate a strong enrichment of non-B DNA conformations relative to interior chromosomal arms. Elevated baseline GC content (often 55–68% GC compared to the 41% human genome average) combines with pronounced strand-specific G/C skew (G4Hunter scores exceeding 1.50 to 2.10) to generate dense clusters of intramolecular G-quadruplex folds. Simultaneously, intronic and intergenic microsatellites within subtelomeric genes—exemplified by AGRN Introns 18–22 on Chr1p36.33 and AXIN1/HBA1 flanking tracts on Chr16p13.3—accumulate extended (TG)_n and (CG)_n dinucleotide arrays that transition into left-handed Z-DNA under negative transcriptional supercoiling.

Telomere Position Effect and Chromatin Looping

In budding yeast (Saccharomyces cerevisiae), Drosophila, and mammalian cells, genes positioned close to telomeres undergo stochastic, heritable transcriptional silencing termed the Telomere Position Effect (TPE). Classical TPE spreads inward for several kilobases as shelterin-recruited histone deacetylases (SIRT6 in mammals; the Sir2–Sir3–Sir4 complex in yeast) deacetylate histone tails and methylate H3K9 (H3K9me3) and H4K20 (H4K20me3), recruiting Heterochromatin Protein 1 (HP1α).

In human chromosomes, where subtelomeres span hundreds of kilobases, long-range three-dimensional chromatin-remodeling and CTCF–cohesin loop extrusion extend this regulation via Telomere Position Effect Over Long Distances (TPE-OLD). Long telomeres loop back to contact subtelomeric loci up to 1–5 Mb away (such as ISG15 on Chr1p36.33 or TERT on Chr5p15.33), maintaining repressive heterochromatin. As somatic cells divide and telomeres shorten below a threshold length, these long-range subtelomeric loops disengage before the onset of a classic DNA-damage senescence response, upregulating subtelomeric genes involved in innate interferon signaling and cellular stress.

TERRA Transcription and R-Loop Maintenance

Initiating within the CpG-rich 29-bp/61-bp repeat promoters of the distal subtelomere, RNA Polymerase II transcribes through the subtelomeric-telomeric junction to synthesize TERRA (UUAGGG_n), a heterogeneous long non-coding RNA ranging from 100 nt to >9 kb. TERRA transcripts remain tethered to chromosome ends both by folding into parallel RNA G-quadruplexes recognized by the shelterin subunit TRF2 (via its basic GAR domain) and by invading the duplex telomeric repeat tract to form three-stranded RNA:DNA hybrids (R-loops). Controlled subtelomeric TERRA R-loops recruit the histone methyltransferase SUV39H1 and the Polycomb Repressive Complex 2 (PRC2) to maintain subtelomeric compaction, and promote homologous-recombination-mediated Alternative Lengthening of Telomeres (ALT) in telomerase-negative sarcomas and glioblastomas.

Antigenic Variation in Eukaryotic Pathogens

Because subtelomeres reside outside the syntenic core where meiotic crossover suppression protects housekeeping loci, natural evolution has co-opted subtelomeric instability as a generator of surface antigen diversity. In the human malaria parasite Plasmodium falciparum, ~60 paralogous var genes encoding P. falciparum Erythrocyte Membrane Protein 1 (PfEMP1) are concentrated within the subtelomeres of all 14 chromosomes. Epigenetic monoallelic exclusion ensures that only a single var locus is transcribed inside an active perinuclear zone while the remaining 59 subtelomeric copies are sequestered in H3K9me3/PfHP1-silenced clusters at the nuclear periphery; frequent ectopic gene conversion between clustered subtelomeric ends generates chimeric PfEMP1 variants that evade host antibody responses.

An analogous subtelomeric archive operates in the African sleeping sickness trypanosome Trypanosoma brucei, which stores more than 1,500 silent Variant Surface Glycoprotein (VSG) genes and pseudogenes inside subtelomeric arrays and minichromosomes, periodically copying a silent subtelomeric VSG cassette into the active Bloodstream Expression Site (BES) via break-induced replication. Likewise, free-living and pathogenic yeasts (Saccharomyces MAL/SUC/FLO families; Candida glabrata EPA adhesins; Pneumocystis MSG surface glycoproteins) concentrate environment-specific metabolic and adhesin gene families within subtelomeres to enable rapid adaptive radiation.

Clinical Subtelomeric Rearrangements and FSHD

In human clinical genetics, cryptic terminal deletions, unbalanced translocations, and inverted duplications involving subtelomeres account for 5% to 7% of unexplained moderate-to-severe intellectual disability and multiple congenital anomaly syndromes. Because high GC content and segmental duplications obscure subtelomeres on conventional G-banded karyotypes (where terminal R-bands appear pale), diagnosis relies on chromosomal microarray (CMA), multiplex ligation-dependent probe amplification (MLPA), or optical genome mapping.

Monosomy 1p36 (deletion of the distal 1.5–5.0 Mb of Chr1p36.33, affecting ~1 in 5,000 live births) removes dosage-sensitive subtelomeric genes including GABRD, KCNAB2, SKI, and PRDM16, resulting in hypotonia, structural heart defects, and refractory seizures. At Chr4q35, Facioscapulohumeral Muscular Dystrophy (FSHD1) arises not from a loss of coding sequence, but from contraction of the subtelomeric D4Z4 macrosatellite array from its normal 11–100 copies down to 1–10 copies on a permissive 4qA haplotype. This subtelomeric repeat contraction relaxes local CpG DNA methylation and H3K9me3 repression, de-repressing the retrogene DUX4 in the terminal D4Z4 unit and triggering myopathic transcriptional toxicity in skeletal muscle.

How do subtelomeres differ from telomeres?

Telomeres are the non-coding, strictly conserved (TTAGGG)_n hexanucleotide repeat arrays (typically 5–15 kb long in humans) bound by the six-protein shelterin complex at the physical ends of linear chromosomes. Subtelomeres are the 100–500 kb transitional regions located immediately inward (centromere-proximal) from the telomere; unlike uniform telomeric repeats, subtelomeres contain protein-coding genes, segmental duplications, non-coding TERRA promoters, and satellite arrays shared among subsets of chromosome arms.

Why are Non-B DNA structures (G-quadruplexes and Z-DNA) enriched in human subtelomeres?

Human subtelomeres undergo frequent double-strand breaks, GC-biased gene conversion during recombination repair, and tandem repeat expansion, driving regional GC content up to 55–65%. This combination of GC-richness, strand-asymmetric purine runs, and intronic (TG)_n/CpG microsatellites creates contiguous sequences that fold into potassium-stabilized G-quadruplexes and supercoil-driven left-handed Z-DNA, which in turn regulate replication fork progression, TERRA transcription, and chromatin looping.

Why do parasitic protozoa store surface antigen genes in subtelomeres?

Subtelomeres cluster together at the nuclear periphery in physical proximity to non-homologous chromosome ends. This spatial clustering allows double-strand breaks within subtelomeric repeats to be repaired via ectopic gene conversion and non-allelic homologous recombination with paralogous donor templates on other chromosomes—rapidly shuffling surface antigen sequences (var in Plasmodium falciparum, VSG in Trypanosoma brucei) without disrupting essential housekeeping genes in the chromosomal core.

1. Riethman, H., Ambrosini, A., Castaneda, C., Finklestein, J., Hu, X. L., Mudunuri, U., Paul, S., & Wei, J. (2004). "Mapping and initial analysis of human subtelomeric sequence assemblies." Genome Research, 14(1), 18–28. https://doi.org/10.1101/gr.1245004
2. Azzalin, C. M., Reichenbach, P., Khoriauli, L., Giulotto, E., & Lingner, J. (2007). "Telomeric repeat–containing RNA and RNA surveillance factors at mammalian chromosome ends." Science, 318(5851), 798–801. https://doi.org/10.1126/science.1147182
3. Robin, J. D., Ludlow, A. T., Batten, K., Magdinier, F., Stadler, G., Wagner, K. R., Shay, J. W., & Wright, W. E. (2014). "Telomere position effect: regulation of gene expression with progressive telomere shortening over long distances." Genes & Development, 28(22), 2464–2476. https://doi.org/10.1101/gad.251041.114
4. Nurk, S., Koren, S., Rhie, A., Rautiainen, M., Bzikadze, A. V., Mikheenko, A., Vollger, M. R., Altemose, N., Uralsky, L., Gershman, A., Aganezov, S., Hoyt, S. J., Diekhans, M., Logsdon, G. A., Alonge, M., Antonarakis, S. E., Borchers, M., Bouffard, G. G., Brooks, S. Y., ... Phillippy, A. M. (2022). "The complete sequence of a human genome." Science, 376(6588), 44–53. https://doi.org/10.1126/science.abj6987