SEARCH · Search PubMed

Results for “Comparative phage genomics”

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.

588 records · Page 9Linked to original sources

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↗

PubMed2026

Integrated miRNA-mRNA profiling reveals candidate regulatory relationships associated with high-fat diet-induced muscle lipid deposition in black seabream (Acanthopagrus schlegelii).

High-fat diets are increasingly used in aquaculture due to their protein-sparing effects; however, the post-transcriptional regulatory mechanisms of fish muscle in response to high-fat diets (HFD) remain unclear. In this study, juvenile black seabream were fed either a normal-fat diet (NFD) or a HFD to investigate the miRNA-mRNA regulatory network associated with diet-induced muscle lipid deposition. Oil Red O staining and biochemical analysis showed that high-fat diet feeding markedly increased lipid droplet accumulation and crude lipid content in muscle, indicating significant induction of muscle lipid deposition. Integrated mRNA and miRNA expression profiling revealed substantial transcriptomic and post-transcriptional responses to high-fat diet challenge. A total of 271 differentially expressed genes were identified, including 120 upregulated and 151 downregulated genes. Through combined target prediction and expression correlation analysis, thirteen candidate inverse miRNA-mRNA relationships were subsequently identified, and RT-qPCR supported the expression patterns of selected miRNAs and mRNAs. These pairs included miR-499-x-dmgdh, miR-499-y-gatm, miR-727-y-ass1, miR-4649-x-foxo4, miR-9129-z-myl7, and several novel miRNA-mediated interactions involving adk, chst11, lypla2, frem2, kcnc4, wars1, bag2, and capn2. Functional analysis suggested that these regulatory pairs were mainly associated with metabolic adaptation, structural remodeling, and cellular stress responses. In particular, gatm, dmgdh, ass1, and adk were associated with energy metabolism-related processes, including pathways previously linked to Ampk regulation, whereas myl7, frem2, and kcnc4 may contribute to muscle structural maintenance and excitability regulation. Overall, this study provides candidate miRNA-mRNA regulatory relationships potentially involved in high-fat diet-induced muscle lipid deposition and adaptive remodeling in black seabream, offering a basis for future functional studies on muscle metabolism and quality regulation in marine fish.

Animals↗

PubMed2026

From commensal to pathobiont: The emergence of virulence-enhanced Escherichia coli in China's food-animal systems - insights with future implications.

A fundamental shift in Escherichia coli epidemiology is being driven by convergence of virulence determinants and antimicrobial resistance within linked human-animal-environment systems. In China, the rapid growth of food-animal production, extensive antimicrobial use, and complex food networks are accelerating the emergence and dissemination of virulence-enhanced E. coli pathobionts. This review synthesizes recent epidemiological, genomics, and outbreak data to characterize China's evolving landscape of food-animal-associated E. coli. We highlight a significant shift from classical pathotypes to hybrid lineages that simultaneously carry virulence factors and last-resort antibiotic resistance determinants, including mcr-1, tet(X4), and blaNDM. These traits disseminate rapidly via plasmid-mediated horizontal gene transfer, facilitating rapid adaptation and enabling cross-sectoral One Health transmission. National surveillance, foodborne outbreak investigations, and whole-genome sequencing data show that food-animal reservoirs are active evolutionary niches that drive pathogen diversity and fitness, rather than serving merely as contamination sources. Whole-genome sequencing also pinpoints high-risk clones (e.g., ST394) and plasmid-mediated co-selection of virulence and AMR. The emergence of hybrid pathotypes (e.g., STEC/ETEC) and AMR-virulence co-selection challenges traditional classification and limits the effectiveness of conventional surveillance approaches. The 2017 colistin ban reduced mcr-1, yet ongoing resistance and emerging tet(X4) demand integrated surveillance. Collectively, these findings call for reconceptualizing E. coli as a dynamic genomic entity embedded within a unified ecological network. Addressing this threat requires an integrated One Health strategy including genomic surveillance, agricultural antimicrobial stewardship, and coordinated food-environment-clinical monitoring to prevent high-risk clone emergence and global spread.

Animals↗

PubMed2026

Norgestrel drives mitochondrial collapse and plasma membrane impairment in Pacific oyster (Crassostrea gigas) sperm by triggering premature acrosome reaction.

The toxic mechanisms of norgestrel (NGT), an emerging marine pollutant, on the sperm from externally fertilized invertebrates remain elusive. This study employed an integrated physiological and multi-omics framework to elucidate how NGT (10 and 1000 ng/L) disrupts acrosome reaction (AR) signaling machinery, thereby impairing the functional integrity of Pacific oyster (Crassostrea gigas, also known as Magallana gigas) sperm. Exposure to NGT triggered a significant, dose-dependent premature AR, characterized by elevated acrosin activity and a loss of acrosomal integrity. Multi-omics integration supports a model in which this premature exocytosis is linked to signaling disturbances, including disruption of calcium signaling and reduced transcript abundance of calmodulin (CaM) and the primary recognition protein zonadhesin (Zan). This signaling interference induced an premature AR, subsequently driving a cascade of bioenergetic and structural failures. At the mitochondrial level, NGT induced abnormal mitochondrial permeability transition pore (mPTP) opening and elevated the transcript levels of antioxidant defense genes (e.g., peroxiredoxin-5, PRDX5). These alterations indicate the occurrence of mitochondrial collapse. Concurrently, scanning electron microscopy verified localized plasma membrane wrinkling and pore formation in sperm. In addition, NGT exposure decreased the transcript abundance of cytoskeleton-related genes, including solute carrier family 26 member 6 (SLC26A6), actin (ACT), and tubulin polymerization promoting protein family member 3 (TPPP3). These molecular changes further disrupted membrane phospholipid homeostasis, as represented by altered glycerophospholipid metabolism. At the same time, cumulative cellular stress was associated with decreased transcript abundance of cytoprotective factors (e.g., baculoviral IAP repeat-containing proteins, birc2) and changes in apoptosis-related genes consistent with activation of a caspase-8-mediated apoptotic programme. In conclusion, NGT, as a representative synthetic progestin, exerts reproductive toxicity by interfering with signaling mediators to induce premature AR, which subsequently exhausts metabolic energy and triggers plasma membrane impairment. These findings provide a critical mechanistic basis for the aquatic ecological risk assessment of synthetic progestins.

Animals↗

PubMed2026

Transcriptomic responses of gill and intestinal tissues in Nile tilapia (Oreochromis niloticus) to bacterial infection following sequential nanoimmersion and hydrogel-based multivalent vaccination.

Bacterial pathogens, including Flavobacterium oreochromis, Aeromonas veronii, Streptococcus agalactiae, and Edwardsiella tarda, represent major infectious threats to Nile tilapia (Oreochromis niloticus). A multivalent vaccination strategy integrating cationic nanoemulsion immersion with oral hydrogel boosters was developed to investigate tissue-specific immune responses at the transcriptomic level. Gill tissues were collected following immersion challenge and intestinal tissues following intraperitoneal injection challenge, reflecting the physiologically relevant infection biology of each pathogen and the mechanistic rationale of each delivery platform. RNA sequencing (RNA-seq) generated high-quality datasets (mapping rate > 81.64%) with strong concordance to quantitative real-time PCR (qRT-PCR) validation (r = 0.83). Comparative transcriptomic analysis revealed distinct yet complementary immune signatures between tissues. Gill transcriptomes were enriched in phagosome, focal adhesion, extracellular matrix-receptor interaction (ECM-receptor interaction), and cytokine-cytokine receptor interaction pathways, accompanied by increased expression of major histocompatibility complex class I/II (MHC class I/II), mannose receptor, αVβ3 integrin, and calnexin, indicating innate activation, enhanced phagocytic capacity, epithelial barrier reinforcement, and adaptive immune coordination. Intestinal transcriptomes showed predominant enrichment of adaptive immune pathways, including the intestinal immune network for immunoglobulin (Ig) production, Forkhead box O (FoxO) signaling, and mitogen-activated protein kinase (MAPK) signaling, with increased expression of T-cell receptor (TCR), inducible T-cell co-stimulator ligand (ICOS-L), C-X-C chemokine receptor type 4 (CXCR4), and polymeric immunoglobulin receptor (pIgR), reflecting T and B cell coordination, lymphocyte trafficking, and mucosal immunoglobulin transport, alongside innate engagement through phagosome pathway enrichment. Shared upregulation of MHC class II, B-cell receptor (BCR) signaling, integrin alpha M (ITGAM), and immunoglobulin-associated components across both tissues suggests coordinated mucosal immune activation through a conserved immune module, warranting direct experimental validation. Collectively, these findings provide transcriptomic evidence that this vaccination strategy elicits an integrated, tissue-specialized immune response, advancing mechanistic understanding of gill and intestinal immunity in vaccine-induced protection of teleost fish.

Animals↗

PubMed2026

Retinal microstructural alterations as early phenotypes of depression in radiogenomics analysis.

BACKGROUND: With the increasing prevalence of depression, there is an urgent clinical need for early screening in depression. The retina offers a promising window for early screening in depression due to its rapid, non-invasive, objective, eye-brain correlated characteristics, but previous research has yielded conflicting alterations in retinal microstructure in depression. METHODS: We screened retinal optical coherence tomography and brain magnetic resonance imaging data in the UK Biobank to enroll 23,225 participants for retinal study of depression occurrence, and 1475 participants for the eye-brain association study. We also used genetic data (ID: ebi-a-GCST90014267 and ukb-d-20,448) from the Integrative Epidemiology Unit Open Genome-Wide Association Study for Mendelian randomization analysis. We used Cox regression to assess the association between retinal microstructure and depression risk, Mendelian randomization to infer causality, and mediation analysis to explore retina-brain pathway association. RESULTS: The Cox regression analysis showed that retinal ganglion cell-inner plexiform layer (GCIPL) thickness remained a significant predictor of depression. The Mendelian randomization analysis indicated a positive statistical association between GCIPL thickness and depression. Moreover, there was a significant positive correlation (all p < 0.001) between the volume of specific depression-related brain regions and the GCIPL thickness. Adjusting for age, sex, and head size, the mediation analysis provided preliminary evidence for a potential anatomical pathway linking retinal GCIPL thickness to depression-related brain regions through primary visual cortex and secondary visual cortex volumes. CONCLUSION: Thickened retinal GCIPL is a potential early phenotype of depression and has a potential association pathway with depression-related brain regions using a radiogenomics approach.

Humans↗

PubMed2026

Quality assessment, prognostic factors, and biomarkers for brain tumor analysis: a comprehensive systematic review.

The brain tumors possess different causative factors and properties, making their diagnosis and treatment difficult. Growth of these cancers usually leads to compression of the adjacent nerves and obstruction of the flow of cerebrospinal fluid, thus leading to increase in intracranial pressure. This affects the working of brain in many ways; thus, the difficulty involved in its treatment. With the improvements in technology in neuroimaging, including Diffusion Tensor Imaging (DTI), Positron Emission Tomography (PET), and multiparametric Magnetic Resonance Imaging (mpMRI), the diagnosis process has become easy. The effectiveness of any form of therapy in such patients depends primarily on their prognosis. While it is a common practice that physicians determine the prognosis of the disease by considering the age of the patient, histological grade of the tumor, and resection status, now this method has become more comprehensive by adding molecular signature and genetic analyses to the list of criteria. Next-generation sequencing (NGS) allows a reliable molecular classification. It increases the level of risk stratification, facilitating the application of therapies tailored to individual patients. Thus, molecular oncology has greatly changed our views on brain tumors' pathology and prognosis while neoadjuvant treatments aim at increasing the survival rate. On the other hand, radiogenomics is a field of study that combines non-invasive imaging phenotypes and genomic information in order to find unique molecular signatures of tumors without collecting samples from tumors. Molecular biomarkers are absolutely essential in the diagnosis of cancer, treatment monitoring, and recurrence of cancer. Advances in liquid biopsy technology, particularly the methods for circulating tumor DNA (ctDNA) and Extracellular Vesicle (EV) based analysis, have enabled the possibility of non-invasive monitoring of the progression of the tumors over time. This review highlights key studies and important scientific works about imaging technologies, biomarkers, and prognostic factors of malignant brain tumors.

Humans↗

PubMed2026

Effects of acute hypoxia followed by reoxygenation on intestinal histomorphology, oxidative stress and hypoxia signaling biomarkers, and microbiota in pikeperch (Sander lucioperca).

In aquatic environments, natural and anthropogenic factors commonly reduce dissolved oxygen (DO) and trigger hypoxia, which threatens the health and survival of aquatic organisms. As an important economic fish species in China, pikeperch (Sander lucioperca) is extremely sensitive to hypoxia. However, there are relatively few reports on how hypoxia and reoxygenation affect its intestinal physiology and microbial community. Three treatment groups were set for pikeperch: normoxia (DO = 8.5 ± 0.5 mg/L), 48 h hypoxia (DO = 2.5 ± 0.1 mg/L), and reoxygenation (48 h hypoxia followed by 6 h reoxygenation at normal DO), to evaluate alterations in intestinal histopathology, tight junction gene expression, oxidative stress, hypoxia signaling molecules and intestinal microbiota composition. The results showed that hypoxia significantly decreased muscularis thickness by approximately 32.5% and reduced the expression of tight junction genes (Occludin, Claudin2, and ZO-2). Moreover, hypoxia significantly increased oxidative stress index levels (GSH-Px, CAT, and MDA), markedly upregulated the expression of Bax, Caspase3, and HIF-1α, while significantly downregulating the expression of Bcl-2, Egln1, and Egln2. Notably, reoxygenation elicited partial compensatory effects against these hypoxia-induced changes. 16S rRNA sequencing analysis revealed that hypoxic stress altered the intestinal microbial community composition of pikeperch and increased its diversity. In the hypoxia group, the abundance of the phylum Bacillota, along with the genera Halomonas and Acinetobacter, was significantly elevated, whereas in the reoxygenation group, the genus Lactobacillus increased approximately 180-fold. The results indicated that hypoxia caused intestinal oxidative damage, cell apoptosis, and intestinal microbiota dysbiosis in pikeperch, while short-term reoxygenation achieved partial recovery from these hypoxia-triggered intestinal injuries. The present research provides valuable references for in-depth exploration of the molecular mechanisms behind the response of pikeperch to acute hypoxia and reoxygenation stress, while also offering a novel perspective to understand the mechanism by which hypoxia impacts intestinal health in fish.

Animals↗

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

Effects of faba bean-based crisping culture on phenotypic characteristics, muscle quality, and serum metabolome in Nile tilapia: Screening biomarkers to assess the degree of crisping.

Feeding Nile tilapia (Oreochromis niloticus) a faba bean-based crisping diet enhances muscle hardness (crispness) and overall flesh quality. However, the underlying mechanisms and reliable biomarkers remain insufficiently defined. This study integrated phenotypic traits, muscle texture, collagen content, serum antioxidant enzyme activities (SOD, CAT, and GSH-Px), MDA levels, and serum metabolomics to understand the determinants of muscle crisping. Fish were assigned to a crisping diet or a control group for 90 days. Individuals in the crisping group were implanted with passive integrated transponder (PIT) tags to enable correlation analyses among phenotypic traits (body weight/length/frame changes), serum indicators (NAM, FAD, and GSH-Px) and muscle hardness. Compared with controls, the crisping diet significantly increased muscle hardness, gumminess, and chewiness, accompanied by elevated collagen content. Antioxidant profiles were altered, with higher activities of serum SOD and CAT, together with elevated MDA levels and reduced GSH-Px activity (P < 0.05). Metabolomic analysis identified 830 differential metabolites (682 upregulated and 148 downregulated), predominantly comprising carboxylic acids and derivatives, glycerophospholipids, and benzene derivatives. Enrichment analysis indicated significant involvement in general metabolic pathways, ATP-binding cassette (ABC) transporters, amino acid biosynthesis, and glycine, serine, and threonine metabolism (P < 0.05). Notably, acetylpyruvate was upregulated in glutathione metabolism, nicotinate and nicotinamide metabolism, and galactose metabolism; pantothenic acid was upregulated in glycine, serine, and threonine metabolism; whereas δ-tocotrienol was downregulated. Correlation analysis revealed weak negative associations between muscle hardness and phenotypic traits (body weight/length/frame changes D5-7, D5-10, D7-8) (P < 0.05). In contrast, serum NAM and FAD were weakly positively correlated with muscle hardness, whereas GSH-Px showed a weak negative correlation (P < 0.05). Collectively, these findings suggest that body weight, body length, frame measurements (D5-7, D5-10, and D7-8), and serum NAM, FAD, and GSH-Px are associated with the degree of muscle crispness in Nile tilapia fed a faba bean-based crisping diet and may serve as candidate biomarkers under these culture conditions.

Animals↗

PubMed2026

Peripheral pain threshold, glycaemic status, and LAMP3 genetic variation: A community-based analysis.

Diabetic polyneuropathy is a common complication of diabetes, yet substantial inter-individual variation in peripheral pain perception suggests underlying genetic influences. This population-based study investigated clinical, metabolic, and genetic determinants of pain threshold using intraepidermal electrical stimulation in 906 participants from the Iwaki Health Promotion Project 2017. Genome-wide association analysis identified 12 loci showing suggestive associations, among which a missense variant in LAMP3 (rs482912) was prioritized as a biologically plausible candidate. Phenotype-stratified analyses showed that individuals carrying the CT or CC genotypes had lower PINT indices than those with the TT genotype, indicating reduced pain thresholds. Notably, the CC genotype retained an association with lower pain threshold using intraepidermal electrical stimulation under conditions of metabolic stress, including impaired glucose tolerance, elevated HbA1c, and obesity, whereas this association was attenuated in the presence of hypertension. Single-cell RNA sequencing analysis of human skin revealed that LAMP3-positive mature dendritic cells, enriched in immunoregulatory molecules, exhibited transcriptional enrichment of inflammatory, antigen-presenting, and nociception-related pathways, including NF-κB, JAK-STAT, cytokine signaling, and neuroimmune sensitization cascades. Autopsy-based skin analysis further demonstrated genotype-associated differences in dermal LAMP3-positive cell infiltration and CD8-positive T-cell abundance, while CD4-positive T-cell abundance and intraepidermal nerve fiber density remained unchanged across genotypes. Taken together, these findings suggest a potential association between LAMP3 variation and individual differences in peripheral pain threshold and provide biological context supporting a role for neuroimmune interactions in early sensory modulation under metabolic stress. Given the suggestive genetic evidence and indirect mechanistic data, these observations should be interpreted as exploratory and hypothesis-generating.

Humans↗

PubMed2026

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↗

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

Genetic evidence for a causal relationship between melatonin metabolism and depression.

To investigate the causal relevance of melatonin metabolism, which provides the biological basis for circulating melatonin levels, to specific depression symptom subtypes, we performed a targeted systematic review of melatonin metabolism pathways in the human brain and liver. Using two-sample Mendelian randomization (MR), we assessed the causal effects of metabolism pathways and/or individual genes on major depressive disorder (MDD) and nine symptom subtypes derived from Patient Health Questionnaire-9 (PHQ-9). Instrumental variables (IVs) were expression quantitative trait loci (eQTL) for eight individual genes, one synthesis route, and three degradation routes. Results were assessed using Bayesian colocalization and phenome-wide association analyses. At the pathway-level, the genetically proxied synthesis-route signal was associated with PHQ-9 Assessment 5 (PHQ9A5, OR: 0.89, 95% CI: 0.85-0.93), but sensitivity analyses suggested this association was primarily driven by TPH1 and may reflect serotonin-related biology. In contrast, higher brain melatonin degradation raised the risk of both PHQ9A1 (OR: 1.03, 95% CI: 1.02-1.04) and PHQ9A7 (OR: 1.03, 95% CI: 1.02-1.03). Within degradation, up-regulation of the kynurenine sub-pathway increased the odds of PHQ9A3 (OR: 1.05, 95% CI: 1.02-1.07), PHQ9A4 (OR = 1.04, 95% CI: 1.02-1.06) and PHQ9A7 (OR: 1.05, 95% CI: 1.02-1.07). Gene-level analyses were largely concordant, except for SULT1A1, whose higher expression was genetically protective for PHQ9A3 but risk-increased for PHQ9A1 and PHQ9A4. Overall, these results demonstrate that melatonin metabolism exerts symptom-specific and pathway-specific causal effects on depression. A stratified view of melatonin's role may help optimize the application of exogenous melatonin supplementation.

Melatonin↗

PubMed2026

Beyond species trees: pervasive gene flow limits phylogenomic resolution in the diversification of Juniperus from the Qinghai-Tibet Plateau.

Understanding how lineages diversify despite persistent ancestral polymorphism and recurrent gene flow remains a central challenge in evolutionary biology. Juniperus distributed across the Qinghai-Tibet Plateau provide an ideal system for addressing this question because repeated geological uplift and climatic oscillations have likely promoted cycles of lineage divergence, range shifts, and secondary contact. Here, we combined approximately 1.08 million genome-wide SNPs from 164 individuals representing thirteen Juniperus lineages with phylogenomic datasets comprising 3,381 nuclear single-copy genes and nearly complete plastomes. We detected extensive phylogenomic discordance and cytonuclear incongruence across genomic datasets. Topology weighting, coalescent simulations, quartet-based tests, and analyses of gene flow and reticulation collectively support the interpretation that these patterns were shaped by the combined effects of prolonged incomplete lineage sorting and gene flow during lineage diversification. Ecological niche analyses further provide a spatial and climatic context in which environmentally similar lineages may have had greater opportunities for secondary contact during historical range shifts. Collectively, our results reveal that the evolutionary history of Qinghai-Tibet Plateau Juniperus is characterized by reticulate diversification rather than strictly bifurcating evolution, and demonstrate how genome-wide discordance can provide biological insights into the evolutionary processes underlying lineage diversification.

Gene Flow↗

PubMed2026

Stage-specific ROMO1 in rheumatoid arthritis: predictive immune insights into the MIF pathway and HLA-DR/IL2RA axis via integrated GWAS, transcriptomic, single-cell, and spatial profiling.

Emerging evidence links reactive oxygen species modulator 1 (ROMO1), a key mitochondrial ROS regulator, to rheumatoid arthritis (RA) pathogenesis. However, its exact mechanism remains elusive given the conflicting evidence about its specific function. We used a four-level integrative framework combining multi-omics data and literature‑supported mechanistic inference. At the genetic level, Mendelian randomization (MR) was performed to explore potential causal relationships between ROMO1, IL2RA, HLA-DR, MIF, and RA risk, followed by differential expression analysis and machine learning-based feature selection to identify key mROS genes. The temporal expression dynamics of ROMO1 were assessed in RA progression. At the cellular and tissue levels, we integrated single-cell RNA sequencing and spatial transcriptomics to map cell-type-specific expression and synovial localization of ROMO1-related immune cells and pathways. Finally, our multi-omics findings were contextualized with literature-supported mechanistic inference. (1) MR results were consistent with a potential protective effect of ROMO1 on RA (OR = 0.52) and its potential regulation of risk factors IL2RA (OR = 0.46) and HLA-DR (OR = 0.40). Conversely, IL2RA (OR = 1.42), HLA-DR (OR = 1.88), and MIF (OR = 1.17) were positively associated with RA risk. Additionally, ROMO1 was identified as a top candidate diagnostic predictor with stage-specific dynamics: downregulated in the early but upregulated in the late/remission stages. (2) Single-cell RNA sequencing showed ROMO1's cell-specific expression in CD14+ HLA-DR+ CD74+ monocytes and CD4+ IL2RA+ T cells. Cell communication analysis further suggested that these cells may participate in MIF pathway regulation. Spatial transcriptomics subsequently identified that ROMO1-related cells localized to synovial pathological regions, with MIF pathway changes correlated with RA progression. (3) Finally, literature-supported mechanistic inference suggests that ROMO1 may modulate mROS levels to promote anti-inflammatory M2 macrophage polarization, which could theoretically contribute to reduced systemic inflammation and the alleviation of multi-organ decline in RA. This integrated multi-omics investigation, supported by literature-based mechanistic inference, suggests ROMO1 as a stage-dependent biomarker candidate and potential immune regulator in RA.

Humans↗

PubMed2026

Construction of an infectious clone of Spodoptera frugiperda densovirus and its biological characteristics.

Densoviruses are highly pathogenic to their insect hosts and have great potential for biocontrol. Spodoptera frugiperda densovirus (SfDV) was isolated from diseased larvae of Spodoptera frugiperda, while its biological functions remain unclear. Herein, we successfully constructed an infectious clone of SfDV. The S. frugiperda larvae transfected with the infectious clone exhibited anorexia, stunted growth, and reduced activity. Histopathological analysis further showed that the epidermis, fat body and trachea were infected instead of muscle and midgut tissues. Transmission electron microscopy (TEM) revealed that numerous virions of about 22 nm were distributed within both the nucleoplasm and cytoplasm of epidermal cells. Moreover, many virions were also found contained within vesicles in the cytoplasm. The replication kinetics of the rescued SfDV (rSfDV) was similar to that of the parental SfDV. The median lethal dose (LD50) and median lethal time (LT50) values of rSfDV were 6.63 × 107 viral genome copies (vgc), 5.23 d, respectively, which were also comparable to those of the parental SfDV. Taken together, the infectious clone of SfDV provides an important tool for further exploring the genome function, pathogenesis, and interactions with its hosts.

Animals↗

PubMed2026

A complete hlyCABD-like RTX operon marks a virulence-associated subset of trh-positive Vibrio parahaemolyticus from Hangzhou Bay, China.

Vibrio parahaemolyticus remains a major cause of seafood-associated gastroenteritis, yet routine surveillance still relies largely on the canonical hemolysin markers thermostable direct hemolysin (tdh) and tdh-related hemolysin (trh). To determine whether this framework overlooks accessory virulence determinants in trh-positive lineages, we analyzed 193 V. parahaemolyticus isolates collected between 2022 and 2025 from clinical, environmental, and seafood-associated sources in the Hangzhou Bay region of China. Serotyping identified 45 serotypes, with O10:K4 predominating among clinical isolates. Both clinical and non-clinical populations showed open pan-genomes, although the non-clinical group carried a larger accessory gene pool. We identified a complete hlyCABD-like RTX operon in 10 trh-positive isolates with T3SS2-associated virulence backgrounds. These RTX-positive isolates were distributed across seven sequence types and three of five phylogenetic groups. This distribution was lineage-restricted but non-clonal. In the representative hybrid-assembled genome, the operon occurred within a mosaic genomic region containing additional virulence- and mobility-associated genes, indicating a composite pathogenicity island-like element. In the tested subset, RTX-positive isolates showed significantly greater hemolytic activity than RTX-negative trh-positive isolates. This significant difference was consistently observed in both plate-based and liquid assays, and within the RTX-positive subset, hlyA expression correlated with hemolytic activity, whereas the trh gene and the tlh (thermolabile hemolysin) gene did not. A complete hlyCABD-like RTX operon therefore identifies a hemolysis-associated subset of trh-positive V. parahaemolyticus and supports its further evaluation as an additional target for food safety surveillance.

Vibrio parahaemolyticus↗