ABSTRACT Metagenomics enables comprehensive exploration of microbial communities but is influenced by library preparation and sequencing technologies, affecting recovery of microbial genomes and proteins. Here, we benchmarked six Illumina-compatible short-read library preparation conditions in triplicate at 2 × 150 bp and 2 × 250 bp read lengths alongside PacBio HiFi long-read sequencing using a composite environmental sample of marine mangrove sediment and terrestrial palm tree soil. Longer short reads (2 × 250 bp) combined with optimal library preparation approaches improved assembly quality, protein detection, and metagenome-assembled genome (MAG) recovery, achieving results approaching those of long-read sequencing. TruSeq libraries at 2 × 250 bp recovered more than sevenfold more unique proteins than the same kit at 2 × 150 bp (811,701 vs 110,108) using the same number of sequencing reads, while recovering a comparable number of high-quality MAGs to PacBio HiFi long-read sequencing (11 vs 18) and surpassing it in protein discovery by almost 10-fold (811,701 vs 87,745) at less than half of the sequencing cost. Furthermore, biosynthetic gene cluster analysis identified 46 biosynthetic gene clusters in TruSeq-250PE assemblies compared to 38 in PacBio HiFi, with several showing no close match in the MIBiG database. Although long reads yield more contiguity and complete genomes, longer short reads offer a cost-effective, scalable alternative for uncovering microbial and functional diversity. These findings provide critical guidance for metagenomic experimental design, demonstrating that strategic selection of library preparation chemistry and sequencing parameters can reveal more unknown microbial information in complex biomes without requiring additional sequencing depth. IMPORTANCE Metagenomic outcomes are strongly influenced by library preparation and sequencing strategies, yet their combined effects in complex environmental samples remain poorly defined. Here, we provide the first direct comparison of Illumina NovaSeq short-read metagenomic sequencing at 2 × 150 bp and 2 × 250 bp across multiple library preparation kits, alongside PacBio HiFi long-read sequencing. We show that sequencing read length and library preparation critically shape assembly quality, protein recovery, and metagenome-assembled genome (MAG) reconstruction. These findings demonstrate that short-read sequencing at 2 × 250 bp, with appropriate library preparation, can match long-read technologies in MAG recovery while substantially surpassing them in protein discovery. With less than half of the sequencing price and a 3.5-fold reduction in cost per gigabase of usable data, this method facilitates more accessible large-scale metagenomic analysis within complex environmental systems.

Metagenomics enables comprehensive exploration of microbial communities but is influenced by library preparation and sequencing technologies, affecting recovery of microbial genomes and proteins. Here, we benchmarked six Illumina short-read library kits at 2×150 bp and 2×250 bp read lengths alongside PacBio HiFi long-read sequencing using a complex environmental sample. Longer short reads (2×250 bp) combined with optimal library preparation approaches notably improved assembly quality, protein detection, and metagenome-assembled genome (MAG) recovery, achieving results similar to those of long-read sequencing. Although long reads yield more contiguous and complete genomes, longer short reads offer a cost-effective, scalable alternative for uncovering microbial and functional diversity. These findings provide critical guidance for metagenomic experimental design, demonstrating the importance of strategic selection of library preparation chemistry and sequencing parameters in revealing unknown microbial information in complex biomes without requiring additional sequencing depth.

In this study, xTitan, a field-deployable, automated, and versatile nucleic acid extraction system was employed to characterize microbial communities in Red Sea-derived samples, including coral colonies, mangrove sediments, and seawater. The use of the xTitan in the field was intended to minimize sample transport bias, obtaining data that may be closer to “ground truth” for microbial diversity. The observed microbial communities from DNA extracted in the field using the xTitan system were compared to DNA extractions performed in a laboratory setting using both xTitan and a standard commercial kit (Qiagen) after approximately 24 h of sample transfer and storage. Microbial community analyses conducted on DNA extracted using the xTitan system and the Qiagen kit yielded similar alpha diversity metric values, with a trend toward higher diversity observed in most samples extracted with the xTitan. The microbial community structure in samples from a Pocillopora verrucosa colony, mangrove sediments, and seawater was affected by the DNA extraction system. In the P. verrucosa colony, 16S rRNA gene sequences affiliated to Endozoicomonas acroporae were preferentially abundant when DNA was extracted in the field with the xTitan system rather than in the lab. In mangrove sediments, significant differences (P-value < 0.05) in beta diversity and functional gene profiles were observed when comparing in-field to in-lab xTitan DNA extracts. In seawater, a pronounced decrease in the relative abundance of cyanobacterial populations was observed when DNA was extracted with both methods after samples were transported to the lab on ice. In addition, hundreds of species in mangrove-associated samples were differentially abundant when DNA was extracted on-site with the xTitan system compared to in-lab extractions. Balneolaceae was one of the most abundant taxa in mangrove sediments and several genera from this family were detected in all replicates across all DNA extraction systems. The usability of different field-deployable instruments for microbial community characterization in marine-derived samples was demonstrated. Moreover, differences in beta diversity were observed when DNA was extracted in-field versus in-lab using the xTitan system, particularly for mangrove-associated samples. These results highlight the value of on-site nucleic acid extraction for enhancing the detection of microbial taxa that can be sensitive to cold storage. This study enabled the testing of the xTitan on Red Sea-derived samples, generating comprehensive information on the effects of DNA extraction systems and transportation of samples on coral and mangrove-associated microbiomes.

We have used multiple sequencing approaches to sequence the genome of a volunteer from Saudi Arabia. We use the resulting data to generate a de novo assembly of the genome, and use different computational approaches to refine the assembly. As a consequence, we provide a contiguous assembly of the complete genome of an individual from Saudi Arabia for all chromosomes except chromosome Y, and label this assembly KSA001. We transferred genome annotations from reference genomes to fully annotate KSA001, and we make all primary sequencing data, the assembly, and the genome annotations freely available in public databases using the FAIR data principles. KSA001 is the first telomere-to-telomere-assembled genome from a Saudi individual that is freely available for any purpose.

Cancer is a multifaceted disease arising from numerous genomic aberrations that have been identified as a result of advancements in sequencing technologies. While next-generation sequencing (NGS), which uses short reads, has transformed cancer research and diagnostics, it is limited by read length. Third-generation sequencing (TGS), led by the Pacific Biosciences and Oxford Nanopore Technologies platforms, employs long-read sequences, which have marked a paradigm shift in cancer research. Cancer genomes often harbour complex events, and TGS, with its ability to span large genomic regions, has facilitated their characterisation, providing a better understanding of how complex rearrangements affect cancer initiation and progression. TGS has also characterised the entire transcriptome of various cancers, revealing cancer-associated isoforms that could serve as biomarkers or therapeutic targets. Furthermore, TGS has advanced cancer research by improving genome assemblies, detecting complex variants, and providing a more complete picture of transcriptomes and epigenomes. This review focuses on TGS and its growing role in cancer research. We investigate its advantages and limitations, providing a rigorous scientific analysis of its use in detecting previously hidden aberrations missed by NGS. This promising technology holds immense potential for both research and clinical applications, with far-reaching implications for cancer diagnosis and treatment.

The tree Eucalyptus camaldulensis is a ubiquitous member of the Eucalyptus genus, which includes several hundred species. Despite the extensive sequencing and assembly of nuclear genomes from various eucalypts, the genus has only one fully annotated and complete mitochondrial genome (mitogenome). Plant mitochondria are characterized by dynamic genomic rearrangements, facilitated by repeat content, a feature that has hindered the assembly of plant mitogenomes. This complexity is evident in the paucity of available mitogenomes. This study, to the best of our knowledge, presents the first E. camaldulensis mitogenome. Our findings suggest the presence of multiple isomeric forms of the E. camaldulensis mitogenome and provide novel insights into minor rearrangements triggered by nested repeat sequences. A comparative sequence analysis of the E. camaldulensis and E. grandis mitogenomes unveils evolutionary changes between the two genomes. A significant divergence is the evolution of a large repeat sequence, which may have contributed to the differences observed between the two genomes. The largest repeat sequences in the E. camaldulensis mitogenome align well with significant yet unexplained structural variations in the E. grandis mitogenome, highlighting the adaptability of repeat sequences in plant mitogenomes.

Colorectal cancer (CRC) ranks as the second leading cause of cancer deaths globally. In recent years, short-read single-cell RNA sequencing (scRNA-seq) has been instrumental in deciphering tumor heterogeneities. However, these studies only enable gene-level quantification but neglect alterations in transcript structures arising from alternative end processing or splicing. In this study, we integrated short- and long-read scRNA-seq of CRC samples to build an isoform-resolution CRC transcriptomic atlas. We identified 394 dysregulated transcript structures in tumor epithelial cells, including 299 resulting from various combinations of splicing events. Second, we characterized genes and isoforms associated with epithelial lineages and subpopulations exhibiting distinct prognoses. Among 31,935 isoforms with novel junctions, 330 were supported by The Cancer Genome Atlas RNA-seq and mass spectrometry data. Finally, we built an algorithm that integrated novel peptides derived from open reading frames of recurrent tumor-specific transcripts with mass spectrometry data and identified recurring neoepitopes that may aid the development of cancer vaccines.

Background Long-read whole genome sequencing (lrWGS) has the potential to address the technical limitations of exome sequencing in ways not possible by short-read WGS. However, its utility in autosomal recessive Mendelian diseases is largely unknown. Methods In a cohort of 34 families in which the suspected autosomal recessive diseases remained undiagnosed by exome sequencing, lrWGS was performed on the Pacific Bioscience Sequel IIe platform. Results Likely causal variants were identified in 13 (38%) of the cohort. These include (1) a homozygous splicing SV in TYMS as a novel candidate gene for lethal neonatal lactic acidosis, (2) a homozygous non-coding SV that we propose impacts STK25 expression and causes a novel neurodevelopmental disorder, (3) a compound heterozygous SV in RP1L1 with complex inheritance pattern in a family with inherited retinal disease, (4) homozygous deep intronic variants in LEMD2 and SNAP91 as novel candidate genes for neurodevelopmental disorders in two families, and (5) a promoter SNV in SLC4A4 causing non-syndromic band keratopathy. Surprisingly, we also encountered causal variants that could have been identified by short-read exome sequencing in 7 families. The latter highlight scenarios that are especially challenging at the interpretation level. Conclusions Our data highlight the continued need to address the interpretation challenges in parallel with efforts to improve the sequencing technology itself. We propose a path forward for the implementation of lrWGS sequencing in the setting of autosomal recessive diseases in a way that maximizes its utility.

Recent years have witnessed a rapid development of sequencing technologies. Fundamental differences and limitations among various platforms impact the time, the cost and the accuracy for sequencing whole genomes. Here we designed a complete de novo plant genome generation workflow that starts from plant tissue samples and produces high-quality draft genomes with relatively modest laboratory and bioinformatic resources within seven days. To optimize our workflow we selected different species of plants which were used to extract high molecular weight DNA, to make PacBio and ONT libraries for sequencing with the Sequel I, Sequel II and GridION platforms. We assembled high-quality draft genomes of two different Eucalyptus species E. rudis , and E. camaldulensis to chromosome level without using additional scaffolding technologies. For the rapid production of de novo genome assembly of plant species we showed that our DNA extraction protocol followed by PacBio high fidelity sequencing, and assembly with new generation assemblers such as hifiasm produce excellent results. Our findings will be a valuable benchmark for groups planning wet- and dry-lab plant genomics research and for high throughput plant genomics initiatives.

Background The consequences of the earth’s daily rotation have led to 24-hour biological rhythms in most organisms. Even parasites have daily rhythms, which, when in synchrony with host rhythms, optimize their fitness. Understanding these rhythms may enable development of novel control strategies that take advantage of rhythmic vulnerabilities. Recent work on blood-dwelling protozoan parasites has revealed daily rhythms in gene expression, physiology, drug sensitivity and the presence of an intrinsic circadian clock. However, similar studies on metazoan parasites are lacking. The aims of this study were to investigate if a metazoan parasite has daily molecular oscillations, whether they give us insight into how these longer-lived organisms can survive host daily cycles over a life-span of many years and to determine whether canonical metazoan circadian clock genes are present and rhythmic. We addressed these questions using the human blood fluke Schistosoma mansoni, that lives in the vasculature for decades and causes the serious neglected tropical disease schistosomiasis. Results Using round-the-clock transcriptomics of male and female adult worms we discovered that ∼2% of its genes followed a daily pattern of expression. Rhythmic processes included a night-time stress response and a day-time metabolic ‘rush hour’. Transcriptional profiles in the female reproductive system were mirrored by daily patterns in egg laying (eggs are the main drivers of the host pathology). Genes cycling with the highest amplitudes include drug targets and a vaccine candidate. These 24hr rhythms may be driven by host rhythms and/or generated by a circadian clock. However, core clock genes are missing and orthologues of secondary clock genes show no 24hr rhythmicity in transcript abundance. Conclusions The daily rhythms identified here reveal temporally compartmentalised internal processes and host interactions over the daily cycle, including processes relevant to within-host survival and between-host transmission. Our findings suggest that if these daily rhythms are generated by an intrinsic circadian clock then the oscillatory mechanism must be distinct from that in other Metazoa. Most importantly, knowing which gene transcripts oscillate at this temporal scale is relevant to functional genomic studies that will lead to the development and delivery of therapeutics against schistosomiasis.

Currently, different sequencing platforms are used to generate plant genomes and no workflow has been properly developed to optimize time, cost, and assembly quality. We present LeafGo, a complete de novo plant genome workflow, that starts from tissue and produces genomes with modest laboratory and bioinformatic resources in approximately 7 days and using one long-read sequencing technology. LeafGo is optimized with ten different plant species, three of which are used to generate high-quality chromosome-level assemblies without any scaffolding technologies. Finally, we report the diploid genomes of Eucalyptus rudis and E. camaldulensis and the allotetraploid genome of Arachis hypogaea.

Background. At least 250 million people worldwide suffer from schistosomiasis, caused by Schistosoma worms. Genome sequences for several Schistosoma species are available, including a high-quality annotated reference for Schistosoma mansoni. There is a pressing need to develop a reliable functional toolkit to translate these data into new biological insights and targets for intervention. CRISPR-Cas9 was recently demonstrated for the first time in S. mansoni, to produce somatic mutations in the omega-1 ( ω1) gene. Methods. We employed CRISPR-Cas9 to introduce somatic mutations in a second gene, SULT-OR, a sulfotransferase expressed in the parasitic stages of S. mansoni, in which mutations confer resistance to the drug oxamniquine. A 262-bp PCR product spanning the region targeted by the gRNA against SULT-OR was amplified, and mutations identified in it by high-throughput sequencing. Results. We found that 0.3-2.0% of aligned reads from CRISPR-Cas9-treated adult worms showed deletions spanning the predicted Cas9 cut site, compared to 0.1-0.2% for sporocysts, while deletions were extremely rare in eggs. The most common deletion observed in adults and sporocysts was a 34 bp-deletion directly upstream of the predicted cut site, but rarer deletions reaching as far as 102 bp upstream of the cut site were also detected. The CRISPR-Cas9-induced deletions, if homozygous, are predicted to cause resistance to oxamniquine by producing frameshifts, ablating SULT-OR transcription, or leading to mRNA degradation via the nonsense-mediated mRNA decay pathway. However, no SULT-OR knock down at the mRNA level was observed, presumably because the cells in which CRISPR-Cas9 did induce mutations represented a small fraction of all cells expressing SULT-OR. Conclusions. Further optimisation of CRISPR-Cas protocols for different developmental stages and particular cell types, including germline cells, will contribute to the generation of a homozygous knock-out in any gene of interest, and in particular the SULT-OR gene to derive an oxamniquine-resistant stable transgenic line.

Abstract Background The consequences of the earth’s daily rotation have led to 24-h biological rhythms in most organisms. Even some parasites are known to have daily rhythms, which, when in synchrony with host rhythms, can optimise their fitness. Understanding these rhythms may enable the development of control strategies that take advantage of rhythmic vulnerabilities. Recent work on protozoan parasites has revealed 24-h rhythms in gene expression, drug sensitivity and the presence of an intrinsic circadian clock; however, similar studies on metazoan parasites are lacking. To address this, we investigated if a metazoan parasite has daily molecular oscillations, whether they reveal how these longer-lived organisms can survive host daily cycles over a lifespan of many years and if animal circadian clock genes are present and rhythmic. We addressed these questions using the human blood fluke Schistosoma mansoni that lives in the vasculature for decades and causes the tropical disease schistosomiasis. Results Using round-the-clock transcriptomics of male and female adult worms collected from experimentally infected mice, we discovered that ~ 2% of its genes followed a daily pattern of expression. Rhythmic processes included a stress response during the host’s active phase and a ‘peak in metabolic activity’ during the host’s resting phase. Transcriptional profiles in the female reproductive system were mirrored by daily patterns in egg laying (eggs are the main drivers of the host pathology). Genes cycling with the highest amplitudes include predicted drug targets and a vaccine candidate. These 24-h rhythms may be driven by host rhythms and/or generated by a circadian clock; however, orthologs of core clock genes are missing and secondary clock genes show no 24-h rhythmicity. Conclusions There are daily rhythms in the transcriptomes of adult S. mansoni, but they appear less pronounced than in other organisms. The rhythms reveal temporally compartmentalised internal processes and host interactions relevant to within-host survival and between-host transmission. Our findings suggest that if these daily rhythms are generated by an intrinsic circadian clock then the oscillatory mechanism must be distinct from that in other animals. We have shown which transcripts oscillate at this temporal scale and this will benefit the development and delivery of treatments against schistosomiasis.

AbstractCRISPR/Cas9-mediated genome editing shows cogent potential for the genetic modification of helminth parasites. Here we report successful gene knock-in (KI) into the genome of the egg ofSchistosoma mansoniby combining CRISPR/Cas9 with single-stranded oligodeoxynucleotides (ssODNs). We edited the acetylcholinesterase (AChE) gene ofS. mansonitargeting two guide RNAs (gRNAs), X5 and X7, located on exon 5 and exon 7 of Smp_154600, respectively. A CRISPR/Cas9-vector encoding gRNA X5 or X7 was introduced by electroporation into eggs recovered from livers of experimentally infected mice. Simultaneously, eggs were transfected with a ssODN donor encoding a stop codon in all six frames, flanked by 50 nt-long 5’- and 3’-homology arms matching the predicted Cas9-catalyzed double stranded break at X5 or X7. Next generation sequencing analysis of reads of amplicon libraries spanning targeted regions revealed that the major modifications induced by CRISPR/Cas9 in the eggs were generated by homology directed repair (HDR). Furthermore, soluble egg antigen from AChE-edited eggs exhibited markedly reduced AChE activity, indicative that programmed Cas9 cleavage mutated the AChE gene. Following injection of AChE-edited schistosome eggs into the tail veins of mice, a significant decrease in circumoval granuloma size was observed in the lungs of the mice. Notably, there was an enhanced Th2 response involving IL-4, −5, −10, and-13 induced by lung cells and splenocytes in mice injected with X5-KI eggs in comparison to control mice injected with unmutated eggs. A Th2-predominant response, with increased levels of IL-4, −13 and GATA3, also was induced by X5 KI eggs in small intestine-draining mesenteric lymph node cells when the gene-edited eggs were introduced into the subserosa of the ileum of the mice. These findings confirmed the potential and the utility of CRISPR/Cas9-mediated genome editing for functional genomics in schistosomes.Author SummarySchistosomiasis is the most devastating of the parasitic helminth diseases. Currently, no vaccines are available for human use and praziquantel is the only available treatment raising considerable concern that drug resistance will develop. A major challenge faced by the schistosomiasis research community is the lack of suitable tools to effectively characterise schistosome gene products as potential new drug and/or vaccine targets. We introduced CRISPR/Cas9 mediated editing intoS. mansonieggs targeting the gene encoding acetylcholinesterase (AChE), a recognized anthelminthic drug target. We found that the major modifications induced by CRISPR/Cas9 in the eggs were generated by homology directed repair (HDR). This platform provides a unique opportunity to generate precise loss-of-function insertions into the schistosome genome. We pre-screened the activity of two guide RNAs of the AChE gene and compared/validated the mutation efficacy using next-generation sequencing analysis at the genomic level and phenotypic modifications at the protein level. That resulted in reduced AChE activity observed in AChE-edited eggs, and decreased lung circumoval granuloma size in mice injected with those edited eggs. The CRISPR/Cas9-genome editing system we established in this study provides a pivotal platform for gene functional studies to identify and test new anti-schistosome intervention targets, which can be extended to the other human schistosome species and other important parasitic helminths.

Despite significant progress, China faces the challenge of re-emerging schistosomiasis transmission in currently controlled areas due, in part, to the presence of a range of animal reservoirs, notably water buffalo and cattle, which can harbor Schistosoma japonicum infections. Environmental, ecological and social-demographic changes in China, shown to affect the distribution of oncomelanid snails, can also impact future schistosomiasis transmission. In light of their importance in the S. japonicum, lifecycle, vaccination has been proposed as a means to reduce the excretion of egg from cattle and buffalo, thereby interrupting transmission from these reservoir hosts to snails. A DNA-based vaccine (SjCTPI) our team developed showed encouraging efficacy against S. japonicum in Chinese water buffaloes. Here we report the results of a double-blind cluster randomized trial aimed at determining the impact of a combination of the SjCTPI bovine vaccine (given as a prime-boost regimen), human mass chemotherapy and snail control on the transmission of S. japonicum in 12 selected administrative villages around the Dongting Lake in Hunan province. The trial confirmed human praziquantel treatment is an effective intervention at the population level. Further, mollusciciding had an indirect ~50% efficacy in reducing human infection rates. Serology showed that the SjCTPI vaccine produced an effective antibody response in vaccinated bovines, resulting in a negative correlation with bovine egg counts observed at all post-vaccination time points. Despite these encouraging outcomes, the effect of the vaccine in preventing human infection was inconclusive. This was likely due to activities undertaken by the China National Schistosomiasis Control Program, notably the treatment, sacrifice or removal of bovines from trial villages, over which we had no control; as a result, the trial design was compromised, reducing power and contaminating outcome measures. This highlights the difficulties in undertaking field trials of this nature and magnitude, particularly over a long period, and emphasizes the importance of mathematical modeling in predicting the potential impact of control intervention measures. A transmission blocking vaccine targeting bovines for the prevention of S. japonicum with the required protective efficacy would be invaluable in tandem with other preventive intervention measures if the goal of eliminating schistosomiasis from China is to become a reality.

We have developed a method based on the Qiagen Genomics-Tip kit for the extraction of HMW DNA from a wide range of plants suitable for third generation sequencing platforms such as PacBio and Oxford Nanopore. This HMW plant DNA was used to make sequencing libraries with high insert N50 (30-40 Kb) for GridION and PromethION Oxford Nanopore machines and Sequel I and Sequel II PacBio machines. This HMW DNA is also suitable for PacBio HiFi libraries. The resulting sequencing data is ideally suited for de novo assembly of plant genomes.

Schistosomiasis is an important parasitic disease, touching roughly 200 million people worldwide. The causative agents are different Schistosoma species. Schistosomes have a complex life cycle, with a freshwater snail as intermediate host. After infection, sporocysts develop inside the snail host and give rise to human dwelling larvae. We present here a detailed step-by-step video instruction in English, French, Spanish and Portuguese that shows how these sporocysts can be manipulated and transferred from one snail to another. This procedure provides a technical basis for different types of ex vivo modifications, such as those used in functional genomics studies.

CRISPR/Cas9-based genome editing has yet to be reported in species of the Platyhelminthes. We tested this approach by targeting omega-1 (ω1) of Schistosoma mansoni as proof of principle. This secreted ribonuclease is crucial for Th2 polarization and granuloma formation. Schistosome eggs were exposed to Cas9 complexed with guide RNA complementary to ω1 by electroporation or by transduction with lentiviral particles. Some eggs were also transfected with a single stranded donor template. Sequences of amplicons from gene-edited parasites exhibited Cas9-catalyzed mutations including homology directed repaired alleles, and other analyses revealed depletion of ω1 transcripts and the ribonuclease. Gene-edited eggs failed to polarize Th2 cytokine responses in macrophage/T-cell co-cultures, while the volume of pulmonary granulomas surrounding ω1-mutated eggs following tail-vein injection into mice was vastly reduced. Knock-out of ω1 and the diminished levels of these cytokines following exposure showcase the novel application of programmed gene editing for functional genomics in schistosomes.