Mechanisms of action and resistance in histone methylation-targeted therapy.

Epigenomes enable the rectification of disordered cancer gene expression, thereby providing new targets for pharmacological interventions. The clinical utility of targeting histone H3 lysine trimethylation (H3K27me3) as an epigenetic hallmark has been demonstrated1-7. However, in actual therapeutic settings, the mechanism by which H3K27me3-targeting therapies exert their effects and the response of tumour cells remain unclear. Here we show the potency and mechanisms of action and resistance of the EZH1-EZH2 dual inhibitor valemetostat in clinical trials of patients with adult T cell leukaemia/lymphoma. Administration of valemetostat reduced tumour size and demonstrated durable clinical response in aggressive lymphomas with multiple genetic mutations. Integrative single-cell analyses showed that valemetostat abolishes the highly condensed chromatin structure formed by the plastic H3K27me3 and neutralizes multiple gene loci, including tumour suppressor genes. Nevertheless, subsequent long-term treatment encounters the emergence of resistant clones with reconstructed aggregate chromatin that closely resemble the pre-dose state. Acquired mutations at the PRC2-compound interface result in the propagation of clones with increased H3K27me3 expression. In patients free of PRC2 mutations, TET2 mutation or elevated DNMT3A expression causes similar chromatin recondensation through de novo DNA methylation in the H3K27me3-associated regions. We identified subpopulations with distinct metabolic and gene translation characteristics implicated in primary susceptibility until the acquisition of the heritable (epi)mutations. Targeting epigenetic drivers and chromatin homeostasis may provide opportunities for further sustained epigenetic cancer therapies.

The non-canonical bivalent gene Wfdc15a controls spermatogenic protease and immune homeostasis.

Male infertility can be caused by chromosomal abnormalities, mutations and epigenetic defects. Epigenetic modifiers pre-program hundreds of spermatogenic genes in spermatogonial stem cells (SSCs) for expression later in spermatids, but it remains mostly unclear whether and how those genes are involved in fertility. Here, we report that Wfdc15a, a WFDC family protease inhibitor pre-programmed by KMT2B, is essential for spermatogenesis. We found that Wfdc15a is a non-canonical bivalent gene carrying both H3K4me3 and facultative H3K9me3 in SSCs, but is later activated along with the loss of H3K9me3 and acquisition of H3K27ac during meiosis. We show that WFDC15A deficiency causes defective spermiogenesis at the beginning of spermatid elongation. Notably, depletion of WFDC15A causes substantial disturbance of the testicular protease-antiprotease network and leads to an orchitis-like inflammatory response associated with TNFα expression in round spermatids. Together, our results reveal a unique epigenetic program regulating innate immunity crucial for fertility.

An improved method for the highly specific detection of transcription start sites.

Precise detection of the transcriptional start site (TSS) is a key for characterizing transcriptional regulation of genes and for annotation of newly sequenced genomes. Here, we describe the development of an improved method, designated 'TSS-seq2.' This method is an iterative improvement of TSS-seq, a previously published enzymatic cap-structure conversion method to detect TSSs in base sequences. By modifying the original procedure, including by introducing split ligation at the key cap-selection step, the yield and the accuracy of the reaction has been substantially improved. For example, TSS-seq2 can be conducted using as little as 5 ng of total RNA with an overall accuracy of 96%; this yield a less-biased and more precise detection of TSS. We then applied TSS-seq2 for TSS analysis of four plant species that had not yet been analyzed by any previous TSS method.

Direct Reprogramming Improves Cardiac Function and Reverses Fibrosis in Chronic Myocardial Infarction.

BACKGROUND: Because adult cardiomyocytes have little regenerative capacity, resident cardiac fibroblasts (CFs) synthesize extracellular matrix after myocardial infarction (MI) to form fibrosis, leading to cardiac dysfunction and heart failure. Therapies that can regenerate the myocardium and reverse fibrosis in chronic MI are lacking. The overexpression of cardiac transcription factors, including Mef2c/Gata4/Tbx5/Hand2 (MGTH), can directly reprogram CFs into induced cardiomyocytes (iCMs) and improve cardiac function under acute MI. However, the ability of in vivo cardiac reprogramming to repair chronic MI with established scars is undetermined. METHODS: We generated a novel Tcf21iCre/reporter/MGTH2A transgenic mouse system in which tamoxifen treatment could induce both MGTH and reporter expression in the resident CFs for cardiac reprogramming and fibroblast lineage tracing. We first tested the efficacy of this transgenic system in vitro and in vivo for acute MI. Next, we analyzed in vivo cardiac reprogramming and fusion events under chronic MI using Tcf21iCre/Tomato/MGTH2A and Tcf21iCre/mTmG/MGTH2A mice, respectively. Microarray and single-cell RNA sequencing were performed to determine the mechanism of cardiac repair by in vivo reprogramming. RESULTS: We confirmed the efficacy of transgenic in vitro and in vivo cardiac reprogramming for acute MI. In chronic MI, in vivo cardiac reprogramming converted ≈2% of resident CFs into iCMs, in which a majority of iCMs were generated by means of bona fide cardiac reprogramming rather than by fusion with cardiomyocytes. Cardiac reprogramming significantly improved myocardial contraction and reduced fibrosis in chronic MI. Microarray analyses revealed that the overexpression of MGTH activated cardiac program and concomitantly suppressed fibroblast and inflammatory signatures in chronic MI. Single-cell RNA sequencing demonstrated that resident CFs consisted of 7 subclusters, in which the profibrotic CF population increased under chronic MI. Cardiac reprogramming suppressed fibroblastic gene expression in chronic MI by means of conversion of profibrotic CFs to a quiescent antifibrotic state. MGTH overexpression induced antifibrotic effects partly by suppression of Meox1, a central regulator of fibroblast activation. CONCLUSIONS: These results demonstrate that cardiac reprogramming could repair chronic MI by means of myocardial regeneration and reduction of fibrosis. These findings present opportunities for the development of new therapies for chronic MI and heart failure.

Positive selective sweeps of epigenetic mutations regulating specialized metabolites in plants.

DNA methylation is an important factor regulating gene expression in organisms. However, whether DNA methylation plays a key role in adaptive evolution is unknown. Here, we show evidence of naturally selected DNA methylation in Arabidopsis thaliana In comparison with single nucleotide polymorphisms, three types of methylation-methylated CGs (mCGs), mCHGs, and mCHHs-contributed highly to variable gene expression levels among an A thaliana population. Such variably expressed genes largely affect a large variation of specialized metabolic quantities. Among the three types of methylations, only mCGs located in promoter regions of genes associated with specialized metabolites show a selective sweep signature in the A thaliana population. Thus, naturally selected mCGs appear to be key mutations that cause the expressional diversity associated with specialized metabolites during plant evolution.

Tsga8 is required for spermatid morphogenesis and male fertility in mice.

During spermatogenesis, intricate gene expression is coordinately regulated by epigenetic modifiers, which are required for differentiation of spermatogonial stem cells (SSCs) contained among undifferentiated spermatogonia. We have previously found that KMT2B conveys H3K4me3 at bivalent and monovalent promoters in undifferentiated spermatogonia. Because these genes are expressed late in spermatogenesis or during embryogenesis, we expect that many of them are potentially programmed by KMT2B for future expression. Here, we show that one of the genes targeted by KMT2B, Tsga8, plays an essential role in spermatid morphogenesis. Loss of Tsga8 in mice leads to male infertility associated with abnormal chromosomal distribution in round spermatids, malformation of elongating spermatid heads and spermiation failure. Tsga8 depletion leads to dysregulation of thousands of genes, including the X-chromosome genes that are reactivated in spermatids, and insufficient nuclear condensation accompanied by reductions of TNP1 and PRM1, key factors for histone-to-protamine transition. Intracytoplasmic sperm injection (ICSI) of spermatids rescued the infertility phenotype, suggesting competency of the spermatid genome for fertilization. Thus, Tsga8 is a KMT2B target that is vitally necessary for spermiogenesis and fertility.

Downregulation of LAMB3 Altered the Carcinogenic Properties of Human Papillomavirus 16-Positive Cervical Cancer Cells.

Nearly all cervical cancer cases are caused by infection with high-risk human papillomavirus (HR-HPV) types. The mechanism of cervical cell transformation is related to the powerful action of viral oncoproteins and cellular gene alterations. Transcriptomic data from cervical cancer and normal cervical cells were utilized to identify upregulated genes and their associated pathways. The laminin subunit beta-3 (LAMB3) mRNAwas overexpressed in cervical cancer and was chosen for functional analysis. The LAMB3 was predominantly expressed in the extracellular region and the plasma membrane, which play a role in protein binding and cell adhesion molecule binding, leading to cell migration and tissue development. LAMB3 was found to be implicated in the pathway in cancer and the PI3K-AKT signaling pathway. LAMB3 knockdown decreased cell migration, invasion, anchorage-dependent and anchorage-independent cell growth and increased the number of apoptotic cells. These effects were linked to a decrease in protein levels involved in the PI3K-AKT signaling pathway and an increase in p53 protein. This study demonstrated that LAMB3 could promote cervical cancer cell migration, invasion and survival.

Metabolic labeling of RNA using multiple ribonucleoside analogs enables the simultaneous evaluation of RNA synthesis and degradation rates.

Gene expression is determined by a balance between RNA synthesis and RNA degradation. To elucidate the underlying regulatory mechanisms and principles of this, simultaneous measurements of RNA synthesis and degradation are required. Here, we report the development of "Dyrec-seq," which uses 4-thiouridine and 5-bromouridine to simultaneously quantify RNA synthesis and degradation rates. Dyrec-seq enabled the quantification of RNA synthesis and degradation rates of 4702 genes in HeLa cells. Functional enrichment analysis showed that the RNA synthesis and degradation rates of genes are actually determined by the genes' biological functions. A comparison of theoretical and experimental analyses revealed that the amount of RNA is determined by the ratio of RNA synthesis to degradation rates, whereas the rapidity of responses to external stimuli is determined only by the degradation rate. This study emphasizes that not only RNA synthesis but also RNA degradation is important in shaping gene expression patterns.

Long-read sequencing for non-small-cell lung cancer genomes.

Here, we report the application of a long-read sequencer, PromethION, for analyzing human cancer genomes. We first conducted whole-genome sequencing on lung cancer cell lines. We found that it is possible to genotype known cancerous mutations, such as point mutations. We also found that long-read sequencing is particularly useful for precisely identifying and characterizing structural aberrations, such as large deletions, gene fusions, and other chromosomal rearrangements. In addition, we identified several medium-sized structural aberrations consisting of complex combinations of local duplications, inversions, and microdeletions. These complex mutations occurred even in key cancer-related genes, such as STK11, NF1, SMARCA4, and PTEN The biological relevance of those mutations was further revealed by epigenome, transcriptome, and protein analyses of the affected signaling pathways. Such structural aberrations were also found in clinical lung adenocarcinoma specimens. Those structural aberrations were unlikely to be reliably detected by conventional short-read sequencing. Therefore, long-read sequencing may contribute to understanding the molecular etiology of patients for whom causative cancerous mutations remain unknown and therapeutic strategies are elusive.

UTX deficiency in neural stem/progenitor cells results in impaired neural development, fetal ventriculomegaly, and postnatal death.

Recent studies have demonstrated that epigenetic modifications are deeply involved in neurogenesis; however, the precise mechanisms remain largely unknown. To determine the role of UTX (also known as KDM6A), a demethylase of histone H3K27, in neural development, we generated Utx-deficient mice in neural stem/progenitor cells (NSPCs). Since Utx is an X chromosome-specific gene, the genotypes are sex-dependent; female mice lose both Utx alleles (UtxΔ/Δ ), and male mice lose one Utx allele yet retain one Uty allele, the counterpart of Utx on the Y chromosome (UtxΔ/Uty ). We found that UtxΔ/Δ mice exhibited fetal ventriculomegaly and died soon after birth. Immunofluorescence staining and EdU labeling revealed a significant increase in NSPCs and a significant decrease in intermediate-progenitor and differentiated neural cells. Molecular analyses revealed the downregulation of pathways related to DNA replication and increased H3K27me3 levels around the transcription start sites in UtxΔ/Δ NSPCs. These results indicate that UTX globally regulates the expression of genes required for proper neural development in NSPCs, and UTX deficiency leads to impaired cell cycle exit, reduced differentiation, and neonatal death. Interestingly, although UtxΔ/Uty mice survived the postnatal period, most died of hydrocephalus, a clinical feature of Kabuki syndrome, a congenital anomaly involving UTX mutations. Our findings provide novel insights into the role of histone modifiers in neural development and suggest that UtxΔ/Uty mice are a potential disease model for Kabuki syndrome.

Long-read whole-genome methylation patterning using enzymatic base conversion and nanopore sequencing.

Long-read whole-genome sequencing analysis of DNA methylation would provide useful information on the chromosomal context of gene expression regulation. Here we describe the development of a method that improves the read length generated by using the bisulfite-sequencing-based approach. In this method, we combined recently developed enzymatic base conversion, where an unmethylated cytosine (C) should be converted to thymine (T), with nanopore sequencing. After methylation-sensitive base conversion, the sequencing library was constructed using long-range polymerase chain reaction. This type of analysis is possible using a minimum of 1 ng genomic DNA, and an N50 read length of 3.4-7.6 kb is achieved. To analyze the produced data, which contained a substantial number of base mismatches due to sequence conversion and an inaccurate base read of the nanopore sequencing, a new analytical pipeline was constructed. To demonstrate the performance of long-read methylation sequencing, breast cancer cell lines and clinical specimens were subjected to analysis, which revealed the chromosomal methylation context of key cancer-related genes, allele-specific methylated genes, and repetitive or deletion regions. This method should convert the intractable specimens for which the amount of available genomic DNA is limited to the tractable targets.

Directed induction of alveolar type I cells derived from pluripotent stem cells via Wnt signaling inhibition.

Alveologenesis is a developmental step involving the expansion of the lung surface area which is essential for gas exchange. The gas exchange process is mediated by alveolar type I (AT1) cells, which are known to be differentiated from alveolar type II (AT2) or bipotent cells. Due to the difficulty of isolating and culturing primary AT1 cells, the mechanism underlying their differentiation is not completely understood. We performed single-cell RNA sequencing (scRNA-seq) of fibroblast-dependent alveolar organoids (FD-AOs), including human induced pluripotent stem cell (hiPSC)-derived epithelial cells and fetal lung fibroblasts, and identified hiPSC-derived AT1 (iAT1) cells. A comparison of the FD-AOs and fibroblast-free alveolar organoids showed that iAT1 cells were mainly present in the FD-AOs. Importantly, the transcriptomes of iAT1 cells were remarkably similar to those of primary AT1 cells. Additionally, XAV-939, a tankyrase inhibitor, increased iAT1 cells in passaged FD-AOs, suggesting that these cells were differentiated from hiPSC-derived AT2 (iAT2) cells through the inhibition of canonical Wnt signaling. Consequently, our scRNA-seq data allowed us to define iAT1 cells and identify FD-AOs as a useful model for investigating the mechanism underlying human AT1 cell differentiation from AT2 cells in vitro.

MCM10 compensates for Myc-induced DNA replication stress in breast cancer stem-like cells.

Cancer stem-like cells (CSCs) induce drug resistance and recurrence of tumors when they experience DNA replication stress. However, the mechanisms underlying DNA replication stress in CSCs and its compensation remain unclear. Here, we demonstrate that upregulated c-Myc expression induces stronger DNA replication stress in patient-derived breast CSCs than in differentiated cancer cells. Our results suggest critical roles for mini-chromosome maintenance protein 10 (MCM10), a firing (activating) factor of DNA replication origins, to compensate for DNA replication stress in CSCs. MCM10 expression is upregulated in CSCs and is maintained by c-Myc. c-Myc-dependent collisions between RNA transcription and DNA replication machinery may occur in nuclei, thereby causing DNA replication stress. MCM10 may activate dormant replication origins close to these collisions to ensure the progression of replication. Moreover, patient-derived breast CSCs were found to be dependent on MCM10 for their maintenance, even after enrichment for CSCs that were resistant to paclitaxel, the standard chemotherapeutic agent. Further, MCM10 depletion decreased the growth of cancer cells, but not of normal cells. Therefore, MCM10 may robustly compensate for DNA replication stress and facilitate genome duplication in cancer cells in the S-phase, which is more pronounced in CSCs. Overall, we provide a preclinical rationale to target the c-Myc-MCM10 axis for preventing drug resistance and recurrence of tumors.

Kmt2b conveys monovalent and bivalent H3K4me3 in mouse spermatogonial stem cells at germline and embryonic promoters.

The mammalian male germline is sustained by a pool of spermatogonial stem cells (SSCs) that can transmit both genetic and epigenetic information to offspring. However, the mechanisms underlying epigenetic transmission remain unclear. The histone methyltransferase Kmt2b is highly expressed in SSCs and is required for the SSC-to-progenitor transition. At the stem-cell stage, Kmt2b catalyzes H3K4me3 at bivalent H3K27me3-marked promoters as well as at promoters of a new class of genes lacking H3K27me3, which we call monovalent. Monovalent genes are mainly activated in late spermatogenesis, whereas most bivalent genes are mainly not expressed until embryonic development. These data suggest that SSCs are epigenetically primed by Kmt2b in two distinguishable ways for the upregulation of gene expression both during the spermatogenic program and through the male germline into the embryo. Because Kmt2b is also the major H3K4 methyltransferase for bivalent promoters in embryonic stem cells, we also propose that Kmt2b has the capacity to prime stem cells epigenetically.

Long-term expansion of alveolar stem cells derived from human iPS cells in organoids.

The stable expansion of tissue-specific stem cells in vitro has contributed to research on several organs. Alveolar epithelial type II (AT2) cells function as tissue stem cells in the lung, but robust models for studying human AT2 cells are lacking. Here we report a method for the efficient generation and long-term expansion of alveolar organoids (AOs) harboring SFTPC+ alveolar stem cells derived from human induced pluripotent stem cells (hiPSCs). hiPSC-derived SFTPC+ cells self-renewed, with transcriptomes and morphology consistent with those of AT2 cells, and were able to differentiate into alveolar epithelial type I (AT1)-like cells. Single-cell RNA-seq of SFTPC+ cells and their progenitors demonstrated that their differentiation process and cellular heterogeneity resembled those of developing AT2 cells in vivo. AOs were applicable to drug toxicology studies recapitulating AT2-cell-specific phenotypes. Our methods can help scientists overcome the limitations of current approaches to the modeling of human alveoli and should be useful for disease modeling and regenerative medicine.

Recent advances in the detection of base modifications using the Nanopore sequencer.

DNA and RNA modifications have important functions, including the regulation of gene expression. Existing methods based on short-read sequencing for the detection of modifications show difficulty in determining the modification patterns of single chromosomes or an entire transcript sequence. Furthermore, the kinds of modifications for which detection methods are available are very limited. The Nanopore sequencer is a single-molecule, long-read sequencer that can directly sequence RNA as well as DNA. Moreover, the Nanopore sequencer detects modifications on long DNA and RNA molecules. In this review, we mainly focus on base modification detection in the DNA and RNA of mammals using the Nanopore sequencer. We summarize current studies of modifications using the Nanopore sequencer, detection tools using statistical tests or machine learning, and applications of this technology, such as analyses of open chromatin, DNA replication, and RNA metabolism.

A genome database for a Japanese population of the larvacean Oikopleura dioica.

The larvacean Oikopleura dioica is a planktonic chordate and is a tunicate that belongs to the closest relatives to vertebrates. Its simple and transparent body, invariant embryonic cell lineages, and short life cycle of 5 days make it a promising model organism for the study of developmental biology. The genome browser OikoBase was established in 2013 using Norwegian O. dioica. However, genome information for other populations is not available, even though many researchers have studied local populations. In the present study, we sequenced using Illumina and PacBio RSII technologies the genome of O. dioica from a southwestern Japanese population that was cultured in our laboratory for 3 years. The genome of Japanese O. dioica was assembled into 576 scaffold sequences with a total length and N50 length of 56.6 and 1.5 Mb, respectively. A total of 18,743 gene models (transcript models) were predicted in the genome assembly, named OSKA2016. In addition, 19,277 non-redundant transcripts were assembled using RNA-seq data. The OSKA2016 has global sequence similarity of only 86.5% when compared with the OikoBase, highlighting the sequence difference between the two far distant O. dioica populations on the globe. The genome assembly, transcript assembly, and transcript models were incorporated into ANISEED (https://www.aniseed.cnrs.fr/) for genome browsing and BLAST searches. Mapping of reads obtained from male- or female-specific genome libraries yielded male-specific scaffolds in the OSKA2016 and revealed that over 2.6 Mb of sequence were included in the male-specific Y-region. The genome and transcriptome resources from two distinct populations will be useful datasets for developmental biology, evolutionary biology, and molecular ecology using this model organism.

A GWAS approach to find SNPs associated with salt removal in rice leaf sheath.

BACKGROUND AND AIMS: The ability for salt removal at the leaf sheath level is considered to be one of the major mechanisms associated with salt tolerance in rice. Thus, understanding the genetic control of the salt removal capacity in leaf sheaths will help improve the molecular breeding of salt-tolerant rice varieties and speed up future varietal development to increase productivity in salt-affected areas. We report a genome-wide association study (GWAS) conducted to find single nucleotide polymorphisms (SNPs) associated with salt removal in leaf sheaths of rice. METHODS: In this study, 296 accessions of a rice (Oryza sativa) diversity panel were used to identify salt removal-related traits and conduct GWAS using 36 901 SNPs. The sheath:blade ratio of Na+ and Cl- concentrations was used to determine the salt removal ability in leaf sheaths. Candidate genes were further narrowed via Gene Ontology and RNA-seq analysis to those whose putative function was likely to be associated with salt transport and were up-regulated in response to salt stress. KEY RESULTS: For the association signals of the Na+ sheath:blade ratio, significant SNPs were found only in the indica sub-population on chromosome 5. Within candidate genes found in the GWAS study, five genes were upregulated and eight genes were downregulated in the internal leaf sheath tissues in the presence of salt stress. CONCLUSIONS: These GWAS data imply that rice accessions in the indica variety group are the main source of genes and alleles associated with Na+ removal in leaf sheaths of rice under salt stress.

A systematic sequencing-based approach for microbial contaminant detection and functional inference.

BACKGROUND: Microbial contamination poses a major difficulty for successful data analysis in biological and biomedical research. Computational approaches utilizing next-generation sequencing (NGS) data offer promising diagnostics to assess the presence of contaminants. However, as host cells are often contaminated by multiple microorganisms, these approaches require careful attention to intra- and interspecies sequence similarities, which have not yet been fully addressed. RESULTS: We present a computational approach that rigorously investigates the genomic origins of sequenced reads, including those mapped to multiple species that have been discarded in previous studies. Through the analysis of large-scale synthetic and public NGS samples, we estimate that 1000-100,000 contaminating microbial reads are detected per million host reads sequenced by RNA-seq. The microbe catalog we established included Cutibacterium as a prevalent contaminant, suggesting that contamination mostly originates from the laboratory environment. Importantly, by applying a systematic method to infer the functional impact of contamination, we revealed that host-contaminant interactions cause profound changes in the host molecular landscapes, as exemplified by changes in inflammatory and apoptotic pathways during Mycoplasma infection of lymphoma cells. CONCLUSIONS: We provide a computational method for profiling microbial contamination on NGS data and suggest that sources of contamination in laboratory reagents and the experimental environment alter the molecular landscape of host cells leading to phenotypic changes. These findings reinforce the concept that precise determination of the origins and functional impacts of contamination is imperative for quality research and illustrate the usefulness of the proposed approach to comprehensively characterize contamination landscapes.

MoMI-G: modular multi-scale integrated genome graph browser.

BACKGROUND: Genome graph is an emerging approach for representing structural variants on genomes with branches. For example, representing structural variants of cancer genomes as a genome graph is more natural than representing such genomes as differences from the linear reference genome. While more and more structural variants are being identified by long-read sequencing, many of them are difficult to visualize using existing structural variants visualization tools. To this end, visualization method for large genome graphs such as human cancer genome graphs is demanded. RESULTS: We developed MOdular Multi-scale Integrated Genome graph browser, MoMI-G, a web-based genome graph browser that can visualize genome graphs with structural variants and supporting evidences such as read alignments, read depth, and annotations. This browser allows more intuitive recognition of large, nested, and potentially more complex structural variations. MoMI-G has view modules for different scales, which allow users to view the whole genome down to nucleotide-level alignments of long reads. Alignments spanning reference alleles and those spanning alternative alleles are shown in the same view. Users can customize the view, if they are not satisfied with the preset views. In addition, MoMI-G has Interval Card Deck, a feature for rapid manual inspection of hundreds of structural variants. Herein, we describe the utility of MoMI-G by using representative examples of large and nested structural variations found in two cell lines, LC-2/ad and CHM1. CONCLUSIONS: Users can inspect complex and large structural variations found by long-read analysis in large genomes such as human genomes more smoothly and more intuitively. In addition, users can easily filter out false positives by manually inspecting hundreds of identified structural variants with supporting long-read alignments and annotations in a short time. SOFTWARE AVAILABILITY: MoMI-G is freely available at https://github.com/MoMI-G/MoMI-G under the MIT license.