SEARCH · PubMed Health

Results for “Comparative transcriptomics”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

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

At least 253 records · Page 14Linked to original sources

PubMed2006

Microarrays reveal distinct gene signatures in the thymus of seropositive and seronegative myasthenia gravis patients and the role of CC chemokine ligand 21 in thymic hyperplasia.

Myasthenia gravis (MG) is an autoimmune disease mainly caused by antiacetylcholine receptor autoantibodies (seropositive (SP) disease) or by Abs against unknown autoantigenic target(s) (seronegative (SN) disease). Thymectomy is usually beneficial although thymic hyperplasia with ectopic germinal centers is mainly observed in SP MG. To understand the role of thymus in the disease process, we compared the thymic transcriptome of non-MG adults to those of SP patients with a low or high degree of hyperplasia or SN patients. Surprisingly, an overexpression of MHC class II, Ig, and B cell marker genes is observed in SP but also SN MG patients. Moreover, we demonstrate an overexpression of CXCL13 in all MG thymuses leading probably to the generalized B cell infiltration. However, we find different chemotactic properties for MG subgroups and, especially, a specific overexpression of CCL21 in hyperplastic thymuses triggering most likely ectopic germinal center development. Besides, SN patients present a peculiar signature with an abnormal expression of genes involved in muscle development and synaptic transmission, but also genes implicated in host response, suggesting that viral infection might be related to SN MG. Altogether, these results underline differential pathogenic mechanisms in the thymus of SP and SN MG and propose new research areas.

Adult↗

PubMed2026

T2T Genome Assembly and Multi-Omics Data Reveal Terrestrial Adaptation and Mucus Biosynthesis in Tropical Leatherleaf Slug (Laevicaulis alte).

Laevichaulis alte is a slug in the order Systellommatophora that evolved from aquatic ancestors and now faces strong challenges from desiccation, respiration on land, and novel pathogens. Its mucus is essential for water retention, locomotion, and defense. To link terrestrial adaptation with mucus biosynthesis, we generated a gap-free genome assembly of L. alte using PacBio HiFi reads, Oxford Nanopore ultra-long reads, and Hi-C data. The genome shows low heterozygosity and holocentromeric chromosomes. Functional metabolomics revealed marked metabolic shifts between L. alte and the closely related aquatic species Peronia verruculata. In L. alte, differential metabolites were enriched in lipid metabolism, immune regulation, and stress response pathways, consistent with life in a dry and microbe-rich terrestrial environment. Comparative genomics and transcriptomics identified candidate genes linked to mucus secretion and physiological adaptation, including VEGF, ASGR2, and COL6A6. Further analyses highlighted the vascular endothelial growth factor (VEGF) gene family as a key regulator connecting angiogenesis, tissue remodeling, and mucus production pathways in L. alte. Together, this gap-free genome and multi-omics dataset establish a molecular framework that links genomic innovation, mucus biology, and terrestrial adaptation in Systellommatophora, and they offer a basis for understanding ecological niche specialization in land molluscs.

Animals↗

PubMed2002

Bacillus subtilis functional genomics: genome-wide analysis of the DegS-DegU regulon by transcriptomics and proteomics.

The DegS-DegU two-component regulatory system of Bacillus subtilis controls various processes that characterize the transition from the exponential to the stationary growth phase, including the induction of extracellular degradative enzymes, expression of late competence genes and down-regulation of the sigma(D) regulon. The degU32(Hy) mutation stabilizes the phosphorylated form of DegU (DegU-P), resulting in overproduction of several extracellular degradative enzymes. In this study, the pleiotropic DegS-DegU regulon was characterized by combining proteomic and transcriptomic approaches. A comparative analysis of wild-type B. subtilis and the degU32(Hy) mutant grown in complex medium was performed during the exponential and in the stationary growth phase. Besides genes already known to be under the control of DegU-P, novel putative members of this regulon were identified. Although the degU32(Hy) mutant is assumed to contain high levels of phosphorylated DegU in the exponential as well as in the stationary growth phase, many genes known to be positively regulated by DegU-P did not show enhanced expression in the mutant strain during exponential growth. This is consistent with the fact that most genes belonging to the DegS-DegU regulon are subject to multiple regulation; this is also reflected in the strong stationary-phase induction of these genes in the mutant strain. As expected, during the exponential growth phase, the sigma(D) regulon was expressed at significantly lower levels in the degU32(Hy) mutant than in the wild type.

Bacillus subtilis↗

PubMed2006

Bdellovibrio: growth and development during the predatory cycle.

Predatory Bdellovibrio enter the periplasm of other Gram-negative bacteria, growing within and consuming them. Unravelling molecular details of this intimate association between bacterial predator and prey is challenging yet fascinating, and might lead to novel antibacterials in the future. Pioneering physiological and biochemical studies described the predatory life of Bdellovibrio in the 1960s and 1970s, later followed by recombinant DNA work in the 1990s, which led to a revival in Bdellovibrio molecular research. This revival continues in the 21st century with the advent of a genome sequence. Now worldwide research is underway on the comparative genomics and transcriptomics of predatory bacteria, and will illuminate the evolutionary adaptations to become predatory, and will hopefully ultimately illuminate how the predatory processes of Bdellovibrio can be employed against pathogenic bacteria and for humankind.

Bacterial Adhesion↗

PubMed2026

Time-dependent effects of rapid-acting antidepressants in iPSC-derived neurons from treatment-resistant depression and healthy volunteers.

Rapid-acting antidepressants like ketamine and serotonergic psychedelics show promise for treatment-resistant depression (TRD), but the molecular mechanisms that contribute to their therapeutic effects remain unclear. Induced pluripotent stem cells (iPSCs) offer a platform to model human cortical neurons and investigate drug effects in a human-relevant system. Here, iPSCs from individuals with TRD and healthy volunteers (HVs) were differentiated into mature cortical-like neurons and treated for six and 24 h with agents being investigated as rapid-acting antidepressants, including (2 R,6 R)-hydroxynorketamine (HNK), psilocybin, lysergic acid diethylamide (LSD), and 2,5-Dimethoxy-4-iodoamphetamine (DOI). Bulk and single-cell RNA sequencing assessed global and cell-type-specific transcriptomic responses. Synaptic proteins were evaluated via Western blotting and immunocytochemistry. To validate translational relevance, transcriptomic results were compared to CSF proteomics from ketamine-treated HVs. Despite differing initial pharmacological targets, overall gene expression across all compounds was highly correlated at matched timepoints compared to vehicle control, suggesting shared downstream effects. Both glutamatergic and serotonergic drugs converged on pathways involving inflammation, mTORC1 signaling, and cellular growth. At the single-cell level, (2 R,6 R)-HNK showed distinct cell-type specific alterations: upregulation in excitatory neurons and concomitant downregulation of inhibitory neuron populations. Differentially expressed genes from (2 R,6 R)-HNK-treated neurons also overlapped with CSF proteomic signatures from ketamine-treated individuals, supporting the model's translational relevance. This study is the first to assess multiple putative rapid-acting antidepressants in parallel using an iPSC-derived neuron model. Both convergent and drug-specific changes in gene expression and pathway enrichment were observed across diverse compounds, supporting the use of human iPSC-derived neurons in antidepressant drug discovery. Clinical Trial Registry: www.clinical trials.gov, NCT02484456.

Journal Article↗

PubMed2003

Redirection of B cell responsiveness by transforming growth factor beta receptor.

The multifunctional transforming growth factor beta receptor (TbetaR) ligand pair plays a central role in the regulation of lymphocyte homeostasis and prevention of autoimmunity. Although the mechanisms underlying the induction of transcriptional modulators by TbetaR have been studied in considerable detail, relatively little is known about the regulatory pathways targeted. To shed light on the mechanisms involved in negative regulation of B cell responses we identified TbetaR-dependent transcriptome changes by comparative gene expression profiling of normal and TbetaR-deficient primary B cells. The data reveal TbetaR-mediated induction of inhibitors of antigen receptor signaling (Ship-1, CD72) as well as inhibitors of the Jak/Stat pathway and signaling by means of Toll-like receptors (SOCS1,3). These inhibitory effects are complemented by induction of antiproliferative transcription factors. In contrast to this inhibition, G protein-coupled receptors such as CXCR4 and agonists mediating Ca2+ flux (inositol trisphosphate receptor subtype 2) are induced by TbetaR, indicating enhancement of the Ca2+ storage/ release system and chemotactic responses. Suppression of proapoptotic genes suggests support of cell survival. Confirming the shift in B cell responsiveness, antigen-receptor-mediated activation of Syk and phospholipase C-gamma2, as well as Stat6 phosphorylation, is inhibited, whereas chemotaxis, Ca2+ release, and cell survival are enhanced in transforming growth factor-beta-sensitive B cells. The data provide a molecular basis for TbetaR-mediated inhibition of B cell responsiveness and indicate that TbetaR maintains homeostasis not only through inhibition of the cell cycle but also by delivering a coherent instructive signal that redirects responsiveness to microenvironmental cues.

Animals↗

PubMed2026

Co-option of ancestral stem regulators drove recurrent evolution of underground storage organs.

Geophytes are plants that produce underground storage organs such as tubers, rhizomes, and bulbs, to facilitate asexual reproduction and withstand a myriad of environmental challenges. While the potato (Solanum tuberosum L.) serves as the primary model for studying tuberization, the genetic mechanisms encoding this trait across diverse angiosperm lineages remains unclear. This study utilized a phylogenomic-transcriptomic approach to compare tuber development across nine tuberizing species with five nontuberizing sister taxa. We identified orthologs of key potato tuberization genes that exhibit similar expression in the stolons or tubers of these distant relatives. In nontuberizing species, these orthologs exhibit distinct expression profiles and are primarily expressed in the stem. This suggests that the independent evolution of tubers across angiosperms resulted from shifts in the expression of preexisting genes that led to their co-option. This process, also known as exaptation, occurs when existing genetic suites are recruited for entirely new biological functions. This mechanism stands in contrast to the repeated loss or gain of genes, which has been associated with the origin of other adaptive plant traits. Furthermore, the co-option of the same genes was observed in species with other stem-derived storage organs, such as rhizomes and runners. These findings reveal a conserved evolutionary model for the development of stem-derived geophyte organs that evolved independently across the flowering plants over the past 160 My.

Plant Tubers↗

PubMed2026

KCFtools: rapid alignment-free method for introgression screening and GWAS using k-mer profiles.

MOTIVATION: In the era of multiple genome references, researchers often align sequencing reads against distinct assemblies or even multiple references simultaneously. This enables applications such as the detection of introgressed segments or highly variable genomic regions, which are especially prevalent in large-genome crop species such as lettuce or wheat. However, these applications come at the cost of increased computational burden, inconsistencies in mapping methods, and reduced reproducibility across studies. To address these limitations, we developed KCFtools, a Java-based toolkit that identifies the presence and absence of k-mers in nonoverlapping genomic or transcriptomic windows by comparing query and reference genomes. This alignment-free approach enables the efficient computation of an identity score for each window, thereby facilitating robust detection of introgressed or variable regions across genomes. RESULTS: We systematically evaluated the performance and accuracy of the k-mer-based method implemented in KCFtools, benchmarking it against conventional single nucleotide variation-based introgression detection pipelines. Our results demonstrate that KCFtools effectively captures introgressed segments and structurally diverse regions, even in species with fragmented or highly divergent reference genomes. In addition, we extended KCFtools to generate genotype matrices from k-mer variation tables. These matrices are compatible with genome-wide association studies software and allow the identification of loci associated with phenotypic traits. We showcase the utility of this approach by detecting known and novel associations for downy mildew resistance in lettuce, underscoring the pipeline's potential for high-resolution, reference-agnostic population genetic analysis. AVAILABILITY AND IMPLEMENTATION: https://github.com/sivasubramanics/kcftools.

Software↗

PubMed2004

Dietary protein-related changes in hepatic transcription correspond to modifications in hepatic protein expression in growing pigs.

In a previous investigation we showed by expression profiling based on transcription analysis using differential display RT-PCR (DDRT-PCR) and real-time RT-PCR that a soy protein diet (SPI) significantly changes the hepatic transcription pattern compared with a casein diet (CAS). The present study was conducted to determine whether the transcriptional modulation is translated into protein expression. The hepatic mRNA abundance of four genes (EP24.16, LC3, NPAP60L, RFC2) that showed diet-related expression in previous DDRT-PCR experiments was analyzed by real-time RT-PCR. Two pigs that showed the most prominent SPI-related changes of transcription and two casein-fed pigs were selected and their hepatic protein pattern was studied comparatively by two-dimensional gel electrophoresis and peptide mass fingerprinting. The two-dimensional protein gel electrophoresis revealed a predominant SPI-associated upregulation of protein expression that corresponded to the results of the mRNA study. Of 380 diet-related protein spots displayed, 215 appeared exclusively or enlarged in the two SPI pigs; 10 of 39 diet-related expressed protein spots extracted could be identified by peptide mass fingerprinting and database search. Compared with the transcriptomics approach, the proteomics approach led in part to the identification of the same diet-associated expressed molecules (plasminogen, trypsin, phospholipase A2, glutathione-S-transferase alpha, retinal binding protein) or at least molecules belonging to the same metabolic pathways (protein and amino acid metabolism, oxidative stress response, lipid metabolism). The present results at the proteome level confirm SPI-related increased oxidative stress response and significant effects on protein biosynthesis already observed at the transcriptome level.

Animals↗

PubMed2001

Gene expression profiling using a novel method: amplified differential gene expression (ADGE).

Amplified differential gene expression (ADGE) is a novel technique, designed to profile gene expression of the whole transcriptome or to compare expression of a set of genes between two samples. ADGE employs hybridization to quadratically amplify the ratio of an expressed gene between control and tester samples before displaying. The subtle structures of adapters and primers are designed for displaying the amplified ratio of an expressed gene between two samples. Four selective nucleotides at the 3' end of primers are used to increase PCR efficiency for targeted molecules and to improve detection of PCR products. Double PCR with the same pair of primers expands the detection range, especially for genes of low abundance. Integration of these steps makes ADGE sensitive and accurate. Application to drug resistant human tumor cell lines showed that ADGE accurately profiled expression levels for induced, repressed or unchanged genes. The qualitative expression patterns for ADGE were verified with RT-PCR.

DNA Primers↗

PubMed2002

Density of points clustering, application to transcriptomic data analysis.

With the increasing amount of data produced by high-throughput technologies in many fields of science, clustering has become an integral step in exploratory data analysis in order to group similar elements into classes. However, many clustering algorithms can only work properly if aided by human expertise. For example, one parameter which is crucial and often manually set is the number of clusters present in the analyzed set. We present a novel stopping rule to find the optimal number of clusters based on the comparison of the density of points inside the clusters and between them. The method is evaluated on synthetic as well as on real transcriptomic data and compared with two current methods. Finally, we illustrate its usefulness in the analysis of the expression profiles of promyelocytic cells before and after treatment with all-trans retinoic acid. Simultaneous clustering for gene regulation and absolute initial expression levels allowed the identification of numerous genes associated with signal transduction revealing the complexity of retinoic acid signaling.

Algorithms↗

PubMed2005

Benchmarking the CATMA microarray. A novel tool for Arabidopsis transcriptome analysis.

Transcript profiling is crucial to study biological systems, and various platforms have been implemented to survey mRNAs at the genome scale. We have assessed the performance of the CATMA microarray designed for Arabidopsis (Arabidopsis thaliana) transcriptome analysis and compared it with the Agilent and Affymetrix commercial platforms. The CATMA array consists of gene-specific sequence tags of 150 to 500 bp, the Agilent (Arabidopsis 2) array of 60mer oligonucleotides, and the Affymetrix gene chip (ATH1) of 25mer oligonucleotide sets. We have matched each probe repertoire with the Arabidopsis genome annotation (The Institute for Genomic Research release 5.0) and determined the correspondence between them. Array performance was analyzed by hybridization with labeled targets derived from eight RNA samples made of shoot total RNA spiked with a calibrated series of 14 control transcripts. CATMA arrays showed the largest dynamic range extending over three to four logs. Agilent and Affymetrix arrays displayed a narrower range, presumably because signal saturation occurred for transcripts at concentrations beyond 1,000 copies per cell. Sensitivity was comparable for all three platforms. For Affymetrix GeneChip data, the RMA software package outperformed Microarray Suite 5.0 for all investigated criteria, confirming that the information provided by the mismatch oligonucleotides has no added value. In addition, taking advantage of replicates in our dataset, we conducted a robust statistical analysis of the platform propensity to yield false positive and false negative differentially expressed genes, and all gave satisfactory results. The results establish the CATMA array as a mature alternative to the Affymetrix and Agilent platforms.

Arabidopsis↗

PubMed2025

Bundle sheath cell-specific expression of chloroplast genes encoding subunits of the NADH dehydrogenase-like complex in maize.

C4 photosynthesis alleviates the limitation caused by the oxygenase activity of Rubisco by partitioning photosynthetic functions between two distinct cell types: bundle sheath cells (BSCs) and mesophyll cells (MCs). These cell types perform different steps of photosynthesis using specialized machinery, accompanied by differential expression of chloroplast genes. To uncover the underlying molecular mechanisms for this differentiation, we isolated BSCs and MCs and compared their chloroplast transcriptomes, focusing on the chloroplast NADH dehydrogenase-like (NDH) complex, which is enriched in BSCs. To investigate whether RNA stabilization contributes to differential gene expression, we analyzed RNA footprints that reflect the binding of pentatricopeptide repeat (PPR) proteins to their RNA targets. We could not detect cell-type-specific accumulation of footprint RNAs. We then focused on transcriptional regulation, specifically on an operon that starts with the rps15 gene. The operon includes six ndh genes and the psaC gene encoding a photosystem I subunit. Transcript levels of all genes in this operon were higher in BSCs than in MCs, suggesting coordinated regulation as a transcriptional unit. Based on the genomic location of the rps15 gene within inverted repeats near the junctions on both sides of the small single copy region, we demonstrated that rps15, through two distinct promoters, is sufficient to drive preferential accumulation of downstream transcripts in BSCs.

Zea mays↗

PubMed2006

In vivo dissection of the Helicobacter pylori Fur regulatory circuit by genome-wide location analysis.

Iron homeostasis is particularly important in pathogenic bacteria, which need to compete with the host for this essential cofactor. In Helicobacter pylori, a causative agent of several gastric pathologies, iron uptake and storage genes are regulated at the transcriptional level by the ferric uptake regulator Fur. The regulatory circuit of Fur has recently come under focus because of an intimate interlink with a broader regulatory network governing metal homeostasis, acidic response, and virulence. To dissect the Fur regulatory circuit and identify in vivo targets of regulation, we developed a genome-wide location analysis protocol which allowed the identification of 200 genomic loci bound by Fur as well as the investigation of the binding efficiency of the protein to these loci in response to iron. Comparative analysis with transcriptomes of wild-type and fur deletion mutant strains allowed the distinction between targets associated with Fur regulation and genes indirectly influenced by the fur mutation. The Fur regulon includes 59 genes, 25 of which appear to be positively regulated. A case study conducted by primer extension analysis of two oppositely regulated genes, hpn2 and flaB, suggests that negative regulation as well as positive regulation occurs at the transcriptional level. Furthermore, the results revealed the existence of 13 Fur targeted loci within polycistronic operons, which were associated with transcript deregulation in the fur mutant strain. This study provides a systematic insight of Fur regulation at the genome-wide level in H. pylori and points to regulatory functions extending beyond the classical Fur repression paradigm.

Bacterial Proteins↗

PubMed2005

Integrative genomic, transcriptional, and proteomic diversity in natural isolates of the human pathogen Burkholderia pseudomallei.

Natural isolates of pathogenic bacteria can exhibit a broad range of phenotypic traits. To investigate the molecular mechanisms contributing to such phenotypic variability, we compared the genomes, transcriptomes, and proteomes of two natural isolates of the gram-negative bacterium Burkholderia pseudomallei, the causative agent of the human disease melioidosis. Significant intrinsic genomic, transcriptional, and proteomic variations were observed between the two strains involving genes of diverse functions. We identified 16 strain-specific regions in the B. pseudomallei K96243 reference genome, and for eight regions their differential presence could be ascribed to either DNA acquisition or loss. A remarkable 43% of the transcriptional differences between the strains could be attributed to genes that were differentially present between K96243 and Bp15682, demonstrating the importance of lateral gene transfer or gene loss events in contributing to pathogen diversity at the gene expression level. Proteins expressed in a strain-specific manner were similarly correlated at the gene expression level, but up to 38% of the global proteomic variation between strains comprised proteins expressed in both strains but associated with strain-specific protein isoforms. Collectively, >65 hypothetical genes were transcriptionally or proteomically expressed, supporting their bona fide biological presence. Our results provide, for the first time, an integrated framework for classifying the repertoire of natural variations existing at distinct molecular levels for an important human pathogen.

Bacterial Proteins↗

PubMed2015

Advanced age is associated with worsened outcomes and a unique genomic response in severely injured patients with hemorrhagic shock.

INTRODUCTION: We wished to characterize the relationship of advanced age to clinical outcomes and to transcriptomic responses after severe blunt traumatic injury with hemorrhagic shock. METHODS: We performed epidemiological, cytokine, and transcriptomic analyses on a prospective, multi-center cohort of 1,928 severely injured patients. RESULTS: We found that there was no difference in injury severity between the aged (age ≥55, n = 533) and young (age <55, n = 1395) cohorts. However, aged patients had more comorbidities. Advanced age was associated with more severe organ failure, infectious complications, ventilator days, and intensive care unit length of stay, as well as, an increased likelihood of being discharged to skilled nursing or long-term care facilities. Additionally, advanced age was an independent predictor of a complicated recovery and 28-day mortality. Acutely after trauma, blood neutrophil genome-wide expression analysis revealed an attenuated transcriptomic response as compared to the young; this attenuated response was supported by the patients' plasma cytokine and chemokine concentrations. Later, these patients demonstrated gene expression changes consistent with simultaneous, persistent pro-inflammatory and immunosuppressive states. CONCLUSIONS: We concluded that advanced age is one of the strongest non-injury related risk factors for poor outcomes after severe trauma with hemorrhagic shock and is associated with an altered and unique peripheral leukocyte genomic response. As the general population's age increases, it will be important to individualize prediction models and therapeutic targets to this high risk cohort.

Adult↗

PubMed2026

Time-Dependent Effects of Rapid-Acting Antidepressants in iPSC-Derived Neurons from Treatment-Resistant Depression and Healthy Volunteers.

UNLABELLED: Rapid-acting antidepressants like ketamine and serotonergic psychedelics show promise for treatment-resistant depression (TRD), but the molecular mechanisms that contribute to their therapeutic effects remain unclear. Induced pluripotent stem cells (iPSCs) offer a platform to model human cortical neurons and investigate drug effects in a human-relevant system. Here, iPSCs from individuals with TRD and healthy volunteers (HVs) were differentiated into mature cortical-like neurons and treated for six and 24 hours with agents being investigated as rapid-acting antidepressants, including (2R,6R)-hydroxynorketamine (HNK), psilocybin, lysergic acid diethylamide (LSD), and 2,5-Dimethoxy-4-iodoamphetamine (DOI). Bulk and single-cell RNA sequencing assessed global and cell-type-specific transcriptomic responses. Synaptic proteins were evaluated via Western blotting and immunocytochemistry. To validate translational relevance, transcriptomic results were compared to CSF proteomics from ketamine-treated HVs. Despite differing initial pharmacological targets, overall gene expression across all compounds was highly correlated at matched timepoints compared to vehicle control, suggesting shared downstream effects. Both glutamatergic and serotonergic drugs converged on pathways involving inflammation, mTORC1 signaling, and cellular growth. At the single-cell level, HNK showed distinct cell-type specific alterations: upregulation in excitatory neurons and concomitant downregulation of inhibitory neuron populations. Differentially expressed genes from HNK-treated neurons also overlapped with CSF proteomic signatures from ketamine-treated individuals, supporting the model's translational relevance. This study is the first to assess multiple putative rapid-acting antidepressants in parallel using an iPSC-derived neuron model. Both convergent and drug-specific changes in gene expression and pathway enrichment were observed across diverse compounds, supporting the use of human iPSC-derived neurons in antidepressant drug discovery. CLINICAL TRIAL REGISTRY: www.clinicaltrials.gov, NCT02484456.

Journal Article↗

PubMed2026

Transcriptomic responses to developmental temperature in two field-collected Spodoptera exigua populations from Korea.

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.

Animals↗