Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “gene family expansion”

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 19 recordsLinked to original sources
PubMed2025

Genome-wide identification and expression profiling of the MADS-box gene family in Lavandula angustifolia.

BACKGROUND: MADS-box genes encode transcription factors critical for plant development, particularly floral organogenesis, flowering time regulation, and adaptation to environmental stresses. Among these, the MIKCC-type genes are pivotal regulators in floral developmental processes. Although the evolutionary diversification and functional dynamics of MADS-box genes have been extensively characterized in model plants such as Arabidopsis thaliana and Oryza sativa, their evolutionary relationships and functional profiles in Lavandula angustifolia, an economically significant aromatic plant, remain poorly understood. RESULTS: Genome-wide analysis identified 173 MADS-box genes in L. angustifolia, categorized into type I (Mα: 26; Mβ: 0; Mγ: 10) and type II (MIKCC: 125; MIKC*: 12) based on phylogenetic comparisons with A. thaliana. The MIKCC subgroup was further subdivided into 12 subclasses, including genes central to the ABCDE model of floral organ specification. Structural analyses revealed distinct conserved motifs and exon-intron configurations specific to each subgroup, indicative of functional divergence. Synteny analysis demonstrated Whole Genome Duplication (WGD) and segmental duplications as major contributors to MIKCC gene family expansion, notably among genes linked to floral organ development. Expression profiling via RNA-seq and quantitative real-time PCR (qPCR) showed type II MADS-box genes exhibited higher expression levels with pronounced tissue-specific and developmental stage-specific expression patterns compared to type I genes. Many type II genes displayed significant associations with floral organogenesis, floral transition, and abiotic stress responses, underscoring their essential roles in reproductive development and environmental adaptability in L. angustifolia. CONCLUSIONS: The identification and comprehensive characterization of 173 MADS-box genes in L. angustifolia highlight the significant expansion of the MIKCC subgroup driven primarily by WGD and segmental duplications. The distinct structural features and specific expression patterns observed provide insights into the functional divergence and complexity of these genes, particularly regarding floral organogenesis and adaptation to environmental stress. This study establishes a robust molecular basis for further functional analysis and genetic improvement of aromatic plants.

MADS Domain Proteins↗
PubMed2025

Comparative genomic analysis of Acer tsinglingense and A. davidii provides insights into nervonic acid biosynthesis, population evolution and genome vulnerability of endangered A. tsinglingense.

Global biodiversity is facing threats from climate change, habitat fragmentation, and anthropogenic activities-pressures that particularly endanger endemic and narrowly distributed species. In this study, the high-quality chromosome-level genomes of two ecologically divergent maples were assembled: the endangered and range-restricted Acer tsinglingense (791.40 Mb) and its widespread congener Acer davidii (1291.99 Mb). Phylogenomic analysis indicates that the two species diverged ~16.3 million years ago, with A. tsinglingense showing notable gene family expansions in secondary metabolite pathways. Notably, the 3-ketoacyl-CoA synthase gene family, which is involved in nervonic acid biosynthesis, underwent significant expansion and tandem duplication in A. tsinglingense, exhibiting high expression in buds. Population genomic analysis revealed that, compared with the widely distributed A. davidii, A. tsinglingense possesses lower genetic diversity, higher harmful mutation load, and signatures of a severe population bottleneck during the Late Pleistocene. Genome-environment association analysis further identified climate-adaptive genomic variations linked to five key environmental factors and projected potential genomic offsets under future climate scenarios. The southern lineage of A. tsinglingense exhibited greater climate sensitivity and genomic vulnerability under strong selective pressures, underscoring its importance as a conservation priority. Our research reveals that metabolic specializations in A. tsinglingense (such as the synthesis of nervonic acid) may confer competitive advantages in specific habitats. However, factors including its restricted distribution, historical population bottlenecks, and accumulated genetic load severely constrain its evolutionary potential to cope with rapid climate change. These findings emphasize the importance of elucidating the genomic basis and mechanisms of endangerment in metabolically specialized and threatened plant species to inform effective conservation strategies.

Genome, Plant↗
PubMed2026

The Key Trichoderma-Induced Gene Encoding a DUF568 Domain-Containing Protein Mediates Defense Responses in Wheat.

Genes encoding DUF568 domain-containing proteins participate in plant stress adaptation. To elucidate the functional role of DUF568 domain-containing genes in Trichoderma-induced wheat defense responses against wheat Fusarium crown rot, we performed a genome-wide identification and characterization of the TaDUF568 gene family in hexaploid wheat (Triticum aestivum L.). In this study, a total of 33 TaDUF568 family genes were systematically identified and characterized at the genome-wide level, exhibiting uneven chromosomal distribution and diverse physicochemical properties. Phylogenetic, structural, and collinearity analyses revealed conserved family characteristics among monocot species. Segmental duplication was verified as the primary driver of gene family expansion. Expression profiling revealed divergent tissue-specific expression patterns among TaDUF568 family members, among which TaDUF568.18 was strongly induced by Trichoderma M2. Subcellular localization assays confirmed that TaDUF568.18 is a plasma membrane-localized protein. Functional validation via stable transgenes demonstrated that overexpression of TaDUF568.18 restricted lesion expansion, improved agronomic traits, and enhanced disease resistance. This study is the first to characterize the wheat DUF568 family and confirm that TaDUF568.18 (annotated as TaAIR12) acts as a positive regulator of Trichoderma-mediated wheat defense, providing a valuable gene resource for wheat disease-resistance breeding.

DUF568↗
PubMed2026

Biocontrol potential and molecular basis of predation in a marine raptorial ciliate.

Predator-prey interactions are widespread across organisms and are key drivers of morphological and behavioral evolution. Despite this, predation remains poorly understood among microbial eukaryotes, mostly due to the absence of a tractable experimental system that allows quantitative, reproducible investigation. This study establishes the marine raptorial ciliate Chaenea vorax as a highly efficient predator, with Rosenzweig-MacArthur model simulations based on predation data showing that only a few dozen individuals can eliminate the vast majority of the facultatively pathogenic ciliate Uronema marinum within 1-2 days, providing a quantitative basis for developing predator-based biocontrol strategies in aquaculture. Genomic analysis shows that C. vorax possesses a highly fragmented macronuclear genome enriched with predation-related pathways, including calcium-mediated contractility, cellular proteolysis, toxin expulsion systems, among others. Transcriptomic profiling during predation events further demonstrates significant upregulation of genes involved in cytoskeletal remodeling, proteolytic activity, and cellular detoxification. Evolutionary analyses suggest that C. vorax has an extremely long evolutionary history, exceptionally high nucleotide diversity even among ciliates, and gene family expansions linked to predatory adaptation. Although the prey possesses certain defensive mechanisms (e.g. trichocysts), these are largely ineffective against short-term predation in closed aquatic environments. These findings provide fundamental insights into the molecular basis of predation in ciliates and suggest the potential utility of C. vorax in biocontrol applications targeting pathogenic ciliates.

Ciliophora↗
PubMed2026

The genome of Thesium ramosoides (Santalales) reveals evolutionary dynamics associated with parasitism and alpine adaptation.

Plant species adapting to complex environments experience contrasting selection pressures that drive the expansion and contraction of different gene families. However, few studies have investigated simultaneous genomic responses to such diverse selective forces. Here, we generate a high-quality genome assembly for the hemiparasitic plant Thesium ramosoides, the first for the largest genus in the Santalales, and explore the genomic basis underlying the evolution of parasitism and alpine adaptation. Unlike many other parasitic plants, the Thesium genome has not undergone additional rounds of whole-genome duplication, making it particularly tractable for studying gene family evolution. Our analyses reveal substantial loss of photosynthesis-related genes and contraction of biotic defense gene families, likely reflecting adaptation to a hemiparasitic lifestyle and reduced pathogen pressure at high altitudes. The absence of key root hair development genes correlates with the degenerate root hair phenotype observed in this species. Furthermore, hallmarks of high-altitude adaptation include the expansion of gene families involved in responses to hypoxia. Notably, expansion of gene families associated with meristem development is consistent with the presence of below-ground crown buds that enable rapid regeneration after mountain fires. Unexpectedly, we detected tandem duplication and diversification of the strigolactone receptor gene D14, which regulates secondary shoot formation, but not of its ancestral paralog KAI2, which mediates seed germination in response to the smoke-derived compound karrikin. This finding suggests divergent signaling mechanisms underlying fire adaptation across different parasitic plant lineages. By integrating time-series transcriptomic data, we propose a post-fire "defense first, repair later, recovery last" model, in which resources are reallocated from immediate defense to rapid repair and ultimately to long-term recovery, to explain the adaptation of T. ramosoides to fire-prone habitats. Our study provides critical insights into the complex and contrasting genomic dynamics that drive adaptation to multiple co-occurring selection pressures.

Genome, Plant↗
PubMed2026

De novo genome assemblies of threatened Asian hornbills (Bucerotidae) reveal declining population trajectories during the late Pleistocene.

BACKGROUND: Asian hornbills are flagship species of the wet tropics that face significant threats from hunting, habitat loss, and fragmentation. Despite being conservation flagships, whole genome information is available for only two of the 32 Asian hornbill species. In this study, we provide the first de novo genome assemblies for four hornbill species (Bucerotidae) in Asia. METHODS: We used a combination of long-read and short-read sequencing data to assemble and annotate de novo hybrid genomes of four species of hornbills. We also assembled and compared mitochondrial genomes of these species. Using a comparative genomics approach, we performed orthology assignment and gene evolution analyses to identify unique gene families in Asian hornbills, gene families that showed significant expansion, their functions and structural variation. Furthermore, using the Pairwise Sequentially Markov Coalescent (PSMC) method, we reconstructed demographic histories of hornbill species to examine changes in their population trajectories in the past. RESULTS: We present hybrid genome assemblies for Great Hornbill (B. bicornis - GH), Rufous-necked Hornbill (A. nipalensis- RNH), Malabar Pied Hornbill (A. coronatus- MPH) and Wreathed Hornbill (R. undulatus- WH). The genome sizes of these hornbills range from 1.1 Gb to 1.3 Gb, with over 95.9% completeness and gene prediction BUSCO. We reported 10,525 orthogroups shared among four Asian hornbill species and identified significant expansion in gene families associated with structural keratin development in Asian hornbills compared to their ancestors. We also provide annotated mitogenomes for each of these species. Furthermore, we found that the WH, a more abundant, widely distributed, and migratory species, showed a higher Ne than the other three hornbill species. However, an overall decline in Ne for all species was recorded during the Pleistocene climatic fluctuations. CONCLUSIONS: We present the first-ever, high-quality reference genomes for the threatened hornbill species from Asia. Hornbills have shown significant expansion in genes involved in structural keratin development. Our results indicate that Pleistocene climatic fluctuations have led to dramatic population declines in all four species. We believe that this study provides robust genomic resources to support future comparative and conservation genomics efforts for hornbills.

Animals↗
PubMed2025

First clinical diagnosis of FAME3 via commercial Long-Read sequencing reveals mosaic repeat expansion in MARCHF6 gene.

Familial Adult Myoclonic Epilepsy type 3 (FAME3) is a rare autosomal dominant disorder characterized by cortical tremor and epilepsy, caused by a noncoding pentanucleotide repeat expansion (TTTTA/TTTCA)n in the MARCHF6 gene. Conventional genetic testing often fails to detect this expansion due to its repetitive structure and intronic location. We evaluated a 61-year-old woman with refractory myoclonic and generalized tonic-clonic seizures, whose prior genetic testing-including exome and genome sequencing-was non-diagnostic. Using PacBio HiFi long-read whole-genome sequencing and the tandem repeat genotyping tool TRGT, we identified a pathogenic MARCHF6 intronic expansion. The proband harbored one allele with 15 TTTTA repeats and a second allele with a compound expansion of 661 TTTTA and 12 TTTCA repeats. Three affected relatives shared similarly expanded alleles, but with increasing repeat size in the latter generations. Importantly, analysis using TRGT-instability revealed repeat mosaicism in all affected individuals, reflected by variability in motif counts across individual sequencing reads. This somatic heterogeneity may contribute to the phenotypic penetrance, variable expressivity and pleiotropism seen in FAME3 disease expression. To our knowledge, this is the first clinical diagnosis of FAME3 using a commercially available long-read sequencing platform, underscoring its diagnostic utility in resolving complex repeat expansion disorders and uncovering biologically relevant mosaicism.

Humans↗
PubMed2026

Evolutionary patterns and repeated adaptive strategies of deep-sea anemones.

Sea anemones occupy the full depth range of the oceans, yet their evolutionary patterns and adaptive strategies to the enigmatic deep sea have remained contentious and poorly resolved. Here, we assemble genomes (n = 13) and transcriptomes for 15 species collected between 432 and 6,000 m and integrate them with publicly available actiniarian data. We find support for a shallow-water origin of Actiniaria through a framework that emphasizes genome-scale changes associated with habitat transitions. Most strikingly, these changes include repeated dismantling of the circadian toolkit across deep-sea lineages. In addition to convergent gene losses in photo- and temperature-regulatory genes, we find that some deep-sea lineages have experienced recurrent loss or pseudogenization of key meiotic genes (e.g., Meiosin, Ythdc2, Spo11, and Mlh3), suggesting reduced meiotic capacity in some lineages. Despite this extensive genomic erosion, deep-sea anemones exhibit molecular tuning: specific amino acid substitutions improve enzyme performance under low-temperature conditions relevant to the deep sea, while selective expansions of gene families related to neural excitability, membrane systems, and other functions may help maintain physiological performance in this environment. Functional assays in yeast indicate enhanced performance of the deep-sea variants at 4°C. These results define a "loss-optimization-innovation" triad that underlies bathymetric adaptations and may apply to other deep-sea fauna worldwide.

Actiniaria↗
PubMed2026

Unraveling the genomic blueprint of the Indian black soldier fly: From genome assembly to evolutionary insights.

The black soldier fly (BSF) (Hermetia illucens) has been renowned for its sustainable bioconversion capabilities, resulting in smart protein production with wide applications in animal feed, bioenergy, and biofertilizer. However, the genetic mechanisms underlying efficient bioconversion and productivity remain poorly understood. To advance strain-specific applications and strengthen genetic resource availability, we present the whole genome sequencing (WGS) data for an Indian isolate of black soldier fly. The assembled genome was 1.46 Gb with a scaffold N50 of 172.7 Mb, and a GC content of 42.6%. Furthermore, 64.17% of genomic sequences were masked as repeated, and 14,317 protein-coding sequences were identified. Variant analysis against the reference genome identified 34.44 million variants (∼33.25 million SNPs and ∼ 1.18 million INDELs), with the majority (99.3%) classified as MODIFIER, 0.54% as LOW impact, 0.14% as MODERATE, and only 0.003% as HIGH impact. Comparative genomic analysis with other related species revealed expansions of gene families in BSF associated with Immune effector (Antimicrobial peptides (AMPs), Lysozymes, and Peptidoglycan Recognition Protein (PGRP) and Detoxification (cytochrome P450 enzymes). Notably, AMPs in the Indian isolate showed enhanced copy number variation in defensin (27) and PGRP (40) compared to reference BSF, suggesting potential regional adaptations to pathogen exposure. Collectively, this genomic data provides an improved resource for evolutionary studies, functional genomics, and targeted genetic improvement of BSF for sustainable bioconversion applications.

Comparative genomics↗
PubMed2026

Revisiting the genome assembly of Lupinus species reveals differential diploidization after a shared whole-genome duplication.

Accurate genome assemblies are essential for comparative genomics, yet Hi-C-guided scaffolding can introduce structural errors that misrepresent chromosome architecture and bias evolutionary inferences. Here, we identified pervasive scaffolding errors-including artificial fusions, internal inversions, and incomplete contig mounting-in 2 previously published Lupinus genomes (L. cosentinii and L. digitatus) using a segmentation method based on long terminal repeat (LTR) retrotransposon density. We reassembled both genomes, producing chromosome-level references of 472.7 Mb (16 chromosomes) and 427.2 Mb (21 chromosomes), with BUSCO completeness >98.5%. Synteny validation and reapplication of LTR profiling confirmed that all prior errors were resolved. Using these corrected genomes together with 4 additional Lupinus species and 2 outgroup legumes, we investigated postpolyploid evolution. Synonymous substitution rate (Ks) analysis revealed a genus-specific whole-genome duplication (WGD) event (Ks = 0.17) shared by all 6 Lupinus species. The proportion of WGD-derived genes varied markedly, from 60% in L. digitatus to only 36% in L. mutabilis, indicating differential diploidization. While all species retained a core set of WGD duplicates enriched in cytoskeleton organization, ion transport, and defense responses, each exhibited lineage-specific functional trajectories: cell wall modification in L. cosentinii and L. digitatus, nitrogen metabolism in L. albus and L. angustifolius, flower development in L. luteus, and stress/lipid metabolism in L. mutabilis. Our corrected assemblies provide optimal references for Lupinus comparative genomics, and our findings demonstrate that a shared WGD event can lead to both conserved and highly divergent postpolyploid fates, likely underpinning adaptive diversification within the genus.

Lupinus↗
PubMed2026

Chromosome-level genome assembly of the bitterling Rhodeus sinensis (Acheilognathidae) reveals genomic signatures associated with its mussel-dependent reproductive system.

Bitterlings (Acheilognathidae) exhibit a unique reproductive strategy characterized by symbiotic embryonic development inside the gill cavities of freshwater unionid mussels. Despite extensive ecological and physiological research on this system, genomic resources for bitterlings have remained limited, hindering comparative and evolutionary studies. Here, we present a high-quality, chromosome-level genome assembly for Rhodeus sinensis, a widely distributed bitterling species in the Korean Peninsula. By combining PacBio Continuous Long Read (CLR) sequencing, Illumina short reads, and Hi-C scaffolding, we generated a 0.77 Gb genome assembly with a scaffold N50 of 30.06 Mb. The final assembly comprises 24 chromosome-scale scaffolds, accounting for 98.3% of the assembled genome, with a BUSCO completeness score of 96.3% against the Actinopterygii_odb10. Comparative genomic analyses identified prominent expansions in gene families associated with alcohol metabolism, lipid catabolism, and oxidative stress responses. These genomic signatures of metabolic rewiring suggest a potential fuel flexibility, which may serve as a critical adaptive mechanism to mitigate the severe hypoxic stress encountered within the host mussel's gill environment. Ultimately, our chromosome-level genome assembly and findings provide a robust genomic foundation, contributing to a deeper understanding of the extreme physiological adaptations and unique life-history evolution within the Acheilognathidae.

Rhodeus sinensis↗
PubMed2025

Genomic signatures of innovation and selection in the extremotolerant yeast Kluyveromyces marxianus.

Extremophiles can be the product of millions of years of evolutionary engineering and refinement. The underlying mechanisms can be quite distinct from the ones operating at earlier stages of trait innovation. In this work, we have developed the compost yeast Kluyveromyces marxianus, which diverged from its closest relative >20 million years ago, as a model for interspecies comparative biology and genomics. We applied a battery of growth assays to species of the Kluyveromyces genus and found that K. marxianus outperformed its relatives in a battery of heat and chemical stress conditions. We then generated and analyzed genomes from across the genus, to find derived genetic features associated with, and potentially causal for, K. marxianus traits. We found robust expansions in gene families in the K. marxianus genome, most notably among genes annotated as transmembrane transporters and in metabolism. In molecular-evolution tests, we identified adaptive protein variants at hundreds of genes, among which plasma membrane transporters were over-represented. Together, these signals enable a model for the molecular mechanisms and evolutionary pressures underlying K. marxianus traits, including gains in transporter function mediating stress resistance, and metabolic variants contributing to its capacity for rapid growth in challenging conditions. Such oligogenic architectures may be the rule rather than the exception in phenotypes that have evolved over long timescales.

Journal Article↗
PubMed2025

The T2T genome assembly of watershield (Brasenia schreberi) unveils genomic insights into aquatic adaptation.

Watershield (Brasenia schreberi), belonging to Cabombaceae within the order Nymphaeales, represents one of the early-diverged angiosperm lineages. This perennial floating leaf freshwater aquatic plant features submerged juvenile leaves enveloped in a thick layer of transparent gelatinous mucilage, aiding in its resistance to aquatic stress. However, the evolutionary history of the mechanisms underlying its specific phenotype remains unclear. In this study, we present the telomere-to-telomere level genome of B. schreberi, unveiling that it underwent two rounds of whole-genome duplications (WGDs) and a recent whole-genome triplication, with the most ancient WGD being shared by Nymphaeaceae. WGD and dispersed duplication significantly contributed to the expansion of gene families, which are primarily associated with environmental adaptation. Additionally, we discovered that mature leaves primarily conduct photosynthesis and may transport nutrients to underwater juvenile leaves for polysaccharide synthesis. We also identified an ancestral broad expression pattern of ABC genes, and the similar expression of anthocyanin biosynthesis genes across all flower organs resulted in entirely purple flowers. Our findings deepen the understanding of the evolution of this specific aquatic plant phenotypes.

Genome, Plant↗
PubMed2026

Comparative genomics and full-length transcriptome profiling of wing morphs in Tetrix grossus (Orthoptera: Tetrigidae).

Wing polymorphism represents a paradigmatic dispersal-reproduction trade-off, yet its molecular basis remains uncharacterised in the phylogenetically distant pygmy grasshoppers (Tetrigidae). Here we integrate comparative genomics across ten orthopteran species with full-length transcriptomics of long-winged (FL) and short-winged (FS) Tetrix grossus. OrthoFinder recovered 118 orthogroups specific to T. grossus. Against a backdrop of pronounced gene-family contraction (36 expansions versus 222 contractions; net -186, mirrored at the ancestral Tetrix node, +37/-140), we identified an ancestral, Tetrix-specific expansion of hormone-regulation (12 genes; fold enrichment 7.93) and lipid/carbohydrate-metabolic families organised into syntenic clusters, alongside 513 positively selected genes enriched for integrin-mediated cell adhesion (6 genes), a process relevant to epithelial and appendage morphogenesis. Full-length transcriptomics of one long-winged (FL) and one short-winged (FS) adult female detected 7530 (FL) and 7515 (FS) expressed genes, with 794 FL- and 776 FS-restricted transcriptome-derived SNP-associated genes. The FL morph was enriched for an EGFR/Ras-Rho developmental-patterning axis and neuromuscular flight genes, whereas the FS morph was enriched for insulin/peptide-hormone response and growth-regulatory loci. Overall, we present genomic resources and testable hypotheses concerning the evolution and regulation of wing morphs in Tetrigidae rather than a validated genetic architecture of wing-morph determination.

Animals↗
PubMed2026

Evolutionary expansion of the NF-Y gene family in bivalves and divergent subunit responses to thermal and pathogenic stress in the noble scallop.

Nuclear factor Y (NF-Y) is a conserved eukaryotic transcription factor complex that specifically interacts with the CCAAT motif. Prior research has demonstrated that this gene family participates in various biological processes, encompassing growth, development, and stress responses, across a broad spectrum of organisms. However, research on the role of the NF-Y family in bivalves remains limited. In this study, we comprehensively identified the NF-Y family in 34 bivalve species, and further investigated its expression in the noble scallop Chlamys nobilis. A total of 296 NF-Y genes were identified and classified into three subfamilies, NF-YA, NF-YB, and NF-YC. Phylogenetic analysis revealed that NF-YA and NF-YC have remained relatively conserved, whereas NF-YB has undergone significant expansion. Additionally, while substantial disparities in gene copy numbers exist across species, the motif composition and exon-intron structures within each subfamily demonstrate notable conservation. Tissue expression profiling revealed distinct expression patterns among CnNF-Y genes, with several members exhibiting relatively high transcript abundance in gonadal tissues. Furthermore, qRT-PCR results demonstrated that CnNF-YA2, CnNF-YB6, and CnNF-YC were significantly and continuously upregulated under heat stress. Conversely, several genes, particularly CnNF-YA2, CnNF-YB3, and CnNF-YB4, exhibited dynamic transcriptional responses to Vibrio parahaemolyticus exposure. These findings enhance our understanding of the evolutionary trajectory and functional diversification of the NF-Y gene family in bivalves, laying a theoretical foundation for future research on thermal adaptation, immune regulation, and molecular breeding in scallops.

Animals↗
PubMed2026

Evolutionary adaptation of CCD4 enzymes in Buddleja alternifolia for crocetin biosynthesis.

INTRODUCTION: Carotenoid cleavage dioxygenase 4 (CCD4) enzymes play central roles in carotenoid turnover and apocarotenoid biosynthesis in plants. Despite their importance, the evolutionary mechanisms underlying diversification of CCD4 catalytic functions remain poorly understood. This study investigated the CCD gene family in Buddleja alternifolia, with particular emphasis on the expansion and functional evolution of the CCD4 subfamily. METHODS: A genome-wide identification and comparative analysis of CCD genes were performed in B. alternifolia. Genomic organization, phylogenetic relationships, and syntenic patterns were analyzed to investigate gene family expansion. Functional characterization of 11 BaCCD4 paralogs was conducted through biochemical assays, while structural analyses were used to identify sequence features associated with differences in substrate cleavage specificity. Gene expression profiling was performed to assess patterns of tissue-specific regulation. RESULTS: Twenty-three CCD genes were identified, including 12 CCD4 paralogs, representing one of the largest CCD4 expansions reported within Lamiales. Syntenic and genomic analyses revealed that recent tandem duplication events, particularly within a CCD4-rich region on chromosome 10, were the primary drivers of this expansion. The presence of pseudogenes in the same region supported an ongoing birth-and-death evolutionary process. Functional analyses demonstrated extensive biochemical diversification among BaCCD4 enzymes despite their high sequence similarity. Several paralogs catalyzed asymmetric carotenoid cleavage leading to citraurin production, whereas two paralogs, KAG8367281 and KAG8375220, exhibited symmetric zeaxanthin cleavage activity, producing crocetin dialdehyde, the direct precursor of crocins. Notably, these crocetin-producing enzymes belonged to closely related paralogous pairs whose counterparts displayed distinct cleavage specificities, indicating rapid neofunctionalization after duplication. Structural analyses suggested that subtle sequence variations, including indels affecting loop regions adjacent to the substrate access channel, may underlie changes in regioselectivity. Expression profiling further revealed tissue-specific expression patterns consistent with functional divergence among paralogs. DISCUSSION: These findings indicate that crocetin-forming activity in B. alternifolia likely evolved through progressive modifications of ancestral CCD4 functions rather than through a single evolutionary event. The remarkable expansion and diversification of the CCD4 subfamily provide evidence for the role of gene duplication and neofunctionalization in shaping carotenoid cleavage specificity. Collectively, this work establishes B. alternifolia as a valuable model for investigating the molecular evolution of CCD4 enzymes and the emergence of specialized apocarotenoid metabolism in plants.

CCD4 evolution↗
PubMed2026

Adaptive genomic evolution and WD40-regulated temporal dynamics of anthocyanins support leaf photoplasticity in Parrotia subaequalis.

BACKGROUND: Parrotia subaequalis, a Tertiary relict endemic to China, plays a significant role in phylogeny and adaptive evolution as a key species in the early differentiation of angiosperms. It has abundant leaf colors and great potential as an ornamental tree. RESULTS: This study assembled the first chromosome-level genome of P. subaequalis (Contig N50 = 2.15 Mb), revealing transposable element proliferation, key paleopolyploid events and dynamic gene family evolution, including the expansion of secondary metabolite transport and synthesis genes (such as WD40, 2OG-FeII_Oxy) and the contraction of gene families related to flower morphogenesis (such as F-box-like, K-box). Through integrative transcriptomics and targeted metabolomics approaches, we further revealed that the color transition of young leaves from red to green was driven by temporal accumulation differences of malvidin-3,5-O-diglucoside, whose biosynthesis is progressively down-regulated during leaf development. WGCNA revealed that a subset of WD40 genes (light-signaling, TTG1/HOS15-like, etc.) coexpresses with anthocyanin biosynthetic genes, like 4CLL9, GT1, in anthocyanin-related modules enriched for auxin signaling and hydrolase activity, suggesting a potential link between WD40 expansion and photoprotective plasticity. Relevant regulatory networks were found to complement the species-specific gene pool related to leaf color regulation. CONCLUSION: This genomic resource of P. subaequalis advanced our understanding of early angiosperm adaptation through neofunctionalized regulatory networks and established a foundation for molecular breeding aimed at enhancing environmental resilience while preserving ornamental traits.

Anthocyanins↗
PubMed2025

Multi-omics analyses provide insights into the molecular basis for salt tolerance of Phyla nodiflora.

The perennial herbaceous plant, Phyla nodiflora (Verbenaceae), which possesses natural resistance to multiple abiotic stresses, is widely used as a pioneer species in island ecological restoration. Due to the lack of information about its genome, the mechanism underlying its tolerance to environmental stresses, such as salinity, is almost entirely unknown. Here, we report on the high-quality genome of P. nodiflora that is 403.07 Mb in size, and which was assembled and anchored onto 18 pseudo-chromosomes. Genomic synteny revealed that P. nodiflora underwent two whole genome duplication events, which promoted the expansion of genes related to environmental adaptation and the biosynthesis of secondary metabolites. An integrated genomic and transcriptomic analysis suggested that salt stress tolerance in P. nodiflora is associated with the expansion and activated expression of genes related to abscisic acid (ABA) homeostasis and signaling. The expansion of ZEP family genes may contribute to the consistent increase in ABA levels under salt stress. Lysine acetylomic analysis revealed that exposure to salt led to widespread protein deacetylation, with these proteins primarily involved in signal transduction, carbohydrate transport and metabolism, and transcription regulation. Deacetylation of glutathione S-transferase increased enzymatic activities in response to salt-induced oxidative stress. Collectively, the genomic, transcriptomic, and lysine acetylomic analyses provide profound insight into the molecular basis of the adaptation of P. nodiflora to salt stress, and will be helpful to engineer salt-tolerant plants for ecological restoration.

Salt Tolerance↗