Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Transcriptomic profiling”

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.

At least 451 records · Page 25Linked to original sources

Thyroid hormone deprivation creates an immunological signature in the mouse liver, involving Kupffer cell presentation as the mouse ages.

PURPOSE: Aging is associated with an increased prevalence of chronic liver diseases suggesting impaired immune and metabolic function. In addition, thyroid hormone (TH) impacts liver physiology and TH deprivation or excess negatively affect organ maintenance. However, whether age-dependent consequences of TH alterations are reflected in a liver-specific adaptation is unknown so far. The present study aimed to characterize the impact of TH deprivation or excess on the liver transcriptome during aging. METHODS: Five- and 21-month-old male C57BL/6 mice were exposed either to chronic TH deprivation or to chronic TH excess and compared to control treatment by microarray-based liver transcriptome analysis. RESULTS: Significant roles of both TH state and age became obvious: Bioinformatic analysis of the liver transcriptome data revealed an age-dependent immune signature by chronic TH deprivation, an age-dependent immune and metabolic signature independent of exogenous TH modulation, as well as an age-dependent metabolic signature by chronic TH excess. Published data of single cell transcriptomic atlas characterizing aging tissues in the mouse were compared with our data and revealed Kupffer cell presentation in the immunological signature by TH deprivation during aging. Literature data for four prominent differentially expressed genes, namely C1qb, C3ar1, Ctss, and Msr1, revealed that the complement system, extracellular matrix remodelling, as well as the proinflammatory phenotype of Kupffer cells are altered by TH deprivation during aging. CONCLUSION: In conclusion, our study illuminates the interplay between TH deprivation, aging, and liver transcriptome signatures, highlighting potential implications for immune function and tissue maintenance, particularly through the modulation of Kupffer cell presentation.

Animals↗

Unveiling the molecular basis of gonadal development: Multi-omics uncovers sex-related genes and steroid pathways in Sinonovacula constricta.

The razor clam Sinonovacula constricta is an economically important cultured mollusk in China, but the molecular mechanism of its gonadal development and sexual differentiation remains unclear. This study integrated gonadal transcriptomic, proteomic, and metabolomic analysis to identify key sex-related molecules. Transcriptome analysis identified 2795 DELs and 6497 DEGs between sexes, including the sex-related genes Fem-1b, Fem-1c, GUCY1B2 and FAT4, as well as a regulatory network of 39 lncRNA-mRNA pairs involving Tektin-4, Ropporin-1, Histone H1, and FoxN4. Proteomic analysis revealed 3217 DEPs: Tektin family members, Ropporin-1 and Tssk proteins were upregulated in the testis, while histone H1 and FAT4 were upregulated in the ovary. Metabolomic analysis detected 409 DEMs, with uridine identified as a potential sex differential marker (upregulated in the ovary), and 23 gonadal development-related DEMs showed sex-specific upregulation. Integrative transcriptome-proteome analysis identified 1543 co-expressed DEGs/DEPs enriched in nucleosome assembly, oxidative phosphorylation, and carbon metabolism, including key sex-related genes AKAP14, Tektin/Tssk families, Histone H1, and FAT4. Transcriptome-metabolome integration identified 32 shared KEGG pathways (e.g., biosynthesis of unsaturated fatty acids, pyrimidine metabolism), while proteome-metabolome integration revealed 5 (positive ion) and 6 (negative ion) co-enriched pathways, with alanine, aspartate and glutamate metabolism and oxidative phosphorylation being functionally relevant to gonadal development. Collectively, these results reveal the molecular basis of gonadal development, highlight critical sex-related genes and steroid metabolic pathways, and provide valuable resources for future reproduction and breeding in S. constricta.

Animals↗

Comparative transcriptome analysis provides insights into dorso-ventral color pattern formation of Holothuria edulis.

Animal body color patterns are highly diverse and play critical roles in camouflage, intraspecific communication, and environmental adaptation. Holothuria edulis, an important echinoderm inhabiting tropical waters, exhibits a typical dorsoventral dichromatism. This unique body color difference represents a key phenotypic trait for its habitat adaptation; however, the core differential genes regulating this trait remain to be elucidated. In this study, comparative transcriptome sequencing was performed on the dorsal and ventral body wall tissues of H. edulis, leading to the identification of a number of differentially expressed genes (DEGs), followed by GO functional annotation and KEGG pathway enrichment analysis. GO enrichment analysis indicated that the DEGs were significantly enriched in functional categories such as extracellular region, peptidase inhibitor activity, and tetrapyrrole binding. KEGG pathway analysis further revealed significant enrichment of protein digestion and absorption, the TNF signaling pathway, and cholesterol metabolism. Notably, the pigmentation-related gene FMO2 was highly expressed in the dorsal body wall tissue, whereas cyp1a1, ZIC1, Slc7a11, WNT-1, and ADAMTS20 were highly expressed in the ventral body wall tissue. This study identified DEGs and enriched pathways associated with dorsoventral body color differences in H. edulis, providing new insights into the molecular regulatory mechanisms underlying body color pattern formation. From the perspective of aquaculture applications, body color is one of the important traits affecting the quality and market value of sea cucumber products. Elucidating the molecular mechanisms of body color variation can provide a scientific basis for molecular marker-assisted breeding of superior sea cucumber variety.

Animals↗

An integrated single-cell and spatial proteotranscriptomics atlas of fibroblast-driven immunoregulation within the human adult oral cavity.

The immunoregulatory architecture of human oral tissues remains poorly defined. We present an integrated single-cell and spatial proteotranscriptomic atlas profiling >250,000 single-cell transcriptomes and >4 million spatially resolved cells across 13 niches. Using our AI-enabled AstroSuite, we defined neighborhoods and interaction modules, revealing peri-epithelial fibroblast-centered hubs enriched in effector cytokines. We harmonized fibroblast subtypes (universal, immune, peri-epithelial, peri-vascular, peri-neural, antigen-presenting cell [APC]-like, stress responsive, and myofibroblasts) with stress-responsive subtypes partitioning between mucosae and glands (type I and II). Spatial multiomics mapped ligand-receptor programs and identified mucosal stress-responsive fibroblasts as putative immunoregulatory hubs. Niche-aware integration of healthy and diseased datasets revealed fibroblast rewiring into inflammatory and reparative niches. Disease neighborhoods exhibited expansion of major histocompatibility complex (MHC)-I+, MHC-II+, and programmed cell death ligand 1 (PD-L1)+ fibroblasts and predicted spatial engagement with T cells at tertiary lymphoid structures. Together, this atlas identifies fibroblasts as central regulators of structural immunity and provides a scalable framework to target stromal-immune interactions across barrier organs.

Journal Article↗

Integrated proteomic network analysis reveals PTPRC as a central hub protein orchestrating co-expression modules and metabolic dysregulation in renal carcinoma: PTPRC protein molecular action.

The occurrence of renal carcinoma is closely related to a variety of molecular mechanisms and metabolic disorders. PTPRC (protein tyrosine phosphatase receptor C), as an important regulatory protein, was studied to reveal the role of PTPRC in renal carcinoma through comprehensive proteomic network analysis, especially its core position in the coordination of co-expression modules and metabolic disorders. This study was the first to download and process multiple publicly available renal cancer transcriptome data to conduct differential gene expression analysis across datasets. Functional enrichment and disease ontology analysis were performed on the transcriptome of renal cancer, and weighted gene co-expression network (WGCNA) was constructed. The results showed that comprehensive principal component analysis revealed significant differences in the transcriptome of renal cancer, and functional annotation revealed specific pathways associated with renal cancer. WGCNA analysis identified tumor-associated co-expression modules, while multi-omics analysis further identified core regulatory networks including PTPRC. As a central hub protein, PTPRC plays an important coordinating role in the co-expression module and metabolic dysregulation of renal carcinoma. This discovery provides a new perspective for understanding the molecular mechanism of kidney cancer.

Humans↗

From stem cells to somites: Revealing genetic and exogenous factors of human embryogenesis.

Stem-cell-based human embryo models offer an ethically tractable platform for studying early human development. This study employs somitoids, three-dimensional models of human somitogenesis, to investigate how transcriptional programs and culture conditions influence somite formation and segmentation. We show that pre-differentiation culture medium impacts the developmental potential of induced pluripotent stem cells (iPSCs), with StemFit medium and Matrigel embedding outperforming mTeSR Plus medium in generating robust somite-like structures. Strikingly, these differences arise despite only subtle changes in transcriptomic and time-resolved proteomic profiles. P300-based proximity labeling also reveals a largely overlapping set of chromatin-associated regulators across iPSC conditions. In somitoids, enhancer-associated profiling highlights factors linked to somitogenesis, including MESP2 and TBX6. Knockout of three identified regulators, BPTF, RBPJ, and CITED2, demonstrate their essential roles in somite formation. Together, these findings highlight how culture conditions and enhancer-associated networks influence early human development and demonstrate somitoids as a scalable system for functional genomics.

Humans↗

Genetic association between epilepsy and gliomas: Insights from Mendelian randomization and single-cell transcriptomic analyses.

BACKGROUND: Seizures are prevalent in glioma patients, especially in those with low-grade gliomas. The interaction between gliomas and epilepsy involves complex biological mechanisms that are not fully understood. METHODS: We collected Genome-Wide Association Study data for epilepsy and gliomas, performed differential expression analysis, and conducted Gene Ontology (GO) enrichment analysis on the identified genes. Single-cell RNA sequencing data (scRNA-seq) from GSE221534 dataset in Gene Expression Omnibus (GEO) were used to analyze cell-cell interactions within glioma samples from patients with and without epilepsy. RESULTS: Mendelian Randomization (MR) analysis revealed significant associations between genetic variants related to epilepsy and glioma risk, suggesting a potential causal relationship, especially in astrocytomas. Differential expression analysis identified epilepsy-related genes that were significantly upregulated in astrocytoma tissues compared to normal brain tissues. GO enrichment analysis indicated that these genes are involved in critical biological processes such as neurogenesis and cellular signaling. The scRNA-seq analysis showed, compared to non-epileptic samples, glioma stem cells, microglia, and NK cells are increased in the core regions of astrocytomas in epileptic patients. Additionally, intercellular communication between tumor cells and other non-tumor cells is markedly enhanced in astrocytoma samples from epileptic patients. CONCLUSION: This study provides evidence of a genetic association between epilepsy and gliomas and elucidates the biological mechanisms through which epilepsy may influence glioma progression.

Humans↗

Early oligodendrocyte dysfunction signature in Alzheimer's disease: Insights from DNA methylomics and transcriptomics.

Much research into the aetiology of Alzheimer's disease (AD) has focused on neuronal cell types, while studies on the contribution of glial cells, particularly oligodendrocytes (OLGs), are only starting to emerge. Altered brain DNA methylation, an epigenetic modification that provides the interplay between genetics and environmental cues to tightly regulate gene expression, is well documented in AD. Yet, cell-type-specific investigations remain limited. Here, we examine the role of DNA methylation and OLGs in AD, and how such changes may impact gene expression. We performed weighted-gene correlation network analysis (WGCNA) on multiple brain omics AD datasets across species: human DNA methylation data from 4 brain regions, human brain single-nuclei RNA sequencing data and mouse brain RNA sequencing data. We compared AD-associated network modules enriched for OLG genes across AD brain regions, as well as with other neurodegenerative disease DNA methylation datasets. We identified a DNA methylation signature associated with AD, enriched for OLGs, and preserved across brain regions representing early and late AD pathology stages. Genes within this signature showed altered expression in AD OLGs, confirming cell-type specificity and relevance to AD. This OLG signature was also preserved in transgenic mice with early Aβ pathology and in other neurodegenerative diseases without Aβ pathology. We reveal a consistent pattern of OLG dysfunction spanning early to late stages of AD, across DNA methylation and gene expression. Our findings highlight OLG-associated DNA methylation changes as important in AD pathogenesis, and possibly in other neurodegenerative diseases, opening new avenues for therapeutic development.

Alzheimer Disease↗

Full-length single-cell spatial transcriptomics reveals spatial and cell-type-specific transcript isoforms in the primate brain.

The primate brain exhibits complex RNA alternative splicing heterogeneity crucial for functional complexity, yet systematic spatial isoform characterization has been lacking. We developed Fullscope-seq, a full-length single-molecule large field-of-view spatial transcriptomics sequencing method at single-cell resolution, based on programmed concatenation cDNA for multiple long-read sequencing platforms. Applying Fullscope-seq to the macaque brain, we uncovered thousands of genes exhibiting differential transcript usage (DTU) across cortical layers, cell types and brain regions. Fullscope-seq resolved hundreds of major isoform switches across distinct brain regions and identified DTUs between superficial and deep cortical layers. Cortical layer-specific DTUs showed cell-composition dependence, whereas regional DTUs were regulated according to both cellular composition and spatial contexts. These isoform variations showed substantial enrichment for neuropsychiatric disorder-associated genes and were conserved across platforms and species. Our study establishes a scalable framework for spatial isoform analysis and provides a resource for understanding transcriptomic diversity in complex tissues.

Animals↗

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↗

Bacterial stress responses lower mRNA-protein level correlations.

Diverse bacterial pathogens have evolved complex regulatory mechanisms to adapt to various environmental stresses during infection. The uncertainty in mRNA-protein levels in response to environmental stressors complicates our understanding of bacterial physiology and their adaptation to stressful environments. To examine this issue, we have integrated transcriptomics and proteomics data on three human bacterial pathogens: Salmonella enterica Typhimurium, Yersinia pseudotuberculosis, and Staphylococcus aureus under 10 infection-relevant stress conditions. We observed positive correlations between mRNA and protein levels, which were decreased under different stress conditions. Essential genes exhibited higher expression levels with lower variation across the conditions and stronger mRNA-protein correlations compared to nonessential genes, highlighting their critical role in bacterial adaptability and survival. Moreover, we identified a substantial number of genes with stress-induced noncorrelating mRNA-protein levels, particularly under conditions triggering strong stress responses. Particularly this level was dramatically lowered for osmotic stress-specific genes affected by impaired translational activity under osmotic stress. Our findings highlight the prevalence of noncorrelating mRNA-protein levels and the potential role of posttranslational modifications in modulating protein levels in response to environmental stressors during infection. This study provides a comprehensive framework for integrating transcriptomics and proteomics data and identifies potential gene products that might significantly impact the ability of diverse bacterial pathogens to adapt to hostile infection environments.

RNA, Messenger↗

Spatial mutual nearest neighbors for spatial transcriptomics data.

MOTIVATION: Mutual nearest neighbors (MNN) is a widely used computational tool to perform batch correction for single-cell RNA-sequencing data. However, in applications such as spatial transcriptomics, it fails to take into account the 2D spatial information. RESULTS: Here, we present spatialMNN, an algorithm that integrates multiple spatial transcriptomic samples and identifies spatial domains. Our approach begins by building a k-nearest neighbors (kNN) graph based on the spatial coordinates, prunes noisy edges, and identifies niches to act as anchor points for each sample. Next, we construct a MNN graph across the samples to identify similar niches. Finally, the spatialMNN graph can be partitioned using existing algorithms, such as the Louvain algorithm to predict spatial domains across the tissue samples. We demonstrate the performance of spatialMNN using large datasets, including one with N = 31 10x Genomics Visium samples. We also evaluate the computing performance of spatialMNN to other popular spatial clustering methods. AVAILABILITY AND IMPLEMENTATION: Our software package is available on GitHub (https://github.com/Pixel-Dream/spatialMNN). The code is available on Zenodo (https://doi.org/10.5281/zenodo.15073963).

Algorithms↗

Enhancing and accelerating cell type deconvolution of large-scale spatial transcriptomics slices with dual network model.

MOTIVATION: Cell type deconvolution deciphers spatial distribution of mRNA transcripts at single cell level by integrating single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics data to infer mixture of cell types of spots in slices. Current algorithms are criticized for neglecting connection between scRNA-seq and spatial transcriptomics data, as well as time-consuming, hampering their application to large-scale datasets. RESULTS: In this study, we propose a joint learning nonnegative matrix factorization algorithm for fast cell type deconvolution (aka jMF2D), which integrates scRNA-seq and spatial transcriptomics data with network models. To bridge scRNA-seq and spatial transcriptomics data, jMF2D jointly learns cell type similarity network to enhance quality of signatures of cell types, thereby promoting accuracy and efficiency of deconvolution. Experiments demonstrate that jMF2D outperforms state-of-the-art baselines in terms of accuracy by saving about 90% running time on various datasets generated by different platforms. Furthermore, it can also facilitates the identification of spatial domains and bio-marker genes, providing an efficient and effective model for analyzing spatial transcriptomics data. AVAILABILITY AND IMPLEMENTATION: The software is coded using python, and is free available for academic https://github.com/xkmaxidian/jMF2D.

Algorithms↗

SpatialRNA: a Python package for easy application of Graph Neural Network models on single-molecule spatial transcriptomics dataset.

SUMMARY: Image-based spatial transcriptomics (iST) deliver gene expression measurements of RNA transcripts in tissue slices with single-molecule resolution and spatial context preserved. Modern Graph Neural Network (GNN) models are promising methods for capturing the complex molecular and cellular phenotypes in tissues at single-transcript and single-cell levels. A key application of GNNs is the detection of spatial domains or niches, that is, groups of molecules and/or cells that collaboratively work together to produce complex phenotypes. Due to the vast number of detected transcripts in (iST) dataset, applying GNNs on RNA molecule graphs is not trivial. We present a Python package, SpatialRNA, for easy (sub)graph generation from tissue samples and provide comprehensive tutorials for convenient and efficient application of Graph Neural Network models under the PyG framework. This highly scalable tool comprehensively segments tissue into spatial domains, aiding in biological interpretation of iST data and its underlying molecular microenvironments. AVAILABILITY AND IMPLEMENTATION: The SpatialRNA package is freely accessible from online repository https://github.com/ruqianl/spatialrna and can be installed via pip. Comprehensive tutorials, guidance on parameter selection, and complete workflows of case studies are available from the documentation website https://ruqianl.github.io/spatialrna_docs/, and uploaded on Zenodo with a DOI 10.5281/zenodo.17339575.

Neural Networks, Computer↗

Unsupervised multiscale clustering of single-cell transcriptomes to identify hierarchical structures of cell subtypes.

BACKGROUND: Cell clustering is an essential step in uncovering cellular architectures in single-cell RNA sequencing (scRNA-seq) data. However, the existing cell clustering approaches are not well designed to dissect complex structures of cellular landscapes at a finer resolution. RESULTS: Here, we develop a multiscale clustering (MSC) approach to construct a sparse cell-cell correlation network for unsupervised identification of de novo cell types and subtypes across multiple resolutions. Based upon simulated silver- and gold-standard data as well as real scRNA-seq data in diseases, MSC demonstrates significantly improved performance compared to established benchmark methods and reveals a biologically meaningful cell hierarchy to facilitate the discovery of novel disease-associated cell subtypes and mechanisms. CONCLUSIONS: We present MSC as a new single-cell multiscale clustering framework as a powerful tool for advancing discoveries in disease-associated cell populations using single-cell sequencing data.

Single-Cell Analysis↗

Dynamic allelic expression in mouse mammary glands across the adult developmental cycle.

The mammary gland, which primarily develops postnatally, undergoes significant changes during pregnancy and lactation to facilitate milk production. Through the generation and analysis of 480 transcriptomes, we provide the most detailed allelic expression map of the mammary gland, cataloguing cell-type-specific expression from ex-vivo purified cell populations over 10 developmental stages, enabling comparative analysis. The work identifies genes involved in the mammary gland cycle, parental-origin-specific and genetic background-specific expression at cellular and temporal resolution, genes associated with human lactation disorders and breast cancer. Genomic imprinting, a mechanism regulating gene expression based on parental origin, is crucial for controlling gene dosage and stem cell potential throughout development. The analysis identified 25 imprinted genes monoallelically expressed in the mammary gland, with several showing allele-specific expression in distinct cell types. No novel imprinted genes were identified and the absence of biallelically expressed imprinted genes suggests that, unlike in brain, selective absence of imprinting does not regulate gene dosage in the mammary gland. This research highlights transcriptional dynamics within mammary gland cells and identifies novel candidate genes potentially significant in the tissue during pregnancy and lactation. Overall, this comprehensive atlas represents a valuable resource for future studies on expression and transcriptional dynamics in mammary cells.

Animals↗

SAGE transcript profiles of normal primary human hepatocytes expressing oncogenic hepatitis B virus X protein.

Hepatitis B virus (HBV) is a major risk factor for hepatocellular carcinoma (HCC). HBV encodes an oncogenic HBx gene that functions as a transcriptional coactivator of multiple cellular genes. To understand the role(s) of HBx in the early genesis of HCC, we systematically analyzed gene expression profiles by serial analysis of gene expression (SAGE) in freshly isolated human primary hepatocytes infected with a replication-defective adenovirus containing HBx. A total of 19,501 sequence tags (representing 1443 unique transcripts) were analyzed, which provide a distribution of a transcriptome characteristic of normal hepatocytes and a profile associated with HBx expression. Examples of the targeted genes were confirmed by the Megarray analysis with a significant correlation between quantitative SAGE and Megarray (r = 0.8, P < 0.005). In HBx-expressing hepatocytes, a total of 57 transcripts (3.9%) were induced, and 46 transcripts (3.3%) were repressed by more than fivefold. Interestingly, most of the HBx-up-regulated transcripts can be clustered into three major classes, including genes that encode ribosomal proteins, transcription factors with zinc-finger motifs, and proteins associated with protein degradation pathways. These results suggest that HBx may function as a major regulator in common cellular pathways that, in turn, regulate protein synthesis, gene transcription, and protein degradation.

Adenoviridae↗

A Functionally Conserved yet Dynamically Evolving Toolkit Underpinning Molluscan Biomineralization: Insights From Shell and Radula.

The molluscan shell and radula constitute pivotal molluscan innovations, each characterized by distinct functions and diverse forms, regulated by the highly specific biomineralization regulatory networks. Despite their paramount importance, the conserved components and adaptive evolutionary processes governing these regulatory networks remain unresolved. To address this knowledge gap, we advocate for the integration of data from less-explored lineages, such as Scaphopoda, as an essential step. This study presents the inaugural comprehensive transcriptome analysis of Pictodentalium vernedei, a representative species of Scaphopoda distinguished by a unique and evolutionarily conserved shell morphology and radula structure. Furthermore, comparative transcriptome/genome analyses are employed to unravel the conservatism and evolutionary innovation of the involved biomineralization regulatory elements. Our findings underscore the central role of secretomes in governing biomineralization processes, and we identified a fundamental set of 26 domains within molluscan secretomes, forming an essential functional protein domain repertoire necessary for the transformation of inorganic ions into biomineralized structures. This core biomineralization toolkit has undergone independent expansion and lineage-specific recruitment, giving rise to novel, modular domain architectures. This may be essential for the functional specialization and morphological diversification of shell and radula structures. These evolutionary processes are driven by the independent co-option of ancient genes and the emergence of novel de novo genes. This comprehensive investigation not only contributes insights into the evolution of molluscan biomineralization structures but also establishes avenues for further scholarly exploration.

Animals↗