PubMed HealthSearch

SEARCH · PubMed Health

Results for “gene evolution”

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

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

At least 19 recordsLinked to original sources

Genome-wide cyclin gene evolution in Arabidopsis and Brassica reveals polyploidization-driven duplication and flowering-time associations.

Cyclin genes are plant cell cycle regulators that play essential roles in growth, development, and reproduction. However, the evolutionary dynamics and genomic organization of cyclin genes across the Brassicaceae family remain poorly understood, particularly in the context of allotetraploid genome evolution. Here, we investigated the diversity, expansion mechanisms, and potential functional diversification of cyclin genes across ten Brassicaceae genomes, including four Arabidopsis and six Brassica species. A total of 1087 cyclin genes representing 23 cyclin types were identified. Comparative genomic analyses revealed that cyclin gene expansion was strongly influenced by polyploidization in Brassica species, with 1845 duplication events involving 1063 genes. Whole-genome duplication was the predominant mechanism driving expansion, while both inter- and intra-genomic duplications contributed to gene retention in tetraploid Brassica species, with the highest duplication frequency observed in Brassica juncea. Across genomes, 120 physical gene clusters were identified, including homogeneous and heterogeneous types. Ortholog analysis between progenitor and allotetraploid species identified 852 orthologous pairs involving 366 genes, indicating extensive conservation following allotetraploid formation. Phylogenetic analysis resolved cyclins into three major clades, while expression-based clustering in Brassica napus grouped genes into four major clusters, suggesting functional diversification. Integration of pan-genomic and flowering-time QTL analyses further identified two cyclin genes, Bna21cycA2 and Bna113cycD4, which contain amino acid polymorphisms and represent putative candidate variations potentially associated with flowering-time variation across multiple genomes. These findings provide new insights into the evolutionary expansion, retention, and potential functional divergence of cyclin genes in Brassicaceae and highlight candidate loci for future functional studies and crop improvement.

Evolution, Molecular

Recent gene duplication and structural remodeling drive rapid lineage-specific gene family evolution in plants.

Gene duplication promotes the generation of novel gene functions and trait diversity across species. Here, we present DupHIST, a computational pipeline that reconstructs the hierarchical timing of gene duplications by integrating maximum likelihood (ML)-based phylogeny with substitution-derived timing via statistical smoothing. Applied to over 4.5 million genes from 114 plant genomes, we successfully inferred duplication histories across nearly 130,000 orthogroups. This large-scale analysis showed that 53.0% of genes arose from recent, lineage-specific duplications, with high concentrations in particular multi-copy families. Among these, NLR, C48, and P450 families exemplified how recently duplicated genes undergo rapid stepwise structural remodeling. This process was primarily driven by small-scale mutations, including insertions, deletions, and frameshifts, that rapidly accumulated shortly after duplication. By resolving the precise duplication order, we reconstructed these architectural changes, thereby enabling both the inference of putative ancestral structures and the exploration of functional diversification arising from structural remodeling. Structure-based clustering further uncovered that recently duplicated, uncharacterized genes retain core domain structures resembling known functional proteins even across phylogenetically distant species lacking sequence homology. Our findings reveal that recent gene duplications and subsequent structural remodeling represent a widespread and lineage-specific force driving rapid diversification of gene families in plants.

Gene duplication history

Gene transfer agents: The ambiguous role of selfless viruses in genetic exchange and bacterial evolution.

Gene transfer agents (GTAs) are genetic elements derived from ancestral bacteriophages that have become domesticated by the host. GTAs are present in diverse prokaryotic organisms, where they can facilitate horizontal gene transfer under certain conditions. Unlike typical bacteriophages, GTAs do not exhibit any preference for the replication or transfer of the genes encoding them; instead, they exhibit a remarkable capacity to package chromosomal, and sometimes extrachromosomal, DNA into virus-like capsids and disseminate it to neighboring cells. Because GTAs resemble defective prophages, identification of novel GTAs is not trivial. The detection of candidates relies on the genetic similarity to known GTAs, which has been fruitful in α-proteobacterial lineages but challenging in more distant bacteria. Here we consider several fundamental questions: What is the true prevalence of GTAs in prokaryote genomes? Given there are high costs for GTA production, what advantage do GTAs provide to the bacterial host to justify their maintenance? How is the bacterial chromosome recognized and processed for inclusion in GTA particles? This article highlights the challenges in comprehensively understanding GTAs' prevalence, function and DNA packaging method. Going forward, broad study of atypical GTAs and use of ecologically relevant conditions are required to uncover their true impact on bacterial chromosome evolution.

Gene Transfer, Horizontal

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

Adaptive Evolution of the PFK Gene Family in Chinese Longsnout Catfish, Leiocassis longirostris.

The Chinese longsnout catfish is a typical carnivorous fish with a relatively weak ability to utilize glucose. However, the genomic basis for its glucose metabolic adaptation remains unclear. In this study, we used comparative genomics methods to systematically analyze the evolutionary characteristics of glucose metabolism-related genes in the Chinese longsnout catfish, focusing on gene family evolution, patterns of expansion and contraction, and selective pressures. The results indicate that glucose metabolism-related genes have undergone significant reshaping during evolution. Genes involved in glucose digestion, absorption, and insulin signaling pathways demonstrate a tendency toward contraction, while those associated with protein and lipid metabolism exhibit expansion. This pattern is consistent with the species' long-term adaptation to a high-protein, high-fat diet. Comparative analysis further revealed that, compared to fish with different dietary habits, certain key genes involved in glycolysis in the Chinese longsnout catfish exhibit a reduction in copy number. Molecular evolutionary analysis showed that key genes involved in glycolysis and gluconeogenesis (including hexokinase 2 (hk2), phosphofructokinase, muscle/platelet (pfkm/p)) exhibit signs of accelerated evolution or positive selection. Notably, the PFK gene family exhibits complex evolutionary characteristics resulting from the combined effects of gene contraction, rapid evolution, and positive selection. In summary, this study reveals the genomic evolutionary basis for the glucose metabolic adaptation of the Chinese longsnout catfish and identifies the PFK gene family as a key candidate for elucidating its unique glucose metabolic characteristics.

Leiocassis longirostris

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

A comprehensive phylogeny of mammalian PRNP gene reveals no influence of prion misfolding propensity on the evolution of this gene.

Prion diseases are invariably fatal neurodegenerative diseases that affect some mammalian species, including humans. These diseases are caused by the misfolding of the cellular prion protein (PrPC) into a pathologic isoform (PrPSc). The prion protein is highly conserved across mammals. However, some species present lower susceptibility to prion diseases than others. This behavior is likely explained by the resistance of these animal species' prion proteins to acquire a pathological conformation. Therefore, the tertiary structure and interspecific variations encoded in the primary structure determine a PrP proneness to misfolding. For this reason, we studied the PRNP gene from a phylogenetic perspective, potentially unveiling evolutionary events related to prion diseases. We generated a database of mammalian PRNP sequences and constructed phylogenetic trees based on nucleotide sequence variations. We aligned 1146 PRNP gene sequences from 901 different mammalian species and built a PRNP gene-based phylogenetic tree. Classical phylogenetic orders tend to maintain their clustering in the PRNP gene tree. Nonetheless, the few differences found may shed some light on potential evolutionary constraints posed by prion disorders. Moreover, this phylogenetic study was combined with an in vitro misfolding study. Protein Misfolding Shaking Amplification (PMSA) was used to evaluate the tendency of many of these proteins to misfold. This comprehensive analysis spanned a wide range of mammalian prion protein sequences and included analysis of different variants with a focus on the human rs1799990 locus (c.385A > G, p.Met129Val). This variant, widely linked to prion disease susceptibility in humans, is explored in the context of its evolutionary origins. All in all, our PRNP gene-based tree, despite showing some topological differences with the reference species tree that could be in some cases related to prion disease susceptibility, is not significantly distinct. Indicating that the proneness of a PrP variant to misfold spontaneously has not shaped the evolution of this gene.

Phylogeny

Diversification, loss, and virulence gains of the major effector AvrStb6 during continental spread of the wheat pathogen Zymoseptoria tritici.

Interactions between plant pathogens and their hosts are highly dynamic and mainly driven by pathogen effectors and plant receptors. Host-pathogen co-evolution can cause rapid diversification or loss of pathogen genes encoding host-exposed proteins. The molecular mechanisms that underpin such sequence dynamics remains poorly investigated at the scale of entire pathogen species. Here, we focus on AvrStb6, a major effector of the global wheat pathogen Zymoseptoria tritici, evolving in response to the cognate receptor Stb6, a resistance widely deployed in wheat. We comprehensively captured effector gene evolution by analyzing a global thousand-genome panel using reference-free sequence analyses. We found that AvrStb6 has diversified into 59 protein isoforms with a strong association to the pathogen spreading to new continents. Across Europe, we found the strongest differentiation of the effector consistent with high rates of Stb6 deployment. The AvrStb6 locus showed also a remarkable diversification in transposable element content with specific expansion patterns across the globe. We detected AvrStb6 gene losses and evidence for transposable element-mediated disruptions. We used virulence datasets of genome-wide association mapping studies to predict virulence changes across the global panel. Genomic predictions suggested marked increases in virulence on Stb6 cultivars concomitant with the spread of the pathogen to Europe and the subsequent spread to further continents. Finally, we genotyped French bread wheat cultivars for Stb6 and monitored resistant cultivar deployment concomitant with AvrStb6 evolution. Taken together, our data provides a comprehensive view of how a rapidly diversifying effector locus can undergo large-scale sequence changes concomitant with gains in virulence on resistant cultivars. The analyses highlight also the need for large-scale pathogen sequencing panels to assess the durability of resistance genes and improve the sustainability of deployment strategies.

Ascomycota

Conserved function of medaka pink-eyed dilution in melanin synthesis and its divergent transcriptional regulation in gonads among vertebrates.

Medaka is emerging as a model organism for the study of vertebrate development and genetics, and its effectiveness in forward genetics should prove equal to that of zebrafish. Here, we identify by positional cloning a gene responsible for the medaka i-3 albino mutant. i-3 larvae have weakly tyrosinase-positive cells but lack strongly positive and dendritic cells, suggesting loss of fully differentiated melanophores. The region surrounding the i-3 locus is syntenic to human 19p13, but a BAC clone covering the i-3 locus contained orthologs located at 15q11-13, including OCA2 (P). Medaka P consists of 842 amino acids and shares approximately 65% identity with mammalian P proteins. The i-3 mutation is a four-base deletion in exon 13, which causes a frameshift and truncation of the protein. We detected medaka P transcripts in melanin-producing eyeballs and (putative) skin melanophores on embryos and an alternatively spliced form in the non-melanin-producing ovary or oocytes. The mouse p is similarly expressed in gonads, but not alternatively spliced. This is the first isolation of nonmammalian P, the functional mechanism of action of which has not yet been elucidated, even in mammals. Further investigation of the functions of P proteins and the regulation of their expression will provide new insight into body color determination and gene evolution.

Amino Acid Sequence

Genome-wide identification, characterization, and expression pattern analysis of the glyoxalase gene family in Phyllostachys pubescens during abiotic stresses.

BACKGROUND: The glyoxalase pathway comprising of three enzymes i.e., glyoxalase I (GLYI), glyoxalase II (GLYII), and glyoxalase III (GLYIII), which play vital role in mitigating abiotic stresses by detoxifying the stress induced cytotoxic metabolite methylglyoxal (MG). Phyllostachys pubescens an ecologically and economically important forest species, plays vital roles in carbon sequestration and climate change mitigation. A genome-wide study was conducted to identify and characterize GLYI, GLYII, and unique DJ-1/GLYIII gene candidates in P. pubescens. The identified members were evaluated based on phylogenetic analysis, gene structure, chromosomal distribution, gene duplication, presence of conserved domain(s) and cis regulatory region. RESULTS: A total of 19 GLYI, 18 GLYII, and 15 GLYIII members were identified, each featuring characteristic domains: glyoxalase, metallo-β-lactamase, and DJ-1/PfpI, respectively. The presence of different cis-elements in the promoter region of the glyoxalase genes gives insights into their role and regulation under hormonal response, developmental processes and stress adaptation. Besides this, stress responsive transcription factors binding sites also dominated the promoter regions of glyoxalase genes. Expression analysis of various glyoxalase genes demonstrated significant variability under different stress conditions, underscoring their potential roles in stress modulation. Significant upregulation of all of the PhGLYI, PhGLYII, and PhGLYIII were observed under cold, drought, heavy metal and salinity stress suggesting their involvement in oxidative stress management, osmotic regulation and remodelling cellular redox homeostasis. Among the glyoxalase genes, PhGLYI-15, PhGLYII-9, and PhGLYIII-3 showed consistent upregulation under various abiotic stresses. CONCLUSIONS: Our findings reveal that glyoxalase genes crucially contribute towards the improvement of cellular osmotic potential in moso bamboo under different abiotic stresses. This study enhances our understanding of glyoxalase genes' evolution and functional roles in plants and opens new avenues for developing stress resilient crop varieties for sustainable agriculture.

Lactoylglutathione Lyase

Fishing for a reelGene: evaluating gene models with evolution and machine learning.

Assembled genomes and their associated annotations have transformed our study of gene function. However, each new annotated assembly generates new gene models. Inconsistencies between annotations likely arise from biological and technical causes, including pseudogene misclassification, transposon activity, and intron retention from sequencing of unspliced transcripts. To evaluate gene model predictions, we developed reelGene, a pipeline of machine learning models focused on (1) transcription boundaries, (2) mRNA integrity, and (3) protein structure. The first two models leverage sequence characteristics and evolutionary conservation across related taxa to learn the grammar of conserved transcription boundaries and mRNA sequences, while the third uses the conserved evolutionary grammar of protein sequences to predict whether a gene can produce a protein. Evaluating 1.8 million transcript models in Zea mays ssp. mays (maize), reelGene classified 28% as incorrectly annotated or non-functional. We find that reelGene classifies 92.2% of genes in the maize proteome and 99.2% of genes within the maize classical gene list as functional. reelGene also provides a way to further investigate genome biology- for instance, reelGene indicates that 10.3% of dispensable genes in B73 are functional, and within retained duplicate genes, reelGene identifies a 30% bias toward the retention of the M1 subgenome when one copy is functional and the other is non-functional. As an annotation-evaluating tool, reelGene is directly applicable to species of the Andropogoneae tribe, including other important crops like sorghum and miscanthus. As a community resource, reelGene has been integrated onto MaizeGDB both as a browser track and as an individual Shiny App, allowing researchers to evaluate gene model accuracy and further investigate genome biology.

Machine Learning

Intra-colony divergence and global allele sharing reflect purifying selection and recombination at the Botryllus histocompatibility factor locus.

Urochordates, the closest relatives of vertebrates, lack adaptive immunity. However, some taxa, such as the colonial species Botryllus schlosseri, provide a unique model for studying innate self/non-self recognition through natural allogeneic transplantation responses. In this species, interactions between colonies are controlled by a highly polymorphic locus, with the Botryllus histocompatibility factor (BHF) being the only gene known to predict tissue fusion or rejection outcomes with complete accuracy. Here, we analyzed full-length BHF alleles from 19 laboratory-born and wild colonies and found that highly divergent alleles tend to coexist within individuals, whereas identical alleles can be shared across continental-scale distances. Despite extensive length variation, evidence of intragenic recombination, and pronounced nucleotide diversity, BHF exhibits limited protein divergence, with 33 alleles encoding only 17 distinct protein variants. Contrary to expectations for polymorphic recognition genes, no evidence of balancing or directional selection was detected. Instead, signatures of purifying selection were observed. We propose that this contrast between nucleotide and protein diversity arises from the combined effects of recombination, human-mediated gene flow, and linkage to nearby loci under balancing selection, while functional constraints maintain protein stability. These findings suggest that extensive protein diversification may not be a universal driver of allorecognition gene evolution.

Animals

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)

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

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

Evolutionary architecture and lineage-specific diversification of Forkhead box transcription factors in Perna viridis.

The Forkhead box (Fox) transcription factors are evolutionarily conserved regulators of development, cell cycle, and apoptosis across metazoans. This study provides the first comprehensive genome-wide analysis of the Fox gene family in the Asian green mussel (Perna viridis). We identified 28 Fox genes distributed across 10 chromosomes. Comparative analysis reveals the absence of the FoxI, FoxQ1, FoxR and FoxS subfamily, consistent with other bivalves and indicative of lineage-specific gene loss during molluscan evolution. Notably, gene duplications in the FoxAB, FoxD, FoxH, FoxN1-4, FoxQ2 and FoxQD subfamilies may reflect functional diversification associated with environmental adaptation. Exon-intron structural variability, including intron loss in several paralogues, suggests structural diversification and potential regulatory variation. Phylogenetic reconstruction confirmed the monophyly of core Fox classes while highlighting divergent expansion patterns in lophotrochozoans. Selection analyses showed strong purifying selection across duplicated Fox paralogs, supporting functional conservation after lineage-specific expansion. Gene Ontology enrichment linked Fox genes to stress response, apoptosis, and transcriptional regulation. By integrating phylogenetic, structural, and transcriptomic analyses, this study provides a genomic framework for understanding Fox gene organisation, evolution, and tissue-associated expression patterns in Perna viridis and establishes a comparative resource for future functional studies in bivalves.

Animals

Comprehensive identification and analysis of clusters of tandemly duplicated genes reveal their contributions to adaptive evolution of green plants.

Tandem gene duplication occurred more frequently compared with the episodic whole-genome duplication (WGD), providing a continuous supply of genetic material for evolutionary innovation and adaptation to changing environments. The rising roles of clusters of tandemly duplicated genes (CTDGs) in the evolution of phenotypic diversity have been unraveled in mammals. However, the content and biological roles of CTDGs remain largely unknown in plants. Here, we comprehensively identified CTDGs in 220 published plant genomes representing major lineages of green plants. The number of CTDGs showed great variation across taxa, ranging from 0 to 6028. The size of CTDGs varied from 2 to 47 genes, with small clusters containing two members predominating. Interestingly, significant expansion of CTDGs was found in early-diverging land plants and is closely associated with the evolution of key traits (e.g., ABA response, plant cuticle, UV-B resistance) required for plants to conquer terrestrial environments. Functional enrichment analysis revealed conserved and specialized functional profiles among different sizes of CTDGs in both Arabidopsis thaliana and the bryophyte Physcomitrium patens. Small CTDGs were enriched in fundamental stress responses, including protein modification, signal transduction, and responses to diverse stress stimuli, while large CTDGs were enriched in more sophisticated processes such as plant hormone biosynthesis and signaling, plant-microbe interactions, and reproductive processes. Expression pattern analyses of CTDGs under different stress conditions in A. thaliana and P. patens revealed that the highest number of CTDGs showed differential expression under drought stress, suggesting important roles of CTDGs in the evolution of desiccation tolerance in early land plants. The results of this study provide new additions to our knowledge about the abundance of CTDGs across green plants and reveal their important contributions to enable plants to overcome stressful environments on land.

Gene Duplication

Positive selection and relaxed purifying selection contribute to rapid evolution of sex-biased genes in green seaweed Ulva.

BACKGROUND: The evolution of differences in gamete size and number between sexes is a cornerstone of sexual selection theories. The green macroalga Ulva, with incipient anisogamy and parthenogenetic gametes, provides a unique system to investigate theoretical predictions regarding the evolutionary pressures that drive the transition from isogamy to anisogamy, particularly in relation to gamete size differentiation and sexual selection. Its minimal gamete dimorphism and facultative parthenogenesis enable a rare window into early evolutionary steps toward anisogamy. RESULTS: By analyzing the expression profiles of sex-biased genes (SBGs) during gametogenesis, we found that SBGs evolve faster than unbiased genes, driven by higher rates of non-synonymous substitution (dN), indicating that SBGs are under stronger selective pressures. Mating type minus-biased genes (mt-BGs) exhibit higher dN/dS values than mating type plus-biased genes (mt+BGs), suggesting stronger selective pressures on mt-BGs, although this difference was not statistically significant (P = 0.08). Using branch-site and RELAX models, we found positive selection and relaxed purifying selection acting on a significant proportion of SBGs, particularly those associated with flagella function. CONCLUSIONS: This study highlights the selective pressures shaping anisogamy and provides insights into the molecular mechanisms underlying its evolution. The faster evolution of SBGs, particularly mt-BGs, and the positive selection on genes associated with motility, such as those related to flagella function, suggest the importance of enhanced gamete motility in the transition to anisogamy. These findings contribute to our understanding of sexual selection and the evolutionary forces that drive the differentiation of gamete size and number between sexes.

Selection, Genetic