Search

Short tandem repeat (STR) expansions are associated with more than 60 genetic disorders. The size and stability of these expansions correlate with the severity and age of onset of the disease. Therefore, being able to accurately detect the absolute length of STRs is important. Current diagnostic assays include laborious lab experiments, including repeat-primed PCR and Southern blotting, that still cannot precisely determine the exact length of very long repeat expansions. Optical genome mapping (OGM) is a cost-effective and easy-to-use alternative to traditional cytogenetic techniques and allows the comprehensive detection of chromosomal aberrations and structural variants >500 bp in length, including insertions, deletions, duplications, inversions, translocations, and copy number variants. Here, we provide methodological guidance for preparing samples and performing OGM as well as running the analysis pipelines and using the specific repeat expansion workflows to determine the exact repeat length of repeat expansions expanded beyond 500 bp. Together these protocols provide all details needed to analyze the length and stability of any repeat expansion with an expected repeat size difference from the expected wild-type allele of >500 bp. © 2024 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Genomic ultra-high-molecular-weight DNA isolation, labeling, and staining Basic Protocol 2: Data generation and genome mapping using the Bionano Saphyr® System Basic Protocol 3: Manual De Novo Assembly workflow Basic Protocol 4: Local guided assembly workflow Basic Protocol 5: EnFocus Fragile X workflow Basic Protocol 6: Molecule distance script workflow., (© 2024 The Author(s). Current Protocols published by Wiley Periodicals LLC.)

Acute leukemia is a particularly problematic collection of hematological cancers, and, while somewhat rare, the survival rate of patients is typically abysmal without bone marrow transplantation. Furthermore, traditional chemotherapies used as standard-of-care for patients cause significant side effects. Understanding the evolution of leukemia to identify novel targets and, therefore, drug treatment regimens is a significant medical need. Genomic rearrangements and other structural variations (SVs) have long been known to be causative and pathogenic in multiple types of cancer, including leukemia. These SVs may be involved in cancer initiation, progression, clonal evolution, and drug resistance, and a better understanding of SVs from individual patients may help guide therapeutic options. Here, we show the utilization of optical genome mapping (OGM) to detect known and novel SVs in the samples of patients with leukemia. Importantly, this technology provides an unprecedented level of granularity and quantitation unavailable to other current techniques and allows for the unbiased detection of novel SVs, which may be relevant to disease pathogenesis and/or drug resistance. Coupled with the chemosensitivities of these samples to FDA-approved oncology drugs, we show how an impartial integrative analysis of these diverse datasets can be used to associate the detected genomic rearrangements with multiple drug sensitivity profiles. Indeed, an insertion in the gene MUSK is shown to be associated with increased sensitivity to the clinically relevant agent Idarubicin, while partial tandem duplication events in the KMT2A gene are related to the efficacy of another frontline treatment, Cytarabine.

Structural variations (SVs) play a key role in the pathogenicity of hematological malignancies. Standard-of-care (SOC) methods such as karyotyping and fluorescence in situ hybridization (FISH), which have been employed globally for the past three decades, have significant limitations in terms of resolution and the number of recurrent aberrations that can be simultaneously assessed, respectively. Next-generation sequencing (NGS)-based technologies are now widely used to detect clinically significant sequence variants but are limited in their ability to accurately detect SVs. Optical genome mapping (OGM) is an emerging technology enabling the genome-wide detection of all classes of SVs at a significantly higher resolution than karyotyping and FISH. OGM requires neither cultured cells nor amplification of DNA, addressing the limitations of culture and amplification biases. This study reports the clinical validation of OGM as a laboratory-developed test (LDT) according to stringent regulatory (CAP/CLIA) guidelines for genome-wide SV detection in different hematological malignancies. In total, 60 cases with hematological malignancies (of various subtypes), 18 controls, and 2 cancer cell lines were used for this study. Ultra-high-molecular-weight DNA was extracted from the samples, fluorescently labeled, and run on the Bionano Saphyr system. A total of 215 datasets, Inc.luding replicates, were generated, and analyzed successfully. Sample data were then analyzed using either disease-specific or pan-cancer-specific BED files to prioritize calls that are known to be diagnostically or prognostically relevant. Sensitivity, specificity, and reproducibility were 100%, 100%, and 96%, respectively. Following the validation, 14 cases and 10 controls were run and analyzed using OGM at three outside laboratories showing reproducibility of 96.4%. OGM found more clinically relevant SVs compared to SOC testing due to its ability to detect all classes of SVs at higher resolution. The results of this validation study demonstrate the superiority of OGM over traditional SOC methods for the detection of SVs for the accurate diagnosis of various hematological malignancies.

The recommended practice for individuals suspected of a genetic etiology for disorders including unexplained developmental delay/intellectual disability (DD/ID), autism spectrum disorders (ASD), and multiple congenital anomalies (MCA) involves a genetic testing workflow including chromosomal microarray (CMA), Fragile-X testing, karyotype analysis, and/or sequencing-based gene panels. Since genomic imbalances are often found to be causative, CMA is recommended as first tier testing for many indications. Optical genome mapping (OGM) is an emerging next generation cytogenomic technique that can detect not only copy number variants (CNVs), triploidy and absence of heterozygosity (AOH) like CMA, but can also define the location of duplications, and detect other structural variants (SVs), including balanced rearrangements and repeat expansions/contractions. This study compares OGM to CMA for clinically reported genomic variants, some of these samples also have structural characterization by fluorescence in situ hybridization (FISH). OGM was performed on IRB approved, de-identified specimens from 55 individuals with genomic abnormalities previously identified by CMA (61 clinically reported abnormalities). SVs identified by OGM were filtered by a control database to remove polymorphic variants and against an established gene list to prioritize clinically relevant findings before comparing with CMA and FISH results. OGM results showed 100% concordance with CMA findings for pathogenic variants and 98% concordant for all pathogenic/likely pathogenic/variants of uncertain significance (VUS), while also providing additional insight into the genomic structure of abnormalities that CMA was unable to provide. OGM demonstrates equivalent performance to CMA for CNV and AOH detection, enhanced by its ability to determine the structure of the genome. This work adds to an increasing body of evidence on the analytical validity and ability to detect clinically relevant abnormalities identified by CMA. Moreover, OGM identifies translocations, structures of duplications and complex CNVs intractable by CMA, yielding additional clinical utility.

Homologous recombination deficiency (HRD) is characterized by the inability of a cell to repair the double-stranded breaks using the homologous recombination repair (HRR) pathway. The deficiency of the HRR pathway results in defective DNA repair, leading to genomic instability and tumorigenesis. The presence of HRD has been found to make tumors sensitive to ICL-inducing platinum-based therapies and poly(adenosine diphosphate [ADP]-ribose) polymerase (PARP) inhibitors (PARPi). However, there are no standardized methods to measure and report HRD phenotypes. Herein, we compare optical genome mapping (OGM), chromosomal microarray (CMA), and a 523-gene NGS panel for HRD score calculations. This retrospective study included the analysis of 196 samples, of which 10 were gliomas, 176 were hematological malignancy samples, and 10 were controls. The 10 gliomas were evaluated with both CMA and OGM, and 30 hematological malignancy samples were evaluated with both the NGS panel and OGM. To verify the scores in a larger cohort, 135 cases were evaluated with the NGS panel and 71 cases with OGM. The HRD scores were calculated using a combination of three HRD signatures that included loss of heterozygosity (LOH), telomeric allelic imbalance (TAI), and large-scale transitions (LST). In the ten glioma cases analyzed with OGM and CMA using the same DNA (to remove any tumor percentage bias), the HRD scores (mean ± SEM) were 13.2 (±4.2) with OGM compared to 3.7 (±1.4) with CMA. In the 30 hematological malignancy cases analyzed with OGM and the 523-gene NGS panel, the HRD scores were 7.6 (±2.2) with OGM compared to 2.6 (±0.8) with the 523-gene NGS panel. OGM detected 70.8% and 66.8% of additional variants that are considered HRD signatures in gliomas and hematological malignancies, respectively. The higher sensitivity of OGM to capture HRD signature variants might enable a more accurate and precise correlation with response to PARPi and platinum-based drugs. This study reveals HRD signatures that are cryptic to current standard of care (SOC) methods used for assessing the HRD phenotype and presents OGM as an attractive alternative with higher resolution and sensitivity to accurately assess the HRD phenotype.

Background: The cancer genome is commonly altered with thousands of structural rearrangements including insertions, deletions, translocation, inversions, duplications, and copy number variations. Thus, structural variant (SV) characterization plays a paramount role in cancer target identification, oncology diagnostics, and personalized medicine. As part of the SEQC2 Consortium effort, the present study established and evaluated a consensus SV call set using a breast cancer reference cell line and matched normal control derived from the same donor, which were used in our companion benchmarking studies as reference samples., Results: We systematically investigated somatic SVs in the reference cancer cell line by comparing to a matched normal cell line using multiple NGS platforms including Illumina short-read, 10X Genomics linked reads, PacBio long reads, Oxford Nanopore long reads, and high-throughput chromosome conformation capture (Hi-C). We established a consensus SV call set of a total of 1788 SVs including 717 deletions, 230 duplications, 551 insertions, 133 inversions, 146 translocations, and 11 breakends for the reference cancer cell line. To independently evaluate and cross-validate the accuracy of our consensus SV call set, we used orthogonal methods including PCR-based validation, Affymetrix arrays, Bionano optical mapping, and identification of fusion genes detected from RNA-seq. We evaluated the strengths and weaknesses of each NGS technology for SV determination, and our findings provide an actionable guide to improve cancer genome SV detection sensitivity and accuracy., Conclusions: A high-confidence consensus SV call set was established for the reference cancer cell line. A large subset of the variants identified was validated by multiple orthogonal methods., (© 2022. The Author(s).)

The Hawaiian monk seal (HMS) is the single extant species of tropical earless seals of the genus Neomonachus. The species survived a severe bottleneck in the late 19th century and experienced subsequent population declines until becoming the subject of a NOAA-led species recovery effort beginning in 1976 when the population was fewer than 1000 animals. Like other recovering species, the Hawaiian monk seal has been reported to have reduced genetic heterogeneity due to the bottleneck and subsequent inbreeding. Here, we report a chromosomal reference assembly for a male animal produced using a variety of methods. The final assembly consisted of 16 autosomes, an X, and portions of the Y chromosomes. We compared variants in this animal to other HMS and to a frequently sequenced human sample, confirming about 12% of the variation seen in man. To confirm that the reference animal was representative of the HMS, we compared his sequence to that of 10 other individuals and noted similarly low variation in all. Variation in the major histocompatibility (MHC) genes was nearly absent compared to the orthologous human loci. Demographic analysis predicts that Hawaiian monk seals have had a long history of small populations preceding the bottleneck, and their current low levels of heterozygosity may indicate specialization to a stable environment. When we compared our reference assembly to that of other species, we observed significant conservation of chromosomal architecture with other pinnipeds, especially other phocids. This reference should be a useful tool for future evolutionary studies as well as the long-term management of this species.

The divergence of chimpanzee and bonobo provides one of the few examples of recent hominid speciation 1,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 apes 1,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.

Genomic structural variants comprise a significant fraction of somatic mutations driving cancer onset and progression. However, such variants are not readily revealed by standard next-generation sequencing. Optical genome mapping (OGM) surpasses short-read sequencing in detecting large (>500 bp) and complex structural variants (SVs) but requires isolation of ultra-high-molecular-weight DNA from the tissue of interest. We have successfully applied a protocol involving a paramagnetic nanobind disc to a wide range of solid tumors. Using as little as 6.5 mg of input tumor tissue, we show successful extraction of high-molecular-weight genomic DNA that provides a high genomic map rate and effective coverage by optical mapping. We demonstrate the system's utility in identifying somatic SVs affecting functional and cancer-related genes for each sample. Duplicate/triplicate analysis of select samples shows intra-sample reliability but also intra-sample heterogeneity. We also demonstrate that simply filtering SVs based on a GRCh38 human control database provides high positive and negative predictive values for true somatic variants. Our results indicate that the solid tissue DNA extraction protocol, OGM and SV analysis can be applied to a wide variety of solid tumors to capture SVs across the entire genome with functional importance in cancer prognosis and treatment.

Inbred animals were historically chosen for genome analysis to circumvent assembly issues caused by haplotype variation but this resulted in a composite of the two genomes. Here we report a haplotype-aware scaffolding and polishing pipeline which was used to create haplotype-resolved, chromosome-level genome assemblies of Angus (taurine) and Brahman (indicine) cattle subspecies from contigs generated by the trio binning method. These assemblies reveal structural and copy number variants that differentiate the subspecies and that variant detection is sensitive to the specific reference genome chosen. Six genes with immune related functions have additional copies in the indicine compared with taurine lineage and an indicus-specific extra copy of fatty acid desaturase is under positive selection. The haplotyped genomes also enable transcripts to be phased to detect allele-specific expression. This work exemplifies the value of haplotype-resolved genomes to better explore evolutionary and functional variations.

The incomplete identification of structural variants (SVs) from whole-genome sequencing data limits studies of human genetic diversity and disease association. Here, we apply a suite of long-read, short-read, strand-specific sequencing technologies, optical mapping, and variant discovery algorithms to comprehensively analyze three trios to define the full spectrum of human genetic variation in a haplotype-resolved manner. We identify 818,054 indel variants (<50 bp) and 27,622 SVs (≥50 bp) per genome. We also discover 156 inversions per genome and 58 of the inversions intersect with the critical regions of recurrent microdeletion and microduplication syndromes. Taken together, our SV callsets represent a three to sevenfold increase in SV detection compared to most standard high-throughput sequencing studies, including those from the 1000 Genomes Project. The methods and the dataset presented serve as a gold standard for the scientific community allowing us to make recommendations for maximizing structural variation sensitivity for future genome sequencing studies.

Human chromosome 15q25 is involved in several disease-associated structural rearrangements, including microdeletions and chromosomal markers with inverted duplications. Using comparative fluorescence in situ hybridization, strand-sequencing, single-molecule, real-time sequencing and Bionano optical mapping analyses, we investigated the organization of the 15q25 region in human and nonhuman primates. We found that two independent inversions occurred in this region after the fission event that gave rise to phylogenetic chromosomes XIV and XV in humans and great apes. One of these inversions is still polymorphic in the human population today and may confer differential susceptibility to 15q25 microdeletions and inverted duplications. The inversion breakpoints map within segmental duplications containing core duplicons of the GOLGA gene family and correspond to the site of an ancestral centromere, which became inactivated about 25 million years ago. The inactivation of this centromere likely released segmental duplications from recombination repression typical of centromeric regions. We hypothesize that this increased the frequency of ectopic recombination creating a hotspot of hominid inversions where dispersed GOLGA core elements now predispose this region to recurrent genomic rearrangements associated with disease., Competing Interests: I have read the journal's policy and the authors of this manuscript have the following competing interests: EEE is on the scientific advisory board (SAB) of DNAnexus, Inc.

In the originally published version of this Article, the sequenced axolotl strain (the homozygous white mutant) was denoted as 'D/D' rather than 'd/d' in Fig. 1a and the accompanying legend, the main text and the Methods section. The original Article has been corrected online.

Genetic studies of human evolution require high-quality contiguous ape genome assemblies that are not guided by the human reference. We coupled long-read sequence assembly and full-length complementary DNA sequencing with a multiplatform scaffolding approach to produce ab initio chimpanzee and orangutan genome assemblies. By comparing these with two long-read de novo human genome assemblies and a gorilla genome assembly, we characterized lineage-specific and shared great ape genetic variation ranging from single- to mega-base pair-sized variants. We identified ~17,000 fixed human-specific structural variants identifying genic and putative regulatory changes that have emerged in humans since divergence from nonhuman apes. Interestingly, these variants are enriched near genes that are down-regulated in human compared to chimpanzee cerebral organoids, particularly in cells analogous to radial glial neural progenitors., (Copyright © 2018 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Salamanders serve as important tetrapod models for developmental, regeneration and evolutionary studies. An extensive molecular toolkit makes the Mexican axolotl (Ambystoma mexicanum) a key representative salamander for molecular investigations. Here we report the sequencing and assembly of the 32-gigabase-pair axolotl genome using an approach that combined long-read sequencing, optical mapping and development of a new genome assembler (MARVEL). We observed a size expansion of introns and intergenic regions, largely attributable to multiplication of long terminal repeat retroelements. We provide evidence that intron size in developmental genes is under constraint and that species-restricted genes may contribute to limb regeneration. The axolotl genome assembly does not contain the essential developmental gene Pax3. However, mutation of the axolotl Pax3 paralogue Pax7 resulted in an axolotl phenotype that was similar to those seen in Pax3 -/- and Pax7 -/- mutant mice. The axolotl genome provides a rich biological resource for developmental and evolutionary studies.

Accurate and contiguous genome assembly is key to a comprehensive understanding of the processes shaping genomic diversity and evolution. Yet, it is frequently constrained by constitutive heterochromatin, usually characterized by highly repetitive DNA. As a key feature of genome architecture associated with centromeric and subtelomeric regions, it locally influences meiotic recombination. In this study, we assess the impact of large tandem repeat arrays on the recombination rate landscape in an avian speciation model, the Eurasian crow. We assembled two high-quality genome references using single-molecule real-time sequencing (long-read assembly [LR]) and single-molecule optical maps (optical map assembly [OM]). A three-way comparison including the published short-read assembly (SR) constructed for the same individual allowed assessing assembly properties and pinpointing misassemblies. By combining information from all three assemblies, we characterized 36 previously unidentified large repetitive regions in the proximity of sequence assembly breakpoints, the majority of which contained complex arrays of a 14-kb satellite repeat or its 1.2-kb subunit. Using whole-genome population resequencing data, we estimated the population-scaled recombination rate (ρ) and found it to be significantly reduced in these regions. These findings are consistent with an effect of low recombination in regions adjacent to centromeric or subtelomeric heterochromatin and add to our understanding of the processes generating widespread heterogeneity in genetic diversity and differentiation along the genome. By combining three different technologies, our results highlight the importance of adding a layer of information on genome structure that is inaccessible to each approach independently., (© 2017 Weissensteiner et al.; Published by Cold Spring Harbor Laboratory Press.)

There is concern that the stresses of inducing pluripotency may lead to deleterious DNA mutations in induced pluripotent stem cell (iPSC) lines, which would compromise their use for cell therapies. Here we report comparative genomic analysis of nine isogenic iPSC lines generated using three reprogramming methods: integrating retroviral vectors, non-integrating Sendai virus and synthetic mRNAs. We used whole-genome sequencing and de novo genome mapping to identify single-nucleotide variants, insertions and deletions, and structural variants. Our results show a moderate number of variants in the iPSCs that were not evident in the parental fibroblasts, which may result from reprogramming. There were only small differences in the total numbers and types of variants among different reprogramming methods. Most importantly, a thorough genomic analysis showed that the variants were generally benign. We conclude that the process of reprogramming is unlikely to introduce variants that would make the cells inappropriate for therapy.