Long-read Sequencing Research Articles

De novo detection of somatic variants in high-quality long-read single-cell RNA sequencing data.

In cancer, genetic and transcriptomic variations generate clonal heterogeneity, leading to treatment resistance. Long-read single-cell RNA sequencing (LR scRNA-seq) has the potential to detect genetic and transcriptomic variations simultaneously. Here, we present LongSom, a computational workflow leveraging high-quality LR scRNA-seq data to call de novo somatic single-nucleotide variants (SNVs), including in mitochondria (mtSNVs), copy-number alterations (CNAs), and gene fusions, to reconstruct the tumor clonal heterogeneity. Before somatic variants calling, LongSom re-annotates marker gene based cell types using cell mutational profiles. LongSom distinguishes somatic SNVs from noise and germline polymorphisms by applying an extensive set of hard filters and statistical tests. Applying LongSom to human ovarian cancer samples, we detected clinically relevant somatic SNVs that were validated against matched DNA samples. Leveraging somatic SNVs and fusions, LongSom found subclones with different predicted treatment outcomes. In summary, LongSom enables de novo variant detection without the need for normal samples, facilitating the study of cancer evolution, clonal heterogeneity, and treatment resistance.

Factors impacting target-enriched long-read sequencing of resistomes and mobilomes.

We investigated the efficiency of target-enriched long-read sequencing (TELSeq) for detecting antimicrobial resistance genes (ARGs) and mobile genetic elements (MGEs) within complex matrices. We aimed to overcome limitations associated with traditional antimicrobial resistance (AMR) detection methods, including short-read shotgun metagenomics, which can lack sensitivity, specificity, and the ability to provide detailed genomic context. By combining biotinylated probe-based enrichment with long-read sequencing, we facilitated the amplification and sequencing of ARGs, eliminating the need for bioinformatic reconstruction. Our experimental design included replicates of human fecal microbiota transplant material, bovine feces, pristine prairie soil, and a mock human gut microbial community, allowing us to examine variables including genomic DNA input and probe set composition. Our findings demonstrated that TELSeq markedly improves the detection rates of ARGs and MGEs compared to traditional sequencing methods, underlining its potential for accurate AMR monitoring. A key insight from our research is the importance of incorporating mobilome profiles to better predict the transferability of ARGs within microbial communities, prompting a recommendation for the use of combined ARG-MGE probe sets for future studies. We also reveal limitations for ARG detection from low-input workflows, and describe the next steps for ongoing protocol refinement to minimize technical variability and expand utility in clinical and public health settings. This effort is part of our broader commitment to advancing methodologies that address the global challenge of AMR.

Best practices for germline variant and DNA methylation analysis of second- and third-generation sequencing data.

This comprehensive review provides insights and suggested strategies for the analysis of germline variants using second- and third-generation sequencing technologies (SGS and TGS). It addresses the critical stages of data processing, starting from alignment and preprocessing to quality control, variant calling, and the removal of artifacts. The document emphasized the importance of meticulous data handling, highlighting advanced methodologies for annotating variants and identifying structural variations and methylated DNA sites. Special attention is given to the inspection of problematic variants, a step that is crucial for ensuring the accuracy of the analysis, particularly in clinical settings where genetic diagnostics can inform patient care. Additionally, the document covers the use of various bioinformatics tools and software that enhance the precision and reliability of these analyses. It outlines best practices for the annotation of variants, including considerations for problematic genetic alterations such as those in the human leukocyte antigen region, runs of homozygosity, and mitochondrial DNA alterations. The document also explores the complexities associated with identifying structural variants and copy number variations, underscoring the challenges posed by these large-scale genomic alterations. The objective is to offer a comprehensive framework for researchers and clinicians, ensuring that genetic analyses conducted with SGS and TGS are both accurate and reproducible. By following these best practices, the document aims to increase the diagnostic accuracy for hereditary diseases, facilitating early diagnosis, prevention, and personalized treatment strategies. This review serves as a valuable resource for both novices and experts in the field, providing insights into the latest advancements and methodologies in genetic analysis. It also aims to encourage the adoption of these practices in diverse research and clinical contexts, promoting consistency and reliability across studies.

The quality and detection limits of mitochondrial heteroplasmy by long read nanopore sequencing

This study evaluates long-read and short-read sequencing for mitochondrial DNA (mtDNA) heteroplasmy detection. 592,315 bootstrapped datasets generated from two single-nucleotide mismatched ultra-deep sequenced mtDNA samples were used to assess basecalling error and accuracy, limit of heteroplasmy detection, and heteroplasmy detection across various coverage depths. Results showed high Phred scores of data with GC-rich sequence bias for long reads. Limit of detection of 12% heteroplasmy was identified, showing strong correlation (R2 ≥ 0.955) with expected heteroplasmy but underreporting tendency of high-level variants. Nanopore sequencing shows potential for direct applicability in mitochondrial diseases diagnostics, but stringent validation processes to ensure diagnostic result quality are required.

ITSxpress version 2: software to rapidly trim internal transcribed spacer sequences with quality scores for amplicon sequencing.

The nuclear ribosomal RNA (rRNA) internal transcribed spacer (ITS) regions are commonly used to identify fungi and other eukaryotic taxa in amplicon sequencing. The highly conserved rRNA regions flanking the ITS are often trimmed before being used for taxonomic assignment. The Python software package ITSxpress rapidly trims single-end or paired-end sequences in FASTQ format for use in amplicon sequence variant clustering methods like DADA2. This new major release of ITSxpress improves the paired-end merging method, simplifies installation of the QIIME 2 ITSxpress plugin, removes major dependencies, adds use cases, and is compatible with newer compression formats. This article discusses the modifications to ITSxpress that improve the output and user experience, leading to a major version increase.IMPORTANCEITSxpress is a sequence trimming method applied to internal transcribed spacer (ITS) amplicon sequences before calling amplicon sequence variants (ASVs). The ITS region is used to understand the composition of eukaryotic microbial communities. ITS sequences provide good taxonomic resolution due to their hypervariability, but are flanked by conserved regions that allow their primers to be more universal. Amplicons generated with such primers contain regions with different evolutionary rates, and trimming these conserved regions results in better taxonomic classification and a more valid set of ASVs. This package can be used for most amplicon sequencing methods including for newer long-read sequencing formats, such as PacBio. ITSxpress can be installed from Bioconda, used as a Docker image, or installed from source code. The package works well with high-performance computing clusters or laptops due to its low-resource requirements.

Integrating short- and full-length 16S rRNA gene sequencing to elucidate microbiome profiles in Pacific white shrimp (Litopenaeus vannamei) ponds.

Despite their immense economic value as a key aquaculture species, the production of Pacific white shrimp (Litopenaeus vannamei) faces significant challenges from intensive farming practices and disease outbreaks. Routine microbial profiling for disease surveillance could be a promising approach to anticipate and control disease outbreaks. To achieve this, accuracy in microbial profiling in shrimp ponds is crucial for enabling targeted action and prevention. Extensive documentation emphasizes that, beyond biological factors (related to the host, diet, or health status during the rearing period), technical elements, including sequencing techniques significantly influence bacterial community profiling. This study investigated the influence of short- and long-read sequencing of 16S rRNA genes on the microbial profiles in shrimp intestines, water, and sediments. The origin of the samples (intestine or environmental) in shrimp culture ponds primarily drove the observed differences in core microbial species. The ecological niches accounted for 56% of bacterial community variations in culture ponds. Both sequencing approaches showed consistent results in identifying higher-rank taxa and assessing alpha and beta diversity. However, at the species level, full-length 16S rRNA gene sequences provided better resolution than V3-V4 sequences. For routine microbial profiling in shrimp culture ponds, our study suggests that short-read sequences were sufficient for determining overall bacterial community.IMPORTANCEThis interdisciplinary study investigated the influence of sequencing techniques on bacterial communities profiling within Pacific white shrimp (Litopenaeus vannamei) ponds. By integrating aquaculture, microbiology, and environmental science, we revealed the role of ecological niches and factors like salinity and pH on microbiota diversity and composition in shrimp intestines, pond water, and sediment. Additionally, we compared the taxonomic resolution using partial versus full-length 16S rRNA gene sequences, highlighting the value of longer amplicons for precise identification of key taxa. These findings provide novel insights into microbial dynamics underlying environmental effects in shrimp aquaculture. Comprehensive characterization of the pond microbiome could lead to management strategies that promote shrimp health and productivity. Furthermore, the potential of a multi-omics approach for integrating complementary data streams to elucidate environment-microbiome-host interactions was highlighted.

A chromosome-level genome assembly of a model conifer plant, the Japanese cedar, Cryptomeria japonica D. Don.

The Japanese cedar (Cryptomeria japonica D. Don) is one of the most important Japanese forest trees, occupying approximately 44% of artificial forests and planted in East Asia, the Azores Archipelago, and certain islands in the Indian Ocean. Although the huge genome of the species (ca. 9 Gbp) with abundant repeat elements may have represented an obstacle for genetic analysis, this species is easily propagated by cutting, flowered by gibberellic acid, transformed by Agrobacterium, and edited by CRISPR/Cas9. These characteristics of C. japonica recommend it as a model conifer species for which reference genome sequences are necessary. Herein, we report the first chromosome-level assembly of C. japonica (2n = 22) using third-generation selfed progeny (estimated homozygosity rate = 0.96). Young leaf tissue was used to extract high molecular weight DNA (> 50kb) for HiFi PacBio long-read sequencing and to construct an Hi-C/Omni-C library for Illumina short-read sequencing. The 29× and 26× genome coverage of HiFi and Illumina reads, respectively, for de novo assembly yielded 2,651 contigs (9.1 Gbp, N50 contig size 12.0 Mbp). Hi-C analysis mapped 97% of the nucleotides on 11 chromosomes. The assembly was verified through comparison with a consensus linkage map comprising 7,781 markers. BUSCO analysis identified ∼ 91% conserved genes. Annotations of genes and comparisons of repeat elements with other Cupressaceae and Pinaceae species provide a fundamental resource for conifer research.

Comparison of Four Different DNA Isolation Methods from MGIT Culture for Long-Read Whole Genome Sequencing of Mycobacterium tuberculosis

Background: Tuberculosis (TB) remains a global health challenge, particularly due to drug resistance and limitations in rapid diagnosis. Next-generation sequencing (NGS), especially long-read whole genome sequencing (WGS), shows promise for rapidly detecting TB and drug resistance, but it requires high-quality DNA, which is difficult to extract from Mycobacterium tuberculosis due to its complex cell wall. Objectives: This study evaluated four DNA isolation methods for extracting pure DNA from M. tuberculosis, aiming to standardize protocols for long-read WGS. Methods: Mycobacterium tuberculosis H37RV colonies were grown in BACTEC MGIT liquid medium. Two pellets were prepared as the initial material for the DNA extraction protocol: Pellets from 1 mL McFarland 2 suspensions and all growing colonies from two MGIT liquid cultures. Four DNA extraction methods were used: The cetyltrimethylammonium bromide (CTAB) method, GeneJET Genomic DNA Purification Kit, Quick-DNA Fecal/Soil Microbe Kit, and Genematrix Tissue/Bacterial DNA Purification Kit, with some modifications. DNA quality was assessed based on concentration, purity, and integrity. Results: Among the tested methods, the Quick-DNA Fecal/Soil Kit yielded approximately 85 ng/mL of DNA and a purity of 1.9 at 260/280 nm from the colonial pellet of two MGIT tubes. However, lower intact DNA [DNA integrity number (DIN) ~ 6.8] was obtained with this kit. The CTAB method provided the highest intact DNA (DIN ~ 9.5), although the purity of the DNA was not sufficient. Conclusions: Based on three repetitions of McF-2 and colonial pellet extractions, the Quick-DNA Fecal/Soil Kit yielded the highest DNA quantity and purity but showed lower integrity compared to other methods, indicating the need for adjustments. A pellet from two MGIT cultures (~ 100 µL) is suitable for long-read WGS with this kit. However, a larger sample size is required to generalize these findings. For effective long-read sequencing of M. tuberculosis, DNA extraction protocols must be optimized to balance yield, fragment size, and purity for accurate sequencing and drug resistance analysis.

DmTGS: Precise Targeted Enrichment Long-Read Sequencing Panel for Tandem Repeat Detection.

Tandem repeats (TRs) are abundant in the human genome and associated with repeat expansion disorders. Our study aimed to develop a tandem repeat panel utilizing targeted long-read sequencing to evaluate known TRs associated with these disorders and assess its clinical utility. We developed a targeted long-read sequencing panel for 70 TR loci, termed dynamic mutation third-generation sequencing (dmTGS), using the PacBio Sequel II platform. We tested 108 samples with suspected repeat expansion disorders and compared the results with conventional molecular methods. For 108 samples, dmTGS achieved an average of 8000 high-fidelity reads per sample, with a mean read length of 4.7 kb and read quality of 99.9%. dmTGS outperformed repeat-primed-PCR and fluorescence amplicon length analysis-PCR in distinguishing expanded from normal alleles and accurately quantifying repeat counts. The method demonstrated high concordance with confirmatory methods (rlinear = 0.991, P < 0.01), and detected mosaicism with sensitivities of 1% for FMR1 CGG premutation and 5% for full mutations. dmTGS successfully identified interruptive motifs in genes that conventional methods had missed. For variable number TRs in the PLIN4 gene, dmTGS identified precise repeat counts and sequence motifs. Screening 57 patients with suspected genetic muscular diseases, dmTGS confirmed repeat expansions in genes such as GIPC1, NOTCH2NLC, NUTM2B-AS1/LOC642361, and DMPK. Additionally, dmTGS detected CCG interruptions in CTG repeats in 8 myotonic dystrophy type 1 patients with detailed characterization. dmTGS accurately detects repeat sizes and interruption motifs associated with repeat expansion disorders and demonstrates superior performance compared to conventional molecular methods.

Identification and annotation of centromeric hypomethylated regions with Centromere Dip Region (CDR)-Finder.

Centromeres are chromosomal regions historically understudied with sequencing technologies due to their repetitive nature and short-read mapping limitations. However, recent improvements in long-read sequencing allowed for the investigation of complex regions of the genome at the sequence and epigenetic levels. Here, we present Centromere Dip Region (CDR)-Finder: a tool to identify regions of hypomethylation within the centromeres of high-quality, contiguous genome assemblies. These regions are typically associated with a unique type of chromatin containing the histone H3 variant CENP-A, which marks the location of the kinetochore. CDR-Finder identifies the CDRs in large and short centromeres and generates a BED file indicating the location of the CDRs within the centromere. It also outputs a plot for visualization, validation, and downstream analysis. CDR-Finder is available at https://github.com/EichlerLab/CDR-Finder.

Fully phased genome assemblies and graph-based genetic variants of the olive flounder, Paralichthys olivaceus

The olive flounder, Paralichthys olivaceus, also known as the Korean halibut, is an economically important flatfish in East Asian countries. Here, we provided four fully phased genome assemblies of two different olive flounder individuals using high-fidelity long-read sequencing and their parental short-read sequencing data. We obtained 42–44 Gb of ~15-kb and ~Q30 high-fidelity long reads, and their assembly quality values were ~53. We annotated ~30 K genes, ~170-Mb repetitive sequences, and ~3 M 5-methylcytosine positions for each genome assembly, and established a graph-based draft pan-genome of the olive flounder. We identified 5 M single-nucleotide variants and 100 K structural variants with their genotype information, where ~13% of the variants were possibly fixed in the two Korean individuals. Based on our chromosome-level genome assembly, we also explored chromosome evolution in the Pleuronectiformes family, as reported earlier. Our high-quality genomic resources will contribute to future genomic selection for accelerating the breeding process of the olive flounder.

MAGqual: a stand-alone pipeline to assess the quality of metagenome-assembled genomes

BackgroundMetagenomics, the whole genome sequencing of microbial communities, has provided insight into complex ecosystems. It has facilitated the discovery of novel microorganisms, explained community interactions and found applications in various fields. Advances in high-throughput and third-generation sequencing technologies have further fuelled its popularity. Nevertheless, managing the vast data produced and addressing variable dataset quality remain ongoing challenges. Another challenge arises from the number of assembly and binning strategies used across studies. Comparing datasets and analysis tools is complex as it requires the quantitative assessment of metagenome quality. The inherent limitations of metagenomic sequencing, which often involves sequencing complex communities, mean community members are challenging to interrogate with traditional culturing methods leading to many lacking reference sequences. MIMAG standards aim to provide a method to assess metagenome quality for comparison but have not been widely adopted.ResultsTo address the need for simple and quick metagenome quality assignation, here we introduce the pipeline MAGqual (Metagenome-Assembled Genome qualifier) and demonstrate its effectiveness at determining metagenomic dataset quality in the context of the MIMAG standards.ConclusionsThe MAGqual pipeline offers an accessible way to evaluate metagenome quality and generate metadata on a large scale. MAGqual is built in Snakemake to ensure readability and scalability, and its open-source nature promotes accessibility, community development, and ease of updates. MAGqual is built in Snakemake, R, and Python and is available under the MIT license on GitHub at https://github.com/ac1513/MAGqual.EB_9RbKd3CYdEq3zCTHKaqVideo

Decoding the Genetic Code: Scientific Exploration of the Telomere-totelomere (T2T) Genome

Abstract: With the development of third-generation sequencing technology, genome assembly has entered a new era. The combination of multiple sequencing methods can combine the strengths of each, resulting in higher-quality assembly results. Telomere-to-Telomere (T2T) genome refers to a zero gap genome that is assembled at the T2T level of one or more chromosomes through the combination of multiple sequencing technologies, such as Pacific Biosciences High-Fidelity (PacBio HiFi), Oxford Nanopore Technologies (ONT) Ultra-long, High-throughput Chromosome Conformation Capture (Hi-C), and others. High-quality reference genomes are the basis for genomics research, and the T2T genome has enabled the exploration of unknown areas of the genome, such as telomeres and centromeres, which has provided a more in-depth direction for research. This review provides a comprehensive overview of the T2T genome, reviewing the development of sequencing technologies and then outlining sequencing strategies, assembly methods, quality assessment processes, and analysis software for the T2T genome with practical examples. Representative T2T genomic species of plants, animals, and microorganisms (e.g., human, Arabidopsis, brewer's yeast, and so on) are presented separately. A summary of the potential and challenges of current T2T genomic applications is provided, along with an outlook for future developments.

Long-Read Sequencing Identifying the Genetic Complexity of Congenital Adrenal Hyperplasia in the Pedigree.

High sequence homology between CYP21A2 and CYP21A1P poses challenges to genetic diagnosis of congenital adrenal hyperplasia (CAH) due to 21-hydroxylase deficiency (21-OHD). Traditional genetic testing is unable to provide an accurate diagnosis due to the genetic complexity of CAH. Deletions, duplications, and recombination breakpoints were precisely identified by long-read sequencing (LRS). This study presented a pregnant woman, a 21-OHD carrier detected by MLPA, and her husband, a normal subject also detected by MLPA. The fetus was suspected of having 21-OHD based on clinical presentations such as enlarged adrenal glands, atypical external genitalia and karyotyping of 46, XX. LRS further identified the fetus as having the most severe salt-wasting (SW) form of 21-OHD with a compound heterozygote genotype. One allele was TNXA/TNXB CH-2, while the other allele was CYP21A1P/CYP21A2 CH-8. LRS precisely determined the genotypes of the fetus's father and grandmother with duplications, which misdiagnosed by MLPA. The multidisciplinary team recommended immediate glucocorticoid and mineralocorticoid treatment for the child after birth to prevent life-threatening adrenal crisis. LRS provides precise diagnosis for family members with CYP21A2 deletion or duplication, improving disease management and preventing potential adrenal crises. When used in pre-pregnancy genetic testing, LRS can indicate high genetic risk and guide the appropriate therapy during pregnancy and immediately after birth.

Identification of complex Plasmodium falciparum genetic backgrounds circulating in Africa: a multicountry genomic epidemiology analysis

The population structure of the malaria parasite Plasmodium falciparum can reveal underlying adaptive evolutionary processes. Selective pressures to maintain complex genetic backgrounds can encourage inbreeding, producing distinct parasite clusters identifiable by population structure analyses. We analysed population structure in 3783 Pfalciparum genomes from 21countries across Africa, provided by the MalariaGEN Pf7 dataset. We used Principal Coordinate Analysis to cluster parasites, identity by descent (IBD) methods to identify genomic regions shared by cluster members, and linkage analyses to establish their co-inheritance patterns. Structural variants were reconstructed by de novo assembly and verified by long-read sequencing. We identified a strongly differentiated cluster of parasites, named AF1, comprising 47 (1·2%) of 3783samples analysed, distributed over 13countries across Africa, at locations over 7000km apart. Members of this cluster share a complex genetic background, consisting of up to 23loci harbouring many highly differentiated variants, rarely observed outside the cluster. IBD analyses revealed common ancestry at these loci, irrespective of sampling location. Outside the shared loci, however, AF1 members appear to outbreed with sympatric parasites. The AF1 differentiated variants comprise structural variations, including a gene conversion involving the dblmsp and dblmsp2 genes, and numerous single nucleotide polymorphisms. Several of the genes harbouring these mutations are functionally related, often involved in interactions with red blood cells including invasion, egress, and erythrocyte antigen export. We propose that AF1 parasites have adapted to some unidentified evolutionary niche, probably involving interactions with host erythrocytes. This adaptation involves a complex compendium of interacting variants that are rarely observed in Africa, which remains mostly intact despite recombination events. The term cryptotype was used to describe a common background interspersed with genomic regions of local origin. Bill & Melinda Gates Foundation.

Genetic Variation in Jamaican Populations of the Coffee Berry Borer, Hypothenemus hampei.

The coffee berry borer (CBB) Hypothenemus hampei was first described in Africa in 1867 and has spread to all major coffee-producing regions worldwide, including Jamaica. Using long-read sequencing, we produced a new high-quality reference genome (172.7 Mb) for the Jamaican strain of the CBB, with 93% of the genome assembled into 14 scaffolds. Whole genome sequencing of pooled samples from different populations across Jamaica showed that the CBB harbors low levels of genetic diversity alongside an excess of low-frequency alleles, indicative of a recent genetic bottleneck. The analyses also showed a recent surge in the activity of transposable elements (TEs), particularly LINE/R1 and LTR/Gypsy elements, within CBB populations. Our findings offer first insights into the evolutionary genomics of CBB populations in Jamaica, highlighting the potential role of TEs in shaping the genome of this important pest species.

The HLA-A Allele, HLA-A*11:01:01:07, Identified by Third-Generation Sequencing.

HLA-A*11:01:01:07 differs from HLA-A*11:01:01:01 by a single nucleotide exchange in intron 1.

Diversity of the Obligate Gut Bacteria Indicates Host-Symbiont Coevolution at the Population Level in the Plataspid Stinkbug Megacopta cribraria.

Ishikawaella is an obligate gut bacterium in stinkbugs that belong to Plataspidae family (Hemiptera: Heteroptera). It is vertically transmitted to newborn nymphs through capsules laid on eggs by maternal stinkbugs. Previous research has established a pattern of strict cospeciation between Plataspidae stinkbugs and Ishikawaella. However, the possibility of host-symbiont coevolution at the population level within Plataspidae stinkbugs has not been thoroughly explored. This study analyzed the samples of Megacopta cribraria from three phylogenetic clades to investigate host-symbiont coevolution in this insect species. We compared data from third-generation sequencing (PacBio), next-generation sequencing (Illumina), and first-generation sequencing (Sanger), and the results indicated that Illumina sequencing most accurately represents the composition of gut bacterial communities. All stinkbug individuals shared a dominant amplicon sequence variant (ASV), which accounted for an average of 65.99% of Ishikawaella sequences (ranging from 58.68% to 87.01%). The top five ASVs (ASV0-ASV4) represented 99.82% of all Ishikawaella sequences. Among these, the number of base substitutions between any two ASVs ranged from 1 to 3, significantly lower than the number of substitutions between the main and minor ASVs. This finding suggests that closely related strains are likely to coexist in the same host. Beta diversity analyses revealed significant differences in Ishikawaella composition among the three phylogenetic clades, providing evidence for host-symbiont coevolution at the population level in Plataspidae stinkbugs.

Genome sequencing and transcriptome analysis provide an insight into the hypoxia resistance of Channa asiatica

Channa asiatica is an economically valuable fish species and excellent model for studying hypoxic tolerance. However, the underlying genetic and molecular mechanisms are poorly understood. In this study, we assembled a high-quality C. asiatica genome (23 chromosomes, totaling 722 Mb) using a combination of Illumina short-read, PacBio long-read, and Hi-C sequencing. Repetitive elements accounted for 28.39%of the C. asiatica genome, and 23,949 protein-coding genes were predicted, with 96.63 % of these functionally annotated. Moreover, a comparative genomic analysis of 12 fish genomes showed that gene families associated with oxygen binding and transport were expanded in C. asiatica. In addition, transcriptome analysis revealed that multiple oxidative stress pathways were activated when C. asiatica was exposed to air. In conclusion, this study provided high-quality genome assembly and transcriptome data, both serving as critical resources for researching the genetic basis of hypoxic tolerance in C. asiatica.

Leveraging the T2T assembly to resolve rare and pathogenic inversions in reference genome gaps.

Chromosomal inversions (INVs) are particularly challenging to detect due to their copy-number neutral state and association with repetitive regions. Inversions represent about 1/20 of all balanced structural chromosome aberrations and can lead to disease by gene disruption or altering regulatory regions of dosage-sensitive genes in cis Short-read genome sequencing (srGS) can only resolve ∼70% of cytogenetically visible inversions referred to clinical diagnostic laboratories, likely due to breakpoints in repetitive regions. Here, we study 12 inversions by long-read genome sequencing (lrGS) (n = 9) or srGS (n = 3) and resolve nine of them. In four cases, the inversion breakpoint region was missing from at least one of the human reference genomes (GRCh37, GRCh38, T2T-CHM13) and a reference agnostic analysis was needed. One of these cases, an INV9 mappable only in de novo assembled lrGS data using T2T-CHM13 disrupts EHMT1 consistent with a Mendelian diagnosis (Kleefstra syndrome 1; MIM#610253). Next, by pairwise comparison between T2T-CHM13, GRCh37, and GRCh38, as well as the chimpanzee and bonobo, we show that hundreds of megabases of sequence are missing from at least one human reference, highlighting that primate genomes contribute to genomic diversity. Aligning population genomic data to these regions indicated that these regions are variable between individuals. Our analysis emphasizes that T2T-CHM13 is necessary to maximize the value of lrGS for optimal inversion detection in clinical diagnostics. These results highlight the importance of leveraging diverse and comprehensive reference genomes to resolve unsolved molecular cases in rare diseases.