Optical Genome Mapping as a Tool to Unveil New Molecular Findings in Hematological Patients with Complex Chromosomal Rearrangements.

Standard cytogenetic techniques (chromosomal banding analysis-CBA, and fluorescence in situ hybridization-FISH) show limits in characterizing complex chromosomal rearrangements and structural variants arising from two or more chromosomal breaks. In this study, we applied optical genome mapping (OGM) to fully characterize two cases of complex chromosomal rearrangements at high resolution. In case 1, an acute myeloid leukemia (AML) patient showing chromothripsis, OGM analysis was fully concordant with classic cytogenetic techniques and helped to better refine chromosomal breakpoints. The OGM results of case 2, a patient with non-Hodgkin lymphoma, were only partially in agreement with previous cytogenetic analyses and helped to better define clonal heterogeneity, overcoming the bias related to clonal selection due to cell culture of cytogenetic techniques. In both cases, OGM analysis led to the identification of molecular markers, helping to define the pathogenesis, classification, and prognosis of the analyzed patients. Despite extensive efforts to study hematologic diseases, standard cytogenetic methods display unsurmountable limits, while OGM is a tool that has the power to overcome these limitations and provide a cytogenetic analysis at higher resolution. As OGM also shows limits in defining regions of a repetitive nature, combining OGM with CBA to obtain a complete cytogenetic characterization would be desirable.

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.

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.

The Application of Optical Genome Mapping (OGM) in Severe Short Stature Caused by Duplication of 15q14q21.3.

(1) Background: Optical genome mapping (OGM) is a novel approach to identifying genomic structural variations with high accuracy and resolution. We report a proband with severe short stature caused by 46, XY, der (16) ins (16;15) (q23; q21.3q14) that was detected by OGM combined with other tests and review the clinical features of patients with duplication within 15q14q21.3; (2) Methods: OGM, whole exon sequencing (WES), copy number variation sequencing (CNV-seq), and karyotyping were used; (3) Results: The proband was a 10.7-year-old boy with a complaint of severe short stature (-3.41SDS) and abnormal gait. He had growth hormone deficiency, lumbar lordosis, and epiphyseal dysplasia of both femurs. WES and CNV-seq showed a 17.27 Mb duplication of chromosome 15, and there was an insertion in chromosome 16 found by karyotyping. Furthermore, OGM revealed that duplication of 15q14q21.3 was inversely inserted into 16q23.1, resulting in two fusion genes. A total of fourteen patients carried the duplication of 15q14q21.3, with thirteen previously reported and one from our center, 42.9% of which were de novo. In addition, neurologic symptoms (71.4%,10/14) were the most common phenotypes; (4) Conclusions: OGM combined with other genetic methods can reveal the genetic etiology of patients with the clinical syndrome, presenting great potential for use in properly diagnosing in the genetic cause of the clinical syndrome.

The AnnotSV webserver in 2023: updated visualization and ranking.

Much of the human genetics variant repertoire is composed of single nucleotide variants (SNV) and small insertion/deletions (indel) but structural variants (SV) remain a major part of our modified DNA. SV detection has often been a complex question to answer either because of the necessity to use different technologies (array CGH, SNP array, Karyotype, Optical Genome Mapping…) to detect each category of SV or to get an appropriate resolution (Whole Genome Sequencing). Thanks to the deluge of pangenomic analysis, Human geneticists are accumulating SV and their interpretation remains time consuming and challenging. The AnnotSV webserver (https://www.lbgi.fr/AnnotSV/) aims at being an efficient tool to (i) annotate and interpret SV potential pathogenicity in the context of human diseases, (ii) recognize potential false positive variants from all the SV identified and (iii) visualize the patient variants repertoire. The most recent developments in the AnnotSV webserver are: (i) updated annotations sources and ranking, (ii) three novel output formats to allow diverse utilization (analysis, pipelines), as well as (iii) two novel user interfaces including an interactive circos view.

DeepOM: single-molecule optical genome mapping via deep learning.

MOTIVATION: Efficient tapping into genomic information from a single microscopic image of an intact DNA molecule is an outstanding challenge and its solution will open new frontiers in molecular diagnostics. Here, a new computational method for optical genome mapping utilizing deep learning is presented, termed DeepOM. Utilization of a convolutional neural network, trained on simulated images of labeled DNA molecules, improves the success rate in the alignment of DNA images to genomic references. RESULTS: The method is evaluated on acquired images of human DNA molecules stretched in nano-channels. The accuracy of the method is benchmarked against state-of-the-art commercial software Bionano Solve. The results show a significant advantage in alignment success rate for molecules shorter than 50 kb. DeepOM improves the yield, sensitivity, and throughput of optical genome mapping experiments in applications of human genomics and microbiology. AVAILABILITY AND IMPLEMENTATION: The source code for the presented method is publicly available at https://github.com/yevgenin/DeepOM.

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.

Apolipoprotein L1 variability is associated with increased risk of renal failure in the Czech population.

BACKGROUND: With stage 5 chronic kidney disease (CKD5) more prevalent in the Czech Republic than in most European countries, genetic susceptibility is potentially implicated. METHODS: In a group of 1489 CKD5 kidney transplantation patients (93% with complete clinical characteristics; mean age 52.0 years, 37% females) and 2559 healthy controls (mean age 49.0 years, 51% females), we examined the prevalence of six APOL1 SNPs (rs73885319, rs71785313, rs13056427, rs136147, rs10854688 and rs9610473) and one newly detected 55-nucleotide insertion/deletion polymorphism. RESULTS: The rs73885319 and rs71785313 variants were monomorphic in the Czech Caucasian population. Genotype frequencies of the three SNPs examined (rs13056427, rs136147 and rs9610473) were almost identical in patients and controls (all P values were between 0.39 and 0.91). Minor homozygotes of rs10854688 were more common between the patients (13.2%) than in controls (10.7%) (OR [95% CI]; 1.32 [1.08-1.64]; P < 0.01). Prevalence of the newly detected 55-bp APOL1 deletion was significantly higher in CKD5 patients (3.0% vs. 1.7%; OR [95% CI]; 1.80 [1.16-2.80]; P < 0.01) compared to controls. Frequencies of some individual APOL1 haplotypes were borderline different between patients and controls. CONCLUSION: We found an association between rs10854688 SNP within the APOL1 gene and end-stage renal disease in the Czech Caucasian population. Further independent studies are required before a conclusive association between the newly detected APOL1 insertion/deletion polymorphism and CKD5 can be confirmed.

Application of Restriction Site-Associated DNA Sequencing (RAD-Seq) for Copy Number Variation and Triploidy Detection in Human.

At present, low-pass whole-genome sequencing (WGS) is frequently used in clinical research and in the screening of copy number variations (CNVs). However, there are still some challenges in the detection of triploids. Restriction site-associated DNA sequencing (RAD-Seq) technology is a reduced-representation genome sequencing technology developed based on next-generation sequencing. Here, we verified whether RAD-Seq could be employed to detect CNVs and triploids. In this study, genomic DNA of 11 samples was extracted employing a routine method and used to build libraries. Five cell lines of known karyotypes and 6 triploid abortion tissue samples were included for RAD-Seq testing. The triploid samples were confirmed by STR analysis and also tested by low-pass WGS. The accuracy and efficiency of detecting CNVs and triploids by RAD-Seq were then assessed, compared with low-pass WGS. In our results, RAD-Seq detected 11 out of 11 (100%) chromosomal abnormalities, including 4 deletions and 1 aneuploidy in the purchased cell lines and all triploid samples. By contrast, these triploids were missed by low-pass WGS. Furthermore, RAD-Seq showed a higher resolution and more accurate allele frequency in the detection of triploids than low-pass WGS. Our study shows that, compared with low-pass WGS, RAD-Seq has relatively higher accuracy in CNV detection at a similar cost and is capable of identifying triploids. Therefore, the application of this technique in medical genetics has a significant potential value.

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.

Validation of Optical Genome Mapping for the Molecular Diagnosis of Facioscapulohumeral Muscular Dystrophy.

The molecular diagnosis of facioscapulohumeral muscular dystrophy (FSHD) relies on detecting contractions of the unique D4Z4 repeat array at the chromosome 4q35 locus in the presence of a permissive 4q35A haplotype. Long, intact DNA molecules are required for accurate sizing of D4Z4 repeats. We validated the use of optical genome mapping to determine size and haplotype of D4Z4 alleles for FSHD analysis. The cohort included 36 unique DNA specimens from fresh blood samples or archived agarose plugs. High-molecular- weight DNA underwent sequence-specific labeling followed by separation and image analysis with data collection on the Saphyr system. D4Z4 allele sizes were calculated and haplotypes determined from the labeling patterns. Each specimen had previous diagnostic testing using restriction enzyme digests with EcoRI, EcoRI/BlnI, XapI, or HindIII, followed by pulsed field gel electrophoresis and Southern blot analysis with appropriate probes. Optical genome mapping detected 4q35 and 10q26 alleles ranging from 1 to 79 D4Z4 repeats and showed strong correlation with Southern blot allele sizing (R2 = 0.95) and haplotyping (133 of 134; 99.4% haplotype match). Analysis of inter-assay and intra-assay runs showed high reproducibility (0.03 to 0.94 %CV). Subsequent optical genome mapping for routine clinical testing from 315 clinical FSHD cases compared favorably with historical result trends. Optical genome mapping is an accurate and highly reproducible method for chromosomal abnormalities associated with FSHD.

A Restriction Endonuclease-Based Assay to Distinguish NANOGP8 Retrogene from Parental NANOG.

NANOG is an embryonic transcription factor, which gets reexpressed in cancer stem or tumor initiating cells. NANOGP8, a retrogene belonging to the NANOG family, is predominantly expressed in cancer cells and shows very high similarity with NANOG both at the nucleotide and at the protein level. The high similarity makes it extremely challenging to distinguish between these two transcription factors. Here we describe a highly efficient restriction endonuclease-based assay, which is performed on cDNA and allows to distinguish NANOGP8 from NANOG. This assay is critical to understand the specific role of NANOGP8 in cancer stemness, which in turn helps to unravel the therapeutic potential of targeting this undruggable transcription factor through gene therapy, for treatment of various cancers.

Generation of restriction endonucleases barcode map to trace SARS-CoV-2 origin and evolution.

Since the first report of SARS-CoV-2 in China in 2019, there has been a huge debate about the origin. In this work, using a different method we aimed to strengthen the observation that no evidence of genetic manipulation has been found by (1) detecting classical restriction site (RS) sequence in human SARS-CoV-2 genomes and (2) comparing them with other recombinant SARS-CoV-like virus created for experimental purposes. Finally, we propose a novel approach consisting in the generation of a restriction endonucleases site map of SARS-CoV-2 and other related coronavirus genomes to be used as a fingerprint to trace the virus evolution.

Development of a multiplex methylation-sensitive restriction enzyme-based SNP typing system for deconvolution of semen-containing mixtures.

The identification of mixed stains has always been a difficult problem in personal identification in the forensic field. In recent years, tissue-specific methylation sites have proven to be very stable biomarkers for distinguishing tissue origin. However, it is still challenging to perform tissue source identification and individual identification simultaneously. In this study, we developed a method that uses tissue-specific methylation markers combined with single-nucleotide polymorphism (SNP) markers to detect semen from mixed biofluids and to identify individuals simultaneously. Semen-specific CpG markers were chosen from the literature and further validated utilizing methylation-sensitive restriction endonuclease (MSRE) combined with PCR technology. The neighboring SNP markers were searched in the flanking sequence of the target CpG within 400 bp, and SNP typing was then carried out through a single-base extension reaction followed by capillary electrophoresis. Eventually, a method of MSRE combined with SNaPshot that could detect 12 compound CpG-SNP markers was developed. Using this system, 10 ng of total DNA and DNA mixture with semen content up to 25% could be typed successfully. Moreover, the cumulative discrimination power of the system in the northern Chinese Han population is 0.9998. This study provides a valuable strategy for forensic practice to perform tissue origin and individual identification from mixed stains simultaneously.

Bio-typing of Mycobacterium avium subspecies paratuberculosis isolates recovered from the Himalayan sheep and goats.

Information on bio-type profile of Mycobacterium avium subspecies paratubeculosis (MAP) in sheep flocks and goat herds of Himalayan region is not reported earlier. The aim of our study was to determine the bio-type of MAP infecting livestock of this region. A total of 71 faecal samples (sheep-57, goats-14) were screened by Ziehl-Neelsen (ZN) staining and IS900 PCR, and then processed for culture on Herrold's egg yolk medium (HEYM) having mycobactin J (MJ). Out of 71 faecal samples, MAP colonies were seen only in four samples (sheep-3 and goat-1). Isolates were confirmed as MAP on the basis of slow growth, acid fastness, MJ dependency, IS900 and IS1311 PCR. All the IS900 and IS1311 PCR positive samples were bio-typed by IS1311 PCR-REA (restriction endonuclease analysis), which confirmed all four isolates as 'bison type.' In IS1311 based phylogeny of MAP isolates by ClustalW method of the MegAlign program of DNASTAR Lasergene software, the four sequences of MAP isolates (NCBI sequence nos. MH988763, MH988765, MH988766 and MH988764) did not show any distinct clustering/grouping pattern. However, these four isolates showed a bit of closeness to the MAP sequences (KC990353.1 and KC990352.1) of 'bison type' isolated from wood bison in Canada. In conclusion, this is the first report on isolation and bio-type profile of MAP infecting sheep and goats of Himalayan region. Study will help in devising prevention and control strategies against spread of MAP infection in livestock population of Himalayan region.

A PCR plus restriction enzyme-based technique for detecting target-enzyme mutations at position Pro-106 in glyphosate-resistant Lolium perenne.

The first step in managing herbicide-resistant weeds is to confirm their resistance status. It is, therefore, crucial to have a rapid, reliable and cost-effective technique to assess samples for herbicide resistance. We designed and evaluated three derived cleaved amplified polymorphic sequence (dCAPS) markers for detecting glyphosate resistance in Lolium perenne. conferred by non-synonymous mutations at codon-106 in the enolpyruvylshikimate-3-phosphate synthase (EPSPS) gene. The dCAPS markers involve amplification of the target region, digestion of the amplified products with restriction enzymes and gel-based visualisation of the digested products. The results showed that all three dCAPS markers could successfully detect mutations at codon-106 in the target enzyme. The dCAPS markers can also inform us of the zygosity state of the resistance allele and was confirmed by sequencing the target region of the EPSPS gene. The markers described here are effective quick tests for the monitoring and evaluation of the target-enzyme mechanism of glyphosate resistance in Lolium perenne.

Relevance of ddRADseq method for species and population delimitation of closely related and widely distributed wolf spiders (Araneae, Lycosidae).

Although species delimitation is often controversial, emerging DNA-based and classical morphology-based methods are rarely compared using large-scale samplings, even less in the case of widely distributed species that have distant, allopatric populations. In the current study, we examined species boundaries within two wolf spider species of the genus Pardosa (Araneae, Lycosidae), P. riparia and P. palustris. Wolf spiders constitute an excellent model for testing the relevance of traditional vs. modern methods in species and population delimitation because several closely related species are distributed over cross-continental geographic ranges. Allopatric populations of the two Pardosa species were sampled across Europe to Far East Russia (latitudinal range > 150°) and several dozen individuals were studied using morphological characters (morphometry of three measures for both sexes, plus five in males only and two in females only), DNA barcoding (COI sequencing) and double-digest restriction site associated DNA sequencing (ddRADseq). The results obtained allow for changing the taxonomic status of two Far East Russian populations to subspecies and ddRADseq proved to be a powerful tool for taxonomic research despite scarce sampling and inherent subjectivity of species delimitation in allopatry. Overall, this study pleads for both multi-criteria and more population-based studies in taxonomy.

A survey of jaagsiekte sheep retrovirus (JSRV) infection in sheep in the three northeastern provinces of China.

Ovine pulmonary adenomatosis (OPA) is caused by jaagsiekte sheep retrovirus (JSRV) and is a chronic, progressive, and infectious neoplastic lung disease in sheep, which causes significant economic losses to the sheep industry. Neither a vaccine nor serological diagnostic methods to detect OPA are available. We performed a JSRV infection survey in sheep using blood samples (n = 1,372) collected in the three northeastern provinces of China (i.e., Inner Mongolia, Heilongjiang, and Jilin) to determine JSRV infection status in sheep herds using a real-time PCR assay targeting the gag gene of JSRV. The ovine endogenous retrovirus sequence was successfully amplified in all sheep samples tested (296 from the Inner Mongolia Autonomous Region, 255 from Jilin province, and 821 from Heilongjiang province). Subsequently, we attempted to distinguish exogenous JSRV (exJSRV) and endogenous JSRV (enJSRV) infections in these JSRV-positive samples using a combination assay that identifies a ScaI restriction site in an amplified 229-bp fragment of the gag gene of JSRV and a "LHMKYXXM" motif in the cytoplasmic tail region of the JSRV envelope protein. The ScaI restriction site is present in all known oncogenic JSRVs but absent in ovine endogenous retroviruses, while the "LHMKYXXM" motif is in all known exJSRVs but not in enJSRVs. Interestingly, one JSRV strain (HH13) from Heilongjiang province contained the "LHMKYXXM" motif but not the ScaI enzyme site. Phylogenetic analysis showed that strain HH13 was closely related to strain enJSRV-21 reported in the USA, indicating that HH13 could be an exogenous virus. Our results provide valuable information for further research on the genetic evolution and pathogenesis of JSRV.

Development of diagnostic assays for differentiation of atypical Aeromonas salmonicida vapA type V and type VI in ballan wrasse (Labrus bergylta, Ascanius).

Aeromonas salmonicida (As) is a highly heterogeneous bacterial species, and strains' host specificity has been reported. Ballan wrasse (Labrus bergylta Ascanius, 1767) is susceptible to atypical As (aAs) vapA type V and type VI in Scotland and Norway. Identification of the bacterium is achieved by culture and molecular techniques; however, the available methods used to distinguish the As types are costly and time-consuming. This paper describes the development of a PCR and a restriction enzyme assay for the detection of aAs vapA type V and type VI in ballan wrasse, respectively. Type V-specific primers were designed on conserved regions of the vapA gene, and the restriction enzyme assay was performed on the PCR products of the hypervariable region of vapA gene for the detection of type VI isolates. Amplification product was produced for type V (254 bp) and restriction bands (368 and 254 bp) for type VI isolates only. In addition, the assays detected type V and type VI isolates in spiked water samples and type V in diagnostic tissue samples. The assays are fast, specific and cost-effective and can be used as specific diagnostic tools for cleaner fish, to detect infectious divergence strains, and to manage and mitigate aAs disease outbreaks through vaccine development.

Poor performance of DNA barcoding and the impact of RAD loci filtering on the species delimitation of an Iberian ant-eating spider.

Genomic data provide unprecedented power for species delimitation. However, current implementations are still time and resource consuming. In addition, bioinformatic processing is contentious and its impact on downstream analyses is insufficiently understood. Here we employ ddRAD sequencing and a thorough sampling for species delimitation in Zodarion styliferum, a widespread Iberian ant-eating spider. We explore the influence of the loci filtering strategy on the downstream phylogenetic analyses, genomic clustering and coalescent species delimitation. We also assess the accuracy of one mitochondrial (COI) and one nuclear (ITS) barcode for fast and inexpensive species delineation in the group. Our genomic data strongly support two morphologically cryptic but ecologically divergent lineages, mainly restricted to the central-eastern and western parts of the Iberian Peninsula, respectively. Larger matrices with more missing data showed increased genomic diversity, supporting that bioinformatic strategies to maximize matrix completion disproportionately exclude loci with the highest mutation rates. Moderate loci filtering gave the best results across analyses: although larger matrices returned concatenated phylogenies with higher support, middle-sized matrices performed better in genetic structure analyses. COI displayed high diversity and a conspicuous barcode gap, revealing 13 mitochondrial lineages. Mitonuclear discordance is consistent with ancestral isolation in multiple groups, probably in glacial refugia, followed by range expansion and secondary contact that produced genomic homogenization. Several apparently (unidirectionally) introgressed specimens further challenge the accuracy of species identification through mitochondrial barcodes in the group. Conversely, ITS failed to separate both lineages of Z. styliferum. This study shows an extreme case of mitonuclear discordance that highlights the limitations of single molecular barcodes for species delimitation, even in presence of distinct barcode gaps, and brings new light on the effects of parameterization on shallow-divergence studies using RAD data.