Nanopore Sequencing Research Articles

The correlation between CpG methylation and gene expression is driven by sequence variants

Gene promoter and enhancer sequences are bound by transcription factors and are depleted of methylated CpG sites (cytosines preceding guanines in DNA). The absence of methylated CpGs in these sequences typically correlates with increased gene expression, indicating a regulatory role for methylation. We used nanopore sequencing to determine haplotype-specific methylation rates of 15.3 million CpG units in 7,179 whole-blood genomes. We identified 189,178 methylation depleted sequences where three or more proximal CpGs were unmethylated on at least one haplotype. A total of 77,789 methylation depleted sequences (~41%) associated with 80,503 cis-acting sequence variants, which we termed allele-specific methylation quantitative trait loci (ASM-QTLs). RNA sequencing of 896 samples from the same blood draws used to perform nanopore sequencing showed that the ASM-QTL, that is, DNA sequence variability, drives most of the correlation found between gene expression and CpG methylation. ASM-QTLs were enriched 40.2-fold (95% confidence interval 32.2, 49.9) among sequence variants associating with hematological traits, demonstrating that ASM-QTLs are important functional units in the noncoding genome.

Integration of CTCF loops, methylome, and transcriptome in differentiating LUHMES as a model for imprinting dynamics of the 15q11-q13 locus in human neurons.

Human cell line models, including the neuronal precursor line LUHMES, are important for investigating developmental transcriptional dynamics within imprinted regions, particularly the 15q11-q13 Angelman (AS) and Prader-Willi (PWS) syndrome locus. AS results from loss of maternal UBE3A in neurons, where the paternal allele is silenced by a convergent antisense transcript UBE3A-ATS, a lncRNA that terminates at PWAR1 in non-neurons. qRT-PCR analysis confirmed the exclusive and progressive increase in UBE3A-ATS in differentiating LUHMES neurons, validating their use for studying UBE3A silencing. Genome-wide transcriptome analyses revealed changes to 11 834 genes during neuronal differentiation, including the upregulation of most genes within the 15q11-q13 locus. To identify dynamic changes in chromatin loops linked to transcriptional activity, we performed a HiChIP validated by 4C, which identified two neuron-specific CTCF loops between MAGEL2-SNRPN and PWAR1-UBE3A. To determine if allele-specific differentially methylated regions (DMR) may be associated with CTCF loop anchors, whole genome long-read nanopore sequencing was performed. We identified a paternally hypomethylated DMR near the SNRPN upstream loop anchor exclusive to neurons and a paternally hypermethylated DMR near the PWAR1 CTCF anchor exclusive to undifferentiated cells, consistent with increases in neuronal transcription. Additionally, DMRs near CTCF loop anchors were observed in both cell types, indicative of allele-specific differences in chromatin loops regulating imprinted transcription. These results provide an integrated view of the 15q11-q13 epigenetic landscape during LUHMES neuronal differentiation, underscoring the complex interplay of transcription, chromatin looping, and DNA methylation. They also provide insights for future therapeutic approaches for AS and PWS.

A cross-sectional survey of Blastocystis sp. and Dientamoeba fragilis in non-human primates and their caregivers in Czech zoos

Intestinal protists in the gut microbiome are increasingly studied, but their basic epidemiology is not well understood. We explored the prevalence, genetic diversity, and potential zoonotic transmission of two protists colonizing the large intestine - Blastocystis sp. and Dientamoeba fragilis - in 37 species of non-human primates (NHPs) and their caregivers in six zoos in the Czech Republic. We analyzed 179 fecal samples (159 from NHPs, 20 from humans) by qPCR. Blastocystis sp. was detected in 54.7% (98/179) of samples, in 24 NHP species and in 57.2% of NHP samples (prevalence ranged between 36 and 80%), and in 35% of human samples (prevalence ranged between 0 and 67%). Using next generation amplicon sequencing, nine Blastocystis subtypes (ST1-ST5, ST7, ST8, and two novel subtypes) were identified. The two new Blastocystis subtypes (named ST47 and ST48) were described using Nanopore sequencing to produce full-length reference sequences of the small subunit ribosomal RNA gene. Some subtypes were shared between NHPs and their caregivers, suggesting potential zoonotic transmission. Mixed subtype colonization was frequently observed, with 52% of sequenced samples containing two or more subtypes. Dientamoeba was found only in NHPs with a prevalence of 6%. This study emphasizes the critical role of molecular diagnostics in epidemiological and transmission studies of these protists and calls for further research to better understand their impact on public health.

Faecal microbial diversity in a cattle herd infected by Mycobacterium avium subsp. paratuberculosis: a possible effect of production status

Mycobacterium avium subsp. paratuberculosis (MAP) causes Johne’s disease, or paratuberculosis (PTB) in ruminants, besides having zoonotic potential. It possibly changes the gut microbiome, but no conclusive data are available yet. This study aimed at investigating the influence of MAP on the faecal microbiome of cattle naturally infected with PTB. In a follow up period of 10 months, PTB status was investigated in a herd of dairy cattle with history of clinical cases. Each animal was tested for MAP infection using serum and milk ELISA for MAP anti-bodies and IS900 real-time PCR and recombinase polymerase amplification assays for MAP DNA in the faeces and milk monthly for 4 successive months, then a last one after 6 months. The faecal samples were subjected to 16S rDNA metagenomic analysis using Oxford Nanopore Sequencing Technology. The microbial content was compared between animal groups based on MAP positivity rate and production status. All animals were MAP positive by one or more tests, but two animals were consistently negative for MAP DNA in the faeces. In all animals, the phyla firmicutes and bacteroidetes were highly enriched with a small contribution of proteobacteria, and increased abundance of the families Oscillospiraceae, Planococcaceae, and Streptococcacaceae was noted. Animals with high MAP positivity rate showed comparable faecal microbial content, although MAP faecal positivity had no significant effect (p > 0.05) on the microbiome. Generally, richness and evenness indices decreased with increasing positivity rate. A significantly different microbial content was found between dry cows and heifers (p < 0.05). Particularly, Oscillospiraceae and Rikenellaceae were enriched in heifers, while Planococcaceae and Streptococcaceae were overrepresented in dry cows. Furthermore, abundance of 72 genera was significantly different between these two groups (p < 0.05). Changes in faecal microbiome composition were notably associated with increasing MAP shedding in the faeces. The present findings suggest a combined influence of the production status and MAP on the cattle faecal microbiome. This possibly correlates with the fate of the infection, the concern in disease control, again remains for further investigations.

Transverse aortic constriction multi-omics analysis uncovers pathophysiological cardiac molecular mechanisms.

Time-course multi-omics data of a murine model of progressive heart failure (HF) induced by transverse aortic constriction (TAC) provide insights into the molecular mechanisms that are causatively involved in contractile failure and structural cardiac remodelling. We employ Illumina-based transcriptomics, Nanopore sequencing and mass spectrometry-based proteomics on samples from the left ventricle (LV) and right ventricle (RV, RNA only) of the heart at 1, 7, 21 and 56 days following TAC and Sham surgery. Here, we present Transverse Aortic COnstriction Multi-omics Analysis (TACOMA), as an interactive web application that integrates and visualizes transcriptomics and proteomics data collected in a TAC time-course experiment. TACOMA enables users to visualize the expression profile of known and novel genes and protein products thereof. Importantly, we capture alternative splicing events by assessing differential transcript and exon usage as well. Co-expression-based clustering algorithms and functional enrichment analysis revealed overrepresented annotations of biological processes and molecular functions at the protein and gene levels. To enhance data integration, TACOMA synchronizes transcriptomics and proteomics profiles, enabling cross-omics comparisons. With TACOMA (https://shiny.dieterichlab.org/app/tacoma), we offer a rich web-based resource to uncover molecular events and biological processes implicated in contractile failure and cardiac hypertrophy. For example, we highlight: (i) changes in metabolic genes and proteins in the time course of hypertrophic growth and contractile impairment; (ii) identification of RNA splicing changes in the expression of Tpm2 isoforms between RV and LV; and (iii) novel transcripts and genes likely contributing to the pathogenesis of HF. We plan to extend these data with additional environmental and genetic models of HF to decipher common and distinct molecular changes in heart diseases of different aetiologies. Database URL: https://shiny.dieterichlab.org/app/tacoma.

Applications of Nanopore sequencing in precision cancer medicine.

Oxford Nanopore Technologies sequencing, also referred to as Nanopore sequencing, stands at the forefront of a revolution in clinical genetics, offering the potential for rapid, long read, and real-time DNA and RNA sequencing. This technology is currently making sequencing more accessible and affordable. In this comprehensive review, we explore its potential regarding precision cancer diagnostics and treatment. We encompass a critical analysis of clinical cases where Nanopore sequencing was successfully applied to identify point mutations, splice variants, gene fusions, epigenetic modifications, non-coding RNAs, and other pivotal biomarkers that defined subsequent treatment strategies. Additionally, we address the challenges of clinical applications of Nanopore sequencing and discuss the current efforts to overcome them.

Targeted nanopore sequencing using clinical specimens for the rapid diagnosis of extrapulmonary tuberculosis

BackgroundThe clinical presentation of extrapulmonary tuberculosis (EPTB) is atypical and it is easily confused with other diseases such as common infections, making prompt diagnosis a great challenge. This study aimed to evaluate the accuracy of targeted nanopore sequencing (TNS) in the diagnosis of EPTB. The diagnostic accuracy of TNS using different types of extrapulmonary specimens was also evaluated.MethodsWe reviewed the clinical data of patients with suspected EPTB for whom TNS was conducted and who were hospitalized at our center. The true positive, false positive, false negative, and true negative values were determined. Indices of diagnostic accuracy were computed, including sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV), and area under the curve (AUC) for TNS and acid-fast bacilli (AFB) culture, and compared with those from clinical diagnosis.Results149 patients were included in the analysis. The overall sensitivity, specificity, PPV, NPV, and AUC of TNS for the diagnosis of EPTB were 86.4%, 87.5%, 97.3%, 55.3%, and 0.87, respectively. For diagnosis by AFB culture, these values were 25.6%, 100.0%, 100.0%, 20.5%, and 0.63, respectively. The most common specimens used were lymph node tissue, cerebrospinal fluid, pleural effusion, and pleural tissue. The diagnostic accuracy of TNS using all types of extrapulmonary specimens was good.ConclusionsTNS demonstrates good diagnostic accuracy in the rapid diagnosis of EPTB and this was true across different types of extrapulmonary specimens.

Chromosome-level genome assembly of the two-spotted spider mite Tetranychus urticae

The two-spotted spider mite, Tetranychus urticae Koch (Acari: Tetranychidae), is a notorious pest in agriculture that has developed resistance to almost all chemical types used for its control. Here, we assembled a chromosome-level genome for the TSSM using Illumina, Nanopore, and Hi-C sequencing technologies. The assembled contigs had a total length of 103.94 Mb with an N50 of 3.46 Mb, with 87.7 Mb of 34 contigs anchored to three chromosomes. The chromosome-level genome assembly had a BUSCO completeness of 94.8%. We identified 15,604 protein-coding genes, with 11,435 genes that could be functionally annotated. The high-quality genome provides invaluable resources for the genetic and evolutionary study of TSSM.

Complete genome sequence of multidrug-resistant Raoultella terrigena strain RT01-5M1 isolated from farmed salmon in Canada.

The genome sequence of multidrug-resistant Raoultella terrigena RT01-5M1 strain isolated from Canadian farmed salmon was determined using Oxford nanopore and Illumina MiSeq sequencers. The assembled chromosome was estimated at 5,699,993 bp in size, with two plasmids, 164,879 bp and 82,046 bp. The chromosome and smaller plasmid contained antimicrobial resistance genes.

High-quality reference genome of cowpea beetle Callosobruchus maculatus

Callosobruchus maculatus is one of the most competitive stored grain pests, which causes a great loss to agricultural economy. However, due to an inadequacy of high-quality reference genome, the molecular mechanisms for olfactory and hypoxic adaptations to stored environments are unknown and require to be revealed urgently, which will contribute to the detection and prevention of the invasive pests C. maculatus. Here, we presented a high-quality chromosome-level genome of C. maculatus based on Illumina, Nanopore and Hi-C sequencing data. The total size was 1.2 Gb, and 65.17% (797.47 Mb) of it was identified to be repeat sequences. Among assembled chromosomes, chromosome 10 was considered the X chromosome according to the evidence of reads coverage and homologous genes among species. The current version of high-quality genome provides preferable data resources for the adaptive evolution research of C. maculatus.

The Jersey Daffodil Project: Integrating nanopore sequencing into classrooms improves STEM skills, scientific identity and career development

Societal Impact StatementThe Jersey Daffodil Project, a secondary school initiative, boosted students' aspirations in science careers. Fieldwork, DNA sequencing and discussions with scientists enriched their learning and aspirations. Surveys revealed a positive shift in students' views towards STEM careers. The project also influenced an entire cohort to pursue biological sciences in 2021. The use of cutting‐edge technology reduced costs and time, showing promise for broader school collaboration. This project not only transforms science education but also inspires students, shaping their futures and fostering a new generation of scientists.

Complete genome sequencing of Agrobacterium leguminum strain EL101 isolated from Khaya senegalensis tree bark using nanopore technology.

We report the complete genome sequence of Agrobacterium leguminum strain EL101, isolated from the tree bark of Khaya senegalensis in Kampar, Perak, Malaysia, obtained using Q20+ Nanopore Sequencing chemistry. The assembled genome has a total length of 5,324,685 bp, comprising a circular chromosome, a linear chromid, and two non-Ti circular plasmids.

Chromosome-Scale Genome Assembly of the Sheep-Biting Louse Bovicola ovis Using Nanopore Sequencing Data and Pore-C Analysis.

Bovicola ovis, commonly known as the sheep-biting louse, is an ectoparasite that adversely affects the sheep industry. Sheep louse infestation lowers the quality of products, including wool and leather, causing a loss of approximately AUD 123M per annum in Australia alone. The lack of a high-quality genome assembly for the sheep-biting louse, as well as any closely related livestock lice, has hindered the development of louse research and management control tools. In this study, we present the assembly of B. ovis with a genome size of ~123 Mbp based on a nanopore long-read sequencing library and Illumina RNA sequencing, complemented with a chromosome-level scaffolding using the Pore-C multiway chromatin contact dataset. Combining multiple alignment and gene prediction tools, a comprehensive annotation on the assembled B. ovis genome was conducted and recalled 11,810 genes as well as other genomic features including orf, ssr, rRNA and tRNA. A manual curation using alignment with the available closely related louse species, Pediculus humanus, increased the number of annotated genes to 16,024. Overall, this study reported critical genetic resources and biological insights for the advancement of sheep louse research and the development of sustainable control strategies in the sheep industry.

Exploring Microorganisms Associated to Acute Febrile Illness and Severe Neurological Disorders of Unknown Origin: A Nanopore Metagenomics Approach.

Acute febrile illness (AFI) and severe neurological disorders (SNDs) often present diagnostic challenges due to their potential origins from a wide range of infectious agents. Nanopore metagenomics is emerging as a powerful tool for identifying the microorganisms potentially responsible for these undiagnosed clinical cases. In this study, we aim to shed light on the etiological agents underlying AFI and SND cases that conventional diagnostic methods have not been able to fully elucidate. Our approach involved analyzing samples from fourteen hospitalized patients using a comprehensive nanopore metagenomic approach. This process included RNA extraction and enrichment using the SMART-9N protocol, followed by nanopore sequencing. Subsequent steps involved quality control, host DNA/cDNA removal, de novo genome assembly, and taxonomic classification. Our findings in AFI cases revealed a spectrum of disease-associated microbes, including Escherichia coli, Streptococcus sp., Human Immunodeficiency Virus 1 (Subtype B), and Human Pegivirus. Similarly, SND cases revealed the presence of pathogens such as Escherichia coli, Clostridium sp., and Dengue virus type 2 (Genotype-II lineage). This study employed a metagenomic analysis method, demonstrating its efficiency and adaptability in pathogen identification. Our investigation successfully identified pathogens likely associated with AFI and SNDs, underscoring the feasibility of retrieving near-complete genomes from RNA viruses. These findings offer promising prospects for advancing our understanding and control of infectious diseases, by facilitating detailed genomic analysis which is critical for developing targeted interventions and therapeutic strategies.

LINE1-mediated epigenetic repression of androgen receptor transcription causes androgen insensitivity syndrome

Androgen insensitivity syndrome (AIS) is a difference of sex development (DSD) characterized by different degrees of undervirilization in individuals with a 46,XY karyotype despite normal to high gonadal testosterone production. Classically, AIS is explained by hemizygous mutations in the X-chromosomal androgen receptor (AR) gene. Nevertheless, the majority of individuals with clinically diagnosed AIS do not carry an AR gene mutation. Here, we present a patient with a 46,XY karyotype, born with undervirilized genitalia, age-appropriate testosterone levels and no uterus, characteristic for AIS. Diagnostic whole exome sequencing (WES) showed a maternally inherited LINE1 (L1) retrotransposon insertion in the 5′ untranslated region (5′UTR) of the AR gene. Long-read nanopore sequencing confirmed this as an insertion of a truncated L1 element of ≈ 2.7 kb and showed an increased DNA methylation at the L1 insertion site in patient-derived genital skin fibroblasts (GSFs) compared to healthy controls. The insertion coincided with reduced AR transcript and protein levels in patient-derived GSFs confirming the clinical diagnosis AIS. Our results underline the relevance of retrotransposons in human disease, and expand the growing list of human diseases associated with them.

Transmission of global clones of NDM-producing Enterobacterales and interspecies spread of IncX3 plasmid harbouring blaNDM-5 in Tokyo

ObjectiveThe aim of this study is to characterise the molecular characteristics of NDM-producing Enterobacterales, which have been on the increase in recent years in Japan, where IMP-producing bacteria are dominant among carbapenemase-producing Enterobacterales. MethodsWe collected 21 strains of NDM-producing Enterobacterales detected between 2015 and 2022 at five hospitals in Tokyo and performed illumina whole genome sequencing. For the seven selected strains, nanopore long-read sequencing was also performed to characterise the plasmids harbouring blaNDM. ResultsFourteen strains were Escherichia coli and all carried blaNDM-5. Among these strains, eight and three were sequence type (ST) 410 and ST167, respectively, and both groups of strains were spread clonally in different hospitals. Two strains of Klebsiella pneumoniae ST147 carrying blaNDM-1 were detected in a hospital, and these strains had also spread clonally. The remainder included Enterobacter hormaechei, Klebsiella quasipneumoniae, Citrobacter amalonaticus, and Klebsiella michiganensis. Plasmid analysis revealed that an identical IncX3 plasmid harbouring blaNDM-5 was shared among four strains of different bacterial species (E. coli, C. amalonaticus, K. michiganensis, and E. hormaechei) detected at the same hospital. In addition, a Klebsiella quasipneumoniae strain detected at a different hospital also carried an IncX3 plasmid with a similar genetic structure. ConclusionsNosocomial spread of multiple multidrug-resistant global clones and transmission of IncX3 plasmids harbouring blaNDM-5 among multiple species were detected as the major pathways of spread of NDM-producing Enterobacterales in Tokyo. Early detection of carriers and measures to prevent nosocomial spread are important to prevent further spread of NDM-producing organisms.

Sequencing a CC239-MRSA-III with a novel composite SCC mec element from Kuwait.

Staphylococcus aureus CC239-MRSA-III is an ancient pandemic strain of hospital-associated, methicillin-resistant S. aureus that spread globally for decades and that still can be found in some parts of the world. In Kuwait, microarray-based surveillance identified from 2019 to 2022 a series of isolates of a hitherto unknown variant of this strain that carried a second set of recombinase genes, ccrA/B-2. To elucidate the structure of its SCCmec element, two isolates were subjected to nanopore sequencing. This revealed, in addition to ccrA/B-2, several SCC-associated genes including speG (spermidine N acetyltransferase) and a gene encoding a large "E-domain containing protein" (dubbed as edcP-SCC). This gene contained three regions consisting of multiple repeating units. In terms of sequence and structure it was similar but not identical to the biofilm-related aap gene from S. epidermidis. A review of published sequences identified edcP-SCC in eighteen genome sequences of S. aureus, S. epidermidis and S. capitis, and frequently it appears in a similar cluster of genes as in the strains sequenced herein. Isolates also carried a prophage with the adhesion factor sasX/sesI and aminoglycoside resistance genes. This is consistent with an affiliation to the "South-East Asian" Clade of CC239. The emergence of edcP-SCC and sasX-positive CC239 strain shows that, against a global trend towards community-associated MRSA, the ancient pandemic CC239 hospital strain still continues to evolve and to cause outbreaks.

Decoding the Epigenetic Landscape: Insights into 5mC and 5hmC Patterns in Mouse Cortical Cell Types.

The DNA modifications, 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC), represent powerful epigenetic regulators of temporal and spatial gene expression. Yet, how the cooperation of these genome-wide, epigenetic marks determine unique transcriptional signatures across different brain cell populations is unclear. Here we applied Nanopore sequencing of native DNA to obtain a complete, genome-wide, single-base resolution atlas of 5mC and 5hmC modifications in neurons, astrocytes and microglia in the mouse cortex (99% genome coverage, 40 million CpG sites). In tandem with RNA sequencing, analysis of 5mC and 5hmC patterns across cell types reveals astrocytes drive uniquely high brain 5hmC levels and support two decades of research regarding methylation patterns, gene expression and alternative splicing, benchmarking this resource. As such, we provide the most comprehensive DNA methylation data in mouse brain as an interactive, online tool (NAM-Me, https://olsenlab.shinyapps.io/NAMME/) to serve as a resource dataset for those interested in the methylome landscape.

Simultaneous de novo calling and phasing of genetic variants at chromosome-scale using NanoStrand-seq

The successful accomplishment of the first telomere-to-telomere human genome assembly, T2T-CHM13, marked a milestone in achieving completeness of the human reference genome. The upcoming era of genome study will focus on fully phased diploid genome assembly, with an emphasis on genetic differences between individual haplotypes. Most existing sequencing approaches only achieved localized haplotype phasing and relied on additional pedigree information for further whole-chromosome scale phasing. The short-read-based Strand-seq method is able to directly phase single nucleotide polymorphisms (SNPs) at whole-chromosome scale but falls short when it comes to phasing structural variations (SVs). To shed light on this issue, we developed a Nanopore sequencing platform-based Strand-seq approach, which we named NanoStrand-seq. This method allowed for de novo SNP calling with high precision (99.52%) and acheived a superior phasing accuracy (0.02% Hamming error rate) at whole-chromosome scale, a level of performance comparable to Strand-seq for haplotype phasing of the GM12878 genome. Importantly, we demonstrated that NanoStrand-seq can efficiently resolve the MHC locus, a highly polymorphic genomic region. Moreover, NanoStrand-seq enabled independent direct calling and phasing of deletions and insertions at whole-chromosome level; when applied to long genomic regions of SNP homozygosity, it outperformed the strategy that combined Strand-seq with bulk long-read sequencing. Finally, we showed that, like Strand-seq, NanoStrand-seq was also applicable to primary cultured cells. Together, here we provided a novel methodology that enabled interrogation of a full spectrum of haplotype-resolved SNPs and SVs at whole-chromosome scale, with broad applications for species with diploid or even potentially polypoid genomes.

Complete mitochondrial genome of Melia azedarach L., reveals two conformations generated by the repeat sequence mediated recombination

Melia azedarach is a species of enormous value of pharmaceutical industries. Although the chloroplast genome of M. azedarach has been explored, the information of mitochondrial genome (Mt genome) remains surprisingly limited. In this study, we used a hybrid assembly strategy of BGI short-reads and Nanopore long-reads to assemble the Mt genome of M. azedarach. The Mt genome of M. azedarach is characterized by two circular chromosomes with 350,142 bp and 290,387 bp in length, respectively, which encodes 35 protein-coding genes (PCGs), 23 tRNA genes, and 3 rRNA genes. A pair of direct repeats (R1 and R2) were associated with genome recombination, resulting in two conformations based on the Sanger sequencing and Oxford Nanopore sequencing. Comparative analysis identified 19 homologous fragments between Mt and chloroplast genome, with the longest fragment of 12,142 bp. The phylogenetic analysis based on PCGs were consist with the latest classification of the Angiosperm Phylogeny Group. Notably, a total of 356 potential RNA editing sites were predicted based on 35 PCGs, and the editing events lead to the formation of the stop codon in the rps10 gene and the start codons in the nad4L and atp9 genes, which were verified by PCR amplification and Sanger sequencing. Taken together, the exploration of M. azedarach gap-free Mt genome provides a new insight into the evolution research and complex mitogenome architecture.