Single-molecule Optical Mapping Research Articles

First-trimester noninvasive prenatal diagnosis of seven facioscapulohumeral muscular dystrophy type 1 families using SNP-based amplicon sequencing: An earlier, rapid and safer way.

The study is to explore the feasibility and value of SNP-based noninvasive prenatal diagnosis (NIPD) for facioscapulohumeral muscular dystrophy type 1 (FSHD1) in early pregnancy weeks. We prospectively collected seven FSHD1 families, with an average gestational age of 8+6. Among these seven couples, there were three affected FSHD1 mothers and four affected fathers. A multiplex-PCR panel comprising 402 amplicons was designed to selective enrich for highly heterozygous SNPs upstream of the DUX4 gene. Risk haplotype was constructed based on familial linkage analysis. Fetal genotypes were accurately inferred through relative haplotype dosage analysis using Bayes Factor. All tests were successfully completed in a single attempt, and no recombination events were detected. NIPD results were provided within a week, which is 4 weeks earlier than karyomapping and 7 weeks earlier than Bionano single-molecule optical mapping (BOM). Ultimately, five FSHD1 fetuses and two normal fetuses were successfully identified, with a 100% concordance rate with karyomapping and BOM. Therefore, SNP-based NIPD for FSHD1 was demonstrated to be feasible and accurate in early weeks of gestation, although the risk of recombination events cannot be completely eliminated. In the future, testing of more cases is still necessary to fully determine the clinical utility.

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.

PLPPR4 haploinsufficiency causes neurodevelopmental disorders by disrupting synaptic plasticity via mTOR signalling.

Phospholipid phosphatase related 4 (PLPPR4), a neuron-specific membrane protein located at the postsynaptic density of glutamatergic synapses, is a putative regulator of neuronal plasticity. However, PLPPR4 dysfunction has not been linked to genetic disorders. In this study, we report three unrelated patients with intellectual disability (ID) or autism spectrum disorder (ASD) who harbour a de novo heterozygous copy number loss of PLPPR4 in 1p21.2p21.3, a heterozygous nonsense mutation in PLPPR4 (NM_014839, c.4C > T, p.Gln2*) and a homozygous splice mutation in PLPPR4 (NM_014839: c.408 + 2 T > C), respectively. Bionano single-molecule optical mapping confirmed PLPPR4 deletion contains no additional pathogenic genes. Our results suggested that the loss of function of PLPPR4 is associated with neurodevelopmental disorders. To test the pathogenesis of PLPPR4, peripheral blood mononuclear cells obtained from the patient with heterozygous deletion of PLPPR4 were induced to specific iPSCs (CHWi001-A) and then differentiated into neurons. The neurons carrying the deletion of PLPPR4 displayed the reduced density of dendritic protrusions, shorter neurites and reduced axon length, suggesting the causal role of PLPPR4 in neurodevelopmental disorders. As the mTOR signalling pathway was essential for regulating the axon maturation and function, we found that mTOR signalling was inhibited with a higher level of p-AKT, p-mTOR and p-ERK1/2, decreased p-PI3K in PLPPR4-iPSCs neurons. Additionally, we found silencing PLPPR4 disturbed the mTOR signalling pathway. Our results suggested PLPPR4 modulates neurodevelopment by affecting the plasticity of neurons via the mTOR signalling pathway.

Complex 4q35 and 10q26 Rearrangements: A Challenge for Molecular Diagnosis of Patients With Facioscapulohumeral Dystrophy.

After clinical evaluation, the molecular diagnosis of type 1 facioscapulohumeral dystrophy (FSHD1) relies in most laboratories on the detection of a shortened D4Z4 array at the 4q35 locus by Southern blotting. In many instances, this molecular diagnosis remains inconclusive and requires additional experiments to determine the number of D4Z4 units or identify somatic mosaicism, 4q-10q translocations, and proximal p13E-11 deletions. These limitations highlight the need for alternative methodologies, illustrated by the recent emergence of novel technologies such as molecular combing (MC), single molecule optical mapping (SMOM), or Oxford Nanopore-based long-read sequencing providing a more comprehensive analysis of 4q and 10q loci. Over the last decade, MC revealed a further increasing complexity in the organization of the 4q and 10q distal regions in patients with FSHD with cis-duplication of D4Z4 arrays in approximately 1%-2% of cases. By using MC, we investigated in our center 2,363 cases for molecular diagnosis of FSHD. We also evaluated whether previously reported cis-duplications might be identified by SMOM using the Bionano EnFocus FSHD 1.0 algorithm. In our cohort of 2,363 samples, we identified 147 individuals carrying an atypical organization of the 4q35 or 10q26 loci. Mosaicism is the most frequent category followed by cis-duplications of the D4Z4 array. We report here chromosomal abnormalities of the 4q35 or 10q26 loci in 54 patients clinically described as FSHD, which are not present in the healthy population. In one-third of the 54 patients, these rearrangements are the only genetic defect suggesting that they might be causative of the disease. By analyzing DNA samples from 3 patients carrying a complex rearrangement of the 4q35 region, we further showed that the SMOM direct assembly of the 4q and 10q alleles failed to reveal these abnormalities and lead to negative results for FSHD molecular diagnosis. This work further highlights the complexity of the 4q and 10q subtelomeric regions and the need of in-depth analyses in a significant number of cases. This work also highlights the complexity of the 4q35 region and interpretation issues with consequences on the molecular diagnosis of patients or genetic counseling.

Dam Assisted Fluorescent Tagging of Chromatin Accessibility (DAFCA) for Optical Genome Mapping in Nanochannel Arrays.

Proteins and enzymes in the cell nucleus require physical access to their DNA target sites in order to perform genomic tasks such as gene activation and transcription. Hence, chromatin accessibility is a central regulator of gene expression, and its genomic profile holds essential information on the cell type and state. We utilized the E. coli Dam methyltransferase in combination with a fluorescent cofactor analogue to generate fluorescent tags in accessible DNA regions within the cell nucleus. The accessible portions of the genome are then detected by single-molecule optical genome mapping in nanochannel arrays. This method allowed us to characterize long-range structural variations and their associated chromatin structure. We show the ability to create whole-genome, allele-specific chromatin accessibility maps composed of long DNA molecules extended in silicon nanochannels.

Detection of cryptic balanced chromosomal rearrangements using high-resolution optical genome mapping

BackgroundChromosomal rearrangements have profound consequences in diverse human genetic diseases. Currently, the detection of balanced chromosomal rearrangements (BCRs) mainly relies on routine cytogenetic G-banded karyotyping. However, cryptic BCRs are hard...

High-resolution breakpoint junction mapping of proximally extended D4Z4 deletions in FSHD1 reveals evidence for a founder effect.

Facioscapulohumeral muscular dystrophy (FSHD) is an inherited myopathy clinically characterized by weakness in the facial, shoulder girdle and upper a muscles. FSHD is caused by chromatin relaxation of the D4Z4 macrosatellite repeat, mostly by a repeat contraction, facilitating ectopic expression of DUX4 in skeletal muscle. Genetic diagnosis for FSHD is generally based on the sizing and haplotyping of the D4Z4 repeat on chromosome 4 by Southern blotting (SB), molecular combing or single-molecule optical mapping, which is usually straight forward but can be complicated by atypical rearrangements of the D4Z4 repeat. One of these rearrangements is a D4Z4 proximally extended deletion (DPED) allele, where not only the D4Z4 repeat is partially deleted, but also sequences immediately proximal to the repeat are lost, which can impede accurate diagnosis in all genetic methods. Previously, we identified several DPED alleles in FSHD and estimated the size of the proximal deletions by a complex pulsed-field gel electrophoresis and SB strategy. Here, using the next-generation sequencing, we have defined the breakpoint junctions of these DPED alleles at the base pair resolution in 12 FSHD families and 4 control individuals facilitating a PCR-based diagnosis of these DPED alleles. Our resultsshow that half of the DPED alleles are derivates of an ancient founder allele. For some DPED alleles, we found that genetic elements are deleted such as DUX4c, FRG2, DBE-T and myogenic enhancers necessitating re-evaluation of their role in FSHD pathogenesis.

Stochastic initiation of DNA replication across the human genome

Wang etal. (2021) comprehensively map DNA replication initiation events across the human genome using single-molecule optical resolution mapping and find that initiation events are randomly distributed across broad initiation zones that are only utilized in a stochastic fashion across a population of cells.

Single-molecule optical genome mapping in nanochannels: multidisciplinarity at the nanoscale.

The human genome contains multiple layers of information that extend beyond the genetic sequence. In fact, identical genetics do not necessarily yield identical phenotypes as evident for the case of two different cell types in the human body. The great variation in structure and function displayed by cells with identical genetic background is attributed to additional genomic information content. This includes large-scale genetic aberrations, as well as diverse epigenetic patterns that are crucial for regulating specific cell functions. These genetic and epigenetic patterns operate in concert in order to maintain specific cellular functions in health and disease. Single-molecule optical genome mapping is a high-throughput genome analysis method that is based on imaging long chromosomal fragments stretched in nanochannel arrays. The access to long DNA molecules coupled with fluorescent tagging of various genomic information presents a unique opportunity to study genetic and epigenetic patterns in the genome at a single-molecule level over large genomic distances. Optical mapping entwines synergistically chemical, physical, and computational advancements, to uncover invaluable biological insights, inaccessible by sequencing technologies. Here we describe the method’s basic principles of operation, and review the various available mechanisms to fluorescently tag genomic information. We present some of the recent biological and clinical impact enabled by optical mapping and present recent approaches for increasing the method’s resolution and accuracy. Finally, we discuss how multiple layers of genomic information may be mapped simultaneously on the same DNA molecule, thus paving the way for characterizing multiple genomic observables on individual DNA molecules.

Single-molecule optical mapping of the distribution of DNA phosphorothioate epigenetics

DNA phosphorothioate (PT) modifications, with the nonbridging phosphate oxygen replaced by sulfur, governed by DndABCDE or SspABCD, are widely distributed in prokaryotes and have a highly unusual feature of occupying only a small portion of available consensus sequences in a genome. Despite the presence of plentiful non-PT-protected consensuses, DNA PT modification is still employed as a recognition tag by the restriction cognate, for example, DndFGH or SspE, to discriminate and destroy PT-lacking foreign DNA. This raises a fundamental question about how PT modifications are distributed along DNA molecules to keep the restriction components in check. Here, we present two single-molecule strategies that take advantage of the nucleophilicity of PT in combination with fluorescent markers for optical mapping of both single- and double-stranded PT modifications across individual DNA molecules. Surprisingly, PT profiles vary markedly from molecule to molecule, with different PT locations and spacing distances between PT pairs, even in the presence of DndFGH or SspE. The results revealed unprecedented PT modification features previously obscured by ensemble averaging, providing novel insights into the riddles regarding unusual target selection by PT modification and restriction components.

Analysis of balanced reciprocal translocations in patients with subfertility using single-molecule optical mapping

Approximately 1% of individuals who carry a balanced reciprocal translocation (BRT) are subfertile. Current karyotyping does not have the resolution to determine whether the breakpoints of the involved chromosomes perturb genes important for fertility. The aim of this study was to apply single-molecule optical mapping (SMOM) to patients presenting for IVF (in vitro fertilization) to ascertain whether the BRT disrupted any genes associated with normal fertility. Nine subfertile patients with different BRTs were recruited for the study. Methyltransferase enzyme DLE1 was used to fluorescently label their genomic DNA samples at the recognition motif CTTAAG. The SMOM was performed on the Bionano platform, and long molecules aligned against the reference genome hg19 to identify the breakpoint regions. Mate-pair and PCR-Sanger sequencing were used to confirm the precise breakpoint sequences. Both breakpoint regions in each of the nine BRTs were finely mapped to small regions of approximately 10Kb, and their positions were consistent with original cytogenetic banding patterns determined by karyotyping. In three BRTs, breakpoints disrupted genes known to be associated with male infertility, namely NUP155 and FNDC3A [46,XY,t(5;13)(p15;q22)], DPY19L1 [46,XY,t(1;7)(p36.3;p15), and BAI3 [46,XY,t(3;6)(p21;q16)]. The SMOM has potential clinical application as a rapid tool to screen patients with BRTs for underlying genetic causes of infertility and other diseases.

The Driver of Extreme Human-Specific Olduvai Repeat Expansion Remains Highly Active in the Human Genome.

Sequences encoding Olduvai protein domains (formerly DUF1220) show the greatest human lineage-specific increase in copy number of any coding region in the genome and have been associated, in a dosage-dependent manner, with brain size, cognitive aptitude, autism, and schizophrenia. Tandem intragenic duplications of a three-domain block, termed the Olduvai triplet, in four NBPF genes in the chromosomal 1q21.1-0.2 region, are primarily responsible for the striking human-specific copy number increase. Interestingly, most of the Olduvai triplets are adjacent to, and transcriptionally coregulated with, three human-specific NOTCH2NL genes that have been shown to promote cortical neurogenesis. Until now, the underlying genomic events that drove the Olduvai hyperamplification in humans have remained unexplained. Here, we show that the presence or absence of an alternative first exon of the Olduvai triplet perfectly discriminates between amplified (58/58) and unamplified (0/12) triplets. We provide sequence and breakpoint analyses that suggest the alternative exon was produced by an nonallelic homologous recombination-based mechanism involving the duplicative transposition of an existing Olduvai exon found in the CON3 domain, which typically occurs at the C-terminal end of NBPF genes. We also provide suggestive in vitro evidence that the alternative exon may promote instability through a putative G-quadraplex (pG4)-based mechanism. Lastly, we use single-molecule optical mapping to characterize the intragenic structural variation observed in NBPF genes in 154 unrelated individuals and 52 related individuals from 16 families and show that the presence of pG4-containing Olduvai triplets is strongly correlated with high levels of Olduvai copy number variation. These results suggest that the same driver of genomic instability that allowed the evolutionarily recent, rapid, and extreme human-specific Olduvai expansion remains highly active in the human genome.

Single-molecule optical mapping enables quantitative measurement of D4Z4 repeats in facioscapulohumeral muscular dystrophy (FSHD)

PurposeFacioscapulohumeral muscular dystrophy (FSHD) is a common adult muscular dystrophy. Over 95% of FSHD cases are associated with contraction of the D4Z4 tandem repeat (~3.3 kb per unit) at 4q35...

Genome of the Komodo dragon reveals adaptations in the cardiovascular and chemosensory systems of monitor lizards

Monitor lizards are unique among ectothermic reptiles in that they have high aerobic capacity and distinctive cardiovascular physiology resembling that of endothermic mammals. Here, we sequence the genome of the Komodo dragon Varanus komodoensis, the largest extant monitor lizard, and generate a high-resolution de novo chromosome-assigned genome assembly for V. komodoensis using a hybrid approach of long-range sequencing and single-molecule optical mapping. Comparing the genome of V. komodoensis with those of related species, we find evidence of positive selection in pathways related to energy metabolism, cardiovascular homoeostasis, and haemostasis. We also show species-specific expansions of a chemoreceptor gene family related to pheromone and kairomone sensing in V. komodoensis and other lizard lineages. Together, these evolutionary signatures of adaptation reveal the genetic underpinnings of the unique Komodo dragon sensory and cardiovascular systems, and suggest that selective pressure altered haemostasis genes to help Komodo dragons evade the anticoagulant effects of their own saliva. The Komodo dragon genome is an important resource for understanding the biology of monitor lizards and reptiles worldwide.

Enrichment of megabase-sized DNA molecules for single-molecule optical mapping and next-generation sequencing

Next-generation sequencing (NGS) has caused a revolution, yet left a gap: long-range genetic information from native, non-amplified DNA fragments is unavailable. It might be obtained by optical mapping of megabase-sized DNA molecules. Frequently only a specific genomic region is of interest, so here we introduce a method for selection and enrichment of megabase-sized DNA molecules intended for single-molecule optical mapping: DNA from a human cell line is digested by the NotI rare-cutting enzyme and size-selected by pulsed-field gel electrophoresis. For demonstration, more than 600 sub-megabase- to megabase-sized DNA molecules were recovered from the gel and analysed by denaturation-renaturation optical mapping. Size-selected molecules from the same gel were sequenced by NGS. The optically mapped molecules and the NGS reads showed enrichment from regions defined by NotI restriction sites. We demonstrate that the unannotated genome can be characterized in a locus-specific manner via molecules partially overlapping with the annotated genome. The method is a promising tool for investigation of structural variants in enriched human genomic regions for both research and diagnostic purposes. Our enrichment method could potentially work with other genomes or target specified regions by applying other genomic editing tools, such as the CRISPR/Cas9 system.

Single-molecule optical genome mapping of a human HapMap and a colorectal cancer cell line.

BackgroundNext-generation sequencing (NGS) technologies have changed our understanding of the variability of the human genome. However, the identification of genome structural variations based on NGS approaches with read lengths of 35–300 bases remains a challenge. Single-molecule optical mapping technologies allow the analysis of DNA molecules of up to 2 Mb and as such are suitable for the identification of large-scale genome structural variations, and for de novo genome assemblies when combined with short-read NGS data. Here we present optical mapping data for two human genomes: the HapMap cell line GM12878 and the colorectal cancer cell line HCT116.FindingsHigh molecular weight DNA was obtained by embedding GM12878 and HCT116 cells, respectively, in agarose plugs, followed by DNA extraction under mild conditions. Genomic DNA was digested with KpnI and 310,000 and 296,000 DNA molecules (≥150 kb and 10 restriction fragments), respectively, were analyzed per cell line using the Argus optical mapping system. Maps were aligned to the human reference by OPTIMA, a new glocal alignment method. Genome coverage of 6.8× and 5.7× was obtained, respectively; 2.9× and 1.7× more than the coverage obtained with previously available software.ConclusionsOptical mapping allows the resolution of large-scale structural variations of the genome, and the scaffold extension of NGS-based de novo assemblies. OPTIMA is an efficient new alignment method; our optical mapping data provide a resource for genome structure analyses of the human HapMap reference cell line GM12878, and the colorectal cancer cell line HCT116.

DNA fluorocode: A single molecule, optical map of DNA with nanometre resolution

We present a new method for single-molecule optical DNA mapping using an exceptionally dense, yet sequence-specific coverage of DNA with a fluorescent probe. The method employs a DNA methyltransferase enzyme to direct the DNA labelling, followed by molecular combing of the DNA onto a polymer-coated surface and subsequent sub-diffraction limit localization of the fluorophores. The result is a ‘DNA fluorocode’; a simple description of the DNA sequence, with a maximum achievable resolution of less than 20 bases, which can be read and analyzed like a barcode. We demonstrate the generation of a fluorocode for genomic DNA from the lambda bacteriophage using a DNA methyltransferase, M.HhaI, to direct fluorescent labels to four-base sequences reading 5′-GCGC-3′. A consensus fluorocode that allows the study of the DNA sequence at the level of an individual labelling site can be generated from a handful of molecules.