SEARCH · Search PubMed
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.
Explore the source record for details and available documents.
INTRODUCTION: Polycystic ovary syndrome (PCOS) is a multifactorial endocrinological disorder with a substantial genetic component. However, the role of genes involved in follicular development and androgen regulation remains incompletely understood, particularly in South Indian populations. This study aimed to evaluate how variations in the ERBB4 and SHBG genes affect PCOS risk. METHODOLOGY: A hospital-based case-control study was conducted among 400 South Indian women, comprising 200 women with PCOS and 200 age-matched healthy controls. Genomic DNA was extracted to study SNPs at ERBB4 (rs2178575 and rs1351592) and SHBG (rs1799941 and rs727428) using ARMS-PCR genotyping. The study compared genotype and allele frequencies between cases and controls while assessing their associations with allelic, homozygous, heterozygous, dominant, recessive, and over-dominant genetic models. Genotyping accuracy was confirmed by re-genotyping and Sanger sequencing of a subset of samples. RESULTS: The ERBB4 rs2178575 polymorphism demonstrated a significant association with PCOS, as the AA genotype and A allele combination increased risk across all three genetic models, including homozygous, recessive, and allelic models. The ERBB4 rs1351592 variant was associated with 3-fold higher risk of PCOS in heterozygous and GC carriers. The SHBG rs1799941 polymorphism showed a significant link to PCOS through its effects on heterozygous and allelic states, whereas rs727428 displayed no significant connection due to its monomorphic distribution. CONCLUSION: These findings suggest that polymorphisms in ERBB4 and SHBG may contribute to PCOS susceptibility in South Indian women in a locus- and model-specific manner, revealing the intricate genetic structure that defines this medical condition.
The beet armyworm, Spodoptera exigua, is a polyphagous insect whose development and seasonal occurrence are strongly influenced by temperature. However, transcriptomic responses to developmental thermal regimes remain insufficiently characterized in field-collected populations. In this study, we compared two Korean field-collected populations of S. exigua: a Haenam population collected in May and initially maintained at 15 ± 1 °C (HN), and a Jeju population collected in July and initially maintained at 27 ± 1 °C (JJ). F1 larvae from each population were reared under three fluctuating developmental temperature regimes: low (15-21 °C), middle (21-27 °C), and high (27-33 °C), followed by RNA-seq analysis. Differential expression analysis revealed population-associated variation in transcriptomic responses across developmental temperatures. HN exhibited a larger number of differentially expressed genes under the high-temperature regime, suggesting stronger transcriptomic sensitivity to elevated developmental temperature. Functional enrichment analyses identified population-associated differences in pathways related to heat response, oxidative metabolism, cytoskeletal organization, cuticle-associated processes, lipid metabolism, and immune-related functions. In JJ, heat-response and cuticle-related expression patterns were more prominent under warmer developmental conditions, whereas HN showed broader changes in stress- and metabolism-associated pathways under high temperature. Overall, this study provides a comparative transcriptomic analysis of two field-collected S. exigua populations under different developmental temperature regimes and identifies RNA-seq-based molecular response patterns associated with population-specific thermal response profiles.
BACKGROUND AIMS: Cell and gene therapy products face unique regulatory challenges due to their biological complexity and the stringent expectations for manufacturing control, analytical testing, and clinical assessment. As the pipeline grows, understanding the drivers of FDA non‑approval is increasingly important for improving first‑cycle success and reducing development delays. In this study, we aimed to characterize the frequency and impact of refuse‑to‑file actions, major amendments, and complete response letters issued for biologics licensing applications; and to identify recurring patterns of deficiencies contributing to delayed approval. METHODS: Publicly available documents were extracted from the FDA website. Data from each approved cell and gene therapy product's regulatory history-including complete responses, major amendments, inspection timing, and cited deficiencies-were compiled and categorized across clinical, quality, and labeling domains. RESULTS: Analysis showed that informational deficiencies were common across modalities. Major amendments occurred in roughly half of applications, and complete responses had the most significant impact. Quality deficiencies appeared in all complete responses and were the predominant barrier, while clinical and labeling issues were less frequent but meaningful when present. CONCLUSION: Overall, these findings highlight the need for proactive FDA engagement, comprehensive readiness, and early inspection preparation to reduce regulatory risk and improve first cycle approval outcomes.
Numerous adaptations to interventions have been included in cognitive behavioral therapy (CBT) for autistic youth. This study examines the degree to which CBT adapted to the social needs of autistic youth confers significant benefit by promoting social responsiveness. A secondary analysis was conducted on a multisite randomized clinical trial (Wood et al. in JAMA Psychiatry 77:474-483, 2020) comparing adapted CBT with standard-of-practice CBT and treatment-as-usual (TA). Autistic youth (N = 167; aged 7-13) with maladaptive and interfering anxiety participated. The adapted CBT (BIACA) uses a modular format, with an emphasis on supplementing common CBT practice elements (reframing and graded exposure) with social skill supports. The primary outcome measure was the Pediatric Anxiety Rating Scale. Social responsiveness was assessed with the Social Responsiveness Scale, Second Edition. Participants' general mental health was assessed as a secondary outcome using the Brief Problem Checklist. Mediation was tested using the SPSS PROCESS macro. Analyses suggested that the effect of adapted CBT on anxiety was mediated by its effects on social responsiveness, with a statistically significant indirect effect. Youth randomly assigned to adapted CBT exhibited better overall mental health at posttreatment compared to those randomized to the other conditions, and this effect was also mediated by improved social responsiveness. CBT adapted to address some of the social needs of autistic youth may enhance mental health outcomes by supporting social responsiveness, perhaps increasing the ease and effectiveness with which some youth can navigate potentially stressful situations such as entering and participating in group activities.
INTRODUCTION: Recent encouraging outcomes with neoadjuvant immune checkpoint inhibitors (ICIs) in mucosal head and neck squamous cell carcinoma (HNSCC) have generated interest in surgical de-escalation. However, the oncologic safety of response-adapted surgery (RAS) and its ability to achieve survival outcomes comparable to baseline-planned surgery (BPS) remain uncertain. METHODS: A systematic search of the PubMed, EMBASE, Cochrane Library, and the Clinical Trials Registry for studies of neoadjuvant ICIs, with or without chemotherapy, in resectable mucosal HNSCC, that explicitly report surgical extent, between 2020-2025 was performed. Two independent reviewers extracted data following PRISMA guidelines. Main outcomes included major pathologic response (MPR), pathologic complete response (pCR), event-free survival (EFS), and overall survival (OS). Study-level proportions were pooled by random effects models. Heterogeneity was assessed by the I2 statistic. RESULTS: The comparative analysis consisted of 4 RAS studies (involving 202 patients) and 11 BPS studies (403 patients). The pooled overall EFS was 83.3% (76.9-88.2) for the former and 82% (75.2-87.2) for the latter (P=.751), and the respective pooled OS was 92.3% (87.5-95.3) and 91.4% (80.3-96.5) (P=.839). The pooled pCR rate was 41.7% (95% CI 5.4-48.4; I2=.0) for RAS and 19.8% (95%CI 13.3-29.6; I2=.62) for BPS (P=.001), while the MPR was not significantly different (59.6%, versus 48.5%, P=.245). RAS was associated with greater organ preservation and reduced need for mandibulectomy and free-flap reconstruction. CONCLUSIONS: RAS following neoadjuvant ICIs in mucosal HNSCC may enable surgical de-escalation with preserved oncologic outcomes and improved function in selected patients. Larger prospective studies are warranted.
BACKGROUND: Aging is associated with progressive declines in skeletal muscle mass (SMM) and substantial interindividual variability in adaptations to resistance training (RT). However, whether training load intensity influences hypertrophic responsiveness in older adults remains unclear. OBJECTIVE: To investigate the effects of two commonly prescribed RT loading schemes (15RM and 10RM) on estimated SMM adaptations and responsiveness in older women. METHODS: Twenty-seven older women (67.7 ± 4.9 years; 65.7 ± 10.6 kg; 154.5 ± 5.8 cm) completed two 8-week RT interventions performed at 15RM and 10RM in a randomized crossover design, separated by a 12-week detraining period. Estimated SMM was assessed using dual-energy X-ray absorptiometry. Responsiveness was operationally defined using the standard error of measurement, and participants were classified as responders (RP) or non-responders (N-RP) according to changes in estimated SMM. RESULTS: Both loading schemes promoted significant and comparable increases in estimated SMM (15RM: +0.62 ± 0.4 kg vs. 10RM: +0.43 ± 0.3 kg). Overall, 61.1% of participants were classified as responders in at least one condition, whereas 38.9% were classified as non-responders. Responsiveness status was consistent across conditions for most participants, with 48.1% classified as responders under both loading schemes and 25.9% as non-responders under both conditions. Only seven participants (25.9%) changed responsiveness status between interventions. No clear pattern indicated a greater likelihood of responsiveness with either loading scheme. CONCLUSIONS: Results suggest that both the 10RM and 15RM RT protocols significantly increase estimated SMM in older women, with no significant differences between training conditions. Although substantial interindividual variability in adaptive responses exists, training load intensity does not appear to meaningfully influence responsiveness status in most individuals. These findings suggest that both loading schemes can be effectively prescribed to promote SMM gains in older women.
We aimed to compare exercise modalities and evaluate dose-response relationships with glycemic control including continuous aerobic exercise (CAE), resistance training (RT), combined exercise (CE), mind-body exercise (MBE), and high-intensity interval training (HIIT) in older adults with type 2 diabetes mellitus (T2DM). Three databases were searched for randomized controlled trials of exercise interventions in older adults with T2DM reporting glycated hemoglobin (HbA1c). Pairwise, Bayesian network, and dose-response meta-analyses were conducted. Compared with control, HIIT demonstrated the largest estimated reduction (MD = -0.95%; 95% CrI -1.45, -0.49), followed by CE (MD = -0.59%; 95% CrI -0.93, -0.25), CAE (MD = -0.46%; 95% CrI -0.69, -0.24), MBE (MD = -0.42%; 95% CrI -0.76, -0.10), and RT (MD = -0.29%; 95% CrI -0.51, -0.08). Dose-response network meta-analyses suggested a non-linear association between overall exercise dose and HbA1c reduction, with maximal estimated benefits at approximately 704 METs-min/week with the 95% CrI excluding zero between 241 and 920 METs-min/week. HIIT demonstrated the steepest estimated dose-response relationship, but with wider credible intervals. Other exercise modalities showed more gradual dose-response patterns across their estimated effective ranges. Our findings suggest that exercise prescription for older adults with T2DM should be individualized according to exercise modality, dose, and health status.
Golden pompano (Trachinotus ovatus) is one of the most economically important marine fish species in China. It is susceptible to low-temperature stress, which significantly challenges its production and supply. Nevertheless, study on the regulatory mechanisms underlying low-temperature stress responses in golden pompano remains limited. Here, we firstly performed a time-series transcriptome analysis to reconstruct dynamic response patterns under low-temperature stress in golden pompano. Transcriptome profiling identified common differentially expressed genes (DEGs), including fos, hlf, and hmgb1, as well as condition-specific DEGs across distinct low-temperature stress groups. Based on cluster analysis, all DEGs were classified into five distinct expression patterns, reflecting diversified regulation of expression in golden pompano during low-temperature stress. Furthermore, condition-specific regulatory modules were explored via weighted gene co-expression network analysis (WGCNA), highlighting that the two module hub genes, serbf2 and lipc, might respond to low-temperature stress by regulating the lipid catabolic process. Subsequently, untargeted metabolomic analysis revealed that glycerophospholipid metabolism was a significantly enriched common pathway, highlighting its crucial role in mediating the response to low-temperature stress. Finally, by integrating transcriptomic and metabolomic analyses, a gene-metabolite interaction network associated with glycerophospholipid metabolism under low-temperature stress was established. These findings underscore the significance of multiple candidate genes and glycerophospholipid metabolism in golden pompano's response to low-temperature stress, thereby laying a solid molecular foundation for the development of low-temperature-tolerant fish strains.
Inhaled antibiotics and DPP-1 inhibitors improve clinical outcomes in bronchiectasis, but whether baseline symptom burden predicts differential treatment responses remains unclear. In this post hoc analysis of two multicenter randomized trials (SAVE-BE, n = 224; TORNASOL, n = 357), we evaluate the association between baseline Quality of Life-Bronchiectasis Respiratory Symptom Scale (QoL-B-RSS) and treatment effects of florensocatib and inhaled tobramycin. In SAVE-BE, florensocatib reduces exacerbation rates versus placebo (relative risk [RR], 0.47; 95% confidence interval [CI], 0.33-0.67; p < 0.0001), with RRs of 0.53 and 0.40 observed in patients with high and low symptom burdens, respectively, but no significant symptomatic improvement. In TORNASOL, tobramycin produces clinically meaningful QoL-B-RSS improvements (exceeding the 8-point cutoff in high-symptom patients) and ameliorates bronchitic symptoms, with greater benefits in those with higher baseline symptom burden. These hypothesis-generating findings suggest that baseline symptom burden may identify differential responses to anti-inflammatory versus anti-infective therapies in bronchiectasis and support its potential as a simple, practical stratification tool to guide personalized treatment.
6PPD-quinone (6PPDQ), an oxidative transformation product of the widely used tire antioxidant 6PPD, is a ubiquitous environmental contaminant with bioaccumulation potential and widespread human exposure. Recent epidemiological evidence indicates a positive association between urinary 6PPDQ levels and colorectal cancer (CRC) risk; however, its biological effects on CRC-related phenotypes and associated molecular responses remain unclear. We integrated bioinformatics analysis, prognostic modeling, molecular docking and dynamics simulations, and in vitro experiments to investigate cellular and molecular responses to 6PPDQ in CRC models. Predicted 6PPDQ targets were intersected with CRC prognosis-related genes from The Cancer Genome Atlas, followed by functional enrichment and LASSO regression to construct a prognostic risk model, with 1-, 3-, and 5-year AUC values of 0.727, 0.754, and 0.778, respectively. Molecular docking and 100-ns molecular dynamics simulations suggested interactions between 6PPDQ and candidate proteins, including CPT2, SHC2, SRMS, and STK35. Functional assays showed that 6PPDQ exposure altered proliferation, wound-closure capacity, and invasion in Caco-2 and HCT116 cells across the nanomolar concentration range, with non-monotonic and cell-line-dependent responses. In contrast, NCM460 cells showed no increase in EdU incorporation at 10 or 100 nM, whereas reduced proliferation at higher concentrations was accompanied by increased LDH release. 6PPDQ also altered the expression of several prognosis-associated candidate genes. These findings identify cellular phenotypes and candidate molecular responses associated with 6PPDQ exposure under the tested in vitro conditions, but do not establish their causal roles or in vivo relevance. Further mechanistic and in vivo studies are required.
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.
The EARLY 2 FACTOR (E2F) transcription factor (TF) family plays a pivotal role in regulating plant development and adaptations to environmental stresses. However, the physiological function of E2Fs in pepper (Capsicum annuum L.) are not well elucidated. In this work, we conduct a comprehensive genome-wide annotation of the E2F family within the Zunla-1 pepper genome and further explore the biological roles of CaDPb in response to drought stress. Through systematic bioinformatics analysis, we identify a total of nine CaE2F genes within the Zunla-1 genome, categorizing them into three distinct subgroups. Additionally, we discover multiple cis-regulatory elements in the CaE2F promoter regions associated with responses to plant hormones and drought stress. Public RNA-seq datasets reveal distinct expression profiles of CaE2F genes across various pepper tissues and their responses to environmental stimuli and plant hormones. Subsequently, the CaDPb gene is further functionally verified in drought response. Our findings indicate that TRV2:CaDPb silenced pepper plants are more sensitivity to drought. Furthermore, we show that CaDPb participates in the regulation of reactive oxygen species (ROS) production, the expression of drought-responsive genes, and the modulation of stomatal aperture. Taken together, our findings provide a comprehensive characterization of E2F genes in pepper and offer insights into the biological function of CaDPb in pepper drought stress response.
BACKGROUND: Gestational diabetes mellitus (GDM) has a high global prevalence and arises from complex interactions between genetic predisposition and environmental factors. GDM is associated with metabolic disturbances and chronic low-grade inflammation, both of which contribute to its pathogenesis. This study aimed to investigate the association between GDM and 135 single-nucleotide polymorphisms (SNPs) across 20 genes related to metabolic traits. METHODS: In this case-control study, 152 pregnant women with GDM and 684 pregnant women with normal glucose tolerance (NGT) who underwent antenatal examination at Siriraj Hospital, Bangkok, were enrolled. Clinical data and blood samples were collected from all participants. Genomic DNA was isolated and subjected to whole-genome sequencing using the DNBSEQ-T7RS high-throughput sequencing platform. Genotype analyses were performed using R software, and haplotype analyses were conducted using the online SNPStats software. RESULTS: After adjusting for maternal age and pre-pregnancy body mass index, polymorphisms in TCF7L2 (rs34872471, rs7901695, rs4506565, rs7903146, rs12243326, and rs12255372), NEDD1 (rs10431408, rs11830756, rs249579, rs249585, and rs4762339), CMIP (rs2306115 and rs201681534), CDKAL1 (rs4710942), GCKR (rs2293572 and rs2293571), and GCK (rs5883890) were significantly associated with the risk of GDM. Haplotype analysis demonstrated that the TCF7L2 rs12243326-rs12255372 CA haplotype was associated with a decreased risk of GDM (OR = 0.44, 95% CI: 0.23-0.81), while the NEDD1 rs249579-rs249585-rs4762339 GGT haplotype was associated with an increased risk of GDM (OR = 1.40, 95% CI: 1.08-1.82). CONCLUSIONS: These findings suggest that genetic variations in TCF7L2, NEDD1, CMIP, CDKAL1, GCK, and GCKR contribute to GDM susceptibility in the Thai population.
Whether estrone (E1)-associated disruption of early reproductive development in fish is accompanied by brain responses in addition to direct gonadal effects remains unclear. Here, juvenile Siniperca chuatsi, a non-model but economically important freshwater species, were exposed for 60 d to 0, 0.01, 0.1, and 1.0 μg/L E1, spanning environmentally reported and elevated concentrations. By integrating waterborne concentration monitoring, histopathology, transcriptomics, and quantitative real-time PCR (qPCR) validation, we evaluated E1-associated changes in brain and gonadal tissues during early reproductive development. Waterborne E1 concentrations remained generally stable throughout the exposure period. At the highest tested concentration (1.0 μg/L), E1 caused neuronal vacuolation and pyknosis in the hypothalamic region and induced distinct ovarian-like structures in the gonads of genetic males. In the brain, cyp19a1, crhr1, and adcy2a were significantly upregulated, whereas egr1 was significantly downregulated, indicating transcriptional changes in genes associated with local estrogen conversion, stress-response/cAMP signaling, and neuronal activity-related regulation within a broader injury/stress-response background. In the gonad, RNA-seq analysis showed significant downregulation of star2, hsd3b1, cyp17a1, and cyp11b and significant upregulation of hsd17b1, suggesting alterations in steroidogenesis-related gene expression at the transcriptomic level. qPCR analysis of selected gonadal candidate genes showed expression directions generally consistent with the RNA-seq results, and these molecular patterns were consistent with the feminized histological phenotype. Together, these results indicate that E1 can disrupt early reproductive development in juvenile S. chuatsi and support a cautious working model in which E1 exposure is accompanied by concurrent brain and gonadal responses. This study provides new evidence for understanding the toxic effects and ecological risk implications of natural estrogen E1 during early fish development.
This research was undertaken to address environmental concerns associated with industrial solid waste and to reduce cement consumption in geotechnical engineering. It specifically investigates the feasibility of using steel slag (SS) and silica fume (SF) as partial substitutes for ordinary Portland cement (OPC) in soil stabilization. The effects of SS, SF, OPC, and initial moisture content on the unconfined compressive strength (UCS) of stabilized soil were investigated through single-factor experiments and response surface methodology (RSM). The results show that SS and SF can synergistically enhance the strength of stabilized soil, although their interaction effect was not statistically significant within the investigated ranges. Compared with soil stabilized solely with OPC, the addition of 18 % SS and 10 % SF reduced OPC consumption by 3 % without compromising strength. Microstructural and compositional analyses further revealed that SS mainly supplied calcium- and silica-bearing components, while SF provided highly reactive silica and micro-filling effects, jointly promoting hydration reactions and improving the compactness of the stabilized soil matrix. As a result, more hydration products were formed in the OPC/SS/SF-stabilized soil than in the OPC-stabilized soil, which contributed to pore filling and strength enhancement. This study provides useful guidance for the sustainable utilization of industrial solid waste and the low-carbon development of soil stabilization materials.
The expansion of aquaculture into inland saline-alkali waters is constrained by high carbonate alkalinity (CA), a severe environmental stressor for crustaceans. However, the systemic molecular mechanisms underlying its lethal toxicity remain poorly understood. In this study, we employed a comparative transcriptomic approach to investigate the tissue-specific responses of Pacific white shrimp, Penaeus vannamei, under acute lethal stress (48 h-LC50). We focused on three functionally distinct organs: the gill, hepatopancreas, and eyestalk. The results revealed a systemic but highly tissue-specific transcriptomic response. The gill, as the primary interface, exhibited severe structural impairment and critical failure of osmoregulation, highlighted by the significant downregulation of delta-1-pyrroline-5-carboxylate synthetase (P5CS). In contrast, the hepatopancreas initiates a profound metabolic trade-off, sacrificing growth-related pathways to bolster a robust antioxidant defense system, as evidenced by the activation of sulfur metabolism and high protein turnover. The eyestalk displayed a striking disconnect between hyperactivated stress signaling pathways (e.g., mTOR/FoxO) and the collapse of its protein secretory machinery, marked by the suppression of the ER translocon component Sec61. Collectively, our findings suggest that lethal alkalinity toxicity in P. vannamei results from systemic collapse driven by a complex interplay of osmoregulatory failure, metabolic trade-offs, and endocrine disruption. This study provides a comprehensive molecular snapshot of an organism at its physiological limit, offering novel insights into the adaptive strategies and ultimate tolerance boundaries of crustaceans in extreme environments.
Long non-coding RNAs (lncRNAs) stand as newly-arisen molecular types that exert regulatory effects, able to operate as competitive endogenous RNAs (ceRNAs) to engage microRNAs (miRNAs) in interaction, resulting in the recovery of target mRNA expression and activity. Increasing evidences indicate that the ceRNA network affects various biological processes in mammals, including development, cellular differentiation, metabolism, immune response, and disease pathogenesis. In teleost fish, the lncRNA-miRNA-mRNA regulatory networks have been reported occasionally. However, up to now, the roles of lncRNAs in the big-belly seahorse (Hippocampus abdominalis) remains unclear. In this study, we reported for the first time, via whole-transcriptome RNA sequencing, the lncRNA mediated ceRNA regulatory network in Vibrio harveyi-infected H. abdominalis. A total of 4197 differentially expressed mRNAs (DE-mRNAs), 1317 DE-lncRNAs, and 183 DE-miRNAs were identified. Furthermore, the crosstalk between miRNAs and lncRNAs as well as between miRNAs and mRNAs was inferred based on the negative correlations between miRNAs and their target lncRNAs/mRNAs. A core immune associated lncRNA-miRNA-mRNA putative regulatory network was thus constructed, comprising 211 lncRNA-miRNA and 224 mRNA-miRNA pairs. In conclusion, our findings provide an integrative overview of the ceRNA regulatory networks on the underlying immune responses to V. harveyi infection in the big-belly seahorse, and offer a solid theoretical foundation for the comparative immunological research of teleost fish.