Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “species differentiation”

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 19 recordsLinked to original sources
PubMed2025

DNABERT-S: pioneering species differentiation with species-aware DNA embeddings.

SUMMARY: We introduce DNABERT-S, a tailored genome model that develops species-aware embeddings to naturally cluster and segregate DNA sequences of different species in the embedding space. Differentiating species from genomic sequences (i.e. DNA and RNA) is vital yet challenging, since many real-world species remain uncharacterized, lacking known genomes for reference. Embedding-based methods are therefore used to differentiate species in an unsupervised manner. DNABERT-S builds upon a pre-trained genome foundation model named DNABERT-2. To encourage effective embeddings to error-prone long-read DNA sequences, we introduce Manifold Instance Mixup (MI-Mix), a contrastive objective that mixes the hidden representations of DNA sequences at randomly selected layers and trains the model to recognize and differentiate these mixed proportions at the output layer. We further enhance it with the proposed Curriculum Contrastive Learning (C2LR) strategy. Empirical results on 28 diverse datasets show DNABERT-S's effectiveness, especially in realistic label-scarce scenarios. For example, it identifies twice more species from a mixture of unlabeled genomic sequences, doubles the Adjusted Rand Index (ARI) in species clustering, and outperforms the top baseline's performance in 10-shot species classification with just a 2-shot training. AVAILABILITY AND IMPLEMENTATION: Model, codes, and data are publically available at https://github.com/MAGICS-LAB/DNABERT_S.

Sequence Analysis, DNA↗
PubMed2024

DNABERT-S: Pioneering Species Differentiation with Species-Aware DNA Embeddings.

We introduce DNABERT-S, a tailored genome model that develops species-aware embeddings to naturally cluster and segregate DNA sequences of different species in the embedding space. Differentiating species from genomic sequences (i.e., DNA and RNA) is vital yet challenging, since many real-world species remain uncharacterized, lacking known genomes for reference. Embedding-based methods are therefore used to differentiate species in an unsupervised manner. DNABERT-S builds upon a pre-trained genome foundation model named DNABERT-2. To encourage effective embeddings to error-prone long-read DNA sequences, we introduce Manifold Instance Mixup (MI-Mix), a contrastive objective that mixes the hidden representations of DNA sequences at randomly selected layers and trains the model to recognize and differentiate these mixed proportions at the output layer. We further enhance it with the proposed Curriculum Contrastive Learning (C2LR) strategy. Empirical results on 23 diverse datasets show DNABERT-S's effectiveness, especially in realistic label-scarce scenarios. For example, it identifies twice more species from a mixture of unlabeled genomic sequences, doubles the Adjusted Rand Index (ARI) in species clustering, and outperforms the top baseline's performance in 10-shot species classification with just a 2-shot training. Model, codes, and data is publicly available at https://github.com/MAGlCS-LAB/DNABERT_S.

Journal Article↗
PubMed2025

Genomic resequencing unravels species differentiation and polyploid origins in the aquatic plant genus Trapa.

Trapa L. is a non-cereal aquatic crop with significant economic and ecological value. However, debates over its classification have caused uncertainties in species differentiation and the mechanisms of polyploid speciation. This study employed whole-genome resequencing together with the fruit morphology of 229 Trapa accessions (153 Asian and 76 North American samples) to elucidate species differentiation and polyploidization events in Trapa. For the species with AA genome and large fruits, clear genetic differentiation was found between two clades with different geographic origins, that is, from the Yangtze River and Amur River basins. The invasive AA species in North America (T. natans) was identified as originating from the Amur River based on genetic and morphological similarities, while all the cultivated accessions were AA species originating from the Yangtze River with severe genetic impoverishment. The separation of the two BB species with small seeds, that is, T. incisa and T. maximowiczii, was strongly supported by both morphological and genetic evidence. For the tetraploids, Asian and North American tetraploids were revealed to have distinct evolutionary origins. Asian allotetraploids (AABB) originated through hybridization between AA diploids from the Yangtze River Basin and BB diploids T. maximowiczii, supported by nuclear and chloroplast evidence. In contrast, the invasive North American tetraploids (T. bispinosa var. iinumai) exhibited an AACC-like genome, suggesting an independent polyploidization involving an unknown "CC" diploid. These findings provide critical insights into Trapa's complex evolutionary history, polyploidizations, and invasive origins, offering a genomic foundation for the conservation and sustainable utilization of the underutilized aquatic crop amid global environmental challenges.

Polyploidy↗
PubMed2026

Revisiting the genome assembly of Lupinus species reveals differential diploidization after a shared whole-genome duplication.

Accurate genome assemblies are essential for comparative genomics, yet Hi-C-guided scaffolding can introduce structural errors that misrepresent chromosome architecture and bias evolutionary inferences. Here, we identified pervasive scaffolding errors-including artificial fusions, internal inversions, and incomplete contig mounting-in 2 previously published Lupinus genomes (L. cosentinii and L. digitatus) using a segmentation method based on long terminal repeat (LTR) retrotransposon density. We reassembled both genomes, producing chromosome-level references of 472.7 Mb (16 chromosomes) and 427.2 Mb (21 chromosomes), with BUSCO completeness >98.5%. Synteny validation and reapplication of LTR profiling confirmed that all prior errors were resolved. Using these corrected genomes together with 4 additional Lupinus species and 2 outgroup legumes, we investigated postpolyploid evolution. Synonymous substitution rate (Ks) analysis revealed a genus-specific whole-genome duplication (WGD) event (Ks = 0.17) shared by all 6 Lupinus species. The proportion of WGD-derived genes varied markedly, from 60% in L. digitatus to only 36% in L. mutabilis, indicating differential diploidization. While all species retained a core set of WGD duplicates enriched in cytoskeleton organization, ion transport, and defense responses, each exhibited lineage-specific functional trajectories: cell wall modification in L. cosentinii and L. digitatus, nitrogen metabolism in L. albus and L. angustifolius, flower development in L. luteus, and stress/lipid metabolism in L. mutabilis. Our corrected assemblies provide optimal references for Lupinus comparative genomics, and our findings demonstrate that a shared WGD event can lead to both conserved and highly divergent postpolyploid fates, likely underpinning adaptive diversification within the genus.

Lupinus↗
PubMed2025

Investigating the miRNA-mRNA interactome of human trabecular meshwork cells treated with TGF-β1 provides insights into the pathogenesis of pseudoexfoliation glaucoma.

Pseudoexfoliation glaucoma is a severe form of secondary open angle glaucoma and is associated with activation of the TGF-β pathway by TGF-β1. MicroRNAs (miRNAs) are small non-coding RNA species that are involved in regulation of mRNA expression and translation. To investigate what glaucomatous changes occur in the trabecular meshwork and how these changes may be regulated by miRNAs, we performed a bioinformatics analysis resulting in a miRNA-mRNA interactome. Primary human trabecular meshwork cells originating from normal donors were treated with TGF-β1 at 5 ng/mL for 24h; total RNA was extracted followed by RNA-Seq and miRNA-Seq. For both mRNA and miRNA species, differential expression was determined using a bioinformatics pipeline consisting of FastQC, STAR, FeatureCounts, edgeR (for miRNA) and DESeq2 (for mRNA). Putative mRNA-miRNA interactions between differentially expressed mRNA and miRNA species were determined using interaction databases miRWalk, miRTarBase, TarBase and TargetScan. To classify mRNA species by function and pathway, gene enrichment was performed using Enrichr. The resulting miRNA-mRNA interactome consisted of 1202 interactions. Some highly connected microRNAs were hsa-let-7e-5p, hsa-miR-20a-5p, hsa-miR-122-5p, and hsa-miR-29c-3p. Most differentially expressed genes were indicated to be regulated by miRNAs. The sub-interactomes of genes involved in specific pseudoexfoliation glaucoma related enrichment terms such as oxidative stress, unfolded protein response, signal molecules and ECM remodelling were determined. This is the first study to present a genome-wide microRNA-mRNA regulatory network for human trabecular meshwork cells treated with TGF-β1 and may serve to generate unbiased hypotheses about regulatory functions and mRNA targets of miRNAs in pseudoexfoliation glaucoma and may help to develop miRNA-based therapeutics.

Humans↗
PubMed2026

Integrated landscape of salivary metagenome and multi-biofluid metabolome characterizes a microbial-metabolic axis in upper gastrointestinal cancer progression.

BACKGROUND: Upper gastrointestinal cancer (UGIC) imposes a major global health burden, yet the stage-specific molecular changes along the microbial-metabolic axis remain limited understood. We aimed to delineate this molecular landscape across UGIC progression and evaluate its potential as non-invasive methods for precision screening. RESULTS: Derived from a multi-center population-based UGIC screening program, we enrolled 420 individuals, stratified into normal, low-grade intraepithelial neoplasia (LGIN), high-grade intraepithelial neoplasia (HGIN), and UGIC (n = 105 per group). Integrated salivary metagenomics and paired salivary/plasma metabolomics were performed to capture local and systemic dysregulation. We uncovered distinct stage-specific divergence during UGIC progression: profound remodeling of the salivary microbiota (104 differential species) and salivary metabolomics (80 differential metabolites) initiated early at the LGIN stage, whereas plasma metabolic dysregulation (40 differential metabolites) peaked significantly later at the HGIN stage. Integrative analysis revealed salivary microbiota related more closely with salivary metabolome than plasma metabolome. Moreover, statistical evidence suggested that dysbiotic salivary microbiota was associated with altered lysine- and tryptophan-related catabolic pathways converging on Acetyl-CoA-related metabolic nodes, supporting a potential metabolic mechanism in precancerous lesions. Finally, the discriminative model integrating metagenomic and metabolomic markers demonstrated promising diagnostic performance in distinguishing these precancerous lesions (LGIN: area under the curve [AUC] = 0.83; HGIN: AUC = 0.77) and UGIC (AUC = 0.76) from normal. CONCLUSION: This study characterizes a stage-specific microbial-metabolic axis that facilitates the comprehensive understanding of UGIC pathogenesis. These multi-biofluid signatures offer a promising non-invasive triage strategy for detecting precancerous lesions and optimizing endoscopic resource allocation. Video Abstract.

Female↗
PubMed2025

Leclercia barmai sp. nov., isolated from worm castings of Eisenia fetida, is a urease-positive, 3-nitropropionic acid and glycerol-consuming bacterium.

A comprehensive polyphasic characterization has validated the unique taxonomic position of a novel bacterium, strain EMC7T, isolated from the worm castings of earthworm, Eisenia fetida, collected from the Centre for Floriculture and Agri-Business Management (COFAM), NBU (26.7072° N, 88.3554° E). Whole-genome sequence of this Gram-stain-negative, facultatively anaerobic, motile, rod-shaped bacterium showed maximum sequence homology with Leclercia adecarboxylata NBRC 102595T, placing it within the genus Leclercia. The genome of EMC7T is 5.03 Mbp with a G + C content of 56.3 mol%. Phylogenetic analyses established its distinctiveness from Leclercia adecarboxylata and Leclercia tamurae. DNA-DNA hybridization (dDDH) value was 23.6%, and the average nucleotide identity (ANI) was 82.1%, both below the thresholds for prokaryotic species differentiation. Predominant fatty acids were C16:0 (29.53%), summed feature 3 (C16:1ω7c/C16:1ω6c, 16.51%), and C18:1ω7c (10.90%). Notably, EMC7T exhibited urease activity and could metabolize 3-nitropropionic acid (3-NPA), glycerol, tellurite, selenate, and selenite, suggesting potential bioremediation applications. Biochemical tests, phenotypic traits, genotypic data, and physiological properties cumulatively differentiated EMC7T from its closest relatives. Based on chemotaxonomic, phenotypic, genomic, and phylogenetic evidence, strain EMC7T represents a novel bacterial species of the genus Leclercia, for which the name Leclercia barmai sp. nov. (type strain EMC7T = MCC 5183T = JCM 36544T) is proposed.

Animals↗
PubMed2026

Sex- and development-specific transcriptomic profiling of venom and silk genes in the wolf spider Pardosa astrigera provides insights into ecological adaptation and predatory strategies.

Spider venom and silk glands are two major secretory systems that contribute to prey capture, defense, and reproduction, but their sex- and development-specific molecular regulation in wandering wolf spiders remains poorly understood. Here, the transcriptome of Pardosa astrigera, an important agricultural natural enemy in China, revealed significant sex- and development-associated molecular differentiation among adult females, adult males, and spiderlings. A total of 100,025 unigenes were obtained, of which 23,852 were functionally annotated, providing a comprehensive transcriptomic resource for this species. Differential expression patterns showed marked variation among groups, with 531, 1792, and 832 DEGs detected in PAF vs PAS, PAM vs PAS, and PAF vs PAM, respectively. These genes were mainly associated with metabolic, oxidation-reduction, cuticle development, MAPK signaling, and lysosome pathways. Fifteen co-expression modules revealed distinct expression patterns. The turquoise, pink, yellow, and red modules were development-related, whereas the blue module was male-biased. Venom- and spidroin-related genes were distributed across multiple modules, suggesting coordinated regulation. Overall, 42 venom peptides, 21 venom proteins, and 11 spidroins were identified. Representative genes showed strongly biased expression, including spiderling-biased U3_Pp1a and U5_Pp1e, female-biased U4_Pp1a, and male-biased SMase D_108750 and PaTuSp_108466. These findings reveal sex- and development-biased expression patterns of venom- and silk-related candidate genes in P. astrigera and may provide molecular insights into ecological adaptation and predatory strategies in wandering wolf spiders.

Animals↗
PubMed2026

An Amplicon Panel for High-Throughput and Low-Cost Genotyping of Yesso Scallop Mizuhopecten yessoensis.

The Yesso scallop Mizuhopecten yessoensis was imported from Japan to western Canada in the late 1980s to establish an economically viable scallop aquaculture industry. Since this time, the industry in Canada has operated with existing genetic diversity within the broodstock, which is considerably limited relative to wild populations. The sector has not been able to realise its full potential in part due to idiopathic hatchery failures and farm stock collapses due to disease outbreaks associated with the intracellular bacterial pathogen Francisella halioticida. To support Yesso scallop production and breeding, here we generate a low-density, genotyping-by-sequencing amplicon panel using single nucleotide polymorphism (SNP) markers that are evenly spaced across the M. yessoensis genome and that show high heterozygosity in Canada and Japan. The panel can also exploit the high genetic polymorphism of the M. yessoensis genome, with de novo SNP calling identifying over 2,500 high quality SNPs within the 579 sequenced amplicons. We demonstrate the utility and versatility of this new genotyping tool for breeding applications including parentage assignment, low density family-based genome-wide association study, trait heritability evaluation to determine potential for genomic selection, and species differentiation (against the weathervane scallop Patinopecten caurinus). We did not find any genomic regions significantly associated with F. halioticida resistance but did identify potential for genomic selection. We could separate the two species based on genotypes, and did not see evidence of a past M. yessoensis x P. caurinus hybridization event within the M. yessoensis breeding population at Vancouver Island University. This low-cost genotyping panel is expected to accelerate selective breeding improvements for M. yessoensis in Canada and elsewhere.

Animals↗
PubMed2026

Development of a Multiplex Polymerase Chain Reaction Assay for Differentiating Three Lactococcus Species Associated With Piscine Lactococcosis.

Piscine lactococcosis is an important bacterial disease of farmed fish. The causative agents, Lactococcus garvieae, Lactococcus petauri and Lactococcus formosensis, are closely related, which complicates species-level identification. We developed a conventional multiplex PCR assay targeting species-specific genes identified by comparative genomic analysis. Average nucleotide identity reassignment of 441 publicly available genome assemblies identified 111 L. garvieae, 255 L. petauri and 75 L. formosensis genomes. Species-specific primers and a tuf-based Lactococcus common control were evaluated using in silico polymerase chain reaction (PCR) against target genomes and 10,460 off-target assemblies representing 474 taxa in 12 genera. Experimental specificity was assessed using six target strains and 19 non-target fish pathogens. Distinct amplicons of 195, 333 and 500 bp were produced for L. garvieae, L. petauri and L. formosensis, respectively, together with a 132-bp control amplicon. No cross-amplification was observed. All target species were detected in mixed-DNA samples and spiked kidney and spleen homogenates from two fish species. Analytical detection limits were estimated based on microscopic cell counts of bacterial suspensions before DNA extraction and ranged from 0.956 to 8.55 cell equivalents per reaction. This assay represents a rapid, low-cost method for differentiating lactococcosis-causing Lactococcus species using standard PCR and agarose gel electrophoresis.

Lactococcus garvieae↗
PubMed2025

Strong but diffuse genetic divergence underlies differentiation in an incipient species of marine stickleback.

Understanding how lineages proceed along the "speciation continuum" and how species boundaries are maintained over time remain central questions in evolutionary biology. Populations early in the speciation process can give us detailed insight into the reproductive barriers that first initiate speciation. In this study, we explore the nature of genomic divergence between two sympatric marine stickleback ecotypes from Atlantic Canada, "whites" and "commons". Males of each ecotype exhibit distinct nuptial colorations, nesting habits, and parental care strategies. Using population genomic analyses of SNPs and copy number variants (CNVs; deletions and duplications) we show that whites and commons consistently form distinct populations. We uncover genomic differentiation in the white ecotype characteristic of an incipient species, showing extremely low genome-wide differentiation (FST) and very recent divergence (~1 kya). Demographic analysis detected very low levels of ongoing gene flow between populations. Our results and prior genomic studies suggest that reproductive isolation is being maintained between ecotypes despite recent evidence that hybridization in nature does occur. Contrary to other systems, we found many small, but dispersed regions of high differentiation throughout the genome rather than explicitly within chromosomal inversions or the sex chromosomes. On chromosomes VII and XVI, we identified CNVs overlapping genes enriched for olfaction, which may play a role in differences in reproductive strategies between ecotypes. Ultimately, our results demonstrate that genome-wide rather than localized differences can underlie the early stages of divergence, and that this pattern is corroborated by both SNPs and CNVs.

Copy Number Variation↗
PubMed2026

Brucellar spondylitis is associated with disturbance in gut microbiota and histamine metabolism associated inflammation.

BACKGROUND: The pathogenesis of brucellar spondylitis (BLS) has traditionally been considered to be primarily limited to local osteoarticular lesions. With the proposal of the "gut-spine axis" concept, the role of intestinal microecological dysbiosis in inflammatory spinal diseases has attracted in an increase of attention. The overactivated inflammatory cytokine network not only mediates bone destruction and intervertebral disc damage, but also forms a bidirectional interaction with gut microbiota dysbiosis through the "gut-spine axis," collectively driving disease progression. However, the inflammatory mechanism by which gut microbiota participates in the pathological process of BLS remains largely unclear. METHODS: This study recruited 20 BLS patients and 20 healthy donors. Multi-omics analysis including metagenomics, untargeted metabolomics, and targeted short-chain fatty acids (SCFAs) analysis, were used to compare the structural differences in gut microbiota between the two groups and screen for signature differential bacterial species. Plasma levels of histamine and histidine decarboxylase were measured by ELISA to clarify the role of differential histidine metabolic pathway in the disease. Additionally, plasma levels of lipopolysaccharide (LPS) and inflammatory cytokines (IL-1β, IL-6, IL-10, IL-17A, TNF-α) were detected by ELISA. The correlation between gut microbiota and inflammatory indicators was further analyzed. RESULTS: Compared to the healthy control group, the α-diversity of the gut microbiota in BLS patients was significantly reduced, with the microbial community structure exhibiting increased homogeneity. Beta diversity analysis revealed significant differences, suggesting that disease progression is associated with an overall imbalance in the gut microbiota and the deterioration of its specific structural composition. At the phylum level, the abundances of Actinomycetota, unclassified_d_Viruses, and Fusobacteriota were significantly increased in the gut microbiota of BLS patients compared to the control group, while the abundances of Bacillota and Pseudomonadota were significantly decreased. Further analysis revealed that, compared to the control group, the generic abundance of Enterococcus was significantly increased, while the proportions of Blautia, Faecalibacterium, Ruminococcus, Agathobacter, Roseburia, Clostridium, Eubacterium, Alistipes and Anaerobutyricum were significantly decreased. At the species level, the abundances of Enterococcus sp and Enterococcus-faecium were increased, whereas Blautia sp, Ruminococcus sp, Faecalibacterium sp, Faecalibacterium prausnitzii, Agathobacter rectalis, Eubacterium sp, Agathobacter sp, and Roseburia sp were decreased. Furthermore, untargeted metabolomics revealed that metabolites were enriched in the histidine metabolic pathway, and the levels of SCFAs including butyrate, isobutyrate, valerate, and 4-methylvalerate in the intestinal contents were reduced in BLS. Functional KEGG profiling revealed that key KOs involved in butyrate synthesis (e.g., K00074, K00172, K01640) and transport were globally downregulated in the patient group, whereas histidine decarboxylase KOs (K01693, K11755, K19787) that convert histidine to pro-inflammatory histamine were significantly enriched. The loss of butyrate-producing symbionts led to SCFAs deficiency and mucosal barrier disruption, creating ecological niches for facultatively anaerobic Enterococcus, which further exacerbated local inflammation via proteolytic fermentation and histamine production. Compared with the control group, BLS patients showed decreased plasma levels of IL-10, while levels of IL-1β, IL-6, IL-17A, and TNF-α were increased, and LPS levels were elevated. In addition, significantly elevated plasma pro-inflammatory LPS levels in patients with BLS suggest disruption of intestinal integrity and permeability. Correlation analysis indicated a close relationship between gut microbiota and inflammation. CONCLUSION: BLS is associated with gut microbiota dysbiosis and alterations in microbial metabolites, which may be linked to inflammatory responses and histamine metabolism. The differential microbial taxa identified in this study could be developed into a stool-based non-invasive diagnostic panel to facilitate early differentiation of BLS from other spinal disorders. Furthermore, restoring gut microbial balance through probiotic supplementation or dietary modulation may represent a promising adjunctive strategy to enhance the efficacy of standard antibiotic therapy and reduce disease recurrence.

Humans↗
PubMed2026

Genome-wide association analysis reveals specialization to hosts and niches in multiple species of the Lactobacillaceae.

The Lactobacillaceae inhabit diverse environments, but the extent of their habitat adaptation remains unclear and the colonization factors unknown. First, we applied multiple machine learning models to determine if we can distinguish strains of the same species isolated from two different habitats based on their gene content. Surprisingly, we show that no species is differentially adapted to the oral cavity versus the human gut, or food versus the human gut, while only Lactobacillus crispatus showed specialization to the human urogenital system versus human gut. We then asked which species of Lactobacillaceae are habitat-specialized and how they could be identified. Using multiple lifestyle predictors incorporated in logistic regression models, we found that Limosilactobacillus reuteri, Ligilactobacillus ruminis, L. salivarius, L. crispatus, and L. mucosae displayed the highest degrees of host specialization. Applying our microbial genome-wide association study tool, aurora, to these species identified genes encoding adhesins and bacteriocins as the strongest and most common adaptation factors. This work establishes a generalizable framework for identifying novel species-habitat pairs with strong evidence of specialization and for uncovering the genomic features underlying within-species host and habitat adaptation.

Humans↗
PubMed2025

Evolutionary legacy of the "living fossil" genus Parrotia (Hamamelidaceae): genomic insights into species divergence and polygenic adaptation.

Despite their long evolutionary history, the genomic basis of adaptation and speciation in "living fossil" plants remain largely unexplored. Parrotia, a Tertiary relict tree genus with two extant species, P. subaequalis and P. persica, exhibits a disjunct distribution between East Asia and West Asia. Here, we present the first chromosome-level assemblies for both species, confirmed their sibling relationship, and dated the speciation event to the early Miocene. The recent proliferation of long-terminal repeat retrotransposons has driven the genome expansion in P. subaequalis. We detected widespread heterogeneous genomic differentiation between species. Extensive signals of divergent selection, local adaptation, and elevated Ka/Ks ratios in Parrotia indicate that this genus has undergone adaptive evolution in distinct refugia, challenging the notion of it as an "evolutionary dead end". Our findings provide new insights into the genomic evolution, environmental adaptation, and speciation of this "living fossil" tree genus.

Genome, Plant↗
PubMed2026

Emergence of Babesia naoakii infection in Indonesian domestic cattle, a new host record in water buffaloes, and characterization of complete mitochondrial protein-coding genes.

Babesia (B.) naoakii, previously referred to as Babesia sp. Mymensingh, is a recently characterized tick-borne haemoprotozoan parasite of cattle. In Indonesia, we first reported its presence in 2022 from clinically affected cattle in Central Java. To investigate the wider epidemiology of this neglected ruminant-associated Babesia species, we surveyed apparently healthy cattle (Bos indicus) and water buffaloes (Bubalus bubalis) across three districts of Java, Indonesia. A PCR assay targeting the B. naoakii-specific apical membrane antigen 1 (ama1) gene detected the parasite occurrence in 34.39% of assessed cattle (87/253; 95% CI: 28.80-40.44%) and 30.77% of water buffaloes (12/39; 95% CI: 18.47-46.52%). These results represent the first record of B. naoakii infection in water buffaloes in the country and confirm that the parasite circulates in subclinically infected bovine hosts. To characterise this apicomplexan parasite further at the molecular level, we assembled in full length the three mitochondrial protein-coding genes (PCGs): cytochrome c oxidase subunits 1 (cox1) and 3 (cox3), as well as cytochrome b (cytb). These genes were reconstructed by next-generation sequencing of blood DNA collected during the acute haemolytic-phase of B. naoakii infection, from calves that subsequently succumbed to the disease in the endemic area. Phylogenetic analyses of the concatenated amino-acid sequences of cox1, cox3, and cytb placed the Indonesian isolates within a well-supported monophyletic clade, distinct from all previously characterised ruminant-associated Babesia species and sister to the Babesia bigemina/Babesia ovata lineage. This placement confirmed species identity and reinforced the genetic distinctiveness of B. naoakii in Indonesia. Notably, although B. naoakii circulates in peripheral blood and mirrors the diagnostic behaviour of the mild pathogen B. bigemina, its clinical impact more closely resembles that of the severe pathogenic B. bovis, particularly in young animals. This diagnostic-clinical discordance highlights the need for B. naoakii-specific molecular surveillance and species-level differentiation in regions of co-endemicity. Given the high prevalence in subclinically B. naoakii-infected adults, the documented severity of babesiosis in calves, and the potential for substantial economic losses, broader epidemiological investigations and species-specific control measures for B. naoakii are urgently performed. The same holds true for future epizootiological investigations of underdiagnosed B. naoakii-infections possibly circulating in Indonesian endemic ruminant bovids such as the banteng (Bos javanicus), the lowland anoa (Bubalus depressicornis) and the tamaraw (Bubalus mindorensis).

Animals↗
PubMed2026

Oxford Nanopore Sequencing of Clinical DNA for Identification and Comparative Genomic Analysis of Erysipelothrix piscisicarius.

The genus Erysipelothrix comprises facultative anaerobic, nonspore-forming, gram-positive bacteria that can cause skin infections and severe diseases such as septicemia and endocarditis in humans. Although E. rhusiopathiae is the primary pathogen, other species may also be involved, necessitating accurate identification. However, 16S rDNA sequencing lacks sufficient resolution to differentiate among Erysipelothrix species. In this study, we used Oxford Nanopore Technology (ONT) to directly sequence low-quality DNA extracted from heart valve tissue of a 66-year-old female patient with a fatal case of septicemia and aortic endocarditis. In contrast to 16S rDNA Illumina sequencing and matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF MS), which incorrectly identified the pathogen as E. rhusiopathiae, direct sequencing via ONT precisely identified E. piscisicarius as the cause of infection. About 1.47 Mb genome was retrieved from nanopore direct sequencing. Within the E. piscisicarius genome, we detected genes associated with virulence. Phylogenetic analysis showed that our strain clustered with a human-derived E. piscisicarius strain from China and swine-derived strains from Brazil. In conclusion, this study demonstrated that ONT can be used to sequence low-quality DNA extracted directly from patient specimens, obtain a draft bacterial genome, and reliably distinguish between pathogenic species.

Aged↗
PubMed2026

Metagenomic profiling of gut microbiome in post-cholecystectomy patients with diarrhea: a nested case-control study.

BACKGROUND: Cholecystectomy can cause diarrhea, with an incidence as high as 57.2%, seriously impacting patient prognosis. To investigate the gut dysbiosis following cholecystectomy and identify microbial biomarkers and functional genomics associated with post-cholecystectomy diarrhea (PCD), we conducted a nested case-control study within a prospective cohort. METHODS: We enrolled a cohort of 160 patients. At follow-up completion, 30 patients who developed PCD were matched with 30 non-PCD (NPCD) controls. 16 S rRNA sequencing was used to analyze gut microbiota structure and diversity (mainly at genus level). Representative fecal samples underwent metagenomic sequencing for species level and genetic differential analysis. RESULTS: The potentially pathogenic bacterial species Coprococcus comes and Blautia sp. were significantly enriched in the gut microbiota of PCD patients, with their abundance positively correlated with the degree of intestinal inflammation. In contrast, the potentially beneficial bacterial species Bacteroides intestinalis and Prevotella copri, known to contribute to lipid metabolism and play a role in modulating gut immunity and suppressing inflammatory responses, were found to be significantly depleted in PCD patients. Further metagenomic functional analysis revealed significant enrichment of pathways related to cell motility, membrane transport, and sulfur metabolism in PCD patients. CONCLUSIONS: This work identified potential beneficial and pathogenic bacterial species associated with the onset of PCD, as well as significantly enriched functional pathways within the intestinal microbiota. These findings provide a scientific basis for elucidating the relationship between PCD and gut microbiota, and provide candidate microbial signatures and functional pathways that may inform future microbiota-targeted strategies, pending external and mechanistic validation.

Humans↗
PubMed2025

Comparative and systems analyses of Leishmania spp. non-coding RNAs through developmental stages.

Leishmania spp. is the etiological agent of leishmaniases, neglected diseases that seek to be eradicated in the coming years. The life cycle of these parasites involves different host and stress environments. In recent years, many studies have shown that several protein-coding genes are directly involved with the development and host interactions. However, little is still known about the role of non-coding RNAs (ncRNAs) in life cycle progression. In this study, we aimed to identify the genomic structure and function of ncRNAs from Leishmania spp. and to get insights into the repertoire of ncRNAs (RNAome) of this protozoan genus. We studied 26 strains corresponding to 16 different species of Leishmania. Our RNAome analysis revealed the presence of several ncRNAs that are shared among different species, allowing us to differentiate between subgenera as well as between species that are canonically related to visceral leishmaniasis. We also studied co-expression relationships between coding genes and ncRNAs which in the amastigote developmental stage for Leishmania braziliensis and Leishmania donovani revealed the presence of miRNA-like transcripts co-expressed with several coding genes involved in starvation, survival and histone modification. This work represents the first effort to characterize the Leishmania ssp. RNAome, supporting further approaches to better understand the role of ncRNAs in gene regulation, infective process, and host-parasite interaction.

RNA, Untranslated↗