Nanopore Long-read Sequencing Research Articles

Validation of Nanopore long-read sequencing to resolve RPGR ORF15 genotypes in individuals with X-linked retinitis pigmentosa.

X-linked retinitis pigmentosa (XLRP) is characterized by progressive vision loss leading to legal blindness in males and a broad severity spectrum in carrier females. Pathogenic alterations of the retinitis pigmentosa GTPase regulator gene (RPGR) are responsible for over 70% of XLRP cases. In the retina, the RPGRORF15 transcript includes a terminal exon, called ORF15, that is altered in the large majority of RPGR-XLRP cases. Unfortunately, due to its highly repetitive sequence, ORF15 represents a considerable challenge in terms of sequencing for molecular diagnostic laboratories. However, in a recent preliminary work Yahya et al. reported a long-read sequencing approach seeming promising. Here, the aim of the study was to validate and integrate this new sequencing strategy in a routine screening workflow. For that purpose, we performed a masked test on 52 genomic DNA samples from male and female individuals carrying 32 different pathogenic ORF15 variations including 20 located in the highly repetitive region of the exon. For the latter, we have obtained a detection rate of 80-85% in males and 60-80% in females after bioinformatic analyses. These numbers raised to 100% for both status after adding a complementary visual inspection of ORF15 long-reads. In accordance with these results, and considering the frequency of ORF15 pathogenic variations in XLRP, we suggest that a long-read screening of ORF15 should be systematically considered before any other sequencing approach in subjects with a diagnosis compatible with XLRP.

P-564 Decoding ovarian response: the predictive power of FSH receptor gene expression in controlled ovarian stimulation

Abstract Study question Can profiling of variation in follicle-stimulating hormone receptor (FSHR) gene expression forecast response to Controlled Ovarian Stimulation (COS) in Assisted Reproductive Techniques? Summary answer Full-length FSHR transcript was the most common mRNA species. Nine distinct transcripts were characterised, but only one (deletion Exon 6) was associated with altered response. What is known already AMH levels and age are considered useful predictors of ovarian response to COS (high AMH and young age typically associated with high yields of oocytes; older ages and low AMH associated with fewer oocytes). However, it is not uncommon to see patients who do not fit this pattern, highlighting the limitations of these parameters in predicting response. Age-independent biomarkers capable of revealing patients at elevated risk of hyperstimulation, or those who might produce fewer than the optimal number of oocytes, are desirable. Alternative transcripts of the FSHR gene have been proposed as an explanation for variation in response to COS. Study design, size, duration FSHR splice variants in mRNA in granulosa cells from 126 patients were examined utilising long-read nanopore sequencing. This strategy permitted the sequencing of mRNA variants without fragmentation, allowing the detection of all splicing events in each mRNA strand. In addition to being categorised based on age and AMH levels, patients were also classified by the expression of FSHR transcript variants and their response to stimulation (high, low, or normal response). Participants/materials, setting, methods Patterns of FSHR mRNA splice variants were analysed in each of 126 patients. This involved extraction of RNA from the granulosa cells, followed by long-read nanopore sequencing, permitting analysis of whole transcripts in a single, contiguous sequence ‘read’. Patients were classified based upon their level of response to COS (low, normal, high), taking into account their age, AMH level and the number of oocytes retrieved. Patterns of FSHR transcript variation were compared to COS response. Main results and the role of chance Full-length FSHR transcript was detected in 97.6% of patients and skipping of Exon 6 was observed in 68%. Five novel transcripts were discovered, including partial insertion of ∼100bp segments of introns 1 or 2, seen in 6.55% of samples. In total, nine distinct FSHR mRNA variants were detected. As expected, increasing age was associated with higher FSH requirement per oocyte collected (p = 0.018), and AMH levels declined with age (2.8± 1.42 in women <38 versus 1.58 ±1.29 in women ≥38; p = 0.011). Consistent with previous analyses, AMH levels and amount of FSH per oocyte retrieved were negatively correlated (p < 0.001). One-third of the patients (33%) in this cohort displayed responses to COS that were not considered to be consistent with expectations based upon age and AMH levels. Expression of the novel FSHR transcript variants was not associated with FSH needed per oocyte retrieved (p = 0.78), nor the proportion of mature oocytes (p = 0.95).The only transcript potentially linked to an altered response to COS was the previously reported mRNA species, lacking exon 6, associated with significantly greater numbers of oocytes produced (p = 0.0087). Interestingly, women ≥38 expressed novel variants more frequently than their younger counterparts (p = 0.018), potentially indicative of an age-related decline in transcriptional fidelity. Limitations, reasons for caution It is possible that the use of controlled ovarian stimulation might alter the dynamics of transcription in cumulus cells. Therefore, it is possible that the spectrum of transcripts seen in this study might differ from those present in unstimulated cycles. Wider implications of the findings Expected associations between age, AMH levels and response to COS were observed. Skipping exon 6 may be associated with a superior response to COS, suggesting FSHR mRNA profiling might have a role in tailoring controlled ovarian stimulation protocols. Novel transcript variants were detected and seen most often in older women. Trial registration number not applicable

Direct nanopore RNA sequencing of umbra-like virus-infected plants reveals long non-coding RNAs, specific cleavage sites, D-RNAs, foldback RNAs, and temporal- and tissue-specific profiles.

The traditional view of plus (+)-strand RNA virus transcriptomes is that infected cells contain a limited variety of viral RNAs, such as full-length (+)-strand genomic RNA(s), (-)-strand replication intermediate(s), 3' co-terminal subgenomic RNA(s), and viral recombinant defective (D)-RNAs. To ascertain the full complement of viral RNAs associated with the simplest plant viruses, long-read direct RNA nanopore sequencing was used to perform transcriptomic analyses of two related umbra-like viruses: citrus yellow vein-associated virus (CY1) from citrus and CY2 from hemp. Analysis of different timepoints/tissues in CY1- and CY2-infected Nicotiana benthamiana plants and CY2-infected hemp revealed: (i) three 5' co-terminal RNAs of 281 nt, 442 nt and 671 nt, each generated by a different mechanism; (ii) D-RNA populations containing the 671 fragment at their 5'ends; (iii) many full-length genomic RNAs and D-RNAs with identical 3'end 61 nt truncations; (iv) virtually all (-)-strand reads missing 3 nt at their 3' termini; (v) (±) foldback RNAs comprising about one-third of all (-)-strand readsand (vi) a higher proportion of full-length gRNAs in roots than in leaves, suggesting that roots may be functioning as a gRNA reservoir. These findings suggest that viral transcriptomes are much more complex than previously thought.

Engineering psychrophilic polymerase for nanopore long-read sequencing.

Unveiling the potential application of psychrophilic polymerases as candidates for polymerase-nanopore long-read sequencing presents a departure from conventional choices such as thermophilic Bacillus stearothermophilus (Bst) renowned for its limitation in temperature and mesophilic Bacillus subtilis phage (phi29) polymerases for limitations in strong exonuclease activity and weak salt tolerance. Exploiting the PB-Bst fusion DNA polymerases from Psychrobacillus (PB) and Bacillus stearothermophilus (Bst), our structural and biochemical analysis reveal a remarkable enhancement in salt tolerance and a concurrent reduction in exonuclease activity, achieved through targeted substitution of a pivotal functional domain. The sulfolobus 7-kDa protein (Sso7d) emerges as a standout fusion domain, imparting significant improvements in PB-Bst processivity. Notably, this study elucidates additional functional sites regulating exonuclease activity (Asp43 and Glu45) and processivity using artificial nucleotides (Glu266, Gln283, Leu334, Glu335, Ser426, and Asp430). By disclosing the intricate dynamics in exonuclease activity, strand displacement, and artificial nucleotide-based processivity at specific functional sites, our findings not only advance the fundamental understanding of psychrophilic polymerases but also provide novel insights into polymerase engineering.

Detection of four isomers of the human cytomegalovirus genome using nanopore long-read sequencing.

Human cytomegalovirus has a linear DNA genome with a total length of approximately 235kb. This large genome is divided into two domains, "Long" and "Short". There are four isomers of the cytomegalovirus genome with different orientations of each domain. To confirm the presence of four types of isomers, it is necessary to identify the sequence of the junction between the domains. However, due to the presence of repeat sequences, it is difficult to determine the junction sequences by next-generation sequencing analysis. To solve this problem, long-read sequencing was performed using the Oxford Nanopore sequencer and the junctions were successfully identified in four isomers in strain Merin and ATCC-2011-3. Nanopore sequencing also revealed the presence of multiple copies of the "a" sequence (a-seq) in the junctions, indicating the diversity of the junction sequences. These results strongly suggest that long-read sequencing using the nanopore sequencer would be beneficial for identifying the complex structure of the cytomegalovirus genome.

Identification of transcriptionally-active human papillomavirus integrants through nanopore sequencing reveals viable targets for gene therapy against cervical cancer.

Integration of the human papillomavirus (HPV) genome into the cellular genome is a key event that leads to constitutive expression of viral oncoprotein E6/E7 and drives the progression of cervical cancer. However, HPV integration patterns differ on a case-by-case basis among related malignancies. Next-generation sequencing technologies still face challenges for interrogating HPV integration sites. In this study, utilizing Nanopore long-read sequencing, we identified 452 and 108 potential integration sites from the cervical cancer cell lines (CaSki and HeLa) and five tissue samples, respectively. Based on long Nanopore chimeric reads, we were able to analyze the methylation status of the HPV long control region (LCR), which controls oncogene E6/E7 expression, and to identify transcriptionally-active integrants among the numerous integrants. As a proof of concept, we identified an active HPV integrant in between RUNX2 and CLIC5 on chromosome 6 in the CaSki cell line, which was supported by ATAC-seq, H3K27Ac ChIP-seq, and RNA-seq analysis. Knockout of the active HPV integrant, by the CRISPR/Cas9 system, dramatically crippled cell proliferation and induced cell senescence. In conclusion, identifying transcriptionally-active HPV integrants with Nanopore sequencing can provide viable targets for gene therapy against HPV-associated cancers.

Accurate parent-of-origin variant assignment for multiple hereditary cancer syndromes using proband-only blood sample analysis.

10516 Background: Predicting which side of the family a germline variant comes from is a critical gap in current clinical practice, vital for risk management, variant curation, and cascade genetic testing. Assignment of any autosomal variant to either parent with 99% accuracy is now possible with only a blood sample from the proband. Parent-of-Origin-Aware genomic analysis (POAga) is achieved by combining methylation and sequence data from Oxford Nanopore long-read sequencing with chromosome-length haplotypes generated from Strand-seq to infer parent of origin of any variant along the length of a chromosome, due to accurate phasing of imprinted differentially methylated regions that occur on each autosome. We sought to validate POAga in common, high penetrant hereditary cancer conditions such as hereditary breast and ovarian cancer (HBOC) and Lynch syndrome, rarer syndromes with parent-of-origin-effects and other genes predisposing to breast cancer and gastrointestinal malignancies that are associated with genes across multiple chromosomes. Methods: Blood samples from carriers of pathogenic variants in ATM, BRCA1, BRCA2, CDH1, MLH1, MSH2, MSH6, PMS2, EPCAM, PALB2, SDHD and SDHAF2, with known parental segregation, are currently being ascertained to determine the analytic validity of POAga in real world samples of differing age, sex, ethnicity, and cancer status. Samples are undergoing whole genome analysis by long and short read sequencing. Parent-of-origin of the pathogenic or likely pathogenic variant or variant of uncertain significance is predicted according to previously described methods (Akbari V, Hanlon VCT, et al. Cell Genom. 2022 Dec 21;3(1):100233.) under an REB approved protocol. Results: To date, analysis is complete for 100 individuals with a total of 107 rare or pathogenic germline variants with known or presumed parental segregation. Germline variants are in SDHD (n = 18), BRCA2 (n = 17), BRCA1 (n = 14), MLH1 (n = 10), PALB2 (n = 9), PMS2 (n = 8), MSH2 (n = 9), MSH6 (n = 7), ATM (n = 5), CDH1 (n = 5), SDHAF2 (n = 1), DICER1 (n = 1), MUTYH (n = 1), RET (n = 1), and EPCAM (n = 1). Of variants able to be assigned a parent-of-origin (n = 104 of 107, 97%), there was complete concordance between the predicted parent-of-origin and known clinical segregation (n = 104, 100%). Conclusions: Results to date support the ability of POAga to accurately infer parent-of-origin of rare or pathogenic variants with known parental segregation using only a blood sample from carriers of diverse hereditary cancer syndromes. Ongoing validation of POAga will continue to test its feasibility in real-world samples and inform its path towards clinically translation. Parent-of-Origin-Aware genomic analysis is a powerful technology that could improve our understanding of hereditary cancer syndromes and transform our ability to conduct genetic cancer risk assessments for patients and families.

Phenomics and genomic features of Enterococcus avium IRMC1622a isolated from a clinical sample of hospitalized patient

BackgroundEnterococcus avium (E. avium) is a Gram-positive nosocomial pathogen that is commonly isolated from the alimentary tract. The objective of this functional genomics study was to identify the resistant genes by analyzing the genome of E. avium IRMC1622a, a type of bacteria found in feces collected from a patient at a Saudi Arabian tertiary hospital. MethodsThe bacterial strain IRMC1622a was identified by 16 S rRNA sequencing as Enterococcus sp. The resistance phenomics were performed using VITEK® 2, and morphological analysis was achieved using a scanning electron microscope (SEM). Finally, the whole bacterial genome of the bacterial strain IRMC1622a was subjected to sequencing during October 2023 using Oxford Nanopore long-read sequencing technology, and mining for resistant genes. ResultsThe results of antimicrobial resistant phenomics indicated that the IRMC1622a strain was sensitive to all tested antimicrobial agents except for erythromycin, and the same result was confirmed by genomic analysis in addition to other classes of antibiotics. SEM showed E. avium IRMC1622a is ovoid shape, in single cells (L 1.2797 ± 0.1490 µm), in pairs (L 1.7333 ± 0.1054 µm), and in chains (L 2.44033 ± 0.1978 µm). The E. avium IRMC1622a genome has 14 (in CARD) antimicrobial resistance genes that were identified with several mechanisms of antimicrobial resistance, such as the efflux pump and conferring antibiotic resistance. The present study revealed that the E. avium IRMC1622a genome contains a high number of genes associated with virulence factors, and 14 matched pathogenic protein families and predicted as human pathogen (probability score 0.855). We report two (ISEnfa4 and ISEfa5) mobile genetic elements for the first time in the E. avium genome. ConclusionsThe study concludes that E. avium IRMC1622a is susceptible to all tested antibacterials except erythromycin. The IRMC1622a has 14 genes encoding antimicrobial resistance mechanisms, including the efflux pump and conferring antibiotic resistance. This could indicate a potential rise in E. avium resistance in healthcare facilities. These observations may raise concerns regarding E. avium resistance in healthcare. We need more research to understand the pathophysiology of E. avium, which leads to hospital-acquired infections.

Unravelling epigenetic mechanisms in Cerastoderma edule genome: a comparison of healthy and neoplastic cockles

Cancer is a multifaceted genetic disease characterized by the acquisition of several essential hallmarks. Notably, certain cancers exhibit horizontal transmissibility, observed across mammalian species and diverse bivalves, the latter referred to as hemic neoplasia. Within this complex landscape, epigenetic mechanisms such as histone modifications and cytosine methylation emerge as fundamental contributors to the pathogenesis of these transmissible cancers. Our study delves into the epigenetic landscape of Cerastoderma edule, focusing on whole-genome methylation and hydroxymethylation profiles in heathy specimens and transmissible neoplasias by means of Nanopore long-read sequencing. Our results unveiled a global hypomethylation in the neoplastic specimens compared to their healthy counterparts, emphasizing the role of DNA methylation in these tumorigenic processes. Furthermore, we verified that intragenic CpG methylation positively correlated with gene expression, emphasizing its role in modulating transcription in healthy and neoplastic cockles, as also highlighted by some up-methylated oncogenic genes. Hydroxymethylation levels were significantly more elevated in the neoplastic samples, particularly within satellites and complex repeats, likely related to structural functions. Additionally, our analysis also revealed distinct methylation and activity patterns in retrotransposons, providing additional insights into bivalve neoplastic processes. Altogether, these findings contribute to understanding the epigenetic dynamics of bivalve neoplasias and shed light on the roles of DNA methylation and hydroxymethylation in tumorigenesis. Understanding these epigenetic alterations holds promise for advancing our broader understanding of cancer epigenetics.

Accurate and comprehensive evaluation of O6-methylguanine-DNA methyltransferase promoter methylation by nanopore sequencing.

The methylation status of the O6-methylguanine-DNA methyltransferase (MGMT) promoter region is essential in evaluating the prognosis and predicting the drug response in patients with glioblastoma. In this study, we evaluated the utility of using nanopore long-read sequencing as a method for assessing methylation levels throughout the MGMT CpG-island, compared its performance to established techniques and demonstrated its clinical applicability. We analysed 165 samples from CNS tumours, focusing on the MGMT CpG-island using nanopore sequencing. Oxford Nanopore Technologies (ONT) MinION and PromethION flow cells were employed for single sample or barcoded assays, guided by a CRISPR/Cas9 protocol, adaptive sampling or as part of a whole genome sequencing assay. Methylation data obtained through nanopore sequencing were compared to results obtained via pyrosequencing and methylation bead arrays. Hierarchical clustering was applied to nanopore sequencing data for patient stratification. Nanopore sequencing displayed a strong correlation (R2 = 0.91) with pyrosequencing results for the four CpGs of MGMT analysed by both methods. The MGMT-STP27 algorithm's classification was effectively reproduced using nanopore data. Unsupervised hierarchical clustering revealed distinct patterns in methylated and unmethylated samples, providing comparable survival prediction capabilities. Nanopore sequencing yielded high-confidence results in a rapid timeframe, typically within hours of sequencing, and extended the analysis to all 98 CpGs of the MGMT CpG-island. This study presents nanopore sequencing as a valid and efficient method for determining MGMT promotor methylation status. It offers a comprehensive view of the MGMT promoter methylation landscape, which enables the identification of potentially clinically relevant subgroups of patients. Further exploration of the clinical implications of patient stratification using nanopore sequencing of MGMT is warranted.

LINE-1 retrotransposons contribute to mouse PV interneuron development

Retrotransposons are mobile DNA sequences duplicated via transcription and reverse transcription of an RNA intermediate. Cis-regulatory elements encoded by retrotransposons can also promote the transcription of adjacent genes. Somatic LINE-1 (L1) retrotransposon insertions have been detected in mammalian neurons. It is, however, unclear whether L1 sequences are mobile in only some neuronal lineages or therein promote neurodevelopmental gene expression. Here we report programmed L1 activation by SOX6, a transcription factor critical for parvalbumin (PV) interneuron development. Mouse PV interneurons permit L1 mobilization in vitro and in vivo, harbor unmethylated L1 promoters and express full-length L1 mRNAs and proteins. Using nanopore long-read sequencing, we identify unmethylated L1s proximal to PV interneuron genes, including a novel L1 promoter-driven Caps2 transcript isoform that enhances neuron morphological complexity in vitro. These data highlight the contribution made by L1 cis-regulatory elements to PV interneuron development and transcriptome diversity, uncovered due to L1 mobility in this milieu.

P29 Long-read sequencing provides unprecedented insight into the skin and wound microbial landscape

Abstract Introduction and aims It is becoming increasingly evident that the microbial inhabitants of our skin play a key role in cutaneous health and disease. However, our understanding of the contribution of the complex skin microbiome to wound healing outcome remains incomplete. This is, in part, due to the historic challenges around traditional sampling and profiling Methods, which often fail to capture the abundance and diversity of the skin microbiome. Recent advances in microbiome sequencing technologies now enable more comprehensive characterization of skin bacterial composition and host-microbe interactions. Thus, the aim of this study was to utilize a novel long-read sequencing platform to provide species-level microbiome characterization and deliver important insight into the drivers of pathological skin repair. Methods Microbial DNA samples were collected from the skin and wounds of healthy volunteers and patients with chronic wounds (n = 40 per group). Long-read nanopore sequencing was used to profile microbiome samples, while RNA-Seq enabled generation of donor transcriptomic profiles. A multifaceted approach was used to integrate these datasets to investigate the complex interplay between host skin pathology and the microbiome. Results We first identified perturbations in the skin microbiome of patients with chronic wounds, which was significantly less diverse than healthy skin and colonized with key bacteria associated with pathogenicity. These alterations in the skin microbiome also correlated with transcriptomic changes, such as upregulation of genes associated with inflammation, stress and poor healing. Interestingly, long-read sequencing data displayed poor correlation with routine clinical microbiology swabs, while specific clinical characteristics (e.g. glycaemic control and wound duration) were strongly linked to wound microbiome composition. Conclusions Collectively, these data provide crucial new insight into the skin and wound microbial landscape, demonstrating a potential role for specific wound bacteria in modulating pathological skin repair.

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.

Simultaneous profiling of chromatin accessibility and DNA methylation in complete plant genomes using long-read sequencing.

Epigenetic regulations, including chromatin accessibility, nucleosome positioningand DNA methylation intricately shape genome function. However, current chromatin profiling techniques relying on short-read sequencing technologies fail to characterise highly repetitive genomic regions and cannot detect multiple chromatin features simultaneously. Here, we performed Simultaneous Accessibility and DNA Methylation Sequencing (SAM-seq) of purified plant nuclei. Thanks to the use of long-read nanopore sequencing, SAM-seq enables high-resolution profiling of m6A-tagged chromatin accessibility together with endogenous cytosine methylation in plants. Analysis of naked genomic DNA revealed significant sequence preference biases of m6A-MTases, controllable through a normalisation step. By applying SAM-seq to Arabidopsis and maize nuclei we obtained fine-grained accessibility and DNA methylation landscapes genome-wide. We uncovered crosstalk between chromatin accessibility and DNA methylation within nucleosomes of genes, TEs, and centromeric repeats. SAM-seq also detects DNA footprints over cis-regulatory regions. Furthermore, using the single-molecule information provided by SAM-seq we identified extensive cellular heterogeneity at chromatin domains with antagonistic chromatin marks, suggesting that bivalency reflects cell-specific regulations. SAM-seq is a powerful approach to simultaneously study multiple epigenetic features over unique and repetitive sequences, opening new opportunities for the investigation of epigenetic mechanisms.

Assembly of continuous high-resolution draft genome sequence of Hemicentrotus pulcherrimus using long-read sequencing.

The update of the draft genome assembly of sea urchin, Hemicentrotus pulcherrimus, which is widely studied in East Asia as a model organism of early development, was performed using Oxford nanopore long-read sequencing. The updated assembly provided ~600-Mb genome sequences divided into 2,163 contigs with N50 = 516 kb. BUSCO completeness score and transcriptome model mapping ratio (TMMR) of the present assembly were obtained as 96.5% and 77.8%, respectively. These results were more continuous with higher resolution than those by the previous version of H. pulcherrimus draft genome, HpulGenome_v1, where the number of scaffolds = 16,251 with a total of ~100 Mb, N50 = 143 kb, BUSCO completeness score = 86.1%, and TMMR = 55.4%. The obtained genome contained 36,055 gene models that were consistent with those in other echinoderms. Additionally, two tandem repeat sequences of early histone gene locus containing 47 copies and 34 copies of all histone genes, and 185 of the homologous sequences of the interspecifically conserved region of the Ars insulator, ArsInsC, were obtained. These results provide further advance for genome-wide research of development, gene regulation, and intranuclear structural dynamics of multicellular organisms using H. pulcherrimus.

Klebsiella pneumoniae sequence type 147: a high-risk clone increasingly associated with plasmids carrying both resistance and virulence elements.

Introduction. The first hybrid resistance/virulence plasmid, combining elements from virulence plasmids described in hypervirulent types of Klebsiella pneumoniae with those from conjugative resistance plasmids, was described in an isolate of sequence type (ST) 147 from 2016. Subsequently, this type has been increasingly associated with these plasmids.Hypothesis or gap statement. The extent of carriage of hybrid virulence/resistance plasmids in nosocomial isolates of K. pneumoniae requires further investigation.Aim. To describe the occurrence of virulence/resistance plasmids among isolates of K. pneumoniae received by the UK reference laboratory, particularly among representatives of ST147, and to compare their sequences.Methodology. Isolates received by the laboratory during 2022 and the first half of 2023 (n=1278) were screened for virulence plasmids by PCR detection of rmpA/rmpA2 and typed by variable-number tandem repeat analysis. Twenty-nine representatives of ST147 (including a single-locus variant) from seven hospital laboratories were subjected to long-read nanopore sequencing using high-accuracy q20 chemistry to provide complete assemblies.Results. rmpA/rmpA2 were detected in 110 isolates, of which 59 belonged to hypervirulent K1-ST23, K2-ST86 and K2-ST65/375. Of the remainder, representatives of ST147 formed the largest group, with 22 rmpA/rmpA2-positive representatives (out of 47 isolates). Representatives were from 19 hospital laboratories, with rmpA/rmpA2-positive isolates from 10. Nanopore sequencing of 29 representatives of ST147 divided them into those with no virulence plasmid (n=12), those with non-New Delhi metallo-β-lactamase (NDM) virulence plasmids (n=6) and those carrying bla NDM-5 (n=9) or bla NDM-1 (n=2) virulence plasmids. These plasmids were of IncFIB(pNDM-Mar)/IncHI1B(pNDM-MAR) replicon types. Most of the non-NDM virulence plasmids were highly similar to the originally described KpvST147L_NDM plasmid. Those carrying bla NDM-5 were highly similar to one another and to previously described plasmids in ST383 and carried an extensive array of resistance genes. Comparison of the fully assembled chromosomes indicated multiple introductions of ST147 in UK hospitals.Conclusion. This study highlights the high proportion of representatives of ST147 that carry IncFIB(pNDM-Mar)/IncHI1B(pNDM-MAR) hybrid resistance virulence plasmids. It is important to be aware of the high probability that representatives of this type carry these plasmids combining resistance and virulence determinants and of the consequent increased risk to patients.

High-resolution and real-time wastewater viral surveillance by Nanopore sequencing

Wastewater genomic sequencing stands as a pivotal complementary tool for viral surveillance in populations. While long-read Nanopore sequencing is a promising platform to provide real-time genomic data, concerns over the sequencing accuracy of the earlier Nanopore versions have somewhat restrained its widespread application in wastewater analysis. Here, we evaluate the latest improved version of Nanopore sequencing (R10.4.1), using SARS-CoV-2 as the model infectious virus, to demonstrate its effectiveness in wastewater viral monitoring. By comparing amplicon lengths of 400 bp and 1200 bp, we revealed that shorter PCR amplification is more suitable for wastewater samples due to viral genome fragmentation. Utilizing mock wastewater samples, we validated the reliability of Nanopore sequencing for variant identification by comparing it with Illumina sequencing results. The strength of Nanopore sequencing in generating real-time genomic data for providing early warning signals was also showcased, indicating that as little as 0.001 Gb of data can provide accurate results for variant prevalence. Our evaluation also identified optimal alteration frequency cutoffs (>50 %) for precise mutation profiling, achieving >99 % precision in detecting single nucleotide variants (SNVs) and insertions/deletions (indels). Monitoring two major wastewater treatment plants in Hong Kong from September 2022 to April 2023, covering over 4.5 million population, we observed a transition in dominant variants from BA.5 to XBB lineages, with XBB.1.5 being the most prevalent variants. Mutation detection also highlighted the potential of wastewater Nanopore sequencing in uncovering novel mutations and revealed links between signature mutations and specific variants. This study not only reveals the environmental implications of Nanopore sequencing in SARS-CoV-2 surveillance but also underscores its potential in broader applications including environmental health monitoring of other epidemic viruses, which could significantly enhance the field of wastewater-based epidemiology.

Combined approaches, including long-read sequencing, address the diagnostic challenge of HYDIN in primary ciliary dyskinesia

Primary ciliary dyskinesia (PCD), a disorder of the motile cilia, is now recognised as an underdiagnosed cause of bronchiectasis. Accurate PCD diagnosis comprises clinical assessment, analysis of cilia and the identification of biallelic variants in one of 50 known PCD-related genes, including HYDIN. HYDIN-related PCD is underdiagnosed due to the presence of a pseudogene, HYDIN2, with 98% sequence homology to HYDIN. This presents a significant challenge for Short-Read Next Generation Sequencing (SR-NGS) and analysis, and many diagnostic PCD gene panels do not include HYDIN. We have used a combined approach of SR-NGS with bioinformatic masking of HYDIN2, and state-of-the-art long-read Nanopore sequencing (LR_NGS), together with analysis of respiratory cilia including transmission electron microscopy and immunofluorescence to address the underdiagnosis of HYDIN as a cause of PCD. Bioinformatic masking of HYDIN2 after SR-NGS facilitated the detection of biallelic HYDIN variants in 15 of 437 families, but compromised the detection of copy number variants. Supplementing testing with LR-NGS detected HYDIN deletions in 2 families, where SR-NGS had detected a single heterozygous HYDIN variant. LR-NGS was also able to confirm true homozygosity in 2 families when parental testing was not possible. Utilising a combined genomic diagnostic approach, biallelic HYDIN variants were detected in 17 families from 242 genetically confirmed PCD cases, comprising 7% of our PCD cohort. This represents the largest reported HYDIN cohort to date and highlights previous underdiagnosis of HYDIN-associated PCD. Moreover this provides further evidence for the utility of LR-NGS in diagnostic testing, particularly for regions of high genomic complexity.

Whole-genome assembly of a hybrid Trypanosoma cruzi strain assembled with Nanopore sequencing alone.

Trypanosoma cruzi is the causative agent of Chagas disease, which causes 10,000 deaths per year. Despite the high mortality associated with Chagas, relatively few parasite genomes have been assembled to date, with genome assemblies unavailable even for some commonly used laboratory strains. This is at least partially due to T. cruzi's highly complex and highly repetitive genome, which defies investigation using traditional short-read sequencing methods. In this study, we have generated a high-quality whole-genome assembly of the hybrid Tulahuen strain, a commercially available type VI strain, using long-read Nanopore sequencing without short-read scaffolding. The assembled genome contains 25% repeat regions, 17% variable multigene family members, and 27% transposable elements (TEs) and is of comparable quality with T. cruzi genome assemblies that utilized both long- and short-read data. Notably, we find that regions with TEs are significantly enriched for multicopy surface proteins, and that surface proteins are, on average, closer to TEs than to other coding regions. This finding suggests that mobile genetic elements such as transposons may drive recombination within surface protein gene families. This work demonstrates the feasibility of Nanopore sequencing to resolve complex regions of T. cruzi genomes, and with these resolved regions, provides support for a possible mechanism for genomic diversification.

Abstract LB185: High resolution telomere measurements in human cancer and aging using long-read nanopore sequencing

Abstract Telomere length is an important biomarker of cellular replicative potential and aging, but existing measurement methods are limited in resolution and accuracy. During human carcinogenesis, telomeres shorten profoundly in the earliest stages that precede tumor invasion. Telomere shortening can either halt cancer progression by inducing senescence and cell death or can promote cancer development by destabilizing chromosomes and generating chromosomal rearrangements. Stabilization of telomeres, typically by the upregulation of telomerase, is required for progression of nearly all human cancers. Here, we deploy digital telomere measurement by nanopore sequencing to understand how distributions of human telomere length change during aging and following carcinogenesis. We measure telomere attrition and de novo elongation with unprecedented resolution in genetically defined populations of human stem cells, in cells from the peripheral blood of healthy individuals and patients with genetic defects in telomere maintenance, cancer patients, and cancer cell lines. We find that human aging is accompanied by a progressive loss of long telomeres and an accumulation of shorter telomeres in white blood cells. In patients with inherited defects in telomere maintenance, the accumulation of short telomeres is more pronounced and correlates with phenotypic severity. We apply machine learning to train a binary classification model that distinguishes healthy individuals from those with telomere biology disorders. Next, we investigated the telomere lengths of a cohort of twenty colorectal tumors with patient-matched healthy colonic epithelia and found that two-thirds of colorectal tumors have telomere length distributions significantly shorter than the surrounding normal tissue. Finally, we explore the utility of telomere length distributions in cancer cell lines as clinically actionable disease biomarkers. This sequencing and bioinformatic pipeline will advance our understanding of telomere maintenance mechanisms and the use of telomere length as a clinical biomarker of aging and disease. Citation Format: Santiago Sanchez, Jessica Gu, Anudeep Golla, William Shomali, Annika Martin, Dirk Hockemeyer, Sharon Savage, Steven Artandi. High resolution telomere measurements in human cancer and aging using long-read nanopore sequencing [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 2 (Late-Breaking, Clinical Trial, and Invited Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(7_Suppl):Abstract nr LB185.