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 757 records · Page 42Linked to original sources

PubMed2006

Transcriptome analysis of Aspergillus fumigatus exposed to voriconazole.

For a comprehensive evaluation of genes that have their expression modulated during exposure of the mycelia to voriconazole, we performed a large-scale analysis of gene expression in Aspergillus fumigatus using a microarray hybridization approach. By comparing the expression of genes between the reference time and after addition of voriconazole (30, 60, 120, and 240 min), we identified 2,271 genes differentially expressed in the wild-type strain. To validate the expression of some of these genes during exposure to voriconazole, we analyzed 13 genes showing higher expression in the presence of voriconazole by real-time RT-PCR. Although the magnitudes of induction differed between the two experimental systems, in about 85% of the cases they were in good agreement with the microarray data. To our knowledge this is the first study of microarray hybridization analysis for a filamentous fungus exposed to an antifungal agent. In our study, we have observed: (i) a decreased mRNA expression of various ergosterol biosynthesis genes; (ii) increased mRNA levels of genes involved in a variety of cell functions, such as transporters, transcription factors, proteins involved in cell metabolism, and hypothetical proteins; and (iii) the involvement of the cyclic AMP-protein kinase signaling pathway in the increased mRNA expression of several of these genes.

Antifungal Agents↗

PubMed2006

Global transcriptomic analysis of Desulfovibrio vulgaris on different electron donors.

Whole-genome microarrays of Desulfovibrio vulgaris were used to determine relative transcript levels in cells grown to exponential or stationary phase on a medium containing either lactate or formate as electron donor. The results showed that 158 and 477 genes were differentially expressed when comparing exponential to stationary phase in lactate- or formate-based media, respectively; and 505 and 355 genes were responsive to the electron donor used at exponential or stationary phase, respectively. Functional analyses suggested that the differentially regulated genes were involved in almost every aspect of cellular metabolism, with genes involved in protein synthesis, carbon, and energy metabolism being the most regulated. The results suggested that HynBA-1 might function as a primary periplasmic hydrogenase responsible for oxidation of H2 linked to the proton gradient in lactate-based medium, while several periplasmic hydrogenases including HynBA-1 and Hyd might carry out this role in formate-based medium. The results also indicated that the alcohol dehydrogenase and heterodisulfide reductase catalyzed pathway for proton gradient formation might be actively functioning for ATP synthesis in D. vulgaris. In addition, hierarchical clustering analysis using expression data across different electron donors and growth phases allowed the identification of the common electron donor independent changes in gene expression specifically associated with the exponential to stationary phase transition, and those specifically associated with the different electron donors independent of growth phase. The study provides the first global description and functional interpretation of transcriptomic response to growth phase and electron donor in D. vulgaris.

Bacterial Proteins↗

PubMed2005

Genetic expression profile of olfactory ensheathing cells is distinct from that of Schwann cells and astrocytes.

Olfactory ensheathing cells (OECs) accompany the axons of olfactory receptor neurons, which regenerate throughout life, from the olfactory mucosa into the olfactory bulb. OECs have shown widely varying efficacy in repairing the injured nervous system. Analysis of the transcriptome of OECs will help in understanding their biology and will provide tools for investigating the mechanisms of their efficacy and interactions with host tissues in lesion models. In this study, we compared the transcriptional profile of cultured OECs with that of Schwann cells (SCs) and astrocytes (ACs), two glial cell types to which OECs have similarities. Two biological replicates of RNA from cultured OECs, SCs, and ACs were hybridized to long oligo rat 5K arrays against a common reference pool of RNA (50% cultured fibroblast RNA and 50% neonatal rat brain RNA). Transcriptional profiles were analyzed by hierarchical clustering, Principal Components Analysis, and the Venn diagram. The three glial cell types had similarly increased or decreased expression of numerous transcripts compared with the reference. However, OECs were distinguishable from both SCs and ACs by a modest number of transcripts, which were significantly enriched or depleted. Furthermore, OECs and SCs were more closely related to each other than to ACs. Expression of selected transcripts not previously characterized in OECs, such as Lyz, Timp2, Gro1 (Cxcl1), Ccl2 (MCP1), Ctgf, and Cebpb, was validated by real-time reverse transcription-polymerase chain reaction (RT-PCR); immunohistochemistry in cultured OECs, SCs, and ACs, and adult tissues was performed to demonstrate their expression at the protein level.

Animals↗

PubMed2001

Molecular neuropathology of transgenic mouse models of Down syndrome.

Down syndrome (DS) is a complex, clinically heterogeneous disorder which shows both impairment of neurodevelopement and the neurodegenerative changes of Alzheimer's disease (AD). The phenotype of DS is caused by triplication of chromosome 21 and transgenic mouse models have been developed, and are being created, that carry single genes and chromosomal segments to excess. For example, transgenic mice containing additional copies of the amyloid precursor protein (APP) gene, have been useful in producing the Abeta deposition characteristic of AD and DS, but not the cytoskeletal changes that are the hallmarks of these human disorders. Such models are useful in replicating aspects of pathogenesis and allow for the testing of therapeutic agents to restore impaired function. Segmental trisomic mouse models, which survive to adulthood and possess three copies of multiple genes responsible for the DS phenotype, such as Ts1Cje and Ts65Dn, have been used to explore aspects of neurodevelopment and neurodegeneration. These animal models show some but not all the pathological, biochemical, and transcriptional changes seen in DS. They also have the advantage of allowing for the testing of therapeutic agents to restore impaired function. Analysis of the transcriptome and proteome of fetal and adult DS indicates that there is a complex relationship between gene dosage, gene and protein expression, and that data from animal models will need to be compared and evaluated in the light of data obtained from DS tissue.

Animals↗

PubMed2026

A critical appraisal of base-resolution m6A profiling techniques.

N6-methyladenosine (m6A) is the most prevalent internal modification in eukaryotic mRNA, influencing RNA fate and gene regulation. Early antibody-based approaches enabled transcriptome-wide profiling but lacked resolution and quantitative accuracy. Newer approaches now achieve base-resolution m6A detection using improved crosslinking, chemical or enzymatic conversion, and single-molecule sequencing. Antibody-free methods provide quantitative stoichiometry from minimal input, while nanopore direct RNA sequencing offers real-time, single-molecule readouts across entire transcriptomes. Collectively, these methods form a versatile toolkit that integrates global mapping with precise site-level analysis, advancing knowledge of context-dependent m6A regulation in physiology and disease. This review compares their principles, strengths, and limitations to guide method selection and highlight how next-generation epitranscriptomic tools are paving the way for clinical and therapeutic applications.

Humans↗

PubMed2026

GRNContext: an interactive web platform for contextualized gene regulatory networks visualization across human cancers.

SUMMARY: While current Gene Regulatory Network (GRN) databases provide comprehensive reference maps of potential interactions between transcription factors and target genes, they do not specify which regulatory interactions are active within specific biological contexts. This limitation is particularly critical in cancer, where transcriptional programs are inherently tissue-specific. To address this gap, we developed GRNContext, an interactive web platform designed for the visualization, exploration, and comparative analysis of gene regulatory networks contextualized across 33 cancer types from The Cancer Genome Atlas (TCGA). Our approach uses the TFLink human reference GRN as a starting point and integrates TCGA transcriptomic profiles to infer cancer-specific regulatory activity. Regulatory relevance was assessed using complementary machine learning and statistical methods, which were unified into a consensus score to prioritize and filter the most relevant candidate regulators for each target gene. By providing both curated context-specific GRNs and a user-friendly platform, GRNContext constitutes a comprehensive and accessible resource that supports mechanistic investigations, hypothesis generation, and translational research focused on transcriptional regulation in cancer. AVAILABILITY AND IMPLEMENTATION: GRNContext is supported by all major browsers and freely available on the web at https://apps.cienciavida.org/grncontext. It is implemented as a client-server web application featuring a FastAPI backend and a React frontend utilizing Cytoscape.js for interactive network visualization, all containerized via Docker for cross-platform compatibility.

Humans↗

PubMed2026

Inheritance of the epigenetic signature and reduced intermuscular bone phenotype acquired via DNA methylation editing of the runx2 b promoter in zebrafish.

The presence of intermuscular bones (IBs) can directly affect the economic value of aquaculture fish. Although genome editing can create IB-free fish by knocking out key IB-related genes, such as runx2b, the associated DNA sequence alterations raise food safety and health concerns, limiting its breeding applications. In this study, we used CRISPR/dCas9-mediated epigenome-editing technology targeting the runx2 b promoter in zebrafish to alter DNA methylation patterns without changing the DNA sequence. Our results showed that higher runx2 b promoter methylation patterns significantly inhibited eGFP mRNA expression levels in the recombinant plasmid. Using the CRISPR/dCas9-Dnmt7 system to enhance methylation of the zebrafish runx2b promoter, we observed a significant decrease in runx2 b mRNA expression levels in the F0 generation. The IBs in the 11 th-16 th muscle segments of the adult F0 fish were significantly shorter compared with the controls. Inbreeding of fish was used to produce F1 and F2 offspring that retained these high promoter methylation levels, along with persistent runx2b expression suppression and IB development inhibition. Transcriptome sequencing analysis suggested that increasing runx2 b promoter methylation levels may synergistically induce additional epigenetic modifications, potentially affecting the PPAR signaling pathway and FoxO transcription factor regulation, which appears to inhibit osteoblast proliferation and differentiation. Overall, this study demonstrates an innovative application of epigenetic editing technology for aquaculture breeding. By precisely regulating the expression patterns of key genes for economically important traits while preserving genomic DNA integrity, this approach provides a theoretical foundation and technical support for improving fish economic traits.

Animals↗

PubMed2004

Proteomic analysis as related to transcriptome data in the lung of chromium(VI)-treated rats.

Assessing the parallelism between transcriptome data and proteome data represents one of the major challenges of post-genomic research. We evaluated the levels of 380 proteins in lung S12 fractions from Sprague-Dawley rats by antibody microarrays. Approximately half of these proteins were detectable under physiological conditions. There was a poor parallelism between mRNA and protein levels for a subset of 84 coinciding or related activities, whose gene expression had previously been investigated by cDNA array. The proportion of detectable proteins was almost twice as high as the proportion of transcriptionally active genes, which reflects the longer half-life of proteins compared to mRNA. Following the local stimulus provided by a short-term and high-dose exposure to sodium dichromate by the intra-tracheal route, 64 additional proteins were detectable in lung S12 fractions, and the correlation between gene expression and protein levels became significant. Sixteen proteins were increased more than twice following chromium(VI) administration. They included ten activities involved in the positive regulation of the cell cycle, three proteins involved in stress response and protein repair, two pro-apoptotic activities, and one protein involved in lipoprotein catabolism. An increase of P53 protein was detected by Western blot in lung nuclear fractions. Post-genomic analyses, highlighting the stimulation of defence mechanisms triggered by DNA damage, contribute to explain the previously reported discrepancy between the ability of chromium(VI) to induce oxidative stress and genotoxic damage in the lung and its failure to induce lung tumors under comparable experimental conditions. The proteome analysis showed a prominent role of apoptosis, counter-balanced by a positive regulation of the cell cycle aimed at replacing lost cells. In conclusion, our results suggest that, under basal conditions, mRNA undergoes a selective inactivation and post-transcriptional regulation resulting in de-coupling of transcriptome data and proteome data. However, this parallelism is re-established when the cell undergoes genotoxic damage.

Animals↗

PubMed2006

Global transcriptome analysis of the responses of a fluoroquinolone-resistant Streptococcus pneumoniae mutant and its parent to ciprofloxacin.

Streptococcus pneumoniae M22 is a multidrug-resistant mutant selected after exposure of capsulated wild-type S. pneumoniae NCTC 7465 (strain M4) to ciprofloxacin. DNA microarray analysis comparing the gene expression profiles of strain M22 with those of strain M4 showed that strain M22 constitutively expressed 22 genes at levels higher than those observed in strain M4 under all conditions studied. These included the genes encoding the enzymes involved in branched-chain amino acid biosynthesis and two genes (patA and patB) with sequences suggestive of ABC transporter proteins. Expression of the patA and patB genes was induced by ciprofloxacin in both strains, but in strain M4 it only reached the levels observed in strain M22 after long incubation with high concentrations of ciprofloxacin. The altered expression profile observed with strain M22 suggested that the mutation or mutations acquired during resistance selection bring the cell into a state in which the expression of critical genes is preemptively altered to correct for the potential effects of ciprofloxacin on gene expression in the parent strain.

Ciprofloxacin↗

PubMed2026

Genomic, transcriptomic, and molecular predictors of response to neoadjuvant therapy in locally advanced rectal cancer: a narrative review.

Total neoadjuvant therapy (TNT) has emerged as a key treatment paradigm for locally advanced rectal cancer, reducing distant metastasis rates and facilitating organ preservation in selected patients. However, treatment response remains heterogeneous, highlighting the need for biomarkers that can guide treatment selection and optimise outcomes. This narrative review synthesises the current evidence regarding tumour-intrinsic genomic biomarkers associated with response to neoadjuvant therapy, encompassing somatic mutations, germline polymorphisms, gene expression profiles, mismatch repair (MMR) status, protein expression, epigenetic markers, and circulating tumour-derived biomarkers across conventional chemoradiotherapy (CRT) and TNT paradigms. Across the reviewed literature, individual somatic mutations, including KRAS, TP53, and BRAF, demonstrated limited reproducibility as predictive biomarkers, although KRAS mutations were recurrently associated with lower pathological complete response (pCR) rates in CRT-era cohorts. Germline polymorphisms in DNA repair (XRCC1) and folate metabolism (MTHFR) genes showed inconsistent associations with treatment response. In contrast, transcriptomic biomarkers demonstrated greater biological coherence, with proliferative, epithelial-mesenchymal transition, and metabolic signatures frequently associated with treatment resistance, while multi-gene classifiers generally outperformed single-gene markers. Among currently available tumour-intrinsic biomarkers, MMR deficiency was the most consistently reported biomarker associated with reduced response to fluoropyrimidine-based regimens, including TNT, although TNT-specific evidence remains comparatively limited. Dynamic circulating tumour DNA (ctDNA) monitoring, particularly ctDNA clearance during or after therapy, was consistently associated with pathological response and long-term oncologic outcomes across reviewed studies, whereas baseline ctDNA levels showed limited predictive value. Overall, the reviewed literature suggests that biomarker research in rectal cancer has evolved from single-gene analyses towards pathway-level and dynamic biomarkers. The integration of transcriptomic signatures, MMR status, and dynamic ctDNA monitoring may represent a promising strategy for personalising neoadjuvant therapy, improving patient selection for organ-preserving approaches, and enhancing oncologic outcomes in locally advanced rectal cancer. Nevertheless, the evidence base remains heterogeneous, and further prospective validation, assay standardisation, and evaluation within contemporary TNT cohorts are required before these biomarkers can be routinely incorporated into clinical decision-making.

Humans↗

PubMed2005

Analysis of methylation-sensitive transcriptome identifies GADD45a as a frequently methylated gene in breast cancer.

Treatment of the breast cancer cell line, MDAMB468 with the DNA methylation inhibitor, 5-azacytidine (5-AzaC) results in growth arrest, whereas the growth of the normal breast epithelial line DU99 (telomerase immortalized) is relatively unaffected. Comparing gene expression profiles of these two lines after 5-AzaC treatment, we identified 36 genes that had relatively low basal levels in MDAMB468 cells compared to the DU99 line and were induced in the cancer cell line but not in the normal breast epithelial line. Of these genes, 33 have associated CpG islands greater than 300 bp in length but only three have been previously described as targets for aberrant methylation in human cancer. Northern blotting for five of these genes (alpha-Catenin, DTR, FYN, GADD45a, and Zyxin) verified the array results. Further analysis of one of these genes, GADD45a, showed that 5-AzaC induced expression in five additional breast cancer cell lines with little or no induction in three additional lines derived from normal breast epithelial cells. The CpG island associated with GADD45a was analysed by bisulfite sequencing, sampling over 100 CpG dinucleotides. We found that four CpG's, located approximately 700 bp upstream of the transcriptional start site are methylated in the majority of breast cancer cell lines and primary tumors but not in DNA from normal breast epithelia or matched lymphocytes from cancer patients. Therefore, this simple method of dynamic transcriptional profiling yielded a series of novel methylation-sensitive genes in breast cancer including the BRCA1 and p53 responsive gene, GADD45a.

Azacitidine↗

PubMed2006

Tale of two metal reducers: comparative proteome analysis of Geobacter sulferreducens PCA and Shewanella oneidensis MR-1.

Geobacter sulfurreducens and S. oneidensis are the subjects of intense research efforts due to their potential applications to bioremediation. The characterization of their proteomes, being done in parallel with the analysis of their genome sequences, transcriptomes, and metabolomes, is providing valuable insights to both their similarities and their differences. A primary target of interest in the proteomes of both of these metal-reducing microbes is the characterization of their c-type cytochromes. The discovery of their full compliment of c-type cytochromes and the description of what growth conditions trigger their expression is central to harnessing their bioremediation potential. Proteome analyses thus far show that both G. sulfurreducens and S. oneidensis share the common location of a majority of their c-type cytochromes in their outer membranes. The c-type cytochromes of G. sulfurreducens, however, appear to be less soluble and therefore more difficult to isolate from the membranes than those expressed by S. oneidensis. The majority of the G. sulfurreducens c-type cytochromes also differ from those of S. oneidensis in that they have higher isoelectric points, most higher than pH 8.0. These characteristics of solubility and isoelectric point could be related and could indicate an underlying functional difference in the strategy for metal reduction between these two microbes. The global proteome results available for G. sulfurreducens and S. oneidensis at the time of this writing are primarily the result of 2DE analysis coupled to protein identification by LC-MS/MS of tryptic peptides from in-gel digests and represent the most abundant proteins detected by Coomassie blue or silver nitrate staining. Currently, several complimentary efforts utilising the 2D-LC-MS/MS approaches are in progress, promising a more complete protein inventory for these microbes in the near future. As these data are added to those already available, the intricate network of metabolic processes, regulation of protein synthesis and protein function, transport of nutrients, and signal transduction will be elucidated. The existing tools of proteomics will be complimented with newer methods such as protein chips and phage display to further characterize these microbial systems. The end result, in the not too distant future, will be predictive models of G. sulfurreducens and S. oneidensis behavior in their natural habitats under a variety of environmental conditions.

Bacterial Proteins↗

PubMed2005

The expression signature of in vitro senescence resembles mouse but not human aging.

BACKGROUND: The biological mechanisms that underlie aging have not yet been fully identified. Senescence, a phenomenon occurring in vitro, limits the number of cell divisions in mammalian cell cultures and has been suggested to contribute to aging. RESULTS: We investigated whether the changes in gene expression that occur during mammalian aging and induction of cellular senescence are similar. We compared changes of gene expression in seven microarray datasets from aging human, mouse and rat, as well as four microarray datasets from senescent cells of man and mouse. The datasets were publicly available or obtained from other laboratories. Correlation measures were used to establish similarities of the expression profiles and gene ontology analyses to identify functional groups of genes that are co-regulated. Robust similarities were established between aging in different species and tissues, indicating that there is an aging transcriptome. Although some cross-species comparisons displayed high correlation, intra-species similarities were more reliable. Similarly, a senescence transcriptome was demonstrated that is conserved across cell types. A similarity between the expression signatures of cellular senescence and aging could be established in mouse, but not in human. CONCLUSION: Our study is the first to use microarray data from several studies and laboratories for dissection of a complex biological phenotype. We demonstrate the presence of a mammalian aging transcriptome, and discuss why similarity between cellular senescence and aging is apparent in aging mice only.

Aging↗

PubMed2006

Reverse transcription using random pentadecamer primers increases yield and quality of resulting cDNA.

Reverse transcription of RNA is an invaluable method for gene expression analysis by real-time PCR or microarray methods. Random primers of varying lengths were compared with respect to their efficiency of priming reverse transcription reactions. The results showed that 15-nucleotide-long random oligonucleotides (pentadecamers) consistently yielded at least 2-fold as much cDNA as did random hexamers using either poly(A) RNA or an amplified version of messenger RNA (aRNA) as a template. The cDNA generated using pentadecamers did not differ in size distribution or the amount of incorporated label compared with cDNA generated with random hexamers. The increased efficiency of priming using random pentadecamers resulted in reverse transcription of >80% of the template aRNA, while random hexamers induced reverse transcription of only 40% of the template aRNA. This suggests a better coverage of the transcriptome when using random pentadecamers over random hexamers. Using the same amount of aRNA as starting material, random pentadecamer-primed reactions resulted in 11-fold more genes being detected in whole transcriptome DNA microarray experiments than random hexamer-primed reactions. The results indicate that random pentadecamers can replace random hexamers in reverse transcription reactions on both poly(A) RNA and amplified RNA, resulting in higher cDNA yields and quality.

DNA Primers↗

PubMed2026

HIF-1 signaling contributes to lenvatinib resistance in patient-derived HCC organoids.

Resistance to lenvatinib remains an important limitation in hepatocellular carcinoma treatment. Six patient-derived organoid lines were established and classified as sensitive or resistant according to ex vivo drug responses, retaining histological and immunophenotypic features of matched parental tumors. Resistant organoids showed unchanged ATP activity, whereas sensitive ones exhibited pronounced morphological changes and reduced ATP activity at higher concentrations. Transcriptome sequencing identified 408 upregulated and 269 downregulated genes in resistant versus sensitive organoids, with HIF-1 signaling among altered pathways. In resistant organoids, lenvatinib increased HIF-1α, ANGPT2, and HK3 mRNA, whereas comparable changes were not detected in sensitive organoids. KC7F2 reduced these transcripts and further decreased ATP activity when combined with lenvatinib. In organoid-derived xenografts, this combination suppressed tumor growth and HIF-1 target expression more than lenvatinib alone, indicating HIF-1 signaling contributes to the resistant phenotype and its inhibition may enhance response.

Drug resistance↗

PubMed2024

Progressive cardiomyopathy with intercalated disc disorganization in a rat model of Becker dystrophy.

Becker muscular dystrophy (BMD) is an X-linked disorder due to in-frame mutations in the DMD gene, leading to a less abundant and truncated dystrophin. BMD is less common and severe than Duchenne muscular dystrophy (DMD) as well as less investigated. To accelerate the search for innovative treatments, we developed a rat model of BMD by deleting the exons 45-47 of the Dmd gene. Here, we report a functional and histopathological evaluation of these rats during their first year of life, compared to DMD and control littermates. BMD rats exhibit moderate damage to locomotor and diaphragmatic muscles but suffer from a progressive cardiomyopathy. Single nuclei RNA-seq analysis of cardiac samples revealed shared transcriptomic abnormalities in BMD and DMD rats and highlighted an altered end-addressing of TMEM65 and Connexin-43 at the intercalated disc, along with electrocardiographic abnormalities. Our study documents the natural history of a translational preclinical model of BMD and reports a cellular mechanism for the cardiac dysfunction in BMD and DMD offering opportunities to further investigate the organization role of dystrophin in intercellular communication.

Animals↗

PubMed2025

Spatial transcriptomics-aided localization for single-cell transcriptomics with STALocator.

Single-cell RNA-sequencing (scRNA-seq) techniques can measure gene expression at single-cell resolution but lack spatial information. Spatial transcriptomics (ST) techniques simultaneously provide gene expression data and spatial information. However, the data quality of the spatial resolution or gene coverage is still much lower than the quality of the single-cell transcriptomics data. To this end, we develop a ST-Aided Locator for single-cell transcriptomics (STALocator) to localize single cells to corresponding ST data. Applications on simulated data showed that STALocator performed better than other localization methods. When applied to the human brain and squamous cell carcinoma data, STALocator could robustly reconstruct the relative spatial organization of critical cell populations. Moreover, STALocator could enhance gene expression patterns for Slide-seqV2 data and predict genome-wide gene expression data for fluorescence in situ hybridization (FISH) and Xenium data, leading to the identification of more spatially variable genes and more biologically relevant Gene Ontology (GO) terms compared with the raw data. A record of this paper's transparent peer review process is included in the supplemental information.

Single-Cell Analysis↗

PubMed2026

Contrasting regulation of protein-coding genes and lncRNA homeologs in allotetraploid Coffea arabica.

A chromosome-level Bourbon assembly revealed that protein-coding homeologs are predominantly co-regulated between subgenomes. In contrast, intergenic lncRNAs display a modest, but statistically consistent bias toward subgenome E across diverse developmental and stress contexts. Coffea arabica is an allotetraploid species derived from natural hybridization between C. canephora and C. eugenioides, which contributed the C and E subgenomes, respectively. This genomic origin poses major challenges for genome assembly, annotation, and the interpretation of gene regulation. In this study, a high-quality genome assembly of C. arabica was generated and annotated, with particular emphasis on identifying protein-coding genes and intergenic long non-coding RNAs (lincRNAs). Homeologous relationships between genes from the C and E subgenomes were established, providing a robust framework to investigate subgenomic conservation and regulatory divergence. Using an extensive collection of publicly available RNA-seq libraries spanning multiple developmental stages, tissues, and environmental conditions, the relative transcriptional contribution of each subgenome was evaluated. On a global scale, gene expression was largely balanced between subgenomes, with no consistent evidence of subgenome dominance. While protein-coding genes showed comparable regulatory behavior across subgenomes, lincRNAs exhibited a more asymmetric expression pattern, suggesting higher subgenome-specific expression that is interpreted here as a consistent directional tendency rather than as evidence of subgenome dominance. Together, these results provide new insights into the regulatory architecture of the C. arabica genome and establish a foundational genomic and transcriptomic resource for future functional studies and crop improvement efforts.

Coffea↗