Human Genome & T2T Milestone

Human Genome & T2T Milestone

For over two decades, the "human reference genome" served as the foundational coordinate system for virtually all genomic analysis — from variant calling to transcript quantification. However, the reference assemblies released by the Human Genome Project and refined through GRCh37 and GRCh38 were fundamentally incomplete. Roughly 8% of the genome, concentrated in highly repetitive regions such as centromeres, satellite arrays, and segmental duplications, could not be assembled or mapped reliably with short-read sequencing technology. These "gaps" were not biologically irrelevant: they contain genes, regulatory elements, and structural variation relevant to health and disease. The Telomere-to-Telomere (T2T) Consortium achieved a landmark in 2022 by publishing the first truly complete human genome assembly, T2T-CHM13. This accomplishment was enabled by long-read sequencing platforms — PacBio HiFi reads and Oxford Nanopore — combined with ultra-long reads and high-fidelity assembly algorithms such as hifiasm. The CHM13 cell line, derived from a complete hydatidiform mole, was strategically chosen because it is essentially homozygous across the genome, dramatically simplifying assembly of a single haploid representation. The T2T assembly added nearly 200 million base pairs of previously unresolved sequence, including fully assembled centromeres and the short arms of acrocentric chromosomes. For bioinformaticians, the T2T milestone has practical consequences. First, reference choice matters: analyses performed against T2T-CHM13 can recover reads that GRCh38 would discard or mis-map, particularly in repetitive and duplicated regions. Second, new annotation resources (such as GENCODE mappings onto T2T) must be used consistently with the reference build — mixing coordinates from different assemblies is a common and serious error. Third, pangenome efforts like the Human Pangenome Reference Consortium build directly on T2T-quality assemblies to represent human diversity more completely than any single reference can. The following Python example demonstrates how to compute basic sequence statistics from a FASTA file, a routine sanity check when working with any reference assembly, including T2T-CHM13: from collections import Counter def fasta_stats(path): """Compute GC content and N-gap statistics for a genome assembly.""" seq = [] with open(path) as fh: for line in fh: if not line.startswith(">"): seq.append(line.strip().upper()) genome = "".join(seq) counts = Counter(genome) total = sum(counts.values()) gc = (counts.get("G", 0) + counts.get("C", 0)) / total * 100 gaps = counts.get("N", 0) print(f"Total length: {total:,} bp") print(f"GC content: {gc:.2f}%") print(f"N bases: {gaps:,} ({gaps/total*100:.4f}%)") fasta_stats("t2t_chm13.fa") A defining feature of T2T-CHM13 is that it contains essentially no unresolved "N" gaps — in contrast to GRCh38, which still carries tens of millions of N characters. The next example shows how to detect and report contiguous gap regions in an assembly, a useful exercise when comparing reference builds: def find_gaps(genome, min_length=10000): """Yield (start, end, length) for gap runs of Ns >= min_length.""" start = None for i, base in enumerate(genome): if base == "N": if start is None: start = i elif start is not None: if i - start >= min_length: yield (start, i, i - start) start = None if start is not None and len(genome) - start >= min_length: yield (start, len(genome), len(genome) - start) # Usage: gaps = list(find_gaps(genome_string)) # GRCh38: hundreds of large gaps; T2T-CHM13: none expected. The T2T-CHM13 assembly (2022) is the first complete, gapless human genome assembly, resolving centromeres, satellite arrays, and segmental duplications that were missing from GRCh38. Long-read sequencing (PacBio HiFi, Oxford Nanopore) combined with modern assemblers like hifiasm made the T2T milestone possible. The CHM13 cell line's near-homozygosity simplified assembly by reducing haplotype complexity to a single haploid representation. Reference build consistency is critical: coordinates, annotations, and variant calls must all correspond to the same assembly (GRCh38 vs. T2T-CHM13). T2T-quality assemblies underpin the Human Pangenome Reference Consortium's effort to build a more inclusive, multi-haplotype reference. Practice exercise: Download the chromosome 1 sequences for both GRCh38 and T2T-CHM13 from the NCBI or Ensembl FTP servers. Using the find_gaps function above, identify and count all gap regions of at least 10 kb in each assembly, and compare their total length. Then compute GC content for the centromeric region (approximately positions 121–128 Mb in T2T-CHM13 chromosome 1) and contrast it with GC content in a gene-rich region of your choice. Write a short summary (150–200 words) explaining how your results illustrate why the T2T assembly matters for studies of repetitive DNA and structural variation.