Long-read Whole-genome Sequencing Research Articles

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.

The interplay of mutagenesis and ecDNA shapes urothelial cancer evolution

Advanced urothelial cancer is a frequently lethal disease characterized by marked genetic heterogeneity1. In this study, we investigated the evolution of genomic signatures caused by endogenous and external mutagenic processes and their interplay with complex structural variants (SVs). We superimposed mutational signatures and phylogenetic analyses of matched serial tumours from patients with urothelial cancer to define the evolutionary dynamics of these processes. We show that APOBEC3-induced mutations are clonal and early, whereas chemotherapy induces mutational bursts of hundreds of late subclonal mutations. Using a genome graph computational tool2, we observed frequent high copy-number circular amplicons characteristic of extrachromosomal DNA (ecDNA)-forming SVs. We characterized the distinct temporal patterns of APOBEC3-induced and chemotherapy-induced mutations within ecDNA-forming SVs, gaining new insights into the timing of these mutagenic processes relative to ecDNA biogenesis. We discovered that most CCND1 amplifications in urothelial cancer arise within circular ecDNA-forming SVs. ecDNA-forming SVs persisted and increased in complexity, incorporating additional DNA segments and contributing to the evolution of treatment resistance. Oxford Nanopore Technologies long-read whole-genome sequencing followed by de novo assembly mapped out CCND1 ecDNA structure. Experimental modelling of CCND1 ecDNA confirmed its role as a driver of treatment resistance. Our findings define fundamental mechanisms that drive urothelial cancer evolution and have important therapeutic implications.

Genomics Insights into Mycolicibacterium Hassiacum Causing Infection in a Cat with Pyogranulomatous Dermatitis and Panniculitis.

Mycolicibacterium hassiacum (homotypic synonym: Mycobacterium hassiacum) represents an ungrouped thermotolerant rapidly growing mycobacteria (RGM) species occasionally associated with infections and disease in humans. In this report, we describe a case of pyogranulomatous dermatitis and panniculitis due to M. hassiacum in an immunocompetent adult cat. To the best of our knowledge, this represents the first report of M. hassiacum infection in animals. We also report the results of the in-depth genome characterization of the isolate using a combined short- and long-read whole-genome sequencing (WGS) approach. We observed the lack of acquired-resistance genes and no evidence of mutations in housekeeping genes associated with resistance to rifampicin and isoniazid. We detected some virulence factors in our isolate, such as some associated with the interaction of mycobacteria with host cells, and the presence of multiple copies of heavy metal resistance genes (arsB, arsR, and arsL/cadL). In conclusion, M. hassiacum should be included among the RGM species associated with feline subcutaneous atypical mycobacteriosis (SAM). A reliable and fast RGM laboratory identification and characterization is important not only for an accurate etiological diagnosis but also for a correct approach to SAM treatment options.

PQEB1: a hospital outbreak plasmid lineage carrying bla KPC-2.

While conducting genomic surveillance of carbapenemase-producing Enterobacteriaceae (CPE) from patient colonisation and clinical infections at Birmingham's Queen Elizabeth Hospital (QE), we identified an N-type plasmid lineage, pQEB1, carrying several antibiotic resistance genes, including the carbapenemase gene bla KPC-2. The pQEB1 lineage is concerning due to its conferral of multidrug resistance, its host range and apparent transmissibility, and its potential for acquiring further resistance genes. Representatives of pQEB1 were found in three sequence types (STs) of Citrobacter freundii, two STs of Enterobacter cloacae, and three species of Klebsiella. Hosts of pQEB1 were isolated from 11 different patients who stayed in various wards throughout the hospital complex over a 13 month period from January 2023 to February 2024. At present, the only representatives of the pQEB1 lineage in GenBank were carried by an Enterobacter hormaechei isolated from a blood sample at the QE in 2016 and a Klebsiella pneumoniae isolated from a urine sample at University Hospitals Coventry and Warwickshire (UHCW) in May 2023. The UHCW patient had been treated at the QE. Long-read whole-genome sequencing was performed on Oxford Nanopore R10.4.1 flow cells, facilitating comparison of complete plasmid sequences. We identified structural variants of pQEB1 and defined the molecular events responsible for them. These have included IS26-mediated inversions and acquisitions of multiple insertion sequences and transposons, including carriers of mercury or arsenic resistance genes. We found that a particular inversion variant of pQEB1 was strongly associated with the QE Liver speciality after appearing in November 2023, but was found in different specialities and wards in January/February 2024. That variant has so far been seen in five different bacterial hosts from six patients, consistent with recent and ongoing inter-host and inter-patient transmission of pQEB1 in this hospital setting.

Reciprocal chromosome translocation t(3;4)(q27;q31.2) with deletion of 3q27 and reduced FBXW7 expression in a patient with developmental delay, hypotonia, and seizures.

Reciprocal chromosomal translocation is one of genomic variations. When cytogenetically de novo reciprocal translocations are identified in patients with some clinical manifestations, the genes in the breakpoints are considered to be related to the clinical features. In this study, we encountered a patient with severe developmental delay, intractable epilepsy, growth failure, distinctive features, and skeletal manifestations. Conventional karyotyping revealed a de novo translocation described as 46,XY,t(3;4)(q27;q31.2). Chromosomal microarray testing detected a 1.25-Mb microdeletion at 3q27.3q28. Although the skeletal manifestations may have been affected by this deletion, the neurological features of this patient were severe and could not be fully explained by this deletion. Since no genomic copy number aberration was detected on chromosome 4, long-read whole-genome sequencing analysis was performed and a precise breakpoint was confirmed. A 460-bp deletion was detected between the two breakpoints; however, no gene was disrupted. FBXW7, the gene responsible for developmental delay, hypotonia, and impaired language, is in the 0.5-Mb telomeric region. Most of the patient's clinical features were considered consistent with symptoms of FBXW7-related disorders, but were more severe. FBXW7 expression in the immortalized lymphoblasts of the patient was reduced compared to that in controls. Based on these findings, we suspect that FBXW7 is affected by downstream position effects of chromosomal translocations. The severe neurological features of the patient may have been affected not only by the 3q27-q28 deletion but also by impaired expression of FBXW7 derived from the breakage of chromosome 4.

A family with neuronal intranuclear inclusion disease with focal segmental glomerulosclerosis.

Neuronal intranuclear inclusion disease (NIID) is a rare neurodegenerative disease caused by the expansion of GGC repeats in the 5'-untranslated region (5'-UTR) of NOTCH2NLC. Although increasing evidence suggests that NIID affects various organs, its association with renal involvement remains unclear. We studied the genetic background of a family with NIID, in which four of five members presented with proteinuria as the initial manifestation. The renal pathology of three patients was diagnosed as focal segmental glomerulosclerosis (FSGS) at a previous hospital. These patients also presented with tremors, retinal degeneration, and episodic neurological events. Finally, one patient exhibited reversible bilateral thalamic high-intensity signal changes on diffusion-weighted imaging during episodic neurological events. Exome sequencing (ES) and nanopore long-read whole-genome sequencing (LR-WGS) were performed on the index case, followed by nanopore target sequencing using Cas9-mediated PCR-free enrichment and methylation analysis. ES revealed no candidate variants; however, nanopore LR-WGS in the index case revealed expansion of short tandem repeats (STR) in NOTCH2NLC. Subsequent nanopore target sequencing using Cas9-mediated PCR-free enrichment showed STR expansion of NOTCH2NLC in an affected sibling and asymptomatic father. Methylation analysis using nanopore data revealed hypermethylation of the expanded allele in the asymptomatic father and partial hypermethylation in a mildly symptomatic sibling, whereas the expanded allele was hypomethylated in the index case. This investigation expands the clinical spectrum of NIID, suggesting that STR expansion of NOTCH2NLC is a cause of renal diseases, including FSGS.

Comprehensive genetic analysis of facioscapulohumeral muscular dystrophy by Nanopore long-read whole-genome sequencing

BackgroundFacioscapulohumeral muscular dystrophy (FSHD) is a high-prevalence autosomal dominant neuromuscular disease characterized by significant clinical and genetic heterogeneity. Genetic diagnosis of FSHD remains a challenge because it cannot be detected by standard sequencing methods and requires a complex diagnosis workflow.MethodsWe developed a comprehensive genetic FSHD detection method based on Oxford Nanopore Technologies (ONT) whole-genome sequencing. Using a case–control design, we applied this procedure to 29 samples and compared the results with those from optical genome mapping (OGM), bisulfite sequencing (BSS), and whole-exome sequencing (WES).ResultsUsing our ONT-based method, we identified 59 haplotypes (35 4qA and 24 4qB) among the 29 samples (including a mosaic sample), as well as the number of D4Z4 repeat units (RUs). The pathogenetic D4Z4 RU contraction identified by our ONT-based method showed 100% concordance with OGM results. The methylation levels of the most distal D4Z4 RU and the double homeobox 4 gene (DUX4) detected by ONT sequencing are highly consistent with the BSS results and showed excellent diagnostic efficiency. Additionally, our ONT-based method provided an independent methylation profile analysis of two permissive 4qA alleles, reflecting a more accurate scenario than traditional BSS. The ONT-based method detected 17 variations in three FSHD2-related genes from nine samples, showing 100% concordance with WES.ConclusionsOur ONT-based FSHD detection method is a comprehensive method for identifying pathogenetic D4Z4 RU contractions, methylation level alterations, allele-specific methylation of two 4qA haplotypes, and variations in FSHD2-related genes, which will all greatly improve genetic testing for FSHD.

Abstract 6247: Comprehensive characterization of cancers with epithelial mesenchymal transition phenotype through long-read sequencing

Abstract Introduction: Pancreatic ductal adenocarcinoma (PDAC) and linitis-plastica type gastric cancer (STAD) are the most representative refractory cancers even in the era of molecularly targeted therapy and immunotherapy. Extensive characterization of these tumors, especially regarding the genomic and epigenomic characteristics of pure cancer cells, has been hampered by the high stromal content in the primary tumor samples. Here we took advantage of cancer cells purified from malignant ascites and established corresponding cell lines to overcome the limitations associated with the analysis and to discover their vulnerabilities. Methods: We developed cell lines from malignant ascitic fluid samples obtained from 60 patients with PDAC and 59 patients with linitis-plastica type STAD. We performed comprehensive multi-omic analyses of these cell lines using long-read whole-genome sequencing (WGS) for complex somatic genomic alterations and for methylation aberration. We also conducted long-read RNA sequencing (iso-seq) for full-length isoform exploration, and chromatin immunoprecipitation coupled with sequencing (ChIP-seq) for enhancer analysis. Complementary short-read whole-genome sequencing and RNA sequencing were also performed. Results: Comprehensive transcriptome analyses identified a cluster with high epithelial-mesenchymal transition (EMT) characteristics, consisting of the majority of the PDAC cell lines and a half of the STAD cell lines. ChIP-seq analysis against H3K27ac revealed the distinct transcriptional circuit that generated EMT cluster-specific super-enhancers governed by SMAD3, independent of cancer type. Most of the transcriptional activity was controlled by the gene methylation. We also identified the differential expression levels of the splicing factors between the EMT-activated and the non-activated groups. Comprehensive isoform analysis by iso-seq revealed that the isoform usage was completely different between the EMT group and the other, especially in genes of the actomyosin pathway which confer characteristic features of the EMT cluster in the cell morphology and cell motility. Using the CHM13 reference genome, long-read WGS analyses revealed the recurrent structural alterations specific to the EMT group such as those involving CDKN2A. We also identified PDAC-specific structural alterations. Conclusion: We identified the distinct molecular characteristics of the EMT-activated group in PDAC and linitis-plastica type STAD which may confer the therapeutic resistance. Our findings provide an interesting insight for the further innovation of the treatment strategy for these intractable cancers. Citation Format: Yosuke Tanaka, Hiroyuki Mano. Comprehensive characterization of cancers with epithelial mesenchymal transition phenotype through long-read sequencing [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 6247.

Abstract 1751: Deep multi-omics analysis of uterine leiomyosarcoma reveals defective HRD signature, a novel gene fusion and the amplification of the poor-prognosis CTHRC1 gene

Abstract Uterine leiomyosarcoma (uLMS) is a rare and aggressive cancer representing approximately 2-5% of all uterine malignancies. The molecular heterogeneity and pathogenesis of uLMS are notwell understood. To unravel the intricate molecular landscape of uLMS, we conducted a multi-omics workflow consisting of 4 high-throughput sequencing technologies; short- and long-readwhole genome sequencing (WGS), RNA-Seq, and Tandem mass tag (TMT) labeling enablingmassively parallel proteome proteomic analyses. Our study cohort consisted of 97 fresh frozentissue samples collected from 68 female patients, with 19 uLMS, 37 leiomyoma (LM), and 41normal myometrium (MM). A deep analysis of somatic variants with clinical impact revealed 46variants present in 90% of samples and two actionable therapeutic targets: IDH1_p.Arg132Cysand KRAS_p.Gly12Val. Also, 80% of the samples presented a chromothripsis signature and 60%presented the SBS3 COSMIC signature related to homologous recombination deficiency (HRD).Analysis of the proliferation score in RNA-Seq data revealed a notable difference among thetumor groups, signifying a higher proliferation index in the uLMS group, suggesting a moreaggressive growth pattern with implications for clinical behavior and disease progression. Afterapplying a customized workflow to identify balanced chromosomal abnormalities of genefusions with clinical impact, the proto-oncogene eukaryotic elongation factor 1 alpha (EEF1A1)emerged as a prominent fusion partner. Eight fusions in 7 samples were identified withtranscriptional overexpression in positive cases. At the proteomic level the enrichment analysisunveiled activated pathways associated with the innate immune system, while suppressedpathways pertained to extracellular matrix (ECM) organization and smooth muscle contraction.With an integrated bioinformatics pipeline, we identified copy number amplification of theCTHRC1 gene in 80% of the samples, accompanied by increased gene expression and proteinlevels. This gene is related to collagen remodeling and degradation and has been demonstratedto be upregulated in some solid tumors and is closely associated with malignancy. Survivalcurve analysis highlighted a significantly poorer prognosis in uLMS exhibiting CTHRC1 geneamplification with potential clinical relevance of CTHRC1 amplification as a prognostic indicatorin uLMS. Taken together these results show the overall suppression of ECM and collagenremodeling in uLMS associated with poor prognosis. Additionally, these data provide the firstevidence that the EEF1A1 gene has a strong selective pressure in uLMS underscoring theindispensable role of precision genomics in advancing precision medicine. Citation Format: Raul Maia Falcao, Jorge Estefano de Souza, William Mathieson, Joseph Woodward Carlson, Tirzah Braz Petta. Deep multi-omics analysis of uterine leiomyosarcoma reveals defective HRD signature, a novel gene fusion and the amplification of the poor-prognosis CTHRC1 gene [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 1751.

Detection of Escherichia coli O157:H7 in Ground Beef Using Long-Read Sequencing.

Foodborne pathogens are a significant cause of illness, and infection with Shiga toxin-producing Escherichia coli (STEC) may lead to life-threatening complications. The current methods to identify STEC in meat involve culture-based, molecular, and proteomic assays and take at least four days to complete. This time could be reduced by using long-read whole-genome sequencing to identify foodborne pathogens. Therefore, the goal of this project was to evaluate the use of long-read sequencing to detect STEC in ground beef. The objectives of the project included establishing optimal sequencing parameters, determining the limit of detection of all STEC virulence genes of interest in pure cultures and spiked ground beef, and evaluating selective sequencing to enhance STEC detection in ground beef. Sequencing libraries were run on the Oxford Nanopore Technologies' MinION sequencer. Optimal sequencing output was obtained using the default parameters in MinKNOW, except for setting the minimum read length to 1 kb. All genes of interest (eae, stx1, stx2, fliC, wzx, wzy, and rrsC) were detected in DNA extracted from STEC pure cultures within 1 h of sequencing, and 30× coverage was obtained within 2 h. All virulence genes were confidently detected in STEC DNA quantities as low as 12.5 ng. In STEC-inoculated ground beef, software-controlled selective sequencing improved virulence gene detection; however, several virulence genes were not detected due to high bovine DNA concentrations in the samples. The growth enrichment of inoculated meat samples in mTSB resulted in a 100-fold increase in virulence gene detection as compared to the unenriched samples. The results of this project suggest that further development of long-read sequencing protocols may result in a faster, less labor-intensive method to detect STEC in ground beef.

In-depth characterization of multidrug-resistant NDM-1 and KPC-3 co-producing Klebsiella pneumoniae bloodstream isolates from Italian hospital patients.

Bloodstream infection (BSI) caused by carbapenem-resistant Klebsiella pneumoniae (KP) poses significant challenges, particularly when the infecting isolate carries multiple antimicrobial resistance (AMR) genes/determinants. This study, employing short- and long-read whole-genome sequencing, characterizes six New Delhi metallo-β-lactamase (NDM) 1 and KP carbapenemase (KPC) 3 co-producing KP isolates, the largest cohort investigated in Europe to date. Five [sequence type (ST) 512] and one (ST11) isolates were recovered from patients who developed BSI from February to August 2022 or February 2023 at two different hospitals in Rome, Italy. Phylogenetic analysis revealed two distinct clusters among ST512 isolates and a separate cluster for the ST11 isolate. Beyond blaNDM-1 and blaKPC-3, various AMR genes, indicative of a multidrug resistance phenotype, including colistin resistance, were found. Each cluster-representative ST512 isolate harbored a blaNDM-1 plasmid (IncC) and a blaKPC-3 plasmid [IncFIB(pQil)/IncFII(K)], while the ST11 isolate harbored a blaNDM-1 plasmid [IncFII(pKPX1)] and a blaKPC-3 plasmid [IncFIB(K)/IncFII(K)]. The blaNDM-1 plasmids carried genes conferring resistance to clinically relevant antimicrobial agents, and the aminoglycoside resistance gene aac(6')-Ib was found on different plasmids. Colistin resistance-associated mgrB/pmrB gene mutations were present in all isolates, and the yersiniabactin-encoding ybt gene was unique to the ST11 isolate. In conclusion, our findings provide insights into the genomic context of blaNDM-1/blaKPC-3 carbapenemase-producing KP isolates.IMPORTANCEThis study underscores the critical role of genomic surveillance as a proactive measure to restrict the spread of carbapenemase-producing KP isolates, especially when key antimicrobial resistance genes, such as blaNDM-1/blaKPC-3, are plasmid borne. In-depth characterization of these isolates may help identify plasmid similarities contributing to their intra-hospital/inter-hospital adaptation and transmission. Despite the lack of data on patient movements, it is possible that carbapenem-resistant isolates were selected to co-produce KP carbapenemase and New Delhi metallo-β-lactamase via plasmid acquisition. Studies employing long-read whole-genome sequencing should be encouraged to address the emergence of KP clones with converging phenotypes of virulence and resistance to last-resort antimicrobial agents.

Plasmid-encoded gene duplications of extended-spectrum β-lactamases in clinical bacterial isolates.

The emergence of extended-spectrum β-lactamase (ESBL)-producing Enterobacteriaceae is an urgent and alarming One Health problem. This study aimed to investigate duplications of plasmid-encoded ESBL genes and their impact on antimicrobial resistance (AMR) phenotypes in clinical and screening isolates. Multi-drug-resistant bacteria from hospitalized patients were collected during routine clinical surveillance from January 2022 to June 2023, and their antimicrobial susceptibility patterns were determined. Genotypes were extracted from long-read whole-genome sequencing data. Furthermore, plasmids and other mobile genetic elements associated with ESBL genes were characterized, and the ESBL genes were correlated to ceftazidime minimal inhibitory concentration (MIC). In total, we identified four cases of plasmid-encoded ESBL gene duplications that match four genetically similar plasmids during the 18-month surveillance period: five Escherichia coli and three Klebsiella pneumoniae isolates. As the ESBL genes were part of transposable elements, the surrounding sequence regions were duplicated as well. In-depth analysis revealed insertion sequence (IS)-mediated transposition mechanisms. Isolates with duplicated ESBL genes exhibited a higher MIC for ceftazidime in comparison to isolates with a single gene copy (3-256 vs. 1.5-32 mg/L, respectively). ESBL gene duplications led to an increased phenotypic resistance against ceftazidime. Our data suggest that ESBL gene duplications by an IS-mediated transposition are a relevant mechanism for how AMR develops in the clinical setting and is part of the microevolution of plasmids.

Long-read sequencing reveals the structural complexity of genomic integration of HPV DNA in cervical cancer cell lines

BackgroundCervical cancer (CC) causes more than 311,000 deaths annually worldwide. The integration of human papillomavirus (HPV) is a crucial genetic event that contributes to cervical carcinogenesis. Despite HPV DNA integration is known to disrupt the genomic architecture of both the host and viral genomes in CC, the complexity of this process remains largely unexplored.ResultsIn this study, we conducted whole-genome sequencing (WGS) at 55-65X coverage utilizing the PacBio long-read sequencing platform in SiHa and HeLa cells, followed by comprehensive analyses of the sequence data to elucidate the complexity of HPV integration. Firstly, our results demonstrated that PacBio long-read sequencing effectively identifies HPV integration breakpoints with comparable accuracy to targeted-capture Next-generation sequencing (NGS) methods. Secondly, we constructed detailed models of complex integrated genome structures that included both the HPV genome and nearby regions of the human genome by utilizing PacBio long-read WGS. Thirdly, our sequencing results revealed the occurrence of a wide variety of genome-wide structural variations (SVs) in SiHa and HeLa cells. Additionally, our analysis further revealed a potential correlation between changes in gene expression levels and SVs on chromosome 13 in the genome of SiHa cells.ConclusionsUsing PacBio long-read sequencing, we have successfully constructed complex models illustrating HPV integrated genome structures in SiHa and HeLa cells. This accomplishment serves as a compelling demonstration of the valuable capabilities of long-read sequencing in detecting and characterizing HPV genomic integration structures within human cells. Furthermore, these findings offer critical insights into the complex process of HPV16 and HPV18 integration and their potential contribution to the development of cervical cancer.

HAT: de novo variant calling for highly accurate short-read and long-read sequencing data.

de novo variants (DNVs) are variants that are present in offspring but not in their parents. DNVs are both important for examining mutation rates as well as in the identification of disease-related variation. While efforts have been made to call DNVs, calling of DNVs is still challenging from parent-child sequenced trio data. We developed Hare And Tortoise (HAT) as an automated DNV detection workflow for highly accurate short-read and long-read sequencing data. Reliable detection of DNVs is important for human genomics and HAT addresses this need. HAT is a computational workflow that begins with aligned read data (i.e. CRAM or BAM) from a parent-child sequenced trio and outputs DNVs. HAT detects high-quality DNVs from Illumina short-read whole-exome sequencing, Illumina short-read whole-genome sequencing, and highly accurate PacBio HiFi long-read whole-genome sequencing data. The quality of these DNVs is high based on a series of quality metrics including number of DNVs per individual, percent of DNVs at CpG sites, and percent of DNVs phased to the paternal chromosome of origin. https://github.com/TNTurnerLab/HAT.

Optical genome mapping identifies a homozygous deletion in the non-coding region of the SCN9A gene in individuals from the same family with congenital insensitivity to pain.

We report an index patient with complete insensitivity to pain and a history of painless fractures, joint hypermobility, and behavioral problems. The index patient descends from a family with notable cases among his maternal relatives, including his aunt and his mother's first cousin, both of whom suffer from congenital insensitivity to pain. The patient had normal results for prior genetic testing: fragile-X syndrome testing, chromosomal microarray analysis, and exome sequencing. Optical genome mapping detected a homozygous deletion affecting the noncoding 5' untranslated region (UTR) and the first non-coding exon of the SCN9A gene in all affected family members, compatible with recessive disease transmission. Pathogenic homozygous loss-of-function variants in the SCN9A gene are associated with impaired pain sensation in humans. Optical genome mapping can thus detect pathogenic structural variants in patients without molecular etiology by standard diagnostic procedures and is a more accessible diagnostic tool than short-read or long-read whole-genome sequencing.

Genotype Analysis of Respiratory Syncytial Virus Before and After the COVID-19 Pandemic Using Whole-Genome Sequencing: A Prospective, Single-Center Study in Korea From 2019 to 2022.

Respiratory syncytial virus (RSV), a highly transmissible virus, is the leading cause of lower respiratory tract infections. We examined molecular changes in the RSV genome before and after the coronavirus disease 2019 (COVID-19) pandemic in Korea, and investigated whether drug-resistant mutations were present. In this prospective, single-center study, RSV-positive respiratory samples were collected between September 2019 and December 2022. Long-read whole-genome sequencing (WGS) was performed, and the presence of known drug-resistant substitutions for palivizumab, nirsevimab, and suptavumab was investigated. Overall, 288 respiratory samples were collected from 276 children. WGS data were available for 133 samples (71 and 62 samples from the pre- and post-pandemic periods, respectively). All RSV-A strains (n = 56) belonged to the GA2.3.5 (ON1) genotype, whereas all RSV-B strains (n = 77) belonged to the GB5.0.5a (BA) genotype. No significant differences in genotypes were observed between the pre- and post-pandemic periods. In addition, no notable mutations related to nirsevimab or palivizumab resistance were detected in the F gene. However, the L172Q and S173L substitutions, which are known to confer resistance to suptavumab, were present in all RSV-B samples. Despite the unprecedented interruption of RSV seasonality, there were no significant molecular changes in circulating RSV strains in Korea related to nirsevimab or palivizumab resistance before and after the COVID-19 pandemic. However, RSV-specific drug-resistance substitutions for suptavumab were identified.

Prolonged carriage of OXA-244-carbapenemase-producing Escherichia coli complicates epidemiological investigations

The rapid increase of OXA-244-producing Escherichia coli, predominantly driven by genetically clustered isolates of sequence type (ST)38, has been observed in at least nine European countries, including Germany. However, the reasons for the spread of OXA-244-producing E. coli remain unclear. Here, we aim to evaluate the possibility of prolonged carriage. We identified a total of six different patients with repeated detection of OXA-244-producing E. coli isolates, which were subjected to both short and long-read whole-genome sequencing (WGS). Besides allelic differences using core genome multilocus sequence typing (cgMLST) analyses, we obtained numbers of single-nucleotide polymorphisms (SNPs) to calculate individual base-pair substitution (BPS) rates. To assess possible re-exposure and risk factors for prolonged carriage, case interviews were conducted. The time between detections ranged from eleven months to more than three years. Initial isolates originated in three+ out of six cases from clinical samples, whereas remaining samples were from screening, mostly in the inpatient setting. As expected, cgMLST analyses showed low numbers of allelic differences between isolates of each case ranging from 1 to 4, whereas numbers of SNPs were between 2 and 99 (mean = 36), thus clearly highlighting the discrepancy between these different bacterial typing approaches. For five out of six cases, observed BPS rates suggest that patients can be colonized with OXA-244-producing E. coli, including ST38 cluster isolates, for extensively long times. Thus, we may have previously missed the epidemiological link between cases because exposure to OXA-244-producing E. coli could have occurred in a time frame, which has not been evaluated in previous investigations. Our results may help to guide future epidemiological investigations as well as to support the interpretation of genetic diversity of OXA-244-producing E. coli, particularly among ST38 cluster isolates.

Advancing animal tuberculosis surveillance using culture-independent long-read whole-genome sequencing.

Animal tuberculosis is a significant infectious disease affecting both livestock and wildlife populations worldwide. Effective disease surveillance and characterization of Mycobacterium bovis (M. bovis) strains are essential for understanding transmission dynamics and implementing control measures. Currently, sequencing of genomic information has relied on culture-based methods, which are time-consuming, resource-demanding, and concerning in terms of biosafety. This study explores the use of culture-independent long-read whole-genome sequencing (WGS) for a better understanding of M. bovis epidemiology in African buffaloes (Syncerus caffer). By comparing two sequencing approaches, we evaluated the efficacy of Illumina WGS performed on culture extracts and culture-independent Oxford Nanopore adaptive sampling (NAS). Our objective was to assess the potential of NAS to detect genomic variants without sample culture. In addition, culture-independent amplicon sequencing, targeting mycobacterial-specific housekeeping and full-length 16S rRNA genes, was applied to investigate the presence of microorganisms, including nontuberculous mycobacteria. The sequencing quality obtained from DNA extracted directly from tissues using NAS is comparable to the sequencing quality of reads generated from culture-derived DNA using both NAS and Illumina technologies. We present a new approach that provides complete and accurate genome sequence reconstruction, culture independently, and using an economically affordable technique.

Genotypic characterization of Streptococcus didelphis causative of fatal infection in white-eared opossums.

Streptococcus didelphis was reported once as related to severe infections in opossums. Thus, we present the first comprehensive whole-genome characterization of clinical S. didelphis strains isolated from white-eared opossums (Didelphis albiventris). Long-read whole-genome sequencing was performed using the MinION platform, which allowed the prediction of several genomic features. We observed that S. didelphis genomes harbor a cluster for streptolysin biosynthesis, and a conserved genomic island with genes involved in transcriptional regulation (arlR) and transmembrane transport (bcrA). Antimicrobial resistance genes for several drug classes were found including beta-lactam, which is the main antimicrobial class used in Streptococcus spp. infections; however no phenotypical resistance was observed. In addition, we predicted the presence of 33 virulence factors in the analyzed genomes. High phylogenetic similarity was observed between clinical and reference strains, yet no clonality was suggested. We also proposed dnaN, gki, pros, and xpt as housekeeping candidates to be used in S. didelphis sequence typing. This is the first whole-genome characterization of S. didelphis, whose data provides important insights on its pathogenicity.

A large pedigree study confirmed the CGG repeat expansion of RILPL1 Is associated with oculopharyngodistal myopathy

BackgroundOculopharyngodistal myopathy (OPDM) is an autosomal dominant adult-onset degenerative muscle disorder characterized by ptosis, ophthalmoplegia and weakness of the facial, pharyngeal and limb muscles. Trinucleotide repeat expansions in non-coding regions of LRP12, G1PC1, NOTCH2NLC and RILPL1 were reported to be the etiologies for OPDM.ResultsIn this study, we performed long-read whole-genome sequencing in a large five-generation family of 156 individuals, including 21 patients diagnosed with typical OPDM. We identified CGG repeat expansions in 5’UTR of RILPL1 gene in all patients we tested while no CGG expansion in unaffected family members. Repeat-primed PCR and fluorescence amplicon length analysis PCR were further confirmed the segregation of CGG expansions in other family members and 1000 normal Chinese controls. Methylation analysis indicated that methylation levels of the RILPL1 gene were unaltered in OPDM patients, which was consistent with previous studies. Our findings provide evidence that RILPL1 is associated OPDM in this large pedigree.ConclusionsOur results identified RILPL1 is the associated the disease in this large pedigree.