Search PubMedSearch

SEARCH · Search PubMed

Results for “Transcriptome sequencing”

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

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

121 records · Page 3Linked to original sources

Unveiling the molecular basis of gonadal development: Multi-omics uncovers sex-related genes and steroid pathways in Sinonovacula constricta.

The razor clam Sinonovacula constricta is an economically important cultured mollusk in China, but the molecular mechanism of its gonadal development and sexual differentiation remains unclear. This study integrated gonadal transcriptomic, proteomic, and metabolomic analysis to identify key sex-related molecules. Transcriptome analysis identified 2795 DELs and 6497 DEGs between sexes, including the sex-related genes Fem-1b, Fem-1c, GUCY1B2 and FAT4, as well as a regulatory network of 39 lncRNA-mRNA pairs involving Tektin-4, Ropporin-1, Histone H1, and FoxN4. Proteomic analysis revealed 3217 DEPs: Tektin family members, Ropporin-1 and Tssk proteins were upregulated in the testis, while histone H1 and FAT4 were upregulated in the ovary. Metabolomic analysis detected 409 DEMs, with uridine identified as a potential sex differential marker (upregulated in the ovary), and 23 gonadal development-related DEMs showed sex-specific upregulation. Integrative transcriptome-proteome analysis identified 1543 co-expressed DEGs/DEPs enriched in nucleosome assembly, oxidative phosphorylation, and carbon metabolism, including key sex-related genes AKAP14, Tektin/Tssk families, Histone H1, and FAT4. Transcriptome-metabolome integration identified 32 shared KEGG pathways (e.g., biosynthesis of unsaturated fatty acids, pyrimidine metabolism), while proteome-metabolome integration revealed 5 (positive ion) and 6 (negative ion) co-enriched pathways, with alanine, aspartate and glutamate metabolism and oxidative phosphorylation being functionally relevant to gonadal development. Collectively, these results reveal the molecular basis of gonadal development, highlight critical sex-related genes and steroid metabolic pathways, and provide valuable resources for future reproduction and breeding in S. constricta.

Animals

Genome-wide analysis of the plant-specific PLATZ gene family in Taraxacum kok-saghyz and its roles in response to drought and salt tolerance.

Abiotic stress severely limits plant growth and productivity. Taraxacum kok-saghyz Rodin (TKS), known for its environmental resilience, represents a valuable resource for identifying stress-tolerant genes to improve stress-adaptive crops. Plant AT-rich protein and zinc-binding protein (PLATZ) transcription factors serve as core regulators of plant growth, developmental processes, and adaptive responses to various stress conditions; however, they remain uncharacterized in TKS. Here, we identified 10 TksPLATZ genes through a whole-genome analysis. Phylogenetically, these genes were grouped into five distinct evolutionary branches. Promoter sequence analysis revealed multiple types of cis-acting regulatory elements that are connected with hormonal signal responses and environmental stress adaptation. Integrated analysis of transcriptome datasets and RT-qPCR validation demonstrated that TksPLATZ genes display tissue-specific expression profiles and show distinct responsive patterns to drought and salt stress treatments. Among them, TksPLATZ1, TksPLATZ2 and TksPLATZ7 were markedly induced under both stressors and were selected for further functional study. We demonstrated that TksPLATZ1, TksPLATZ2 and TksPLATZ7 localize to the cell nucleus and act as transcriptional activators and repressors, respectively. Phenotypic data from overexpression experiments in plants confirm that heterologous expression of TksPLATZ1, TksPLATZ2, and TksPLATZ7 enhances the tolerance of Arabidopsis to salt and osmotic stress. These findings provide valuable genetic resources for improving plant tolerance to environmental stresses.

Salt Tolerance

Distinct cell morphotypes of Aureobasidium melanogenum ZN exhibit differential functional profiles in promoting maize growth.

Black yeast-like fungi of the genus Aureobasidium exhibit morphological plasticity, but whether distinct cellular states within the same genetic background are associated with different plant growth-promoting functions remains unclear. Here, yeast-like cells (YL), swollen cells (SC), and chlamydospores (CH) of Aureobasidium melanogenum ZN were characterized. YL was associated mainly with siderophore production and laccase activity, SC with extracellular polysaccharide accumulation, and CH with phosphate mobilization and higher ammonia and IAA production. Whole-genome and comparative genomic analyses revealed a shared repertoire related to nutrient acquisition, auxin-associated metabolism, extracellular oxidation, and carbohydrate remodeling, with expansions in nutrient- and cell-surface-related gene families. Transcriptomic and metabolomic analyses showed distinct deployment of these capacities, with CH exhibiting broad reprogramming of tryptophan-associated, nitrogen, phosphate, central-carbon, and amino-acid metabolism. In maize, CH at the optimal inoculation concentration of 105 CFU·mL-1 produced the strongest growth promotion, increasing plant height, dry biomass, root length, root surface area, and root volume by 58.6%, 365.1%, 191.0%, 194.3%, and 222.4%, respectively. Consistent with this pronounced growth phenotype, maize root transcriptomics showed coordinated CH-induced responses involving root development, nutrient transport, redox regulation, and root-interface remodeling. Root-zone tracking showed greater short-term stability and persistence of CH. These findings identify cellular state as an important functional dimension of Aureobasidium-plant interactions and provide a basis for developing fungal inoculants with defined beneficial cellular states.

Zea mays

Integrated analysis uncovers exogenous induction and molecular regulation of erinacine A accumulation in Hericium erinaceus.

Erinacine A, a cyathane-type diterpenoid mainly from Hericium erinaceus mycelia, exhibits prominent neurotrophic and neuroprotective activities, making it a promising candidate for managing neurodegenerative diseases. However, its low abundance and unclear genetic regulatory mechanisms hinder its application as a nutraceutical. This study aimed to decipher its regulatory mechanisms and enhance production. Four exogenous inducers were screened, with salicylic acid (SA) and ergosterol (ERG) significantly increasing erinacine A content by 62.21% and 146.70% at 20 days, respectively. Transcriptome and WGCNA of inducer-treated sample identified darkorange and magenta modules associated with erinacine A biosynthesis, with the eri gene cluster enriched in the darkorange module and eriG and eriF as hub genes. Forward genetic analysis via QTL mapping of the HeD127 dikaryon population revealed significant phenotypic variation in erinacine A content (0.341-13.085 mg/g) and identified two loci (erA-1 and erA-2) explaining 18.63% of phenotypic variation. Integrating these forward and reverse genetic analyses revealed that salicylic acid and ergosterol synergistically regulate core carbon metabolic pathways to augment acetyl-CoA supply for the mevalonate pathway, suppressed competitive metabolism, enhanced diterpene skeleton construction and structural modification. These results deepen our understanding of the genetic and molecular basis governing accumulation of erinacine A, and facilitate its application in neuroprotective pharmaceuticals.

Diterpenes

Long-term heat exposure reshapes muscle molecular regulation and enhances thermal tolerance in Clarias fuscus.

Rapid fluctuations in water temperature driven by global warming have become a major abiotic stressor affecting muscle function in teleost fish. This study examined the effects of long-term thermal conditions on heat tolerance in Clarias fuscus. Fish were maintained for 90 days at either a normal temperature group (NT, 26 °C) or a high-temperature group (HT, 34 °C). Subsequently, muscle histology, and transcriptomic profiles were observed following acute high-temperature exposure (34 °C) and after temperature recovery (26 °C). Histological analysis showed that fish from the NT under acute high-temperature stress exhibited severe muscle damage (atrophy, myofilament disruption, and myolysis), whereas fish from the HT displayed markedly reduced lesions. RNA-seq profiling revealed 5769 differentially expressed genes (DEGs) in the NT and 3292 DEGs in the HT following acute temperature challenges. Functional enrichment indicated that, in the HT, modulation of key cell cycle regulators (e.g., ccna, ccnb, cdk1, cdk2) contributed to alleviating muscle damage caused by temperature fluctuations. In the NT, genes associated with ribosome biogenesis (e.g., nop56, riok2, riok1) were up-regulated and then down-regulated during temperature fluctuation, whereas p53 in the cell cycle pathway showed the opposite expression pattern. These findings demonstrate long-term heat exposure reshapes molecular expression and regulatory mechanisms in the muscle of C. fuscus, thereby enhancing thermal tolerance and adaptability, and providing a theoretical basis for breeding heat-resistant, high-quality aquaculture strains.

Animals

Physiological and molecular responses of coelomocytes to low- to mid-frequency acoustic exposure in the sea urchin Strongylocentrotus intermedius.

Underwater noise is a widespread environmental pollutant in marine ecosystems, yet the effects of low- to mid-frequency acoustic exposure on immune physiology and molecular responses in echinoderms remain unclear. In this study, the sea urchin Strongylocentrotus intermedius was exposed to continuous pure-tone acoustic stimulation at 80, 125, 250, 500, 750, and 1000 Hz for 3 h. Results showed that acoustic exposure significantly affected redox homeostasis, energy metabolism, and immune function in S. intermedius coelomocytes: the antioxidant system and glutathione redox balance were altered, as indicated by increased superoxide dismutase (SOD) and catalase (CAT) activities, elevated reduced glutathione (GSH) content, and a higher GSH/GSSG ratio; glycolysis-related enzyme activities were enhanced, with increased pyruvate kinase (PK) activity under 125-500 Hz exposure and elevated hexokinase (HK) activity at 250 Hz; and immune function was impaired, as shown by increased coelomocyte mortality, reduced phagocytic activity, and inhibited acid phosphatase (ACP) and alkaline phosphatase (AKP) activities, whereas respiratory burst activity showed no significant change. Among all treatments, 250 Hz induced the most pronounced physiological responses. Transcriptomic analysis of coelomocytes from the 250 Hz group identified 663 differentially expressed genes, including 537 upregulated and 126 downregulated genes, mainly enriched in pathways related to apoptosis, phagosome, lysosome, glutathione metabolism, arachidonic acid metabolism, and carbohydrate metabolism. These findings indicate that low- to mid-frequency acoustic exposure can act as a physiological and molecular stressor to S. intermedius coelomocytes by affecting redox homeostasis, enhancing energy metabolism, and suppressing immune effector processes, with 250 Hz showing the strongest effect under the present exposure conditions. This study provides experimental evidence for evaluating the potential biological effects of low- to mid-frequency acoustic exposure on benthic echinoderms.

Animals

Correlative analysis of endogenous miRNA expression profiles underlying brown planthopper adaptation to resistant rice.

The brown planthopper (Nilaparvata lugens Stål, BPH) is a major insect pest threatening global rice production. However, the molecular mechanisms underlying the adaptation of BPH populations with different virulence levels to resistant rice cultivars remain poorly understood. MicroRNAs (miRNAs), as key post-transcriptional regulators, play critical roles in host adaptation in herbivorous insects. In this study, we analyzed the miRNA expression profiles of a high-virulent population (IR56p) and a low-virulence population (TN1p) after feeding on susceptible (TN1) and resistant (IR56) rice cultivars. Our findings reveal distinct miRNA-mediated regulatory strategies employed by the two populations. The IR56p population showed downregulation of miRNAs including miR-10, miR-124, and miR-316, showing an inverse correlation with increased expression of predicted target genes involved in detoxification (carboxylesterase, UDP-glycosyltransferase) and effector function (calmodulin). In contrast, several miRNAs highly expressed in IR56p, including miR-307, miR-317, and miR-275, were predicted to target rice genes associated with hormone signaling, cell wall biosynthesis, and oxidative homeostasis, suggesting a possible but unproven inter-species regulatory role that requires functional validation. Collectively, these descriptive and correlative findings provide hypothesis generating insights into insect-plant coevolution and identifies candidate molecular targets for future functional validation and RNA interference-based pest management strategies.

Animals

Diagnostic and clinical utility of exome sequencing and chromosomal microarray in children with GDD/iD: a meta-analysis.

BACKGROUND: Global developmental delay/intellectual disability (GDD/ID) is among the most common neurodevelopmental disorders, with up to half of cases are attributed to genetic factors. Chromosome microarray (CMA) has traditionally been the primary genetic test for idiopathic GDD/ID. However, whole exome sequencing (WES) and whole genome sequencing (WGS) have recently emerged, substantially increasing diagnostic yields in these populations. METHODS: We conducted a comprehensive literature search of PubMed, Scopus, EMBASE, and the Cochrane Library from inception to April 29, 2025. Studies reporting the diagnostic utility of these tests in children with GDD/ID were included and analyzed. RESULTS: A total of 102 studies, comprising 55,752 children, were reviewed. The pooled diagnostic yield of WES was 0.37 (95% CI: 0.33-0.41; I2 = 93%), significantly higher than that of CMA at 0.19 (95% CI: 0.16-0.21; I2 = 95%). Subgroup analyses showed that WES yielded significantly higher diagnostic rates than CMA in both same-sample comparisons (OR = 2.27, 95% CI: 1.08-4.78) and different-sample comparisons (OR = 1.65, 95% CI: 1.15-2.37). Only one study evaluated WGS, reporting a diagnostic yield of 0.27. Meta-regression revealed a significant association between CMA diagnostic yield and the proportion of male participants (p&#x2009;<&#x2009;0.01), but not with WES. No significant difference in diagnostic utility was observed between isolated GDD/ID and GDD/ID with comorbidities. CONCLUSION: In children with unexplained GDD/ID, WES demonstrates superior diagnostic and clinical utility compared to CMA. Incorporating WES as a first-line investigation in the diagnostic evaluation of GDD/ID may be warranted.

Humans

Whole genome sequencing of unusual Hepatitis C virus subtypes and drug resistance analysis during direct-acting antiviral therapy in India.

INTRODUCTION AND OBJECTIVES: Pangenotypic direct-acting antivirals (DAA) are effective against highly prevalent Hepatitis C virus (HCV) subtypes, but have been clinically validated almost exclusively in high-income countries. Unusual HCV subtypes may carry natural polymorphisms, potentially impacting DAA susceptibility. We conducted full-genome characterization and resistance analysis of unusual HCV subtypes in patients receiving DAA treatment. PATIENTS AND METHODS: In this prospective hospital-based study, eligible patients were screened for anti-HCV antibodies and active infection was confirmed by diagnostic 5'NCR-based HCV RNA detection. Genotyping was performed by core region sequencing, and viral load quantified by real-time PCR. For whole genome sequencing, multiplex primers were designed using alignments of global reference sequences. Sequencing was carried out using the Oxford Nanopore Technology platform. Phylogenetic analysis used multiple sequence alignment and the HCV-GLUE resource for resistance-associated substitution (RAS) analysis. RESULTS: Predominant genotype was genotype 3 in 64.3% (n = 45); genotype 6 in 21.4% (n = 15); and genotype 1 in 14.2% (n = 10). Unusual HCV subtype 6xa was detected in two patients and showed no NS5A resistance mutations. One genotype 3b patient relapsed at 24 weeks post-DAA treatment completion and carried NS5A resistance-associated substitutions 30 K and 31 M both at baseline and at relapse, conferring high-level resistance to NS5A inhibitors. CONCLUSION: This is the first report from India of whole genome sequencing of HCV subtype 6xa. The identification of NS5A resistance mutations in the 3b relapse case underscores challenges for global HCV elimination strategies.

Humans

Evaluation of one-step amplicon-based targeted enrichment for SARS-CoV-2 whole-genome sequencing using the Midnight amplicon scheme.

Genomic surveillance proved invaluable during the COVID-19 pandemic for tracking SARS-CoV-2 variants and guiding outbreak responses, underscoring the ongoing need to reduce whole-genome sequencing (WGS) costs and improve workflow efficiency to ensure accessibility in resource limited settings. Here, we evaluated a one-step reverse transcription polymerase chain reaction (RT-PCR) approach using the Midnight V2 primer scheme for targeted amplification of the SARS-CoV-2 genome, assessed its compatibility with Illumina sequencing, and compared its performance to a well-established two-step method. Initially, we determined optimal RT-PCR reaction conditions using the Midnight V2 primer panel for the one-step RT-PCR kit and scaled reaction volumes for both RT-PCR and library preparation. Clinical specimens (n&#x202f;=&#x202f;53) that had undergone routine WGS for surveillance purposes using the established two-step RT-PCR method were compared using the one-step RT-PCR assay. For samples with genome completeness greater than 70%, both methods gave comparable results with similar sequence coverage and 100% concordance for lineage assignment. Further investigation revealed a higher percentage of reads aligning to the SARS-CoV-2 genome with a greater depth of coverage using the one-step method compared to the two-step method. Finally, analysis of scaled one-step and library reaction volumes revealed significant cost savings for samples undergoing WGS. Overall, the results presented here verify the accuracy and reproducibility of one-step targeted amplification and offer an efficient and cost-effective workflow for routine SARS-CoV-2 genomic surveillance.

Humans

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

Routine methods misidentify Serratia spp.: Limitations of MALDI-TOF MS revealed by whole-genome sequencing.

Accurate species-level identification within the genus Serratia remains challenging due to extensive phenotypic overlap and high genomic relatedness among closely related and recently described taxa. This study presents an evaluation of routine and genome-based identification approaches applied to clinical Serratia isolates, integrating phenotypic assays, MALDI-TOF MS (Bruker Daltonics), 16S rRNA gene sequencing, and Whole-Genome Sequencing (WGS). A total of 103 isolates collected from a teaching hospital were analyzed. WGS was performed on a subset of isolates. Conventional biochemical methods classified all isolates as Serratia marcescens, whereas MALDI-TOF MS identified 60.1% as S. marcescens, 11.6% as S. ureilytica, and 28.1% just at the genus level. Peak analysis from MALDI-TOF MS revealed specific peaks associated with S. marcescens and S. ureilytica, but limited discriminatory power. WGS of six isolates initially identified as S. ureilytica by MALDI-TOF MS revealed reclassification as Serratia sarumanii (n = 5) and Serratia montpellierensis (n = 1), supported by Average Nucleotide Identity (ANI), Average Amino Acid Identity (AAI), and Digital DNA-DNA Hybridization (dDDH) thresholds. In contrast, 16S rRNA analysis showed limited species-level resolution. Phylogenomic and SNP-based analyses confirmed these classifications with strong support. Overall, this study underscores the critical role of high-resolution genomic approaches for precise species identification and highlights the need for continuous expansion and curation of MALDI-TOF MS reference databases to support reliable clinical diagnostics and epidemiological surveillance of emerging Serratia species.

Spectrometry, Mass, Matrix-Assisted Laser Desorpti

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

Conserved host-exclusive oligonucleotide motifs enriched in pathogenic genes of human oncogenic viruses.

Comparative viral genomics can reveal sequence-level constraints influencing virus-host interactions. Relative minimal absent words (rMAWs) are short oligonucleotide motifs present in viral genomes but completely absent from the host, potentially reflecting selective pressures related to host adaptation and immune evasion. Using the EAGLE algorithm and the GRCh38 human reference genome, we systematically screened for prevalent rMAWs (prMAWs) across six major human oncogenic viruses: Epstein-Barr virus (EBV), hepatitis B virus (HBV), hepatitis C virus (HCV), human papillomavirus (HPV), human T-cell leukemia virus type 1 (HTLV-1), and human herpesvirus 8/Kaposi's sarcoma-associated herpesvirus (HHV-8/KSHV). highly conserved 11- and 12-bp prMAWs were identified in EBV, HBV, HTLV-1, and HHV-8/KSHV, with sequence prevalences ranging from 91.5% to 97.9%. Conversely, no short prMAWs were detected in HCV or HPV, likely reflecting differences in genome architecture, mutation rates, and long-term host adaptation to the human host. Importantly, the identified host-exclusive motifs exhibited non-random genomic distribution and were preferentially embedded within viral genes central to replication, persistence, immune modulation, and oncogenesis, including EBNA-1 (EBV), HBx (HBV), Tax-associated regions (HTLV-1), and lytic replication genes of HHV-8/KSHV. Notably, all detected prMAWs were enriched in GC nucleotides and exhibited marked CpG over-representation, suggesting sequence constraints associated with epigenetic regulation and viral persistence. Collectively, these highly conserved, host-exclusive signatures offer promising, candidates for sequence-directed approaches in the diagnosis, monitoring, and investigation of virus-associated cancers.

Humans

A conserved distal-tail helical extension defines a tailspike attachment architecture in Gram-negative siphophages.

Rapid growth of bacteriophage genome collections has outpaced functional annotation of tail-tip proteins, limiting comparative analysis of host-recognition structures. Starting from a shared distal-tail gene organization in the Salmonella phages 9NA and Jersey, I developed a morphogenetic bioinformatic framework integrating gene synteny, sequence comparison, profile hidden Markov model (HMM) screening, structural evidence, structure-aware searching, and AlphaFold modeling. Comparison with the experimentally characterized lambda and Sf11 tail assemblies identified a predominantly alpha-helical C-terminal extension of the distal-tail (DT) protein associated with tailspike attachment, termed the distal-tail helical extension (DT-helix). Screening 541,986 proteins from 5167 complete NCBI RefSeq tailed-phage genomes, followed by evidence-based evaluation of sequence, genomic context, and structural architecture, identified 165 curated DT-helical-extension-associated phages. Their DT proteins segregated into six sequence groups. In the four principal multi-member groups, cognate tailspikes showed group-specific conservation in proximal N-terminal regions but substantially greater downstream diversity, consistent with sequence constraint at the DT-tailspike attachment boundary. A complementary ProstT5/Foldseek search supported the established groups but revealed no convincing additional highly divergent family. Together with the experimentally characterized Sf11 attachment interface, these findings define a recurrent morphogenetic architecture linking conserved distal-tail scaffolds to more variable receptor-binding proteins across siphophages infecting Gram-negative bacteria. Although universal exchangeability is not established, the identified scaffold-receptor-binding boundaries provide a framework for molecular characterization and rational phage engineering. Accession-level information for the 165 curated phages is available through PhageTailDB.

Viral Tail Proteins

Mitochondrial DNA diversity in Ecuadorian populations: Recurrence of variant 16136 within haplogroup B2.

The identification of lineage-defining variants, frequently found in the coding region of mitochondrial DNA (mtDNA), is essential for refining haplogroup classification. Most mtDNA studies in South American populations have focused on the control region (CR), which has provided important insights into population structure and maternal lineage origins, although information needed for more robust phylogenetic resolution has been neglected. This study investigates the maternal genetic structure of Ecuadorian populations by combining CR and whole mitogenome analyses. Sequences from the mtDNA CR were obtained from 461 individuals (253 Mestizos and 208 Native Americans), while complete mitogenomes were sequenced for 127 individuals to improve phylogenetic resolution by identifying lineage-defining variants present in coding region. Most mtDNA haplogroups in the two population groups analyzed were of Native American origin (A2, B2, B4, C1, D1, D4), with significant differences in the distribution of specific lineages between them. Among Mestizos, African haplogroups (all within the L branches) and Eurasian haplogroups (H, K, R, U) were detected at low frequencies, whereas no African lineages were observed among Native Americans. The results obtained highlighted a heterogeneity within Ecuadorian populations that must be considered when developing mtDNA haplotype databases for forensic purposes. Whole mitogenome sequences enabled the identification of variants that refined haplogroup classifications, provided a more accurate reconstruction of the maternal genetic diversity, and improve the discrimination between Native American and Asian maternal lineages within haplogroup B4b.

Humans

Comparative genomic epidemiology of food- and patient-derived diarrheagenic Escherichia coli from sentinel surveillance in Southeast China.

Diarrheagenic Escherichia coli (DEC) remains an important foodborne pathogen, yet long-term comparative genomic surveillance data jointly characterizing food-derived and patient-derived isolates remain limited. This surveillance-based comparative study integrated antimicrobial susceptibility testing and whole-genome sequencing to characterize diarrheagenic Escherichia coli isolates recovered from food and patient sources in Lishui, Southeast China, during 2018-2025, with emphasis on occurrence, resistance profiles, genomic backgrounds, and plasmid replicon-associated features. Antimicrobial susceptibility testing was performed for 258 selected isolates, and whole-genome sequencing was conducted for a curated analytical subset of 204 isolates. The sequenced subset was used for diversity-oriented comparative genomic analysis rather than for unbiased prevalence estimation of the entire DEC collection. EAEC predominated in both sources, although food-associated occurrence was heterogeneous across categories, with the highest recovery rate observed in raw meat. Patient-derived isolates showed a broader overall resistance burden, whereas food-derived isolates retained substantial resistance to tetracycline, chloramphenicol, and florfenicol. Phylogenetic analysis showed partial overlap in genomic backgrounds between food-derived and patient-derived isolates, while representative resistance determinants displayed both broadly distributed and lineage-enriched patterns. Replicon-based plasmid profiling identified 42 plasmid types, including 12 detected in both sources, with IncF-related replicons predominating among these shared profiles. Several food-derived isolates carried multiple plasmid replicon types that were also observed in patient-derived isolates. Overall, food-derived and patient-derived DEC showed partial overlap in genomic backgrounds, resistance determinants, and replicon-defined plasmid profiles within this surveillance setting, while retaining source-associated heterogeneity. These findings should be interpreted as surveillance-based comparative evidence rather than as evidence of direct source attribution or transmission.

Humans

Metatranscriptomic analysis of viral sequences associated with Culex nigripalpus at an Alabama aquaculture site.

Mosquitoes associated with aquaculture habitats can harbor diverse viruses, yet the viromes of many locally abundant species remain poorly characterized. At an aquaculture-associated site in Auburn, Alabama, we surveyed mosquito populations and found Culex nigripalpus to be the dominant species collected. To characterize viruses associated with this mosquito, we performed RNA-seq on pooled female Cx. nigripalpus and compared complementary bioinformatic workflows for viral detection and genome recovery. One workflow removed host-associated reads by mapping to the closest available mosquito reference genome prior to assembly, whereas a second workflow used fully de novo assembly and viral database annotation. Additional protein-level filtering, cross-workflow comparison, and comparison of Trinity and rnaSPAdes assemblies were used to prioritize well-supported viral candidates. Across the original analyses, 16 submitted accessions corresponding to 12 collapsed virus/name groups were recovered, including Merida virus, Hubei mosquito virus 5, Zhejiang mosquito virus, Hubei virga-like virus 3, Rinkaby virus, Elemess virus, Qingnian mosquito virus, Serbia narna-like virus 2, XiangYun narna-levi-like virus 8, Ecclesville picorna-like virus, and baculovirus-like fragments. Several candidates were supported across multiple workflows, while others were recovered only under specific analytical conditions, indicating that candidate recovery was influenced by assembly and filtering choices. Selected viral contigs were independently supported by RT-PCR amplification. Overall, these results provide a first characterization of viral sequences associated with Cx. nigripalpus from an Alabama aquaculture-associated site and show that comparison across assembly and filtering strategies helped prioritize the most consistently supported viral candidates.

Animals