Ordered restriction maps of Saccharomyces cerevisiae chromosomes constructed by optical mapping

Diagnostic and Prognostic/Therapeutic Significance of Comprehensive Analysis of Bone and Soft Tissue Tumors Using Optical Genome Mapping and Next-Generation Sequencing

Detecting somatic structural variants (SVs), copy number variants (CNVs), and mutations in bone and soft tissue tumors is essential for accurately diagnosing, treating, and prognosticating outcomes. Optical genome mapping (OGM) holds promise to yield useful data on SVs and CNVs but requires fresh or snap-frozen tissues. This study aimed to evaluate the clinical utility of data from OGM compared with current standard-of-care cytogenetic testing. We evaluated 60 consecutive specimens from bone and soft tissue tumors using OGM and karyotyping, fluorescence in situ hybridization, gene fusion assays, and deep next-generation sequencing. OGM accurately identified diagnostic SVs/CNVs previously detected by karyotyping and fluorescence in situ hybridization (specificity = 100%). OGM identified diagnostic and pathogenic SVs/CNVs (∼23% of cases) undetected by karyotyping (cryptic/submicroscopic). OGM allowed the detection and further characterization of complex structural rearrangements including chromoanagenesis (27% of cases) and complex 3- to 6-way translocations (15% of cases). In addition to identifying 321 SVs and CNVs among cases with chromoanagenesis events, OGM identified approximately 9 SVs and 12 CNVs per sample. A combination of OGM and deep next-generation sequencing data identified diagnostic, disease-associated, and pathogenic SVs, CNVs, and mutations in ∼98% of the cases. Our cohort contained the most extensive collection of bone and soft tissue tumors profiled by OGM. OGM had excellent concordance with standard-of-care cytogenetic testing, detecting and assigning high-resolution genome-wide genomic abnormalities with higher sensitivity than routine testing. This is the first and largest study to provide insights into the clinical utility of combined OGM and deep sequencing for the pathologic diagnosis and potential prognostication of bone and soft tissue tumors in routine clinical practice.

Molecular diagnosis and preimplantation genetic testing for chromosome 1q21.1 recurrent microduplication

As the development of molecular diagnostic methods, a large number of clinically relevant or disease-related copy number variations (CNVs) could be detected, and the demand for genetic counselling and clinical treatment is also increasing. For patients with pathogenic or likely pathogenic CNVs, preimplantation genetic testing (PGT) could provide a feasible path to prevent the inheritance of the genetic disorder in the offspring. In this study, we included a couple with 1q21.1 recurrent microduplication to conduct molecular diagnosis and PGT clinical application. The optical genome mapping (OGM) successfully verified the orientation and location of the microduplication, which further proved OGM as a promising approach for chromosomal anomalies detection with high resolutions. In PGT application, linkage-analysis-based PGT and high resolution PGT-A were simultaneously conducted for the pedigree and all the embryos. The results were consistent between linkage analysis and high resolution aneuploid analysis in the targeted region. One embryo that was absent of paternal 1q21.1q21.2 duplication was selected for further transplantation. This successful clinical practice in this study shed light for future molecular diagnosis and PGT application in tandem microduplications.

Sequence-Specific Minor Groove Binders in Labeling and Single-Molecule Analysis of DNA

The ability to address specific sequences within DNA is of tremendous interest in biotechnology and biomedicine. Various technologies have been established over the past few decades, such as nicking enzymes and methyltransferase-directed sequence-specific labeling, transcription activator-like effector nucleases (TALENs), the CRISPR-Cas9 system, and polyamides of heterocycles as sequence-specific DNA minor groove binders. Pyrrole-imidazole polyamides have been reported to recognize predetermined DNA sequences, and some successful attempts have demonstrated their potential in regulating gene expression. However, few studies on single-molecule labeling and analysis of DNA have been explored, particularly at single-targeting sites. In this study, we rationally designed and synthesized a set of functional minor groove binders, varying in structures, sequence information addressed, and methods of dye introduction. Their potential for sequence-specific labeling and single-molecule DNA analysis was evaluated through chromatographic and on-surface optical assays. First results indicated that, while they yielded excellent imaging output, the labeling specificity of the hairpin polyamides for single-molecule use was hindered by single-mismatch sites. To address this issue, and in an unprecedented approach, we devised a competitive binding strategy that utilizes ethidium bromide as a nonspecific binder to competitively block the mismatches, significantly enhancing the labeling specificity. These findings provide valuable insights into the use of hairpin polyamides as sequence-programmable and enzyme-free DNA labelers in the field of sequence-specific DNA recognition and modification.

Sequencing and Optical Genome Mapping for the Adventurous Chemist

This review provides a comprehensive overview of the chemistries and workflows of the sequencing methods that have been or are currently commercially available, providing a very brief historical introduction to each method. The main optical genome mapping approaches are introduced in the same manner, although only a subset of these are or have ever been commercially available. The review comes with a deck of slides containing all of the figures for ease of access and consultation.

ct-DNA compaction by nanoparticles formed by silica and gemini surfactants having hydroxyl group substituted spacers: In vitro, in vivo, and ex vivo gene uptake to cancer cells

Hybrid nanoparticles formed by Silica (SiO2) coated with cationic gemini surfactants with variable hydroxyl group substituted spacers, 12-4(OH)-12,2Br- and 12-4(OH)2-12,2Br- have shown a great extent of compaction of calf thymus DNA (ct-DNA) compared to conventional counterpart cationic surfactant, dodecyl trimethylammonium bromide (DTAB). Study shows not only the hydrophobicity of the spacer but also the hydrogen bonding interactions between the hydroxyl group substituted spacer and DNA have a great role in DNA compaction. 12-4(OH)2-12,2Br- is more efficient in compacting ct-DNA compared to 12-4(OH)-12,2Br- due to the stronger binding of the former with ct-DNA than the latter. While 12-4(OH)-12,2Br- makes 50 % ct-DNA compaction at its 0.63 μM concentration in the presence of SiO2 nanoparticles, the same % of compaction can be achieved at a concentration as low as 0.25 μM of 12-4(OH)2-12,2Br-. However, DTAB makes 50 % ct-DNA compaction at a concentration as high as 7.00 μM under the same condition. Therefore, the present systems address the very common challenge, i.e., cytotoxicity due to cationic surfactants. The system of 12-4(OH)2-12,2Br- coated SiO2 nanoparticles displays the maximum cell viability (≥90 %), causing the least cell death in the mouse fibroblast cells (NIH3T3) cell lines compared to the cell viability of ≤80 % for DTAB. 12-4(OH)2-12,2Br- coated SiO2 nanoparticles system has presented excellent in vitro cellular uptake of genes on mouse mammary gland adenocarcinoma (4T1) cells after incubating for 3 h and 6 h. In vivo study shows that 12-4(OH)2-12,2Br- coated SiO2 nanoparticles system takes the highest amount of ct-DNA in cells and tumors in a time-dependent manner. The ex vivo studies using different organs of the mice demonstrate that the tumor sites in the breast of the mice are most affected by these formulations. Cytotoxicity assays and cellular uptake studies suggest that the present systems can be used for potential applications for gene delivery and oncological therapies.

Cytogenetics in the management of hematological malignancies: An overview of alternative technologies for cytogenetic characterization

Genomic characterization is an essential part of the clinical management of hematological malignancies for diagnostic, prognostic and therapeutic purposes. Although CBA and FISH are still the gold standard in hematology for the detection of CNA and SV, some alternative technologies are intended to complement their deficiencies or even replace them in the more or less near future. In this article, we provide a technological overview of these alternatives. CMA is the historical and well established technique for the high-resolution detection of CNA. For SV detection, there are emerging techniques based on the study of chromatin conformation and more established ones such as RTMLPA for the detection of fusion transcripts and RNA-seq to reveal the molecular consequences of SV. Comprehensive techniques that detect both CNA and SV are the most interesting because they provide all the information in a single examination. Among these, OGM is a promising emerging higher-solution technique that offers a complete solution at a contained cost, at the expense of a relatively low throughput per machine. WGS remains the most adaptable solution, with long-read approaches enabling very high-resolution detection of CAs, but requiring a heavy bioinformatics installation and at a still high cost. However, the development of high-resolution genome-wide detection techniques for CAs allows for a much better description of chromoanagenesis. Therefore, we have included in this review an update on the various existing mechanisms and their consequences and implications, especially prognostic, in hematological malignancies.

Scanning Electron Microscopy Imaging of Large DNA Molecules Using a Metal-Free Electro-Stain Composed of DNA-Binding Proteins and Synthetic Polymers

This paper presents the first scanning electron microscopy (SEM)-based DNA imaging in biological samples. This novel approach incorporates a metal-free electro-stain reagent, formulated by combining DNA-binding proteins and synthetic polymers to enhance the visibility of 2-nm-thick DNA under SEM. Notably, DNA molecules stain with proteins and polymers appear as dark lines under SEM. The resulting DNA images exhibit a thickness of 15.0±4.0 nm. As SEM is the primary platform, it integrates seamlessly with various chemically functionalized large surfaces with the aid of microfluidic devices. The approach allows high-resolution imaging of various DNA structures including linear, circular, single-stranded DNA and RNA, originating from nuclear and mitochondrial genomes. Furthermore, quantum dots are successfully visualized as bright labels that are sequence-specifically incorporated into DNA molecules, which highlights the potential for SEM-based optical DNA mapping. In conclusion, DNA imaging using SEM with the novel electro-stain offers electron microscopic resolution with the ease of optical microscopy.

Optical Genome Mapping for Applications in Repeat Expansion Disorders

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.

Sequencing and characterizing short tandem repeats in the human genome

Short tandem repeats (STRs) are highly polymorphic sequences throughout the human genome that are composed of repeated copies of a 1-6-bp motif. Over 1 million variable STR loci are known, some of which regulate gene expression and influence complex traits, such as height. Moreover, variants in at least 60 STR loci cause genetic disorders, including Huntington disease and fragile X syndrome. Accurately identifying and genotyping STR variants is challenging, in particular mapping short reads to repetitive regions and inferring expanded repeat lengths. Recent advances in sequencing technology and computational tools for STR genotyping from sequencing data promise to help overcome this challenge and solve genetically unresolved cases and the 'missing heritability' of polygenic traits. Here, we compare STR genotyping methods, analytical tools and their applications to understand the effect of STR variation on health and disease. We identify emergent opportunities to refine genotyping and quality-control approaches as well as to integrate STRs into variant-calling workflows and large cohort analyses.

Methods of Detection and Mechanisms of Origin of Complex Structural Genome Variations

Based on classical karyotyping, structural genome variations (SVs) have generally been considered to be either "simple" (with one or two breakpoints) or "complex" (with more than two breakpoints). Studying the breakpoints of SVs at nucleotide resolution revealed additional, subtle structural variations, such that even "simple" SVs turned out to be "complex." Genome-wide sequencing methods, such as fosmid and paired-end mapping, short-read and long-read whole genome sequencing, and single-molecule optical mapping, also indicated that the number of SVs per individual was considerably larger than expected from karyotyping and high-resolution chromosomal array-based studies. Interestingly, SVs were detected in studies of cohorts of individuals without clinical phenotypes. The common denominator of all SVs appears to be a failure to accurately repair DNA double-strand breaks (DSBs) or to halt cell cycle progression if DSBs persist. This review discusses the various DSB response mechanisms during the mitotic cell cycle and during meiosis and their regulation. Emphasis is given to the molecular mechanisms involved in the formation of translocations, deletions, duplications, and inversions during or shortly after meiosis I. Recently, CRISPR-Cas9 studies have provided unexpected insights into the formation of translocations and chromothripsis by both breakage-fusion-bridge and micronucleus-dependent mechanisms.

Identification of cryptic balanced translocations in couples with unexplained recurrent pregnancy loss based upon embryonic PGT-A results

PURPOSE

The goal of this study is to determine whether any balanced translocation (BT) had been missed by previous karyotyping in patients with unexplained recurrent pregnancy loss (uRPL).

METHODS

This case series included 48 uRPL-affected couples with normal karyotypes. The embryos from these couples have all undergone preimplantation testing for aneuploidies (PGT-A). Based on the PGT-A's results, 48 couples could be categorized into two groups: 17 couples whose multiple embryos were detected with similar structural variations (SVs, segmental/complete) and 31 couples without such findings but who did not develop any euploid embryo despite at least three high-quality blastocysts being tested. The peripheral blood sample of each partner was then collected for mate-pair sequencing (MPseq) to determine whether any of them were BT carriers.

RESULTS

MPseq analyses identified 13 BTs in the 17 couples whose multiple embryos had similar SVs detected (13/17, 76.47%) and three BTs in the 31 couples without euploid embryo obtained (3/31, 9.7%). Among the 16 MPseq-identified BTs, six were missed due to the limited resolution of G-banding karyotyping analysis, and the rest were mostly owing to the similar banding patterns and/or comparable sizes shared by the two segments exchanged.

CONCLUSION

A normal karyotype does not eliminate the possibility of carrying BT for couples with uRPL. The use of PGT-A allows us to perceive the "carrier couples" missed by karyotyping analysis, providing an increased risk of finding cryptic BTs if similar SVs are always detected on two chromosomes among multiple embryos. Nonetheless, certain carriers with translocated segments of sub-resolution may still go unnoticed.

Optical Genome Mapping: A Machine-Based Platform in Cytogenomics

Optical genome mapping (OGM) has generated excitement following decades of research and development. Now, commercially available technical platforms have been used to compare various other cytogenetic and cytogenomic technologies, including karyotype, microarrays, and DNA sequencing, with impressive results. In this chapter, using OGM as a case study, we advocate for a new trend in future cytogenomics, emphasizing the power of machine automation to deliver higher-quality cytogenomic data. By briefly discussing OGM, along with its major advantages and limitations, we underscore the importance of karyotype-based genomic research, from both a theoretical framework and a new technology perspective. We also call for the encouragement of further technological platform development for the future of cytogenetics and cytogenomics.

Design of optimal labeling patterns for optical genome mapping via information theory

MOTIVATION

Optical genome mapping (OGM) is a technique that extracts partial genomic information from optically imaged and linearized DNA fragments containing fluorescently labeled short sequence patterns. This information can be used for various genomic analyses and applications, such as the detection of structural variations and copy-number variations, epigenomic profiling, and microbial species identification. Currently, the choice of labeled patterns is based on the available biochemical methods and is not necessarily optimized for the application.

RESULTS

In this work, we develop a model of OGM based on information theory, which enables the design of optimal labeling patterns for specific applications and target organism genomes. We validated the model through experimental OGM on human DNA and simulations on bacterial DNA. Our model predicts up to 10-fold improved accuracy by optimal choice of labeling patterns, which may guide future development of OGM biochemical labeling methods and significantly improve its accuracy and yield for applications such as epigenomic profiling and cultivation-free pathogen identification in clinical samples.

AVAILABILITY AND IMPLEMENTATION

https://github.com/yevgenin/PatternCode.

Telomere-to-telomere assembly of diploid chromosomes with Verkko

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 high-fidelity 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 by integrating 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.

Optical Genome Mapping for Oncology Applications

Optical genome mapping (OGM) is a next-generation cytogenomic technology that has the potential to replace standard-of-care technologies used in the genetic workup of various malignancies. The ability to detect various classes of structural variations that include copy number variations, deletions, duplications, balanced and unbalanced events (insertions, inversions, and translocation) and complex genomic rearrangements in a single assay and analysis demonstrates the utility of the technology in tumor research and clinical application. Herein, we provide the methodological details for performing OGM and pre- and post-analytical quality control (QC) checks and describe critical steps that should be performed with caution, probable causes for specific QC failures, and potential method modifications that could be implemented as part of troubleshooting. The protocol description and troubleshooting guide should help new and current users of the technology to improve or troubleshoot the problems (if any) in their workflow. © 2023 Wiley Periodicals LLC. Basic Protocol: Optical genome mapping.

Microfluidic delivery of cutting enzymes for fragmentation of surface-adsorbed DNA molecules

We describe a method for fragmenting, in-situ, surface-adsorbed and immobilized DNAs on polymethylmethacrylate(PMMA)-coated silicon substrates using microfluidic delivery of the cutting enzyme DNase I. Soft lithography is used to produce silicone elastomer (Sylgard 184) gratings which form microfluidic channels for delivery of the enzyme. Bovine serum albumin (BSA) is used to reduce DNase I adsorption to the walls of the microchannels and enable diffusion of the cutting enzyme to a distance of 10mm. Due to the DNAs being immobilized, the fragment order is maintained on the surface. Possible methods of preserving the order for application to sequencing are discussed.

Genomics in the long-read sequencing era

Long-read sequencing (LRS) technologies have provided extremely powerful tools to explore genomes. While in the early years these methods suffered technical limitations, they have recently made significant progress in terms of read length, throughput, and accuracy and bioinformatics tools have strongly improved. Here, we aim to review the current status of LRS technologies, the development of novel methods, and the impact on genomics research. We will explore the most impactful recent findings made possible by these technologies focusing on high-resolution sequencing of genomes and transcriptomes and the direct detection of DNA and RNA modifications. We will also discuss how LRS methods promise a more comprehensive understanding of human genetic variation, transcriptomics, and epigenetics for the coming years.

Role of Gemini Surfactants with Variable Spacers and SiO2 Nanoparticles in ct-DNA Compaction and Applications toward In Vitro/In Vivo Gene Delivery

Compaction of calf thymus DNA (ct-DNA) by two cationic gemini surfactants, 12-4-12 and 12-8-12, in the absence and presence of negatively charged SiO₂ nanoparticles (NPs) (∼100 nm) has been explored using various techniques. 12-8-12 having a longer hydrophobic spacer induces a greater extent of ct-DNA compaction than 12-4-12, which becomes more efficient with SiO₂ NPs. While 50% ct-DNA compaction in the presence of SiO₂ NPs occurs at ∼77 nM of 12-8-12 and ∼130 nM of 12-4-12, but a conventional counterpart surfactant, DTAB, does it at its concentration as high as ∼7 μM. Time-resolved fluorescence anisotropy measurements show changes in the rotational dynamics of a fluorescent probe, DAPI, and helix segments in the condensed DNA. Fluorescence lifetime data and ethidium bromide exclusion assays reveal the binding sites of surfactants to ct-DNA. 12-8-12 with SiO₂ NPs has shown the highest cell viability (≥90%) and least cell death in the human embryonic kidney (HEK) 293 cell lines in contrast to the cell viability of ≤80% for DTAB. These results show that 12-8-12 with SiO₂ NPs has the highest time and dose-dependent cytotoxicity compared to 12-8-12 and 12-4-12 in the murine breast cancer 4T1 cell line. Fluorescence microscopy and flow cytometry are performed for in vitro cellular uptake of YOYO-1-labeled ct-DNA with surfactants and SiO₂ NPs using 4T1 cells after 3 and 6 h incubations. The in vivo tumor accumulation studies are carried out using a real-time in vivo imaging system after intravenous injection of the samples into 4T1 tumor-bearing mice. 12-8-12 with SiO₂ has delivered the highest amount of ct-DNA in cells and tumors in a time-dependent manner. Thus, the application of a gemini surfactant with a hydrophobic spacer and SiO₂ NPs in compacting and delivering ct-DNA to the tumor is proven, warranting its further exploration in nucleic acid therapy for cancer treatment.

Sequence-specific DNA labelling for fluorescence microscopy

The preservation of nucleus structure during microscopy imaging is a top priority for understanding chromatin organization, genome dynamics, and gene expression regulation. In this review, we summarize the sequence-specific DNA labelling methods that can be used for imaging in fixed and/or living cells without harsh treatment and DNA denaturation: (i) hairpin polyamides, (ii) triplex-forming oligonucleotides, (iii) dCas9 proteins, (iv) transcription activator-like effectors (TALEs) and (v) DNA methyltransferases (MTases). All these techniques are capable of identifying repetitive DNA loci and robust probes are available for telomeres and centromeres, but visualizing single-copy sequences is still challenging. In our futuristic vision, we see gradual replacement of the historically important fluorescence in situ hybridization (FISH) by less invasive and non-destructive methods compatible with live cell imaging. Combined with super-resolution fluorescence microscopy, these methods will open the possibility to look into unperturbed structure and dynamics of chromatin in living cells, tissues and whole organisms.

Optical Genome Mapping for a Patient with a Congenital Disorder and Chromosomal Translocation

We performed optical genome mapping (OGM), a newly developed cytogenetic technique, for a patient with a disorder of sex development (DSD) and a 46,XX,t(9;11)(p22;p13) karyotype. The results of OGM were validated using other methods. OGM detected a 9;11 reciprocal translocation and successfully mapped its breakpoints to small regions of 0.9-12.3 kb. OGM identified 46 additional small structural variants, only three of which were detected by array-based comparative genomic hybridization. OGM suggested the presence of complex rearrangements on chromosome 10; however, these variants appeared to be artifacts. The 9;11 translocation was unlikely to be associated with DSD, while the pathogenicity of the other structural variants remained unknown. These results indicate that OGM is a powerful tool for detecting and characterizing chromosomal structural variations, although the current methods of OGM data analyses need to be improved.

Population genomic analysis provides evidence of the past success and future potential of South China tiger captive conservation

BACKGROUND

Among six extant tiger subspecies, the South China tiger (Panthera tigris amoyensis) once was widely distributed but is now the rarest one and extinct in the wild. All living South China tigers are descendants of only two male and four female wild-caught tigers and they survive solely in zoos after 60 years of effective conservation efforts. Inbreeding depression and hybridization with other tiger subspecies were believed to have occurred within the small, captive South China tiger population. It is therefore urgently needed to examine the genomic landscape of existing genetic variation among the South China tigers.

RESULTS

In this study, we assembled a high-quality chromosome-level genome using long-read sequences and re-sequenced 29 high-depth genomes of the South China tigers. By combining and comparing our data with the other 40 genomes of six tiger subspecies, we identified two significantly differentiated genomic lineages among the South China tigers, which harbored some rare genetic variants introgressed from other tiger subspecies and thus maintained a moderate genetic diversity. We noticed that the South China tiger had higher F_ROH values for longer runs of homozygosity (ROH > 1 Mb), an indication of recent inbreeding/founder events. We also observed that the South China tiger had the least frequent homozygous genotypes of both high- and moderate-impact deleterious mutations, and lower mutation loads than both Amur and Sumatran tigers. Altogether, our analyses indicated an effective genetic purging of deleterious mutations in homozygous states from the South China tiger, following its population contraction with a controlled increase in inbreeding based on its pedigree records.

CONCLUSIONS

The identification of two unique founder/genomic lineages coupled with active genetic purging of deleterious mutations in homozygous states and the genomic resources generated in our study pave the way for a genomics-informed conservation, following the real-time monitoring and rational exchange of reproductive South China tigers among zoos.

Resequencing of a Pekin duck breeding population provides insights into the genomic response to short-term artificial selection

BACKGROUND

Short-term, intense artificial selection drives fast phenotypic changes in domestic animals and leaves imprints on their genomes. However, the genetic basis of this selection response is poorly understood. To better address this, we employed the Pekin duck Z2 pure line, in which the breast muscle weight was increased nearly 3-fold after 10 generations of breeding. We denovo assembled a high-quality reference genome of a female Pekin duck of this line (GCA_003850225.1) and identified 8.60 million genetic variants in 119 individuals among 10 generations of the breeding population.

RESULTS

We identified 53 selected regions between the first and tenth generations, and 93.8% of the identified variations were enriched in regulatory and noncoding regions. Integrating the selection signatures and genome-wide association approach, we found that 2 regions covering 0.36 Mb containing UTP25 and FBRSL1 were most likely to contribute to breast muscle weight improvement. The major allele frequencies of these 2 loci increased gradually with each generation following the same trend. Additionally, we found that a copy number variation region containing the entire EXOC4 gene could explain 1.9% of the variance in breast muscle weight, indicating that the nervous system may play a role in economic trait improvement.

CONCLUSIONS

Our study not only provides insights into genomic dynamics under intense artificial selection but also provides resources for genomics-enabled improvements in duck breeding.

Optical Genome Mapping for Cytogenetic Diagnostics in AML

The classification and risk stratification of acute myeloid leukemia (AML) is based on reliable genetic diagnostics. A broad and expanding variety of relevant aberrations are structural variants beyond single-nucleotide variants. Optical Genome Mapping is an unbiased, genome-wide, amplification-free method for the detection of structural variants. In this review, the current knowledge of Optical Genome Mapping (OGM) with regard to diagnostics in hematological malignancies in general, and AML in specific, is summarized. Furthermore, this review focuses on the ability of OGM to expand the use of cytogenetic diagnostics in AML and perhaps even replace older techniques such as chromosomal-banding analysis, fluorescence in situ hybridization, or copy number variation microarrays. Finally, OGM is compared to amplification-based techniques and a brief outlook for future directions is given.

Era of gapless plant genomes: innovations in sequencing and mapping technologies revolutionize genomics and breeding

Whole-genome sequencing and assembly have revolutionized plant genetics and molecular biology over the last two decades. However, significant shortcomings in first- and second-generation technology resulted in imperfect reference genomes: numerous and large gaps of low quality or undeterminable sequence in areas of highly repetitive DNA along with limited chromosomal phasing restricted the ability of researchers to characterize regulatory noncoding elements and genic regions that underwent recent duplication events. Recently, advances in long-read sequencing have resulted in the first gapless, telomere-to-telomere (T2T) assemblies of plant genomes. This leap forward has the potential to increase the speed and confidence of genomics and molecular experimentation while reducing costs for the research community.

Testing for Facioscapulohumeral Muscular Dystrophy with Optical Genome Mapping

The introduction of optical genome mapping has improved time constraints and a lack of specificity from previous methodologies when performing genome-wide analyses of samples. Optical genome mapping allows for the detection of structural variations, aberrations, and functionality traits from a single stained molecule of DNA. Though the preparation time is increased compared to previously utilized visualization techniques, optical genome mapping significantly reduces the time needed for analysis. Specifically, individual disease pipelines have been developed to rapidly analyze prepared samples. One of these diseases, Facioscapulohumeral Muscular Dystrophy (FSHD), is detected through quantification of the D4Z4 repeat array on chromosome 4q35. Optical genome mapping, with the ability to enumerate the repeats of the D4Z4 array, has demonstrated the capability to precisely diagnose FSHD. In this protocol, the preparation of samples and subsequent loading and analysis in an optical genome mapping system is discussed for the detection and analysis of FSHD. These methods should prove highly useful in FSHD analyses and beyond with the development of further disease analysis pipelines within the instrument. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Genomic DNA isolation, labeling, and staining Basic Protocol 2: Mapping and analysis with the Bionano Saphyr® system.

Detection of Transgene Location in the CYP2A13/2B6/2F1-transgenic Mouse Model using Optical Genome Mapping Technology

Most transgenic mouse models are generated through random integration of the transgene. The location of the transgene provides valuable information for assessing potential effects of the transgenesis on the host and for designing genotyping protocols that can amplify across the integration site, but it is challenging to identify. Here, we report the successful utility of optical genome mapping technology to identify the transgene insertion site in a CYP2A13/2B6/2F1-transgenic mouse model, which produces three human cytochrome P450 (P450) enzymes (CYP2A13, CYP2B6, and CYP2F1) that are encoded by neighboring genes on human chromosome 19. These enzymes metabolize many drugs, respiratory toxicants, and chemical carcinogens. Initial efforts to identify candidate insertion sites by whole genome sequencing was unsuccessful, apparently because the transgene is located in a region of the mouse genome that contains highly repetitive sequences. Subsequent utility of the optical genome mapping approach, which compares genome-wide marker distribution between the transgenic mouse genome and a reference mouse (GRCm38) or human (GRCh38) genome, localized the insertion site to mouse chromosome 14, between two marker positions at 4451324 base pair and 4485032 base pair. A transgene-mouse genome junction sequence was further identified through long-polymerase chain reaction amplification and DNA sequencing at GRCm38 Chr.14:4484726. The transgene insertion (∼2.4 megabase pair) contained 5-7 copies of the human transgenes, which replaced a 26.9-33.4 kilobase pair mouse genomic region, including exons 1-4 of Gm3182, a predicted and highly redundant gene. Finally, the sequencing results enabled the design of a new genotyping protocol that can distinguish between hemizygous and homozygous CYP2A13/2B6/2F1-transgenic mice. SIGNIFICANCE STATEMENT: This study characterizes the genomic structure of, and provides a new genotyping method for, a transgenic mouse model that expresses three human P450 enzymes, CYP2A13, CYP2B6, and CYP2F1, that are important in xenobiotic metabolism and toxicity. The demonstrated success in applying the optical genome mapping technology for identification of transgene insertion sites should encourage others to do the same for other transgenic models generated through random integration, including most of the currently available human P450 transgenic mouse models.

Flow Sorting-Assisted Optical Mapping

Optical mapping-a technique that visualizes short sequence motives along DNA molecules of hundred kilobases to megabase in size-has found an important place in genome research. It is widely used to facilitate genome sequence assemblies and analyses of genome structural variations. Application of the technique is conditional on availability of highly pure ultra-long high-molecular-weight DNA (uHMW DNA), which is challenging to achieve in plants due to the presence of the cell wall, chloroplasts, and secondary metabolites, just as a high content of polysaccharides and DNA nucleases in some species. These obstacles can be overcome by employment of flow cytometry, enabling a fast and highly efficient purification of cell nuclei or metaphase chromosomes, which are afterward embedded in agarose plugs and used to isolate the uHMW DNA in situ. Here, we provide a detailed protocol for the flow sorting-assisted uHMW DNA preparation that has been successfully used to construct whole-genome as well as chromosomal optical maps for 20 plant species from several plant families.

Evaluation of optical genome mapping for detecting chromosomal translocation in clinical cytogenetics

Background

Balanced reciprocal translocation is one of the most common chromosomal abnormalities in humans that may lead to infertility, recurrent pregnancy loss, or having children with physical or mental abnormalities. Karyotyping and FISH are traditional detection approaches with a low resolution. Bionano optical genome mapping (OGM) developed in recent years can be used to analyze chromosomal abnormalities at a higher resolution, providing the possibility of more in‐depth analyses of balanced chromosome translocations.

Methods

To evaluate the feasibility of OGM to detect chromosome balanced translocations, 10 genetic outpatients were collected and detected simultaneously by karyotype analysis, FISH, CNV‐seq, and Bionano OGM in this study.

Results

The results showed that the karyotypes of the patients were detected by karyotype analysis, FISH, and Bionano OGM, but one patient with karyotype t(Y,19) was not correctly detected by OGM. There were not find any chromosome abnormality by CNV‐seq. More importantly, OGM allowed the location of the mutation to the gene level, which is important for aiding diagnoses, compared to karyotype analysis, and FISH.

Conclusions

This study shows that OGM can be a high adjunctive diagnostic method for detecting balanced chromosome translocations, but the accuracy and precision of OGM detecting mutations need to be gradually improved in telomere and centromere regions.

RegScaf: a regression approach to scaffolding

MOTIVATION

Crucial to the correctness of a genome assembly is the accuracy of the underlying scaffolds that specify the orders and orientations of contigs together with the gap distances between contigs. The current methods construct scaffolds based on the alignments of 'linking' reads against contigs. We found that some 'optimal' alignments are mistaken due to factors such as the contig boundary effect, particularly in the presence of repeats. Occasionally, the incorrect alignments can even overwhelm the correct ones. The detection of the incorrect linking information is challenging in any existing methods.

RESULTS

In this study, we present a novel scaffolding method RegScaf. It first examines the distribution of distances between contigs from read alignment by the kernel density. When multiple modes are shown in a density, orientation-supported links are grouped into clusters, each of which defines a linking distance corresponding to a mode. The linear model parameterizes contigs by their positions on the genome; then each linking distance between a pair of contigs is taken as an observation on the difference of their positions. The parameters are estimated by minimizing a global loss function, which is a version of trimmed sum of squares. The least trimmed squares estimate has such a high breakdown value that it can automatically remove the mistaken linking distances. The results on both synthetic and real datasets demonstrate that RegScaf outperforms some popular scaffolders, especially in the accuracy of gap estimates by substantially reducing extremely abnormal errors. Its strength in resolving repeat regions is exemplified by a real case. Its adaptability to large genomes and TGS long reads is validated as well.

AVAILABILITY AND IMPLEMENTATION

RegScaf is publicly available at https://github.com/lemontealala/RegScaf.git.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Clinical whole-genome sequencing in cancer diagnosis

Characterizing the genomic landscape of cancers is a routine part of clinical care that began with the discovery of the Philadelphia chromosome and has since coevolved with genomic technologies. Genomic analysis of tumors at the nucleotide level using DNA sequencing has revolutionized the understanding of cancer biology and identified new molecular drivers of disease that have led to therapeutic advances and improved patient outcomes. However, the application of next-generation sequencing in the clinical laboratory has generally been limited until very recently to targeted analysis of selected genes. Recent technological innovations and reductions in sequencing costs are now able to deliver the long-promised goal of tumor whole-genome sequencing as a practical clinical assay.

Optical Genome Mapping as a Diagnostic Tool in Pediatric Acute Myeloid Leukemia

Pediatric AML is characterized by numerous genetic aberrations (chromosomal translocations, deletions, insertions) impacting its classification for risk of treatment failure. Aberrations are described by classical cytogenetic procedures (karyotyping, FISH), which harbor limitations (low resolution, need for cell cultivation, cost-intensiveness, experienced staff required). Optical Genome Mapping (OGM) is an emerging chip-based DNA technique combining high resolution (~500 bp) with a relatively short turnaround time. Twenty-four pediatric patients with AML, bi-lineage leukemia, and mixed-phenotype acute leukemia were analyzed by OGM, and the results were compared with cytogenetics. Results were discrepant in 17/24 (70%) cases, including 32 previously unknown alterations called by OGM only. One newly detected deletion and two translocations were validated by primer walking, breakpoint-spanning PCR, and DNA sequencing. As an added benefit, in two cases, OGM identified a new minimal residual disease (MRD) marker. Comparing impact on risk stratification in de novo AML, 19/20 (95%) cases had concordant results while only OGM unraveled another high-risk aberration. Thus, OGM considerably expands the methodological spectrum to optimize the diagnosis of pediatric AML via the identification of new aberrations. Results will contribute to a better understanding of leukemogenesis in pediatric AML. In addition, aberrations identified by OGM may provide markers for MRD monitoring.

Progress in Methods for Copy Number Variation Profiling

Copy number variations (CNVs) are the predominant class of structural genomic variations involved in the processes of evolutionary adaptation, genomic disorders, and disease progression. Compared with single-nucleotide variants, there have been challenges associated with the detection of CNVs owing to their diverse sizes. However, the field has seen significant progress in the past 20–30 years. This has been made possible due to the rapid development of molecular diagnostic methods which ensure a more detailed view of the genome structure, further complemented by recent advances in computational methods. Here, we review the major approaches that have been used to routinely detect CNVs, ranging from cytogenetics to the latest sequencing technologies, and then cover their specific features.

From Genome Sequencing to CRISPR-Based Genome Editing for Climate-Resilient Forest Trees

Due to the economic and ecological importance of forest trees, modern breeding and genetic manipulation of forest trees have become increasingly prevalent. The CRISPR-based technology provides a versatile, powerful, and widely accepted tool for analyzing gene function and precise genetic modification in virtually any species but remains largely unexplored in forest species. Rapidly accumulating genetic and genomic resources for forest trees enabled the identification of numerous genes and biological processes that are associated with important traits such as wood quality, drought, or pest resistance, facilitating the selection of suitable gene editing targets. Here, we introduce and discuss the latest progress, opportunities, and challenges of genome sequencing and editing for improving forest sustainability.

Applications of Optical Mapping for Plant Genome Assembly and Structural Variation Detection

Optical mapping plays an important role in plant genomics, particularly in plant genome assembly and large-scale structural variation detection. While DNA sequencing provides base-by-base nucleotide information, optical mapping shows the physical locations of selected enzyme restriction sites in a genome. The long single-molecule maps produced by optical mapping make it a useful auxiliary technique to DNA sequencing, which generally cannot span large and complex genomic regions. Although optical mapping, therefore, offers unique advantages to researchers, there are few dedicated tools to assist in optical mapping analyses. In this chapter, we present runBNG2, a successor of runBNG to help optical-mapping data analysis for diverse datasets.

Combining dense and sparse labeling in optical DNA mapping

Optical DNA mapping (ODM) is based on fluorescent labeling, stretching and imaging of single DNA molecules to obtain sequence-specific fluorescence profiles, DNA barcodes. These barcodes can be mapped to theoretical counterparts obtained from DNA reference sequences, which in turn allow for DNA identification in complex samples and for detecting structural changes in individual DNA molecules. There are several types of DNA labeling schemes for ODM and for each labeling type one or several types of match scoring methods are used. By combining the information from multiple labeling schemes one can potentially improve mapping confidence; however, combining match scores from different labeling assays has not been implemented yet. In this study, we introduce two theoretical methods for dealing with analysis of DNA molecules with multiple label types. In our first method, we convert the alignment scores, given as output from the different assays, into p-values using carefully crafted null models. We then combine the p-values for different label types using standard methods to obtain a combined match score and an associated combined p-value. In the second method, we use a block bootstrap approach to check for the uniqueness of a match to a database for all barcodes matching with a combined p-value below a predefined threshold. For obtaining experimental dual-labeled DNA barcodes, we introduce a novel assay where we cut plasmid DNA molecules from bacteria with restriction enzymes and the cut sites serve as sequence-specific markers, which together with barcodes obtained using the established competitive binding labeling method, form a dual-labeled barcode. All experimental data in this study originates from this assay, but we point out that our theoretical framework can be used to combine data from all kinds of available optical DNA mapping assays. We test our multiple labeling frameworks on barcodes from two different plasmids and synthetically generated barcodes (combined competitive-binding- and nick-labeling). It is demonstrated that by simultaneously using the information from all label types, we can substantially increase the significance when we match experimental barcodes to a database consisting of theoretical barcodes for all sequenced plasmids.

DNA bridges: A novel platform for single-molecule sequencing and other DNA-protein interaction applications

DNA molecular combing is a technique that stretches thousands of long individual DNA molecules (up to 10 Mbp) into a parallel configuration on surface. It has previously been proposed to sequence these molecules by synthesis. However, this approach poses two critical challenges: 1-Combed DNA molecules are overstretched and therefore a nonoptimal substrate for polymerase extension. 2-The combing surface sterically impedes full enzymatic access to the DNA backbone. Here, we introduce a novel approach that attaches thousands of molecules to a removable surface, with a tunable stretching factor. Next, we dissolve portions of the surface, leaving the DNA molecules suspended as 'bridges'. We demonstrate that the suspended molecules are enzymatically accessible, and we have used an enzyme to incorporate labeled nucleotides, as predicted by the specific molecular sequence. Our results suggest that this novel platform is a promising candidate to achieve high-throughput sequencing of Mbp-long molecules, which could have additional genomic applications, such as the study of other protein-DNA interactions.

Filling gaps of genome scaffolds via probabilistic searching optical maps against assembly graph

BACKGROUND

Optical maps record locations of specific enzyme recognition sites within long genome fragments. This long-distance information enables aligning genome assembly contigs onto optical maps and ordering contigs into scaffolds. The generated scaffolds, however, often contain a large amount of gaps. To fill these gaps, a feasible way is to search genome assembly graph for the best-matching contig paths that connect boundary contigs of gaps. The combination of searching and evaluation procedures might be "searching followed by evaluation", which is infeasible for long gaps, or "searching by evaluation", which heavily relies on heuristics and thus usually yields unreliable contig paths.

RESULTS

We here report an accurate and efficient approach to filling gaps of genome scaffolds with aids of optical maps. Using simulated data from 12 species and real data from 3 species, we demonstrate the successful application of our approach in gap filling with improved accuracy and completeness of genome scaffolds.

CONCLUSION

Our approach applies a sequential Bayesian updating technique to measure the similarity between optical maps and candidate contig paths. Using this similarity to guide path searching, our approach achieves higher accuracy than the existing "searching by evaluation" strategy that relies on heuristics. Furthermore, unlike the "searching followed by evaluation" strategy enumerating all possible paths, our approach prunes the unlikely sub-paths and extends the highly-probable ones only, thus significantly increasing searching efficiency.

Signal-based optical map alignment

In genomics, optical mapping technology provides long-range contiguity information to improve genome sequence assemblies and detect structural variation. Originally a laborious manual process, Bionano Genomics platforms now offer high-throughput, automated optical mapping based on chips packed with nanochannels through which unwound DNA is guided and the fluorescent DNA backbone and specific restriction sites are recorded. Although the raw image data obtained is of high quality, the processing and assembly software accompanying the platforms is closed source and does not seem to make full use of data, labeling approximately half of the measured signals as unusable. Here we introduce two new software tools, independent of Bionano Genomics software, to extract and process molecules from raw images (OptiScan) and to perform molecule-to-molecule and molecule-to-reference alignments using a novel signal-based approach (OptiMap). We demonstrate that the molecules detected by OptiScan can yield better assemblies, and that the approach taken by OptiMap results in higher use of molecules from the raw data. These tools lay the foundation for a suite of open-source methods to process and analyze high-throughput optical mapping data. The Python implementations of the OptiTools are publicly available through http://www.bif.wur.nl/.

Expanding the conservation genomics toolbox: Incorporating structural variants to enhance genomic studies for species of conservation concern

Structural variants (SVs) are large rearrangements (>50 bp) within the genome that impact gene function and the content and structure of chromosomes. As a result, SVs are a significant source of functional genomic variation, that is, variation at genomic regions underpinning phenotype differences, that can have large effects on individual and population fitness. While there are increasing opportunities to investigate functional genomic variation in threatened species via single nucleotide polymorphism (SNP) data sets, SVs remain understudied despite their potential influence on fitness traits of conservation interest. In this future-focused Opinion, we contend that characterizing SVs offers the conservation genomics community an exciting opportunity to complement SNP-based approaches to enhance species recovery. We also leverage the existing literature-predominantly in human health, agriculture and ecoevolutionary biology-to identify approaches for readily characterizing SVs and consider how integrating these into the conservation genomics toolbox may transform the way we manage some of the world's most threatened species.

Optical genome mapping enables constitutional chromosomal aberration detection

Chromosomal aberrations including structural variations (SVs) are a major cause of human genetic diseases. Their detection in clinical routine still relies on standard cytogenetics. Drawbacks of these tests are a very low resolution (karyotyping) and the inability to detect balanced SVs or indicate the genomic localization and orientation of duplicated segments or insertions (copy number variant [CNV] microarrays). Here, we investigated the ability of optical genome mapping (OGM) to detect known constitutional chromosomal aberrations. Ultra-high-molecular-weight DNA was isolated from 85 blood or cultured cells and processed via OGM. A de novo genome assembly was performed followed by structural variant and CNV calling and annotation, and results were compared to known aberrations from standard-of-care tests (karyotype, FISH, and/or CNV microarray). In total, we analyzed 99 chromosomal aberrations, including seven aneuploidies, 19 deletions, 20 duplications, 34 translocations, six inversions, two insertions, six isochromosomes, one ring chromosome, and four complex rearrangements. Several of these variants encompass complex regions of the human genome involved in repeat-mediated microdeletion/microduplication syndromes. High-resolution OGM reached 100% concordance compared to standard assays for all aberrations with non-centromeric breakpoints. This proof-of-principle study demonstrates the ability of OGM to detect nearly all types of chromosomal aberrations. We also suggest suited filtering strategies to prioritize clinically relevant aberrations and discuss future improvements. These results highlight the potential for OGM to provide a cost-effective and easy-to-use alternative that would allow comprehensive detection of chromosomal aberrations and structural variants, which could give rise to an era of "next-generation cytogenetics."

Next-generation cytogenetics: Comprehensive assessment of 52 hematological malignancy genomes by optical genome mapping

Somatic structural variants (SVs) are important drivers of cancer development and progression. In a diagnostic set-up, especially for hematological malignancies, the comprehensive analysis of all SVs in a given sample still requires a combination of cytogenetic techniques, including karyotyping, FISH, and CNV microarrays. We hypothesize that the combination of these classical approaches could be replaced by optical genome mapping (OGM). Samples from 52 individuals with a clinical diagnosis of a hematological malignancy, divided into simple (<5 aberrations, n = 36) and complex (≥5 aberrations, n = 16) cases, were processed for OGM, reaching on average: 283-fold genome coverage. OGM called a total of 918 high-confidence SVs per sample, of which, on average, 13 were rare and >100 kb. In addition, on average, 73 CNVs were called per sample, of which six were >5 Mb. For the 36 simple cases, all clinically reported aberrations were detected, including deletions, insertions, inversions, aneuploidies, and translocations. For the 16 complex cases, results were largely concordant between standard-of-care and OGM, but OGM often revealed higher complexity than previously recognized. Detailed technical comparison with standard-of-care tests showed high analytical validity of OGM, resulting in a sensitivity of 100% and a positive predictive value of >80%. Importantly, OGM resulted in a more complete assessment than any previous single test and most likely reported the most accurate underlying genomic architecture (e.g., for complex translocations, chromoanagenesis, and marker chromosomes). In conclusion, the excellent concordance of OGM with diagnostic standard assays demonstrates its potential to replace classical cytogenetic tests as well as to rapidly map novel leukemia drivers.

Multicolor Whole-Genome Mapping in Nanochannels for Genetic Analysis

Analysis of structural variations (SVs) is important to understand mutations underlying genetic disorders and pathogenic conditions. However, characterizing SVs using short-read, high-throughput sequencing technology is difficult. Although long-read sequencing technologies are being increasingly employed in characterizing SVs, their low throughput and high costs discourage widespread adoption. Sequence motif-based optical mapping in nanochannels is useful in whole-genome mapping and SV detection, but it is not possible to precisely locate the breakpoints or estimate the copy numbers. We present here a universal multicolor mapping strategy in nanochannels combining conventional sequence-motif labeling system with Cas9-mediated target-specific labeling of any 20-base sequences (20mers) to create custom labels and detect new features. The sequence motifs are labeled with green fluorophores and the 20mers are labeled with red fluorophores. Using this strategy, it is possible to not only detect the SVs but also utilize custom labels to interrogate the features not accessible to motif-labeling, locate breakpoints, and precisely estimate copy numbers of genomic repeats. We validated our approach by quantifying the D4Z4 copy numbers, a known biomarker for facioscapulohumeral muscular dystrophy (FSHD) and estimating the telomere length, a clinical biomarker for assessing disease risk factors in aging-related diseases and malignant cancers. We also demonstrate the application of our methodology in discovering transposable long non-interspersed Elements 1 (LINE-1) insertions across the whole genome.

Optical genome mapping, a promising alternative to gold standard cytogenetic approaches in a series of acute lymphoblastic leukemias

Acute lymphoblastic leukemias (ALL) are characterized by a large number of cytogenetic abnormalities of clinical interest that require the use of several complementary techniques. Optical genome mapping (OGM) is based on analysis of ultra-high molecular weight DNA molecules that provides a high-resolution genome-wide analysis highlighting copy number and structural anomalies, including balanced translocations. We compared OGM to standard techniques (karyotyping, fluorescent in situ hybridization, single nucleotide polymorphism-array and reverse transcription multiplex ligation-dependent probe amplification) in 10 selected B or T-ALL. Eighty abnormalities were found using standard techniques of which 72 (90%) were correctly detected using OGM. Eight discrepancies were identified, while 12 additional anomalies were found by OGM. Among the discrepancies, four were detected in raw data but not retained because of filtering issues. However, four were truly missed, either because of a low variant allele frequency or because of a low coverage of some regions. Of the additional anomalies revealed by OGM, seven were confirmed by another technique, some of which are recurrent in ALL such as LMO2-TRA and MYC-TRB fusions. Despite false positive anomalies due to background noise and a case of inter-sample contamination secondarily identified, the OGM technology was relatively simple to use with little practice. Thus, OGM represents a promising alternative to cytogenetic techniques currently performed for ALL characterization. It enables a time and cost effective analysis allowing identification of complex cytogenetic events, including those currently inaccessible to standard techniques.

Fast and efficient Rmap assembly using the Bi-labelled de Bruijn graph

Genome wide optical maps are high resolution restriction maps that give a unique numeric representation to a genome. They are produced by assembling hundreds of thousands of single molecule optical maps, which are called Rmaps. Unfortunately, there are very few choices for assembling Rmap data. There exists only one publicly-available non-proprietary method for assembly and one proprietary software that is available via an executable. Furthermore, the publicly-available method, by Valouev et al. (Proc Natl Acad Sci USA 103(43):15770-15775, 2006), follows the overlap-layout-consensus (OLC) paradigm, and therefore, is unable to scale for relatively large genomes. The algorithm behind the proprietary method, Bionano Genomics' Solve, is largely unknown. In this paper, we extend the definition of bi-labels in the paired de Bruijn graph to the context of optical mapping data, and present the first de Bruijn graph based method for Rmap assembly. We implement our approach, which we refer to as RMAPPER, and compare its performance against the assembler of Valouev et al. (Proc Natl Acad Sci USA 103(43):15770-15775, 2006) and Solve by Bionano Genomics on data from three genomes: E. coli, human, and climbing perch fish (Anabas Testudineus). Our method was able to successfully run on all three genomes. The method of Valouev et al. (Proc Natl Acad Sci USA 103(43):15770-15775, 2006) only successfully ran on E. coli. Moreover, on the human genome RMAPPER was at least 130 times faster than Bionano Solve, used five times less memory and produced the highest genome fraction with zero mis-assemblies. Our software, RMAPPER is written in C++ and is publicly available under GNU General Public License at https://github.com/kingufl/Rmapper .

FaNDOM: Fast nested distance-based seeding of optical maps

Optical mapping (OM) provides single-molecule readouts of fluorescently labeled sequence motifs on long fragments of DNA, resolved to nucleotide-level coordinates. With the advent of microfluidic technologies for analysis of DNA molecules, it is possible to inexpensively generate long OM data ( > 150 kbp) at high coverage. In addition to scaffolding for de novo assembly, OM data can be aligned to a reference genome for identification of genomic structural variants. We introduce FaNDOM (Fast Nested Distance Seeding of Optical Maps)-an optical map alignment tool that greatly reduces the search space of the alignment process. On four benchmark human datasets, FaNDOM was significantly (4-14×) faster than competing tools while maintaining comparable sensitivity and specificity. We used FaNDOM to map variants in three cancer cell lines and identified many biologically interesting structural variants, including deletions, duplications, gene fusions and gene-disrupting rearrangements. FaNDOM is publicly available at https://github.com/jluebeck/FaNDOM.

Comparison of structural variants detected by optical mapping with long-read next-generation sequencing

MOTIVATION

Recent studies have shown the potential of using long-read whole-genome sequencing (WGS) approaches and optical mapping (OM) for the detection of clinically relevant structural variants (SVs) in cancer research. Three main long-read WGS platforms are currently in use: Pacific Biosciences (PacBio), Oxford Nanopore Technologies (ONT) and 10x Genomics. Recently, whole-genome OM technology (Bionano Genomics) has been introduced into human diagnostics. Questions remain about the accuracy of these long-read sequencing platforms, how comparable/interchangeable they are when searching for SVs and to what extent they can be replaced or supplemented by OM. Moreover, no tool can effectively compare SVs obtained by OM and WGS.

RESULTS

This study compared optical maps of the breast cancer cell line SKBR3 with AnnotSV outputs from WGS platforms. For this purpose, a software tool with comparative and filtering features was developed. The majority of SVs up to a 50 kbp distance variance threshold found by OM were confirmed by all WGS platforms, and 99% of translocations and 80% of deletions found by OM were confirmed by both PacBio and ONT, with ∼70% being confirmed by 10x Genomics in combination with PacBio and/or ONT. Interestingly, long deletions (>100 kbp) were detected only by 10x Genomics. Regarding insertions, ∼72% was confirmed by PacBio and ONT, but none by 10x Genomics. Inversions and duplications detected by OM were not detected by WGS. Moreover, the tool enabled the confirmation of SVs that overlapped in the same gene(s) and was applied to the filtering of disease-associated SVs.

AVAILABILITY

https://github.com/novosadt/om-annotsv-svc.

Research advances in and prospects of ornamental plant genomics

The term 'ornamental plant' refers to all plants with ornamental value, which generally have beautiful flowers or special plant architectures. China is rich in ornamental plant resources and known as the "mother of gardens". Genomics is the science of studying genomes and is useful for carrying out research on genome evolution, genomic variations, gene regulation, and important biological mechanisms based on detailed genome sequence information. Due to the diversity of ornamental plants and high sequencing costs, the progress of genome research on ornamental plants has been slow for a long time. With the emergence of new sequencing technologies and a reduction in costs since the whole-genome sequencing of the first ornamental plant (Prunus mume) was completed in 2012, whole-genome sequencing of more than 69 ornamental plants has been completed in <10 years. In this review, whole-genome sequencing and resequencing of ornamental plants will be discussed. We provide analysis with regard to basic data from whole-genome studies of important ornamental plants, the regulation of important ornamental traits, and application prospects.

Construction of Whole Genomes from Scaffolds Using Single Cell Strand-Seq Data

Accurate reference genome sequences provide the foundation for modern molecular biology and genomics as the interpretation of sequence data to study evolution, gene expression, and epigenetics depends heavily on the quality of the genome assembly used for its alignment. Correctly organising sequenced fragments such as contigs and scaffolds in relation to each other is a critical and often challenging step in the construction of robust genome references. We previously identified misoriented regions in the mouse and human reference assemblies using Strand-seq, a single cell sequencing technique that preserves DNA directionality Here we demonstrate the ability of Strand-seq to build and correct full-length chromosomes by identifying which scaffolds belong to the same chromosome and determining their correct order and orientation, without the need for overlapping sequences. We demonstrate that Strand-seq exquisitely maps assembly fragments into large related groups and chromosome-sized clusters without using new assembly data. Using template strand inheritance as a bi-allelic marker, we employ genetic mapping principles to cluster scaffolds that are derived from the same chromosome and order them within the chromosome based solely on directionality of DNA strand inheritance. We prove the utility of our approach by generating improved genome assemblies for several model organisms including the ferret, pig, Xenopus, zebrafish, Tasmanian devil and the Guinea pig.