Search

Here the Human Pangenome Reference Consortium presents a first draft of the human pangenome reference. The pangenome contains 47 phased, diploid assemblies from a cohort of genetically diverse individuals1. These assemblies cover more than 99% of the expected sequence in each genome and are more than 99% accurate at the structural and base pair levels. Based on alignments of the assemblies, we generate a draft pangenome that captures known variants and haplotypes and reveals new alleles at structurally complex loci. We also add 119 million base pairs of euchromatic polymorphic sequences and 1,115 gene duplications relative to the existing reference GRCh38. Roughly 90 million of the additional base pairs are derived from structural variation. Using our draft pangenome to analyse short-read data reduced small variant discovery errors by 34% and increased the number of structural variants detected per haplotype by 104% compared with GRCh38-based workflows, which enabled the typing of the vast majority of structural variant alleles per sample.

Most large-scale genetic studies of autism have focused on the discovery of genes by proving an enrichment of de novo mutations (DNMs) in autism probands or characterizing polygenic risk based on the association of common variants. We present evidence in support of an oligogenic model where two or more ultrarare mutations of more modest effect are preferentially transmitted to children with autism. Such private gene-disruptive mutations are enriched in families where there are multiple affected individuals, emerged two or three generations ago, and map to genes not previously associated with autism. Although no single gene has reached statistical significance, this class of variation should be considered along with genetic and nongenetic factors to better explain the etiology of this complex trait.

Purpose: Genome-wide sequencing is increasingly being performed during pregnancy to identify the genetic cause of congenital anomalies. The interpretation of prenatally identified variants can be challenging and is hampered by our often limited knowledge of prenatal phenotypes. To better delineate the prenatal phenotype of Coffin-Siris syndrome (CSS), we collected clinical data from patients with a prenatal phenotype and a pathogenic variant in one of the CSS-associated genes. Methods: Clinical data was collected through an extensive web-based survey. Results: We included 44 patients with a variant in a CSS-associated gene and a prenatal phenotype; 9 of these patients have been reported before. Prenatal anomalies that were frequently observed in our cohort include hydrocephalus, agenesis of the corpus callosum, hypoplastic left heart syndrome, persistent left vena cava, diaphragmatic hernia, renal agenesis, and intrauterine growth restriction. Anal anomalies were frequently identified after birth in patients with ARID1A variants (6/14, 43%). Interestingly, pathogenic ARID1A variants were much more frequently identified in the current prenatal cohort (16/44, 36%) than in postnatal CSS cohorts (5%-9%). Conclusion: Our data shed new light on the prenatal phenotype of patients with pathogenic variants in CSS genes. (C) 2022 American College of Medical Genetics and Genomics. Published by Elsevier Inc. All rights reserved.

A critical step in preserving protein homeostasis is the recognition, binding, unfolding, and translocation of protein substrates by six AAA-ATPase proteasome subunits (ATPase-associated with various cellular activities) termed PSMC1-6, which are required for degradation of proteins by 26 S proteasomes. Here, we identified 15 de novo missense variants in the PSMC3 gene encoding the AAA-ATPase proteasome subunit PSMC3/Rpt5 in 23 unrelated heterozygous patients with an autosomal dominant form of neurodevelopmental delay and intellectual disability. Expression of PSMC3 variants in mouse neuronal cultures led to altered dendrite development, and deletion of the PSMC3 fly ortholog Rpt5 impaired reversal learning capabilities in fruit flies. Structural modeling as well as proteomic and transcriptomic analyses of T cells derived from patients with PSMC3 variants implicated the PSMC3 variants in proteasome dysfunction through disruption of substrate translocation, induction of proteotoxic stress, and alterations in proteins controlling developmental and innate immune programs. The proteostatic perturbations in T cells from patients with PSMC3 variants correlated with a dysregulation in type I interferon (IFN) signaling in these T cells, which could be blocked by inhibition of the intracellular stress sensor protein kinase R (PKR). These results suggest that proteotoxic stress activated PKR in patient-derived T cells, resulting in a type I IFN response. The potential relationship among proteosome dysfunction, type I IFN production, and neurodevelopment suggests new directions in our understanding of pathogenesis in some neurodevelopmental disorders.

We completely sequenced and assembled all centromeres from a second human genome and used two reference sets to benchmark genetic, epigenetic, and evolutionary variation within centromeres from a diversity panel of humans and apes. We find that centromere single-nucleotide variation can increase by up to 4.1-fold relative to other genomic regions, with the caveat that up to 45.8% of centromeric sequence, on average, cannot be reliably aligned with current methods due to the emergence of new alpha-satellite higher-order repeat (HOR) structures and two to threefold differences in the length of the centromeres. The extent to which this occurs differs depending on the chromosome and haplotype. Comparing the two sets of complete human centromeres, we find that eight harbor distinctly different alpha-satellite HOR array structures and four contain novel alpha-satellite HOR variants in high abundance. DNA methylation and CENP-A chromatin immunoprecipitation experiments show that 26% of the centromeres differ in their kinetochore position by at least 500 kbp - a property not readily associated with novel alpha-satellite HORs. 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 nearly complete turnover of alpha-satellite HORs, but with idiosyncratic changes in structure characteristic to each species. Phylogenetic reconstruction of human haplotypes supports limited to no recombination between the p- and q-arms of human chromosomes and reveals that novel alpha-satellite HORs share a monophyletic origin, providing a strategy to estimate the rate of saltatory amplification and mutation of human centromeric DNA.

Advances in long-read sequencing (LRS) technology continue to make whole-genome sequencing more complete, affordable, and accurate. LRS provides significant advantages over short-read sequencing approaches, including phasedde novogenome assembly, access to previously excluded genomic regions, and discovery of more complex structural variants (SVs) associated with disease. Limitations remain with respect to cost, scalability, and platform-dependent read accuracy and the tradeoffs between sequence coverage and sensitivity of variant discovery are important experimental considerations for the application of LRS. We compare the genetic variant calling precision and recall of Oxford Nanopore Technologies (ONT) and PacBio HiFi platforms over a range of sequence coverages. For read-based applications, LRS sensitivity begins to plateau around 12-fold coverage with a majority of variants called with reasonable accuracy (F1 score above 0.5), and both platforms perform well for SV detection. Genome assembly increases variant calling precision and recall of SVs and indels in HiFi datasets with HiFi outperforming ONT in quality as measured by the F1 score of assembly-based variant callsets. While both technologies continue to evolve, our work offers guidance to design cost-effective experimental strategies that do not compromise on discovering novel biology.

SummaryThe focal attachment of the kinetochore to the centromere core is essential for genome maintenance, yet the highly repetitive nature of human centromeres limits our understanding of their chromatin organization. We demonstrate that single-molecule chromatin fiber sequencing can uniquely resolve chromatin organization within centromeres at single-molecule and single-nucleotide resolution. We find that the centromere core contains a dichotomous chromatin organization not found elsewhere in the genome, which is characterized by highly accessible chromatin patches heterogeneously punctuated amongst tightly compacted nucleosome arrays. These highly accessible chromatin patches correspond to sites of kinetochore attachment, and clustered CENP-B occupancy within these patches phase nucleosome arrays to the alpha-satellite repeat. This dichotomous chromatin organization is conserved among humans despite the marked divergence of the underlying alpha-satellite organization and is similarly conserved in gibbon centromeres that lack alpha-satellite repeats, indicating that functional conservation within centromeres is mediated at the level of chromatin, not DNA sequence.HighlightsDichotomous accessible and compacted chromatin (dichromatin) marks centromere coresHighly accessible chromatin patches punctuate sites of kinetochore attachmentDichromatin can form irrespective of CENP-B occupancyConservation within centromeres is mediated at the level of chromatin, not DNA

To better understand the pattern of primate genome structural variation, we sequenced and assembled using multiple long-read sequencing technologies the genomes of eight nonhuman primate species, including New World monkeys (owl monkey and marmoset), Old World monkey (macaque), Asian apes (orangutan and gibbon), and African ape lineages (gorilla, bonobo, and chimpanzee). Compared to the human genome, we identified 1,338,997 lineage-specific fixed structural variants (SVs) disrupting 1,561 protein-coding genes and 136,932 regulatory elements, including the most complete set of human-specific fixed differences. Across 50 million years of primate evolution, we estimate that 819.47 Mbp or ~27% of the genome has been affected by SVs based on analysis of these primate lineages. We identify 1,607 structurally divergent regions (SDRs) wherein recurrent structural variation contributes to creating SV hotspots where genes are recurrently lost (CARDs,ABCD7,OLAH) and new lineage-specific genes are generated (e.g.,CKAP2,NEK5) and have become targets of rapid chromosomal diversification and positive selection (e.g.,RGPDs). High-fidelity long-read sequencing has made these dynamic regions of the genome accessible for sequence-level analyses within and between primate species for the first time.

The genetic mechanisms underlying the expansion in size and complexity of the human brain remains poorly understood. L1 retrotransposons are a source of divergent genetic information in hominoid genomes, but their importance in physiological functions and their contribution to human brain evolution is largely unknown. Using multi-omic profiling we here demonstrate that L1-promoters are dynamically active in the developing and adult human brain. L1s generate hundreds of developmentally regulated and cell-type specific transcripts, many which are co-opted as chimeric transcripts or regulatory RNAs. One L1-derived lncRNA, LINC01876, is a human-specific transcript expressed exclusively during brain development. CRISPRi-silencing of LINC01876 results in reduced size of cerebral organoids and premature differentiation of neural progenitors, implicating L1s in human-specific developmental processes. In summary, our results demonstrate that L1-derived transcripts provide a previously undescribed layer of primate- and human-specific transcriptome complexity that contributes to the functional diversification of the human brain.

Rare genomic disorders (RGDs) confer elevated risk for neurodevelopmental psychiatric disorders. In this era of intense genomics discoveries, the landscape of RGDs is rapidly evolving. However, there has not been comparable progress to date in scalable, harmonized phenotyping methods. As a result, beyond associations with categorical diagnoses, the effects on dimensional traits remain unclear for many RGDs. The nature and specificity of RGD effects on cognitive and behavioral traits is an area of intense investigation: RGDs are frequently associated with more than one psychiatric condition, and those studied to date affect, to varying degrees, a broad range of developmental and cognitive functions. Although many RGDs have large effects, phenotypic expression is typically influenced by additional genomic and environmental factors. There is emerging evidence that using polygenic risk scores in individuals with RGDs offers opportunities to refine prediction, thus allowing for the identification of those at greatest risk of psychiatric illness. However, translation into the clinic is hindered by roadblocks, which include limited genetic testing in clinical psychiatry, and the lack of guidelines for following individuals with RGDs, who are at high risk of developing psychiatric symptoms. The Genes to Mental Health Network (G2MH) is a newly funded National Institute of Mental Health initiative that will collect, share, and analyze large-scale data sets combining genomics and dimensional measures of psychopathology spanning diverse populations and geography. The authors present here the most recent understanding of the effects of RGDs on dimensional behavioral traits and risk for psychiatric conditions and discuss strategies that will be pursued within the G2MH network, as well as how expected results can be translated into clinical practice to improve patient outcomes.

Typical genotyping workflows map reads to a reference genome before identifying genetic variants. Generating such alignments introduces reference biases and comes with substantial computational burden. Furthermore, short-read lengths limit the ability to characterize repetitive genomic regions, which are particularly challenging for fast k-mer-based genotypers. In the present study, we propose a new algorithm, PanGenie, that leverages a haplotype-resolved pangenome reference together with k-mer counts from short-read sequencing data to genotype a wide spectrum of genetic variation—a process we refer to as genome inference. Compared with mapping-based approaches, PanGenie is more than 4 times faster at 30-fold coverage and achieves better genotype concordances for almost all variant types and coverages tested. Improvements are especially pronounced for large insertions (≥50 bp) and variants in repetitive regions, enabling the inclusion of these classes of variants in genome-wide association studies. PanGenie efficiently leverages the increasing amount of haplotype-resolved assemblies to unravel the functional impact of previously inaccessible variants while being faster compared with alignment-based workflows.

Studies of de novo mutation (DNM) have typically excluded some of the most repetitive and complex regions of the genome because these regions cannot be unambiguously mapped with short-read sequencing data. To better understand the genome-wide pattern of DNM, we generated long-read sequence data from an autism parent-child quad with an affected female where no pathogenic variant had been discovered in short-read Illumina sequence data. We deeply sequenced all four individuals by using three sequencing platforms (Illumina, Oxford Nanopore, and Pacific Biosciences) and three complementary technologies (Strand-seq, optical mapping, and 10X Genomics). Using long-read sequencing, we initially discovered and validated 171 DNMs across two children—a 20% increase in the number of de novo single-nucleotide variants (SNVs) and indels when compared to short-read callsets. The number of DNMs further increased by 5% when considering a more complete human reference (T2T-CHM13) because of the recovery of events in regions absent from GRCh38 (e.g., three DNMs in heterochromatic satellites). In total, we validated 195 de novo germline mutations and 23 potential post-zygotic mosaic mutations across both children; the overall true substitution rate based on this integrated callset is at least 1.41 × 10(−8) substitutions per nucleotide per generation. We also identified six de novo insertions and deletions in tandem repeats, two of which represent structural variants. We demonstrate that long-read sequencing and assembly, especially when combined with a more complete reference genome, increases the number of DNMs by >25% compared to previous studies, providing a more complete catalog of DNM compared to short-read data alone.

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.

Complex structural variants (CSVs) are genomic alterations that have more than two breakpoints and are considered as the simultaneous occurrence of simple structural variants. However, detecting the compounded mutational signals of CSVs is challenging through a commonly used model-match strategy. As a result, there has been limited progress for CSV discovery compared with simple structural variants. Here, we systematically analyzed the multi-breakpoint connection feature of CSVs, and proposed Mako, utilizing a bottom-up guided model-free strategy, to detect CSVs from paired-end short-read sequencing. Specifically, we implemented a graph-based pattern growth approach, where the graph depicts potential breakpoint connections, and pattern growth enables CSV detection without pre-defined models. Comprehensive evaluations on both simulated and real datasets revealed that Mako outperformed other algorithms. Notably, validation rates of CSVs on real data based on experimental and computational validations as well as manual inspections are around 70%, where the medians of experimental and computational breakpoint shift are 13 bp and 26 bp, respectively. Moreover, the Mako CSV subgraph effectively characterized the breakpoint connections of a CSV event and uncovered a total of 15 CSV types, including two novel types of adjacent segment swap and tandem dispersed duplication. Further analysis of these CSVs also revealed the impact of sequence homology on the formation of CSVs. Mako is publicly available at https://github.com/xjtu-omics/Mako.

The first telomere-to-telomere (T2T) human genome assembly (T2T-CHM13) release was a milestone in human genomics. The T2T-CHM13 genome assembly extends our understanding of telomeres, centromeres, segmental duplication, and other complex regions. The current human genome reference (GRCh38) has been widely used in various human genomic studies. However, the large-scale genomic differences between these two important genome assemblies are not characterized in detail yet. Here, we identify 590 discrepant regions (∼226 Mbp) in total. In addition to the previously reported ‘non-syntenic’ regions, we identify 67 additional large-scale discrepant regions and precisely categorize them into four structural types with a newly developed website tool (SynPlotter). The discrepant regions (∼20.4 Mbp) excluding telomeric and centromeric regions are highly structurally polymorphic in humans, where copy number variation are likely associated with various human disease and disease susceptibility, such as immune and neurodevelopmental disorders. The analyses of a newly identified discrepant region—theKLRCgene cluster—shows that the depletion ofKLRC2by a single deletion event is associated with natural killer cell differentiation in ∼20% of humans. Meanwhile, the rapid amino acid replacements within KLRC3 is consistent with the action of natural selection during primate evolution. Our study furthers our understanding of the large-scale structural variation differences between these two crucial human reference genomes and future interpretation of studies of human genetic variation.

CHD8, a major autism gene, functions in chromatin remodelling and has various roles involving several biological pathways. Therefore, unsurprisingly, previous studies have shown that intellectual developmental disorder with autism and macrocephaly (IDDAM), the syndrome caused by pathogenic variants in CHD8, consists of a broad range of phenotypic abnormalities. We collected and reviewed 106 individuals with IDDAM, including 36 individuals not previously published, thus enabling thorough genotype–phenotype analyses, involving the CHD8 mutation spectrum, characterization of the CHD8 DNA methylation episignature, and the systematic analysis of phenotypes collected in Human Phenotype Ontology (HPO). We identified 29 unique nonsense, 25 frameshift, 24 missense, and 12 splice site variants. Furthermore, two unique inframe deletions, one larger deletion (exons 26–28), and one translocation were observed. Methylation analysis was performed for 13 patients, 11 of which showed the previously established episignature for IDDAM (85%) associated with CHD8 haploinsufficiency, one analysis was inconclusive, and one showing a possible gain-of-function signature instead of the expected haploinsufficiency signature was observed. Consistent with previous studies, phenotypical abnormalities affected multiple organ systems. Many neurological abnormalities, like intellectual disability (68%) and hypotonia (29%) were observed, as well as a wide variety of behavioural abnormalities (88%). Most frequently observed behavioural problems included autism spectrum disorder (76%), short attention span (32%), abnormal social behaviour (31%), sleep disturbance (29%) and impaired social interactions (28%). Furthermore, abnormalities in the digestive (53%), musculoskeletal (79%) and genitourinary systems (18%) were noted. Although no significant difference in severity was observed between males and females, individuals with a missense variant were less severely affected. Our study provides an extensive review of all phenotypic abnormalities in patients with IDDAM and provides clinical recommendations, which will be of significant value to individuals with a pathogenic variant in CHD8, their families, and clinicians as it gives a more refined insight into the clinical and molecular spectrum of IDDAM, which is essential for accurate care and counselling.

The human Y chromosome has been notoriously difficult to sequence and assemble because of its complex repeat structure including long palindromes, tandem repeats, and segmental duplications. 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 finished. 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, revealing the complete ampliconic structures ofTSPY, DAZ, andRBMY; 42 additional protein-coding genes, mostly from theTSPYgene family; and an alternating pattern of human satellite 1 and 3 blocks in the heterochromatic Yq12 region. We have combined T2T-Y with a prior assembly of the CHM13 genome and mapped available population variation, clinical variants, and functional genomics data to produce a complete and comprehensive reference sequence for all 24 human chromosomes.

During mammalian development, differences in chromatin state coincide with cellular differentiation and reflect changes in the gene regulatory landscape1. In the developing brain, cell fate specification and topographic identity are important for defining cell identity2 and confer selective vulnerabilities to neurodevelopmental disorders3. Here, to identify cell-type-specific chromatin accessibility patterns in the developing human brain, we used a single-cell assay for transposase accessibility by sequencing (scATAC-seq) in primary tissue samples from the human forebrain. We applied unbiased analyses to identify genomic loci that undergo extensive cell-type- and brain-region-specific changes in accessibility during neurogenesis, and an integrative analysis to predict cell-type-specific candidate regulatory elements. We found that cerebral organoids recapitulate most putative cell-type-specific enhancer accessibility patterns but lack many cell-type-specific open chromatin regions that are found in vivo. Systematic comparison of chromatin accessibility across brain regions revealed unexpected diversity among neural progenitor cells in the cerebral cortex and implicated retinoic acid signalling in the specification of neuronal lineage identity in the prefrontal cortex. Together, our results reveal the important contribution of chromatin state to the emerging patterns of cell type diversity and cell fate specification and provide a blueprint for evaluating the fidelity and robustness of cerebral organoids as a model for cortical development., Analysis of chromatin state at a single-cell level in samples of developing human forebrain demonstrate both cell-type-specific and region-specific changes during neurogenesis.

Most genetic studies consider autism spectrum disorder (ASD) and developmental disorder (DD) separately despite overwhelming comorbidity and shared genetic etiology. Here, we analyzed de novo variants (DNVs) from 15,560 ASD (6,557 from SPARK) and 31,052 DD trios independently and also combined as broader neurodevelopmental disorders (NDDs) using three models. We identify 615 NDD candidate genes (false discovery rate [FDR] < 0.05) supported by ≥1 models, including 138 reaching Bonferroni exome-wide significance ( P < 3.64e–7) in all models. The genes group into five functional networks associating with different brain developmental lineages based on single-cell nuclei transcriptomic data. We find no evidence for ASD-specific genes in contrast to 18 genes significantly enriched for DD. There are 53 genes that show mutational bias, including enrichments for missense ( n = 41) or truncating ( n = 12) DNVs. We also find 10 genes with evidence of male- or female-bias enrichment, including 4 X chromosome genes with significant female burden ( DDX3X , MECP2 , WDR45 , and HDAC8) . This large-scale integrative analysis identifies candidates and functional subsets of NDD genes.

Autism spectrum disorder (ASD) represents a group of neurodevelopmental phenotypes with a strong genetic component. An excess of likely gene-disruptive (LGD) mutations in GIGYF1 was implicated in ASD. Here, we report that GIGYF1 is the second-most mutated gene among known ASD high-confidence risk genes. We investigated the inheritance of 46 GIGYF1 LGD variants, including the highly recurrent mutation c.333del:p.L111Rfs*234. Inherited GIGYF1 heterozygous LGD variants were 1.8 times more common than de novo mutations. Among individuals with ASD, cognitive impairments were less likely in those with GIGYF1 LGD variants relative to those with other high-confidence gene mutations. Using a Gigyf1 conditional KO mouse model, we showed that haploinsufficiency in the developing brain led to social impairments without significant cognitive impairments. In contrast, homozygous mice showed more severe social disability as well as cognitive impairments. Gigyf1 deficiency in mice led to a reduction in the number of upper-layer cortical neurons, accompanied by a decrease in proliferation and increase in differentiation of neural progenitor cells. We showed that GIGYF1 regulated the recycling of IGF-1R to the cell surface. KO of GIGYF1 led to a decreased level of IGF-1R on the cell surface, disrupting the IGF-1R/ERK signaling pathway. In summary, our findings show that GIGYF1 is a regulator of IGF-1R recycling. Haploinsufficiency of GIGYF1 was associated with autistic behavior, likely through interference with IGF-1R/ERK signaling pathway.

Rapid genetic testing in the critical care setting enables targeted evaluations, directs therapies, and helps families and care providers make informed decisions about goals of care. We tested whether we could perform ultra-rapid assessment of genetic risk for a Mendelian condition, based on information from an affected sibling, in a newborn via whole-genome sequencing using the Oxford Nanopore platform. By optimization of the DNA extraction and library preparation steps paired with targeted analysis, we were able to demonstrate within three hours of birth that the newborn was neither affected nor a carrier for variants underlying acrodermatitis enteropathica. This proof-of-concept experiment demonstrates how prior knowledge of familial variants can be used to rapidly evaluate an individual at-risk for a genetic disease.

There are more than 55,000 variable number tandem repeats (VNTRs) in the human genome, notable for both their striking polymorphism and mutability. Despite their role in human evolution and genomic variation, they have yet to be studied collectively and in detail, partially owing to their large size, variability, and predominant location in noncoding regions. Here, we examine 467 VNTRs that are human-specific expansions, unique to one location in the genome, and not associated with retrotransposons. We leverage publicly available long-read genomes, including from the Human Genome Structural Variant Consortium, to ascertain the exact nucleotide composition of these VNTRs and compare their composition of alleles. We then confirm repeat unit composition in more than 3000 short-read samples from the 1000 Genomes Project. Our analysis reveals that these VNTRs contain highly structured repeat motif organization, modified by frequent deletion and duplication events. Although overall VNTR compositions tend to remain similar between 1000 Genomes Project superpopulations, we describe a notable exception with substantial differences in repeat composition (in PCBP3), as well as several VNTRs that are significantly different in length between superpopulations (in ART1, PROP1, DYNC2I1, and LOC102723906). We also observe that most of these VNTRs are expanded in archaic human genomes, yet remain stable in length between single generations. Collectively, our findings indicate that repeat motif variability, repeat composition, and repeat length are all informative modalities to consider when characterizing VNTRs and their contribution to genomic variation.

Summary The visualization and analysis of genomic repeats is typically accomplished using dot plots; however, the emergence of telomere-to-telomere assemblies with multi-megabase repeats requires new visualization strategies. Here, we introduce StainedGlass, which can generate publication-quality figures and interactive visualizations that depict the identity and orientation of multi-megabase tandem repeat structures at a genome-wide scale. The tool can rapidly reveal higher-order structures and improve the inference of evolutionary history for some of the most complex regions of genomes. Availability and implementation StainedGlass is implemented using Snakemake and available open source under the MIT license at https://mrvollger.github.io/StainedGlass/. Supplementary information Supplementary data are available at Bioinformatics online.

Stress granules (SGs) are cytoplasmic assemblies in response to a variety of stressors. We report a new neurodevelopmental disorder (NDD) with common features of language problems, intellectual disability, and behavioral issues caused by de novo likely gene-disruptive variants in UBAP2L , which encodes an essential regulator of SG assembly. Ubap2l haploinsufficiency in mouse led to social and cognitive impairments accompanied by disrupted neurogenesis and reduced SG formation during early brain development. On the basis of data from 40,853 individuals with NDDs, we report a nominally significant excess of de novo variants within 29 genes that are not implicated in NDDs, including 3 essential genes ( G3BP1 , G3BP2 , and UBAP2L ) in the core SG interaction network. We validated that NDD-related de novo variants in newly implicated and known NDD genes, such as CAPRIN1 , disrupt the interaction of the core SG network and interfere with SG formation. Together, our findings suggest the common SG pathology in NDDs.

Immunoglobulins (IGs), crucial components of the adaptive immune system, are encoded by three genomic loci. However, the complexity of the IG loci severely limits the effective use of short read sequencing, limiting our knowledge of population diversity in these loci. We leveraged existing long read whole-genome sequencing (WGS) data, fosmid technology, and IG targeted single-molecule, real-time (SMRT) long-read sequencing (IG-Cap) to create haplotype-resolved assemblies of the IG Lambda (IGL) locus from 6 ethnically diverse individuals. In addition, we generated 10 diploid assemblies of IGL from a diverse cohort of individuals utilizing IG-cap. From these 16 individuals, we identified significant allelic diversity, including 37 novel IGLV alleles. In addition, we observed highly elevated single nucleotide variation (SNV) in IGLV genes relative to IGL intergenic and genomic background SNV density. By comparing SNV calls between our high quality assemblies and existing short read datasets from the same individuals, we show a high propensity for false-positives in the short read datasets. Finally, for the first time, we nucleotide-resolved common 5-10 Kb duplications in the IGLC region that contain functional IGLJ and IGLC genes. Together these data represent a significant advancement in our understanding of genetic variation and population diversity in the IGL locus.

Motivation Highly contiguous de novo phased diploid genome assemblies are now feasible for large numbers of species and individuals. Methods are needed to validate assembly accuracy and detect misassemblies with orthologous sequencing data to allow for confident downstream analyses. Results We developed GAVISUNK, an open-source pipeline that detects misassemblies and produces a set of reliable regions genome-wide by assessing concordance of distances between unique k-mers in Pacific Biosciences high-fidelity assemblies and raw Oxford Nanopore Technologies reads. Availability and implementation GAVISUNK is available at https://github.com/pdishuck/GAVISUNK. Supplementary information Supplementary data are available at Bioinformatics online.

The Telomere-to-Telomere consortium recently assembled the first truly complete sequence of a human genome. To resolve the most complex repeats, this project relied on manual integration of ultra-long Oxford Nanopore sequencing reads with a high-resolution assembly graph built from long, accurate PacBio HiFi reads. We have improved and automated this strategy in Verkko, an iterative, graph-based pipeline for assembling complete, diploid genomes. Verkko begins with a multiplex de Bruijn graph built from long, accurate reads and progressively simplifies this graph via the integration of ultra-long reads and haplotype-specific markers. The result is a phased, diploid assembly of both haplotypes, with many chromosomes automatically assembled from telomere to telomere. Running Verkko on the HG002 human genome resulted in 20 of 46 diploid chromosomes assembled without gaps at 99.9997% accuracy. The complete assembly of diploid genomes is a critical step towards the construction of comprehensive pangenome databases and chromosome-scale comparative genomics.

MotivationHighly contiguous de novo genome assemblies are now feasible for large numbers of species and individuals. Methods are needed to validate assembly accuracy and detect misassemblies with orthologous sequencing data to allow for confident downstream analyses.ResultsWe developed GAVISUNK, an open-source pipeline that detects misassemblies and produces a set of reliable regions genome-wide by assessing concordance of distances between unique k-mers in Pacific Biosciences high-fidelity (HiFi) assemblies and raw Oxford Nanopore Technologies reads.AvailabilityGAVISUNK is available at https://github.com/pdishuck/GAVISUNK.Contacteee@gs.washington.edu

The divergence of chimpanzee and bonobo provides one of the few examples of recent hominid speciation1,2. Here we describe a fully annotated, high-quality bonobo genome assembly, which was constructed without guidance from reference genomes by applying a multiplatform genomics approach. We generate a bonobo genome assembly in which more than 98% of genes are completely annotated and 99% of the gaps are closed, including the resolution of about half of the segmental duplications and almost all of the full-length mobile elements. We compare the bonobo genome to those of other great apes1,3–5 and identify more than 5,569 fixed structural variants that specifically distinguish the bonobo and chimpanzee lineages. We focus on genes that have been lost, changed in structure or expanded in the last few million years of bonobo evolution. We produce a high-resolution map of incomplete lineage sorting and estimate that around 5.1% of the human genome is genetically closer to chimpanzee or bonobo and that more than 36.5% of the genome shows incomplete lineage sorting if we consider a deeper phylogeny including gorilla and orangutan. We also show that 26% of the segments of incomplete lineage sorting between human and chimpanzee or human and bonobo are non-randomly distributed and that genes within these clustered segments show significant excess of amino acid replacement compared to the rest of the genome., A high-quality bonobo genome assembly provides insights into incomplete lineage sorting in hominids and its relevance to gene evolution and the genetic relationship among living hominids.

Virtually all genome sequencing efforts in national biobanks, complex and Mendelian disease programs, and medical genetic initiatives are reliant upon short-read whole-genome sequencing (srWGS), which presents challenges for the detection of structural variants (SVs) relative to emerging long-read WGS (lrWGS) technologies. Given this ubiquity of srWGS in large-scale genomics initiatives, we sought to establish expectations for routine SV detection from this data type by comparison with lrWGS assembly, as well as to quantify the genomic properties and added value of SVs uniquely accessible to each technology. Analyses from the Human Genome Structural Variation Consortium (HGSVC) of three families captured ~11,000 SVs per genome from srWGS and ~25,000 SVs per genome from lrWGS assembly. Detection power and precision for SV discovery varied dramatically by genomic context and variant class: 9.7% of the current GRCh38 reference is defined by segmental duplication (SD) and simple repeat (SR), yet 91.4% of deletions that were specifically discovered by lrWGS localized to these regions. Across the remaining 90.3% of reference sequence, we observed extremely high (93.8%) concordance between technologies for deletions in these datasets. In contrast, lrWGS was superior for detection of insertions across all genomic contexts. Given that non-SD/SR sequences encompass 95.9% of currently annotated disease-associated exons, improved sensitivity from lrWGS to discover novel pathogenic deletions in these currently interpretable genomic regions is likely to be incremental. However, these analyses highlight the considerable added value of assembly-based lrWGS to create new catalogs of insertions and transposable elements, as well as disease-associated repeat expansions in genomic sequences that were previously recalcitrant to routine assessment.

SUMMARYDiverse inbred mouse strains are among the foremost models for biomedical research, yet genome characterization of many strains has been fundamentally lacking in comparison to human genomics research. In particular, the discovery and cataloging of structural variants is incomplete, limiting the discovery of potentially causative alleles for phenotypic variation across individuals. Here, we utilized long-read sequencing to resolve genome-wide structural variants (SVs, variants ≥ 50 bp) in 20 genetically distinct inbred mice. We report 413,758 site-specific SVs that affect 13% (356 Mbp) of the current mouse reference assembly, including 510 previously unannotated variants which alter coding sequences. We find that 39% of SVs are attributed to transposable element (TE) variation accounting for 75% of bases altered by SV. We then utilized this callset to investigate the impact of TE heterogeneity on mouse embryonic stem cells (mESCs), and find multiple TE classes that influence chromatin accessibility across loci. We also identify strain-specific transcription start sites originating in polymorphic TEs that modify gene expression. Our work provides the first long-read based analysis of mouse SVs and illustrates that previously unresolved TEs underlie epigenetic and transcriptome differences in mESCs.

Since the introduction of next‐generation sequencing, an increasing number of disorders have been discovered to have genetic etiology. To address diverse clinical questions and coordinate research activities that arise with the identification of these rare disorders, we developed the Human Disease Genes website series (HDG website series): an international digital library that records detailed information on the clinical phenotype of novel genetic variants in the human genome (https://humandiseasegenes.info/). Each gene website is moderated by a dedicated team of clinicians and researchers, focused on specific genes, and provides up‐to‐date—including unpublished—clinical information. The HDG website series is expanding rapidly with 424 genes currently adopted by 325 moderators from across the globe. On average, a gene website has detailed phenotypic information of 14.4 patients. There are multiple examples of added value, one being the ARID1B gene website, which was recently utilized in research to collect clinical information of 81 new patients. Additionally, several gene websites have more data available than currently published in the literature. In conclusion, the HDG website series provides an easily accessible, open and up‐to‐date clinical data resource for patients with pathogenic variants of individual genes. This is a valuable resource not only for clinicians dealing with rare genetic disorders such as developmental delay and autism, but other professionals working in diagnostics and basic research. Since the HDG website series is a dynamic platform, its data also include the phenotype of yet unpublished patients curated by professionals providing higher quality clinical detail to improve management of these rare disorders.

Unlike copy number variants (CNVs), inversions remain an underexplored genetic variation class. By integrating multiple genomic technologies, we discover 729 inversions in 41 human genomes. Approximately 85% of inversions

BackgroundWerner syndrome (WS) is an autosomal recessive progeroid syndrome caused by variants inWRN. The International Registry of Werner Syndrome has identified biallelic pathogenic variants in 179/188 cases of classical WS. In the remaining nine cases, only one heterozygous pathogenic variant has been identified.MethodsTargeted long-read sequencing (T-LRS) on an Oxford Nanopore platform was used to search for a second pathogenic variant inWRN. Previously, T-LRS was successfully used to identify missing variants and analyse complex rearrangements.ResultsWe identified a second pathogenic variant in eight of nine unsolved WS cases. In five cases, T-LRS identified intronic splice variants that were confirmed by either RT-PCR or exon trapping to affect splicing; in one case, T-LRS identified a 339 kbp deletion, and in two cases, pathogenic missense variants. Phasing of long reads predicted all newly identified variants were on a different haplotype than the previously known variant. Finally, in one case, RT-PCR previously identified skipping of exon 20; however, T-LRS did not detect a pathogenic DNA sequence variant.ConclusionT-LRS is an effective method for identifying missing pathogenic variants. Although limitations with computational prediction algorithms can hinder the interpretation of variants, T-LRS is particularly effective in identifying intronic variants.

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.

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.

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.