Sequence Reads Research Articles

Deep learning revealed the distribution and evolution patterns for invertible promoters across bacterial lineages.

Invertible promoters (invertons) are crucial regulatory elements in bacteria, facilitating gene expression changes under stress. Despite their importance, their prevalence and the range of regulated gene functions are largely unknown. We introduced DeepInverton, a deep learning model that identifies invertons across a broad phylogenetic spectrum without using sequencing reads. By analyzing 68 733 bacterial genomes and 9382 metagenomes, we have uncovered over 200 000 nonredundant invertons and have also highlighted their abundance in pathogens. Additionally, we identified a post-Cambrian Explosion increase of invertons, paralleling species diversification. Furthermore, we revealed that invertons regulate diverse functions, including antimicrobial resistance and biofilm formation, underscoring their role in environmental adaptation. Notably, the majority of inverton identifications by DeepInverton have been confirmed by the in vitro experiments. The comprehensive inverton profiles have deepened our understanding of invertons at pan-genome and pan-metagenome scales, enabling a broad spectrum of applications in microbial ecology and synthetic biology.

Alternative transcripts recode human genes to express overlapping, frameshifted microproteins.

Overlapping genes were thought to be essentially absent from the human genome until the discovery of abundant, frameshifted internal open reading frames (iORFs) nested within annotated protein coding sequences. However, it is currently unclear how many functional human iORFs exist and how they are expressed. We demonstrate that, in hundreds of cases, alternative transcript variants that bypass the start codon of annotated coding sequences (CDSs) can recode a human gene to express the iORF-encoded microprotein. While many human genes generate such non- coding alternative transcripts, they are poorly annotated. Here we develope a new analysis pipeline enabling the assignment of translated human iORFs to alternative transcripts, and provide long- read sequencing and molecular validation of their expression in dozens of cases. Finally, we demonstrate that a conserved DEDD2 iORF switches the function of this gene from pro- to anti- apoptotic. This work thus demonstrates that alternative transcript variants can broadly reprogram human genes to express frameshifted iORFs, revealing new levels of complexity in the human transcriptome and proteome.

Canine Multiple System Degeneration Associated with Sequence Variants in SERAC1.

Canine multiple system degeneration (CMSD) is an early onset, progressive movement disorder affecting Kerry Blue Terriers and Chinese Crested dogs. The associated pathologic lesions include degeneration of the cerebellum, caudate nucleus, and substantia nigra. CMSD is inherited as an autosomal recessive trait in both dog breeds. Previous linkage mapping localized the CMSD locus to a 15 MB region on canine chromosome 1 (CFA1). Next-generation sequencing was used to generate whole-genome sequences from the DNA of an affected dog from each breed. The resulting sequence reads were aligned to the NCBI canine reference genome (build 3.1). Among the homozygous sequence variants within the CFA1 target region, a nonsense variant in exon 15 of SERAC1 was identified in the affected Kerry Blue Terrier, while in the Chinese Crested dog, a 4 bp deletion in the SERAC1 exon 4 acceptor splice site was found. RT-PCR showed that this deletion resulted in exon 4 skipping. Genotyping of large cohorts of Kerry Blue Terriers and Chinese Crested dogs for the respective breed-specific SERAC1 variants showed complete concordance between genotype and disease phenotype. Genotype-phenotype concordance was also observed in offspring generated by cross breeding between SERAC1-heterozygous Kerry Blue Terrier and Chinese Crested dogs, with only the compound heterozygotes exhibiting the disease phenotype, further confirming the recessive inheritance of CMSD. Variants in human SERAC1 are associated with disorders with a range of ages of disease onset and patterns of clinical signs, but that are all characterized by movement abnormalities similar to those of the dogs with CMSD. Canine CMSD could serve as a valuable model to elucidate the mechanisms underlying SERAC1-deficiency disorders and to evaluate potential therapeutic interventions.

A first report of rotavirus B from Zambian pigs leading to the discovery of a novel VP4 genotype P[9]

BackgroundRotavirus B (RVB) causes diarrhea in humans and pigs. Although various RVB strains were identified in humans and various animals globally, little is known about the epidemiology RVB infection in Africa. In this study, we attempted to examine the prevalence of RVB infection in pig populations in Zambia.MethodsMetagenomic analyses were conducted on pig feces collected in Zambia to detect double stranded RNA viruses, including RVB. To clarify the prevalence of RVB infection in pig populations in Zambia, 147 fecal samples were screened for the RVB detection by RT-qPCR. Full genome sequence of a detected RVB was determined by Sanger sequencing and genetically analyzed.ResultsThe metagenomic analyses revealed that RVB sequence reads and contigs of RVB were detected from one fecal sample collected from pigs in Zambia. RT-qPCR screening detected RVB genomes in 36.7% (54/147) of fecal samples. Among 54 positive samples, 13 were positive in non-diarrheal samples (n = 48, 27.1%) and 41 in diarrheal samples (n = 99, 41.4%). Genetic analyses demonstrated that all the segments of ZP18-18, except for VP4, had high nucleotide sequence identities (80.6–92.6%) with all other known RVB strains detected in pigs. In contrast, the VP4 sequence of ZP18-18 was highly divergent from other RVB strains (< 64.6% identities) and formed a distinct lineage in the phylogenetic tree. Notably, the VP8 subunit of the VP4 showed remarkably low amino acid identities (33.3%) to those of known RVB strains, indicating that the VP8 subunit of ZP18-18 was unique among RVB strains. According to the whole genome classification for RVB, ZP18-18 was assigned to a genotype constellation, G18-P[9]-I12-R4-C4-M4-A8-N10-T5-E4-H7 with the newly established VP4 genotype P[9].ConclusionsThis current study updates the geographical distribution and the genetic diversity of RVB. Given the lack of information regarding RVB in Africa, further RVB surveillance is required to assess the potential risk to humans and animals.

Estimating and correcting index hopping misassignments in single-cell RNA-seq data.

Index hopping causes read assignment errors in data from multiplexed sequencing libraries. This issue has become more prevalent with the widespread use of high-capacity sequencers and highly multiplexed single-cell RNA sequencing (scRNA-seq) libraries. We conducted deep, plate-based scRNA-seq on a mixed population of mouse skin cells. Analysis of transcriptomes from 1152 cells identified four distinct cell types. To estimate the error rate in sample assignment due to index hopping, we employed differential expression analysis to identify signature genes that were highly and specifically expressed in each cell type. We quantified the proportion of misassigned reads by examining the detection rates of signature genes in other cell types. Remarkably, regardless of gene expression levels, we estimated that 0.65% of reads per gene were assigned to incorrect cell across our data. To computationally compensate for index hopping, we developed a simple correction method wherein, for each gene, 0.65% of the library's average expression level was subtracted from the expression of each cell. This correction had notable effects on transcriptome analyses, including increased cell-cell clustering distance and alterations in intermediate state assignments of cell differentiation. These findings underscore the potential impact of index hopping on experimental results. In conclusion, we devised a straightforward method to estimate and correct for the index hopping rate by quantifying misassigned genes in distinct cell types within an scRNA-seq library. This approach can be applied to any barcoded, multiplexed scRNA-seq library containing cells with distinct expression profiles, allowing for correction of the expression matrix before conducting biological analysis.

Novel variants and genotype-phenotype correlation in a multicentre cohort of GNE myopathy in China

BackgroundGlcNAc2-epimerase (GNE) myopathy is a rare autosomal recessive disorder caused by pathogenic variants in the GNE gene, which is essential for the sialic acid biosynthesis pathway.ObjectiveThis multi-centre study aimed to...

Candida albicans isolates contain frequent heterozygous structural variants and transposable elements within genes and centromeres.

The human fungal pathogen Candida albicans poses a significant burden on global health, causing high rates of mortality and antifungal drug resistance. C. albicans is a heterozygous diploid organism that reproduces asexually. Structural variants (SVs) are an important source of genomic rearrangement, particularly in species that lack sexual recombination. To comprehensively investigate SVs across clinical isolates of C. albicans, we conducted long read sequencing and genome-wide SV analysis in three distantly related clinical isolates. Our work included a new, comprehensive analysis of transposable element (TE) composition, location and diversity. SVs and TEs are frequently close to coding sequences and many SVs are heterozygous, suggesting that SVs might impact gene and allele-specific expression. Most SVs are uniquely present in only one clinical isolate, indicating that SVs represent a significant source of intra-species genetic variation. We identified multiple, distinct SVs at the centromeres of Chromosome 4 and Chromosome 5, including inversions and transposon polymorphisms. These two chromosomes are often aneuploid in drug resistant clinical isolates, and can form isochromosome structures with breakpoints near the centromere. Further screening of 100 clinical isolates confirmed the widespread presence of centromeric SVs in C. albicans, often appearing in a heterozygous state, indicating that SVs are contributing to centromere evolution in C. albicans Together, these findings highlight that SVs and TEs are common across diverse clinical isolates of C. albicans and that the centromeres of this organism are important sites of genome rearrangement.

Revisiting the Briggs Ancient DNA Damage Model: A Fast Maximum Likelihood Method to Estimate Post-Mortem Damage.

One essential initial step in the analysis of ancient DNA is to authenticate that the DNA sequencing reads are actually from ancient DNA. This is done by assessing if the reads exhibit typical characteristics of post-mortem damage (PMD), including cytosine deamination and nicks. We present a novel statistical method implemented in a fast multithreaded programme, ngsBriggs that enables rapid quantification of PMD by estimation of the Briggs ancient damage model parameters (Briggs parameters). Using a multinomial model with maximum likelihood fit, ngsBriggs accurately estimates the parameters of the Briggs model, quantifying the PMD signal from single and double-stranded DNA regions. We extend the original Briggs model to capture PMD signals for contemporary sequencing platforms and show that ngsBriggs accurately estimates the Briggs parameters across a variety of contamination levels. Classification of reads into ancient or modern reads, for the purpose of decontamination, is significantly more accurate using ngsBriggs than using other methods available. Furthermore, ngsBriggs is substantially faster than other state-of-the-art methods. ngsBriggs offers a practical and accurate method for researchers seeking to authenticate ancient DNA and improve the quality of their data.

Chromosome Structural Rearrangements in Invasive Haplodiploid Ambrosia Beetles Revealed by the Genomes of Euwallacea fornicatus (Eichhoff) and Euwallacea similis (Ferrari) (Coleoptera, Curculionidae, Scolytinae).

Bark and ambrosia beetles are among the most ecologically and economically damaging introduced plant pests worldwide. Life history traits including polyphagy, haplodiploidy, inbreeding polygyny, and symbiosis with fungi contribute to their dispersal and impact. Species vary in their interactions with host trees, with many attacking stressed or recently dead trees, such as the globally distributed Euwallacea similis (Ferrari). Other species, like the Polyphagous Shot Hole Borer Euwallacea fornicatus (Eichhoff), can attack over 680 host plants and is causing considerable economic damage in several countries. Despite their notoriety, publicly accessible genomic resources for Euwallacea Hopkins species are scarce, hampering our understanding of their invasive capabilities as well as modern control measures, surveillance, and management. Using a combination of long and short read sequencing platforms, we assembled and annotated high quality (BUSCO > 98% complete) pseudo-chromosome-level genomes for these species. Comparative macrosynteny analysis identified an increased number of pseudo-chromosome scaffolds in the haplodiploid inbreeding species of Euwallacea compared to diploid outbred species, due to fission events. This suggests that life history traits can impact chromosome structure. Further, the genome of E. fornicatus had a higher relative proportion of repetitive elements, up to 17% more, than E. similis. Metagenomic assembly pipelines identified microbiota associated with both species including Fusarium fungal symbionts and a novel Wolbachia strain. These novel genomes of haplodiploid inbreeding species will contribute to the understanding of how life history traits are related to their evolution and to the management of these invasive pests.

Transcriptional variation and RNA polymorphism among different Lentinula edodes (Berk.) Pegler strains.

Lentinula edodes is a commercially important mushroom known for its nutritional and therapeutic values. However, the molecular mechanisms underlying the distinct nutritional and physiological attributes of various L. edodes strains are not well understood. This study focused on three Lentinula strains (DMRO-356, DMRO-623, and DMRO-388s) with different nutritional and productivity profiles. Illumina sequencing was used to perform a whole-transcriptome analysis, conducting 100-base pair paired-end sequencing of total messenger RNA (mRNA) in duplicate, resulting in 28-48 million sequencing reads per strain. After rigorous data filtering, over 99% of high-quality reads were retained, and more than 95% were aligned to the Lentinula genome. Differential gene expression analyses identified 2210 differentially expressed genes between DMRO-356 and DMRO-623, 862 between DMRO-356 and DMRO-388s, and 2212 between DMRO-623 and DMRO-388s. Significant genetic variations were found among the strains, including 7753 single nucleotide polymorphisms (SNPs) in DMRO-356 versus DMRO-623 and 4080 SNPs in DMRO-356 versus DMRO-388s. Additionally, 349 insertions/deletions (InDels) were found in DMRO-356/DMRO-623 and 218 in DMRO-356/DMRO-388 s. Non-synonymous SNPs, which alter amino acid compositions, were analyzed, showing a preference for polar over charged amino acids. These differentially expressed genes were associated with various nutritional and developmental processes, highlighting the importance of genetic variations in shaping amino acid composition and potentially affecting protein function. This study is the first comprehensive exploration of transcriptional differences among Lentinula strains available for its cultivation, providing valuable insights to enhance mushroom quality and productivity. © 2024 Society of Chemical Industry.

5-hydroxymethylcytosine sequencing in plasma cell-free DNA identifies unique epigenomic features in prostate cancer patients resistant to androgen deprivation therapies.

Currently, no biomarkers are available to identify resistance to androgen-deprivation therapies (ADT) in men with hormone-naive prostate cancer. Since 5-hydroxymethylcytosines (5hmC) in gene body are associated with gene activation, in this study, we evaluated whether 5hmC signatures in cell-free DNA (cfDNA) predicts early resistance to ADT. We collected a total of 139 serial plasma samples from 55 prostate cancer patients receiving ADT at three time points including baseline (prior to initiating ADT, N=55), 3-month (after initiating ADT, N=55), and disease progression (N=15) within 24 months or 24-month if no progression was detected (N=14). To quantify 5hmC abundance across the genome, we used selective chemical labeling sequencing and mapped sequence reads to individual genes. Differential methylation analysis in baseline samples identified significant 5hmC difference in 1,642 of 23,433 genes between patients with and without progression (false discovery rate, FDR<0.1). Patients with disease progression showed significant 5hmC enrichments in multiple hallmark gene sets with androgen responses as top enriched gene set (FDR=1.19E-13). Interestingly, this enrichment was driven by a subgroup of patients featuring a significant 5hmC hypermethylation in the gene sets involving AR , FOXA1 and GRHL2 . To quantify overall activities of these gene sets, we developed a gene set activity scoring algorithm and observed significant association of high activity scores with poor progression-free survival (P<0.05). Longitudinal analysis showed that the high activity scores were significantly reduced after 3-months of initiating ADT (P<0.0001) but returned to higher levels when the disease was progressed (P<0.05). This study demonstrates that 5hmC-based activity scores from gene sets involved in AR , FOXA1 and GRHL2 may be used as biomarkers to determine early treatment resistance, monitor disease progression, and potentially identify patients who would benefit from upfront treatment intensification.

Genomic analyses of the southern and eastern yellowjacket wasps (Hymenoptera: Vespidae) reveal evolutionary signatures of social life

Abstract Insects have evolved remarkably complex social systems. Social wasps are particularly noteworthy because they display gradations in social behaviors. Here, we sequence the genomes of two highly diverged Vespula wasps, V. squamosa and V. maculifrons Buysson (Hymenoptera: Vespidae), to gain greater insight into the evolution of sociality. Both V. squamosa and V. maculifrons are social wasps that live in large colonies characterized by distinct queen and worker castes. However, V. squamosa is a facultative social parasite, and V. maculifrons is its frequent host. We found that the genomes of both species were ~200 Mbp in size, similar to the genome sizes of congeneric species. Analyses of gene expression from members of different castes and developmental stages revealed similarities in expression patterns among immature life stages. We also found evidence of DNA methylation within the genome of both species by directly analyzing DNA sequence reads. Moreover, genes that were highly and uniformly expressed were also relatively highly methylated. We further uncovered evidence of differences in patterns of molecular evolution in the two taxa, consistent with V. squamosa exhibiting alterations in evolutionary pressures associated with its facultatively parasitic or polygyne life history. Finally, rates of gene evolution were correlated with variation in gene expression between castes and developmental stages, as expected if more highly expressed genes were subject to stronger levels of selection. Overall, this study expands our understanding of how social behavior relates to genome evolution in insects.

Sequence Flow: interactive web application for visualizing partial order alignments

BackgroundMultiple sequence alignment (MSA) has proven extremely useful in computational biology, especially in inferring evolutionary relationships via phylogenetic analysis and providing insight into protein structure and function. An alternative to the standard MSA model is partial order alignment (POA), in which aligned sequences are represented as paths in a graph rather than rows in a matrix. While the POA model has proven useful in several applications (e.g. sequencing reads assembly and pangenome structure exploration), we lack efficient visualization tools that could highlight its advantages.ResultsWe propose Sequence Flow – a web application designed to address the above problem. Sequence Flow presents the POA as a Sankey diagram, a kind of graph visualisation typically used for graphs representing flowcharts. Sequence Flow enables interactive alignment exploration, including fragment selection, highlighting a selected group of sequences, modification of the position of graph nodes, structure simplification etc. After adjustment, the visualization can be saved as a high-quality graphic file. Thanks to the use of SanKEY.js – a JavaScript library for creating Sankey diagrams, designed specifically to visualize POAs, Sequence Flow provides satisfactory performance even with large alignments.ConclusionsWe provide Sankey diagram-based POA visualization tools for both end users (Sequence Flow) and bioinformatic software developers (SanKEY.js). Sequence Flow webservice is available at https://sequenceflow.mimuw.edu.pl/. The source code for SanKEY.js is available at https://github.com/Krzysiekzd/SanKEY.js and for Sequence Flow at https://github.com/Krzysiekzd/SequenceFlow.

Comprehensive Analysis of Complex Copy Number Variations in the Plasmodium falciparum Genome via Long-Read Sequencing: Implications for Antimalarial Resistance and Broad Genomic Insights

Abstract Introduction/Objective The genome of Plasmodium falciparum, characterized by a high AT content and a propensity for copy number variations (CNVs), presents a model for understanding how genomes adapt to stress. Our research focused on the prevalence and distribution of CNVs highlights long-read sequencing as a tool for detecting complex structural variants in other models. Methods/Case Report P. falciparum parasites were cultured and treated with DSM1 and aphidicolin that induce replication stress. Long-read sequencing was conducted using the Oxford Nanopore Technology Minion, with an emphasis on an optimizing DNA extraction method for preserving intact chromosomes. Visualization of CNVs, including tandem and inverted duplications, was achieved through a custom R Shiny application that allowed both quantitative and qualitative assessments. Results (if a Case Study enter NA) Long read sequencing identified an enrichment of multiple CNV types in treated samples, with a significant enrichment of inverted duplications as well as other types of structural variants. Aphidicolin and DSM1-treated samples demonstrated significant accumulation of these variations compared to controls (2.0 and 5.6-fold, respectively with a p value of &amp;lt; 0.0001 in both treatment groups). These occurred mainly in core genomic regions that do not contain variable gene copies. An increase in the number of inverted duplications may be related to stress response to the treatment and may be a potential adaptive mechanism of the parasite. Our manual single-read visualization technique made it possible to visualize the structure of various CNV types, which provides unique insight on CNVs as they arise in single genomes. Conclusion We hypothesize that inverted duplications represent stalled replication forks, and we are currently assessing their frequency in replicating and non-replicating cells. This study highlights the role of long-read sequencing in detecting complex genomic structures in P. falciparum, which may influence the parasite’s resistance to drug treatment. These findings not only advance our understanding of malaria pathology but also aid our comprehension of how cancer cells might react to drug-induced stress and the mechanisms behind CNV formation in neoplastic cells.

Direct genome sequencing of respiratory viruses from low viral load clinical specimens using the target capture sequencing technology.

The use of metagenomic next-generation sequencing technology to obtain complete viral genome sequences directly from clinical samples with low viral load remains challenging-especially in the case of respiratory viruses-due to the low copy number of viral versus host genomes. To overcome this limitation, target capture sequencing for the enrichment of specific genomes has been developed and applied for direct genome sequencing of viruses. However, as the efficiency of enrichment varies depending on the probes, the type of clinical sample, etc., validation is essential before target capture sequencing can be applied to clinical diagnostics. In this study, we evaluated the utility of target capture sequencing with a comprehensive viral probe panel for clinical respiratory specimens collected from patients diagnosed with SARS-CoV-2 or influenza type A. We focused on clinical specimens containing low copy numbers of viral genomes. Target capture sequencing yielded approximately 180- and 2,000-fold higher read counts of SARS-CoV-2 and influenza A virus, respectively, than metagenomic sequencing when the RNA extracted from specimens contained 59.3 copies/µL of SARS-CoV-2 or 625.1 copies/µL of influenza A virus. In addition, the target capture sequencing identified sequence reads in all SARS-CoV-2- or influenza type A-positive specimens with <26 RNA copies/µL, some of which also yielded >70% of the full-length genomes of SARS-CoV-2 or influenza A virus. Furthermore, the target capture sequencing using comprehensive probes identified co-infections with viruses other than SARS-CoV-2, suggesting that this approach will not only detect a wide range of viruses but also contribute to epidemiological studies.IMPORTANCETarget capture sequencing has been developed and applied for direct genome sequencing of viruses in clinical specimens to overcome the low detection sensitivity of metagenomic next-generation sequencing. In this study, we evaluated the utility of target capture sequencing with a comprehensive viral probe panel for clinical respiratory specimens collected from patients diagnosed with SARS-CoV-2 or influenza type A, focusing on clinical specimens containing low copy numbers of viral genomes. Our results showed that the target capture sequencing yielded dramatically higher read counts than metagenomic sequencing for both viruses. Furthermore, the target capture sequencing using comprehensive probes identified co-infections with other viruses, suggesting that this approach will not only detect a wide range of viruses but also contribute to epidemiological studies.

Metagenomic Analysis of Bacterial Communities in Bee Bread in Türkiye

This study aims to investigate the bacterial community structure in bee bread samples collected from 10 provinces of Türkiye using next-generation sequencing (NGS) and metagenomic analysis. Bacterial genomic DNA was extracted and sequenced using Illumina MiSeq. Bioinformatic analysis involved quality assessment, OTU classification, principal coordinate analysis (PCoA), and diversity index calculations. Heatmap and PCoA were utilized to explore the impact of locality and ecological zones on microbial diversity. Metagenomic analysis of 12 bee bread samples revealed 276,583 high-quality sequencing reads. The dominant bacterial phyla identified were Proteobacteria, Actinobacteria, Cyanobacteria, and Firmicutes. At the genus level, Streptomyces, Streptococcus, Bacillus, and Synechococcus were the most abundant, with Streptomyces and Bacillus playing key roles in the fermentation process of bee bread. The Shannon diversity index (ranging from 2.92 to 4.26) and Simpson's index (0.83 to 0.95) indicated high species diversity and relative abundance in bee bread. The study underscores the need for locality-specific approaches in beekeeping management and highlights the potential significance of beneficial bacterial taxa, particularly those involved in fermentation, in contributing to the nutritional and health properties of bee bread. These findings provide a foundation for further research on the microbial dynamics that support bee colony health.

Emergence and Rapid Diagnosis of Talaromyces marneffei Infections in Renal Transplant Recipients by Next-Generation Sequencing

In the last few years, next-generation sequencing (NGS) has emerged as a technology for laboratory diagnosis of many culture-negative infections and slow-growing microorganisms. In this study, we describe the use of metagenomic NGS (mNGS) for rapid diagnosis of T. marneffei infection in a 37-year-old renal transplant recipient who presented with chronic pneumonia syndrome. Bronchoalveolar lavage for mNGS was positive for T. marneffei sequence reads. Prolonged incubation of the bronchoalveolar lavage revealed T. marneffei colonies after 6 days of incubation. Analysis of 23 cases of T. marneffei infections in renal transplant recipients from the literature revealed that the number of cases ranged from 1 to 4 cases per five years from 1990 to 2020; but increased rapidly to 9 cases from 2021 to 2023, with 7 of them diagnosed by NGS. Twenty of the 23 cases were from T. marneffei-endemic areas [southern part of mainland China (n = 9); Hong Kong (n = 4); northeastern India (n = 2); Indonesia (n = 1) and Taiwan (n = 4)]. For the 3 patients from non-T. marneffei-endemic areas [United Kingdom (n = 2) and Australia (n = 1)], they had travel histories to China and Vietnam respectively. The time taken for diagnosis by mNGS [median 1 (range 1 to 2) day] was significantly shorter than that for fungal culture [median 6 (range 3 to 15) days] (P = 0.002). mNGS is useful for picking up more cases of T. marneffei infections in renal transplant recipients as well as providing a rapid diagnosis. Talaromycosis is an emerging fungal infection in renal transplant recipients.

Kallisto, bustools and kb-python for quantifying bulk, single-cell and single-nucleus RNA-seq.

The term 'RNA-seq' refers to a collection of assays based on sequencing experiments that involve quantifying RNA species from bulk tissue, single cells or single nuclei. The kallisto, bustools and kb-python programs are free, open-source software tools for performing this analysis that together can produce gene expression quantification from raw sequencing reads. The quantifications can be individualized for multiple cells, multiple samples or both. Additionally, these tools allow gene expression values to be classified as originating from nascent RNA species or mature RNA species, making this workflow amenable to both cell-based and nucleus-based assays. This protocol describes in detail how to use kallisto and bustools in conjunction with a wrapper, kb-python, to preprocess RNA-seq data. Execution of this protocol requires basic familiarity with a command line environment. With this protocol, quantification of a moderately sized RNA-seq dataset can be completed within minutes.

Transformer-based representation learning and multiple-instance learning for cancer diagnosis exclusively from raw sequencing fragments of bisulfite-treated plasmacell-free DNA.

Early cancer diagnosis from bisulfite-treated cell-free DNA (cfDNA) fragments requires tedious data analytical procedures. Here, we present a deep-learning-based approach for early cancer interception and diagnosis (DECIDIA) that can achieve accurate cancer diagnosis exclusively from bisulfite-treated cfDNA sequencing fragments. DECIDIA relies on transformer-based representation learning of DNA fragments and weakly supervised multiple-instance learning for classification. We systematically evaluate the performance of DECIDIA for cancer diagnosis and cancer type prediction on a curated dataset of 5389 samples that consist of colorectal cancer (CRC; n = 1574), hepatocellular cell carcinoma (HCC; n = 1181), lung cancer (n = 654), and non-cancer control (n = 1980). DECIDIA achieved an area under the receiver operating curve (AUROC) of 0.980 (95% CI, 0.976-0.984) in 10-fold cross-validation settings on the CRC dataset by differentiating cancer patients from cancer-free controls, outperforming benchmarked methods that are based on methylation intensities. Noticeably, DECIDIA achieved an AUROC of 0.910 (95% CI, 0.896-0.924) on the externally independent HCC testing set in distinguishing HCC patients from cancer-free controls, although there was no HCC data used in model development. In the settings of cancer-type classification, we observed that DECIDIA achieved a micro-average AUROC of 0.963 (95% CI, 0.960-0.966) and an overall accuracy of 82.8% (95% CI, 81.8-83.9). In addition, we distilled four sequence signatures from the raw sequencing reads that exhibited differential patterns in cancer versus control and among different cancer types. Our approach represents a new paradigm towards eliminating the tedious data analytical procedures for liquid biopsy that uses bisulfite-treated cfDNA methylome.

HyLight: Strain aware assembly of low coverage metagenomes

Different strains of identical species can vary substantially in terms of their spectrum of biomedically relevant phenotypes. Reconstructing the genomes of microbial communities at the level of their strains poses significant challenges, because sequencing errors can obscure strain-specific variants. Next-generation sequencing (NGS) reads are too short to resolve complex genomic regions. Third-generation sequencing (TGS) reads, although longer, are prone to higher error rates or substantially more expensive. Limiting TGS coverage to reduce costs compromises the accuracy of the assemblies. This explains why prior approaches agree on losses in strain awareness, accuracy, tendentially excessive costs, or combinations thereof. We introduce HyLight, a metagenome assembly approach that addresses these challenges by implementing the complementary strengths of TGS and NGS data. HyLight employs strain-resolved overlap graphs (OG) to accurately reconstruct individual strains within microbial communities. Our experiments demonstrate that HyLight produces strain-aware and contiguous assemblies at minimal error content, while significantly reducing costs because utilizing low-coverage TGS data. HyLight achieves an average improvement of 19.05% in preserving strain identity and demonstrates near-complete strain awareness across diverse datasets. In summary, HyLight offers considerable advances in metagenome assembly, insofar as it delivers significantly enhanced strain awareness, contiguity, and accuracy without the typical compromises observed in existing approaches.