SEARCH · Search PubMed

Results for “Genome-wide characterisation”

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.

279 records · Page 2Linked to original sources

PubMed2026

Exploratory identification and cellular functional characterization of ppiabl as a candidate gene associated with growth traits in Paralichthys olivaceus.

The Japanese flounder (Paralichthys olivaceus) is an important mariculture species. However, the genetic mechanisms underlying its growth traits remain poorly understood. To explore the genetic basis of growth variation, whole-genome resequencing was performed in a cultured cohort of 60 individuals, followed by exploratory genome-wide association analysis and candidate-gene prioritization. The results revealed heritability estimates of 0.40 for body weight and 0.24 for body length, with substantial overlap in associated loci between the two traits. Exploratory association and variant-annotation analyses prioritized ppiabl, which carries a nonconservative missense variant, as a candidate gene for further investigation. Tissue expression analysis showed that ppiabl was highly expressed in muscle tissue. This gene encodes a protein belonging to the conserved peptidyl-prolyl cis-trans isomerase family. In Japanese flounder primary muscle cells, ppiabl knockdown was associated with altered expression of growth-related genes and an increased G1-phase fraction, whereas overexpression produced changes in the opposite direction. In line with this, fast-growing individuals were found to have significantly larger muscle fiber areas than slow-growing ones. These findings suggest that ppiabl may be involved in muscle-related cellular processes associated with growth variation in Japanese flounder, although its contribution to whole-animal growth requires further validation. Overall, this exploratory study prioritizes ppiabl as a candidate gene potentially associated with growth-related cellular processes in Japanese flounder, although validation in larger independent populations and in vivo models is required.

Animals↗

PubMed2026

Cross-tissue multi-omics integration highlights BPHL and mitochondrial targets in Alzheimer's disease.

BACKGROUND: Mitochondrial dysfunction is a hallmark of Alzheimer's disease (AD), yet specific molecular targets remain to be fully characterized. METHODS: A summary-data-based Mendelian randomization (SMR) framework integrated AD genome-wide association study (GWAS) statistics (39,918 cases) with blood DNA methylation quantitative trait loci (mQTL), gene expression (eQTL), and protein (pQTL) data for 1136 mitochondria-related genes. Associations were assessed using Bayesian colocalization and HEIDI testing. Tissue relevance was evaluated in four brain regions (hippocampus, amygdala, cortex, frontal cortex) using GTEx and external transcriptomic datasets. RESULTS: Screening identified eight candidates supported across blood mQTL and eQTL layers. Stepwise central nervous system (CNS) evaluation singled out biphenyl hydrolase-like (BPHL) as the consistent candidate. Higher genetically predicted BPHL expression was associated with reduced AD risk across the hippocampus (OR=0.920, 95% CI 0.873-0.970), amygdala (OR=0.925, 95%CI 0.880-0.973), cortex (OR=0.943, 95% CI 0.908-0.978), and frontal cortex (OR=0.938, 95%CI 0.901-0.976). These findings aligned with protein-protein interactions connecting BPHL to respiratory complexes and lower BPHL expression in independent AD brains. Functional enrichment converged on oxidative phosphorylation pathways. CONCLUSIONS: By integrating multi-omics data with tissue-specific validation, this study nominates BPHL as a consistent protective candidate in the brain. These findings provide genetic support for mitochondrial molecular perturbations in AD, offering insights for future validation.

Alzheimer Disease↗

PubMed2026

Prioritizing Parkinson's disease risk-associated mitochondrial candidate genes via multi-omics integrative analysis.

BACKGROUND: Mitochondrial dysfunction has been implicated in Parkinson's disease (PD), but the genetically regulated mitochondrial genes associated with PD risk remain incompletely defined. METHODS: We conducted a summary-data-based genetic epidemiology study integrating summary-based Mendelian randomization (SMR), Heterogeneity in dependent instruments (HEIDI) filtering, and Bayesian colocalization to prioritize mitochondrial-related molecular features associated with PD risk. Mitochondrial-related genes were defined using MitoCarta3.0. Genetically predicted gene expression and plasma protein abundance were evaluated using expression quantitative trait loci (eQTL) data from eQTLGen and GTEx v8, and protein quantitative trait loci (pQTL) data was assessed using International Parkinson's Disease Genomics Consortium (IPDGC) as the discovery genome-wide association study (GWAS) and FinnGen as the replication dataset. Prespecified QTL analyses were interpreted using FDR correction, HEIDI filtering, and colocalization support. DNA methylation QTL analysis, mitochondrial phenotype MR, and single-nucleus RNA-seq analysis were performed as complementary analyses. RESULTS: In the primary eQTL analysis, higher genetically predicted TTC19 expression was associated with lower PD risk (OR = 0.80, 95% CI: 0.74-0.87, PPH4 = 0.80), whereas higher MALSU1 expression was associated with increased PD risk (OR = 2.21, 95% CI: 1.59-3.06, PPH4 = 0.96). Both associations survived FDR correction, passed HEIDI filtering, and showed colocalization support. GTEx whole-blood data supported the direction of the TTC19 association. No mitochondrial protein reached significance after FDR correction and colocalization filtering in the primary pQTL analysis. Complementary methylation analysis highlighted cg06270993 as an exploratory regulatory signal for MALSU1. CONCLUSIONS: This MR-colocalization study prioritizes TTC19 and MALSU1 as genetically supported mitochondrial-related candidate genes associated with PD risk. Further validation is required to define their functional roles in PD pathogenesis.

Humans↗

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

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

Integrative genomic and transcriptomic analyses identify key regulators of skin pigmentation in Larimichthys crocea.

The yellow body coloration of large yellow croaker (Larimichthys crocea) constitutes a crucial economic trait, yet its underlying genetic regulatory mechanisms remain poorly understood. This study systematically elucidated the molecular basis of body color variation by integrating genome resequencing and skin transcriptome analyses, combined with the contextual analysis of key pigmentation-related genes and phenotypic histological validation. 200 phenotyped individuals (including yellow-selected lines, F1 progeny, and normal control groups, all derived from a well-characterized aquaculture stock) identified 39 significantly associated SNPs (-log₁₀(P) ≥ 6), mapping to multiple candidate genes. These genes were significantly enriched in pathways related to pigment deposition (GO:0033059), melanosome organization (GO:0032438), melanogenesis, and tyrosine metabolism. Cross-developmental stage transcriptome analysis revealed 2395 differentially expressed genes (DEGs). Multi-omics integration identified eight overlapping candidate genes, including tyrp1, slc45a2, oca2, and dgat2, among which tyrp1 was prioritized for in-depth validation based on its core regulatory role in eumelanin synthesis, significant SNP association signal, and consistent downregulation in transcriptomic data. Experimental validation demonstrated that the g.895C > T mutation in exon 2 of tyrp1b was strongly significantly associated with the yellow phenotype: the frequency of mutant genotypes (TT/CT) reached 92.86%in the yellow-selected group, whereas the control group exclusively exhibited the wild-type genotype (CC). qPCR confirmed significantly downregulated tyrp1b expression in the skin of yellow individuals, consistent with the transcriptome trend. Histological and stereomicroscopic observations of skin tissues further validated the physiological basis of the yellow phenotype, revealing a significant reduction in melanophore number and abnormal melanosome morphology in yellow-phenotype individuals, accompanied by increased xanthophore density. These results suggest that tyrp1b mutation is strongly associated with the yellow phenotype. However, the presence of a wild-type CC individual in the yellow group indicates that this mutation is not strictly required for yellow coloration, suggesting that other genetic or environmental factors may also contribute to the phenotype, Additionally, downregulation of the carotenoid metabolism gene bco2 coupled with upregulation of xdh, together with the functional changes of slc45a2 and oca2, may synergistically promote xanthophore pigment deposition, contributing to the yellow phenotype. As melanin synthesis in large yellow croaker relies on the conserved tyrosinase pathway and transporter proteins, mutations in associated genes (tyrp1b, slc45a2, oca2) represent a primary underlying cause for the loss of melanin-based coloration and transition to a yellow phenotype in L. crocea. These findings provide key molecular targets and a theoretical foundation for molecular breeding of body color in this species, and also enrich the understanding of xanthism regulatory mechanisms in teleosts.

Animals↗

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

Introgression shapes the genomic conflict landscape of Malus, providing evidence for a reticulate backbone in a woody crop lineage.

Phylogenomic discordance is widespread across plants, but its evolutionary significance is often obscured when conflict is treated primarily as analytical noise rather than as evidence of underlying processes. In woody lineages in particular, incomplete lineage sorting, introgression, and genome duplication can interact over long timescales to produce complex genomic histories that are not adequately summarized by a strictly bifurcating tree. Here, we use Malus as a model woody genus to investigate how these processes structure conflict across a genus-scale, accession-based phylogenomic framework. Using broad taxon sampling, hundreds of nuclear loci, plastid genomes, and genome-wide SNP summaries, we reconstruct a robust nuclear backbone for sampled Malus lineages and evaluate where discordance is concentrated and which processes best explain it. Nuclear analyses resolve eight major clades, whereas conflict is non-random and localized to recurrent hotspots rather than evenly distributed across the tree. Cytonuclear discordance is similarly concentrated, especially around Clade H, represented by sampled accessions of M. tschonoskii, where localized plastid-nuclear disagreement is consistent with candidate plastid capture or organellar introgression. Multiple complementary analyses further indicate that the strongest conflict is not explained by ILS alone, but instead reflects lineage-structured introgression, while polyploid complexes represent additional localized sources of evolutionary complexity. Together, these results provide evidence for a reticulate genomic backbone in Malus and show how integrating nuclear, plastid, and genome-wide conflict analyses can help distinguish background discordance from process-specific signals in woody plant radiations. Several lineage-level reticulation hypotheses identified here should now be tested with broader population-level sampling and curated reference accessions.

Malus↗

PubMed2026

Clinical performance of two lithium disilicate CAD/CAM materials in posterior Class II inlay restorations: A 48-month randomised split-mouth clinical trial.

OBJECTIVES: To compare the clinical performance of Amber Mill (AM) and IPS e.max CAD (EM) lithium disilicate computer-aided design/computer-aided manufacturing (CAD/CAM) materials in posterior Class II inlay restorations and characterise their baseline properties. METHODS: Thirty-four adults received paired AM and EM posterior Class II inlays (68 restorations) in a triple-blind randomised split-mouth trial followed for 48 months. Restorations were evaluated at baseline and annually using revised World Dental Federation (FDI) criteria, with fracture and retention as the primary endpoint. Baseline characterisation included flexural strength, shear bond strength, translucency parameter, and scanning electron microscopy. McNemar, Wilcoxon signed-rank, Friedman, one-way analysis of variance, Tukey post hoc, and inter-rater agreement analyses were used. RESULTS: At 48 months, 18 paired participants were available for primary analysis. Failures occurred in 2 of 18 AM restorations and in 3 of 18 EM restorations, corresponding to success rates of 88.9% and 83.3%, respectively, with no significant between-material difference (McNemar p = 1.000). No catastrophic bulk ceramic fracture was observed. Secondary FDI scores remained mostly within the clinically acceptable range; marginal staining deteriorated over time in both groups (p < .001) without significant between-material differences. Baseline material testing showed significant material- and translucency-dependent differences in flexural strength, shear bond strength, and translucency. CONCLUSIONS: Within the limitations of the 48-month follow-up and the tested Class II inlay indication, AM showed clinical performance comparable to EM. Observed clinical complications were related to retention or marginal/interface behaviour. CLINICAL SIGNIFICANCE: For posterior Class II lithium disilicate CAD/CAM inlays, medium-term complications were mainly retention/interface-related, suggesting adhesive-interface durability may be as important as baseline ceramic strength.

Humans↗

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

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

Characteristics of post-exercise responders versus non-responders following aerobic or isometric exercise in physically inactive adults of African and South Asian descent with high-normal blood pressure or grade I hypertension.

OBJECTIVE: To investigate interindividual variability in post-exercise hypotension (PEH) and to characterise cardiovascular and autonomic differences between responders and non-responders following aerobic and isometric exercise in adults of African and South Asian descent with elevated blood pressure (BP). METHODS: Physically inactive adults of African and South Asian descent living in Suriname (18-65 years) with high-normal BP or grade I hypertension participated in a randomised controlled crossover trial. In this randomised cross-over trial, 47 adults (50.1 ± 10.8 years; 38% male) with high-normal blood pressure or grade I hypertension completed three conditions: aerobic exercise (30 min at 40-60% heart rate reserve), isometric handgrip exercise, and a non-exercise control. Ambulatory BP was assessed over 24 h. PEH was defined as the net effect: (post-exercise - pre-exercise) - (post-control - pre-control). Participants were classified as responders if daytime BP decreased ≥5 mmHg. Arterial stiffness, cardiac, and autonomic parameters were assessed. RESULTS: Following aerobic exercise, 46% of participants were classified as systolic responders compared with 28% after isometric exercise. No baseline differences were observed in demographic or clinical characteristics between responders and non-responders, suggesting that PEH variability may reflect underlying physiological rather than clinical differences. Aerobic responders demonstrated greater reductions in aortic augmentation index (-19.4% vs. -10.9%, p = 0.05), larger increases in stroke volume (+8.1 vs. -5.3 mL, p = 0.05) and cardiac output (+1.54 ± 1.89 vs. +0.58 ± 1.60 L/min, p = 0.009), and more favourable autonomic recovery. Among all variables, only the change in cardiac output was associated with PEH magnitude (r = -0.46, p = 0.006). No consistent physiological differences were observed following isometric exercise. CONCLUSION: PEH following aerobic exercise is characterised by a distinct responder phenotype associated with greater reductions in aortic augmentation index and favourable cardiac adaptations. These findings highlight substantial interindividual variability in BP responses and support the need for individualised exercise strategies in hypertension management.

Adolescent↗

PubMed2026

Comparative genomic and proteomic analysis reveals orthogroup structured evolution of tick protease inhibitors.

Protease inhibitors (PIs) play central roles in regulating endogenous proteolysis and host-parasite interactions in ticks. However, the evolutionary architecture underlying their diversification across tick lineages remains insufficiently resolved. Here, we performed a genome-wide comparative analysis of predicted proteomes from 14 tick species to systematically characterize PI repertoires. In total, 4931 putative PIs were identified and grouped into 20 families using the MEROPS classification system. Further, PI families such as Antistasin, WAP-type, and Pacifastin, which have not previously been systematically reported in tick genomes, were classified. Orthogroup inference demonstrated that PI expansion is structured at the level of evolutionary lineages rather than uniformly across families. By stratifying orthogroups according to duplication burden and taxonomic conservation, we identified a broadly conserved single-copy core under strong purifying selection. Motif level analysis of serpin reactive center loops further revealed conservation of inhibitory specificity within single copy orthogroups and diversification of key functional residues in duplication-associated lineages. Integration of secretion prediction and tissue-resolved proteomics from Hyalomma anatolicum and Rhipicephalus microplus demonstrated that evolutionary stratification is reflected at the protein level. Together, these findings provide an orthogroup-resolved evolutionary framework linking duplication dynamics, molecular evolution, and tissue-level protein deployment. This integrative approach offers a systematic basis for prioritizing conserved and diversified PI lineages for future functional and anti-tick intervention studies.

Animals↗

PubMed2026

Genomic identification and functional characterization of the nuclear receptor gene family in relation to sex determination and gonad development in the Pacific oyster (Crassostrea gigas).

Nuclear receptors (NRs) are a large superfamily of transcription factors that control a wide range of physiological processes by modulating the expression of downstream target genes. Numerous studies have confirmed that NR family members play critical and conserved roles in sex determination and gonadal development across metazoans. However, in mollusks, systematic characterization of NRs and their potential functions in gonadal regulation remain largely unexplored. In this study, 46 NR gene family members in the Pacific oyster (Crassostrea gigas) were identified and assigned to eight subfamilies. All NR family members contain at least one of the two core domains (DNA-binding domain, DBD; ligand-binding domain, LBD), and conserved exon-intron structures were observed within the same subgroup, indicating their evolutionary conservation. Furthermore, expression profiling revealed high expression of CgNR2F, CgNR5A1-1, and CgNR0B1 in undifferentiated gonads, suggesting their potential involvement in sex determination. CgNR1A and CgNR2E5 were specifically expressed in female gonads and exhibited female-biased expression patterns, indicating a putative role in ovarian development. Moreover, CgNR3A and CgNR3B showed high expression levels during the undifferentiated stage and early male development stage, implying their possible participation in male gonadal development and gametogenesis. These results expand the understanding of the NR gene family in C. gigas and help elucidate the potential functions of NR genes in sex determination and gonadal development.

Animals↗

PubMed2026

Identification and characterization of G protein-coupled receptors in the nocturnal halictid bee Megalopta genalis.

G protein-coupled receptors (GPCRs) are one of the largest families of membrane proteins in insects, regulating vision, neural signal transduction, and various physiological behaviors. Megalopta genalis exhibits a unique facultatively eusocial lifestyle and possesses adaptations for nocturnal activity; however, its GPCR family has not yet been systematically characterized. In this study, we performed genome-wide identification, phylogenetic analysis, and expression profiling of GPCRs in M. genalis by integrating genomic annotation and transcriptomic analysis. The results showed that a total of 99 GPCRs were identified in the genome of M. genalis, which were classified into four major families. Here, we show that M. genalis has undergone lineage-specific GPCR repertoire remodeling, marked by the expansion of novel orphan receptors and the systematic loss of multiple receptor subtypes, such as the neuropeptide receptors MIP-R and NPFR. Moreover, opsins have formed a diverse array of combinations and non-GPCR odorant receptors have undergone significant expansion via tandem duplication. Together, these features may represent part of the molecular repertoire associated with the adaptation of M. genalis to a nocturnal lifestyle. Furthermore, transcriptomic analysis revealed distinct spatiotemporal expression divergence within each of the Mth/Mthl and Fz GPCR families, suggesting functional specialization across development and adult tissues. This study provides the first systematic identification and initial functional characterization of GPCRs in M. genalis, revealing an evolutionary pattern characterized by the coexistence of contraction and expansion within the GPCR family. These findings lay a foundation for further studies aimed at elucidating the roles of these GPCRs in regulating M. genalis physiology and behavior.

Animals↗

PubMed2026

ChIP-seq profiling identifies diapause-regulated H3K27me3 targets in the fat body of Culex pipiens.

Culex pipiens, a principal vector of significant arboviruses, survives winter through diapause, a hormonally controlled inactive phase that enhances endurance under severe cold circumstances. Recent data suggests that epigenetic processes, namely histone post-translational modifications (hPTMs), play a crucial role in regulating seasonal dormancy. Prior studies from our laboratory indicated a decrease in the methylation of Histone 3 (H3K27me3) in diapausing fat body tissue, associated with elevated expression of the histone demethylase UTX. Nonetheless, the precise genomic areas impacted by these chromatin alterations remained unidentified. We used chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-seq) to delineate the genome-wide distribution of H3K27me3 across fat body chromatin in diapausing (D) and non-diapausing (ND) female Cx. pipiens. Notably, the higher signal at transcription start sites (TSSs) reflects localized redistribution rather than a global decrease, as diapausing fat bodies retain less H3K27me3 overall but concentrate it at promoters. To investigate the functional significance of these chromatin alterations, we confirmed a number of target loci via ChIP-qPCR and assessed gene expression with qRT-PCR. We identified many critical genes that were markedly increased in diapausing mosquitoes, exhibiting an inverse relation to H3K27me3 enrichment. Our data demonstrates different H3K27me3 chromatin landscapes between diapausing and non-diapausing Cx. pipiens, corroborating a hypothesis of selective, locus-specific repression in the non-diapause state and its targeted removal during diapause to permit activation of dormancy-associated genes. These results suggest that chromatin remodeling is a core driver of the diapause switch.

Animals↗

PubMed2026

Characterization and functional insights of histone deacetylases in bivalves: implications for temperature and immune response in Chlamys nobilis.

Histone deacetylases serve as pivotal epigenetic regulators that modulate chromatin remodeling and gene transcription, playing critical roles in immune defense and environmental stress responses in aquatic organisms. However, the evolutionary characteristics and functional roles of the HDAC family in bivalves remain poorly understood. In this study, genome-wide identification of the HDAC family across 30 bivalve species yielded 558 HDAC genes. Phylogenetic reconstruction categorized these genes into four conserved groups and revealed a unique, bivalve-specific SIRT8 clade. Using the noble scallop Chlamys nobilis as a representative model, expression profiling revealed distinct expression patterns among CnHDAC members. Class I and most Class III members were predominantly expressed in the gonads, while Class II members were enriched in immune-related tissues, implying their potential involvement in bivalve immunity. Upon temperature stress, CnHDAC1/2, CnHDAC11-1, CnHDAC11-2, CnSIRT2-1, CnSIRT4, CnSIRT6, and CnSIRT8-3 were significantly induced, highlighting their critical roles in temperature adaptation. Upon Vibrio exposure, CnHDAC1/2, CnHDAC8, CnSIRT4, and CnSIRT6 were upregulated, while CnHDAC4/5/7/9, CnHDAC6/10, CnSIRT2-2, CnSIRT5, CnSIRT7, and CnSIRT8-3 were downregulated, suggesting a coordinated epigenetic regulatory mechanism underlying host immune defense. In conclusion, this study systematically elucidates the evolutionary landscape of the HDAC family and underscores its potential involvement in environmental resilience and host immunity, providing a theoretical basis for the breeding of disease-resistant and stress-tolerant aquaculture bivalves.

Animals↗

PubMed2026

The effects of cold temperature on the development, microbiome, and transcriptome of the sea anemone Nematostella vectensis.

Thermal conditions impact essentially all aspects of the physiology for ectotherms. While the effects of high temperatures have been widely studied, cold temperature effects on aquatic invertebrates and their microbial communities have been poorly characterized. To determine the diverse effects of exposure to cold temperatures, we assessed acute and long-term impacts of ecologically relevant low temperatures on the development, microbiome, and gene expression of the sea anemone Nematostella vectensis. Two hours post fertilization, embryos were exposed to temperatures from 4°C to 35°C and development rate to the juvenile stage was quantified. We found temperature impacts the development rate of embryos, where lower temperatures extended development time and resulted in mortality below 10°C. For both microbiome and host transcriptomic responses, anemones were held at 20°C, 10°C, and 0°C and compared at 24 hours and 7 days. Extended exposures to colder temperatures caused restructuring of the host-associated microbiome, with the loss of common taxonomic groups from the class Bacteroidia and Bacilli. Lastly, cold stress induced significant changes in gene expression, which were more pronounced at the 10°C than 0°C but showed little change over time in each temperature. Interestingly, expression of genes associated with innate immunity were among the most differentially expressed genes including heat shock proteins and innate immune genes providing a potential host-imposed mechanism to explain the shift in the microbiome. Overall, cold temperatures have broad effects on many facets of this sea anemone and its microbial community and indicate the importance of cold temperature events when characterizing how ectotherms acclimate to thermal variation.

Animals↗