SEARCH · Search PubMed

Results for “single-cell RNA sequencing”

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

SourceResearch connectionFrom yearThrough year

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

374 records · Page 4Linked to original sources

PubMed2026

Dynamics of antibiotic resistance genes co-occurrence with pathogenic and non-pathogenic bacteria throughout wastewater treatment processes.

Wastewater treatment plants (WWTPs) are recognized hotspots for antibiotic resistance genes (ARGs) and pathogenic bacteria. Despite advancements in treatment technologies, the persistence of ARGs and pathogenic bacteria remains a concern. In this study, we analyzed the dynamic changes in ARGs and bacterial communities throughout the treatment processes within an anaerobic-anoxic-oxic (AAO) WWTP over one week by using HT-qPCR coupled with 16S rRNA gene amplicon sequencing. The connectedness index, based on network analysis, showed that the dynamics of ARGs and mobile genetic elements (MGEs) were more strongly associated with potentially pathogenic bacteria than with non-pathogenic bacteria, suggesting that ARG immigration and dissemination in the WWTP were likely driven by potentially pathogenic taxa. The AAO treatment significantly reduced ARGs in final effluent (EF) (∼64 %) and residual sludge (RS) (∼81 %); however, potential hosts of ARGs such as Comamonas testosteroni and Clostridioides difficile persisted with minimal changes in relative abundance and remained detectable in EF and RS. Notably, the abundance of ARGs was lower in RS than in EF, and source tracking analysis identified influent as the primary source of ARGs and potentially pathogenic taxa in EF, underscoring the greater health risks associated with effluent discharge.

Wastewater↗

PubMed2026

Integrative multi-omics reveals a fibroblast-centered, ZFHX3-prioritized regulatory framework linking sick sinus syndrome and atrial fibrillation.

OBJECTIVE: To define shared genetic and multi-scale mechanisms underlying comorbidity between sick sinus syndrome (SSS) and atrial fibrillation (AF). METHODS: We integrated genome-wide association study (GWAS) summary statistics for SSS and AF with Genotype-Tissue Expression (GTEx) expression and splicing quantitative trait loci (eQTL/sQTL), atrial single-cell and spatial transcriptomics, and epigenomics. We identified trait-relevant tissues and pathways, prioritized shared cell types, quantified genome-wide and local genetic sharing, detected joint loci by cross-trait meta-analysis, and linked loci to regulatory programs via colocalization and cell-prioritized co-expression networks. RESULTS: Both traits showed strongest enrichment in cardiac tissue, especially Heart Atrial Appendage. Fibroblasts from the left atrial appendage were consistently prioritized as the key shared cell population. SSS and AF displayed significant positive genome-wide genetic correlation, with multiple locally shared regions, including six major loci. Cross-trait meta-analysis identified eight joint-phenotype SNPs implicating four susceptibility genes. ZFHX3 was the leading tissue-cell-gene candidate, acting as a hub in fibroblast co-expression modules and colocalizing with cardiac regulatory signals. CONCLUSION: Shared liability for SSS and AF is highly tissue- and cell-specific, converging on regulatory networks in atrial appendage fibroblasts, with ZFHX3 serving as a central mechanistic and biomarker node.

Humans↗

PubMed2026

Genome-wide SNP data support species boundaries in sympatric Polylepis Ruiz & Pav. (Rosaceae) species from Bolivia and Ecuador.

Species delimitation in the South American genus Polylepis is notoriously challenging due to high morphological similarity and phenotypic plasticity, likely driven by hybridization and gene flow. Previous phylogenetic studies suggested that genetic structure aligns more strongly with geography than with taxonomy, questioning existing species concepts and hampering conservation efforts. We used double-digest RAD sequencing (ddRADseq) to generate genome-wide SNP data for 11 Polylepis species sampled across multiple localities in Bolivia and Ecuador. Population genetic analyses, phylogenetic inference, and network approaches were combined to assess whether genetic structure aligns more closely with taxonomy or geography. Morphologically defined species formed largely cohesive genetic lineages across regions, with species identity explaining substantially more genetic variation than locality. While localized admixture and reticulation were detected among closely related taxa, widespread species showed strong genetic cohesion and clear separation from congeners. Our results indicate that the sampled Polylepis species from Bolivia and Ecuador maintain distinct genetic identities despite localized signals consistent with gene flow. This genome-wide support for current taxonomy highlights Polylepis as a valuable model for studying speciation under gene flow and indicates that multiple geographic sampling will be essential in reconstructing a robust phylogeny of the genus, with important implications for conservation planning in Andean montane forests.

Bolivia↗

PubMed2026

Responding to a protracted tuberculosis outbreak: lessons from multiple rounds of investigation in a Chinese boarding school.

PURPOSE: This study analysed a multi-semester pulmonary tuberculosis (PTB) cluster outbreak in a Chinese boarding school to provide evidence for future epidemic control. METHODS: Contacts were screened via symptoms, infection tests and chest radiography. Screening expanded progressively from close contacts to same-floor contacts, then all students and staff. Whole-genome sequencing (WGS) with single nucleotide polymorphism (SNP) and bioinformatics analysis was used for lineage classification, transmission clustering (≤12 SNPs defining a cluster) and drug resistance prediction. RESULTS: From 2020 to 2022, 20 students were diagnosed with PTB, half laboratory-confirmed. Most cases clustered in class 16 and were epidemiologically linked to the primary case (case 0), who had household PTB exposure. Case 0 and case 1 had diagnostic delays exceeding 3 and 6 months, respectively. WGS of five isolates (case 1, 3, 4, 9 and 10) collected over three semesters showed all belonged to lineage 2 and differed by ≤12 SNPs, confirming the same transmission chain. The infection rate in class 16 (46.34%) was significantly higher than other case classes (19.05%) and classes without cases (8.27%) (χ2 = 61.169, p < 0.001). No new cases were detected during a one-year follow-up of students involved in the outbreak after the final round of screening, nor among household contacts of all cases followed up to the present. CONCLUSIONS: Lack of entry health examinations facilitated the outbreak. Delayed diagnosis, incomplete contact screening and absence of preventive treatment led to cross-semester persistence. The infection rate disparity confirms class 16 as the outbreak epicentre. Improving community case management, extending contact follow-up and enhancing cluster outbreak measures are recommended to prevent future outbreaks.

Humans↗

PubMed2026

Cytonuclear conflict and reticulate evolution in the Morelloid clade (Solanum, Solanaceae): Insights from genome skimming and network Phylogenomics.

The Morelloid clade (black nightshades) is one of the most strongly supported clades within the megadiverse Solanum genus. It comprises 76 globally distributed, non-spiny herbaceous and suffrutescent species. While often erroneously considered poisonous weeds, several species are economically important as orphan crops. The clade is closely related to tomato and potato but, due to a lack of focused breeding efforts, remains a putative reservoir of genetic diversity for crop improvement. Despite this potential, we lack fundamental knowledge on the evolution of the Morelloid clade. The group includes polyploid species with unknown parental origins-likely reflecting reticulate processes such as hybridization, introgression, and associated backcrossing events. Prior analyses have been unable to disentangle these processes, leaving the mechanisms underlying reticulate evolution in the Morelloid clade poorly understood. Here, we use genome skimming to produce a well-supported maximum likelihood plastid phylogeny from complete circularized plastomes and a coalescent-based species tree from combined Angiosperms353 and conserved ortholog set nuclear markers. Our dataset, composed of previously published data and deep genome skimming from herbarium samples, spans 26 Morelloid species. To investigate phylogenetic discordance, we used a nuclear phylogenetic network, multispecies coalescent simulations, a fused rooted nuclear chloroplast tree, and quantification of nuclear gene tree concordance. We show that incongruence between nuclear and plastid trees is pervasive and cannot be explained by incomplete lineage sorting alone. Instead, our results demonstrate that events consistent with repeated chloroplast capture have shaped the reticulate evolutionary history of the clade, especially among African polyploid and Pan-American diploid lineages.

Phylogeny↗

PubMed2026

Enrichment of Lysobacter in a long-term organically managed agricultural field with low soilborne disease incidence.

Disease-suppressive soils, in which soilborne pathogens are naturally suppressed, offer a promising model for sustainable crop protection, particularly in organic farming systems where chemical disease control options are limited. Although disease suppression in these soils is considered to rely on biological control, the underlying mechanisms remain poorly understood. In this study, we investigated soil from a long-term organically managed field in Shiga Prefecture, Japan, where soilborne disease incidence has remained consistently low, to identify bacterial community features potentially associated with this field. The 16S rRNA gene amplicon sequencing indicated that this soil harbored a bacterial community distinct from those of nearby agricultural soils. Following the application of organic compounds, the genus Lysobacter, a taxon with known antagonistic activity against plant pathogens, was markedly enriched in response to proteinaceous organic inputs. This enrichment was consistent across sampling times and specific to certain proteinaceous organic inputs, whereas minimal effects were observed on chitin, N-acetyl-d-glucosamine, or cysteine. Broader soil surveys indicated that Lysobacter enrichment was not strictly associated with whether soils had been managed under organic or conventional farming practices. Stepwise multiple regression analysis identified 10 co-occurring bacterial genera that were strongly associated with Lysobacter abundance. These findings highlight condition-dependent Lysobacter enrichment as a characteristic microbial response to proteinaceous organic amendments in this low-disease-incidence field and provide microbial insights that may inform microbiome-based strategies for sustainable soil management.

Lysobacter↗

PubMed2026

Nanopore-based epigenomic profiling reveals the absence of widespread CpG methylation in the African swine fever virus genome.

DNA methylation is a critical epigenetic mechanism implicated in regulating replication and transcription in DNA viruses. However, the epigenetic landscape of African swine fever virus (ASFV), a large double-stranded DNA virus infecting pigs, remains controversial. Here, we systematically profiled the DNA methylome of the first ASFV strain isolated in Hong Kong (HK_NT_202103) using Oxford Nanopore Technologies (ONT) R10.4.1 sequencing. We employed a paired design: native whole-genome sequencing (WGS) against a methylation-free whole-genome amplification (WGA) control. Using conservative thresholds, we found no evidence of 5-methylcytosine (5mC), especially typical CpG methylation, across the viral genome. Importantly, clear CpG methylation signals were successfully detected in the host genome from WGS data, confirming the functionality of the workflow to detect 5mC at CG sites. While widespread 5mC seems absent, a small number of putative N6-methyladenine (6mA) loci were identified. A specific 6mA candidate exhibited raw ionic current disruptions and gene-level intersection with another ASFV isolate (CAS19-01/2019), although it lacked single-base consensus across different methylation callers or between the two isolates. Although our biological findings are restricted to a single isolate under specific experimental conditions, this study introduces a novel, highly rigorous ONT framework for viral epigenomics research. Furthermore, the absence of ASFV CpG methylation indicates that host CpG-depletion remains a viable strategy for viral metagenomic enrichment. Ultimately, our work offers a critical methodological baseline for ASFV surveillance and highlights the necessity of targeted experimental validation for rare viral modifications.

African Swine Fever Virus↗

PubMed2026

Sea urchin co-culture boosts abalone growth by reducing environmental stress and remodeling gut microbiota.

Biofouling and microenvironmental deterioration are major bottlenecks restricting the intensive aquaculture of Pacific abalone (Haliotis discus hannai). While co-culturing offers an eco-friendly mitigation strategy, the underlying mechanisms promoting abalone growth remain poorly understood. This study evaluated the growth performance of H. d. hannai co-cultured with varying densities of the sea urchin (Strongylocentrotus intermedius). By employing transcriptome and 16S rRNA sequencing of the abalone gut, we investigated the synergistic responses of host gene expression and gut microbiota. Compared with the monoculture group, the co-culture groups showed significantly less biofouling and greater growth of abalone, with the co-culture (n = 15) exhibiting the best outcomes. Transcriptomic analysis revealed 1444, 760, and 508 DEGs in G5, G10, and G15, respectively, compared with G0. These DEGs were significantly enriched in metabolic pathways, including glycolysis and sterol metabolism, indicating a shift in intestinal energy metabolism from stress defense toward growth under co-culture conditions. Gut microbiota profiling identified Proteobacteria and Firmicutes as the dominant phyla, with specific functional taxa (e.g., Psychrilyobacter and Akkermansia) enriched in a density-dependent manner. Furthermore, correlation analysis demonstrated that growth traits positively correlated with growth-promoting taxa (e.g., the unclassified AB1 lineage), but negatively correlated with potentially opportunistic taxa (e.g., Tabrizicola). These findings provide insights into a potential synergistic mechanism of "environmental stress alleviation-metabolic reprogramming-microecological remodeling" driving abalone growth, providing a theoretical foundation for optimizing co-culture systems and developing growth-associated biomarkers.

Animals↗

PubMed2025

New insights into soil amendment: Impact of humic acid on typical antibiotic resistance in agricultural soil.

Humic acid (HA) addition can improve agricultural soil, but little is known about how it affects the soil resistome. In this study, we used selective agar plate combined with quantitative PCR (qPCR) and 16S rRNA gene sequencing to investigate how HA influences antibiotic resistant bacteria (ARB) and antibiotic resistant genes (ARGs) in soil contaminated with erythromycin and kanamycin. 0.1 % HA reduced the abundance of culturable erythromycin-resistant bacteria (ERB), while promoting the growth of kanamycin-resistant bacteria (KRB). Lysinibacillus and Paenibacillus were the dominant genera in ERB and KRB, respectively, governing the changes in their abundances. At this concentration, the Lysinibacillus abundance in ERB decreased from 96.74 % to 70.57 %. Meanwhile, that of Paenibacillus in KRB increased from 33.40 % to 77.44 %. The copy number of ermF decreased after HA addition, while that of ermB increased. Furthermore, 0.1 % HA significantly reduced the copy number and relative abundance of aadA1 and aac(6')-Ib (aka aacA4)-03 in the soil. Changes in these two types of ARB and ARGs were primarily driven by shifts in the microbial community structure. Soil physicochemical properties, particularly increased organic matter (OM), altered the absolute abundance of ermB. Meanwhile, changes in intI1 abundance determined the risk associated with aadA1 and aac(6')-Ib (aka aacA4)-03. These findings emphasize the dual role of HA in the dissemination of antibiotic resistance in agricultural soils and highlight the necessity of considering dose-dependent effects when applying HA as a soil amendment.

Soil Microbiology↗

PubMed2026

Whole-exome characterization of host genetic variation in HIV-associated genes across the high-prevalence Mizo population, Northeast India.

BACKGROUND: The Mizoram state of Northeast India has one of the highest HIV prevalence rates in Asia, yet the host genetic factors influencing HIV susceptibility in this Tibeto-Burman population remain uncharacterised. METHODS: We performed whole-exome sequencing using Illumina NovaSeq 6000, mean coverage 100X on 76 HIV-negative Mizo individuals. Variants were called using GATK HaplotypeCaller v4.3 against GRCh38p14, annotated with ANNOVAR, and filtered using hard-quality thresholds (QD ≥ 2, SOR ≤ 3, MQ ≥ 40, DP ≥ 10, GQ ≥ 20). The allele frequencies were compared against gnomAD v2.1.1 population databases. Hardy-Weinberg equilibrium was assessed using the Wigginton exact test with Bonferroni correction. RESULTS: Post-quality filtering resulted in 12,011 sample-variants across 2,821 unique positions from 36 HIV-associated loci (33 protein-coding genes, 2 chemokine ligands, and 3 lncRNA targets). Of these, 784 observations (51 unique positions) were high-impact nonsynonymous or loss-of-function variants. ADAR rs2229857 (p.K384R, NM_015840) was the most frequently observed variant (Mizo carrier frequency = 0.895; 95% CI: 0.806-0.946). CXCR1 rs16858808 (p.R335C) showed the greatest population enrichment (Mizo carrier frequency = 0.197; 95% CI: 0.123-0.300; 7.65-fold carrier-frequency enrichment versus gnomAD South Asian; CADD = 15.60). Sixteen of 20 tested variants deviated from Hardy-Weinberg equilibrium after Bonferroni correction (p < 0.0025), predominantly showing excess homozygosity consistent with the endogamous Mizo population. The protective variant CCR5-Δ32 was absent in all the 76 individuals tested. CONCLUSION: This first whole-exome characterization of HIV host genes in the Mizo population identifies CXCR1 rs16858808 as the most population-enriched functional variant and reveals a pervasive endogamy signature. These findings provide a population-specific genetic framework for future HIV susceptibility studies and ART pharmacogenomics research.

Humans↗

PubMed2026

Can't see the forest for the trees: The influence of marker type on inferred phylogenetic relationships in a cosmopolitan bat genus.

Fine-resolution information on species relationships and biological diversity is critically needed to guide conservation efforts amidst rapid environmental changes. Systematics, which forms the foundation of this knowledge, has been revolutionized by phylogenomics, utilizing genome-scale datasets. However, the use of diverse marker types, non-comparable taxon sampling, and outgroup selection can lead to conflicting phylogenetic hypotheses. These inconsistencies complicate study comparisons and hinder our ability to assess marker-specific impacts on phylogenetic resolution. The phylogenetic reconstruction of the bat genus Myotis, encompassing over 140 species and characterized by a rapid radiation in the last 20 million years, has been particularly influenced by these challenges. Achieving phylogenetic resolution in Myotis is particularly complex due to subtle interspecific differences in both morphological and molecular traits. Mitochondrial and nuclear markers often produce discordant trees, influenced by hybridization, introgression, and methodological variations. In this study, we employed a consistent taxonomic sample set of 44 Myotis taxa to evaluate the impact of five different genetic marker types on phylogenetic reconstruction. We observed significant discordance between topologies derived from conserved nuclear and mitochondrial markers and found that transposable elements were inadequate for resolving relationships across the entire genus. Our results also clarify the placement of previously problematic taxa within the genus. These findings emphasize the importance of aligning genetic marker choice with specific phylogenetic questions and highlight the influence of taxonomic and methodological variation on phylogenomic outcomes. This work provides a framework for improving phylogenetic inference in rapidly radiating groups and enhances our understanding of evolutionary history in Myotis.

Animals↗

PubMed2026

Transcriptomic insights into thermal stress reveal physiological trade-off between thermal stress adaptation and reproductive investment in Spodoptera litura.

Spodoptera litura, a highly polyphagous lepidopteran pest, poses a major threat to agricultural productivity due to its remarkable adaptability to diverse environmental conditions. Although heat stress is known to trigger transcriptional reprogramming in insects, the molecular mechanisms underlying thermal stress responses in S. litura remain poorly understood. In the present study, fourth-instar larvae were exposed to acute heat stress (44 °C) and compared with control conditions (27 ± 1 °C) to investigate heat-induced transcriptional alterations affecting physiology and reproduction. High-quality RNA-Seq data achieved more than 80% mapping efficiency, with a total of 15,782 transcripts were identified. Transcriptome analysis of S. litura larvae showed 323 differentially expressed genes (DEGs), of which 262 genes were significantly upregulated and 61 were downregulated in heat-stressed larvae compared to the control group. The DEGs were associated with stress response, reproduction, signalling, proteostasis, detoxification, oxidative stress, metabolism, development, and chromatin regulation. Heat shock proteins genes, including HSP70, HSP90, and HSP27, together with co-chaperones such as TRET-1, STIP1, and Starvin, were strongly upregulated, indicating enhanced cellular protection against protein damage and oxidative stress under heat stress. Conversely, key reproductive and cell cycle-related genes, including BARR, CAPD2, FEO, CDK2 and MORULA, were significantly downregulated, suggesting reproductive impairment and developmental arrest. RT-qPCR validation corroborated the RNA-Seq findings, demonstrating a heat-induced physiological trade-off that prioritizes survival over reproduction. Consistent with these molecular responses, heat-stressed insects exhibited marked reproductive impairment, including significant reductions in gonadosomatic index, eupyrene sperm bundle count, mating frequency, mating success, female calling behaviour, copulation duration, fecundity, and egg fertility. Collectively, these findings provide comprehensive insights into the molecular basis of thermal adaptation in S. litura and demonstrate that acute heat stress compromises reproductive fitness while activating conserved stress-response pathways that promote short-term survival.

Animals↗

PubMed2026

Simultaneously PYCR-1 and ALH-6 inhibition exacerbates 6-PPD quinone toxicity via disrupting proline and glutamate metabolisms and activating insulin signals in Caenorhabditis elegans.

Glutamate synthesized from the proline can serve as a precursor for key intermediate metabolites of citric acid cycle. Recently, we observed reduced glutamate content and expression of alh-6 controlling glutamate synthesis by 6-PPD quinone (6-PPDQ) in Caenorhabditis elegans. However, possible effect of 6-PPDQ on proline synthesis and the association with 6-PPDQ toxicity induction remain unclear. After 0.1-10 μg/L 6-PPDQ exposure, proline content was further reduced, and expression of pycr-1 governing proline biosynthesis was decreased. In 6-PPDQ exposed nematodes, RNA interference (RNAi) of pycr-1 decreased α-ketoglutarate content, enhanced mitochondrial dysfunction, reduced nicotinamide adenine dinucleotide (NADH) and reduced flavine adenine dinucleotide (FADH₂) contents, inhibited mitochondrial complex I/II activities, and decreased expressions of gas-1 and mev-1. Moreover, compared to single RNAi, double RNAi of pycr-1 and alh-6 exacerbated the 6-PPDQ toxicity in reducing α-ketoglutarate, NADH, and FADH₂ contents, and suppressing mitochondrial complex I/II activities and gas-1 and mev-1 expressions. Additionally, double RNAi of pycr-1 and alh-6 intensified toxicity of 6-PPDQ on longevity and caused upregulation of insulin ligand and receptor genes and downregulation of daf-16 and its targeted genes in 6-PPDQ exposed nematodes. Furthermore, after 6-PPDQ exposure, daf-16 RNAi suppressed pycr-1 and alh-6 expressions, suggesting formation of a regulatory feedback loop between pycr-1/alh-6 and daf-16. Our findings highlight involvement of disrupted proline and glutamate metabolisms in 6-PPDQ-induced mitochondrial dysfunction and reduced longevity.

Animals↗

PubMed2026

Comparison of paralog identification methods and their impact on species tree topologies in target capture phylogenomics within the Sindora clade (Detarioideae: Leguminosae).

Target capture is a common method of generating high throughput DNA sequencing data for phylogenetic reconstruction of species relationships, for which single copy genes are usually most informative. However, a pervasive problem with target capture is that putatively single copy genes may in fact be paralogs resulting from gene duplication, which are problematic for phylogenetic inference because their evolutionary history may differ from the divergence history of species. Here, we use as a case study a target enrichment dataset of 88 species of Detarioideae (Leguminosae) with a focus on the Sindora clade to examine approaches for handling paralogs, including the built-in paralog handling functions in HybPiper and CAPTUS, plus subsequent steps using Putative Paralog Detection and the tree-based Yang & Smith orthology inference approach. We compare the paralogs flagged using these methods and verify their performance with BLAST mapping against a reference genome sequence of Sindora glabra, and then subsequently compare the species tree topologies produced across these methods. Our comparisons of paralogs flagged across the Sindora clade show that the Putative Paralog Detection pipeline was the most accurate in identifying paralogs in terms of its similarity to the BLAST mapping, followed by the built-in paralog identification function of CAPTUS. However, the results we recovered for the Detarioideae subfamily suggest that the largest differences in species tree topology resulted from the use of paralog-filtered alignments (such as with the Putative Paralog Detection pipeline and the Yang & Smith orthology inference approaches) rather than just by removing the sequences of identified paralogous genes. This was the true for HybPiper-assembled datasets but was not seen in CAPTUS-assembled datasets. In all comparisons, the topological differences caused by different paralog handling methods tended to be confined to clades where processes such as hybridisation and introgression are prevalent. Our study provides a roadmap to establish the best approach to identify, eliminate or separate paralogs in the absence of a chromosomally contiguous reference genome for a study group, and highlights the importance of careful data inspection and processing in addition to understanding the extent of paralogy and paralog characteristics (e.g. sequence divergence between copies) for their study group.

Phylogeny↗

PubMed2026

Antennal transcriptome analysis of chemosensory proteins in the raspberry weevil, Aegorhinus superciliosus (Coleoptera: Curculionidae).

Aegorhinus superciliosus (Coleoptera: Curculionidae) is a polyphagous pest of economic importance in southern Chile, the chemical ecology of which remains poorly characterized. Across insect species, chemosensory proteins, including odorant receptors (ORs), gustatory receptors (GRs), ionotropic receptors (IRs), odorant-binding proteins (OBPs), chemosensory proteins (CSPs), and sensory neuron membrane proteins (SNMPs), mediate the detection of chemical cues involved in host selection, reproduction, and other ecologically relevant behaviors. In this study, the antennal transcriptome of adult A. superciliosus was sequenced and analyzed using a de novo RNA-seq approach. Three independent biological replicates per sex were used for RNA-seq, and the same number of independent biological replicates was used for RT-qPCR validation; sequencing yielded 147,409,936 high-quality reads after quality filtering. A total of 112 candidate chemosensory genes were identified, comprising 43 ORs, 34 OBPs, 10 CSPs, 18 IRs, 5 GRs, and 2 SNMPs. Phylogenetic analyses assigned these candidate proteins to established clades, providing a comparative framework for functional inference for ORs and OBPs. Sex- and tissue-biased expression analyses revealed that several ORs, including AsupOR4, AsupOR19, and AsupOBP13, exhibit antennal enrichment and sex-specific expression patterns. Notably, AsupOR19 and AsupOBP13 displayed strong female-biased expression. In addition, transcripts of selected ORs and OBPs were detected in non-antennal tissues, such as the rostrum and legs, suggesting potential functional versatility beyond canonical olfaction. Together, these findings represent the first molecular identification of the chemosensory repertoire of A. superciliosus. This study establishes a foundation for reverse chemical ecology approaches aimed at identifying behaviorally active volatile organic compounds (VOCs) toward environmentally sustainable strategies for integrated pest management.

Animals↗

PubMed2026

First identification and molecular subtyping of Blastocystis spp. in donkeys in Aksaray province, Türkiye.

Blastocystis is a common intestinal protist worldwide that can infect humans and animals. Although its molecular epidemiology in Türkiye is mostly focused primarily on humans and livestock, equids have received limited attention despite their traditional roles and frequent contact with humans and other animals in rural environments. This study aimed to determine the molecular prevalence and subtype (ST) distribution of Blastocystis spp. in donkeys in Aksaray Province, providing the first molecular data on donkeys in Türkiye. A total of 182 fresh fecal samples were collected from donkeys in nine villages within Aksaray province. Genomic DNA was extracted, and the small subunit ribosomal RNA (SSU rRNA) gene fragment of Blastocystis spp. was amplified via PCR analysis. Positive isolates were sequenced bidirectionally for identification and subsequent phylogenetic analysis of Blastocystis in donkeys. The overall molecular prevalence of Blastocystis spp. in donkeys was 4.4% (8/182). The infection rate was higher in young donkeys (under 3 years old; 8.33%) than in adults (3 years or older; 2.46%). However, this difference was not statistically significant. Sequence analysis of the positive PCR products revealed the presence of one known livestock-specific subtype, ST10. Phylogenetic analysis showed that the ST10 isolates characterized in this study clustered with isolates identified from different hosts. This study provides the first molecular data on Blastocystis presence in donkeys in Türkiye. The exclusive detection of ST10 suggests potential cross-species transmission, likely facilitated by the traditional practice of co-housing donkeys with other animals in confined barns. These findings indicate that donkeys may contribute to Blastocystis transmission, underscoring the importance of a "One Health" approach in future epidemiological surveillance.

Animals↗

PubMed2026

Uce-based phylogeny and classification of Megachilini.

The generic-level classification of the bee tribe Megachilini (Megachilidae) has remained controversial due to poor phylogenetic resolution at the base of the group, particularly among the brood parasitic genera and the numerous dauber ("Chalicodoma s. l.") lineages. We present a phylogenomic analysis of Megachilini based on ultraconserved elements (UCEs), sampling 52 ingroup taxa with emphasis on the dauber lineages. We also present a combined UCE + six-gene analysis to improve taxon coverage, resulting in a dataset with 127 ingroup taxa. Maximum likelihood, coalescent, and Bayesian analyses of multiple UCE matrices recover largely congruent topologies with substantially improved support relative to previous studies. Our results strongly support the monophyly of Megachilini, the early divergence of Noteriades and Gronoceras, and a single origin of brood parasitism. All remaining non-parasitic Megachilini form a moderately supported clade sister to the brood parasitic lineage. The leafcutter bees are monophyletic and nested within dauber lineages. Several major dauber clades are consistently recovered, including an exclusively Australian clade corresponding to the Hackeriapis group of subgenera, while several recognized subgenera are paraphyletic. The lineage known as Morphella, previously placed in synonymy with the subgenus Callomegachile, was not closely related to that subgenus and is here treated as a valid subgenus. Divergence-time analyses place the crown age of Megachilini in the late Eocene to early Oligocene, with major extant lineages diversifying during the Miocene. Limited morphological diagnosability of several clades indicates that splitting non-parasitic lineages into numerous genera would result in an impractical classification that would widen the gap between taxonomists and non-specialists and exacerbate the taxonomic impediment in bees. We therefore advocate retaining a single genus Megachile for non-parasitic Megachilini (excluding Noteriades and Gronoceras), as the classification best supported by phylogenomic evidence and most robust to future taxon sampling.

Animals↗

PubMed2026

Microbial allies in a cotton pest: A descriptive account of associated microbiota dynamics in Dysdercus cingulatus across development.

BACKGROUND: Hemipteran insects harbour several symbiotic partners, mainly bacteria, which play pivotal roles for hosts like dietary provision, support overall physiology, xenobiotic degradation and manipulate/regulate behaviour. Most of these symbionts usually reside and operate from the digestive tracts of the animals. Cotton is one of the major cash crops in India and Dysdercus cingulatus (D. cingulatus) though a secondary pest, is causing significant destruction of cotton bolls, poor lint quality and reduce oil content of seeds. Premature opening of cotton bolls often leads to bacterial and fungal infections, thus resulting in extensive economic loss worldwide. D. cingulatus is a hemimetabolous insect that comprises of developmental stages like egg, nymph (5 instar stages), and adult. The present work explored the ontogeny specific diversity in the associated microbiota and predicted their probable functional inputs in D. cingulatus. RESULTS: The data obtained using 16S rRNA gene sequencing (NovaSeq 6000) revealed presence of members of Proteobacteria (65.83%), Firmicutes (24%), Actinobacteria (10%) phyla throughout the ontogeny of D. cingulatus. Highest alpha diversity of these symbiotic bacteria was recorded in the third instar nymphs in contrast to rest of the developmental stages. Among all the observed genera, Stenotrophomonas, Hungatella and Glutamicibacter were predominant from egg to adult stages. MicFunPred, a tool used for predicting the probable functional inputs of these symbionts, hinted at their probable stage specific contribution in crucial biochemical pathways such as polyketide biosynthesis, ascorbate/aldarate metabolism, pentose phosphate and glyoxylate cycles, steroid hormone and peptidoglycan biosynthesis, and glycolysis/pyruvate metabolism. CONCLUSIONS: The primary investigations on the ontogenetic composition and diversity of associated microbiota, suggest dynamic shifts in D. cingulatus, concurrent with their probable functions/roles in the host development and metabolism. To the best of our knowledge, this is the first report on symbiotic microbiota variation across the developmental stages of D. cingulatus that provides preliminary descriptive observations that may guide future functional and experimental investigations into microbiota-based pest management.

Animals↗