SEARCH · Search PubMed

Results for “gene expression evolution”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

SourceResearch connectionFrom yearThrough year

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

At least 19 recordsLinked to original sources

PubMed2018

Convergent evolution of gene expression in two high-toothed stickleback populations.

Changes in developmental gene regulatory networks enable evolved changes in morphology. These changes can be in cis regulatory elements that act in an allele-specific manner, or changes to the overall trans regulatory environment that interacts with cis regulatory sequences. Here we address several questions about the evolution of gene expression accompanying a convergently evolved constructive morphological trait, increases in tooth number in two independently derived freshwater populations of threespine stickleback fish (Gasterosteus aculeatus). Are convergently evolved cis and/or trans changes in gene expression associated with convergently evolved morphological evolution? Do cis or trans regulatory changes contribute more to gene expression changes accompanying an evolved morphological gain trait? Transcriptome data from dental tissue of ancestral low-toothed and two independently derived high-toothed stickleback populations revealed significantly shared gene expression changes that have convergently evolved in the two high-toothed populations. Comparing cis and trans regulatory changes using phased gene expression data from F1 hybrids, we found that trans regulatory changes were predominant and more likely to be shared among both high-toothed populations. In contrast, while cis regulatory changes have evolved in both high-toothed populations, overall these changes were distinct and not shared among high-toothed populations. Together these data suggest that a convergently evolved trait can occur through genetically distinct regulatory changes that converge on similar trans regulatory environments.

Alleles↗

PubMed2026

Highlights from the 14th International Conference for Plant Mitochondrial Biology: Current Trends and Future Directions.

Plant mitochondrial biology is undergoing a rapid transformation driven by advances in genomics, structural biology, quantitative imaging, and genome engineering. Once focused primarily on respiration and bioenergetics, the field now encompasses diverse areas including genome evolution, gene expression, organelle dynamics, stress signaling, metabolism, and biotechnology. The 14th International Conference for Plant Mitochondrial Biology (ICPMB), held in Kagoshima, Japan, from 18-22 May 2026 (Fig. 1), brought together researchers to discuss recent advances across these rapidly expanding research areas. This meeting report summarizes the major scientific advances presented at ICPMB 2026 and highlights emerging directions that are defining the future of plant mitochondrial biology.

Cytoplasmic male sterility (CMS)↗

PubMed2025

Transposable elements drive evolution and perturb gene expression in Brassica rapa and B. oleracea.

Transposable elements (TEs) significantly influence genomic diversity and gene regulation in plants. Brassica rapa and B. oleracea, with their distinct domestication histories, offer excellent models to explore TE dynamics. Here, we developed a refined TE classification method and systematically analyzed TEs across 12 B. rapa and B. oleracea genomes, identifying 1878 TE families. Approximately half (49.5%) of these TE families were shared between the two species, reflecting a common evolutionary origin, whereas species-specific expansions, particularly among long-terminal repeat (LTR) retrotransposons, underscore their roles in genomic differentiation. We notably characterized a heat-responsive Ty1-copia family (Copia0035) in B. oleracea roots, distinguished by low GC content and the absence of CG and CHG methylation motifs, sharing regulatory similarities with the Arabidopsis heat-induced ONSEN element. Syntenic analyses of gene-TE associations highlighted significant intraspecies TE insertion variability, with more accession-specific insertions in B. rapa and more conserved insertions, often associated with distinct morphotypes in B. oleracea. Gene ontology enrichment indicated TE involvement in developmental, reproductive, and stress response pathways. Transcriptome analysis across diverse accessions revealed that genes proximal to TEs, particularly those regulating floral development and flowering time, exhibit increased expression variability. These findings advance our understanding of TE-mediated genome evolution in Brassica species and underscore their potential utility in breeding and genome engineering strategies for crop improvement.

DNA Transposable Elements↗

PubMed2026

Evolution and Expression Divergence of Legume PAL Genes Suggest Associations with Drought Response and Root Nodule Development.

Comparative genomic analyses provide insight into the mechanisms underlying gene-family evolution and crop adaptation. Here, we used the legume phenylalanine ammonia-lyase (PAL) gene family as a model and integrated pan-genomic, phylogenetic, molecular evolutionary, duplication-mode, and transcriptomic analyses, while developing GFtool for gene family identification. Across 45 genomes, we identified 302 PAL genes and classified them into five Groups. Groups 1-3 represented ancient lineages shared with outgroups, whereas Groups 4 and 5 were legume-specific. Molecular-clock analyses placed the divergence of Group 2 near the Paleocene-Eocene transition, while Groups 4 and 5 diversified from the middle Eocene to the early Oligocene. WGD/segmental duplication broadly contributed to PAL copy-number expansion, whereas tandem duplication was enriched in Group 5 of Papilionoideae. Group 2 genes showed drought-induced expression, whereas Group 5 genes were associated with early root nodule development. GFtool provides a scalable framework for gene-family studies.

Fabaceae↗

PubMed2024

Stochastic modeling of single-cell gene expression adaptation reveals non-genomic contribution to evolution of tumor subclones.

Cancer progression is an evolutionary process driven by the selection of cells adapted to gain growth advantage. We present a formal study on the adaptation of gene expression in subclonal evolution. We model evolutionary changes in gene expression as stochastic Ornstein-Uhlenbeck processes, jointly leveraging the evolutionary history of subclones and single-cell expression data. Applying our model to sublines derived from single cells of a mouse melanoma revealed that sublines with distinct phenotypes are underlined by different patterns of gene expression adaptation, indicating non-genetic mechanisms of cancer evolution. Sublines previously observed to be resistant to anti-CTLA4 treatment showed adaptive expression of genes related to invasion and non-canonical Wnt signaling, whereas sublines that responded to treatment showed adaptive expression of genes related to proliferation and canonical Wnt signaling. Our results suggest that clonal phenotypes emerge as the result of specific adaptivity patterns of gene expression. A record of this paper's transparent peer review process is included in the supplemental information.

Animals↗

PubMed2025

Genome-Wide Mining of lncRNAs Reveals Their Potential Regulatory Role in the Evolution of Viviparity.

Reproduction in vertebrates usually involves egg-laying (oviparity) or live-bearing (viviparity). Oviparity is the ancestral trait from which viviparity has independently evolved more than 100 times in squamate reptiles. This transition involves a series of physiological and structural changes, including the degeneration of eggshell and the evolution of a placenta and differences in the temporal and spatial expression patterns of some functional genes that drive the structural transformation. Long non-coding RNAs (lncRNAs) play important roles in the regulation of gene expression, yet it remains unclear whether they participate in gene expression shifts during the transition from oviparity to viviparity, and if so how. Therefore, we employ deep mining to identify novel lncRNAs of a closely related oviparous-viviparous pair of lizards (Phrynocephalus przewalskii and P. vlangalii). We construct cis- and trans-regulatory networks between lncRNAs and target genes using the transcriptomic data of oviduct or uteri tissues across reproductive periods. Results show that lncRNAs that regulate eggshell gland developmental genes in the oviparous lizard are lost or less expressed in the viviparous lizard. A number of lncRNAs involved in the regulation of placental development and embryo attachment in viviparous species have no orthologs in oviparous species, and others show little or no expression. Accordingly, lncRNAs may play important regulatory roles in the physiological and structural changes in the transition from oviparity to viviparity. These results open doors to the further elucidation of genetic regulatory networks.

Animals↗

PubMed2025

The evolutionary trajectories and gene regulatory roles of nuclear-integrated plastid DNA: clues for enhancing environmental adaptation in Caryophyllales.

Environmental stimuli can induce the transfer of chloroplast DNA to the nuclear genome, resulting in nuclear-integrated plastid DNAs (NUPTs). However, their role in plant adaptability remains unclear. Species within the Caryophyllales order, known for their adaptation to extreme environments, provide an ideal model for studying the evolutionary dynamics and functions of NUPTs. In this study, we analyzed NUPTs in 24 Caryophyllales species to investigate their evolution and regulatory roles in gene expression, particularly in response to environmental stimuli. We found significant interspecies variation in NUPT abundance, ranging from 566 insertions in Amaranthus cruentus to 3585 in Beta vulgaris, with sizes spanning from 100 bp to over 100 kb. Approximately 62% of NUPTs were inserted within the last 20 million years, while some species exhibit insertion peaks dating back 49 million years. NUPT presence/absence polymorphisms in six related species suggest that NUPT insertions and deletions are dynamic processes influenced by phylogeny. NUPTs predominantly integrate into intergenic regions but also insert into genes and promoters, with certain regions acting as hotspots. Notably, NUPTs introduce numerous environmental-responsive cis-acting elements in promoter regions. Genes with NUPT insertions in their promoters are significantly enriched for functions related to environmental response. Further luciferase assays in Spinacia oleracea demonstrated that NUPT insertions can regulate the expression of genes related to environmental responses, indicating their potential role in adaptive evolution. Overall, our study provides insights into NUPT evolution and their influence on gene function and plant adaptability to environmental stimuli.

Plastids↗

PubMed2025

The release of sexual conflict after sex loss is associated with evolutionary changes in gene expression.

Sexual conflict can arise because males and females, while sharing most of their genome, can have different phenotypic optima. Sexually dimorphic gene expression may help reduce conflict, but the expression of many genes may remain sub-optimal owing to unresolved tensions between the sexes. Asexual lineages lack such conflict, making them relevant models for understanding the extent to which sexual conflict influences gene expression. We investigate the evolution of sexual conflict subsequent to sex loss by contrasting the gene expression patterns of sexual and asexual lineages in the pea aphid Acyrthosiphon pisum. Although asexual lineages of this aphid produce a small number of males in autumn, their mating opportunities are limited because of geographic isolation between sexual and asexual lineages. Therefore, gene expression in parthenogenetic females of asexual lineages is no longer constrained by that of other morphs. We found that the expression of genes in males from asexual lineages tended towards the parthenogenetic female optimum, in agreement with theoretical predictions. Surprisingly, males and parthenogenetic females of asexual lineages overexpressed genes normally found in the ovaries and testes of sexual morphs. These changes in gene expression in asexual lineages may arise from the relaxation of selection or the dysregulation of gene networks otherwise used in sexual lineages.

Animals↗

PubMed2026

Expression patterns of risk genes associated with three evolutionarily relevant syndromes in rhesus macaque and human brains.

Depressive disorder (DD), Alzheimer's disease (AD), and schizophrenia (SZ) are evolutionarily relevant traits that disrupt neural networks supporting affect and cognition. While genome-wide association studies have identified risk-related genes for these diseases, how the expression of these genes compare across species remains unclear. In this study, we examined the spatial and temporal expression of ~2000 disease-associated genes in human and rhesus macaque brains. Distinct cross-species signatures emerged. DD-linked genes showed broad cortical-subcortical expression in humans but were confined to subcortical regions in macaques. The divergent subset was enriched for neuron differentiation, migration, synaptic signalling, and cognition. SZ-linked genes were expressed across cortical-subcortical-cerebellar structures in humans. AD-linked genes showed postnatal cortical-hippocampal macaque expression, and broader cortical-subcortical human expression. Cross-species spatial comparisons revealed a significant negative correlation for DD genes, suggesting a broader spatial distribution of DD-related gene expression in humans, extending to distributed emotion-cognition networks, compared to affective hubs in macaques. SZ genes exhibited a similar, though non-significant, negative trend, while AD genes showed a weak, non-significant correlation, indicating an absence of systematic expression shifts. Together, evolutionary shifts in gene expression may have shaped emotional and cognitive functions in humans, and susceptibility to psychiatric and neurodegenerative disorders.

Animals↗

PubMed2024

A mammalian tripartite enhancer cluster controls hypothalamic Pomc expression, food intake, and body weight.

Food intake and energy balance are tightly regulated by a group of hypothalamic arcuate neurons expressing the proopiomelanocortin (POMC) gene. In mammals, arcuate-specific POMC expression is driven by two cis-acting transcriptional enhancers known as nPE1 and nPE2. Because mutant mice lacking these two enhancers still showed hypothalamic Pomc mRNA, we searched for additional elements contributing to arcuate Pomc expression. By combining molecular evolution with reporter gene expression in transgenic zebrafish and mice, here, we identified a mammalian arcuate-specific Pomc enhancer that we named nPE3, carrying several binding sites also present in nPE1 and nPE2 for transcription factors known to activate neuronal Pomc expression, such as ISL1, NKX2.1, and ERα. We found that nPE3 originated in the lineage leading to placental mammals and remained under purifying selection in all mammalian orders, although it was lost in Simiiformes (monkeys, apes, and humans) following a unique segmental deletion event. Interestingly, ablation of nPE3 from the mouse genome led to a drastic reduction (>70%) in hypothalamic Pomc mRNA during development and only moderate (<33%) in adult mice. Comparison between double (nPE1 and nPE2) and triple (nPE1, nPE2, and nPE3) enhancer mutants revealed the relative contribution of nPE3 to hypothalamic Pomc expression and its importance in the control of food intake and adiposity in male and female mice. Altogether, these results demonstrate that nPE3 integrates a tripartite cluster of partially redundant enhancers that originated upon a triple convergent evolutionary process in mammals and that is critical for hypothalamic Pomc expression and body weight homeostasis.

Animals↗

PubMed2026

A high-quality draft genome assembly of Johnsongrass illuminates relationships between polyploidization, crop-wild hybridization, and reproductive biology.

Johnsongrass [Sorghum halepense (L.) Pers.] is an allopolyploid, rhizomatous, perennial grass species and one of the most troublesome weeds in global agriculture. We assembled the first Johnsongrass genome to clarify poorly understood genetic factors influencing variable rates of crop-wild hybridization with cultivated sorghum [S. bicolor (L.) Moench]. The draft genome assembly has a total size of 3.26 Gb and BUSCO completeness of 95.3%. We also report the first evolutionary analysis of INHIBITION OF ALIEN POLLEN (IAP), the only known cross-(in)compatibility locus in the genus. Our results reveal an evolutionary history of genome instability, including the loss of distinct parental subgenomes, and suggest that Nebraska accession 'J-37,' the genome donor, is a segmental allotetraploid that may function as a diploid or aneuploid during meiosis. Genome instability could explain observations of variable ploidies in Johnsongrass and facilitate ongoing hybridization with sorghum where gamete ploidies and IAP alleles match. Given this information, we provide a suggested research framework for studying evolution and gene expression in the Sorghum genus where crop-wild hybridization occurs and for predicting the potential for hybridization between specific crossing partners. Collectively, this work will bolster efforts to study and manage reproductive biology in other crop-wild polyploid complexes.

Sorghum↗

PubMed2026

Heterokairic Genes and the Eco-Evo-Devo of Timing.

Concepts of developmental timing have traditionally been framed under heterochrony as evolved (genetically based) differences in timing, while environmentally induced shifts in timing within genotypes have been treated more loosely. In this article, heterokairy is presented as plasticity in the timing of developmental events, and the term "heterokairic genes" is proposed for environmentally modulated heterochronic genes that underlie this plasticity. Evidence from nematodes, insects, plants, and vertebrates is assembled, with emphasis placed on systems where environmental cues are relayed through endocrine or metabolic pathways to known timing modules/genes. On this basis, a distinction is drawn between validated heterokairic genes, supported by direct mechanistic data, and a broader set of candidates inferred from gene-environment interactions in developmental timing. The eco-evolutionary consequences of such genes are considered, and experimental and genomic strategies for their identification are outlined. It is argued that heterokairic genes provide a useful bridge between environmental variation, developmental mechanisms, and evolutionary change in timing.

Animals↗

PubMed2004

Genomic regionality in rates of evolution is not explained by clustering of genes of comparable expression profile.

In mammalian genomes, linked genes show similar rates of evolution, both at fourfold degenerate synonymous sites (K4) and at nonsynonymous sites (KA). Although it has been suggested that the local similarity in the synonymous substitution rate is an artifact caused by the inclusion of disparately evolving gene pairs, we demonstrate here that this is not the case: after removal of disparately evolving genes, both (1) linked genes and (2) introns from the same gene have more similar silent substitution rates than expected by chance. What causes the local similarity in both synonymous and nonsynonymous substitution rates? One class of hypotheses argues that both may be related to the observed clustering of genes of comparable expression profile. We investigate these hypotheses using substitution rates from both human-mouse and mouse-rat comparisons, and employing three different methods to assay expression parameters. Although we confirm a negative correlation of expression breadth with both K4 and KA, we find no evidence that clustering of similarly expressed genes explains the clustering of genes of comparable substitution rates. If gene expression is not responsible, what about other causes? At least in the human-mouse comparison, the local similarity in KA can be explained by the covariation of KA and K4. As regards K4, our results appear consistent with the notion that local similarity is due to processes associated with meiotic recombination.

Animals↗

PubMed2026

Genome-wide characterization of the bZIP gene family in Rattus norvegicus and expression profiling analysis during brain development.

BACKGROUND: The brown rat (Rattus norvegicus) serves as a cornerstone model organism in biomedical research, particularly for understanding physiological homeostasis and stress responses. The basic leucine zipper (bZIP) transcription factor family is a pivotal regulatory network involved in growth, organogenesis, and neurodevelopment. Despite its importance, a systematic characterization of the bZIP gene family in rats has remained elusive. RESULTS: In this study, we performed a genome-wide identification of 61 RnbZIP genes, which were categorized into 10 distinct subfamilies based on phylogenetic relationships and chromosomal localization. Structural analysis revealed conserved motif arrangements within subfamilies, while collinearity analysis identified significant gene duplication events-predominantly tandem and segmental duplications-that have driven the evolutionary expansion of the RnbZIP family. Quantitative analysis showed that members within the same subfamily shared 45%-92% sequence similarity (calculated using the BLOSUM62 scoring matrix), and all duplicated gene pairs underwent strong purifying selection (Ka/Ks < 1). Comparative genomics across seven rodent species further underscored the evolutionary conservation and divergence of these factors. Expression profiling across diverse organs and brain developmental stages indicated that RnbZIP genes exhibit high tissue specificity. Notably, 10 candidate genes, including RnbZIP01, RnbZIP02, and RnbZIP08, demonstrated dynamic expression patterns during brain maturation, suggesting their essential roles in neurodevelopmental processes. CONCLUSIONS: Our findings provide a comprehensive structural and evolutionary framework for the RnbZIP gene family, highlighting their potential regulatory functions in rat organogenesis and brain development. This study establishes a valuable resource for further functional characterization of specific bZIP members in mammalian neurological systems.

Animals↗

PubMed2026

Comprehensive identification and evolutionary analysis of the Wnt gene family in bivalves: Insights into the larval development of the noble scallop Chlamys nobilis.

The Wnt gene family regulates fundamental developmental processes in metazoans, but its evolutionary composition and developmental deployment in bivalves remain largely unresolved. Here, we performed a comparative genomic analysis of Wnt genes in 19 bivalve species and examined developmental expression profiles in the noble scallop Chlamys nobilis, with Crassostrea gigas and Chlamys farreri used for cross-species comparison. A total of 235 Wnt genes were identified and assigned to 12 subfamilies. No reliable Wnt3 ortholog was detected in any analyzed bivalve, supporting the view that Wnt3 loss occurred early during lophotrochozoan evolution rather than representing a lineage-specific absence. Most Wnt proteins retained the conserved WNT domain, indicating strong structural conservation, whereas lineage-specific copy-number variation and gene loss were observed among species. C. farreri and C. gigas each retained 12 Wnt genes and lacked Wnt3, whereas C. nobilis lacked Wnt3, Wnt7, and Wnt16. Developmental transcriptome analysis and RT-qPCR revealed clear stage-specific expression patterns. In C. gigas, Wnt2/10/A were highly expressed during earlydevelopment and peaked around the D-shaped larval stage, while Wnt8 and Wnt11 showed distinct stage-specific peaks. By contrast, Wnt1/5/6/9 were more active during later larval development or juvenile formation. These results provide a comparative framework for bivalve Wnt evolution and identify candidate Wnt genes potentially involved in larval development and aquaculture-relevant developmental transitions.

Animals↗

PubMed2026

Molecular Evolution and Expression Analysis of the ADH Gene Family in Apple Bud Mutants.

Alcohol dehydrogenase (ADH) catalyzes the reduction of aldehydes to alcohols, key precursor substrates for volatile ester biosynthesis, which determines the characteristic aroma of apple fruit. However, a comprehensive genome-wide investigation of the ADH gene family in apple has been lacking. In this study, we systematically identified ADH genes in the apple genome using integrated bioinformatics approaches, including phylogenetic analysis, synteny evaluation, promoter cis-element prediction, codon usage bias assessment, and protein interaction network modeling. Expression patterns were examined through transcriptomic data and validated by RT-qPCR analysis across different organs and among 'Red Delicious' and its four bud mutant lines. We identified 44 ADH genes, with 12 forming a prominent cluster on chromosome 1. RT-qPCR analysis revealed that MdADH20 was dramatically upregulated in the 'Red Chief' mutant (relative expression of 59.38), suggesting its pivotal role. Phylogenetic analysis revealed a close evolutionary relationship with wild strawberry. The encoded proteins were generally stable and predominantly localized to the cytoplasm. Promoter analysis showed enrichment of growth/development-related and ARE elements, while codon usage analysis identified AGA, GCU, GUU, and CUU as preferred codons. Protein interaction prediction suggested MdADH19 and MdADH20 as hub proteins. Expression profiling and RT-qPCR further identified MdADH20 as a core candidate gene, characterized by its stable and high expression, particularly in the 'Red Delicious' mutant. Its central position in the predicted protein-protein interaction network suggests a potential regulatory role in the aroma biosynthesis pathway of apple fruit. This study provides the first systematic genome-wide characterization of the apple ADH gene family, establishing a theoretical groundwork for deciphering aroma biosynthesis mechanisms and offering potential target genes for flavor improvement through bud mutation breeding strategies.

ADH gene family↗

PubMed2011

Identification of transcriptome SNPs between Xiphophorus lines and species for assessing allele specific gene expression within F₁ interspecies hybrids.

Variations in gene expression are essential for the evolution of novel phenotypes and for speciation. Studying allelic specific gene expression (ASGE) within interspecies hybrids provides a unique opportunity to reveal underlying mechanisms of genetic variation. Using Xiphophorus interspecies hybrid fishes and high-throughput next generation sequencing technology, we were able to assess variations between two closely related vertebrate species, Xiphophorus maculatus and Xiphophorus couchianus, and their F(1) interspecies hybrids. We constructed transcriptome-wide SNP polymorphism sets between two highly inbred X. maculatus lines (JP 163 A and B), and between X. maculatus and a second species, X. couchianus. The X. maculatus JP 163 A and B parental lines have been separated in the laboratory for ≈70 years and we were able to identify SNPs at a resolution of 1 SNP per 49 kb of transcriptome. In contrast, SNP polymorphisms between X. couchianus and X. maculatus species, which diverged ≈5-10 million years ago, were identified about every 700 bp. Using 6524 transcripts with identified SNPs between the two parental species (X. maculatus and X. couchianus), we mapped RNA-seq reads to determine ASGE within F(1) interspecies hybrids. We developed an in silico X. couchianus transcriptome by replacing 90,788 SNP bases for X. maculatus transcriptome with the consensus X. couchianus SNP bases and provide evidence that this procedure overcomes read mapping biases. Employment of the in silico reference transcriptome and tolerating 5 mismatches during read mapping allow direct assessment of ASGE in the F(1) interspecies hybrids. Overall, these results show that Xiphophorus is a tractable vertebrate experimental model to investigate how genetic variations that occur during speciation may affect gene interactions and the regulation of gene expression.

Alleles↗

PubMed2025

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↗