The variation and evolution of complete human centromeres.

Human centromeres have been traditionally very difficult to sequence and assemble owing to their repetitive nature and large size1. As a result, patterns of human centromeric variation and models for their evolution and function remain incomplete, despite centromeres being among the most rapidly mutating regions2,3. Here, using long-read sequencing, we completely sequenced and assembled all centromeres from a second human genome and compared it to the finished reference genome4,5. We find that the two sets of centromeres show at least a 4.1-fold increase in single-nucleotide variation when compared with their unique flanks and vary up to 3-fold in size. Moreover, we find that 45.8% of centromeric sequence cannot be reliably aligned using standard methods owing to the emergence of new α-satellite higher-order repeats (HORs). DNA methylation and CENP-A chromatin immunoprecipitation experiments show that 26% of the centromeres differ in their kinetochore position by >500 kb. To understand evolutionary change, we selected six chromosomes and sequenced and assembled 31 orthologous centromeres from the common chimpanzee, orangutan and macaque genomes. Comparative analyses reveal a nearly complete turnover of α-satellite HORs, with characteristic idiosyncratic changes in α-satellite HORs for each species. Phylogenetic reconstruction of human haplotypes supports limited to no recombination between the short (p) and long (q) arms across centromeres and reveals that novel α-satellite HORs share a monophyletic origin, providing a strategy to estimate the rate of saltatory amplification and mutation of human centromeric DNA.

Recombination between heterologous human acrocentric chromosomes.

The short arms of the human acrocentric chromosomes 13, 14, 15, 21 and 22 (SAACs) share large homologous regions, including ribosomal DNA repeats and extended segmental duplications1,2. Although the resolution of these regions in the first complete assembly of a human genome-the Telomere-to-Telomere Consortium's CHM13 assembly (T2T-CHM13)-provided a model of their homology3, it remained unclear whether these patterns were ancestral or maintained by ongoing recombination exchange. Here we show that acrocentric chromosomes contain pseudo-homologous regions (PHRs) indicative of recombination between non-homologous sequences. Utilizing an all-to-all comparison of the human pangenome from the Human Pangenome Reference Consortium4 (HPRC), we find that contigs from all of the SAACs form a community. A variation graph5 constructed from centromere-spanning acrocentric contigs indicates the presence of regions in which most contigs appear nearly identical between heterologous acrocentric chromosomes in T2T-CHM13. Except on chromosome 15, we observe faster decay of linkage disequilibrium in the pseudo-homologous regions than in the corresponding short and long arms, indicating higher rates of recombination6,7. The pseudo-homologous regions include sequences that have previously been shown to lie at the breakpoint of Robertsonian translocations8, and their arrangement is compatible with crossover in inverted duplications on chromosomes 13, 14 and 21. The ubiquity of signals of recombination between heterologous acrocentric chromosomes seen in the HPRC draft pangenome suggests that these shared sequences form the basis for recurrent Robertsonian translocations, providing sequence and population-based confirmation of hypotheses first developed from cytogenetic studies 50 years ago9.

The complete sequence of a human Y chromosome.

The human Y chromosome has been notoriously difficult to sequence and assemble because of its complex repeat structure that includes long palindromes, tandem repeats and segmental duplications1-3. As a result, more than half of the Y chromosome is missing from the GRCh38 reference sequence and it remains the last human chromosome to be finished4,5. Here, the Telomere-to-Telomere (T2T) consortium presents the complete 62,460,029-base-pair sequence of a human Y chromosome from the HG002 genome (T2T-Y) that corrects multiple errors in GRCh38-Y and adds over 30 million base pairs of sequence to the reference, showing the complete ampliconic structures of gene families TSPY, DAZ and RBMY; 41 additional protein-coding genes, mostly from the TSPY family; and an alternating pattern of human satellite 1 and 3 blocks in the heterochromatic Yq12 region. We have combined T2T-Y with a previous assembly of the CHM13 genome4 and mapped available population variation, clinical variants and functional genomics data to produce a complete and comprehensive reference sequence for all 24 human chromosomes.

Semi-automated assembly of high-quality diploid human reference genomes.

The current human reference genome, GRCh38, represents over 20 years of effort to generate a high-quality assembly, which has benefitted society1,2. However, it still has many gaps and errors, and does not represent a biological genome as it is a blend of multiple individuals3,4. Recently, a high-quality telomere-to-telomere reference, CHM13, was generated with the latest long-read technologies, but it was derived from a hydatidiform mole cell line with a nearly homozygous genome5. To address these limitations, the Human Pangenome Reference Consortium formed with the goal of creating high-quality, cost-effective, diploid genome assemblies for a pangenome reference that represents human genetic diversity6. Here, in our first scientific report, we determined which combination of current genome sequencing and assembly approaches yield the most complete and accurate diploid genome assembly with minimal manual curation. Approaches that used highly accurate long reads and parent-child data with graph-based haplotype phasing during assembly outperformed those that did not. Developing a combination of the top-performing methods, we generated our first high-quality diploid reference assembly, containing only approximately four gaps per chromosome on average, with most chromosomes within ±1% of the length of CHM13. Nearly 48% of protein-coding genes have non-synonymous amino acid changes between haplotypes, and centromeric regions showed the highest diversity. Our findings serve as a foundation for assembling near-complete diploid human genomes at scale for a pangenome reference to capture global genetic variation from single nucleotides to structural rearrangements.

Cell-cycle-coupled structural oscillation of centromeric nucleosomes in yeast.

The centromere is a specialized chromosomal structure that regulates chromosome segregation. Centromeres are marked by a histone H3 variant. In budding yeast, the histone H3 variant Cse4 is present in a single centromeric nucleosome. Experimental evidence supports several different models for the structure of centromeric nucleosomes. To investigate Cse4 copy number in live yeast, we developed a method coupling fluorescence correlation spectroscopy and calibrated imaging. We find that centromeric nucleosomes have one copy of Cse4 during most of the cell cycle, whereas two copies are detected at anaphase. The proposal of an anaphase-coupled structural change is supported by Cse4-Cse4 interactions, incorporation of Cse4, and the absence of Scm3 in anaphase. Nucleosome reconstitution and ChIP suggests both Cse4 structures contain H2A/H2B. The increase in Cse4 intensity and deposition at anaphase are also observed in Candida albicans. Our experimental evidence supports a cell-cycle-coupled oscillation of centromeric nucleosome structure in yeast.

The structure, function and evolution of a complete human chromosome 8.

The complete assembly of each human chromosome is essential for understanding human biology and evolution1,2. Here we use complementary long-read sequencing technologies to complete the linear assembly of human chromosome 8. Our assembly resolves the sequence of five previously long-standing gaps, including a 2.08-Mb centromeric α-satellite array, a 644-kb copy number polymorphism in the ß-defensin gene cluster that is important for disease risk, and an 863-kb variable number tandem repeat at chromosome 8q21.2 that can function as a neocentromere. We show that the centromeric α-satellite array is generally methylated except for a 73-kb hypomethylated region of diverse higher-order α-satellites enriched with CENP-A nucleosomes, consistent with the location of the kinetochore. In addition, we confirm the overall organization and methylation pattern of the centromere in a diploid human genome. Using a dual long-read sequencing approach, we complete high-quality draft assemblies of the orthologous centromere from chromosome 8 in chimpanzee, orangutan and macaque to reconstruct its evolutionary history. Comparative and phylogenetic analyses show that the higher-order α-satellite structure evolved in the great ape ancestor with a layered symmetry, in which more ancient higher-order repeats locate peripherally to monomeric α-satellites. We estimate that the mutation rate of centromeric satellite DNA is accelerated by more than 2.2-fold compared to the unique portions of the genome, and this acceleration extends into the flanking sequence.

Regulators of rDNA array morphology in fission yeast.

Nucleolar morphology is a well-established indicator of ribosome biogenesis activity that has served as the foundation of many screens investigating ribosome production. Missing from this field of study is a broad-scale investigation of the regulation of ribosomal DNA morphology, despite the essential role of rRNA gene transcription in modulating ribosome output. We hypothesized that the morphology of rDNA arrays reflects ribosome biogenesis activity. We established GapR-GFP, a prokaryotic DNA-binding protein that recognizes transcriptionally-induced overtwisted DNA, as a live visual fluorescent marker for quantitative analysis of rDNA organization in Schizosaccharomyces pombe. We found that the morphology-which we refer to as spatial organization-of the rDNA arrays is dynamic throughout the cell cycle, under glucose starvation, RNA pol I inhibition, and TOR activation. Screening the haploid S. pombe Bioneer deletion collection for spatial organization phenotypes revealed large ribosomal protein (RPL) gene deletions that alter rDNA organization. Further work revealed RPL gene deletion mutants with altered rDNA organization also demonstrate resistance to the TOR inhibitor Torin1. A genetic analysis of signaling pathways essential for this resistance phenotype implicated many factors including a conserved MAPK, Pmk1, previously linked to extracellular stress responses. We propose RPL gene deletion triggers altered rDNA morphology due to compensatory changes in ribosome biogenesis via multiple signaling pathways, and we further suggest compensatory responses may contribute to human diseases such as ribosomopathies. Altogether, GapR-GFP is a powerful tool for live visual reporting on rDNA morphology under myriad conditions.

A working model for the formation of Robertsonian chromosomes.

Robertsonian chromosomes form by fusion of two chromosomes that have centromeres located near their ends, known as acrocentric or telocentric chromosomes. This fusion creates a new metacentric chromosome and is a major mechanism of karyotype evolution and speciation. Robertsonian chromosomes are common in nature and were first described in grasshoppers by the zoologist W. R. B. Robertson more than 100 years ago. They have since been observed in many species, including catfish, sheep, butterflies, bats, bovids, rodents and humans, and are the most common chromosomal change in mammals. Robertsonian translocations are particularly rampant in the house mouse, Mus musculus domesticus, where they exhibit meiotic drive and create reproductive isolation. Recent progress has been made in understanding how Robertsonian chromosomes form in the human genome, highlighting some of the fundamental principles of how and why these types of fusion events occur so frequently. Consequences of these fusions include infertility and Down's syndrome. In this Hypothesis, I postulate that the conditions that allow these fusions to form are threefold: (1) sequence homology on non-homologous chromosomes, often in the form of repetitive DNA; (2) recombination initiation during meiosis; and (3) physical proximity of the homologous sequences in three-dimensional space. This Hypothesis highlights the latest progress in understanding human Robertsonian translocations within the context of the broader literature on Robertsonian chromosomes.

Myc promotes polyploidy in murine trophoblast cells and suppresses senescence.

The placenta is essential for reproductive success. The murine placenta includes polyploid giant cells that are crucial for its function. Polyploidy occurs broadly in nature but its regulators and significance in the placenta are unknown. We have discovered that many murine placental cell types are polyploid and have identified factors that license polyploidy using single-cell RNA sequencing. Myc is a key regulator of polyploidy and placental development, and is required for multiple rounds of DNA replication, likely via endocycles, in trophoblast giant cells. Furthermore, MYC supports the expression of DNA replication and nucleotide biosynthesis genes along with ribosomal RNA. Increased DNA damage and senescence occur in trophoblast giant cells without Myc, accompanied by senescence in the neighboring maternal decidua. These data reveal Myc is essential for polyploidy to support normal placental development, thereby preventing premature senescence. Our study, combined with available literature, suggests that Myc is an evolutionarily conserved regulator of polyploidy.

Telomere-to-telomere assembly of a complete human X chromosome.

After two decades of improvements, the current human reference genome (GRCh38) is the most accurate and complete vertebrate genome ever produced. However, no single chromosome has been finished end to end, and hundreds of unresolved gaps persist1,2. Here we present a human genome assembly that surpasses the continuity of GRCh382, along with a gapless, telomere-to-telomere assembly of a human chromosome. This was enabled by high-coverage, ultra-long-read nanopore sequencing of the complete hydatidiform mole CHM13 genome, combined with complementary technologies for quality improvement and validation. Focusing our efforts on the human X chromosome3, we reconstructed the centromeric satellite DNA array (approximately 3.1 Mb) and closed the 29 remaining gaps in the current reference, including new sequences from the human pseudoautosomal regions and from cancer-testis ampliconic gene families (CT-X and GAGE). These sequences will be integrated into future human reference genome releases. In addition, the complete chromosome X, combined with the ultra-long nanopore data, allowed us to map methylation patterns across complex tandem repeats and satellite arrays. Our results demonstrate that finishing the entire human genome is now within reach, and the data presented here will facilitate ongoing efforts to complete the other human chromosomes.

Regulators of the cohesin network.

Chromosome cohesion is the term used to describe the cellular process in which sister chromatids are held together from the time of their replication until the time of separation at the metaphase to anaphase transition. In this capacity, chromosome cohesion, especially at centromeric regions, is essential for chromosome segregation. However, cohesion of noncentromeric DNA sequences has been shown to occur during double-strand break (DSB) repair and the transcriptional regulation of genes. Cohesion for the purposes of accurate chromosome segregation, DSB repair, and gene regulation are all achieved through a similar network of proteins, but cohesion for each purpose may be regulated differently. In this review, we focus on recent developments regarding the regulation of this multipurpose network for tying DNA sequences together. In particular, regulation via effectors and posttranslational modifications are reviewed. A picture is emerging in which complex regulatory networks are capable of differential regulation of cohesion in various contexts.

Maternal Smc3 protects the integrity of the zygotic genome through DNA replication and mitosis.

Aneuploidy is frequently observed in oocytes and early embryos, begging the question of how genome integrity is monitored and preserved during this crucial period. SMC3 is a subunit of the cohesin complex that supports genome integrity, but its role in maintaining the genome during this window of mammalian development is unknown. We discovered that, although depletion of Smc3 following meiotic S phase in mouse oocytes allowed accurate meiotic chromosome segregation, adult females were infertile. We provide evidence that DNA lesions accumulated following S phase in SMC3-deficient zygotes, followed by mitosis with lagging chromosomes, elongated spindles, micronuclei, and arrest at the two-cell stage. Remarkably, although centromeric cohesion was defective, the dosage of SMC3 was sufficient to enable embryogenesis in juvenile mutant females. Our findings suggest that, despite previous reports of aneuploidy in early embryos, chromosome missegregation in zygotes halts embryogenesis at the two-cell stage. Smc3 is a maternal gene with essential functions in the repair of spontaneous damage associated with DNA replication and subsequent chromosome segregation in zygotes, making cohesin a key protector of the zygotic genome.

Integrative structure and functional anatomy of a nuclear pore complex.

Nuclear pore complexes play central roles as gatekeepers of RNA and protein transport between the cytoplasm and nucleoplasm. However, their large size and dynamic nature have impeded a full structural and functional elucidation. Here we determined the structure of the entire 552-protein nuclear pore complex of the yeast Saccharomyces cerevisiae at sub-nanometre precision by satisfying a wide range of data relating to the molecular arrangement of its constituents. The nuclear pore complex incorporates sturdy diagonal columns and connector cables attached to these columns, imbuing the structure with strength and flexibility. These cables also tie together all other elements of the nuclear pore complex, including membrane-interacting regions, outer rings and RNA-processing platforms. Inwardly directed anchors create a high density of transport factor-docking Phe-Gly repeats in the central channel, organized into distinct functional units. This integrative structure enables us to rationalize the architecture, transport mechanism and evolutionary origins of the nuclear pore complex.

Taking cohesin and condensin in context.

Structural maintenance of chromosome (SMC) protein complexes, including cohesin and condensin, are increasingly being recognized for their important role in cancer and development, making it critical that we understand how these evolutionarily conserved multi-subunit protein complexes associate with and organize the genome. We review adaptor proteins for SMC complexes and how these adaptors may capture SMC complexes following loop extrusion to provide a framework for chromosome organization.

Ribosomal DNA instability and genome adaptability.

Ribosomes are large, multi-subunit ribonucleoprotein complexes, essential for protein synthesis. To meet the high cellular demand for ribosomes, all eukaryotes have numerous copies of ribosomal DNA (rDNA) genes that encode ribosomal RNA (rRNA), usually far in excess of the requirement for ribosome biogenesis. In all eukaryotes studied, rDNA genes are arranged in one or more clusters of tandem repeats localized to nucleoli. The tandem arrangement of repeats, combined with the high rates of transcription at the rDNA loci, and the difficulty of replicating repetitive sequences make the rDNA inherently unstable and particularly susceptible to large variations in repeat copy number. Despite mounting evidence suggesting extra-ribosomal functions of the rDNA, its repetitive nature has excluded it from traditional sequencing-based studies. However, more recently, several studies have revealed the unique potential of the rDNA to act as a "canary in the coalmine," being particularly sensitive to genomic stresses and acting as a source of adaptive response. Here, we review evidence uncovering mechanisms of regulation of instability and copy number variation at the rDNA and their role in adaptation to the environment, which could serve to understand the basic principles governing the behavior of other tandem repeats and their role in shaping the genome.

Ribosomal DNA-connecting ribosome biogenesis and chromosome biology.

Ribosomal DNA, the topic of this special issue, has long fascinated biologists. The RNA products of the ribosomal DNA are the ribosomal RNAs that are part of the ribosome. In this special issue, we focus on the sequence, molecular organization, repair, stability, copy number, and peculiar genetics of this region of the genome. The locus can impact not only the translational capability of cells, but also genome organization, stability and integrity, providing a link between translation and chromosome biology.

Ribosomal DNA and the nucleolus in the context of genome organization.

The nucleolus constitutes a prominent nuclear compartment, a membraneless organelle that was first documented in the 1830s. The fact that specific chromosomal regions were present in the nucleolus was recognized by Barbara McClintock in the 1930s, and these regions were termed nucleolar organizing regions, or NORs. The primary function of ribosomal DNA (rDNA) is to produce RNA components of ribosomes. Yet, ribosomal DNA also plays a pivotal role in nuclear organization by assembling the nucleolus. This review is focused on the rDNA and associated proteins in the context of genome organization. Recent advances in understanding chromatin organization suggest that chromosomes are organized into topological domains by a DNA loop extrusion process. We discuss the perspective that rDNA may also be organized in topological domains constrained by structural maintenance of chromosome protein complexes such as cohesin and condensin. Moreover, biophysical studies indicate that the nucleolar compartment may be formed by active processes as well as phase separation, a perspective that lends further insight into nucleolar organization. The application of the latest perspectives and technologies to this organelle help further elucidate its role in nuclear structure and function.

Chromatin assembly factor-1 (CAF-1) chaperone regulates Cse4 deposition into chromatin in budding yeast.

Correct localization of the centromeric histone variant CenH3/CENP-A/Cse4 is an important part of faithful chromosome segregation. Mislocalization of CenH3 could affect chromosome segregation, DNA replication and transcription. CENP-A is often overexpressed and mislocalized in cancer genomes, but the underlying mechanisms are not understood. One major regulator of Cse4 deposition is Psh1, an E3 ubiquitin ligase that controls levels of Cse4 to prevent deposition into non-centromeric regions. We present evidence that Chromatin assembly factor-1 (CAF-1), an evolutionarily conserved histone H3/H4 chaperone with subunits shown previously to interact with CenH3 in flies and human cells, regulates Cse4 deposition in budding yeast. yCAF-1 interacts with Cse4 and can assemble Cse4 nucleosomes in vitro. Loss of yCAF-1 dramatically reduces the amount of Cse4 deposited into chromatin genome-wide when Cse4 is overexpressed. The incorporation of Cse4 genome-wide may have multifactorial effects on growth and gene expression. Loss of yCAF-1 can rescue growth defects and some changes in gene expression associated with Cse4 deposition that occur in the absence of Psh1-mediated proteolysis. Incorporation of Cse4 into promoter nucleosomes at transcriptionally active genes depends on yCAF-1. Overall our findings suggest CAF-1 can act as a CenH3 chaperone, regulating levels and incorporation of CenH3 in chromatin.

DNA replication stress restricts ribosomal DNA copy number.

Ribosomal RNAs (rRNAs) in budding yeast are encoded by ~100-200 repeats of a 9.1kb sequence arranged in tandem on chromosome XII, the ribosomal DNA (rDNA) locus. Copy number of rDNA repeat units in eukaryotic cells is maintained far in excess of the requirement for ribosome biogenesis. Despite the importance of the repeats for both ribosomal and non-ribosomal functions, it is currently not known how "normal" copy number is determined or maintained. To identify essential genes involved in the maintenance of rDNA copy number, we developed a droplet digital PCR based assay to measure rDNA copy number in yeast and used it to screen a yeast conditional temperature-sensitive mutant collection of essential genes. Our screen revealed that low rDNA copy number is associated with compromised DNA replication. Further, subculturing yeast under two separate conditions of DNA replication stress selected for a contraction of the rDNA array independent of the replication fork blocking protein, Fob1. Interestingly, cells with a contracted array grew better than their counterparts with normal copy number under conditions of DNA replication stress. Our data indicate that DNA replication stresses select for a smaller rDNA array. We speculate that this liberates scarce replication factors for use by the rest of the genome, which in turn helps cells complete DNA replication and continue to propagate. Interestingly, tumors from mini chromosome maintenance 2 (MCM2)-deficient mice also show a loss of rDNA repeats. Our data suggest that a reduction in rDNA copy number may indicate a history of DNA replication stress, and that rDNA array size could serve as a diagnostic marker for replication stress. Taken together, these data begin to suggest the selective pressures that combine to yield a "normal" rDNA copy number.

Ribosomal DNA copy number loss and sequence variation in cancer.

Ribosomal DNA is one of the most variable regions in the human genome with respect to copy number. Despite the importance of rDNA for cellular function, we know virtually nothing about what governs its copy number, stability, and sequence in the mammalian genome due to challenges associated with mapping and analysis. We applied computational and droplet digital PCR approaches to measure rDNA copy number in normal and cancer states in human and mouse genomes. We find that copy number and sequence can change in cancer genomes. Counterintuitively, human cancer genomes show a loss of copies, accompanied by global copy number co-variation. The sequence can also be more variable in the cancer genome. Cancer genomes with lower copies have mutational evidence of mTOR hyperactivity. The PTEN phosphatase is a tumor suppressor that is critical for genome stability and a negative regulator of the mTOR kinase pathway. Surprisingly, but consistent with the human cancer genomes, hematopoietic cancer stem cells from a Pten-/- mouse model for leukemia have lower rDNA copy number than normal tissue, despite increased proliferation, rRNA production, and protein synthesis. Loss of copies occurs early and is associated with hypersensitivity to DNA damage. Therefore, copy loss is a recurrent feature in cancers associated with mTOR activation. Ribosomal DNA copy number may be a simple and useful indicator of whether a cancer will be sensitive to DNA damaging treatments.