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.

Genomic insights into metabolic flux in hummingbirds.

Hummingbirds are very well adapted to sustain efficient and rapid metabolic shifts. They oxidize ingested nectar to directly fuel flight when foraging but have to switch to oxidizing stored lipids derived from ingested sugars during the night or long-distance migratory flights. Understanding how this organism moderates energy turnover is hampered by a lack of information regarding how relevant enzymes differ in sequence, expression, and regulation. To explore these questions, we generated a chromosome-scale genome assembly of the ruby-throated hummingbird (A. colubris) using a combination of long- and short-read sequencing, scaffolding it using existing assemblies. We then used hybrid long- and short-read RNA sequencing of liver and muscle tissue in fasted and fed metabolic states for a comprehensive transcriptome assembly and annotation. Our genomic and transcriptomic data found positive selection of key metabolic genes in nectivorous avian species and deletion of critical genes (SLC2A4, GCK) involved in glucostasis in other vertebrates. We found expression of a fructose-specific version of SLC2A5 putatively in place of insulin-sensitive SLC2A5, with predicted protein models suggesting affinity for both fructose and glucose. Alternative isoforms may even act to sequester fructose to preclude limitations from transport in metabolism. Finally, we identified differentially expressed genes from fasted and fed hummingbirds, suggesting key pathways for the rapid metabolic switch hummingbirds undergo.

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.

Chromosome-specific telomere lengths and the minimal functional telomere revealed by nanopore sequencing.

We developed a method to tag telomeres and measure telomere length by nanopore sequencing in the yeast S. cerevisiae Nanopore allows long-read sequencing through the telomere, through the subtelomere, and into unique chromosomal sequence, enabling assignment of telomere length to a specific chromosome end. We observed chromosome end-specific telomere lengths that were stable over 120 cell divisions. These stable chromosome-specific telomere lengths may be explained by slow clonal variation or may represent a new biological mechanism that maintains equilibrium unique to each chromosome end. We examined the role of RIF1 and TEL1 in telomere length regulation and found that TEL1 is epistatic to RIF1 at most telomeres, consistent with the literature. However, at telomeres that lack subtelomeric Y' sequences, tel1Δ rif1Δ double mutants had a very small, but significant, increase in telomere length compared with the tel1Δ single mutant, suggesting an influence of Y' elements on telomere length regulation. We sequenced telomeres in a telomerase-null mutant (est2Δ) and found the minimal telomere length to be ∼75 bp. In these est2Δ mutants, there were apparent telomere recombination events at individual telomeres before the generation of survivors, and these events were significantly reduced in est2Δ rad52Δ double mutants. The rate of telomere shortening in the absence of telomerase was similar across all chromosome ends at ∼5 bp per generation. This new method gives quantitative, high-resolution telomere length measurement at each individual chromosome end and suggests possible new biological mechanisms regulating telomere length.

Simultaneous profiling of chromatin accessibility and methylation on human cell lines with nanopore sequencing.

Probing epigenetic features on DNA has tremendous potential to advance our understanding of the phased epigenome. In this study, we use nanopore sequencing to evaluate CpG methylation and chromatin accessibility simultaneously on long strands of DNA by applying GpC methyltransferase to exogenously label open chromatin. We performed nanopore sequencing of nucleosome occupancy and methylome (nanoNOMe) on four human cell lines (GM12878, MCF-10A, MCF-7 and MDA-MB-231). The single-molecule resolution allows footprinting of protein and nucleosome binding, and determination of the combinatorial promoter epigenetic signature on individual molecules. Long-read sequencing makes it possible to robustly assign reads to haplotypes, allowing us to generate a fully phased human epigenome, consisting of chromosome-level allele-specific profiles of CpG methylation and chromatin accessibility. We further apply this to a breast cancer model to evaluate differential methylation and accessibility between cancerous and noncancerous cells.

Clinical validation of a cell-free DNA fragmentome assay for augmentation of lung cancer early detection.

Lung cancer screening via annual low-dose computed tomography (LDCT) has poor adoption. We conducted a prospective case-control study among 958 individuals eligible for lung cancer screening to develop a blood-based lung cancer detection test that when positive is followed by an LDCT. Changes in genome-wide cell-free DNA (cfDNA) fragmentation profiles (fragmentomes) in peripheral blood reflected genomic and chromatin characteristics of lung cancer. We applied machine learning to fragmentome features to identify individuals who were more or less likely to have lung cancer. We trained the classifier using 576 cases and controls from study samples, and then validated it in a held-out group of 382 cases and controls. The validation demonstrated high sensitivity for lung cancer, and consistency across demographic groups and comorbid conditions. Applying test performance to the screening eligible population in a five-year model with modest utilization assumptions suggested the potential to prevent thousands of lung cancer deaths.

Epigenetic patterns in a complete human genome.

The completion of a telomere-to-telomere human reference genome, T2T-CHM13, has resolved complex regions of the genome, including repetitive and homologous regions. Here, we present a high-resolution epigenetic study of previously unresolved sequences, representing entire acrocentric chromosome short arms, gene family expansions, and a diverse collection of repeat classes. This resource precisely maps CpG methylation (32.28 million CpGs), DNA accessibility, and short-read datasets (166,058 previously unresolved chromatin immunoprecipitation sequencing peaks) to provide evidence of activity across previously unidentified or corrected genes and reveals clinically relevant paralog-specific regulation. Probing CpG methylation across human centromeres from six diverse individuals generated an estimate of variability in kinetochore localization. This analysis provides a framework with which to investigate the most elusive regions of the human genome, granting insights into epigenetic regulation.

Segmental duplications and their variation in a complete human genome.

Despite their importance in disease and evolution, highly identical segmental duplications (SDs) are among the last regions of the human reference genome (GRCh38) to be fully sequenced. Using a complete telomere-to-telomere human genome (T2T-CHM13), we present a comprehensive view of human SD organization. SDs account for nearly one-third of the additional sequence, increasing the genome-wide estimate from 5.4 to 7.0% [218 million base pairs (Mbp)]. An analysis of 268 human genomes shows that 91% of the previously unresolved T2T-CHM13 SD sequence (68.3 Mbp) better represents human copy number variation. Comparing long-read assemblies from human (n = 12) and nonhuman primate (n = 5) genomes, we systematically reconstruct the evolution and structural haplotype diversity of biomedically relevant and duplicated genes. This analysis reveals patterns of structural heterozygosity and evolutionary differences in SD organization between humans and other primates.

From telomere to telomere: The transcriptional and epigenetic state of human repeat elements.

Mobile elements and repetitive genomic regions are sources of lineage-specific genomic innovation and uniquely fingerprint individual genomes. Comprehensive analyses of such repeat elements, including those found in more complex regions of the genome, require a complete, linear genome assembly. We present a de novo repeat discovery and annotation of the T2T-CHM13 human reference genome. We identified previously unknown satellite arrays, expanded the catalog of variants and families for repeats and mobile elements, characterized classes of complex composite repeats, and located retroelement transduction events. We detected nascent transcription and delineated CpG methylation profiles to define the structure of transcriptionally active retroelements in humans, including those in centromeres. These data expand our insight into the diversity, distribution, and evolution of repetitive regions that have shaped the human genome.

Complete genomic and epigenetic maps of human centromeres.

Existing human genome assemblies have almost entirely excluded repetitive sequences within and near centromeres, limiting our understanding of their organization, evolution, and functions, which include facilitating proper chromosome segregation. Now, a complete, telomere-to-telomere human genome assembly (T2T-CHM13) has enabled us to comprehensively characterize pericentromeric and centromeric repeats, which constitute 6.2% of the genome (189.9 megabases). Detailed maps of these regions revealed multimegabase structural rearrangements, including in active centromeric repeat arrays. Analysis of centromere-associated sequences uncovered a strong relationship between the position of the centromere and the evolution of the surrounding DNA through layered repeat expansions. Furthermore, comparisons of chromosome X centromeres across a diverse panel of individuals illuminated high degrees of structural, epigenetic, and sequence variation in these complex and rapidly evolving regions.

The complete sequence of a human genome.

Since its initial release in 2000, the human reference genome has covered only the euchromatic fraction of the genome, leaving important heterochromatic regions unfinished. Addressing the remaining 8% of the genome, the Telomere-to-Telomere (T2T) Consortium presents a complete 3.055 billion-base pair sequence of a human genome, T2T-CHM13, that includes gapless assemblies for all chromosomes except Y, corrects errors in the prior references, and introduces nearly 200 million base pairs of sequence containing 1956 gene predictions, 99 of which are predicted to be protein coding. The completed regions include all centromeric satellite arrays, recent segmental duplications, and the short arms of all five acrocentric chromosomes, unlocking these complex regions of the genome to variational and functional studies.

De novo genome assembly of the tobacco hornworm moth (Manduca sexta).

The tobacco hornworm, Manduca sexta, is a lepidopteran insect that is used extensively as a model system for studying insect biology, development, neuroscience, and immunity. However, current studies rely on the highly fragmented reference genome Msex_1.0, which was created using now-outdated technologies and is hindered by a variety of deficiencies and inaccuracies. We present a new reference genome for M. sexta, JHU_Msex_v1.0, applying a combination of modern technologies in a de novo assembly to increase continuity, accuracy, and completeness. The assembly is 470 Mb and is ∼20× more continuous than the original assembly, with scaffold N50 > 14 Mb. We annotated the assembly by lifting over existing annotations and supplementing with additional supporting RNA-based data for a total of 25,256 genes. The new reference assembly is accessible in annotated form for public use. We demonstrate that improved continuity of the M. sexta genome improves resequencing studies and benefits future research on M. sexta as a model organism.

Epigenetically regulated digital signaling defines epithelial innate immunity at the tissue level.

To prevent damage to the host or its commensal microbiota, epithelial tissues must match the intensity of the immune response to the severity of a biological threat. Toll-like receptors allow epithelial cells to identify microbe associated molecular patterns. However, the mechanisms that mitigate biological noise in single cells to ensure quantitatively appropriate responses remain unclear. Here we address this question using single cell and single molecule approaches in mammary epithelial cells and primary organoids. We find that epithelial tissues respond to bacterial microbe associated molecular patterns by activating a subset of cells in an all-or-nothing (i.e. digital) manner. The maximum fraction of responsive cells is regulated by a bimodal epigenetic switch that licenses the TLR2 promoter for transcription across multiple generations. This mechanism confers a flexible memory of inflammatory events as well as unique spatio-temporal control of epithelial tissue-level immune responses. We propose that epigenetic licensing in individual cells allows for long-term, quantitative fine-tuning of population-level responses.

D-Serine Potently Drives Ligand-Binding Domain Closure in the Ionotropic Glutamate Receptor GluD2.

Despite their classification as ionotropic glutamate receptors, GluD receptors are not functional ligand-gated ion channels and do not bind glutamate. GluD2 receptors bind D-serine and coordinate transsynaptic complexes that regulate synaptic plasticity. Instead of opening the ion channel pore, mechanical tension produced from closure of GluD2 ligand-binding domains (LBDs) drives conformational rearrangements for non-ionotropic signaling. We report computed conformational free energy landscapes for the GluD2 LBD in apo and D-serine-bound states. Unexpectedly, the conformational free energy associated with GluD2 LBD closure upon D-serine binding is greater than that for AMPA, NMDA, and kainate receptor LBDs upon agonist binding. This excludes insufficient force generation as an explanation for lack of ion channel activity in GluD2 receptors and suggests that non-ionotropic conformational rearrangements do more work than pore opening. We also report free energy landscapes for GluD2 LBD harboring a neurodegenerative mutation and demonstrate selective stabilization of closed conformations in the apo state.

Investigating the potential use of an ionic liquid (1-Butyl-1-methylpyrrolidinium bis(trifluoromethylsulfonyl)imide) as an anti-fungal treatment against the amphibian chytrid fungus, Batrachochytrium dendrobatidis.

The disease chytridiomycosis, caused by the pathogenic chytrid fungus, Batrachochytrium dendrobatidis (Bd), has contributed to global amphibian declines. Bd infects the keratinized epidermal tissue in amphibians and causes hyperkeratosis and excessive skin shedding. In individuals of susceptible species, the regulatory function of the amphibian's skin is disrupted resulting in an electrolyte depletion, osmotic imbalance, and eventually death. Safe and effective treatments for chytridiomycosis are urgently needed to control chytrid fungal infections and stabilize populations of endangered amphibian species in captivity and in the wild. Currently, the most widely used anti-Bd treatment is itraconazole. Preparations of itraconazole formulated for amphibian use has proved effective, but treatment involves short baths over seven to ten days, a process which is logistically challenging, stressful, and causes long-term health effects. Here, we explore a novel anti-fungal therapeutic using a single application of the ionic liquid, 1-Butyl-1-methylpyrrolidinium bis(trifluoromethylsulfonyl)imide (BMP-NTf2), for the treatment of chytridiomycosis. BMP-NTf2 was found be effective at killing Bd in vitro at low concentrations (1:1000 dilution). We tested BMP-NTf2 in vivo on two amphibian species, one that is relatively tolerant of chytridiomycosis (Pseudacris regilla) and one that is highly susceptible (Dendrobates tinctorius). A toxicity trial revealed a surprising interaction between Bd infection status and the impact of BMP-NTf2 on D. tinctorius survival. Uninfected D. tinctorius tolerated BMP-NTf2 (mean ± SE; 96.01 ± 9.00 µl/g), such that only 1 out of 30 frogs died following treatment (at a dose of 156.95 µL/g), whereas, a lower dose (mean ± SE; 97.45 ± 3.52 µL/g) was not tolerated by Bd-infected D. tinctorius, where 15 of 23 frogs died shortly upon BMP-NTf2 application. Those that tolerated the BMP-NTf2 application did not exhibit Bd clearance. Thus, BMP-NTf2 application, under the conditions tested here, is not a suitable option for clearing Bd infection in D. tinctorius. However, different results were obtained for P. regilla. Two topical applications of BMP-NTf2 on Bd-infected P. regilla (using a lower BMP-NTf2 dose than on D. tinctorius, mean ± SE; 9.42 ± 1.43 µL/g) reduced Bd growth, although the effect was lower than that obtained by daily doses of itracanozole (50% frogs exhibited complete clearance on day 16 vs. 100% for itracanozole). Our findings suggest that BMP-NTf2 has the potential to treat Bd infection, however the effect depends on several parameters. Further optimization of dose and schedule are needed before BMP-NTf2 can be considered as a safe and effective alternative to more conventional antifungal agents, such as itraconazole.