Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “nanopore sequencing”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8Linked to original sources

Comprehensive benchmarking of somatic structural variant detection at ultra-low allele fractions.

Postzygotic mosaicism gives rise to somatic structural variants (SVs) at ultra-low variant allele fractions (VAFs), which pose challenges for detection due to the high-coverage sequencing required and noise introduced by sequencing artifacts. Although somatic SV detection has been extensively studied in cancer, these studies are not directly applicable to the study of tissue mosaicism, as they rely on matched normals, target higher VAF ranges, and are enriched for different types of SVs. We present comprehensive benchmark data and best practices for non-cancer somatic SV detection. We created a synthetic mosaic sample by combining six HapMap individuals at varying proportions, generating allele fractions as low as 0.25%. This sample was sequenced to ~2,300x total coverage using Illumina, PacBio, and Nanopore technologies across multiple sequencing centers. A high-confidence benchmark SV set containing over 21,000 pseudo-somatic insertions and deletions ≥50bp was derived from haplotype-resolved assemblies. We evaluated 12 SV discovery pipelines and identified caller-specific strengths and sequencing platform-specific shortcomings. We find that short read-based approaches show reduced recall for insertions and repeat-associated SVs, whereas long-read sequencing achieves high accuracy throughout the genome, increasing linearly with coverage. The best algorithm's sensitivity exceeded 80% for VAFs ≥4% and 15% for VAFs of 0.5-1% with 60x coverage. The publicly available benchmarking data and comparative analysis of current methods provide a foundation for robust discovery of SV mosaicism in non-cancer tissues..

Journal Article↗

Degradation of ribosomal RNA during Plasmodium falciparum gametocytogenesis.

The life cycle of Plasmodium falciparum is characterized by complex regulatory changes that allow adaptation of the parasites to different environmental conditions, which are especially pronounced during transmission between the mammalian host and the insect vector. Previous studies have shown that P. falciparum uses three types of ribosomal RNAs (rRNA A, S1, and S2 types) at different stages of its life cycle. We used Oxford Nanopore Technologies direct RNA sequencing to investigate the dynamics of rRNA usage throughout the parasite's intraerythrocytic development, as well as in salivary gland sporozoites. Our study revealed a preponderance of A-type rRNAs during the intraerythrocytic cycle and gametocytogenesis, while S-type rRNAs slowly increase in abundance in mosquito stages starting 3 days post-infection. Salivary gland sporozoites showed an even proportion of all rRNA types. By examining the length distributions of rRNA molecules, we detected an extensive and specific degradation of rRNAs during gametocytogenesis, starting in stage II gametocytes and continuing until the final stages of gametocyte development. We hypothesize that rRNA degradation may be linked to the global translational repression and metabolic quiescence described in stage V gametocytes, similar to mechanisms observed in bacterial and eukaryotic stress responses.IMPORTANCEOur study uses Oxford Nanopore direct RNA sequencing of tightly synchronized blood-stage Plasmodium falciparum parasites to investigate the expression of ribosomal RNAs during asexual and sexual development. P. falciparum utilizes distinct types of rRNA during its development. However, due to the challenges of differentiating these highly similar molecules, their regulation and the mechanism underlying the switch between rRNA types remain unclear. We observe significant rRNA degradation in mature gametocytes, leading us to propose that this potentially leads to a reduced number of functional ribosomes when parasites become quiescent and translation is repressed.

Plasmodium falciparum↗

Exploring biosynthetic potential of the endophytic Penicillium turbatum BLH34 using whole-genome sequence analysis and molecular networking.

An in-depth genomic and metabolomic investigation was conducted on the endophytic fungus Penicillium turbatum BLH34, isolated from Macleaya cordata. Hybrid sequencing (Illumina-Nanopore) generated a high-quality 27.9 Mb genome (GC 48.6%) encoding 9798 proteins, with functional annotation linking 5350 genes to the NCBI non-redundant database and 3404 to KEGG pathways. AntiSMASH analysis uncovered 35 biosynthetic gene clusters (BGCs), 23 of which lacked homology to known pathways, highlighting BLH34's potential for novel metabolite discovery. Molecular networking (GNPS) and LC-MS/MS identified 19 specialised metabolites, including antimicrobial polyketides. Bioassays demonstrated potent inhibition against Staphylococcus aureus (36 mm), Bacillus subtilis (28 mm) and Escherichia coli (24 mm), underscoring its pharmaceutical relevance.

Penicillium↗

Radiation-resistant and desiccation-tolerant bacteria from the Chavara-Neendakara high background radiation area, india: phenotypic characterisation and genomic insights.

Radiation-resistant microorganisms that survive high doses of ionising radiation serve as valuable models for understanding stress adaptation; however, the genomic determinants underlying extreme radiation tolerance in bacteria from natural environments with high background radiation remain insufficiently characterised. Bacterial isolates from the Chavara-Neendakara HBRA (Kerala, India) were evaluated for desiccation tolerance, and the desiccation-resistant isolates were subsequently exposed to gamma irradiation (1-10 kGy) using a 60Co source. Isolates were identified through 16S rRNA sequencing, morphologically characterised by FE-SEM, and screened for antibiotic susceptibility. The highly radiation-resistant strain underwent whole-genome sequencing via Oxford Nanopore Technology, with De novo assembly, polishing, and genome annotation. Four bacterial isolates (Micrococcaceae and Paenibacillaceae) exhibited D10 values of 1-7 kGy, including one multidrug-resistant strain; no endospores were observed in the Paenibacillus isolate under the tested conditions. Paenibacillus sp. HBRA004 survived 10 kGy gamma radiation, exceeding all previously reported HBRA isolates by over fourfold. Its 5.0 Mbp genome (GC = 48.27%, ≥ 99% completeness) encodes five mechanistically independent DNA repair pathways; homologous recombination (recA, recN, radA), base excision repair (mutM, mutY, mutT), mismatch repair (mutL, mutS), nucleotide excision repair (uvrA, uvrB, uvrD), and non-homologous end joining (ku, ligD), alongside a redundant antioxidant network comprising triple-copy Fe/Mn-family superoxide dismutases and ahpC peroxiredoxin. A thioredoxin system (trxA, trxB, msrA) and manganese uptake via mntH may contribute to further layers of ROS defence. Their specific contribution to the HBRA004 phenotype remains to be experimentally and comparatively validated. These findings represent the first genomically characterised 10 kGy-resistant bacterial isolate from the Chavara-Neendakara HBRA, establishing a new benchmark for radiation tolerance within this ecologically significant environment. Pathway depth, gene copy amplification, and Mn/Fe homeostasis appear to be candidate mechanisms contributing to high-level radiation tolerance, consistent with patterns in other radiation-resistant taxa, though their contribution requires functional validation.

India↗

Genomic signatures of host-range divergence in the generalist Beauveria bassiana and the specialist Beauveria brongniartii.

Entomopathogenic fungi of the genus Beauveria are widely used biological control agents that infect diverse insect hosts and can also associate with plants as rhizosphere colonizers and endophytes. Within this genus, Beauveria bassiana is a cosmopolitan generalist, whereas Beauveria brongniartii exhibits a narrower host range, primarily targeting soil-dwelling coleopteran larvae with limited evidence of plant colonization. To explore genomic differentiation associated with this ecological divergence, the commercially exploited B. brongniartii strain BIPESCO2 and B. bassiana ATHUM 4946 were sequenced using Oxford Nanopore technology, followed by comparative genomic analyses across multiple strains. Orthology identified species-specific gene families, although overall genome architecture and core gene content were highly conserved. The CAZyme repertoires were nearly identical, indicating retention of a versatile enzymatic toolkit supporting plant association, saprotrophy, and insect pathogenicity. In contrast, biosynthetic gene clusters displayed substantial variation, including structural remodeling of Beauveria-specific virulence-associated clusters and expansion of type I polyketide synthase clusters in B. brongniartii. Effector prediction revealed a conserved core of largely uncharacterized proteins alongside species-specific orthogroups enriched in adhesion-, immunity-, and cuticle-interaction domains. Together, these findings indicate that host-range divergence in Beauveria is associated with compartmentalized genomic differentiation, particularly in secondary metabolism and a limited subset of lineage-specific virulence factors, rather than in the conserved core infection machinery.

Beauveria↗

Novel bacterial hosts and mobile genetic structure of tet(X) variants in tetracycline-contaminated aquatic environment uncovered by culture and long-read metagenomics.

Clinically important tigecycline (3rd-generation tetracycline) resistance tet(X) variants were inferred to have evolutionarily originated from environmental bacteria, and have been recognized among environment, human and animals. However, genetic basis for environmental proliferation and dissemination of tet(X) variants remains ambiguous. This study profiled tet(X) variants at gene, contig, isolate, and community levels in environmental community subjected to long-term stepwise increasing oxytetracycline (1st-generation tetracycline) or tigecycline pressure using long-term microcosm experiments, quantitative PCR, bacterial isolation, whole-genome sequencing, and Nanopore-based long-read metagenomics. We confirmed that both oxytetracycline and tigecycline enriched the abundance of tetracycline resistance genes especially oxytetracycline-enriched tet(X3). Unexpectedly diverse bacterial hosts and genetic structure of tet(X)-positive mobile elements in the environment microbiome were identified using bacterial isolation and long-read Nanopore metagenomics. Pseudomonas defluvii was first reported to carry tet(X3) in the chromosome, forming IS26-tet(X3)-res-ISCR2 circular intermediate to transfer between different DNA molecules. Database mining revealed similar mobile segments have prevailed among animal-derived Acinetobacter species. Unlike the widely reported ISCR2-mediated transfer of tet(X6), we identified a novel mobile multidrug transposon TnAs3 where tet(X6) and class 1 integron co-transferred as its passenger region. Mobile tet(X2)-ere(D)-aadS-erm(F)-blaOXA-347 segment was annotated in Runella, and co-occurrences of tet(X2) and ere(D), aadS, blaOXA-347 were also found in Flavobacterium, Arsenicibacter, Chryseobacterium and Pedobacter. Overall, tetracycline-contaminated aquatic microbiome harboured diverse mobile tet(X)-positive segments which have not yet been acquired by clinical pathogens, and thus served as the genetic pool of tet(X) variants together with indigenous bacterial hosts, especially the newly reported Pseudomonas defluvii. Reducing pollution of older-generation tetracyclines would be a proactive way to mitigate environmental evolution and possible clinical effects of tet(X) variants.

Metagenomics↗

The planktonic microbiome of the Great Barrier Reef.

Large genome databases have markedly improved our understanding of marine microorganisms1-5. Although these resources have focused on prokaryotes, genomes from many dominant marine lineages, such as Pelagibacter and Prochlorococcus, are conspicuously underrepresented. Here we present the Great Barrier Reef Microbial Genomes Database (GBR-MGD), comprising 5,283 prokaryotic genomes obtained from Great Barrier Reef seawater samples using Nanopore and Illumina sequencing, including a collection of high-quality genomes of underrepresented groups. We show that standard short-read assemblies miss these populations owing to a combination of strain heterogeneity and low-GC-percentage sequencing bias. The GBR-MGD also comprises 20 chromosome-level picoeukaryote and 808,585 viral genomes, including a newly described clade of marine Crassvirales. We demonstrate the utility of the GBR-MGD to identify indicator taxa that can reliably predict the effects of reef management practices, such as the establishment of marine protected zones.

Bacteria↗

Leveraging ONT move table values for signal aware variant calling.

Oxford Nanopore Technologies (ONT) sequencing enables long-range haplotype phasing and contiguous genome assembly but still exhibits elevated error rates that challenge small variant calling, particularly for insertions and deletions (Indels). While raw electrical signals contain rich information, existing signal-aware methods require computationally intensive processing of large signal files. Here, we present Clair3 v2, a method that leverages the ONT move table-a lightweight byproduct of basecalling that maps signal events to nucleotide positions-to improve variant calling accuracy. Clair3 v2 builds upon Clair3 and integrates signal-level dwelling time to significantly enhance variant calling performance. We also propose a genome position based circular buffer to incorporate dwelling time with minimal computational overhead. Benchmarking across six Genome in a Bottle samples demonstrates substantial improvements in variant calling accuracy. With HAC basecalling, Clair3 v2 achieves a mean SNP F1-score of 97.69% at 10 × depth (compared to 96.45% for baseline Clair3), and Indel F1 scores improved from 64.27% to 76.70%, while gains persisted at higher depths. The benefits were most pronounced for longer Indels and in complex genomic regions, where Indel F1 scores in long homopolymer regions improved from 14.3% to 45.2%. Benchmark results across various basecalling modes, samples, and coverage settings outperformed Clair3 baselines and other methods, including DeepVariant and Dorado Variant, and demonstrate the significant benefits of Clair3 v2. Furthermore, Clair3 v2 incurs negligible runtime compared to standard Clair3, making it practical for routine use.

Sequence Analysis, DNA↗

Molecular landscape of methicillin-resistant Staphylococcus aureus strains in clinical infections from hospitals in Lagos, Nigeria.

BACKGROUND AND OBJECTIVES: Multidrug-resistant Staphylococcus aureus (MRSA) accounts for a significant proportion of antimicrobial resistance (AMR)-associated infections worldwide. This study investigated the molecular profile of MRSA in Nigeria, providing valuable genomic data to fill existing knowledge gaps and highlighting its importance in the context of the global AMR crisis. METHODS: A total of 107 isolates were obtained from patient samples, including wound swabs/pus (65 isolates, 60.7%), blood cultures (16 isolates, 15%), urine/urinary catheter (8 isolates, 7.5%) and other sources. Species identification was performed using MALDI-TOF, and antimicrobial susceptibility testing was performed using the VITEK®2 system. Genomic DNA was extracted and subjected to whole-genome sequencing using short-read Illumina technology. In addition, a subset of isolates underwent long-read sequencing using Oxford Nanopore technology. RESULTS: Among the 107 isolates, 63 (59%) were identified as MRSA, with 58 (92%) carrying the mecA gene. The MRSA isolates exhibited high resistance to non-β-lactam antibiotics, particularly trimethoprim/sulfamethoxazole (95.3%), erythromycin (76.6%), gentamicin (71.4%) and quinolones (69.8%). The most prevalent MRSA belonged to the Bengal Bay clone [t657/ST772/Staphylococcal Cassette Chromosome mec (SCCmec) V(5C2)/Panton-Valentine leukocidin (PVL) + methicillin-susceptible Staphylococcus aureus (MRSA)], followed by t4690/ST152/SCCmec Vc(5C2&5)/PVL + MRSA and ST8 (t008, n = 1; t064, n = 4)/SCCmec Vc(5C2&5). Phylogenetic analysis suggests both community/associated transmission and possible importation of strains. CONCLUSIONS: This study highlights the significant burden of MRSA in Nigeria, with the high-risk Bengal Bay MRSA clone as the most common strain. The widespread resistance to non-β-lactam antibiotics underscores the urgent need for enhanced surveillance, infection control and antibiotic stewardship to mitigate its spread.

Journal Article↗

The UTRs of Leishmania donovani vary in length and are enriched in potential regulatory structures.

Leishmania spp. regulate gene expression largely post-transcriptionally, yet untranslated regions (UTRs) remain poorly delineated. We generated high-quality genome and transcriptome datasets for Leishmania donovani strain 1S2D (Ld1S) by combining PacBio HiFi de novo assembly with Oxford Nanopore direct RNA sequencing of promastigotes and axenic amastigotes. The genome assembly consists of 65 scaffolds totaling ~33.3 Mb. Structural comparisons to LdBPK282A1 revealed numerous rearrangements, including some reshuffling genes among polycistronic transcription units and validated by polycistronic reads from RNA sequencing. Promastigote and amastigote RNA sequencing produced 469,010 and 46,729 monocistronic reads containing a spliced-leader and a polyA tail sequences, defining 8,479 transcripts and supporting 7,415 of the 7,969 annotated protein coding genes, as well as 604 putative long non-coding RNAs. We annotated UTRs for 4,921 genes and observed that putative RNA G-quadruplexes were markedly enriched in UTRs. We also noted that 31.9% and 11.5% were expressed into multiple isoforms in promastigotes and amastigotes, respectively. Collectively, these data provide a comprehensive annotation of L. donovani genes and their UTRs and reveal widespread and stage-specific UTR length polymorphisms, and, overall, points to an important role of 3' UTR in post-transcriptional regulation in L. donovani.

Journal Article↗

Long-read transcriptomics corrects Trichomonas vaginalis intron annotations and refines transcript-end features.

BACKGROUND: Trichomonas vaginalis causes the most prevalent non-viral sexually transmitted infection worldwide. Despite its large genome (181.5 Mb; 36,310 predicted protein-coding genes in NYU_TvagG3_2), intron annotations remain limited and inconsistently validated. A recent short-read RNA-seq study reported 63 putative active introns, but short reads can misassign splice boundaries and cannot resolve complete transcript structures. METHODS: We integrated Oxford Nanopore direct RNA sequencing (DRS), ONT cDNA long-read sequencing, and Illumina RNA-seq to refine intron annotations, transcript-end features, and UTR boundaries in T. vaginalis. Candidate introns were validated by targeted PCR and Sanger sequencing, and representative splicing events were further assessed using public SRA datasets. RESULTS: Starting from 31 historically annotated introns, motif-guided long-read screening and orthogonal validation identified 17 additional validated introns, increasing the curated set to 48 confirmed introns. Among these 17 events, three were previously unrecognized in the current NYU_TvagG3_2 reference annotation. We also corrected five reported loci, including two false-positive introns, two splice-coordinate misannotations, and one gene-sequence error. DRS further supported transcript termination site mapping, UAAA polyadenylation-signal profiling relative to poly(A) addition sites, and single-molecule poly(A)-tail estimation. StringTie mixed-mode assemblies provided updated UTR boundaries for intron-bearing transcripts and transcripts without curated introns. CONCLUSIONS: This study provides a rigorously validated, long-read-refined resource of intron annotations, UTR boundaries, and UAAA-guided transcript-end features for T. vaginalis, together with a reproducible workflow for non-model protists. These refinements improve the current reference annotation and support future studies of functional genomics, parasite biology, pathogenesis, and diagnostic development.

Trichomonas vaginalis↗

A MAGIBU-based model for pediatric and juvenile CNS tumors: an in-house epigenetic decision-support framework compared with online DNA methylation classifiers.

Background: DNA methylation profiling is a tool that provides key support for central nervous system (CNS) tumor classification. However, diagnostically ambiguous pediatric cases may result in discordant outputs across classifiers. We developed MAGIBU, a cross-platform, projection-based framework that embeds individual methylomes into a fixed CNS reference landscape, ranking diagnostic entities by local epigenetic proximity to support clinician-led integrative diagnosis. Methods: As a proof-of-concept, we evaluated MAGIBU in eight morphologically challenging pediatric/juvenile CNS tumors with unresolved diagnoses after institutional and central pathology review. To establish a benchmark in the absence of a definitive histopathological ground truth, a consensus epigenetic reference was defined a priori for cases showing concordant results between the Heidelberg CNS Tumor Methylation Classifier and Methylscape Analysis. Comparisons were also performed with Epigenomic Digital Pathology (EpiDiP). To validate MAGIBU beyond this discovery cohort, performance was assessed at the family level across the CNS methylation spectrum (n = 678, 28 methylation families), on non-array platforms (whole-genome bisulfite sequencing and Oxford Nanopore), and in a focused analysis of the low-grade glioma and diffuse midline glioma compartment across four independent cohorts (n = 670). Results: In the discovery cohort, MAGIBU achieved high concordance with the consensus reference (Cohen's κ = 0.855), outperforming EpiDiP (κ = 0.278), which frequently placed low-grade tumors in proximity to higher-grade reference regions. Conclusions: MAGIBU provides a stable, quantitative differential diagnosis framework that mitigates the limitations of rigid categorical assignments. By leveraging a distance-based proximity metric, it offers a transparent decision-support tool that integrates effectively with clinical, radiological, and molecular data. While performance is inherently dependent on reference atlas composition, MAGIBU represents a robust complementary approach for the diagnostic workup of ambiguous CNS tumors.

Brain↗

Sequence-specific detection of individual DNA strands using engineered nanopores.

We describe biosensor elements that are capable of identifying individual DNA strands with single-base resolution. Each biosensor element consists of an individual DNA oligonucleotide covalently attached within the lumen of the alpha-hemolysin (alphaHL) pore to form a "DNA-nanopore". The binding of single-stranded DNA (ssDNA) molecules to the tethered DNA strand causes changes in the ionic current flowing through a nanopore. On the basis of DNA duplex lifetimes, the DNA-nanopores are able to discriminate between individual DNA strands up to 30 nucleotides in length differing by a single base substitution. This was exemplified by the detection of a drug resistance-conferring mutation in the reverse transcriptase gene of HIV. In addition, the approach was used to sequence a complete codon in an individual DNA strand tethered to a nanopore.

Base Pair Mismatch↗

Chromosome-level genome assembly and annotation of Pterygoplichthys pardalis.

Suckermouth catfishes, with their evolved powerful features, have become notorious invasive species, causing significant damage to aquatic ecosystems. However, the lack of high-quality genomes severely restricts research on this group within the field. In this study, we de novo assembled the chromosome-level genome assembly of Pterygoplichthys pardalis using multiple platforms of sequencing data, including Illumina short reads, Nanopore long reads, and Hi-C sequencing reads, resulting in a 1.51 Gb genome assembly. Multiple evaluations, including read mapping ratio (98.52%), transcript mapping ratio (99.61%), conserved BUSCO gene set (98.8%), and N50 score (49.47 Mb), indicated the high continuity and accuracy of the genome assembly we generated. Genome annotation found that 0.97 Gb of genome sequences are repetitive sequences, accounting for 64.47% of the genome assembly. Further, 23,859 protein-coding genes were successfully predicted, 92.92% of which could be annotated in functional databases. This high-quality genome assembly of P. pardalis provides a valuable resource for understanding the genetic underpinnings of P. pardalis's invasive success and offers critical data for future fisheries research and management.

Animals↗

Toward single molecule DNA sequencing: direct identification of ribonucleoside and deoxyribonucleoside 5'-monophosphates by using an engineered protein nanopore equipped with a molecular adapter.

Individual nucleic acid molecules might be sequenced by the identification of nucleoside 5'-monophosphates as they are released by processive exonucleases. Here, we show that single molecule detection with a modified protein nanopore can be used to identify ribonucleoside and 2'-deoxyribonucleoside 5'-monophosphates, thereby taking a step along this path. Distinct levels of current block are observed for each of the four members of a set of nucleoside 5'-monophosphates when the molecules bind within a mutant alpha-hemolysin pore, (M113R)(7), equipped with the molecular adapter heptakis-(6-deoxy-6-amino)-beta-cyclodextrin. While our results compare favorably with alternative approaches, further work will be required to improve the accuracy of identification of the nucleic acid bases, to feed each released nucleotide into the pore, and to ensure that every nucleotide is captured by the adapter.

DNA↗

Population-scale detection of methylation outliers from long-read genome sequencing.

BACKGROUND: Aberrant DNA methylation can mediate the functional effects of rare genetic variation and contribute to imprinting disorders, repeat expansion diseases, and other pathogenic regulatory mechanisms. Long-read sequencing technologies now enable genome-wide detection of CpG methylation alongside genetic variation from a single assay. However, methods for systematic identification and interpretation of methylation outliers from long-read sequencing data remain limited. METHODS: We developed METAFORA, a computational workflow for detecting methylation outlier regions from PacBio and Oxford Nanopore long-read sequencing data. METAFORA constructs population-level methylation references, segments the genome into correlated CpG blocks, infers technical and biological sources of variation through hidden factor estimation, models uncertainty due to variable depth sequencing, and computes covariate-adjusted methylation outlier scores for individual samples. We applied METAFORA across large long-read sequencing cohorts and integrated methylation outliers with multi-omic data. METAFORA is implemented as a snakemake workflow available at https://github.com/tjense25/METAFORA. RESULTS: METAFORA identified methylation outlier regions associated with rare structural variants, tandem repeat expansions, and imprinting abnormalities. We found outlier regions were enriched for molecular outliers across transcriptomic and chromatin accessibility datasets, supporting their functional relevance in gene regulation. In a representative case, METAFORA identified an imprinting defect affecting the GNAS locus associated with an STX16 deletion. CONCLUSIONS: METAFORA enables scalable detection and interpretation of methylation outliers from long-read sequencing data and provides a framework for integrating epigenetic outliers with genomic and multi-omic analyses. These approaches may improve interpretation of rare regulatory variation and support discovery of clinically relevant epigenetic abnormalities in genomic medicine.

DNA methylation↗

The chromosome-level genome of Stylosanthes guianensis provides insights into genome evolution and environmental adaptation.

Stylosanthes guianensis is a leguminous forage crop of significant economic importance, primarily distributed in tropical and subtropical regions. It exhibits strong adaptability to various stresses, yet the genetic basis underlying this trait remains unclear. In this study, we constructed the first chromosome-scale reference genome of S. guianensis using a combination of Nanopore and Hi-C sequencing technologies. The assembled genome size is 1254 Mb, with 10 pseudochromosomes. Using Nanopore full-length transcriptome data, we generated high-quality transcript-level gene annotations, identifying 36 585 gene models and 110 601 transcripts. The repetitive sequences in S. guianensis account for 79.16% of the genome, with the extensive expansion of Gypsy elements in long terminal repeats contributing to its genome size enlargement. Comparative genomic and transcriptomic analyses revealed that flavonoid metabolism plays a pivotal role in stress adaptation, providing new insights into the genetic basis of stress tolerance. Additionally, we generated whole-genome methylation profiles under cold treatment and control conditions, offering valuable data for future epigenomic research. These findings provide essential molecular resources for understanding stress resilience in S. guianensis and advancing its molecular breeding.

Genome, Plant↗